High-copy transposons from a pathogen give rise to a conserved sRNA family with a novel host immunity target.

Small RNAs (sRNAs) are involved in gene silencing in multiple ways including through cross-kingdom transfers from parasites to their hosts. Little is known about the evolutionary mechanisms enabling eukaryotic microbes to evolve functional mimics of host small regulatory RNAs. Here, we describe the identification and functional characterization of SINE_sRNA1, a sRNA family derived from highly abundant SINE retrotransposons in the genome of the wheat powdery mildew pathogen. SINE_sRNA1 is encoded by a sequence motif that is conserved in multiple SINE families and corresponds to a functional plant miRNA mimic targeting Tae_AP1, a wheat gene encoding an aspartic protease only found in monocots. Tae_AP1 has a novel function enhancing both pattern-triggered immunity (PTI) and effector-triggered immunity (ETI) thus contributing to the cross activation of plant defenses. We conclude that SINE_sRNA1 and Tae_AP1 are functional innovations suggesting the contribution of transposons to the evolutionary arms race between a parasite and its host.

The Zymoseptoria tritici Avirulence Factor AvrStb6 Accumulates in Hyphae Close to Stomata and Triggers a Wheat Defense Response Hindering Fungal Penetration.

Zymoseptoria tritici, the causal agent of Septoria tritici blotch, is one of Europe's most damaging wheat pathogens, causing significant economic losses. Genetic resistance is a common strategy to control the disease, Stb6 being a resistance gene used for more than 100 years in Europe. This study investigates the molecular mechanisms underlying Stb6-mediated resistance. Utilizing confocal microscopy imaging, we determined that Z. tritici epiphytic hyphae mainly accumulate the corresponding avirulence factor AvrStb6 in close proximity to stomata. Consequently, the progression of AvrStb6-expressing avirulent strains is hampered during penetration. The fungal growth inhibition co-occurs with a transcriptional reprogramming in wheat characterized by an induction of immune responses, genes involved in stomatal regulation, and cell wall-related genes. Overall, we shed light on the gene-for-gene resistance mechanisms in the wheat-Z. tritici pathosystem at the cytological and transcriptomic level, and our results highlight that stomatal penetration is a critical process for pathogenicity and resistance. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Botryosphaeriaceae gene machinery: Correlation between diversity and virulence.

The Botryosphaeriaceae family comprises numerous fungal pathogens capable of causing economically meaningful diseases in a wide range of crops. Many of its members can live as endophytes and turn into aggressive pathogens following the onset of environmental stress events. Their ability to cause disease may rely on the production of a broad set of effectors, such as cell wall-degrading enzymes, secondary metabolites, and peptidases. Here, we conducted comparative analyses of 41 genomes representing six Botryosphaeriaceae genera to provide insights into the genetic features linked to pathogenicity and virulence. We show that these Botryosphaeriaceae genomes possess a large diversity of carbohydrate-active enzymes (CAZymes; 128 families) and peptidases (45 families). Botryosphaeria, Neofusicoccum, and Lasiodiplodia presented the highest number of genes encoding CAZymes involved in the degradation of the plant cell wall components. The genus Botryosphaeria also exhibited the highest abundance of secreted CAZymes and peptidases. Generally, the secondary metabolites gene cluster profile was consistent in the Botryosphaeriaceae family, except for Diplodia and Neoscytalidium. At the strain level, Neofusicoccum parvum NpBt67 stood out among all the Botryosphaeriaceae genomes, presenting a higher number of secretome constituents. In contrast, the Diplodia strains showed the lowest richness of the pathogenicity- and virulence-related genes, which may correlate with their low virulence reported in previous studies. Overall, these results contribute to a better understanding of the mechanisms underlying pathogenicity and virulence in remarkable Botryosphaeriaceae species. Our results also support that Botryosphaeriaceae species could be used as an interesting biotechnological tool for lignocellulose fractionation and bioeconomy.

A survey of lineage-specific genes in Triticeae reveals de novo gene evolution from genomic raw material.

Diploid plant genomes typically contain ~35,000 genes, almost all belonging to highly conserved gene families. Only a small fraction are lineage-specific, which are found in only one or few closely related species. Little is known about how genes arise de novo in plant genomes and how often this occurs; however, they are believed to be important for plants diversification and adaptation. We developed a pipeline to identify lineage-specific genes in Triticeae, using newly available genome assemblies of wheat, barley, and rye. Applying a set of stringent criteria, we identified 5942 candidate Triticeae-specific genes (TSGs), of which 2337 were validated as protein-coding genes in wheat. Differential gene expression analyses revealed that stress-induced wheat TSGs are strongly enriched in putative secreted proteins. Some were previously described to be involved in Triticeae non-host resistance and cold response. Additionally, we show that 1079 TSGs have sequence homology to transposable elements (TEs), ~68% of them deriving from regulatory non-coding regions of Gypsy retrotransposons. Most importantly, we demonstrate that these TSGs are enriched in transmembrane domains and are among the most highly expressed wheat genes overall. To summarize, we conclude that de novo gene formation is relatively rare and that Triticeae probably possess ~779 lineage-specific genes per haploid genome. TSGs, which respond to pathogen and environmental stresses, may be interesting candidates for future targeted resistance breeding in Triticeae. Finally, we propose that non-coding regions of TEs might provide important genetic raw material for the functional innovation of TM domains and the evolution of novel secreted proteins.

Ancient variation of the AvrPm17 gene in powdery mildew limits the effectiveness of the introgressed rye Pm17 resistance gene in wheat.

Introgressions of chromosomal segments from related species into wheat are important sources of resistance against fungal diseases. The durability and effectiveness of introgressed resistance genes upon agricultural deployment is highly variable-a phenomenon that remains poorly understood, as the corresponding fungal avirulence genes are largely unknown. Until its breakdown, the Pm17 resistance gene introgressed from rye to wheat provided broad resistance against powdery mildew (Blumeria graminis). Here, we used quantitative trait locus (QTL) mapping to identify the corresponding wheat mildew avirulence effector AvrPm17. It is encoded by two paralogous genes that exhibit signatures of reoccurring gene conversion events and are members of a mildew sublineage specific effector cluster. Extensive haplovariant mining in wheat mildew and related sublineages identified several ancient virulent AvrPm17 variants that were present as standing genetic variation in wheat powdery mildew prior to the Pm17 introgression, thereby paving the way for the rapid breakdown of the Pm17 resistance. QTL mapping in mildew identified a second genetic component likely corresponding to an additional resistance gene present on the 1AL.1RS translocation carrying Pm17. This gene remained previously undetected due to suppressed recombination within the introgressed rye chromosomal segment. We conclude that the initial effectiveness of 1AL.1RS was based on simultaneous introgression of two genetically linked resistance genes. Our results demonstrate the relevance of pathogen-based genetic approaches to disentangling complex resistance loci in wheat. We propose that identification and monitoring of avirulence gene diversity in pathogen populations become an integral part of introgression breeding to ensure effective and durable resistance in wheat.

Wheat Pm4 resistance to powdery mildew is controlled by alternative splice variants encoding chimeric proteins.

Crop breeding for resistance to pathogens largely relies on genes encoding receptors that confer race-specific immunity. Here, we report the identification of the wheat Pm4 race-specific resistance gene to powdery mildew. Pm4 encodes a putative chimeric protein of a serine/threonine kinase and multiple C2 domains and transmembrane regions, a unique domain architecture among known resistance proteins. Pm4 undergoes constitutive alternative splicing, generating two isoforms with different protein domain topologies that are both essential for resistance function. Both isoforms interact and localize to the endoplasmatic reticulum when co-expressed. Pm4 reveals additional diversity of immune receptor architecture to be explored for breeding and suggests an endoplasmatic reticulum-based molecular mechanism of Pm4-mediated race-specific resistance.

Chromosome-scale genome assembly provides insights into rye biology, evolution and agronomic potential.

Rye (Secale cereale L.) is an exceptionally climate-resilient cereal crop, used extensively to produce improved wheat varieties via introgressive hybridization and possessing the entire repertoire of genes necessary to enable hybrid breeding. Rye is allogamous and only recently domesticated, thus giving cultivated ryes access to a diverse and exploitable wild gene pool. To further enhance the agronomic potential of rye, we produced a chromosome-scale annotated assembly of the 7.9-gigabase rye genome and extensively validated its quality by using a suite of molecular genetic resources. We demonstrate applications of this resource with a broad range of investigations. We present findings on cultivated rye's incomplete genetic isolation from wild relatives, mechanisms of genome structural evolution, pathogen resistance, low-temperature tolerance, fertility control systems for hybrid breeding and the yield benefits of rye-wheat introgressions.

Identification of specificity-defining amino acids of the wheat immune receptor Pm2 and powdery mildew effector AvrPm2.

Plant nucleotide-binding leucine-rich repeat receptors (NLRs) act as intracellular sensors for pathogen-derived effector proteins and trigger an immune response, frequently resulting in the hypersensitive cell death response (HR) of the infected host cell. The wheat (Triticum aestivum) NLR Pm2 confers resistance against the fungal pathogen Blumeria graminis f. sp. tritici (Bgt) if the isolate contains the specific RNase-like effector AvrPm2. We identified and isolated seven new Pm2 alleles (Pm2e-i) in the wheat D-genome ancestor Aegilops tauschii and two new natural AvrPm2 haplotypes from Bgt. Upon transient co-expression in Nicotiana benthamiana, we observed a variant-specific HR of the Pm2 variants Pm2a and Pm2i towards AvrPm2 or its homolog from the AvrPm2 effector family, BgtE-5843, respectively. Through the introduction of naturally occurring non-synonymous single nucleotide polymorphisms and structure-guided mutations, we identified single amino acids in both the wheat NLR Pm2 and the fungal effector proteins AvrPm2 and BgtE-5843 responsible for the variant-specific HR of the Pm2 variants. Exchanging these amino acids led to a modified HR of the Pm2-AvrPm2 interaction and allowed the identification of the effector head epitope, a 20-amino-acid long unit of AvrPm2 involved in the HR. Swapping of the AvrPm2 head epitope to the non-HR-triggering AvrPm2 family member BgtE-5846 led to gain of a HR by Pm2a. Our study presents a molecular approach to identify crucial effector surface structures involved in the HR and demonstrates that natural and induced diversity in an immune receptor and its corresponding effectors can provide the basis for understanding and modifying NLR-effector specificity.

Domestication of High-Copy Transposons Underlays the Wheat Small RNA Response to an Obligate Pathogen.

Plant genomes have evolved several evolutionary mechanisms to tolerate and make use of transposable elements (TEs). Of these, transposon domestication into cis-regulatory and microRNA (miRNA) sequences is proposed to contribute to abiotic/biotic stress adaptation in plants. The wheat genome is derived at 85% from TEs, and contains thousands of miniature inverted-repeat transposable elements (MITEs), whose sequences are particularly prone for domestication into miRNA precursors. In this study, we investigate the contribution of TEs to the wheat small RNA immune response to the lineage-specific, obligate powdery mildew pathogen. We show that MITEs of the Mariner superfamily contribute the largest diversity of miRNAs to the wheat immune response. In particular, MITE precursors of miRNAs are wide-spread over the wheat genome, and highly conserved copies are found in the Lr34 and QPm.tut-4A mildew resistance loci. Our work suggests that transposon domestication is an important evolutionary force driving miRNA functional innovation in wheat immunity.

The AvrPm3-Pm3 effector-NLR interactions control both race-specific resistance and host-specificity of cereal mildews on wheat.

The wheat Pm3 resistance gene against the powdery mildew pathogen occurs as an allelic series encoding functionally different immune receptors which induce resistance upon recognition of isolate-specific avirulence (AVR) effectors from the pathogen. Here, we describe the identification of five effector proteins from the mildew pathogens of wheat, rye, and the wild grass Dactylis glomerata, specifically recognized by the PM3B, PM3C and PM3D receptors. Together with the earlier identified AVRPM3A2/F2, the recognized AVRs of PM3B/C, (AVRPM3B2/C2), and PM3D (AVRPM3D3) belong to a large group of proteins with low sequence homology but predicted structural similarities. AvrPm3b2/c2 and AvrPm3d3 are conserved in all tested isolates of wheat and rye mildew, and non-host infection assays demonstrate that Pm3b, Pm3c, and Pm3d are also restricting the growth of rye mildew on wheat. Furthermore, divergent AVR homologues from non-adapted rye and Dactylis mildews are recognized by PM3B, PM3C, or PM3D, demonstrating their involvement in host specificity.

Fungal resistance mediated by maize wall-associated kinase ZmWAK-RLK1 correlates with reduced benzoxazinoid content.

Wall-associated kinases (WAKs) have recently been identified as major components of fungal and bacterial disease resistance in several cereal crop species. However, the molecular mechanisms of WAK-mediated resistance remain largely unknown. Here, we investigated the function of the maize gene ZmWAK-RLK1 (Htn1) that confers quantitative resistance to northern corn leaf blight (NCLB) caused by the hemibiotrophic fungal pathogen Exserohilum turcicum. ZmWAK-RLK1 was found to localize to the plasma membrane and its presence resulted in a modification of the infection process by reducing pathogen penetration into host tissues. A large-scale transcriptome analysis of near-isogenic lines (NILs) differing for ZmWAK-RLK1 revealed that several differentially expressed genes are involved in the biosynthesis of the secondary metabolites benzoxazinoids (BXs). The contents of several BXs including DIM2 BOA-Glc were significantly lower when ZmWAK-RLK1 is present. DIM2 BOA-Glc concentration was significantly elevated in ZmWAK-RLK1 mutants with compromised NCLB resistance. Maize mutants that were affected in overall BXs biosynthesis or content of DIM2 BOA-Glc showed increased NCLB resistance. We conclude that Htn1-mediated NCLB resistance is associated with a reduction of BX secondary metabolites. These findings suggest a link between WAK-mediated quantitative disease resistance and changes in biochemical fluxes starting with indole-3-glycerol phosphate.

A chromosome-scale genome assembly reveals a highly dynamic effector repertoire of wheat powdery mildew.

Blumeria graminis f. sp. tritici (B.g. tritici) is the causal agent of the wheat powdery mildew disease. The highly fragmented B.g. tritici genome available so far has prevented a systematic analysis of effector genes that are known to be involved in host adaptation. To study the diversity and evolution of effector genes we produced a chromosome-scale assembly of the B.g. tritici genome. The genome assembly and annotation was achieved by combining long-read sequencing with high-density genetic mapping, bacterial artificial chromosome fingerprinting and transcriptomics. We found that the 166.6 Mb B.g. tritici genome encodes 844 candidate effector genes, over 40% more than previously reported. Candidate effector genes have characteristic local genomic organization such as gene clustering and enrichment for recombination-active regions and certain transposable element families. A large group of 412 candidate effector genes shows high plasticity in terms of copy number variation in a global set of 36 isolates and of transcription levels. Our data suggest that copy number variation and transcriptional flexibility are the main drivers for adaptation in B.g. tritici. The high repeat content may play a role in providing a genomic environment that allows rapid evolution of effector genes with selection as the driving force.

Non-parent of Origin Expression of Numerous Effector Genes Indicates a Role of Gene Regulation in Host Adaption of the Hybrid Triticale Powdery Mildew Pathogen.

Powdery mildew is an important disease of cereals. It is caused by one species, Blumeria graminis, which is divided into formae speciales each of which is highly specialized to one host. Recently, a new form capable of growing on triticale (B.g. triticale) has emerged through hybridization between wheat and rye mildews (B.g. tritici and B.g. secalis, respectively). In this work, we used RNA sequencing to study the molecular basis of host adaptation in B.g. triticale. We analyzed gene expression in three B.g. tritici isolates, two B.g. secalis isolates and two B.g. triticale isolates and identified a core set of putative effector genes that are highly expressed in all formae speciales. We also found that the genes differentially expressed between isolates of the same form as well as between different formae speciales were enriched in putative effectors. Their coding genes belong to several families including some which contain known members of mildew avirulence (Avr) and suppressor (Svr) genes. Based on these findings we propose that effectors play an important role in host adaptation that is mechanistically based on Avr-Resistance gene-Svr interactions. We also found that gene expression in the B.g. triticale hybrid is mostly conserved with the parent-of-origin, but some genes inherited from B.g. tritici showed a B.g. secalis-like expression. Finally, we identified 11 unambiguous cases of putative effector genes with hybrid-specific, non-parent of origin gene expression, and we propose that they are possible determinants of host specialization in triticale mildew. These data suggest that altered expression of multiple effector genes, in particular Avr and Svr related factors, might play a role in mildew host adaptation based on hybridization.

Distinct domains of the AVRPM3A2/F2 avirulence protein from wheat powdery mildew are involved in immune receptor recognition and putative effector function.

Recognition of the AVRPM3A2/F2 avirulence protein from powdery mildew by the wheat PM3A/F immune receptor induces a hypersensitive response after co-expression in Nicotiana benthamiana. The molecular determinants of this interaction and how they shape natural AvrPm3a2/f2 allelic diversity are unknown. We sequenced the AvrPm3a2/f2 gene in a worldwide collection of 272 mildew isolates. Using the natural polymorphisms of AvrPm3a2/f2 as well as sequence information from related gene family members, we tested 85 single-residue-altered AVRPM3A2/F2 variants with PM3A, PM3F and PM3FL456P/Y458H (modified for improved signaling) in Nicotiana benthamiana for effects on recognition. An intact AvrPm3a2/f2 gene was found in all analyzed isolates and the protein variant recognized by PM3A/F occurred globally at high frequencies. Single-residue alterations in AVRPM3A2/F2 mostly disrupted, but occasionally enhanced, the recognition response by PM3A, PM3F and PM3FL456P/Y458H . Residues enhancing hypersensitive responses constituted a protein domain separate from both naturally occurring polymorphisms and positively selected residues of the gene family. These results demonstrate the utility of using gene family sequence diversity to screen residues for their role in recognition. This approach identified a putative interaction surface in AVRPM3A2/F2 not polymorphic in natural alleles. We conclude that molecular mechanisms besides recognition drive AvrPm3a2/f2 diversification.

Cereal powdery mildew effectors: a complex toolbox for an obligate pathogen.

Cereal powdery mildews are major pathogens of cultivated monocot crops, and all are obligate biotrophic fungi that can only grow and reproduce on living hosts. This lifestyle is combined with extreme host specialization where every mildew subspecies (referred to as forma specialis) can only infect one plant species. Recently there has been much progress in our understanding of the possible roles effectors play in this complex host-pathogen interaction. Here, we review current knowledge on the origin, evolution, and mode of action of cereal mildew effectors, with a particular focus on recent advances in the identification of bona fide effectors and avirulence effector proteins from wheat and barley powdery mildews.

Rapid turnover of effectors in grass powdery mildew (Blumeria graminis).

BACKGROUND: Grass powdery mildew (Blumeria graminis, Ascomycota) is a major pathogen of cereal crops and has become a model organism for obligate biotrophic fungal pathogens of plants. The sequenced genomes of two formae speciales (ff.spp.), B.g. hordei and B.g. tritici (pathogens of barley and wheat), were found to be enriched in candidate effector genes (CEGs). Similar to other filamentous pathogens, CEGs in B. graminis are under positive selection. Additionally, effectors are more likely to have presence-absence polymorphisms than other genes among different strains. RESULTS: Here we identified effectors in the genomes of three additional host-specific lineages of B. graminis (B.g. poae, B.g. avenae and B.g. infecting Lolium) which diverged between 24 and 5 million years ago (Mya). We found that most CEGs in B. graminis are clustered in families and that most families are present in both reference genomes (B.g. hordei and B.g. tritici) and in the genomes of all three newly annotated lineages. We identified conserved protein domains including a novel lipid binding domain. The phylogenetic analysis showed that frequent gene duplications and losses shaped the diversity of the effector repertoires of the different lineages through their evolutionary history. We observed several lineage-specific expansions where large clades of CEGs originated in only one lineage from a single gene through repeated gene duplications. When we applied a birth-death model we found that the turnover rate (the rate at which genes are deleted and duplicated) of CEG families is much higher than for non-CEG families. The analysis of genomic context revealed that the immediate surroundings of CEGs are enriched in transposable elements (TE) which could play a role in the duplication and deletion of CEGs. CONCLUSIONS: The CEG repertoires of related pathogens diverged dramatically in short evolutionary times because of rapid turnover and of positive selection fixing non-synonymous mutations. While signatures of positive selection on effector sequences are the expected outcome of the evolutionary "arms race" between pathogen and plant immune system, it is more difficult to infer the mechanisms and evolutionary forces that maintained an extreme turnover rate in CEG families of B. graminis for several millions of years.

AvrPm2 encodes an RNase-like avirulence effector which is conserved in the two different specialized forms of wheat and rye powdery mildew fungus.

There is a large diversity of genetically defined resistance genes in bread wheat against the powdery mildew pathogen Blumeria graminis (B. g.) f. sp. tritici. Many confer race-specific resistance to this pathogen, but until now only the mildew avirulence gene AvrPm3a2/f2 that is recognized by Pm3a/f was known molecularly. We performed map-based cloning and genome-wide association studies to isolate a candidate for the mildew avirulence gene AvrPm2. We then used transient expression assays in Nicotiana benthamiana to demonstrate specific and strong recognition of AvrPm2 by Pm2. The virulent AvrPm2 allele arose from a conserved 12 kb deletion, while there is no protein sequence diversity in the gene pool of avirulent B. g. tritici isolates. We found one polymorphic AvrPm2 allele in B. g. triticale and one orthologue in B. g. secalis and both are recognized by Pm2. AvrPm2 belongs to a small gene family encoding structurally conserved RNase-like effectors, including Avra13 from B. g. hordei, the cognate Avr of the barley resistance gene Mla13. These results demonstrate the conservation of functional avirulence genes in two cereal powdery mildews specialized on different hosts, thus providing a possible explanation for successful introgression of resistance genes from rye or other grass relatives to wheat.

Hybridization of powdery mildew strains gives rise to pathogens on novel agricultural crop species.

Throughout the history of agriculture, many new crop species (polyploids or artificial hybrids) have been introduced to diversify products or to increase yield. However, little is known about how these new crops influence the evolution of new pathogens and diseases. Triticale is an artificial hybrid of wheat and rye, and it was resistant to the fungal pathogen powdery mildew (Blumeria graminis) until 2001 (refs. 1,2,3). We sequenced and compared the genomes of 46 powdery mildew isolates covering several formae speciales. We found that B. graminis f. sp. triticale, which grows on triticale and wheat, is a hybrid between wheat powdery mildew (B. graminis f. sp. tritici) and mildew specialized on rye (B. graminis f. sp. secalis). Our data show that the hybrid of the two mildews specialized on two different hosts can infect the hybrid plant species originating from those two hosts. We conclude that hybridization between mildews specialized on different species is a mechanism of adaptation to new crops introduced by agriculture.

Multiple Avirulence Loci and Allele-Specific Effector Recognition Control the Pm3 Race-Specific Resistance of Wheat to Powdery Mildew.

In cereals, several mildew resistance genes occur as large allelic series; for example, in wheat (Triticum aestivum and Triticum turgidum), 17 functional Pm3 alleles confer agronomically important race-specific resistance to powdery mildew (Blumeria graminis). The molecular basis of race specificity has been characterized in wheat, but little is known about the corresponding avirulence genes in powdery mildew. Here, we dissected the genetics of avirulence for six Pm3 alleles and found that three major Avr loci affect avirulence, with a common locus_1 involved in all AvrPm3-Pm3 interactions. We cloned the effector gene AvrPm3(a2/f2) from locus_2, which is recognized by the Pm3a and Pm3f alleles. Induction of a Pm3 allele-dependent hypersensitive response in transient assays in Nicotiana benthamiana and in wheat demonstrated specificity. Gene expression analysis of Bcg1 (encoded by locus_1) and AvrPm3 (a2/f2) revealed significant differences between isolates, indicating that in addition to protein polymorphisms, expression levels play a role in avirulence. We propose a model for race specificity involving three components: an allele-specific avirulence effector, a resistance gene allele, and a pathogen-encoded suppressor of avirulence. Thus, whereas a genetically simple allelic series controls specificity in the plant host, recognition on the pathogen side is more complex, allowing flexible evolutionary responses and adaptation to resistance genes.