Bulked Segregant RNA-Seq Reveals Distinct Expression Profiling in Chinese Wheat Cultivar Jimai 23 Responding to Powdery Mildew

Analysis of Defense-Related Gene Expression and Leaf Metabolome in Wheat During the Early Infection Stages of Blumeria graminis f. sp. tritici

Blumeria graminis f. sp. tritici (Bgt) is an obligate biotrophic fungal pathogen responsible for powdery mildew in bread wheat (Triticum aestivum). Upon Bgt infection, the wheat plant activates basal defense mechanisms, namely PAMP-triggered immunity, in the leaves during the first few days. Understanding this early stage of quantitative resistance is crucial for developing new breeding tools and evaluating plant resistance inducers for sustainable agricultural practices. In this sense, we used a combination of transcriptomic and metabolomic approaches to analyze the early steps of the interaction between Bgt and the moderately susceptible wheat cultivar Pakito. Bgt infection resulted in an increasing expression of genes encoding pathogenesis-related (PR) proteins (PR1, PR4, PR5, and PR8) known to target the pathogen, during the first 48 h postinoculation. Moreover, RT-qPCR and metabolomic analyses pointed out the importance of the phenylpropanoid pathway in quantitative resistance against Bgt. Among metabolites linked to this pathway, hydroxycinnamic acid amides containing agmatine and putrescine as amine components accumulated from the second to the fourth day after inoculation. This suggests their involvement in quantitative resistance via cross-linking processes in cell walls for reinforcement, which is supported by the up-regulation of PAL (phenylalanine ammonia-lyase), PR15 (oxalate oxidase) and POX (peroxidase) after inoculation. Finally, pipecolic acid, which is considered a signal involved in systemic acquired resistance, accumulated after inoculation. These new insights lead to a better understanding of basal defense in wheat leaves after Bgt infection.

Identification of a Wheat Powdery Mildew Dominant Resistance Gene in the Pm5 Locus for High-Throughput Marker-Assisted Selection

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), poses a severe threat to wheat yield and quality worldwide. Rapid identification and the accurate transference of effective resistance genes are important to the development of resistant cultivars and the sustainable control of this disease. In the present study, the wheat line AL11 exhibited high levels of resistance to powdery mildew at both the seedling and adult plant stages. Genetic analysis of the AL11 × 'Shixin 733' mapping population revealed that its resistance was controlled by a single dominant gene, tentatively designated PmAL11. Using bulked segregant RNA-Seq and molecular marker analysis, PmAL11 was mapped to the Pm5 locus on chromosome 7B where it cosegregated with the functional marker Pm5e-KASP. Sequence alignment analysis revealed that the Pm5e-homologous sequence in AL11 was identical to the reported recessive gene Pm5e in wheat landrace 'Fuzhuang 30'. It appears that PmAL11 was most probably Pm5e, but it was mediated by a dominant inheritance pattern, so it should provide a valuable resistance resource for both genetic study and wheat breeding. To efficiently use and trace PmAL11 in breeding, a new kompetitive allele-specific PCR marker AL11-K2488 that cosegregated with this gene was developed and confirmed to be applicable in the different wheat backgrounds, thus promoting its use in the marker-assisted selection of PmAL11.

Screening of Candidate Genes Associated with Brown Stripe Resistance in Sugarcane via BSR-seq Analysis

Sugarcane brown stripe (SBS), caused by the fungal pathogen Helminthosporium stenospilum, is one of the most serious threats to sugarcane production. However, its outbreaks and epidemics require suitable climatic conditions, resulting in the inefficient improvement of the SBS resistance by phenotype selection. The sugarcane F₁ population of SBS-resistant YT93-159 × SBS-susceptible ROC22 was used for constructing the bulks. Bulked segregant RNA-seq (BSR-seq) was then performed on the parents YT93-159 (T01) and ROC22 (T02), and the opposite bulks of 30 SBS-susceptible individuals mixed bulk (T03) and 30 SBS-resistant individuals mixed bulk (T04) collected from 287 F₁ individuals. A total of 170.00 Gb of clean data containing 297,921 SNPs and 70,426 genes were obtained. Differentially expressed genes (DEGs) analysis suggested that 7787 and 5911 DEGs were identified in the parents (T01 vs. T02) and two mixed bulks (T03 vs. T04), respectively. In addition, 25,363 high-quality and credible SNPs were obtained using the genome analysis toolkit GATK for SNP calling. Subsequently, six candidate regions with a total length of 8.72 Mb, which were located in the chromosomes 4B and 7C of sugarcane wild species Saccharum spontaneum, were identified, and 279 genes associated with SBS-resistance were annotated by ED algorithm and ΔSNP-index. Furthermore, the expression profiles of candidate genes were verified by quantitative real-time PCR (qRT-PCR) analysis, and the results showed that eight genes (LRR-RLK, DHAR1, WRKY7, RLK1, BLH4, AK3, CRK34, and NDA2) and seven genes (WRKY31, CIPK2, CKA1, CDPK6, PFK4, CBL2, and PR2) of the 20 tested genes were significantly up-regulated in YT93-159 and ROC22, respectively. Finally, a potential molecular mechanism of sugarcane response to H. stenospilum infection is illustrate that the activations of ROS signaling, MAPK cascade signaling, Ca²⁺ signaling, ABA signaling, and the ASA-GSH cycle jointly promote the SBS resistance in sugarcane. This study provides abundant gene resources for the SBS resistance breeding in sugarcane.

Genetic identification of SNP markers and candidate genes associated with sugarcane smut resistance using BSR-Seq

Sugarcane smut caused by Sporisorium scitamineum is one of the most severe fungal diseases worldwide. In this study, a cross was made between a smut-resistant variety YT93-159 and a smut-susceptible variety ROC22, and 312 progenies were obtained. Two bulks of progenies were then constructed, one consisted of 27 highly smut resistant progenies and the other 24 smut susceptible progenies. Total RNAs of the progenies of each bulk, were pooled and subject to bulked segregant RNA-sequence analysis (BSR-Seq). A total of 164.44 Gb clean data containing 2,341,449 SNPs and 64,999 genes were obtained, 7,295 of which were differentially expressed genes (DEGs). These DEGs were mainly enriched in stress-related metabolic pathways, including carbon metabolism, phenylalanine metabolism, plant hormone signal transduction, glutathione metabolism, and plant-pathogen interactions. Besides, 45,946 high-quality, credible SNPs, a 1.27 Mb region at Saccharum spontaneum chromosome Chr5B (68,904,827 to 70,172,982), and 129 candidate genes were identified to be associated with smut resistance. Among them, twenty-four genes, either encoding key enzymes involved in signaling pathways or being transcription factors, were found to be very closely associated with stress resistance. RT-qPCR analysis demonstrated that they played a positive role in smut resistance. Finally, a potential molecular mechanism of sugarcane and S. scitamineum interaction is depicted that activations of MAPK cascade signaling, ROS signaling, Ca²⁺ signaling, and PAL metabolic pathway and initiation of the glyoxalase system jointly promote the resistance to S. scitamineum in sugarcane. This study provides potential SNP markers and candidate gene resources for smut resistance breeding in sugarcane.

Identification of the powdery mildew resistance gene in wheat breeding line Yannong 99102-06188 via bulked segregant exome capture sequencing

Powdery mildew of wheat (Triticum aestivum), caused by Blumeria graminis f.sp. tritici (Bgt), is a destructive disease that seriously threatens the yield and quality of its host. Identifying resistance genes is the most attractive and effective strategy for developing disease-resistant cultivars and controlling this disease. In this study, a wheat breeding line Yannong 99102-06188 (YN99102), an elite derivative line from the same breeding process as the famous wheat cultivar Yannong 999, showed high resistance to powdery mildew at the whole growth stages. Genetic analysis was carried out using Bgt isolate E09 and a population of YN99102 crossed with a susceptible parent Jinhe 13-205 (JH13-205). The result indicated that a single recessive gene, tentatively designated pmYN99102, conferred seedling resistance to the Bgt isolate E09. Using bulked segregant exome capture sequencing (BSE-Seq), pmYN99102 was physically located to a ~33.7 Mb (691.0-724.7 Mb) interval on the chromosome arm 2BL, and this interval was further locked in a 1.5 cM genetic interval using molecular markers, which was aligned to a 9.0 Mb physical interval (699.2-708.2 Mb). Based on the analysis of physical location, origin, resistant spectrum, and inherited pattern, pmYN99102 differed from those of the reported powdery mildew (Pm) resistance genes on 2BL, suggesting pmYN99102 is most likely a new Pm gene/allele in the targeted interval. To transfer pmYN99102 to different genetic backgrounds using marker-assisted selection (MAS), 18 closely linked markers were tested for their availability in different genetic backgrounds for MAS, and all markers expect for YTU103-97 can be used in MAS for tracking pmYN99102 when it transferred into those susceptible cultivars.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

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.

Identification of Rf Genes in Hexaploid Wheat (Triticumaestivum L.) by RNA-Seq and Paralog Analyses

Among the natural mechanisms used for wheat hybrid breeding, the most desirable is the system combining the cytoplasmic male sterility (cms) of the female parent with the fertility-restoring genes (Rf) of the male parent. The objective of this study was to identify Rf candidate genes in the wheat genome on the basis of transcriptome sequencing (RNA-seq) and paralog analysis data. Total RNA was isolated from the anthers of two fertility-restorer (Primépi and Patras) and two non-restorer (Astoria and Grana) varieties at the tetrad and late uninucleate microspore stages. Of 36,912 differentially expressed genes (DEGs), 21 encoding domains in known fertility-restoring proteins were selected. To enrich the pool of Rf candidates, 52 paralogs (PAGs) of the 21 selected DEGs were included in the analyses. The expression profiles of most of the DEGs and PAGs determined bioinformatically were as expected (i.e., they were overexpressed in at least one fertility-restorer variety). However, these results were only partially consistent with the quantitative real-time PCR data. The DEG and PAG promoters included cis-regulatory elements common among PPR-encoding genes. On the basis of the obtained results, we designated seven genes as Rf candidate genes, six of which were identified for the first time in this study.

Bulked Segregant RNA-Seq Provides Distinctive Expression Profile Against Powdery Mildew in the Wheat Genotype YD588

Wheat powdery mildew, caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), is a destructive disease leading to huge yield losses in production. Host resistance can greatly contribute to the control of the disease. To explore potential genes related to the powdery mildew (Pm) resistance, in this study, we used a resistant genotype YD588 to investigate the potential resistance components and profiled its expression in response to powdery mildew infection. Genetic analysis showed that a single dominant gene, tentatively designated PmYD588, conferred resistance to powdery mildew in YD588. Using bulked segregant RNA-Seq (BSR-Seq) and single nucleotide polymorphism (SNP) association analysis, two high-confidence candidate regions were detected in the chromosome arm 2B, spanning 453,752,054-506,356,791 and 584,117,809-664,221,850 bp, respectively. To confirm the candidate region, molecular markers were developed using the BSR-Seq data and mapped PmYD588 to an interval of 4.2 cM by using the markers YTU588-004 and YTU588-008. The physical position was subsequently locked into the interval of 647.1-656.0 Mb, which was different from those of Pm6, Pm33, Pm51, Pm52, Pm63, Pm64, PmQ, PmKN0816, MlZec1, and MlAB10 on the same chromosome arm in its position, suggesting that it is most likely a new Pm gene. To explore the potential regulatory genes of the R gene, 2,973 differentially expressed genes (DEGs) between the parents and bulks were analyzed using gene ontology (GO), clusters of orthologous group (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Based on the data, we selected 23 potential regulated genes in the enriched pathway of plant-pathogen interaction and detected their temporal expression patterns using an additional set of wheat samples and time-course analysis postinoculation with Bgt. As a result, six disease-related genes showed distinctive expression profiles after Bgt invasion and can serve as key candidates for the dissection of resistance mechanisms and improvement of durable resistance to wheat powdery mildew.