Prediction and analysis of three gene families related to leaf rust (Puccinia triticina) resistance in wheat (Triticum aestivum L.)

Leaf rust (Puccinia recondita f. sp. secalis) triggers substantial changes in rye (Secale cereale L.) at the transcriptome and metabolome levels

BACKGROUND

Rye (Secale cereale L.) is a cereal crop highly tolerant to environmental stresses, including abiotic and biotic stresses (e.g., fungal diseases). Among these fungal diseases, leaf rust (LR) is a major threat to rye production. Despite extensive research, the genetic basis of the rye immune response to LR remains unclear.

RESULTS

An RNA-seq analysis was conducted to examine the immune response of three unrelated rye inbred lines (D33, D39, and L318) infected with compatible and incompatible Puccinia recondita f. sp. secalis (Prs) isolates. In total, 877 unique differentially expressed genes (DEGs) were identified at 20 and 36 h post-treatment (hpt). Most of the DEGs were up-regulated. Two lines (D39 and L318) had more up-regulated genes than down-regulated genes, whereas the opposite trend was observed for line D33. The functional classification of the DEGs helped identify the largest gene groups regulated by LR. Notably, these groups included several DEGs encoding cytochrome P450, receptor-like kinases, methylesterases, pathogenesis-related protein-1, xyloglucan endotransglucosylases/hydrolases, and peroxidases. The metabolomic response was highly conserved among the genotypes, with line D33 displaying the most genotype-specific changes in secondary metabolites. The effect of pathogen compatibility on metabolomic changes was less than the effects of the time-points and genotypes. Accordingly, the secondary metabolome of rye is altered by the recognition of the pathogen rather than by a successful infection. The results of the enrichment analysis of the DEGs and differentially accumulated metabolites (DAMs) reflected the involvement of phenylpropanoid and diterpenoid biosynthesis as well as thiamine metabolism in the rye immune response.

CONCLUSION

Our work provides novel insights into the genetic and metabolic responses of rye to LR. Numerous immune response-related DEGs and DAMs were identified, thereby clarifying the mechanisms underlying the rye response to compatible and incompatible Prs isolates during the early stages of LR development. The integration of transcriptomic and metabolomic analyses elucidated the contributions of phenylpropanoid biosynthesis and flavonoid pathways to the rye immune response to Prs. This combined analysis of omics data provides valuable insights relevant for future research conducted to enhance rye resistance to LR.

Genome-wide analysis of transmembrane 9 superfamily genes in wheat (Triticum aestivum) and their expression in the roots under nitrogen limitation and Bacillus amyloliquefaciens PDR1 treatment conditions

Introduction

Transmembrane 9 superfamily (TM9SF) proteins play significant roles in plant physiology. However, these proteins are poorly characterized in wheat (Triticum aestivum). The present study aimed at the genome-wide analysis of putative wheat TM9SF (TraesTM9SF) proteins and their potential involvement in response to nitrogen limitation and Bacillus amyloliquefaciens PDR1 treatments.

Methods

TraesTM9SF genes were retrieved from the wheat genome, and their physiochemical properties, alignment, phylogenetic, motif structure, cis-regulatory element, synteny, protein-protein interaction (PPI), and transcription factor (TF) prediction analyses were performed. Transcriptome sequencing and quantitative real-time polymerase reaction (qRT-PCR) were performed to detect gene expression in roots under single or combined treatments with nitrogen limitation and B. amyloliquefaciens PDR1.

Results and discussion

Forty-seven TraesTM9SF genes were identified in the wheat genome, highlighting the significance of these genes in wheat. TraesTM9SF genes were absent on some wheat chromosomes and were unevenly distributed on the other chromosomes, indicating that potential regulatory functions and evolutionary events may have shaped the TraesTM9SF gene family. Fifty-four cis-regulatory elements, including light-response, hormone response, biotic/abiotic stress, and development cis-regulatory elements, were present in the TraesTM9SF promoter regions. No duplication of TraesTM9SF genes in the wheat genome was recorded, and 177 TFs were predicted to target the 47 TraesTM9SF genes in a complex regulatory network. These findings offer valued data for predicting the putative functions of uncharacterized TM9SF genes. Moreover, transcriptome analysis and validation by qRT-PCR indicated that the TraesTM9SF genes are expressed in the root system of wheat and are potentially involved in the response of this plant to single or combined treatments with nitrogen limitation and B. amyloliquefaciens PDR1, suggesting their functional roles in plant growth, development, and stress responses.

Conclusion

These findings may be vital in further investigation of the function and biological applications of TM9SF genes in wheat.

Genome-Wide Association Study Reveals Novel Powdery Mildew Resistance Loci in Bread Wheat

Powdery mildew (PM), caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), significantly threatens global bread wheat production. Although the use of resistant cultivars is an effective strategy for managing PM, currently available wheat cultivars lack sufficient levels of resistance. To tackle this challenge, we conducted a comprehensive genome-wide association study (GWAS) using a diverse panel of 286 bread wheat genotypes. Over three consecutive years (2020-2021, 2021-2022, and 2022-2023), these genotypes were extensively evaluated for PM severity under field conditions following inoculation with virulent Bgt isolates. The panel was previously genotyped using the Illumina 90K Infinium iSelect assay to obtain genome-wide single-nucleotide polymorphism (SNP) marker coverage. By applying FarmCPU, a multilocus mixed model, we identified a total of 113 marker-trait associations (MTAs) located on chromosomes 1A, 1B, 2B, 3A, 3B, 4A, 4B, 5A, 5B, 6B, 7A, and 7B at a significance level of p ≤ 0.001. Notably, four novel MTAs on chromosome 6B were consistently detected in 2020-2021 and 2021-2022. Furthermore, within the confidence intervals of the identified SNPs, we identified 96 candidate genes belonging to different proteins including 12 disease resistance/host-pathogen interaction-related protein families. Among these, protein kinases, leucine-rich repeats, and zinc finger proteins were of particular interest due to their potential roles in PM resistance. These identified loci can serve as targets for breeding programs aimed at developing disease-resistant wheat cultivars.

Identification and characteristics of wheat Lr orthologs in three rye inbred lines

The genetic background of the immune response of rye to leaf rust (LR), although extensively studied, is still not well understood. The recent publication of the genome of rye line Lo7 and the development of efficient transcriptomic methods has aided the search for genes that confer resistance to this disease. In this study, we investigated the potential role of rye orthologs of wheat Lr genes (Lr1, Lr10, Lr21, Lr22a, and RGA2/T10rga2-1A) in the LR seedling-stage resistance of inbred rye lines D33, D39, and L318. Bioinformatics analysis uncovered numerous Lr orthologs in the Lo7 genome, namely, 14 ScLr1, 15 ScRga2, and 2 ScLr21 paralogs, and 1 each of ScLr10 and ScLr22a genes. The paralogs of ScLr1, ScRga2, and ScLr21 were structurally different from one another and their wheat counterparts. According to an RNA sequencing analysis, only four wheat Lr gene orthologs identified in the Lo7 genome (ScLr1_3, ScLr1_4, ScLr1_8, and ScRga2_6) were differentially expressed; all four were downregulated after infection with compatible or incompatible isolates of Puccinia recondita f. sp. secalis (Prs). Using a more precise tool, RT-qPCR, we found that two genes were upregulated at 20 h post-infection, namely, ScLr1_4 and ScLr1_8 in lines D33 and D39, respectively, both of which have been found to be resistant to LR under field conditions and after treatment with a semi-compatible Prs strain. We were unable to discern any universal pattern of gene expression after Prs infection; on the contrary, all detected relationships were plant genotype-, Prs isolate-, or time-specific. Nevertheless, at least some Lr orthologs in rye (namely, ScLr1_3 ScLr1_4, ScLr1_8, and ScRga2_6), even though mainly downregulated, may play an important role in the response of rye to LR.

Exploring genomic regions involved in bread wheat resistance to leaf rust at seedling/adult stages by using GWAS analysis

BACKGROUND

Global wheat productivity is seriously challenged by a range of rust pathogens, especially leaf rust derived from Puccinia triticina. Since the most efficient approach to control leaf rust is genetic resistance, many efforts have been made to uncover resistance genes; however, it demands an ongoing exploration for effective resistance sources because of the advent of novel virulent races. Thus, the current study was focused on detecting leaf rust resistance-related genomic loci against the P. triticina prevalent races by GWAS in a set of Iranian cultivars and landraces.

RESULTS

Evaluation of 320 Iranian bread wheat cultivars and landraces against four prevalent rust pathotypes of P. triticina (LR-99-2, LR-98-12, LR-98-22, and LR-97-12) indicated the diversity in wheat accessions responses to P. triticina. From GWAS results, 80 leaf rust resistance QTLs were located in the surrounding known QTLs/genes on almost chromosomes, except for 1D, 3D, 4D, and 7D. Of these, six MTAs (rs20781/rs20782 associated with resistance to LR-97-12; rs49543/rs52026 for LR-98-22; rs44885/rs44886 for LR-98-22/LR-98-1/LR-99-2) were found on genomic regions where no resistance genes previously reported, suggesting new loci conferring resistance to leaf rust. The GBLUP genomic prediction model appeared better than RR-BLUP and BRR, reflecting that GBLUP is a potent model for genomic selection in wheat accessions.

CONCLUSIONS

Overall, the newly identified MTAs as well as the highly resistant accessions in the recent work provide an opportunity towards improving leaf rust resistance.

Harnessing genetic resistance to rusts in wheat and integrated rust management methods to develop more durable resistant cultivars

Wheat is one of the most important staple foods on earth. Leaf rust, stem rust and stripe rust, caused by Puccini triticina, Puccinia f. sp. graminis and Puccinia f. sp. striiformis, respectively, continue to threaten wheat production worldwide. Utilization of resistant cultivars is the most effective and chemical-free strategy to control rust diseases. Convectional and molecular biology techniques identified more than 200 resistance genes and their associated markers from common wheat and wheat wild relatives, which can be used by breeders in resistance breeding programmes. However, there is continuous emergence of new races of rust pathogens with novel degrees of virulence, thus rendering wheat resistance genes ineffective. An integration of genomic selection, genome editing, molecular breeding and marker-assisted selection, and phenotypic evaluations is required in developing high quality wheat varieties with resistance to multiple pathogens. Although host genotype resistance and application of fungicides are the most generally utilized approaches for controlling wheat rusts, effective agronomic methods are required to reduce disease management costs and increase wheat production sustainability. This review gives a critical overview of the current knowledge of rust resistance, particularly race-specific and non-race specific resistance, the role of pathogenesis-related proteins, non-coding RNAs, and transcription factors in rust resistance, and the molecular basis of interactions between wheat and rust pathogens. It will also discuss the new advances on how integrated rust management methods can assist in developing more durable resistant cultivars in these pathosystems.

Diagnostic accuracy of genetic markers for identification of the Lr46/Yr29 “slow rusting” locus in wheat (Triticum aestivum L.)

Wheat leaf rust, caused by fungal pathogen Puccinia triticina Erikss, annually contributes to production losses as high as 40% in susceptible varieties and remains as one of the most damaging diseases of wheat worldwide. Currently, one of the major challenges of wheat geneticists and breeders is to accumulate major genes for durability of rust resistance called “slow rusting” genes using marker-assisted selection (MAS). Until now, eight genes (Lr34/Yr18, Lr46/Yr29, Lr67/Yr46, Lr68, Lr74, Lr75, Lr77, and Lr78) conferring resistance against multiple fungal pathogens have been identified in wheat gene pool and the molecular markers were developed for them. In MAS practice, it is a common problem that cultivars exhibiting desirable marker genotypes may not necessarily have the targeted genes or alleles and vice versa, which is known as “false positives.” The aim of this study was to compare the available four markers: Xwmc44, Xgwm259, Xbarc80, and csLV46G22 markers (not published yet), for the identification of the Lr46/Yr29 loci in 73 genotypes of wheat, which were reported as sources of various “slow rusting” genes, including 60 with confirmed Lr46/Yr29 gene, reported in the literature. This research revealed that csLV46G22 together with Xwmc44 is most suitable for the identification of resistance allele of the Lr46/Yr29 gene; however, there is a need to clone the Lr46/Yr29 loci to identify and verify the allelic variation of the gene and the function.

Identification of leaf rust resistance genes Lr34 and Lr46 in common wheat (Triticum aestivum L. ssp. aestivum) lines of different origin using multiplex PCR

Leaf rust caused by the fungus Puccinia recondita f. sp. tritici is one of the most dangerous diseases of common wheat. Infections caused by fungal pathogens reduce the quantity and quality of yields of many cereal species. The most effective method to limit plant infection is to use cultivars that show rust resistance. Genetically conditioned horizontal-type resistance (racial-nonspecific) is a desirable trait because it is characterized by more stable expression compared to major (R) genes that induce racially specific resistance, often overcome by pathogens. Horizontal resistance is conditioned by the presence of slow rust genes, which include genes Lr34 and Lr46. This study aimed to identify markers linked to both genes in 64 common wheat lines and to develop multiplex PCR reaction conditions that were applied to identify both genes simultaneously. The degree of infestation of the analyzed lines was also assessed in field conditions during the growing season of 2017 and 2018. Simple sequence repeat anchored-polymerase chain reaction (SSR-PCR) marker csLV was identified during analysis in line PHR 4947. The presence of a specific sequence has also been confirmed in multiplex PCR analyses. In addition to gene Lr34, gene Lr46 was identified in this genotype. Lines PHR 4947 and PHR 4819 were characterized by the highest leaf rust resistance in field conditions. During STS-PCR analyses, the marker wmc44 of gene Lr46 was identified in most of the analyzed lines. This marker was not present in the following genotypes: PHR 4670, PHR 4800, PHR 4859, PHR 4907, PHR 4922, PHR 4949, PHR 4957, PHR 4995, and PHR 4997. The presence of a specific sequence has also been confirmed in multiplex PCR analyses. Genotypes carrying the markers of the analyzed gene showed good resistance to leaf rust in field conditions in both 2017 and 2018. Research has demonstrated that marker assisted selection (MAS) and multiplex PCR techniques are excellent tools for selecting genotypes resistant to leaf rust.

Molecular mapping of a new recessive wheat leaf rust resistance gene originating from Triticum spelta

TSD276-2, a wheat genetic stock derived from the cross Agra Local/T. spelta 276 showed broad spectrum resistance against leaf rust pathogen. Genetic analysis was undertaken using F₁, F₂, F_2:3 and BC₁F₁ generations derived from the cross TSD276-2/Agra Local. The results revealed a single recessive gene for leaf rust resistance, tentatively named as LrTs_276-2, in TSD276-2. Molecular mapping of leaf rust resistance gene LrTs_276-2 in TSD276-2 was done using SNP-based PCR and SSR markers. For Bulked Segregant Analysis (BSA), two bulks viz. resistant bulk and susceptible bulk, and the parents TSD276-2 and Agra Local were genotyped for SNPs using AFFYMETRIX 35K Wheat Breeders' AXIOM array. T. spelta 276 was also genotyped and used as a check. BSA indicated that the gene for leaf rust resistance in TSD276-2 is located on chromosome arm 1DS. Putatively linked SNPs on chromosome arm 1DS were converted into PCR-based markers. Polymorphic SSR markers on chromosome arm 1DS were also identified. Final linkage map was constructed using one SNP-based PCR and three SSR markers. The rust reaction and chromosomal location suggest that LrTs_276-2 is a new leaf rust resistance gene which may be useful in broadening the genetic base of leaf rust resistance in wheat.

A Meta-Analysis of Quantitative Trait Loci Associated with Multiple Disease Resistance in Rice (Oryza sativa L.)

Rice blast, sheath blight and bacterial leaf blight are major rice diseases found worldwide. The development of resistant cultivars is generally perceived as the most effective way to combat these diseases. Plant disease resistance is a polygenic trait where a combinatorial effect of major and minor genes affects this trait. To locate the source of this trait, various quantitative trait loci (QTL) mapping studies have been performed in the past two decades. However, investigating the congruency between the reported QTL is a daunting task due to the heterogeneity amongst the QTLs studied. Hence, the aim of our study is to integrate the reported QTLs for resistance against rice blast, sheath blight and bacterial leaf blight and objectively analyze and consolidate the location of QTL clusters in the chromosomes, reducing the QTL intervals and thus identifying candidate genes within the selected meta-QTL. A total of twenty-seven studies for resistance QTLs to rice blast (8), sheath blight (15) and bacterial leaf blight (4) was compiled for QTL projection and analyses. Cumulatively, 333 QTLs associated with rice blast (114), sheath blight (151) and bacterial leaf blight (68) resistance were compiled, where 303 QTLs could be projected onto a consensus map saturated with 7633 loci. Meta-QTL analysis on 294 QTLs yielded 48 meta-QTLs, where QTLs with membership probability lower than 60% were excluded, reducing the number of QTLs within the meta-QTL to 274. Further, three meta-QTL regions (MQTL2.5, MQTL8.1 and MQTL9.1) were selected for functional analysis on the basis that MQTL2.5 harbors the highest number of QTLs; meanwhile, MQTL8.1 and MQTL9.1 have QTLs associated with all three diseases mentioned above. The functional analysis allows for determination of enriched gene ontology and resistance gene analogs (RGAs) and other defense-related genes. To summarize, MQTL2.5, MQTL8.1 and MQTL9.1 have a considerable number of R-genes that account for 10.21%, 4.08% and 6.42% of the total genes found in these meta-QTLs, respectively. Defense genes constitute around 3.70%, 8.16% and 6.42% of the total number of genes in MQTL2.5, MQTL8.1 and MQTL9.1, respectively. This frequency is higher than the total frequency of defense genes in the rice genome, which is 0.0096% (167 defense genes/17,272 total genes). The integration of the QTLs facilitates the identification of QTL hotspots for rice blast, sheath blight and bacterial blight resistance with reduced intervals, which helps to reduce linkage drag in breeding. The candidate genes within the promising regions could be utilized for improvement through genetical engineering.

The progress of leaf rust research in wheat

Leaf rust (also called brown rust) in wheat, caused by fungal pathogen Puccinia triticina Erikss. (Pt) is one of the major constraints in wheat production worldwide. Pt is widespread with diverse population structure and undergoes rapid evolution to produce new virulent races against resistant cultivars that are regularly developed to provide resistance against the prevailing races of the pathogen. Occasionally, the disease may also take the shape of an epidemic in some wheat-growing areas causing major economic losses. In the recent past, substantial progress has been made in characterizing the sources of leaf rust resistance including non-host resistance (NHR). Progress has also been made in elucidating the population biology of Pt and the mechanisms of wheat-Pt interaction. So far, ∼80 leaf rust resistance genes (Lr genes) have been identified and characterized; some of them have also been used for the development of resistant wheat cultivars. It has also been shown that a gene-for-gene relationship exists between individual wheat Lr genes and the corresponding Pt Avr genes so that no Lr gene can provide resistance unless the prevailing race of the pathogen carries the corresponding Avr gene. Several Lr genes have also been cloned and their products characterized, although no Avr gene corresponding a specific Lr gene has so far been identified. However, several candidate effectors for Pt have been identified and functionally characterized using genome-wide analyses, transcriptomics, RNA sequencing, bimolecular fluorescence complementation (BiFC), virus-induced gene silencing (VIGS), transient expression and other approaches. This review summarizes available information on different aspects of the pathogen Pt, genetics/genomics of leaf rust resistance in wheat including cloning and characterization of Lr genes and epigenetic regulation of disease resistance.

A Comparative Transcriptomic and Proteomic Analysis of Hexaploid Wheat's Responses to Colonization by Bacillus velezensis and Gaeumannomyces graminis, Both Separately and Combined

Tritrophic interactions involving a biocontrol agent, a pathogen, and a plant have been analyzed predominantly from the perspective of the biocontrol agent. To explore the adaptive strategies of wheat in response to beneficial, pathogenic, and combined microorganisms, we performed the first comprehensive transcriptomic, proteomic, and biochemical analysis in wheat roots after exposure to Bacillus velezensis CC09, Gaeumannomyces graminis var. tritici, and their combined colonization, respectively. The transcriptional or translational programming of wheat roots inoculated with beneficial B. velezensis showed mild alterations compared with that of pathogenic G. graminis var. tritici. However, the combination of B. velezensis and G. graminis var. tritici activated a larger transcriptional or translational program than for each single microorganism, although the gene expression pattern was similar to that of individual infection by G. graminis var. tritici, suggesting a prioritization of defense against G. graminis var. tritici infection. Surprisingly, pathogen-associated molecular pattern-triggered immunity and effector-triggered immunity made wheat pretreated with B. velezensis more sensitive to subsequent G. graminis var. tritici infection. Additionally, B. velezensis triggered a salicylic acid (SA)-dependent mode of induced systemic resistance that resembles pathogen-induced systemic acquired resistance. Wheat plants mainly depend on SA-mediated resistance, and not that mediated by jasmonic acid (JA), against the necrotrophic pathogen G. graminis var. tritici. Moreover, SA-JA interactions resulted in antagonistic effects regardless of the type of microorganisms in wheat. Further enhancement of SA-dependent defense responses such as lignification to the combined infection was shown to reduce the level of induced JA-dependent defense against subsequent infection with G. graminis var. tritici. Altogether, our results demonstrate how the hexaploid monocot wheat responds to beneficial or pathogenic microorganisms and prolongs the onset of take-all disease through modulation of cell reprogramming and signaling events.

Inferring defense-related gene families in Arabidopsis and wheat

BACKGROUND

A large number of disease resistance genes or QTLs in crop plants are identified through conventional genetics and genomic tools, but their functional or molecular characterization remains costly, labor-intensive and inaccurate largely due to the lack of deep sequencing of large and complex genomes of many important crops such as allohexaploid wheat (Triticum aestivum L.). On the other hand, gene annotation and relevant genomic resources for disease resistance and other defense-related traits are more abundant in model plant Arabidopsis (Arabidopsis thaliana). The objectives of this study are (i) to infer homology of defense-related genes in Arabidopsis and wheat and (ii) to classify these homologous genes into different gene families.

RESULTS

We employed three bioinformatics and genomics approaches to identifying candidate genes known to affect plant defense and to classifying these protein-coding genes into different gene families in Arabidopsis. These approaches predicted up to 1790 candidate genes in 11 gene families for Arabidopsis defense to biotic stresses. The 11 gene families included ABC, NLR and START, the three families that are already known to confer rust resistance in wheat, and eight new families. The distributions of predicted SNPs for individual rust resistance genes were highly skewed towards specific gene families, including eight one-to-one uniquely matched pairs: Lr21-NLR, Lr34-ABC, Lr37-START, Sr2-Cupin, Yr24-Transcription factor, Yr26-Transporter, Yr36-Kinase and Yr53-Kinase. Two of these pairs, Lr21-NLR and Lr34-ABC, are expected because Lr21 and Lr34 are well known to confer race-specific and race-nonspecific resistance to leaf rust (Puccinia triticina) and they encode NLR and ABC proteins.

CONCLUSIONS

Our inference of 11 known and new gene families enhances current understanding of functional diversity with defense-related genes in genomes of model plant Arabidopsis and cereal crop wheat. Our comparative genomic analysis of Arabidopsis and wheat genomes is complementary to the conventional map-based or marker-based approaches for identification of genes or QTLs for rust resistance genes in wheat and other cereals. Race-specific and race-nonspecific candidate genes predicted by our study may be further tested and combined in breeding for durable resistance to wheat rusts and other pathogens.