Bulked Segregant RNA Sequencing Research Articles

Broad-Spectrum Resistance to Leaf Rust in the Argentinean Wheat Cultivar “Klein Proteo” Is Controlled by LrKP Located on Chromosome 2BS

Wheat leaf rust, caused by Puccinia triticina, is a severe fungal disease threatening global wheat production. The rational application of genetic loci controlling wheat resistance to leaf rust in breeding practice is still the best choice for disease control. A previous study indicated that the Argentinean wheat cultivar “Klein Proteo” might carry leaf rust resistance (Lr) genes Lr3a and Lr10, as well as an unknown Lr gene. In this study, seedlings of “Klein Proteo” showed high resistance to all the 20 Pt pathotypes isolated in China. Using bulked segregant RNA sequencing (BSR-seq) and developed CAPS markers, the single-dominant gene LrKP was initially mapped to a 114–168 Mb region on chromosome 2BS. Using gene-specific primers of a previously cloned chromosome 2BS-located Lr13 gene, we found that “Klein Proteo” also carried the Lr13 gene. Moreover, the expression of Lr13 in the resistant bulk was significantly higher than that in the susceptible bulk. Nevertheless, “Klein Proteo” showed a much broader and higher resistance compared with the near isogenic line and “ZhouMai 22” carrying Lr13. In conclusion, the wheat cultivar “Klein Proteo” showed great potential in the genetic improvement of wheat resistance to leaf rust in China and the genetic bases controlling the broad-spectrum resistance were initially revealed.

Genetically Dissecting the Novel Powdery Mildew Resistance Gene in Wheat Breeding Line PBDH1607.

Powdery mildew is one of the most destructive diseases in wheat production. Identifying novel resistance genes and deploying them in new cultivars is the most effective approach to minimize wheat losses caused by powdery mildew. In this study, wheat breeding line PBDH1607 showed high resistance to powdery mildew at both the seedling and adult plant stages. Genetic analysis of the seedling data demonstrated that the resistance was controlled by a single dominant gene, tentatively designated PmPBDH. The ΔSNP index based on bulked segregant RNA sequencing indicated that PmPBDH was associated with an interval of about 30.8 Mb (713.5 to 744.3 Mb) on chromosome arm 4AL. Using newly developed markers, we mapped PmPBDH to a 3.2-cM interval covering 7.1 Mb (719,055,516 to 726,215,121 bp). This interval differed from those of Pm61 (717,963,176 to 719,260,469 bp), MlIW30 (732,769,506 to 732,790,522 bp), and MlNSF10 (729,275,816 to 731,365,462 bp) reported on the same chromosome arm. PmPBDH also differed from Pm61, MlIW30, and MlNSF10 by its response spectrum, origin, or inheritance mode, suggesting that PmPBDH should be a new Pm gene. In the candidate interval, five genes were found to be associated with PmPBDH via time course gene expression analysis, and thus they are candidate genes of PmPBDH. Six closely linked markers, including two kompetitive allele-specific PCR markers, were confirmed to be applicable for tracking PmPBDH in marker-assisted breeding.

Molecular Mapping of a Recessive Gene for Stripe Rust Resistance at the YrCf75 Locus Using Bulked Segregant Analysis Combined with Single Nucleotide Polymorphism Genotyping Arrays and Bulked Segregant RNA-Sequencing.

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most important diseases in wheat worldwide. Planting resistant varieties is the most economical, effective, and environment-friendly measure to control wheat stripe rust. Changfeng 75, a Chinese winter wheat variety, shows high stripe rust resistance in both seedling and adult-plant stages. The seedling tests of F1, F2, and F2:3 populations derived from Mingxian 169/Changfeng 75 inoculated with Chinese predominant Pst race CYR34 showed that the stripe rust resistance of Changfeng 75 was controlled by a single recessive gene. The locus was temporarily designated as YrCf75. Bulked segregant analysis (BSA) combined with the wheat 660K single-nucleotide polymorphism (SNP) array and bulked segregant RNA-sequencing indicated that the proportion of polymorphic SNPs on wheat chromosome 2A was the highest, which suggested that YrCf75 was likely located on chromosome 2A. Two hundred and twenty-five Kompetitive allele-specific PCR (KASP) and 75 simple sequence repeat (SSR) markers on chromosome 2A were used to map YrCf75 using the BSA approach. Linkage analysis indicated that 31 KASP markers and one SSR marker were linked to YrCf75, and the genetic distances of the two closest flanking KASP markers, AX-1110060462 and AX-111004763, were 1.2 and 2.7 cM, respectively. YrCf75 was located on wheat chromosome 2AL. The molecular detection, resistance specificity, and chromosome location showed that YrCf75 is likely a new gene that is different from the known stripe rust resistance genes (Yr1 and Yr32) on wheat chromosome 2AL.

Mutations in the miRNA165/166 binding site of the HB2 gene result in pleiotropic effects on morphological traits in wheat

Leaf, spike, stem, and root morphologies are key factors that determine crop growth, development, and productivity. Multiple genes that control these morphological traits have been identified in Arabidopsis, rice, maize, and other plant species. However, little is known about the genomic regions and genes associated with morphological traits in wheat. Here, we identified the ethyl methanesulfonate-derived mutant wheat line M133 that displays multiple morphological changes that include upward-curled leaves, paired spikelets, dwarfism, and delayed heading. Using bulked segregant RNA sequencing (BSR-seq) and a high-resolution genetic map, we identified TraesCS1D02G155200 (HB-D2) as a potential candidate gene. HB-D2 encodes a class III homeodomain-leucine zipper (HD-ZIP III) transcription factor, and the mutation was located in the miRNA165/166 complementary site, resulting in a resistant allele designated rHb-D2. The relative expression of rHb2 in the mutant plants was significantly higher (P < 0.01) than in plants homozygous for the WT allele. Independent resistant mutations that disrupt the miRNA165/166 complementary sites in the A- (rHb-A2) and B-genome (rHb-B2) homoeologs showed similar phenotypic alterations, but the relative intensity of the effects was different. Transgenic plants expressing rHb-D2 gene driven by the maize UBIQUITIN (UBI) promoter showed similar phenotypes to the rHb-D2 mutant. These results confirmed that HB-D2 is the causal gene responsible for the mutant phenotypes. Finally, a survey of 1397 wheat accessions showed that the complementary sites for miRNA165/166 in all three HB2 homoeologs are highly conserved. Our results suggest that HB2 plays an important role in regulating growth and development in wheat.

Fine mapping and analysis of candidate genes for qFT7.1, a major quantitative trait locus controlling flowering time in Brassica rapa L.

qFT7.1, a major QTL for flowering time in Brassica rapa was fine-mapped to chromosome A07 in a 56.4-kb interval, in which the most likely candidate gene is BraA07g018240.3C. In Brassica rapa, flowering time (FT) is an important agronomic trait that affects the yield, quality, and adaption. FT is a complicated trait that is regulated by many genes and is affected greatly by the environment. In this study, a chromosome segment substitution line (CSSL), CSSL16, was selected that showed later flowering than the recurrent parent, a rapid-cycling inbred line of B. rapa (RcBr). Using Bulked Segregant RNA sequencing, we identified a late flowering quantitative trait locus (QTL), designated as qFT7.1, on chromosome A07, based on a secondary-F2 population derived from the cross between CSSL16 and RcBr. qFT7.1 was further validated by conventional QTL mapping. This QTL explained 39.9% (logarithm of odds = 32.2) of the phenotypic variations and was fine mapped to a 56.4-kb interval using recombinant analysis. Expression analysis suggested that BraA07g018240.3C, which is homologous to ATC (encoding Arabidopsis thaliana CENTRORADIALIS homologue), a gene for delayed flowering in Arabidopsis, as the most promising candidate gene. Sequence analysis demonstrated that two synonymous mutations existed in the coding region and numerous bases replacements existed in promoter region between BraA07g018240.3C from CSSL16 and RcBr. The results will increase our knowledge related to the molecular mechanism of late flowering in B. rapa and lays a solid foundation for the breeding of late bolting B. rapa.

Conjunctive Analyses of Bulk Segregant Analysis Sequencing and Bulk Segregant RNA Sequencing to Identify Candidate Genes Controlling Spikelet Sterility of Foxtail Millet.

Foxtail millet has gradually become a model gramineous C4 crop owing to its short growth period and small genome. Research on the development of its spikelets is not only directly related to the yield and economic value of foxtail millet but also can provide a reference for studying the fertility of other C4 crops. In this study, a hybrid population containing 200 offspring was constructed from the Xinong8852 and An15 parental lines, and two extreme trait populations were constructed from the F2 generation for analysis of the spikelet sterility. The F2 population conformed to a 3:1 Mendelian segregation ratio, and it was thus concluded that this trait is likely controlled by a single recessive gene. Bulk segregant analysis sequencing (BSA-Seq) was used to determine the candidate regions and candidate genes related to the development of foxtail millet spikelets. Additionally, the functional analysis of differentially expressed genes in populations with different traits was conducted by bulk segregant RNA sequencing (BSR-Seq). Finally, conjunctive analysis of BSA-Seq and BSR-Seq results, combined with biological information analysis, revealed six genes on chromosome VII that were ultimately identified as candidate genes controlling foxtail millet spikelet development. This study provides a new reference for research on foxtail millet sterility and lays a solid foundation for the examination of fertility in other gramineous crops.

Bulked Segregant RNA Sequencing Revealed Difference Between Virulent and Avirulent Brown Planthoppers.

The brown planthopper (Nilaparvata lugens Stål, BPH) is one of the most devastating insect pests of rice (Oryza sativa L.), but BPH populations have varying degrees of virulence to rice varieties carrying different resistance genes. To help efforts to characterize these variations we applied bulked segregant RNA sequencing (BSR-seq) to identify differentially expressed genes (DEGs) and genetic loci associated with BPH virulence to YHY15 rice plants carrying the resistance gene Bph15. BPHs that are highly virulent or avirulent to these plants were selected from an F2 population to form two contrasting bulks, and BSR-seq identified 751 DEGs between the bulks. Genes associated with carbohydrate, amino acid and nucleotide metabolism, the endocrine system, and signal transduction were upregulated in the avirulent insects when they fed on these plants. The results also indicated that shifts in lipid metabolism and digestive system pathways were crucial for the virulent BPHs’ adaptation to the resistant rice. We identified 24 single-nucleotide polymorphisms (SNPs) in 21 genes linked with BPH virulence. Possible roles of genes apparently linked to BPH virulence are discussed. Our results provide potentially valuable information for further studies of BPH virulence mechanisms and development of robust control strategies.

Map-Based Cloning and Characterization of Br-dyp1, a Gene Conferring Dark Yellow Petal Color Trait in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Flower color is an important trait in Brassica species. However, genes responsible for the dark yellow flower trait in Chinese cabbage have not been reported. In this study, we identified a dark-yellow-flowered Chinese cabbage line SD369. Genetic analysis indicated that the dark yellow flower trait in SD369 was controlled by a single recessive locus, Br-dyp1 (dark yellow petal color 1 in Brassica rapa). Using bulked segregant RNA sequencing and kompetitive allele-specific PCR assays, Br-dyp1 was fine-mapped to an interval of 53.6 kb on chromosome A09. Functional annotation analysis, expression analysis, and sequence variation analysis revealed that Bra037130 (BraA09.ZEP), which encodes a zeaxanthin epoxidase, was the most likely candidate gene for Br-dyp1. Carotenoid profile analysis suggested that Bra037130 (BraA09.ZEP) might participate in the epoxidation from zeaxanthin to violaxanthin. The 679 bp insertion in dark yellow petal caused premature stop codon, thus caused the loss-of-function of the enzyme zeaxanthin epoxidase (ZEP), which disturbed the carotenoid metabolism, and caused the increased accumulation of total carotenoid, and finally converted the flower color from yellow to dark yellow. Comparative transcriptome analysis also showed that the “carotenoid biosynthesis” pathway was significantly enriched, and genes involved in carotenoid degradation and abscisic acid biosynthesis and metabolism were significantly downregulated. Furthermore, we developed and validated the functional marker Br-dyp1-InDel for Br-dyp1. Overall, these results provide insight into the molecular basis of carotenoid-based flower coloration in B. rapa and reveal valuable information for marker-assisted selection breeding in Chinese cabbage.

Identification of a Pm4 Allele as a Powdery Mildew Resistance Gene in Wheat Line Xiaomaomai.

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most destructive foliar diseases of wheat. In this study, we combined the bulked segregant RNA sequencing (BSR-seq) and comparative genomics analysis to localize the powdery mildew resistance gene in Chinese landrace Xiaomaomai. Genetic analysis of F1 plants from a crossing of Xiaomaomai × Lumai23 and the derived F2 population suggests that a single recessive gene, designated as pmXMM, confers the resistance in this germplasm. A genetic linkage map was constructed using the newly developed SNP markers and pmXMM was mapped to the distal end of chromosome 2AL. The two flanking markers 2AL15 and 2AL34 were closely linked to pmXMM at the genetic distance of 3.9 cM and 1.4 cM, respectively. Using the diagnostic primers of Pm4, we confirmed that Xiaomaomai carries a Pm4 allele and the gene function was further validated by the virus-induced gene silencing (VIGS). In addition, we systematically analyzed pmXMM in comparison with the other Pm4 alleles. The results suggest that pmXMM is identical to Pm4d and Pm4e at sequence level. Pm4b is also not different from Pm4c according to their genome/amino acid sequences. Only a few nucleotide variances were detected between pmXMM and Pm4a/b, which indicate the haplotype variation of the Pm4 gene.

Quantitative Trait Loci Mapping and Development of KASP Marker Smut Screening Assay Using High-Density Genetic Map and Bulked Segregant RNA Sequencing in Sugarcane (Saccharum spp.).

Sugarcane is one of the most important industrial crops globally. It is the second largest source of bioethanol, and a major crop for biomass-derived electricity and sugar worldwide. Smut, caused by Sporisorium scitamineum, is a major sugarcane disease in many countries, and is managed by smut-resistant varieties. In China, smut remains the single largest constraint for sugarcane production, and consequently it impacts the value of sugarcane as an energy feedstock. Quantitative trait loci (QTLs) associated with smut resistance and linked diagnostic markers are valuable tools for smut resistance breeding. Here, we developed an F1 population (192 progeny) by crossing two sugarcane varieties with contrasting smut resistance and used for genome-wide single nucleotide polymorphism (SNP) discovery and mapping, using a high-throughput genotyping method called “specific locus amplified fragment sequencing (SLAF-seq) and bulked-segregant RNA sequencing (BSR-seq). SLAF-seq generated 148,500 polymorphic SNP markers. Using SNP and previously identified SSR markers, an integrated genetic map with an average 1.96 cM marker interval was produced. With this genetic map and smut resistance scores of the F1 individuals from four crop years, 21 major QTLs were mapped, with a phenotypic variance explanation (PVE) > 8.0%. Among them, 10 QTLs were stable (repeatable) with PVEs ranging from 8.0 to 81.7%. Further, four QTLs were detected based on BSR-seq analysis. aligning major QTLs with the genome of a sugarcane progenitor Saccharum spontaneum, six markers were found co-localized. Markers located in QTLs and functional annotation of BSR-seq-derived unigenes helped identify four disease resistance candidate genes located in major QTLs. 77 SNPs from major QTLs were then converted to Kompetitive Allele-Specific PCR (KASP) markers, of which five were highly significantly linked to smut resistance. The co-localized QTLs, candidate resistance genes, and KASP markers identified in this study provide practically useful tools for marker-assisted sugarcane smut resistance breeding.

Fine Mapping of qd1, a Dominant Gene that Regulates Stem Elongation in Bread Wheat.

Stem elongation is a critical phase for yield determination and, as a major trait, is targeted for manipulation for improvement in bread wheat (Triticum aestivum L.). In a previous study, we characterized a mutant showing rapid stem elongation but with no effect on plant height at maturity. The present study aimed to finely map the underlying mutated gene, qd1, in this mutant. By analyzing an F2 segregating population consisting of 606 individuals, we found that the qd1 gene behaved in a dominant manner. Moreover, by using the bulked segregant RNA sequencing (BSR-seq)-based linkage analysis method, we initially mapped the qd1 gene to a 13.55 Mb region on chromosome 4B (from 15.41 to 28.96 Mb). This result was further confirmed in F2 and BC3F2 segregating populations. Furthermore, by using transcriptome sequencing data, we developed 14 Kompetitive Allele-Specific PCR (KASP) markers and then mapped the qd1 gene to a smaller and more precise 5.08 Mb interval from 26.80 to 31.88 Mb. To develop additional markers to finely map the qd1 gene, a total of 4,481 single-nucleotide polymorphisms (SNPs) within the 5.08 Mb interval were screened, and 25 KASP markers were developed based on 10x-depth genome resequencing data from both wild-type (WT) and mutant plants. The qd1 gene was finally mapped to a 1.33 Mb interval from 28.86 to 30.19 Mb on chromosome 4B. Four candidate genes were identified in this region. Among them, the expression pattern of only TraesCS4B02G042300 in the stems was concurrent with the stem development of the mutant and WT. The qd1 gene could be used in conjunction with molecular markers to manipulate stem development in the future.

Molecular Characterization of All-Stage and Adult-Plant Resistance Loci Against Powdery Mildew in Winter Wheat Cultivar Liangxing 99 Using BSR-Seq Technology.

Winter wheat cultivar Liangxing 99, which carries gene Pm52, is resistant to powdery mildew at both seedling and adult-plant stages. An F2:6 recombinant inbred line population from cross Liangxing 99 × Zhongzuo 9504 was phenotyped with Blumeria graminis f. sp. tritici isolate Bgt27 at the adult-plant stage in four field tests and the seedling stage in a greenhouse test. The analysis of bulk segregant RNA sequencing (BSR-Seq) identified a single-nucleotide polymorphism-enriched locus, Qaprpm.caas.2B, on chromosome 2BL in the same genomic interval of Pm52 associated with the all-stage resistance (ASR) and Qaprpm.caas.7A on chromosome 7AL associated with the adult-plant resistance (APR) against the disease. Qaprpm.caas.2B was detected in a 1.3 cM genetic interval between markers Xicscl726 and XicsK128 in which Pm52 was placed with a range of logarithm of odds (LOD) values from 28.1 to 34.6, and the phenotype variations explained in terms of maximum disease severity (MDS) ranged from 45 to 52%. The LOD peak of Qaprpm.caas.7A was localized in a 4.6 cM interval between markers XicsK7A8 and XicsK7A26 and explained the phenotypic variation of MDS ranging from 13 to 16%. The results of this study confirmed Pm52 for ASR and identified Qaprpm.caas.7A for APR to powdery mildew in Liangxing 99.

TEAseq-based identification of 35,696 Dissociation insertional mutations facilitates functional genomic studies in maize

In plants, transposable element (TE)-triggered mutants are important resources for functional genomic studies. However, conventional approaches for genome-wide identification of TE insertion sites are costly and laborious. This study developed a novel, rapid, and high-throughput TE insertion site identification workflow based on next-generation sequencing and named it Transposable Element Amplicon Sequencing (TEAseq). Using TEAseq, we systemically profiled the Dissociation (Ds) insertion sites in 1606 independent Ds insertional mutants in advanced backcross generation using K17 as background. The Ac-containing individuals were excluded for getting rid of the potential somatic insertions. We characterized 35,696 germinal Ds insertions tagging 10,323 genes, representing approximately 23.3% of the total genes in the maize genome. The insertion sites were presented in chromosomal hotspots around the ancestral Ds loci, and insertions occurred preferentially in gene body regions. Furthermore, we mapped a loss-of-function AGL2 gene using bulked segregant RNA-sequencing assay and proved that AGL2 is essential for seed development. We additionally established an open-access database named MEILAM for easy access to Ds insertional mutations. Overall, our results have provided an efficient workflow for TE insertion identification and rich sequence-indexed mutant resources for maize functional genomic studies.

RNA-seq bulked segregant analysis combined with KASP genotyping rapidly identified PmCH7087 as responsible for powdery mildew resistance in wheat.

Powdery mildew causes considerable yield losses in common wheat (Triticum aestivum L.) production. Mapping and cloning powdery mildew-resistant quantitative trait loci can benefit stable yield production by facilitating the breeding of resistant varieties. In this study, we used the powdery mildew resistance introgression line 'CH7087' (harboring the resistance gene PmCH7087) and developed a large F2 population and a corresponding F2:3 segregation population containing 2,000 family lines for molecular mapping of PmCH7087. Genetic analysis demonstrated that the resistance phenotype was controlled by a single dominant gene. According to the performance exhibited by the F2:3 lines, 50 resistant lines and 50 susceptible lines without phenotype segregation were chosen for pooling and bulked segregant RNA sequencing (BSR-Seq) analysis. A region spanning 42.77 Mb was identified, and genotyping of an additional 183 F2:3 lines with extreme phenotypes using 20 kompetitive allele-specific polymerase chain reaction (KASP) markers in the BSR-Seq mapping regions confirmed this region and narrowed it to 9.68 Mb, in which 45 genes were identified and annotated. Five of these transcripts harbored nonsynonymous single nucleotide polymorphisms between the two parents, with the transcripts of TraesCS2B01G302800 being involved in signal transduction. Furthermore, TraesCS2B01G302800.2 was annotated as the closest homologue of serine/threonine-protein kinase PBS1, a typical participant in the plant disease immune response, indicating that TraesCS2B01G302800 was the candidate gene of PmCH7087. Our results may facilitate future research attempting to improve powdery mildew resistance in wheat and to identify candidate genes for further verification and gene cloning.

Mapping of wheat stripe rust resistance gene Yr041133 by BSR-Seq analysis

Puccinia striiformis Westend. f. sp. tritici (Pst) pathotype CYR34 is widely virulent and prevalent in China. Here, we report identification of a strpie rust resistance (Yr) gene, designated Yr041133, in winter wheat line 041133. This line produced a hypersensitive reaction to CYR34 and conferred resistance to 13 other pathotypes. Resistance to CYR34 in line 041133 was controlled by a single dominant gene. Bulked segregant RNA sequencing (BSR-Seq) was performed on a pair of RNA bulks generated by pooling resistant and susceptible recombinant inbred lines. Yr041133 was mapped to a 1.7 cM genetic interval on the chromosome arm 7BL that corresponded to a 0.8 Mb physical interval (608.9–609.7 Mb) in the Chinese Spring reference genome. Based on its unique physical location Yr041133 differred from the other Yr genes on this chromosome arm.

A predicted NEDD8 conjugating enzyme gene identified as a Capsicum candidate Rf gene using bulk segregant RNA sequencing

Cytoplasmic male sterility (CMS) is an important tool for producing F1 hybrids, which can exhibit heterosis. The companion system, restorer-of-fertility (Rf), is poorly understood at the molecular level and would be valuable in producing restorer lines for hybrid seed production. The identity of the Rf gene in Capsicum (pepper) is currently unclear. In this study, using bulked segregant RNA sequencing (BSR-seq), a strong candidate Rf gene, Capana06g002866, which is annotated as a NEDD8 conjugating enzyme E2, was identified. Capana06g002866 has an ORF of 555 bp in length encoding 184 amino acids; it can be cloned from F1 plants from the hybridization of the CMS line 8A and restorer line R1 but is not found in CMS line 8A. With qRT-PCR validation, Capana06g002866 was found to be upregulated in restorer accessions compared to sterile accessions. The relative expression in flower buds increased with the developmental stage in F1 plants, while the expression was very low in all flower bud stages of the CMS lines. These results provide new insights into the Rf gene in pepper and will be useful for other crops utilizing the CMS system.

Identification of a locus associated with genic male sterility in maize via EMS mutagenesis and bulked-segregant RNA-seq

Genic male sterility (GMS) is one of the most important resources for exploiting heterosis in crop breeding, so that identifying genomic loci regulating GMS is desirable. However, many regulatory genes controlling GMS have not yet been characterized in maize, owing partly to a lack of genetic materials. We generated a recessive male-sterile maize mutant in the Jing 724 genetic background via ethyl methanesulfonate treatment, and found the male sterility to be due to a single gene mutation. Bulk-segregant RNA sequencing of three replicates indicated that one genomic region located at the end of chromosome 4 was associated with the observed mutant phenotype. Among genes with nonsynonymous mutations, Zm00001d053895 (bHLH51) showed abolished expression in the sterile bulks and was annotated as a bHLH transcription factor orthologous to Arabidopsis AMS, suggesting an association with the male sterility of the mutant. Kompetitive Allele-Specific PCR assays further validated the exclusive correlation of male sterility with the single C-to-T mutation in the fifth exon. The new maize mutant and the potential SNP locus provide novel genetic material for investigating the molecular mechanism underlying tapetal development and may facilitate the improvement of hybrid production systems.

A rare single nucleotide variant in Pm5e confers powdery mildew resistance in common wheat.

Powdery mildew poses severe threats to wheat production. The most sustainable way to control this disease is through planting resistant cultivars. We report the map-based cloning of the powdery mildew resistance allele Pm5e from a Chinese wheat landrace. We applied a two-step bulked segregant RNA sequencing (BSR-Seq) approach in developing tightly linked or co-segregating markers to Pm5e. The first BSR-Seq used phenotypically contrasting bulks of recombinant inbred lines (RILs) to identify Pm5e-linked markers. The second BSR-Seq utilized bulks of genetic recombinants screened from a fine-mapping population to precisely quantify the associated genomic variation in the mapping interval, and identified the Pm5e candidate genes. The function of Pm5e was validated by transgenic assay, loss-of-function mutants and haplotype association analysis. Pm5e encodes a nucleotide-binding domain leucine-rich-repeat-containing (NLR) protein. A rare nonsynonymous single nucleotide variant (SNV) within the C-terminal leucine rich repeat (LRR) domain is responsible for the gain of powdery mildew resistance function of Pm5e, an allele endemic to wheat landraces of Shaanxi province of China. Results from this study demonstrate the value of landraces in discovering useful genes for modern wheat breeding. The key SNV associated with powdery mildew resistance will be useful for marker-assisted selection of Pm5e in wheat breeding programs.

Two Clubroot-Resistance Genes, Rcr3 and Rcr9wa, Mapped in Brassica rapa Using Bulk Segregant RNA Sequencing.

Genetic resistance is widely used to manage clubroot (Plasmodiophora brassicae) in brassica crops, but new pathotypes have recently been identified on canola (Brassica napus) on the Canadian prairies. Resistance effective against both the most prevalent pathotype (3H, based on the Canadian Clubroot Differential system) and the new pathotypes is needed. BC1 plants of Brassica rapa from a cross of line 96-6990-2 (clubroot resistance originating from turnip cultivar ‘Waaslander’) and a susceptible doubled-haploid line, ACDC, exhibited a 1:1 segregation for resistance against pathotypes 3H and 5X. A resistance gene designated as Rcr3 was mapped initially based on the percentage of polymorphic variants using bulked segregant RNA sequencing (BSR-Seq) and further mapped using Kompetitive Allele Specific PCR. DNA variants were identified by assembling short reads against a reference genome of B. rapa. Rcr3 was mapped into chromosome A08. It was flanked by single nucleotide polymorphisms (SNP) markers (A90_A08_SNP_M12 and M16) between 10.00 and 10.23 Mb, in an interval of 231.6 Kb. There were 32 genes in the Rcr3 interval. Three genes (Bra020951, Bra020974, and Bra020979) were annotated with disease resistance mechanisms, which are potential candidates for Rcr3. Another resistance gene, designated as Rcr9wa, for resistance to pathotype 5X was mapped, with the flanking markers (A90_A08_SNP_M28 and M79) between 10.85 and 11.17 Mb using the SNP sites identified through BSR-Seq for Rcr3. There were 44 genes in the Rcr9wa interval, three of which (Bra020827, Bra020828, Bra020814) were annotated as immune-system-process related genes, which are potential candidates for Rcr9wa.

Identification of two recessive etiolation genes (py1, py2) in pakchoi (Brassica rapa L. ssp. chinensis)

BackgroundLeaf color is a major agronomic trait, which has a strong influence on crop yields. Isolating leaf color mutants can represent valuable materials for research in chlorophyll (Chl) biosynthesis and metabolism regulation.ResultsIn this study, we identified a stably inherited yellow leaf mutant derived from ‘Huaguan’ pakchoi variety via isolated microspore culture and designated as pylm. This mutant displayed yellow leaves after germination. Its etiolated phenotype was nonlethal and stable during the whole growth period. Its growth was weak and its hypocotyls were markedly elongated. Genetic analysis revealed that two recessive nuclear genes, named py1 and py2, are responsible for the etiolation phenotype. Bulked segregant RNA sequencing (BSR-Seq) showed that py1 and py2 were mapped on chromosomes A09 and A07, respectively. The genes were single Mendelian factors in F3:4 populations based on a 3:1 phenotypic segregation ratio. The py1 was localized to a 258.3-kb interval on a 34-gene genome. The differentially expressed gene BraA09004189 was detected in the py1 mapping region and regulated heme catabolism. One single-nucleotide polymorphism (SNP) of BraA09004189 occurred in pylm. A candidate gene-specific SNP marker in 1520 F3:4 yellow-colored individuals co-segregated with py1. For py2, 1860 recessive homozygous F3:4 individuals were investigated and localized py2 to a 4.4-kb interval. Of the five genes in this region, BraA07001774 was predicted as a candidate for py2. It encoded an embryo defective 1187 and a phosphotransferase related to chlorophyll deficiency and hypocotyl elongation. One SNP of BraA07001774 occurred in pylm. It caused a single amino acid mutation from Asp to Asn. According to quantitative real-time polymerase chain reaction (qRT-PCR), BraA07001774 was downregulated in pylm.ConclusionsOur study identified a Chl deficiency mutant pylm in pakchoi. Two recessive nuclear genes named py1 and py2 had a significant effect on etiolation. Candidate genes regulating etiolation were identified as BraA09004189 and BraA07001774, respectively. These findings will elucidate chlorophyll metabolism and the molecular mechanisms of the gene interactions controlling pakchoi etiolation.