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

Population genomics and molecular epidemiology of wheat powdery mildew in Europe

Agricultural diseases are a major threat to sustainable food production. Yet, for many pathogens we know exceptionally little about their epidemiological and population dynamics, and this knowledge gap is slowing the development of efficient control strategies. Here we study the population genomics and molecular epidemiology of wheat powdery mildew, a disease caused by the biotrophic fungus Blumeria graminis forma specialis tritici (Bgt). We sampled Bgt across two consecutive years, 2022 and 2023, and compiled a genomic dataset of 415 Bgt isolates from 22 countries in Europe and surrounding regions. We identified a single epidemic unit in the north of Europe, consisting of a highly homogeneous population. Conversely, the south of Europe hosts smaller local populations which are less interconnected. In addition, we show that the population structure can be largely predicted by the prevalent wind patterns. We identified several loci that were under selection in the recent past, including fungicide targets and avirulence genes. Some of these loci are common between populations, while others are not, suggesting different local selective pressures. We reconstructed the evolutionary history of one of these loci, AvrPm17, coding for an effector recognized by the wheat receptor Pm17. We found evidence for a soft sweep on standing genetic variation. Multiple AvrPm17 haplotypes, which can partially escape recognition by Pm17, spread rapidly throughout the continent upon its introduction in the early 2000s. We also identified a new virulent variant, which emerged more recently and can evade Pm17 resistance altogether. Overall, we highlight the potential of genomic surveillance in resolving the evolutionary and epidemiological dynamics of agricultural pathogens, as well as in guiding control strategies.

The dynamics of fungal genome organization and its impact on host adaptation and antifungal resistance

Fungi are a diverse kingdom characterized by remarkable genomic plasticity that facilitates pathogenicity and adaptation to adverse environmental conditions. In this review, we delve into the dynamic organization of fungal genomes and its implications for host adaptation and antifungal resistance. We examine key features and the heterogeneity of genomes across different fungal species, including but not limited to their chromosome content, DNA composition, distribution and arrangement of their content across chromosomes, and other major traits. We further highlight how this variability in genomic traits influences their virulence and adaptation to adverse conditions. Fungal genomes exhibit large variations in size, gene content, and structural features, such as the abundance of transposable elements (TEs), compartmentalization into gene-rich and TE-rich regions, and the presence or absence of dispensable chromosomes. Genomic structural variations are equally diverse in fungi, ranging from whole-chromosome duplications that may enhance tolerance to antifungal compounds, to targeted deletion of effector encoding genes that may promote virulence. Finally, the often-overlooked fungal mitochondrial genomes can also affect virulence and resistance to fungicides. Such and other features of fungal genome organization are reviewed and discussed in the context of host-microbe interactions and antifungal resistance.

A powdery mildew core effector protein targets the host endosome tethering complexes HOPS and CORVET in barley

Powdery mildew fungi are serious pathogens affecting many plant species. Their genomes encode extensive repertoires of secreted effector proteins that suppress host immunity. Here, we revised and analyzed the candidate secreted effector protein (CSEP) effectome of the powdery mildew fungus, Blumeria hordei (Bh). We identified seven putative effectors that are broadly conserved in powdery mildew species, suggesting that they are core effectors of these phytopathogens. We showed that one of these effectors, CSEP0214, interacts with the barley (Hordeum vulgare) vacuolar protein-sorting 18 (VPS18) protein, a shared component of the class C core vacuole/endosome tethering (CORVET) and homotypic fusion and protein-sorting (HOPS) endosomal tethering complexes that mediate fusion of early endosomes and multivesicular bodies, respectively, with the central vacuole. Overexpression of CSEP0214 and knockdown of either VPS18, HOPS-specific VPS41, or CORVET-specific VPS8 blocked the vacuolar pathway and the accumulation of the fluorescent vacuolar marker protein (SP)-RFP-AFVY in the endoplasmic reticulum. Moreover, CSEP0214 inhibited the interaction between VPS18 and VPS16, which are both shared components of CORVET as well as HOPS. Additionally, introducing CSEP0214 into barley leaf cells blocked the hypersensitive cell death response associated with resistance gene-mediated immunity, indicating that endomembrane trafficking is required for this process. CSEP0214 expression also prevented callose deposition in cell wall appositions at attack sites and encasements of fungal infection structures. Our results indicate that the powdery mildew core effector CSEP0214 is an essential suppressor of plant immunity.

Genomic Surveillance and Molecular Evolution of Fungicide Resistance in European Populations of Wheat Powdery Mildew

Fungicides are used in agriculture to manage fungal infections and maintain crop yield and quality. In Europe, their application on cereals increased drastically starting from the mid 1970s, contributing to a significant improvement in yields. However, extensive usage has led to the rapid evolution of resistant pathogen populations within just a few years of fungicide deployment. Here we focus on wheat powdery mildew, a disease caused by the ascomycete fungus Blumeria graminis forma specialis tritici (Bgt). Previous research on Bgt documented the emergence of resistance to different fungicides and identified various resistance mechanisms. Yet, the frequency, distribution and evolutionary dynamics of fungicide resistance in Bgt populations remain largely unexplored. In this study, we leveraged extensive sampling and whole-genome sequencing of Bgt populations in Europe and the Mediterranean to investigate the population genetics and molecular epidemiology of fungicide resistance towards five major fungicide classes. We analysed gene sequences and copy number variation of eight known fungicide target genes in 415 Bgt isolates sampled between 1980 and 2023. We observed that mutations conferring resistance to various fungicides increased in frequency over time and had distinct geographic distributions, probably due to diverse deployment of fungicides across different regions. For demethylation inhibitor fungicides, we identified multiple independent events of resistance emergence with distinct mutational profiles, and we tracked their rapid spread in the last decades. Overall, we revealed the evolutionary and epidemiological dynamics of fungicide resistance mutations in European Bgt populations. These results underscore the potential of genomic surveillance and population genetics to enhance our understanding of fungicide resistance.

From Lesions to Lessons: Two Decades of Filamentous Plant Pathogen Genomics

Many filamentous microorganisms, such as fungi and oomycetes, have evolved the ability to colonize plants and cause devastating crop diseases. Coevolutionary conflicts with their hosts have shaped the genomes of these plant pathogens. Over the past 20 years, genomics and genomics-enabled technologies have revealed remarkable diversity in genome size, architecture, and gene regulatory mechanisms. Technical and conceptual advances continue to provide novel insights into evolutionary dynamics, diversification of distinct genomic compartments, and facilitated molecular disease diagnostics. In this review, we discuss how genomics has advanced our understanding of genome organization and plant-pathogen coevolution and provide a perspective on future developments in the field. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Technical Advances Drive the Molecular Understanding of Effectors from Wheat and Barley Powdery Mildew Fungi

Pathogens manipulate host physiology through the secretion of virulence factors (effectors) to invade and proliferate on the host. The molecular functions of effectors inside plant hosts have been of interest in the field of molecular plant-microbe interactions. Obligate biotrophic pathogens, such as rusts and powdery mildews, cannot proliferate outside of plant hosts. In addition to the inhibition of the plant's immune components, these pathogens are under particular pressure to extract nutrients efficiently from the host. Understanding the molecular basis of infections mediated by obligate biotrophic pathogens is significant because of their impact in modern agriculture. In addition, powdery mildews serve as excellent models for obligate biotrophic cereal pathogens. Here, we summarize the current knowledge on the effectorome of the barley and wheat powdery mildews and putative molecular virulence functions of effectors. We emphasize the availability of comprehensive genomic, transcriptomic, and proteomic resources and discuss the methodological approaches used for identifying candidate effectors, assessing effector virulence traits, and identifying effector targets in the host. We highlight established and more recently employed methodologies, discuss limitations, and suggest additional strategies. We identify open questions and discuss how addressing them with currently available resources will enhance our understanding of Blumeria candidates for secretor effector proteins. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Recombination landscape and karyotypic variations revealed by linkage mapping in the grapevine downy mildew pathogen Plasmopara viticola

Plasmopara viticola, the causal agent of grapevine downy mildew, is a biotrophic oomycete engaged in a tight coevolutionary relationship with its host. Rapid adaptation of the pathogen is favored by annual sexual reproduction that generates genotypic diversity. With the aim of studying the recombination landscape across the P. viticola genome, we generated 2 half-sibling F1 progenies (N = 189 and 162). Using targeted SNP sequencing, between 1,405 and 1,894 markers were included in parental linkage maps, and a consensus map was obtained by integrating 4,509 markers. The reference genome could be assembled into 17 pseudochromosomes, anchoring 88% of its physical length. We observed a strong collinearity between parental genomes and extensive synteny with the downy mildew Peronospora effusa. In the consensus map, the median recombination rate was 13.8 cM/Mb. The local recombination rate was highly variable along chromosomes, and recombination was suppressed in putative centromeric regions. Recombination rate was found negatively correlated with repeats' coverage and positively correlated with gene coverage. However, genes encoding secreted proteins and putative effectors were underrepresented in highly recombining regions. In both progenies, about 5% of the individuals presented karyotypic anomalies. Aneuploidies and triploidies almost exclusively originated from the male-transmitted chromosomes. Triploids resulted from fertilization by diploid gametes, but also from dispermy. Obligatory sexual reproduction each year may explain the lower level of karyotypic variation in P. viticola compared to other oomycetes. The linkage maps will be useful to guide future de novo chromosome-scale assemblies of P. viticola genomes and to perform forward genetics.

Avirulence depletion assay: Combining R gene-mediated selection with bulk sequencing for rapid avirulence gene identification in wheat powdery mildew

Wheat production is threatened by multiple fungal pathogens, such as the wheat powdery mildew fungus (Blumeria graminis f. sp. tritici, Bgt). Wheat resistance breeding frequently relies on the use of resistance (R) genes that encode diverse immune receptors which detect specific avirulence (AVR) effectors and subsequently induce an immune response. While R gene cloning has accelerated recently, AVR identification in many pathogens including Bgt lags behind, preventing pathogen-informed deployment of resistance sources. Here we describe a new "avirulence depletion (AD) assay" for rapid identification of AVR genes in Bgt. This assay relies on the selection of a segregating, haploid F1 progeny population on a resistant host, followed by bulk sequencing, thereby allowing rapid avirulence candidate gene identification with high mapping resolution. In a proof-of-concept experiment we mapped the AVR component of the wheat immune receptor Pm3a to a 25 kb genomic interval in Bgt harboring a single effector, the previously described AvrPm3a2/f2. Subsequently, we applied the AD assay to map the unknown AVR effector recognized by the Pm60 immune receptor. We show that AvrPm60 is encoded by three tandemly arrayed, nearly identical effector genes that trigger an immune response upon co-expression with Pm60 and its alleles Pm60a and Pm60b. We furthermore provide evidence that Pm60 outperforms Pm60a and Pm60b through more efficient recognition of AvrPm60 effectors, suggesting it should be prioritized for wheat breeding. Finally, we show that virulence towards Pm60 is caused by simultaneous deletion of all AvrPm60 gene paralogs and that isolates lacking AvrPm60 are especially prevalent in the US thereby limiting the potential of Pm60 in this region. The AD assay is a powerful new tool for rapid and inexpensive AVR identification in Bgt with the potential to contribute to pathogen-informed breeding decisions for the use of novel R genes and regionally tailored gene deployment.

Recombination Rate Variation in Social Insects: An Adaptive Perspective

Social insects have the highest rates of meiotic recombination among Metazoa, but there is considerable variation within the Hymenoptera. We synthesize the literature to investigate several hypotheses for these elevated recombination rates. We reexamine the long-standing Red Queen hypothesis, considering how social aspects of immunity could lead to increases in recombination. We examine the possibility of positive feedback between gene duplication and recombination rate in the context of caste specialization. We introduce a novel hypothesis that recombination rate may be driven up by direct selection on recombination activity in response to increases in lifespan. Finally, we find that the role of population size in recombination rate evolution remains opaque, despite the long-standing popularity of this hypothesis. Moreover, our review emphasizes how the varied life histories of social insect species provide an effective framework for advancing a broader understanding of adaptively driven variation in recombination rates.

A fungal plant pathogen overcomes mlo-mediated broad-spectrum disease resistance by rapid gene loss

Hosts and pathogens typically engage in a coevolutionary arms race. This also applies to phytopathogenic powdery mildew fungi, which can rapidly overcome plant resistance and perform host jumps. Using experimental evolution, we show that the powdery mildew pathogen Blumeria hordei is capable of breaking the agriculturally important broad-spectrum resistance conditioned by barley loss-of-function mlo mutants. Partial mlo virulence of evolved B. hordei isolates is correlated with a distinctive pattern of adaptive mutations, including small-sized (c. 8-40 kb) deletions, of which one is linked to the de novo insertion of a transposable element. Occurrence of the mutations is associated with a transcriptional induction of effector protein-encoding genes that is absent in mlo-avirulent isolates on mlo mutant plants. The detected mutational spectrum comprises the same loci in at least two independently isolated mlo-virulent isolates, indicating convergent multigenic evolution. The mutational events emerged in part early (within the first five asexual generations) during experimental evolution, likely generating a founder population in which incipient mlo virulence was later stabilized by additional events. This work highlights the rapid dynamic genome evolution of an obligate biotrophic plant pathogen with a transposon-enriched genome.

A kinase fusion protein from Aegilops longissima confers resistance to wheat powdery mildew

Many disease resistance genes have been introgressed into wheat from its wild relatives. However, reduced recombination within the introgressed segments hinders the cloning of the introgressed genes. Here, we have cloned the powdery mildew resistance gene Pm13, which is introgressed into wheat from Aegilops longissima, using a method that combines physical mapping with radiation-induced chromosomal aberrations and transcriptome sequencing analysis of ethyl methanesulfonate (EMS)-induced loss-of-function mutants. Pm13 encodes a kinase fusion protein, designated MLKL-K, with an N-terminal domain of mixed lineage kinase domain-like protein (MLKL_NTD domain) and a C-terminal serine/threonine kinase domain bridged by a brace. The resistance function of Pm13 is validated through transient and stable transgenic complementation assays. Transient over-expression analyses in Nicotiana benthamiana leaves and wheat protoplasts reveal that the fragment Brace-Kinase122-476 of MLKL-K is capable of inducing cell death, which is dependent on a functional kinase domain and the three α-helices in the brace region close to the N-terminus of the kinase domain.

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, an sRNA family derived from highly abundant short interspersed nuclear element (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 microRNA (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), thereby 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. [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.

Wheat zinc finger protein TaZF interacts with both the powdery mildew AvrPm2 protein and the corresponding wheat Pm2a immune receptor

Plant defense responses to pathogens are induced after direct or indirect perception of effector proteins or their activity on host proteins. In fungal-plant interactions, relatively little is known about whether, in addition to avirulence effectors and immune receptors, other proteins contribute to specific recognition. The nucleotide-binding leucine-rich repeat (NLR) immune receptor Pm2a in wheat recognizes the fungal powdery mildew effector AvrPm2. We found that the predicted wheat zinc finger TaZF interacts with both the fungal avirulence protein AvrPm2 and the wheat NLR Pm2a. We further demonstrated that the virulent AvrPm2-H2 variant does not interact with TaZF. TaZF silencing in wheat resulted in a reduction but not a loss of Pm2a-mediated powdery mildew resistance. Interaction studies showed that the leucine-rich repeat domain of Pm2a is the mediator of the interaction with TaZF. TaZF recruits both Pm2a and AvrPm2 from the cytosol to the nucleus, resulting in nuclear localization of Pm2a, TaZF, and AvrPm2 in wheat. We propose that TaZF acts as a facilitator of Pm2a-dependent AvrPm2 effector recognition. Our findings highlight the importance of identifying effector host targets for characterization of NLR-mediated effector recognition.

Long-term and rapid evolution in powdery mildew fungi

The powdery mildew fungi (Erysiphaceae) are globally distributed plant pathogens with a range of more than 10,000 plant hosts. In this review, we discuss the long- and short-term evolution of these obligate biotrophic fungi and outline their diversity with respect to morphology, lifestyle, and host range. We highlight their remarkable ability to rapidly overcome plant immunity, evolve fungicide resistance, and broaden their host range, for example, through adaptation and hybridization. Recent advances in genomics and proteomics, particularly in cereal powdery mildews (genus Blumeria), provided first insights into mechanisms of genomic adaptation in these fungi. Transposable elements play key roles in shaping their genomes, where even close relatives exhibit diversified patterns of recent and ongoing transposon activity. These transposons are ubiquitously distributed in the powdery mildew genomes, resulting in a highly adaptive genome architecture lacking obvious regions of conserved gene space. Transposons can also be neofunctionalized to encode novel virulence factors, particularly candidate secreted effector proteins, which may undermine the plant immune system. In cereals like barley and wheat, some of these effectors are recognized by plant immune receptors encoded by resistance genes with numerous allelic variants. These effectors determine incompatibility ("avirulence") and evolve rapidly through sequence diversification and copy number variation. Altogether, powdery mildew fungi possess plastic genomes that enable their fast evolutionary adaptation towards overcoming plant immunity, host barriers, and chemical stress such as fungicides, foreshadowing future outbreaks, host range shifts and expansions as well as potential pandemics by these pathogens.

Powdery mildew effectors AVRA1 and BEC1016 target the ER J-domain protein HvERdj3B required for immunity in barley

The barley powdery mildew fungus, Blumeria hordei (Bh), secretes hundreds of candidate secreted effector proteins (CSEPs) to facilitate pathogen infection and colonization. One of these, CSEP0008, is directly recognized by the barley nucleotide-binding leucine-rich-repeat (NLR) receptor MLA1 and therefore is designated AVRA1. Here, we show that AVRA1 and the sequence-unrelated Bh effector BEC1016 (CSEP0491) suppress immunity in barley. We used yeast two-hybrid next-generation interaction screens (Y2H-NGIS), followed by binary Y2H and in planta protein-protein interactions studies, and identified a common barley target of AVRA1 and BEC1016, the endoplasmic reticulum (ER)-localized J-domain protein HvERdj3B. Silencing of this ER quality control (ERQC) protein increased Bh penetration. HvERdj3B is ER luminal, and we showed using split GFP that AVRA1 and BEC1016 translocate into the ER signal peptide-independently. Overexpression of the two effectors impeded trafficking of a vacuolar marker through the ER; silencing of HvERdj3B also exhibited this same cellular phenotype, coinciding with the effectors targeting this ERQC component. Together, these results suggest that the barley innate immunity, preventing Bh entry into epidermal cells, requires ERQC. Here, the J-domain protein HvERdj3B appears to be essential and can be regulated by AVRA1 and BEC1016. Plant disease resistance often occurs upon direct or indirect recognition of pathogen effectors by host NLR receptors. Previous work has shown that AVRA1 is directly recognized in the cytosol by the immune receptor MLA1. We speculate that the AVRA1 J-domain target being inside the ER, where it is inapproachable by NLRs, has forced the plant to evolve this challenging direct recognition.

Deciphering the Genomic Landscape and Virulence Mechanisms of the Wheat Powdery Mildew Pathogen Blumeria graminis f. sp. tritici Wtn1: Insights from Integrated Genome Assembly and Conidial Transcriptomics

A high-quality genome sequence from an Indian isolate of Blumeria graminis f. sp. tritici Wtn1, a persistent threat in wheat farming, was obtained using a hybrid method. The assembly of over 9.24 million DNA-sequence reads resulted in 93 contigs, totaling a 140.61 Mb genome size, potentially encoding 8480 genes. Notably, more than 73.80% of the genome, spanning approximately 102.14 Mb, comprises retro-elements, LTR elements, and P elements, influencing evolution and adaptation significantly. The phylogenomic analysis placed B. graminis f. sp. tritici Wtn1 in a distinct monocot-infecting clade. A total of 583 tRNA anticodon sequences were identified from the whole genome of the native virulent strain B. graminis f. sp. tritici, which comprises distinct genome features with high counts of tRNA anticodons for leucine (70), cysteine (61), alanine (58), and arginine (45), with only two stop codons (Opal and Ochre) present and the absence of the Amber stop codon. Comparative InterProScan analysis unveiled "shared and unique" proteins in B. graminis f. sp. tritici Wtn1. Identified were 7707 protein-encoding genes, annotated to different categories such as 805 effectors, 156 CAZymes, 6102 orthologous proteins, and 3180 distinct protein families (PFAMs). Among the effectors, genes like Avra10, Avrk1, Bcg-7, BEC1005, CSEP0105, CSEP0162, BEC1016, BEC1040, and HopI1 closely linked to pathogenesis and virulence were recognized. Transcriptome analysis highlighted abundant proteins associated with RNA processing and modification, post-translational modification, protein turnover, chaperones, and signal transduction. Examining the Environmental Information Processing Pathways in B. graminis f. sp. tritici Wtn1 revealed 393 genes across 33 signal transduction pathways. The key pathways included yeast MAPK signaling (53 genes), mTOR signaling (38 genes), PI3K-Akt signaling (23 genes), and AMPK signaling (21 genes). Additionally, pathways like FoxO, Phosphatidylinositol, the two-component system, and Ras signaling showed significant gene representation, each with 15-16 genes, key SNPs, and Indels in specific chromosomes highlighting their relevance to environmental responses and pathotype evolution. The SNP and InDel analysis resulted in about 3.56 million variants, including 3.45 million SNPs, 5050 insertions, and 5651 deletions within the whole genome of B. graminis f. sp. tritici Wtn1. These comprehensive genome and transcriptome datasets serve as crucial resources for understanding the pathogenicity, virulence effectors, retro-elements, and evolutionary origins of B. graminis f. sp. tritici Wtn1, aiding in developing robust strategies for the effective management of wheat powdery mildew.

Segmental duplications drive the evolution of accessory regions in a major crop pathogen

Many pathogens evolved compartmentalized genomes with conserved core and variable accessory regions (ARs) that carry effector genes mediating virulence. The fungal plant pathogen Fusarium oxysporum has such ARs, often spanning entire chromosomes. The presence of specific ARs influences the host range, and horizontal transfer of ARs can modify the pathogenicity of the receiving strain. However, how these ARs evolve in strains that infect the same host remains largely unknown. We defined the pan-genome of 69 diverse F. oxysporum strains that cause Fusarium wilt of banana, a significant constraint to global banana production, and analyzed the diversity and evolution of the ARs. Accessory regions in F. oxysporum strains infecting the same banana cultivar are highly diverse, and we could not identify any shared genomic regions and in planta-induced effectors. We demonstrate that segmental duplications drive the evolution of ARs. Furthermore, we show that recent segmental duplications specifically in accessory chromosomes cause the expansion of ARs in F. oxysporum. Taken together, we conclude that extensive recent duplications drive the evolution of ARs in F. oxysporum, which contribute to the evolution of virulence.

Mutagenesis of Wheat Powdery Mildew Reveals a Single Gene Controlling Both NLR and Tandem Kinase-Mediated Immunity

Blumeria graminis f. sp. tritici (Bgt) is a globally important fungal wheat pathogen. Some wheat genotypes contain powdery mildew resistance (Pm) genes encoding immune receptors that recognize specific fungal-secreted effector proteins, defined as avirulence (Avr) factors. Identifying Avr factors is vital for understanding the mechanisms, functioning, and durability of wheat resistance. Here, we present AvrXpose, an approach to identify Avr genes in Bgt by generating gain-of-virulence mutants on Pm genes. We first identified six Bgt mutants with gain of virulence on Pm3b and Pm3c. They all had point mutations, deletions or insertions of transposable elements within the corresponding AvrPm3b2/c2 gene or its promoter region. We further selected six mutants on Pm3a, aiming to identify the yet unknown AvrPm3a3 recognized by Pm3a, in addition to the previously described AvrPm3a2/f2. Surprisingly, Pm3a virulence in the obtained mutants was always accompanied by an additional gain of virulence on the unrelated tandem kinase resistance gene WTK4. No virulence toward 11 additional R genes tested was observed, indicating that the gain of virulence was specific for Pm3a and WTK4. Several independently obtained Pm3a-WTK4 mutants have mutations in Bgt-646, a gene encoding a putative, nonsecreted ankyrin repeat-containing protein. Gene expression analysis suggests that Bgt-646 regulates a subset of effector genes. We conclude that Bgt-646 is a common factor required for avirulence on both a specific nucleotide-binding leucine-rich repeat and a WTK immune receptor. Our findings suggest that, beyond effectors, another type of pathogen protein can control the race-specific interaction between powdery mildew and wheat. [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.

Two functional CC-NBS-LRR proteins from rye chromosome 6RS confer differential age-related powdery mildew resistance to wheat

Rye (Secale cereale), a valuable relative of wheat, contains abundant powdery mildew resistance (Pm) genes. Using physical mapping, transcriptome sequencing, barley stripe mosaic virus-induced gene silencing, ethyl methane sulfonate mutagenesis, and stable transformation, we isolated and validated two coiled-coil, nucleotide-binding site and leucine-rich repeat (CC-NBS-LRR) alleles, PmTR1 and PmTR3, located on rye chromosome 6RS from different triticale lines. PmTR1 confers age-related resistance starting from the three-leaf stage, whereas its allele, PmTR3, confers typical all-stage resistance, which may be associated with their differential gene expression patterns. Overexpression in Nicotiana benthamiana showed that the CC, CC-NBS, and CC-LRR fragments of PMTR1 induce cell death, whereas in PMTR3 the CC and full-length fragments perform this function. Luciferase complementation imaging and pull-down assays revealed distinct interaction activities between the CC and NBS fragments. Our study elucidates two novel rye-derived Pm genes and their derivative germplasm resources and provides novel insights into the mechanism of age-related resistance, which can aid the improvement of resistance against wheat powdery mildew.

Long noncoding RNAs emerge from transposon-derived antisense sequences and may contribute to infection stage-specific transposon regulation in a fungal phytopathogen

BACKGROUND

The genome of the obligate biotrophic phytopathogenic barley powdery mildew fungus Blumeria hordei is inflated due to highly abundant and possibly active transposable elements (TEs). In the absence of the otherwise common repeat-induced point mutation transposon defense mechanism, noncoding RNAs could be key for regulating the activity of TEs and coding genes during the pathogenic life cycle.

RESULTS

We performed time-course whole-transcriptome shotgun sequencing (RNA-seq) of total RNA derived from infected barley leaf epidermis at various stages of fungal pathogenesis and observed significant transcript accumulation and time point-dependent regulation of TEs in B. hordei. Using a manually curated consensus database of 344 TEs, we discovered phased small RNAs mapping to 104 consensus transposons, suggesting that RNA interference contributes significantly to their regulation. Further, we identified 5,127 long noncoding RNAs (lncRNAs) genome-wide in B. hordei, of which 823 originated from the antisense strand of a TE. Co-expression network analysis of lncRNAs, TEs, and coding genes throughout the asexual life cycle of B. hordei points at extensive positive and negative co-regulation of lncRNAs, subsets of TEs and coding genes.

CONCLUSIONS

Our work suggests that similar to mammals and plants, fungal lncRNAs support the dynamic modulation of transcript levels, including TEs, during pivotal stages of host infection. The lncRNAs may support transcriptional diversity and plasticity amid loss of coding genes in powdery mildew fungi and may give rise to novel regulatory elements and virulence peptides, thus representing key drivers of rapid evolutionary adaptation to promote pathogenicity and overcome host defense.

Virulence Associations and Global Context of AvrStb6 Genetic Diversity in Iranian Populations of Zymoseptoria tritici

Managing pathogen damage in wheat production is important for sustaining yields. Fungal plant pathogen genomes encode many small secreted proteins acting as effectors that play key roles in the successful colonization of host tissue and triggering host defenses. AvrStb6 is the first described Zymoseptoria tritici avirulence effector, which triggers Stb6-mediated immunity in the wheat host in a gene-for-gene manner. Evasion of major resistance factors such as Stb6 challenges deployment decisions on wheat cultivars. In this study, we analyzed the evolution of the AvrStb6 effector in Iranian isolates of Z. tritici. In total, 78 isolates were isolated and purified from 30 infected wheat specimens collected from the East Azerbaijan and Ardabil provinces of Iran. The pathogenicity of all isolates was evaluated on the susceptible wheat cultivar 'Tajan'. A subset of 40 isolates were also tested for pathogenicity on the resistant cultivar 'Shafir' carrying Stb6. Genetic diversity at the AvrStb6 locus was analyzed for 14 isolates covering the breadth of the observed disease severity. The AvrStb6 sequence variation was high, with virulent isolates carrying highly diverse AvrStb6 haplotypes. In an analysis including more than 1,000 additional AvrStb6 sequences from a global set of isolates, we found that virulent isolates carried AvrStb6 haplotypes either clustering with known virulent haplotypes on different continents or constituting previously unknown haplotypes. Furthermore, we found that AvrStb6 variants from avirulent isolates clustered with known avirulent genotypes from Europe. Our study highlights the relevance of AvrStb6 for Z. tritici virulence and the exceptional global diversity patterns of this effector.

Wheat powdery mildew resistance: from gene identification to immunity deployment

Powdery mildew is one of the most devastating diseases on wheat and is caused by the obligate biotrophic phytopathogen Blumeria graminis f. sp. tritici (Bgt). Due to the complexity of the large genome of wheat and its close relatives, the identification of powdery mildew resistance genes had been hampered for a long time until recent progress in large-scale sequencing, genomics, and rapid gene isolation techniques. Here, we describe and summarize the current advances in wheat powdery mildew resistance, emphasizing the most recent discoveries about the identification of genes conferring powdery mildew resistance and the similarity, diversity and molecular function of those genes. Multilayered resistance to powdery mildew in wheat could be used for counteracting Bgt, including durable, broad spectrum but partial resistance, as well as race-specific and mostly complete resistance mediated by nucleotide-binding and leucine rich repeat domain (NLR) proteins. In addition to the above mentioned layers, manipulation of susceptibility (S) and negative regulator genes may represent another layer that can be used for durable and broad-spectrum resistance in wheat. We propose that it is promising to develop effective and durable strategies to combat powdery mildew in wheat by simultaneous deployment of multilayered immunity.

An effector of Erysiphe necator translocates to chloroplasts and plasma membrane to suppress host immunity in grapevine

The powdery mildew (Erysiphe necator) is a prevalent pathogen hampering grapevine growth in the vineyard. An arsenal of candidate secreted effector proteins (CSEPs) was encoded in the E. necator genome, but it is largely unclear what role CSEPs plays during the E. necator infection. In the present study, we identified a secreted effector CSEP080 of E. necator, which was located in plant chloroplasts and plasma membrane. Transient expressing CSEP080 promotes plant photosynthesis and inhibits INF1-induced cell death in tobacco leaves. We found that CSEP080 was a necessary effector for the E. necator pathogenicity, which interacted with grapevine chloroplast protein VviB6f (cytochrome b6-f complex iron-sulfur subunit), affecting plant photosynthesis. Transient silencing VviB6f increased the plant hydrogen peroxide production, and the plant resistance to powdery mildew. In addition, CSEP080 manipulated the VviPE (pectinesterase) to promote pectin degradation. Our results demonstrated the molecular mechanisms that an effector of E. necator translocates to host chloroplasts and plasma membrane, which suppresses with the grapevine immunity system by targeting the chloroplast protein VviB6f to suppress hydrogen peroxide accumulation and manipulating VviPE to promote pectin degradation.

Spray-induced gene silencing to identify powdery mildew gene targets and processes for powdery mildew control

Spray-induced gene silencing (SIGS) is an emerging tool for crop pest protection. It utilizes exogenously applied double-stranded RNA to specifically reduce pest target gene expression using endogenous RNA interference machinery. In this study, SIGS methods were developed and optimized for powdery mildew fungi, which are widespread obligate biotrophic fungi that infect agricultural crops, using the known azole-fungicide target cytochrome P450 51 (CYP51) in the Golovinomyces orontii-Arabidopsis thaliana pathosystem. Additional screening resulted in the identification of conserved gene targets and processes important to powdery mildew proliferation: apoptosis-antagonizing transcription factor in essential cellular metabolism and stress response; lipid catabolism genes lipase a, lipase 1, and acetyl-CoA oxidase in energy production; and genes involved in manipulation of the plant host via abscisic acid metabolism (9-cis-epoxycarotenoid dioxygenase, xanthoxin dehydrogenase, and a putative abscisic acid G-protein coupled receptor) and secretion of the effector protein, effector candidate 2. Powdery mildew is the dominant disease impacting grapes and extensive powdery mildew resistance to applied fungicides has been reported. We therefore developed SIGS for the Erysiphe necator-Vitis vinifera system and tested six successful targets identified using the G. orontii-A. thaliana system. For all targets tested, a similar reduction in powdery mildew disease was observed between systems. This indicates screening of broadly conserved targets in the G. orontii-A. thaliana pathosystem identifies targets and processes for the successful control of other powdery mildew fungi. The efficacy of SIGS on powdery mildew fungi makes SIGS an exciting prospect for commercial powdery mildew control.

A chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

Structural polymorphisms within a common powdery mildew effector scaffold as a driver of coevolution with cereal immune receptors

In plants, host-pathogen coevolution often manifests in reciprocal, adaptive genetic changes through variations in host nucleotide-binding leucine-rich repeat immune receptors (NLRs) and virulence-promoting pathogen effectors. In grass powdery mildew (PM) fungi, an extreme expansion of a RNase-like effector family, termed RALPH, dominates the effector repertoire, with some members recognized as avirulence (AVR) effectors by cereal NLR receptors. We report the structures of the sequence-unrelated barley PM effectors AVRA6, AVRA7, and allelic AVRA10/AVRA22 variants, which are detected by highly sequence-related barley NLRs MLA6, MLA7, MLA10, and MLA22 and of wheat PM AVRPM2 detected by the unrelated wheat NLR PM2. The AVR effectors adopt a common scaffold, which is shared with the RNase T1/F1 family. We found striking variations in the number, position, and length of individual structural elements between RALPH AVRs, which is associated with a differentiation of RALPH effector subfamilies. We show that all RALPH AVRs tested have lost nuclease and synthetase activities of the RNase T1/F1 family and lack significant binding to RNA, implying that their virulence activities are associated with neo-functionalization events. Structure-guided mutagenesis identified six AVRA6 residues that are sufficient to turn a sequence-diverged member of the same RALPH subfamily into an effector specifically detected by MLA6. Similar structure-guided information for AVRA10 and AVRA22 indicates that MLA receptors detect largely distinct effector surface patches. Thus, coupling of sequence and structural polymorphisms within the RALPH scaffold of PMs facilitated escape from NLR recognition and potential acquisition of diverse virulence functions.

Two pathogen loci determine Blumeria graminis f. sp. tritici virulence to wheat resistance gene Pm1a

Blumeria graminis f. sp. tritici (Bgt) is a globally important fungal pathogen of wheat that can rapidly evolve to defeat wheat powdery mildew (Pm) resistance genes. Despite periodic regional deployment of the Pm1a resistance gene in US wheat production, Bgt strains that overcome Pm1a have been notably nonpersistent in the United States, while on other continents, they are more widely established. A genome-wide association study (GWAS) was conducted to map sequence variants associated with Pm1a virulence in 216 Bgt isolates from six countries, including the United States. A virulence variant apparently unique to Bgt isolates from the United States was detected in the previously mapped gene AvrPm1a (BgtE-5612) on Bgt chromosome 6; an in vitro growth assay suggested no fitness reduction associated with this variant. A gene on Bgt chromosome 8, Bgt-51526, was shown to function as a second determinant of Pm1a virulence, and despite < 30% amino acid identity, BGT-51526 and BGTE-5612 were predicted to share > 85% of their secondary structure. A co-expression study in Nicotiana benthamiana showed that BGTE-5612 and BGT-51526 each produce a PM1A-dependent hypersensitive response. More than one member of a B. graminis effector family can be recognized by a single wheat immune receptor, and a two-gene model is necessary to explain virulence to Pm1a.

The broad use of the Pm8 resistance gene in wheat resulted in hypermutation of the AvrPm8 gene in the powdery mildew pathogen

BACKGROUND

Worldwide wheat production is under constant threat by fast-evolving fungal pathogens. In the last decades, wheat breeding for disease resistance heavily relied on the introgression of chromosomal segments from related species as genetic sources of new resistance. The Pm8 resistance gene against the powdery mildew disease has been introgressed from rye into wheat as part of a large 1BL.1RS chromosomal translocation encompassing multiple disease resistance genes and yield components. Due to its high agronomic value, this translocation has seen continuous global use since the 1960s on large growth areas, even after Pm8 resistance was overcome by the powdery mildew pathogen. The long-term use of Pm8 at a global scale provided the unique opportunity to study the consequences of such extensive resistance gene application on pathogen evolution.

RESULTS

Using genome-wide association studies in a population of wheat mildew isolates, we identified the avirulence effector AvrPm8 specifically recognized by Pm8. Haplovariant mining in a global mildew population covering all major wheat growing areas of the world revealed 17 virulent haplotypes of the AvrPm8 gene that grouped into two functional categories. The first one comprised amino acid polymorphisms at a single position along the AvrPm8 protein, which we confirmed to be crucial for the recognition by Pm8. The second category consisted of numerous destructive mutations to the AvrPm8 open reading frame such as disruptions of the start codon, gene truncations, gene deletions, and interference with mRNA splicing. With the exception of a single, likely ancient, gain-of-virulence mutation found in mildew isolates around the world, all AvrPm8 virulence haplotypes were found in geographically restricted regions, indicating that they occurred recently as a consequence of the frequent Pm8 use.

CONCLUSIONS

In this study, we show that the broad and prolonged use of the Pm8 gene in wheat production worldwide resulted in a multitude of gain-of-virulence mechanisms affecting the AvrPm8 gene in the wheat powdery mildew pathogen. Based on our findings, we conclude that both standing genetic variation as well as locally occurring new mutations contributed to the global breakdown of the Pm8 resistance gene introgression.

Pathogenicity Variation in Two Genomes of Cercospora Species Causing Gray Leaf Spot in Maize

The gray leaf spots caused by Cercospora spp. severely affect the yield and quality of maize. However, the evolutionary relation and pathogenicity variation between species of the Cercospora genus is largely unknown. In this study, we constructed high-quality reference genomes by nanopore sequencing two Cercospora species, namely, C. zeae-maydis and C. zeina, with differing pathogenicity, collected from northeast (Liaoning [LN]) and southeast (Yunnan [YN]) China, respectively. The genome size of C. zeae-maydis-LN is 45.08 Mb, containing 10,839 annotated genes, whereas that of Cercospora zeina-YN is 42.18 Mb, containing 10,867 annotated genes, of which approximately 86.58% are common in the two species. The difference in their genome size is largely attributed to increased long terminal repeat retrotransposons of 3.8 Mb in total length in C. zeae-maydis-LN. There are 41 and 30 carbohydrate-binding gene subfamilies identified in C. zeae-maydis-LN and C. zeina-YN, respectively. A higher number of carbohydrate-binding families found in C. zeae-maydis-LN, and its unique CBM4, CBM37, and CBM66, in particular, may contribute to variation in pathogenicity between the two species, as the carbohydrate-binding genes are known to encode cell wall-degrading enzymes. Moreover, there are 114 and 107 effectors predicted, with 47 and 46 having unique potential pathogenicity in C. zeae-maydis-LN and C. zeina-YN, respectively. Of eight effectors randomly selected for pathogenic testing, five were found to inhibit cell apoptosis induced by Bcl-2-associated X. Taken together, our results provide genomic insights into variation in pathogenicity between C. zeae-maydis and C. zeina. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Whole-Genome Metalloproteases in the Wheat Sharp Eyespot Pathogen Rhizoctonia cerealis and a Role in Fungal Virulence

Rhizoctonia cerealis is the causal agent of sharp eyespot, a devastating disease of cereal crops including wheat. Several metalloproteases have been implicated in pathogenic virulence, but little is known about whole-genome metalloproteases in R. cerealis. In this study, a total of 116 metalloproteases-encoding genes were identified and characterized from the R. cerealis Rc207 genome. The gene expression profiles and phylogenetic relationship of 11 MEP36/fungalysin metalloproteases were examined during the fungal infection to wheat, and function of an upregulated secretory MEP36 named RcFL1 was validated. Of 11 MEP36 family metalloproteases, ten, except RcFL5, were predicted to be secreted proteins and nine encoding genes, but not RcFL5 and RcFL2, were expressed during the R. cerealis infection process. Phylogenetic analysis suggested that MEP36 metalloproteases in R. cerealis were closely related to those of Rhizoctonia solani but were remote to those of Bipolaris sorokiniana, Fusarium graminearum, F. pseudograminearum, and Pyricularia oryzae. Expression of RcFL1 was significantly upregulated during the infection process and induced plant cell death in wheat to promote the virulence of the pathogen. The MEP36 domain was necessary for the activities of RcFL1. Furthermore, RcFL1 could repress the expression of wheat genes coding for the chitin elicitor receptor kinase TaCERK1 and chitinases. These results suggest that this MEP36 metalloprotease RcFL1 may function as a virulence factor of R. cerealis through inhibiting host chitin-triggered immunity and chitinases. This study provides insights on pathogenic mechanisms of R. cerealis. RcFL1 likely is an important gene resource for improving resistance of wheat to R. cerealis through host-induced gene silencing strategy.

Global genomic analyses of wheat powdery mildew reveal association of pathogen spread with historical human migration and trade

The fungus Blumeria graminis f. sp. tritici causes wheat powdery mildew disease. Here, we study its spread and evolution by analyzing a global sample of 172 mildew genomes. Our analyses show that B.g. tritici emerged in the Fertile Crescent during wheat domestication. After it spread throughout Eurasia, colonization brought it to America, where it hybridized with unknown grass mildew species. Recent trade brought USA strains to Japan, and European strains to China. In both places, they hybridized with local ancestral strains. Thus, although mildew spreads by wind regionally, our results indicate that humans drove its global spread throughout history and that mildew rapidly evolved through hybridization.

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

Domesticated and wild wheat relatives provide an important source of new immune receptors for wheat resistance breeding against fungal pathogens. The durability of these resistance genes is variable and difficult to predict, yet it is crucial for effective resistance breeding. We identified a fungal effector protein recognized by an immune receptor introgressed from rye to wheat. We found that variants of the effector allowing the fungus to overcome the resistance are ancient. They were already present in the wheat powdery mildew gene pool before the introgression of the immune receptor and are therefore responsible for the rapid resistance breakdown. Our study demonstrates that the effort to identify durable resistance genes cannot be dissociated from studies of pathogen avirulence genes.

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.

Beyond Nuclear Ribosomal DNA Sequences: Evolution, Taxonomy, and Closest Known Saprobic Relatives of Powdery Mildew Fungi (Erysiphaceae) Inferred From Their First Comprehensive Genome-Scale Phylogenetic Analyses

Powdery mildew fungi (Erysiphaceae), common obligate biotrophic pathogens of many plants, including important agricultural and horticultural crops, represent a monophyletic lineage within the Ascomycota. Within the Erysiphaceae, molecular phylogenetic relationships and DNA-based species and genera delimitations were up to now mostly based on nuclear ribosomal DNA (nrDNA) phylogenies. This is the first comprehensive genome-scale phylogenetic analysis of this group using 751 single-copy orthologous sequences extracted from 24 selected powdery mildew genomes and 14 additional genomes from Helotiales, the fungal order that includes the Erysiphaceae. Representative genomes of all powdery mildew species with publicly available whole-genome sequencing (WGS) data that were of sufficient quality were included in the analyses. The 24 powdery mildew genomes included in the analysis represented 17 species belonging to eight out of 19 genera recognized within the Erysiphaceae. The epiphytic genera, all but one represented by multiple genomes, belonged each to distinct, well-supported lineages. Three hemiendophytic genera, each represented by a single genome, together formed the hemiendophytic lineage. Out of the 14 other taxa from the Helotiales, Arachnopeziza araneosa, a saprobic species, was the only taxon that grouped together with the 24 genome-sequenced powdery mildew fungi in a monophyletic clade. The close phylogenetic relationship between the Erysiphaceae and Arachnopeziza was revealed earlier by a phylogenomic study of the Leotiomycetes. Further analyses of powdery mildew and Arachnopeziza genomes may discover signatures of the evolutionary processes that have led to obligate biotrophy from a saprobic way of life. A separate phylogeny was produced using the 18S, 5.8S, and 28S nrDNA sequences of the same set of powdery mildew specimens and compared to the genome-scale phylogeny. The nrDNA phylogeny was largely congruent to the phylogeny produced using 751 orthologs. This part of the study has revealed multiple contamination and other quality issues in some powdery mildew genomes. We recommend that the presence of 28S, internal transcribed spacer (ITS), and 18S nrDNA sequences in powdery mildew WGS datasets that are identical to those determined by Sanger sequencing should be used to assess the quality of assemblies, in addition to the commonly used Benchmarking Universal Single-Copy Orthologs (BUSCO) values.

First Draft Genome Assemblies of Pleochaeta shiraiana and Phyllactinia moricola, Two Tree-Parasitic Powdery Mildew Fungi with Hemiendophytic Mycelia

Powdery mildew fungi (Erysiphaceae) are widespread obligate biotrophic plant pathogens. Thus, applying genetic and omics approaches to study these fungi remains a major challenge, particularly for species with hemiendophytic mycelium. These belong to a distinct phylogenetic lineage within the family Erysiphaceae. To date, only a single draft genome assembly is available for this clade, obtained for Leveillula taurica. Here, we generated the first draft genome assemblies of Pleochaeta shiraiana and Phyllactinia moricola, two tree-parasitic powdery mildew species with hemiendophytic mycelium, representing two genera that have not yet been investigated with genomics tools. The Pleochaeta shiraiana assembly was 96,769,103 bp in length and consisted of 14,447 scaffolds, and the Phyllactinia moricola assembly was 180,382,532 bp in length on 45,569 scaffolds. Together with the draft genome of L. taurica, these resources will be pivotal for understanding the molecular basis of the lifestyle of these fungi, which is unique within the family Erysiphaceae.

Host Adaptation Through Hybridization: Genome Analysis of Triticale Powdery Mildew Reveals Unique Combination of Lineage-Specific Effectors

The emergence of new fungal pathogens through hybridization represents a serious challenge for agriculture. Hybridization between the wheat mildew (Blumeria graminis f. sp. tritici) and rye mildew (B. graminis f. sp. secalis) pathogens has led to the emergence of a new mildew form (B. graminis f. sp. triticale) growing on triticale, a man-made amphiploid crop derived from crossing rye and wheat, which was originally resistant to the powdery mildew disease. The identification of the genetic basis of host adaptation in triticale mildew has been hampered by the lack of a reference genome. Here, we report the 141.4-Mb reference assembly of triticale mildew isolate THUN-12 derived from long-read sequencing and genetic map-based scaffolding. All 11 triticale mildew chromosomes were assembled from telomere-to-telomere and revealed that 19.7% of the hybrid genome was inherited from the rye mildew parental lineage. We identified lineage-specific regions in the hybrid, inherited from the rye or wheat mildew parental lineages, that harbor numerous bona fide candidate effectors. We propose that the combination of lineage-specific effectors in the hybrid genome is crucial for host adaptation, allowing the fungus to simultaneously circumvent the immune systems contributed by wheat and rye in the triticale crop. In line with this, we demonstrate the functional transfer of the SvrPm3 effector from wheat to triticale mildew, a virulence effector that specifically suppresses resistance of the wheat Pm3 allelic series. This transfer is the likely underlying cause for the observed poor effectiveness of several Pm3 alleles against triticale mildew and exemplifies the negative implications of pathogen hybridizations on resistance breeding.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

An automated and combinative method for the predictive ranking of candidate effector proteins of fungal plant pathogens

Fungal plant-pathogens promote infection of their hosts through the release of 'effectors'-a broad class of cytotoxic or virulence-promoting molecules. Effectors may be recognised by resistance or sensitivity receptors in the host, which can determine disease outcomes. Accurate prediction of effectors remains a major challenge in plant pathology, but if achieved will facilitate rapid improvements to host disease resistance. This study presents a novel tool and pipeline for the ranking of predicted effector candidates-Predector-which interfaces with multiple software tools and methods, aggregates disparate features that are relevant to fungal effector proteins, and applies a pairwise learning to rank approach. Predector outperformed a typical combination of secretion and effector prediction methods in terms of ranking performance when applied to a curated set of confirmed effectors derived from multiple species. We present Predector ( https://github.com/ccdmb/predector ) as a useful tool for the ranking of predicted effector candidates, which also aggregates and reports additional supporting information relevant to effector and secretome prediction in a simple, efficient, and reproducible manner.

The Powdery Mildew Effector CSEP0027 Interacts With Barley Catalase to Regulate Host Immunity

Powdery mildew is one of the most important fungal pathogen diseases. The genome of barley mildew fungus, Blumeria graminis f. sp. hordei (Bgh), encodes a large number of candidate secreted effector proteins (CSEPs). So far, the function and mechanism of most CSEPs remain largely unknown. Here, we identify a Bgh effector CSEP0027, a member of family 41, triggering cell death in Nicotiana benthamiana. CSEP0027 contains a functional signal peptide (SP), verified by yeast secretion assay. We show that CSEP0027 promotes Bgh virulence in barley infection using transient gene expression and host-induced gene silencing (HIGS). Barley catalase HvCAT1 is identified as a CSEP0027 interactor by yeast two-hybrid (Y2H) screening, and the interaction is verified in yeast, in vitro and in vivo. The coexpression of CSEP0027 and HvCAT1 in barley cells results in altered localization of HvCAT1 from the peroxisome to the nucleus. Barley stripe mosaic virus (BSMV)-silencing and transiently-induced gene silencing (TIGS) assays reveal that HvCAT1 is required for barley immunity against Bgh. We propose that CSEP0027 interacts with barley HvCAT1 to regulate the host immunity and likely reactive oxygen species (ROS) homeostasis to promote fungal virulence during barley infection.

A devil's bargain with transposable elements in plant pathogens

Transposable elements (TEs) spread in genomes through self-copying mechanisms and are a major cause of genome expansions. Plant pathogens have finely tuned the expression of virulence factors to rely on epigenetic control targeted at nearby TEs. Stress experienced during the plant infection process leads to derepression of TEs and concurrently allows the expression of virulence factors. We argue that the derepression of TEs elements causes an evolutionary conflict by favoring TEs that can be reactivated. Active TEs and recent genome size expansions indicate that plant pathogens could face long-term consequences from the short-term benefit of fine-tuning the infection process. Hence, encoding key virulence factors close to TEs under epigenetic control constitutes a devil's bargain for pathogens.

Powdery Mildew Resistance Phenotypes of Wheat Gene Bank Accessions

Powdery mildew (Blumeria graminis f. sp. tritici) is a common pathogen of bread wheat (Triticum aestivum L.), and genetic resistance is an effective and environmentally friendly method to reduce its adverse impact. The introgression of novel genes from wheat progenitors and related species can increase the diversity of disease resistance and accumulation of minor genes to improve the crop's resistance durability. To accomplish these two actions, host genotypes without major resistances should be preferably used. Therefore, the main aim of this study was to carry out seedling tests to detect such resistances in a set of wheat accessions from the Czech gene bank and to group the cultivars according to their phenotype. Ear progenies of 448 selected cultivars originating from 33 countries were inoculated with three isolates of the pathogen. Twenty-eight cultivars were heterogeneous, and 110 cultivars showed resistance to at least one isolate. Fifty-nine cultivars, mostly from Northwest Europe, were resistant to all three isolates were more than three times more frequently recorded in spring than in winter cultivars. Results will facilitate a rational and practical approach preferably using the set of cultivars without major resistances for both mentioned methods of breeding wheat cultivars resistant to powdery mildew.

Tactics of host manipulation by intracellular effectors from plant pathogenic fungi

Fungal pathogens can secrete hundreds of effectors, some of which are known to promote host susceptibility. This biological complexity, together with the lack of genetic tools in some fungi, presents a substantial challenge to develop a broad picture of the mechanisms these pathogens use for host manipulation. Nevertheless, recent advances in understanding individual effector functions are beginning to flesh out our view of fungal pathogenesis. This review discusses some of the latest findings that illustrate how effectors from diverse species use similar strategies to modulate plant physiology to their advantage. We also summarize recent breakthroughs in the identification of effectors from challenging systems, like obligate biotrophs, and emerging concepts such as the 'iceberg model' to explain how the activation of plant immunity can be turned off by effectors with suppressive activity.

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.

A Hybrid Genome Assembly Resource for Podosphaera xanthii, the Main Causal Agent of Powdery Mildew Disease in Cucurbits

Podosphaera xanthii is the main causal agent of powdery mildew in cucurbits and, arguably, the most important fungal pathogen of cucurbit crops. Here, we present the first reference genome assembly for P. xanthii. We performed a hybrid genome assembly, using reads from Illumina NextSeq550 and PacBio Sequel S3. The short and long reads were assembled into 1,727 scaffolds with an N₅₀ size of 163,173 bp, resulting in a 142-Mb genome size. The combination of homology-based and ab initio predictions allowed the prediction of 14,911 complete genes. Repetitive sequences comprised 76.2% of the genome. Our P. xanthii genome assembly improves considerably the molecular resources for research on P. xanthii-cucurbit interactions and provides new opportunities for further genomics, transcriptomics, and evolutionary studies in powdery mildew fungi.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A highly differentiated region of wheat chromosome 7AL encodes a Pm1a immune receptor that recognizes its corresponding AvrPm1a effector from Blumeria graminis

Pm1a, the first powdery mildew resistance gene described in wheat, is part of a complex resistance (R) gene cluster located in a distal region of chromosome 7AL that has suppressed genetic recombination. A nucleotide-binding, leucine-rich repeat (NLR) immune receptor gene was isolated using mutagenesis and R gene enrichment sequencing (MutRenSeq). Stable transformation confirmed Pm1a identity which induced a strong resistance phenotype in transgenic plants upon challenge with avirulent Blumeria graminis (wheat powdery mildew) pathogens. A high-density genetic map of a B. graminis family segregating for Pm1a avirulence combined with pathogen genome resequencing and RNA sequencing (RNAseq) identified AvrPm1a effector gene candidates. In planta expression identified an effector, with an N terminal Y/FxC motif, that induced a strong hypersensitive response when co-expressed with Pm1a in Nicotiana benthamiana. Single chromosome enrichment sequencing (ChromSeq) and assembly of chromosome 7A suggested that suppressed recombination around the Pm1a region was due to a rearrangement involving chromosomes 7A, 7B and 7D. The cloning of Pm1a and its identification in a highly rearranged region of chromosome 7A provides insight into the role of chromosomal rearrangements in the evolution of this complex resistance cluster.

The leucine-rich repeats in allelic barley MLA immune receptors define specificity towards sequence-unrelated powdery mildew avirulence effectors with a predicted common RNase-like fold

Nucleotide-binding domain leucine-rich repeat-containing receptors (NLRs) in plants can detect avirulence (AVR) effectors of pathogenic microbes. The Mildew locus a (Mla) NLR gene has been shown to confer resistance against diverse fungal pathogens in cereal crops. In barley, Mla has undergone allelic diversification in the host population and confers isolate-specific immunity against the powdery mildew-causing fungal pathogen Blumeria graminis forma specialis hordei (Bgh). We previously isolated the Bgh effectors AVRA1, AVRA7, AVRA9, AVRA13, and allelic AVRA10/AVRA22, which are recognized by matching MLA1, MLA7, MLA9, MLA13, MLA10 and MLA22, respectively. Here, we extend our knowledge of the Bgh effector repertoire by isolating the AVRA6 effector, which belongs to the family of catalytically inactive RNase-Like Proteins expressed in Haustoria (RALPHs). Using structural prediction, we also identified RNase-like folds in AVRA1, AVRA7, AVRA10/AVRA22, and AVRA13, suggesting that allelic MLA recognition specificities could detect structurally related avirulence effectors. To better understand the mechanism underlying the recognition of effectors by MLAs, we deployed chimeric MLA1 and MLA6, as well as chimeric MLA10 and MLA22 receptors in plant co-expression assays, which showed that the recognition specificity for AVRA1 and AVRA6 as well as allelic AVRA10 and AVRA22 is largely determined by the receptors' C-terminal leucine-rich repeats (LRRs). The design of avirulence effector hybrids allowed us to identify four specific AVRA10 and five specific AVRA22 aa residues that are necessary to confer MLA10- and MLA22-specific recognition, respectively. This suggests that the MLA LRR mediates isolate-specific recognition of structurally related AVRA effectors. Thus, functional diversification of multi-allelic MLA receptors may be driven by a common structural effector scaffold, which could be facilitated by proliferation of the RALPH effector family in the pathogen genome.

Genome-Wide Analyses of Repeat-Induced Point Mutations in the Ascomycota

The Repeat-Induced Point (RIP) mutation pathway is a fungus-specific genome defense mechanism that mitigates the deleterious consequences of repeated genomic regions and transposable elements (TEs). RIP mutates targeted sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence and extent of RIP with a sliding-window approach. Using genome-wide RIP data and two sets of control groups, the association between RIP, TEs, and GC content were contrasted in organisms capable and incapable of RIP. Based on these data, we then set out to determine the extent and occurrence of RIP in 58 representatives of the Ascomycota. The findings were summarized by placing each of the fungi investigated in one of six categories based on the extent of genome-wide RIP. In silico RIP analyses, using a sliding-window approach with stringent RIP parameters, implemented simultaneously within the same genetic context, on high quality genome assemblies, yielded superior results in determining the genome-wide RIP among the Ascomycota. Most Ascomycota had RIP and these mutations were particularly widespread among classes of the Pezizomycotina, including the early diverging Orbiliomycetes and the Pezizomycetes. The most extreme cases of RIP were limited to representatives of the Dothideomycetes and Sordariomycetes. By contrast, the genomes of the Taphrinomycotina and Saccharomycotina contained no detectable evidence of RIP. Also, recent losses in RIP combined with controlled TE proliferation in the Pezizomycotina subphyla may promote substantial genome enlargement as well as the formation of sub-genomic compartments. These findings have broadened our understanding of the taxonomic range and extent of RIP in Ascomycota and how this pathway affects the genomes of fungi harboring it.

Comparative Transcriptome Analysis of Wheat Lines in the Field Reveals Multiple Essential Biochemical Pathways Suppressed by Obligate Pathogens

Mildew and rust are the most devastating cereal pathogens, and in wheat they can cause up to 50% yield loss every year. Wheat lines containing resistance genes are used to effectively control fungal diseases, but the molecular mechanisms underlying the interaction between wheat and its fungal pathogens are poorly understood. Here, we used RNA sequencing (RNA-Seq) to compare the transcriptomic landscape of susceptible and resistant wheat lines to identify genes and pathways that are targeted by obligate biotrophic fungal pathogens. The five lines differed in the expression of thousands of genes under infection as well as control conditions. Generally, mixed infection with powdery mildew and leaf rust resulted in downregulation of numerous genes in susceptible lines. Interestingly, transcriptomic comparison between the nearly isogenic lines Thatcher and Thatcher-Lr34 identified 753 genes that are uniquely downregulated in the susceptible line upon infection. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis, revealed the suppression of six major biochemical pathways, namely nuclear transport, alternative splicing, DNA damage response, ubiquitin-mediated proteolysis, phosphoinositol signaling, and photosynthesis. We conclude that powdery mildew and leaf rust evade the wheat defense system by suppression of programmed cell death (PCD) and responses to cellular damage. Considering the broad range of the induced changes, we propose that the pathogen targets "master regulators" at critical steps in the respective pathways. Identification of these wheat genes targeted by the pathogen could inspire new directions for future wheat breeding.

Effector Biology of Biotrophic Plant Fungal Pathogens: Current Advances and Future Prospects

The interaction of fungal pathogens with their host requires a novel invading mechanism and the presence of various virulence-associated components responsible for promoting the infection. The small secretory proteins, explicitly known as effector proteins, are one of the prime mechanisms of host manipulation utilized by the pathogen to disarm the host. Several effector proteins are known to translocate from fungus to the plant cell for host manipulation. Many fungal effectors have been identified using genomic, transcriptomic, and bioinformatics approaches. Most of the effector proteins are devoid of any conserved signatures, and their prediction based on sequence homology is very challenging, therefore by combining the sequence consensus based upon machine learning features, multiple tools have also been developed for predicting apoplastic and cytoplasmic effectors. Various post-genomics approaches like transcriptomics of virulent isolates have also been utilized for identifying active consortia of effectors. Significant progress has been made in understanding biotrophic effectors; however, most of it is underway due to their complex interaction with host and complicated recognition and signaling networks. This review discusses advances, and challenges in effector identification and highlighted various features of the potential effector proteins and approaches for understanding their genetics and strategies for regulation.

Ultraviolet Mutagenesis Coupled with Next-Generation Sequencing as a Method for Functional Interrogation of Powdery Mildew Genomes

Powdery mildews are obligate biotrophic fungal pathogens causing important diseases of plants worldwide. Very little is known about the requirements for their pathogenicity at the molecular level. This is largely due to the inability to culture these organisms in vitro or to modify them genetically. Here, we describe a mutagenesis procedure based on ultraviolet (UV) irradiation to accumulate mutations in the haploid genome of the barley powdery mildew pathogen Blumeria graminis f. sp. hordei. Exposure of B. graminis f. sp. hordei conidia to different durations of UV-C radiation (10 s to 12 min) resulted in a reduced number of macroscopically visible fungal colonies. B. graminis f. sp. hordei colony number was negatively correlated with exposure time and the total number of consecutive cycles of UV irradiation. Dark incubation following UV exposure further reduced fungal viability, implying that photoreactivation is an important component of DNA repair in B. graminis f. sp. hordei. After several rounds of UV mutagenesis, we selected two mutant isolates in addition to the parental B. graminis f. sp. hordei K1 isolate for whole-genome resequencing. By combining automated prediction of sequence variants and their manual validation, we identified unique UV-induced mutations in the genomes of the two isolates. Most of these mutations were in the up- or downstream regions of genes or in the intergenic space. Some of the variants detected in genes led to predicted missense mutations. As an additional insight, our bioinformatic analyses revealed a complex population structure within supposedly clonal B. graminis f. sp. hordei isolates.

The Many Questions about Mini Chromosomes in Colletotrichum spp

Many fungal pathogens carry accessory regions in their genome, which are not required for vegetative fitness. Often, although not always, these regions occur as relatively small chromosomes in different species. Such mini chromosomes appear to be a typical feature of many filamentous plant pathogens. Since these regions often carry genes coding for effectors or toxin-producing enzymes, they may be directly related to virulence of the respective pathogen. In this review, we outline the situation of small accessory chromosomes in the genus Colletotrichum, which accounts for ecologically important plant diseases. We summarize which species carry accessory chromosomes, their gene content, and chromosomal makeup. We discuss the large variation in size and number even between different isolates of the same species, their potential roles in host range, and possible mechanisms for intra- and interspecies exchange of these interesting genetic elements.