Near-isogenic Lines Research Articles

Gene expression of near-isogenic lines (NILs) carrying blast resistance genes Pi9 and Pi54 in the background of rice cultivar Mushk Budji.

Kashmir valley, India is a homeland to rice landraces like Zag, Nunbeoul, Qadirbeigh, Kawkadur, Kamad, Mushk Budji, etc., generally characterized by short grains, aroma, earliness and cold tolerance. Mushk Budji is a commercially important speciality rice known for its taste and aroma, nonetheless, is extremely vulnerable to blast disease. Through the use of the marker-assisted backcrossing (MABC) approach, a set of 24 Near-isogenic lines (NILs) was created, and the lines with the highest background genome recovery were chosen. The expression analysis was carried out for the component genes and other eight pathway genes related to blast resistance. The major blast resistance genes Pi9 (from IRBL-9W) and Pi54 (from DHMAS 70Q 164-1b) were incorporated following simultaneous-but-step-wise MABC. The NILs harbouring genes Pi9 + Pi54, Pi9 and Pi54 expressed resistance to isolate (Mo-nwi-kash-32) under controlled and natural field conditions. The loci controlling ETI (effector triggered immunity) included the gene Pi9 and showed 61.18 and 60.27 fold change in relative gene expression in Pi54 + Pi9 and Pi9 carrying NILs against RP Mushk Budji. Pi54 was up regulated and showed 41 and 21 fold change in relative gene expression for NIL-Pi54 + Pi9 and NIL-Pi54, respectively. Among the pathway genes, LOC_Os01g60600 (WRKY 108) recorded 8 and 7.5 fold up regulation in Pi9 and Pi54 NILs. The NILs showed recurrent parent genome recovery (RPG) per cent of 81.67 to 92.54 and were on par in performance to recurrent parent Mushk Budji. The lines were utilized to study the expression of the loci controlling WRKYs, peroxidases and chitinases that confer overall ETI response.

Fine mapping of a Fusarium crown rot resistant locus on chromosome arm 6HL in barley by exploiting near isogenic lines, transcriptome profiling, and a large near isogenic line-derived population

Key messageThis study reported validation and fine mapping of a Fusarium crown rot resistant locus on chromosome arm 6HL in barley using near isogenic lines, transcriptome sequences, and a large near isogenic line-derived population.Fusarium crown rot (FCR), caused by Fusarium pseudograminearum, is a chronic and serious disease affecting cereal production in semi-arid regions globally. The increasing prevalence of this disease in recent years is attributed to the widespread adoption of minimum tillage and stubble retention practices. In the study reported here, we generated eight pairs of near isogenic lines (NILs) targeting a putative QTL (Qcrs.caf-6H) conferring FCR resistance in barley. Assessing the NILs confirmed the large effect of this locus. Aimed to develop markers that can be reliably used in incorporating this resistant allele into breeding programs and identify candidate genes, transcriptomic analyses were conducted against three of the NIL pairs and a large NIL-derived population consisting of 1085 F7 recombinant inbred lines generated. By analyzing the transcriptomic data and the fine mapping population, Qcrs.caf-6H was delineated into an interval of 0.9 cM covering a physical distance of ~ 547 kb. Six markers co-segregating with this locus were developed. Based on differential gene expression and SNP variations between the two isolines among the three NIL pairs, candidate genes underlying the resistance at this locus were detected. These results would improve the efficiency of incorporating the targeted locus into barley breeding programs and facilitate the cloning of causal gene(s) responsible for the resistance.

Transcriptome profiling reveals the expression and regulation of genes associated with Fusarium wilt resistance in chickpea (Cicer arietinum L.).

Fusarium wilt (FW) is one of the most significant biotic stresses limiting chickpea production worldwide. To dissect the molecular mechanism of FW resistance in chickpea, comparative transcriptome analyses of contrasting resistance sources of chickpea genotypes under control and Fusarium oxysporum f. sp. ciceris (Foc) inoculated conditions were performed. The high-throughput transcriptome sequencing generated about 1137 million sequencing reads from 24 samples representing two resistant genotypes, two susceptible genotypes, and two near-isogenic lines under control and stress conditions at two-time points (7th- and 12th-day post-inoculation). The analysis identified 5182 differentially expressed genes (DEGs) between different combinations of chickpea genotypes. Functional annotation of these genes indicated their involvement in various biological processes such as defense response, cell wall biogenesis, secondary metabolism, and disease resistance. A significant number (382) of transcription factor encoding genes exhibited differential expression patterns under stress. Further, a considerable number of the identified DEGs (287) co-localized with previously reported quantitative trait locus for FW resistance. Several resistance/susceptibility-related genes, such as SERINE/THREONINE PROTEIN KINASE, DIRIGENT, and MLO exhibiting contrasting expression patterns in resistant and susceptible genotypes upon Foc inoculation, were identified. The results presented in the study provide valuable insights into the transcriptional dynamics associated with FW stress response in chickpea and provide candidate genes for the development of disease-resistant chickpea cultivars.

Comparative transcriptome profiling of near isogenic lines PBW343 and FLW29 to unravel defense related genes and pathways contributing to stripe rust resistance in wheat.

Stripe rust (Sr), caused by Puccinia striiformis f. sp. tritici (Pst), is the most devastating disease that poses serious threat to the wheat-growing nations across the globe. Developing resistant cultivars is the most challenging aspect in wheat breeding. The function of resistance genes (R genes) and the mechanisms by which they influence plant-host interactions are poorly understood. In the present investigation, comparative transcriptome analysis was carried out by involving two near-isogenic lines (NILs) PBW343 and FLW29. The seedlings of both the genotypes were inoculated with Pst pathotype 46S119. In total, 1106 differentially expressed genes (DEGs) were identified at early stage of infection (12 hpi), whereas expressions of 877 and 1737 DEGs were observed at later stages (48 and 72 hpi) in FLW29. The identified DEGs were comprised of defense-related genes including putative R genes, 7 WRKY transcriptional factors, calcium, and hormonal signaling associated genes. Moreover, pathways involved in signaling of receptor kinases, G protein, and light showed higher expression in resistant cultivar and were common across different time points. Quantitative real-time PCR was used to further confirm the transcriptional expression of eight critical genes involved in plant defense mechanism against stripe rust. The information about genes are likely to improve our knowledge of the genetic mechanism that controls the stripe rust resistance in wheat, and data on resistance response-linked genes and pathways will be a significant resource for future research.

Mapping of the bs5 and bs6 non-race-specific recessive resistances against bacterial spot of pepper.

Bacterial spot caused by Xanthomonas euvesicatoria is a major disease of pepper (Capsicum annuum L.) in warm and humid production environments. Use of genetically resistant cultivars is an effective approach to manage bacterial spot. Two recessive resistance genes, bs5 and bs6, confer non-race-specific resistance against bacterial spot. The objective of our study was to map these two loci in the pepper genome. We used a genotyping-by-sequencing approach to initially map the position of the two resistances. Segregating populations for bs5 and bs6 were developed by crossing susceptible Early CalWonder (ECW) with near-isogenic lines ECW50R (bs5 introgression) or ECW60R (bs6 introgression). Following fine-mapping, bs5 was delimited to a ~535 Kbp interval on chromosome 3, and bs6 to a ~666 Kbp interval in chromosome 6. We identified 14 and 8 candidate resistance genes for bs5 and bs6, respectively, based on predicted protein coding polymorphisms between ECW and the corresponding resistant parent. This research enhances marker-assisted selection of bs5 and bs6 in breeding programs and is a crucial step towards elucidating the molecular mechanisms underlying the resistances.

Proteomic analysis of near-isogenic lines reveals key biomarkers on wheat chromosome 4B conferring drought tolerance.

Drought is a major constraint for wheat production that is receiving increased attention due to global climate change. This study conducted isobaric tags for relative and absolute quantitation proteomic analysis on near-isogenic lines to shed light on the underlying mechanism of qDSI.4B.1 quantitative trait loci (QTL) on the short arm of chromosome 4B conferring drought tolerance in wheat. Comparing tolerant with susceptible isolines, 41 differentially expressed proteins were identified to be responsible for drought tolerance with a p-value of < 0.05 and fold change >1.3 or <0.7. These proteins were mainly enriched in hydrogen peroxide metabolic activity, reactive oxygen species metabolic activity, photosynthetic activity, intracellular protein transport, cellular macromolecule localization, and response to oxidative stress. Prediction of protein interactions and pathways analysis revealed the interaction between transcription, translation, protein export, photosynthesis, and carbohydrate metabolism as the most important pathways responsible for drought tolerance. The five proteins, including 30S ribosomal protein S15, SRP54 domain-containing protein, auxin-repressed protein, serine hydroxymethyltransferase, and an uncharacterized protein with encoding genes on 4BS, were suggested as candidate proteins responsible for drought tolerance in qDSI.4B.1 QTL. The gene coding SRP54 protein was also one of the differentially expressed genes in our previous transcriptomic study.

Comparison of gene expression changes in two wheat varieties with different phenotype to strip rust using RNA-Seq analysis

The fungus Puccinia striiformis f. sp. tritici (Pst) is an important threat to wheat production because it can cause wheat stripe rust. The present study aimed to identify new stripe rust resistance genes and to provide a theoretical and practical basis for breeding wheat varieties with broad spectrum, stable, and durable resistance. Wheat leaves inoculated with wheat stripe rust fungus Chinese yellow rust 34 were collected at different time points for transcriptomic analysis based on the wheat stripe rust susceptible varieties AVOCET S (AVS) and AVSYr15NIL [near-isogenic line (NIL) derived from AVS]. The results showed that the number of upregulated genes in the two varieties was 294, 364, 398, and 604, and the number of downregulated genes was 520, 178, 570, and 345 on the 1st, 3rd, 5th, and 7th days post inoculation, respectively. Gene Ontology and Kyoto Encyclopedia of Gene and Genomes enrichment analyses found enrichment of differentially expressed genes in the peroxisome proliferators-activated receptor signaling pathways, plant–pathogen interaction, and styrene acrylic acid biosynthesis that encoded protein kinases, signal transduction, transcription factors, and functional protein components. Differentially expressed genes were randomly selected for quantitative reverse transcription PCR analysis, and the change trend was the same as in the transcriptome data. The results of this study suggest that genes in AVSYr15NIL related to the stripe rust response could be valuable for understanding the mechanisms involved in stripe rust resistance.

Comparison of cotton (Gossypium hirsutum) hybrids and their non-transgenic in India

It is still not clear whether differential response of Bt cotton and conventional non-Bt cotton genotypes in terms of growth, phenology and yield ability are hypothetical or realistic. This study investigated the response of three Bt cotton hybrids in comparison to their conventional non-transgenic version (near-isogenic lines) along with two locally adapted genotypes for growth, yield attributes, phenological characteristics and rooting behavior. Therefore, two experiments were conducted during summer season of year 2016 and 2017 to evaluate differences among 3 popular Bt cotton hybrids along with their non-Bt version. First field experiment was conducted as a pot study wherein shoot and root parameters of 3 Bt cotton hybrids (Ankur 3028, MRC 7017, Bioseed 6588) were compared with their non-Bt version along with 2 local checks (i.e.,hirsutum cotton Cv. F2228 and arboreum cottonCv. FDK124). Second field experiment evaluated growth, phenological characteristics, yield attributes and productivity potential of above mentioned Bt and non-Bt hybrids. Findings of pot experiment could not establish any significant difference among studied Bt and non-Bt hybrids for root and shoot characteristics. However, arboreum cotton FDK124 exhibited higher root length than all genotypes. Results of second experiment revealed significant differences for phenological characteristics, wherein Bt version of all the tested hybrids exhibited earliness for flowering initiation, 50 % flowering and days to maturity over their non-Bt version. Seed cotton yield (SCY) in all hybrids possessing Bt gene was significantly higher than their non-Bt version primarily due to higher bolls plant-1. Studies conclude that introgression of Bt gene could not alter rooting parameters of evaluated cotton genotypes but significantly affected phenological stages by inducing earliness and higher SCY in Bt versions due to improved bolls plant-1. Keywords: Bt cotton; Near-isogenic lines; Phenology; Root dry matter; Seed cotton yield

Tissues and mechanisms associated with Verticillium wilt resistance in tomato using bi-grafted near-isogenic lines.

Host resistance is the primary means to control Verticillium dahliae, a soil-borne pathogen causing major losses on a broad range of plants, including tomato. The tissues and mechanisms responsible for resistance remain obscure. In the field, resistant tomato used as rootstocks does not confer resistance. Here, we created bi-grafted plants with near-isogenic lines (NILs) exhibiting (Ve1) or lacking (ve1) resistance to V. dahliae race 1. Ten days after inoculation, scion and rootstock tissues were subjected to differential gene expression and co-expression network analyses. Symptoms only developed in susceptible scions regardless of the rootstock. Infection caused more dramatic alteration of tomato gene expression in susceptible compared to resistant tissues, including pathogen receptor, signaling pathway, PR protein, and cell wall modification genes. Differences were observed between scions and rootstocks, primarily related to physiological processes in these tissues. Gene expression in scions was influenced by the rootstock genotype. A few genes were associated with the Ve1 genotype, which was independent of infection or tissue type. Several were physically clustered, some near the Ve1 locus on chromosome 9. Transcripts mapped to V. dahliae were dominated by secreted candidate effector proteins. These findings advance knowledge of molecular mechanisms underlying the tomato-V. dahliae interaction.

Genetic loci regulating the concentrations of anthocyanins and proanthocyanidins in the pericarps of purple and red rice.

The pigmented flavonoids, anthocyanins and proanthocyanidins, have health promoting properties. Previous work determined that the genes Pb and Rc turn on and off the biosynthesis of anthocyanins (purple) and proanthocyanidins (red), respectively. Not yet known is how the concentrations of these pigmented flavonoids are regulated in grain pericarps. Quantitative trait locus (QTL) analysis in a population of rice (Oryza sativa L.) F5 recombinant inbred lines from white pericarp "IR36ae" x red+purple pericarp "242" revealed three QTLs associated with grain concentrations of anthocyanins (TAC) or proanthocyanidins (PA). Both TAC and PA independently mapped to a 1.5 Mb QTL region on chromosome 3 between RM3400 (at 15.8 Mb) and RM15123 (17.3 Mb), named qPR3. Across 2 years, qPR3 explained 36.3% of variance in TAC and 35.8% in PA variance not attributable to Pb or Rc. The qPR3 region encompasses Kala3, a MYB transcription factor previously known to regulate purple grain characteristics. Study of PbPbRcrc progeny showed that TAC of RcRc near isogenic lines (NILs) was 2.1-4.5x that of rcrc. Similarly, study of PbPbRcRc NILs, which had 70% higher PA than pbpbRcRc NILs, revealed a mutual enhancement, not a trade-off between these compounds that share precursors. This suggests that Pb and Rc upregulate genes in a shared pathway as they activate TAC and PA synthesis, respectively. This study provides molecular markers for facilitating marker-assisted selection of qPR3, qPR5, and qPR7 to enhance grain concentrations of pigmented flavonoids and documented that stacking Rc and Pb genes further increases both flavonoid compounds.

A systematic dissection in oilseed rape provides insights into the genetic architecture and molecular mechanism of yield heterosis.

Heterosis refers to the better performance of cross progeny compared with inbred parents, and its utilization contributes greatly to agricultural production. Several hypotheses have been proposed to explain heterosis mainly including dominance, over-dominance (or pseudo-overdominance) and epistasis. However, systematic dissection and verification of these hypotheses are rarely documented. Here, comparison of heterosis level across different traits showed that the strong heterosis of composite traits (such as yield) could be attributed to the multiplicative effects of moderate heterosis of component traits, whether at the genome or locus level. Yield heterosis was regulated by a complex trait-QTL network that was characterized by obvious centre-periphery structure, hub QTL, complex up/downstream and positive/negative feedback relationships. More importantly, we showed that better-parent heterosis on yield could be produced in a cross of two near-isogenic lines by the pyramiding and complementation of two major heterotic QTL showing partial-dominance on yield components. The causal gene (BnaA9.CYP78A9) of QC14 was identified, and its heterotic effect results from the heterozygous status of a CACTA-like transposable element in its upstream regulatory region, which led to partial dominance at expression and auxin levels, thus resulting in non-additive expression of downstream responsive genes involved in cell cycle and proliferation, eventually leading to the heterosis of cell number. Taken together, the results at the phenotypic, genetic and molecular levels were highly consistent, which demonstrated that the pyramiding effect of heterotic QTL and the multiplicative effect of individual component traits could well explain substantial parts of yield heterosis in oilseed rape. These results provide in-depth insights into the genetic architecture and molecular mechanism of yield heterosis.

Multiplexed Host-Induced Gene Silencing of Aspergillus flavus Genes Confers Aflatoxin Resistance in Groundnut.

Aflatoxins are immunosuppressive and carcinogenic secondary metabolites, produced by the filamentous ascomycete Aspergillus flavus, that are hazardous to animal and human health. In this study, we show that multiplexed host-induced gene silencing (HIGS) of Aspergillus flavus genes essential for fungal sporulation and aflatoxin production (nsdC, veA, aflR, and aflM) confers enhanced resistance to Aspergillus infection and aflatoxin contamination in groundnut (<20 ppb). Comparative proteomic analysis of contrasting groundnut genotypes (WT and near-isogenic HIGS lines) supported a better understanding of the molecular processes underlying the induced resistance and identified several groundnut metabolites that might play a significant role in resistance to Aspergillus infection and aflatoxin contamination. Fungal differentiation and pathogenicity proteins, including calmodulin, transcriptional activator-HacA, kynurenine 3-monooxygenase 2, VeA, VelC, and several aflatoxin pathway biosynthetic enzymes, were downregulated in Aspergillus infecting the HIGS lines. Additionally, in the resistant HIGS lines, a number of host resistance proteins associated with fatty acid metabolism were strongly induced, including phosphatidylinositol phosphate kinase, lysophosphatidic acyltransferase-5, palmitoyl-monogalactosyldiacylglycerol Δ-7 desaturase, ceramide kinase-related protein, sphingolipid Δ-8 desaturase, and phospholipase-D. Combined, this knowledge can be used for groundnut pre-breeding and breeding programs to provide a safe and secure food supply.

Transcriptomic profiling of near-isogenic lines reveals candidate genes for a significant locus conferring metribuzin resistance in wheat

BackgroundWeeds reduce wheat yields in dryland farming systems. Herbicides such as metribuzin are commonly used to control weeds. However, wheat has a narrow safety margin against metribuzin. Standing crops such as wheat with weeds in the same field can also be killed by the same dose of metribuzin. Therefore, it is important to identify metribuzin resistance genes and understand the resistance mechanism in wheat for sustainable crop production. A previous study identified a significant metribuzin resistance wheat QTL, Qsns.uwa.4 A.2, explaining 69% of the phenotypic variance for metribuzin resistance.ResultsTwo NIL pairs with the most contrasting performance in the metribuzin treatment and different in genetic backgrounds were compared using RNA sequence analysis, identifying nine candidate genes underlying Qsns.uwa.4 A.2 responsible for metribuzin resistance. Quantitative RT-qPCR further validated the candidate genes, with TraesCS4A03G1099000 (nitrate excretion transporter), TraesCS4A03G1181300 (aspartyl protease), and TraesCS4A03G0741300 (glycine-rich proteins) identified as key factors for metribuzin resistance.ConclusionIdentified markers and key candidate genes can be used for selecting metribuzin resistance in wheat.

Don't panic any more: MORE PANICLES 3 increases rice yield at elevated CO2 levels.

Don't panic any more: MORE PANICLES 3 increases rice yield at elevated CO2 levels.

Long coleoptile genotype and soil texture interactions determine establishment success and early growth parameters of wheat sown at depth

Long coleoptile (LC) - wheat genotypes sown deeper (> 10 cm) than current practice (< 5 cm) enables earlier deep sowing into subsurface soil moisture. The LC genotype shows promise to broaden sowing windows and remove reliance on the shifting seasonal break for crop establishment. The suitability of this sowing strategy to the highly variable soils of the southern Australian wheatbelt warrants investigation to identify suitable environments for potential establishment benefits. We conducted controlled environment room (CER) and field experiments using pairs of near-isogenic-lines (NILs) that differed genetically only for coleoptile length. The experiments investigated the suitability of deep sown (12 cm) LC wheats in soils with different texture (sand versus loam). The CER experiment identified that LC genotypes exhibited greater emergence and early seedling growth, producing greater shoot and root biomass, when deep sown (12 cm) in loam. While LC and short coleoptile (SC) genotypes produced similar shoot growth when sown at depth in sand, LC produced greater root biomass reflected in more seminal roots. In the field on sand, LC and SC had similar emergence and early shoot growth characteristics at all sowing depths (5, 8 and 10 cm). However, the heavy-textured, loam soil reduced seedling emergence and restricted leaf growth with deep sowing (12 cm) in the SC. Our results indicate that the LC trait was more critical for plant emergence and early growth parameters (shoot and root biomass) from depth in a heavier textured loam compared with a sand. By identifying environments that will likely benefit most from deep sown LC, we are one step closer to adapting crop production to shifting rainfall patterns in southern Australia.

Introgression and Environmental Effects on Gene Expression and Aroma Volatiles As Biomarkers in a Melon Near-Isogenic Line

Introgression and Environmental Effects on Gene Expression and Aroma Volatiles As Biomarkers in a Melon Near-Isogenic Line

Mapping and validation of quantitative trait loci associated with dorsal aleurone thickness in rice (Oryza sativa)

Mapping of QTLs for dorsal aleurone thickness (DAT) was performed using chromosome segment substitution lines in rice. Three QTLs, qDAT3.1, qDAT3.2, and qDAT7.1, were detected in multiple environments. As a specified endosperm cell type, the aleurone has an abundance of various nutrients. Increasing the number of aleurone layers is a practicable way of developing highly nutritious cereals. Identifying genes that can increase aleurone thickness is useful for the breeding of aleurone traits to improve the nutritional and health values of rice. Here, we found that iodine staining could efficiently distinguish the aleurone layers, which revealed great variation of the aleurone thickness in rice, especially at the dorsal side of the seed. Therefore, we used a population of chromosome segmental substitution lines (CSSLs) derived from Koshihikari and Nona Bokra for quantitative trait locus (QTL) analysis of the dorsal aleurone thickness (DAT). Three QTLs, qDAT3.1, qDAT3.2, and qDAT7.1, were detected in multiple seasons. Among these, qDAT3.2 colocalizes with Hd6 and Hd16, two QTLs previously identified to regulate the heading date of Koshihikari, explaining the negative correlation between the DAT and days to heading (DTH) in rice. We also provide evidence that early-heading ensures the filling of rice seed under a relatively high temperature to promote aleurone thickening. qDAT7.1, the most stable QTL expressed in different environments, functions independently from heading date. Although Nona Bokra has a lower DAT, its qDAT7.1 allele significantly increased DAT in rice, which was further validated using two near-isogenic lines (NILs). These findings pave the way for further gene cloning of aleurone-related QTLs and may aid the development of highly nutritious rice.

Divergent selection of KNR6 maximizes grain production by balancing the flowering-time adaptation and ear size in maize.

Divergent selection of KNR6 maximizes grain production by balancing the flowering-time adaptation and ear size in maize.

Gene B5 in Cotton Confers High and Broad Resistance to Bacterial Blight and Conditions High Amounts of Sesquiterpenoid Phytoalexins.

Bacterial blight resistance gene B5 has received little attention since it was first described in 1950. A near-isogenic line (NIL) of Gossypium hirsutum cotton, AcB5, was generated in an otherwise bacterial-blight-susceptible 'Acala 44' background. The introgressed locus B5 in AcB5 conferred strong and broad-spectrum resistance to bacterial blight. Segregation patterns of test crosses under Oklahoma field conditions indicated that AcB5 is likely homozygous for resistance at two loci with partial dominance gene action. In controlled-environment conditions, two of the four copies of B5 were required for effective resistance. Contrary to expectations of gene-for-gene theory, AcB5 conferred high resistance toward isogenic strains of Xanthomonas citri subsp. malvacearum carrying cloned avirulence genes avrB4, avrb7, avrBIn, avrB101, and avrB102, respectively, and weaker resistance toward the strain carrying cloned avrb6. The hypothesis that each B gene, in the absence of a polygenic complex, triggers sesquiterpenoid phytoalexin production was tested by measurement of cadalene and lacinilene phytoalexins during resistant responses in five NILs carrying different B genes, four other lines carrying multiple resistance genes, as well as susceptible Ac44E. Phytoalexin production was an obvious, but variable, response in all nine resistant lines. AcB5 accumulated an order of magnitude more of all four phytoalexins than any of the other resistant NILs. Its total levels were comparable to those detected in OK1.2, a highly resistant line that possesses several B genes in a polygenic background.

Transcriptome analysis of two near-isogenic lines of bell pepper (Capsicum annuum) infected with bell pepper endornavirus and pepper mild mottle virus.

Transcriptome analysis of two near-isogenic lines of bell pepper (Capsicum annuum) infected with bell pepper endornavirus and pepper mild mottle virus.