Coevolution between a family of parasite virulence effectors and a class of LINE-1 retrotransposons

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.

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.

The Venturia inaequalis effector repertoire is dominated by expanded families with predicted structural similarity, but unrelated sequence, to avirulence proteins from other plant-pathogenic fungi

Background

Scab, caused by the biotrophic fungus Venturia inaequalis, is the most economically important disease of apples worldwide. During infection, V. inaequalis occupies the subcuticular environment, where it secretes virulence factors, termed effectors, to promote host colonization. Consistent with other plant-pathogenic fungi, many of these effectors are expected to be non-enzymatic proteins, some of which can be recognized by corresponding host resistance proteins to activate plant defences, thus acting as avirulence determinants. To develop durable control strategies against scab, a better understanding of the roles that these effector proteins play in promoting subcuticular growth by V. inaequalis, as well as in activating, suppressing, or circumventing resistance protein-mediated defences in apple, is required.

Results

We generated the first comprehensive RNA-seq transcriptome of V. inaequalis during colonization of apple. Analysis of this transcriptome revealed five temporal waves of gene expression that peaked during early, mid, or mid-late infection. While the number of genes encoding secreted, non-enzymatic proteinaceous effector candidates (ECs) varied in each wave, most belonged to waves that peaked in expression during mid-late infection. Spectral clustering based on sequence similarity determined that the majority of ECs belonged to expanded protein families. To gain insights into function, the tertiary structures of ECs were predicted using AlphaFold2. Strikingly, despite an absence of sequence similarity, many ECs were predicted to have structural similarity to avirulence proteins from other plant-pathogenic fungi, including members of the MAX, LARS, ToxA and FOLD effector families. In addition, several other ECs, including an EC family with sequence similarity to the AvrLm6 avirulence effector from Leptosphaeria maculans, were predicted to adopt a KP6-like fold. Thus, proteins with a KP6-like fold represent another structural family of effectors shared among plant-pathogenic fungi.

Conclusions

Our study reveals the transcriptomic profile underpinning subcuticular growth by V. inaequalis and provides an enriched list of ECs that can be investigated for roles in virulence and avirulence. Furthermore, our study supports the idea that numerous sequence-unrelated effectors across plant-pathogenic fungi share common structural folds. In doing so, our study gives weight to the hypothesis that many fungal effectors evolved from ancestral genes through duplication, followed by sequence diversification, to produce sequence-unrelated but structurally similar proteins.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01442-9.

Transcriptional activity and epigenetic regulation of transposable elements in the symbiotic fungus Rhizophagus irregularis

Arbuscular mycorrhizal (AM) fungi form mutualistic relationships with most land plant species. AM fungi have long been considered as ancient asexuals. Long-term clonal evolution would be remarkable for a eukaryotic lineage and suggests the importance of alternative mechanisms to promote genetic variability facilitating adaptation. Here, we assessed the potential of transposable elements for generating such genomic diversity. The dynamic expression of TEs during Rhizophagus irregularis spore development suggests ongoing TE activity. We find Mutator-like elements located near genes belonging to highly expanded gene families. Whole-genome epigenomic profiling of R. irregularis provides direct evidence of DNA methylation and small RNA production occurring at TE loci. Our results support a model in which TE activity shapes the genome, while DNA methylation and small RNA-mediated silencing keep their overproliferation in check. We propose that a well-controlled TE activity directly contributes to genome evolution in AM fungi.

Variable genome evolution in fungi after transposon-mediated amplification of a housekeeping gene

Background

Transposable elements (TEs) can be key drivers of evolution, but the mechanisms and scope of how they impact gene and genome function are largely unknown. Previous analyses revealed that TE-mediated gene amplifications can have variable effects on fungal genomes, from inactivation of function to production of multiple active copies. For example, a DNA methyltransferase gene in the wheat pathogen Zymoseptoria tritici (synonym Mycosphaerella graminicola) was amplified to tens of copies, all of which were inactivated by Repeat-Induced Point mutation (RIP) including the original, resulting in loss of cytosine methylation. In another wheat pathogen, Pyrenophora tritici-repentis, a histone H3 gene was amplified to tens of copies with little evidence of RIP, leading to many potentially active copies. To further test the effects of transposon-aided gene amplifications on genome evolution and architecture, the repetitive fraction of the significantly expanded genome of the banana pathogen, Pseudocercospora fijiensis, was analyzed in greater detail.

Results

These analyses identified a housekeeping gene, histone H3, which was captured and amplified to hundreds of copies by a hAT DNA transposon, all of which were inactivated by RIP, except for the original. In P. fijiensis the original H3 gene probably was not protected from RIP, but most likely was maintained intact due to strong purifying selection. Comparative analyses revealed that a similar event occurred in five additional genomes representing the fungal genera Cercospora, Pseudocercospora and Sphaerulina.

Conclusions

These results indicate that the interplay of TEs and RIP can result in different and unpredictable fates of amplified genes, with variable effects on gene and genome evolution.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0177-0) contains supplementary material, which is available to authorized users.

A barley powdery mildew fungus non-autonomous retrotransposon encodes a peptide that supports penetration success on barley

Pathogens overcome plant immunity by means of secreted effectors. Host effector targets often act in pathogen defense, but might also support fungal accommodation or nutrition. The barley ROP GTPase HvRACB is involved in accommodation of fungal haustoria of the powdery mildew fungus Blumeria graminis f.sp. hordei (Bgh) in barley epidermal cells. We found that HvRACB interacts with the ROP-interactive peptide 1 (ROPIP1) that is encoded on the active non-long terminal repeat retroelement Eg-R1 of Bgh. Overexpression of ROPIP1 in barley epidermal cells and host-induced post-transcriptional gene silencing (HIGS) of ROPIP1 suggested that ROPIP1 is involved in virulence of Bgh. Bimolecular fluorescence complementation and co-localization supported that ROPIP1 can interact with activated HvRACB in planta. We show that ROPIP1 is expressed by Bgh on barley and translocated into the cytoplasm of infected barley cells. ROPIP1 is recruited to microtubules upon co-expression of MICROTUBULE ASSOCIATED ROP GTPase ACTIVATING PROTEIN (HvMAGAP1) and can destabilize cortical microtubules. The data suggest that Bgh ROPIP targets HvRACB and manipulates host cell microtubule organization for facilitated host cell entry. This points to a possible neo-functionalization of retroelement-derived transcripts for the evolution of a pathogen virulence effector.

Novel jack-in-the-box effector of the barley powdery mildew pathogen?

This article comments on: Nottensteiner M, Zechmann B, McCollum C, Hückelhoven R. 2018. A barley powdery mildew fungus non-autonomous retrotransposon encodes a peptide that supports penetration success on barley. Journal of Experimental Botany 69, 3745–3758.

Comparative analysis of the predicted secretomes of Rosaceae scab pathogens Venturia inaequalis and V. pirina reveals expanded effector families and putative determinants of host range

BACKGROUND

Fungal plant pathogens belonging to the genus Venturia cause damaging scab diseases of members of the Rosaceae. In terms of economic impact, the most important of these are V. inaequalis, which infects apple, and V. pirina, which is a pathogen of European pear. Given that Venturia fungi colonise the sub-cuticular space without penetrating plant cells, it is assumed that effectors that contribute to virulence and determination of host range will be secreted into this plant-pathogen interface. Thus the predicted secretomes of a range of isolates of Venturia with distinct host-ranges were interrogated to reveal putative proteins involved in virulence and pathogenicity.

RESULTS

Genomes of Venturia pirina (one European pear scab isolate) and Venturia inaequalis (three apple scab, and one loquat scab, isolates) were sequenced and the predicted secretomes of each isolate identified. RNA-Seq was conducted on the apple-specific V. inaequalis isolate Vi1 (in vitro and infected apple leaves) to highlight virulence and pathogenicity components of the secretome. Genes encoding over 600 small secreted proteins (candidate effectors) were identified, most of which are novel to Venturia, with expansion of putative effector families a feature of the genus. Numerous genes with similarity to Leptosphaeria maculans AvrLm6 and the Verticillium spp. Ave1 were identified. Candidates for avirulence effectors with cognate resistance genes involved in race-cultivar specificity were identified, as were putative proteins involved in host-species determination. Candidate effectors were found, on average, to be in regions of relatively low gene-density and in closer proximity to repeats (e.g. transposable elements), compared with core eukaryotic genes.

CONCLUSIONS

Comparative secretomics has revealed candidate effectors from Venturia fungal plant pathogens that attack pome fruit. Effectors that are putative determinants of host range were identified; both those that may be involved in race-cultivar and host-species specificity. Since many of the effector candidates are in close proximity to repetitive sequences this may point to a possible mechanism for the effector gene family expansion observed and a route to diversification via transposition and repeat-induced point mutation.

De novo Analysis of the Epiphytic Transcriptome of the Cucurbit Powdery Mildew Fungus Podosphaera xanthii and Identification of Candidate Secreted Effector Proteins

The cucurbit powdery mildew fungus Podosphaera xanthii is a major limiting factor for cucurbit production worldwide. Despite the fungus’s agronomic and economic importance, very little is known about fundamental aspects of P. xanthii biology, such as obligate biotrophy or pathogenesis. To design more durable control strategies, genomic information about P. xanthii is needed. Powdery mildews are fungal pathogens with large genomes compared with those of other fungi, which contain vast amounts of repetitive DNA sequences, much of which is composed of retrotransposons. To reduce genome complexity, in this work we aimed to obtain and analyse the epiphytic transcriptome of P. xanthii as a starting point for genomic research. Total RNA was isolated from epiphytic fungal material, and the corresponding cDNA library was sequenced using a 454 GS FLX platform. Over 676,562 reads were obtained and assembled into 37,241 contigs. Annotation data identified 8,798 putative genes with different orthologues. As described for other powdery mildew fungi, a similar set of missing core ascomycete genes was found, which may explain obligate biotrophy. To gain insight into the plant-pathogen relationships, special attention was focused on the analysis of the secretome. After this analysis, 137 putative secreted proteins were identified, including 53 candidate secreted effector proteins (CSEPs). Consistent with a putative role in pathogenesis, the expression profile observed for some of these CSEPs showed expression maxima at the beginning of the infection process at 24 h after inoculation, when the primary appressoria are mostly formed. Our data mark the onset of genomics research into this very important pathogen of cucurbits and shed some light on the intimate relationship between this pathogen and its host plant.

Avirulence Genes in Cereal Powdery Mildews: The Gene-for-Gene Hypothesis 2.0

The gene-for-gene hypothesis states that for each gene controlling resistance in the host, there is a corresponding, specific gene controlling avirulence in the pathogen. Allelic series of the cereal mildew resistance genes Pm3 and Mla provide an excellent system for genetic and molecular analysis of resistance specificity. Despite this opportunity for molecular research, avirulence genes in mildews remain underexplored. Earlier work in barley powdery mildew (B.g. hordei) has shown that the reaction to some Mla resistance alleles is controlled by multiple genes. Similarly, several genes are involved in the specific interaction of wheat mildew (B.g. tritici) with the Pm3 allelic series. We found that two mildew genes control avirulence on Pm3f: one gene is involved in recognition by the resistance protein as demonstrated by functional studies in wheat and the heterologous host Nicotiana benthamiana. A second gene is a suppressor, and resistance is only observed in mildew genotypes combining the inactive suppressor and the recognized Avr. We propose that such suppressor/avirulence gene combinations provide the basis of specificity in mildews. Depending on the particular gene combinations in a mildew race, different genes will be genetically identified as the "avirulence" gene. Additionally, the observation of two LINE retrotransposon-encoded avirulence genes in B.g. hordei further suggests that the control of avirulence in mildew is more complex than a canonical gene-for-gene interaction. To fully understand the mildew-cereal interactions, more knowledge on avirulence determinants is needed and we propose ways how this can be achieved based on recent advances in the field.

A Large Family of AvrLm6-like Genes in the Apple and Pear Scab Pathogens, Venturia inaequalis and Venturia pirina

Venturia inaequalis and V. pirina are Dothideomycete fungi that cause apple scab and pear scab disease, respectively. Whole genome sequencing of V. inaequalis and V. pirina isolates has revealed predicted proteins with sequence similarity to AvrLm6, a Leptosphaeria maculans effector that triggers a resistance response in Brassica napus and B. juncea carrying the resistance gene, Rlm6. AvrLm6-like genes are present as large families (>15 members) in all sequenced strains of V. inaequalis and V. pirina, while in L. maculans, only AvrLm6 and a single paralog have been identified. The Venturia AvrLm6-like genes are located in gene-poor regions of the genomes, and mostly in close proximity to transposable elements, which may explain the expansion of these gene families. An AvrLm6-like gene from V. inaequalis with the highest sequence identity to AvrLm6 was unable to trigger a resistance response in Rlm6-carrying B. juncea. RNA-seq and qRT-PCR gene expression analyses, of in planta- and in vitro-grown V. inaequalis, has revealed that many of the AvrLm6-like genes are expressed during infection. An AvrLm6 homolog from V. inaequalis that is up-regulated during infection was shown (using an eYFP-fusion protein construct) to be localized to the sub-cuticular stroma during biotrophic infection of apple hypocotyls.

Evolution of the EKA family of powdery mildew avirulence-effector genes from the ORF 1 of a LINE retrotransposon

Background

The Avrk1 and Avra10 avirulence (AVR) genes encode effectors that increase the pathogenicity of the fungus Blumeria graminis f.sp. hordei (Bgh), the powdery mildew pathogen, in susceptible barley plants. In resistant barley, MLK1 and MLA10 resistance proteins recognize the presence of AVRK1 and AVRA10, eliciting the hypersensitive response typical of gene for gene interactions. Avrk1 and Avra10 have more than 1350 homologues in Bgh genome, forming the EKA (Effectors homologous to Avrk1 and Avra10) gene family.

Results

We tested the hypothesis that the EKA family originated from degenerate copies of Class I LINE retrotransposons by analysing the EKA family in the genome of Bgh isolate DH14 with bioinformatic tools specially developed for the analysis of Transposable Elements (TE) in genomes. The Class I LINE retrotransposon copies homologous to Avrk1 and Avra10 represent 6.5 % of the Bgh annotated genome and, among them, we identified 293 AVR/effector candidate genes. We also experimentally identified peptides that indicated the translation of several predicted proteins from EKA family members, which had higher relative abundance in haustoria than in hyphae.

Conclusions

Our analyses indicate that Avrk1 and Avra10 have evolved from part of the ORF1 gene of Class I LINE retrotransposons. The co-option of Avra10 and Avrk1 as effectors from truncated copies of retrotransposons explains the huge number of homologues in Bgh genome that could act as dynamic reservoirs from which new effector genes may evolve. These data provide further evidence for recruitment of retrotransposons in the evolution of new biological functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2185-x) contains supplementary material, which is available to authorized users.

The landscape of transposable elements in the finished genome of the fungal wheat pathogen Mycosphaerella graminicola

Background

In addition to gene identification and annotation, repetitive sequence analysis has become an integral part of genome sequencing projects. Identification of repeats is important not only because it improves gene prediction, but also because of the role that repetitive sequences play in determining the structure and evolution of genes and genomes. Several methods using different repeat-finding strategies are available for whole-genome repeat sequence analysis. Four independent approaches were used to identify and characterize the repetitive fraction of the Mycosphaerella graminicola (synonym Zymoseptoria tritici) genome. This ascomycete fungus is a wheat pathogen and its finished genome comprises 21 chromosomes, eight of which can be lost with no obvious effects on fitness so are dispensable.

Results

Using a combination of four repeat-finding methods, at least 17% of the M. graminicola genome was estimated to be repetitive. Class I transposable elements, that amplify via an RNA intermediate, account for about 70% of the total repetitive content in the M. graminicola genome. The dispensable chromosomes had a higher percentage of repetitive elements as compared to the core chromosomes. Distribution of repeats across the chromosomes also varied, with at least six chromosomes showing a non-random distribution of repetitive elements. Repeat families showed transition mutations and a CpA → TpA dinucleotide bias, indicating the presence of a repeat-induced point mutation (RIP)-like mechanism in M. graminicola. One gene family and two repeat families specific to subtelomeres also were identified in the M. graminicola genome. A total of 78 putative clusters of nested elements was found in the M. graminicola genome. Several genes with putative roles in pathogenicity were found associated with these nested repeat clusters. This analysis of the transposable element content in the finished M. graminicola genome resulted in a thorough and highly curated database of repetitive sequences.

Conclusions

This comprehensive analysis will serve as a scaffold to address additional biological questions regarding the origin and fate of transposable elements in fungi. Future analyses of the distribution of repetitive sequences in M. graminicola also will be able to provide insights into the association of repeats with genes and their potential role in gene and genome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1132) contains supplementary material, which is available to authorized users.

Adaptive genomic structural variation in the grape powdery mildew pathogen, Erysiphe necator

BACKGROUND

Powdery mildew, caused by the obligate biotrophic fungus Erysiphe necator, is an economically important disease of grapevines worldwide. Large quantities of fungicides are used for its control, accelerating the incidence of fungicide-resistance. Copy number variations (CNVs) are unbalanced changes in the structure of the genome that have been associated with complex traits. In addition to providing the first description of the large and highly repetitive genome of E. necator, this study describes the impact of genomic structural variation on fungicide resistance in Erysiphe necator.

RESULTS

A shotgun approach was applied to sequence and assemble the genome of five E. necator isolates, and RNA-seq and comparative genomics were used to predict and annotate protein-coding genes. Our results show that the E. necator genome is exceptionally large and repetitive and suggest that transposable elements are responsible for genome expansion. Frequent structural variations were found between isolates and included copy number variation in EnCYP51, the target of the commonly used sterol demethylase inhibitor (DMI) fungicides. A panel of 89 additional E. necator isolates collected from diverse vineyard sites was screened for copy number variation in the EnCYP51 gene and for presence/absence of a point mutation (Y136F) known to result in higher fungicide tolerance. We show that an increase in EnCYP51 copy number is significantly more likely to be detected in isolates collected from fungicide-treated vineyards. Increased EnCYP51 copy numbers were detected with the Y136F allele, suggesting that an increase in copy number becomes advantageous only after the fungicide-tolerant allele is acquired. We also show that EnCYP51 copy number influences expression in a gene-dose dependent manner and correlates with fungal growth in the presence of a DMI fungicide.

CONCLUSIONS

Taken together our results show that CNV can be adaptive in the development of resistance to fungicides by providing increasing quantitative protection in a gene-dosage dependent manner. The results of this work not only demonstrate the effectiveness of using genomics to dissect complex traits in organisms with very limited molecular information, but also may have broader implications for understanding genomic dynamics in response to strong selective pressure in other pathogens with similar genome architectures.

In silico analysis of the core signaling proteome from the barley powdery mildew pathogen (Blumeria graminis f.sp. hordei)

BACKGROUND

Compared to other ascomycetes, the barley powdery mildew pathogen Blumeria graminis f.sp. hordei (Bgh) has a large genome (ca. 120 Mbp) that harbors a relatively small number of protein-coding genes (ca. 6500). This genomic assemblage is thought to be the result of numerous gene losses, which likely represent an evolutionary adaptation to a parasitic lifestyle in close association with its host plant, barley (Hordeum vulgare). Approximately 8% of the Bgh genes are predicted to encode virulence effectors that are secreted into host tissue and/or cells to promote pathogenesis; the remaining proteome is largely uncharacterized at present.

RESULTS

We provide a comparative analysis of the conceptual Bgh proteome, with an emphasis on proteins with known roles in fungal development and pathogenicity, for example heterotrimeric G proteins and G protein coupled receptors; components of calcium and cAMP signaling; small monomeric GTPases; mitogen-activated protein cascades and transcription factors. The predicted Bgh proteome lacks a number of proteins that are otherwise conserved in filamentous fungi, including two proteins that are required for the formation of anastomoses (somatic hyphal connections). By contrast, apart from minor modifications, all major canonical signaling pathways are retained in Bgh. A family of kinases that preferentially occur in pathogenic species of the fungal clade Leotiomyceta is unusually expanded in Bgh and its close relative, Blumeria graminis f.sp. tritici.

CONCLUSIONS

Our analysis reveals characteristic features of the proteome of a fungal phytopathogen that occupies an extreme habitat: the living plant cell.

Interaction of a Blumeria graminis f. sp. hordei effector candidate with a barley ARF-GAP suggests that host vesicle trafficking is a fungal pathogenicity target

Filamentous phytopathogens, such as fungi and oomycetes, secrete effector proteins to establish successful interactions with their plant hosts. In contrast with oomycetes, little is known about effector functions in true fungi. We used a bioinformatics pipeline to identify Blumeria effector candidates (BECs) from the obligate biotrophic barley powdery mildew pathogen, Blumeria graminis f. sp. hordei (Bgh). BEC1-BEC5 are expressed at different time points during barley infection. BEC1, BEC2 and BEC4 have orthologues in the Arabidopsis thaliana-infecting powdery mildew fungus Golovinomyces orontii. Arabidopsis lines stably expressing the G. orontii BEC2 orthologue, GoEC2, are more susceptible to infection with the non-adapted fungus Erysiphe pisi, suggesting that GoEC2 contributes to powdery mildew virulence. For BEC3 and BEC4, we identified thiopurine methyltransferase, a ubiquitin-conjugating enzyme, and an ADP ribosylation factor-GTPase-activating protein (ARF-GAP) as potential host targets. Arabidopsis knockout lines of the respective HvARF-GAP orthologue (AtAGD5) allowed higher entry levels of E. pisi, but exhibited elevated resistance to the oomycete Hyaloperonospora arabidopsidis. We hypothesize that ARF-GAP proteins are conserved targets of powdery and downy mildew effectors, and we speculate that BEC4 might interfere with defence-associated host vesicle trafficking.

An immunity-triggering effector from the Barley smut fungus Ustilago hordei resides in an Ustilaginaceae-specific cluster bearing signs of transposable element-assisted evolution

The basidiomycete smut fungus Ustilago hordei was previously shown to comprise isolates that are avirulent on various barley host cultivars. Through genetic crosses we had revealed that a dominant avirulence locus UhAvr1 which triggers immunity in barley cultivar Hannchen harboring resistance gene Ruh1, resided within an 80-kb region. DNA sequence analysis of this genetically delimited region uncovered the presence of 7 candidate secreted effector proteins. Sequence comparison of their coding sequences among virulent and avirulent parental and field isolates could not distinguish UhAvr1 candidates. Systematic deletion and complementation analyses revealed that UhAvr1 is UHOR_10022 which codes for a small effector protein of 171 amino acids with a predicted 19 amino acid signal peptide. Virulence in the parental isolate is caused by the insertion of a fragment of 5.5 kb with similarity to a common U. hordei transposable element (TE), interrupting the promoter of UhAvr1 and thereby changing expression and hence recognition of UhAVR1p. This rearrangement is likely caused by activities of TEs and variation is seen among isolates. Using GFP-chimeric constructs we show that UhAvr1 is induced only in mated dikaryotic hyphae upon sensing and infecting barley coleoptile cells. When infecting Hannchen, UhAVR1p causes local callose deposition and the production of reactive oxygen species and necrosis indicative of the immune response. UhAvr1 does not contribute significantly to overall virulence. UhAvr1 is located in a cluster of ten effectors with several paralogs and over 50% of TEs. This cluster is syntenous with clusters in closely-related U. maydis and Sporisorium reilianum. In these corn-infecting species, these clusters harbor however more and further diversified homologous effector families but very few TEs. This increased variability may have resulted from past selection pressure by resistance genes since U. maydis is not known to trigger immunity in its corn host.

Transposable element dynamics among asymbiotic and ectomycorrhizal Amanita fungi

Transposable elements (TEs) are ubiquitous inhabitants of eukaryotic genomes and their proliferation and dispersal shape genome architectures and diversity. Nevertheless, TE dynamics are often explored for one species at a time and are rarely considered in ecological contexts. Recent work with plant pathogens suggests a link between symbiosis and TE abundance. The genomes of pathogenic fungi appear to house an increased abundance of TEs, and TEs are frequently associated with the genes involved in symbiosis. To investigate whether this pattern is general, and relevant to mutualistic plant-fungal symbioses, we sequenced the genomes of related asymbiotic (AS) and ectomycorrhizal (ECM) Amanita fungi. Using methods developed to interrogate both assembled and unassembled sequences, we characterized and quantified TEs across three AS and three ECM species, including the AS outgroup Volvariella volvacea. The ECM genomes are characterized by abundant numbers of TEs, an especially prominent feature of unassembled sequencing libraries. Increased TE activity in ECM species is also supported by phylogenetic analysis of the three most abundant TE superfamilies; phylogenies revealed many radiations within contemporary ECM species. However, the AS species Amanita thiersii also houses extensive amplifications of elements, highlighting the influence of additional evolutionary parameters on TE abundance. Our analyses provide further evidence for a link between symbiotic associations among plants and fungi, and increased TE activity, while highlighting the importance individual species’ natural histories may have in shaping genome architecture.

Formae speciales of cereal powdery mildew: close or distant relatives?

Powdery mildew is an important disease of cereals, affecting both grain yield and end-use quality. The causal agent of powdery mildew on cereals, Blumeria graminis, has been classified into eight formae speciales (ff.spp.), infecting crops and wild grasses. Advances in research on host specificity and resistance, and on pathogen phylogeny and origins, have brought aspects of the subspecific classification system of B. graminis into ff.spp. into question, because it is based on adaptation to certain hosts rather than strict host specialization. Cereals therefore cannot be considered as typical non-hosts to non-adapted ff.spp. We introduce the term 'non-adapted resistance' of cereals to inappropriate ff.spp. of B. graminis, which involves both pathogen-associated molecular pattern-triggered immunity (PTI) and effector-triggered immunity (ETI). There is no clear distinction between the mechanisms of resistance to adapted and non-adapted ff.spp. Molecular evolutionary data suggest that the taxonomic grouping of B. graminis into different ff.spp. is not consistent with the phylogeny of the fungus. Imprecise estimates of mutation rates and the lack of genetic variation in introduced populations may explain the uncertainty with regard to divergence times, in the Miocene or Holocene epochs, of ff.spp. of B. graminis which infect cereal crop species. We propose that most evidence favours divergence in the Holocene, during the course of early agriculture. We also propose that the forma specialis concept should be retained for B. graminis pathogenic on cultivated cereals to include clades of the fungus which are strongly specialized to these hosts, i.e. ff.spp. hordei, secalis and tritici, as well as avenae from cultivated A. sativa, and that the forma specialis concept should no longer be applied to B. graminis from most wild grasses.

Small RNAs: a new paradigm in plant-microbe interactions

A never-ending arms race drives coevolution between pathogens and hosts. In plants, pathogen attacks invoke multiple layers of host immune responses. Many pathogens deliver effector proteins into host cells to suppress host immunity, and many plants have evolved resistance proteins to recognize effectors and trigger robust resistance. Here, we discuss findings on noncoding small RNAs (sRNAs) from plants and pathogens, which regulate host immunity and pathogen virulence. Recent discoveries have unveiled the role of noncoding sRNAs from eukaryotic pathogens and bacteria in pathogenicity in both plant and animal hosts. The discovery of fungal sRNAs that are delivered into host cells to suppress plant immunity added sRNAs to the list of pathogen effectors. Similar to protein effector genes, many of these sRNAs are generated from transposable element (TE) regions, which are likely to contribute to rapidly evolving virulence and host adaptation. We also discuss RNA silencing that occurs between organisms.

Epigenetic control of effectors in plant pathogens

Plant pathogens display impressive versatility in adapting to host immune systems. Pathogen effector proteins facilitate disease but can become avirulence (Avr) factors when the host acquires discrete recognition capabilities that trigger immunity. The mechanisms that lead to changes to pathogen Avr factors that enable escape from host immunity are diverse, and include epigenetic switches that allow for reuse or recycling of effectors. This perspective outlines possibilities of how epigenetic control of Avr effector gene expression may have arisen and persisted in filamentous plant pathogens, and how it presents special problems for diagnosis and detection of specific pathogen strains or pathotypes.

Mosaic genome structure of the barley powdery mildew pathogen and conservation of transcriptional programs in divergent hosts

Barley powdery mildew, Blumeria graminis f. sp. hordei (Bgh), is an obligate biotrophic ascomycete fungal pathogen that can grow and reproduce only on living cells of wild or domesticated barley (Hordeum sp.). Domestication and deployment of resistant barley cultivars by humans selected for amplification of Bgh isolates with different virulence combinations. We sequenced the genomes of two European Bgh isolates, A6 and K1, for comparative analysis with the reference genome of isolate DH14. This revealed a mosaic genome structure consisting of large isolate-specific DNA blocks with either high or low SNP densities. Some of the highly polymorphic blocks likely accumulated SNPs for over 10,000 years, well before the domestication of barley. These isolate-specific blocks of alternating monomorphic and polymorphic regions imply an exceptionally large standing genetic variation in the Bgh population and might be generated and maintained by rare outbreeding and frequent clonal reproduction. RNA-sequencing experiments with isolates A6 and K1 during four early stages of compatible and incompatible interactions on leaves of partially immunocompromised Arabidopsis mutants revealed a conserved Bgh transcriptional program during pathogenesis compared with the natural host barley despite ~200 million years of reproductive isolation of these hosts. Transcripts encoding candidate-secreted effector proteins are massively induced in successive waves. A specific decrease in candidate-secreted effector protein transcript abundance in the incompatible interaction follows extensive transcriptional reprogramming of the host transcriptome and coincides with the onset of localized host cell death, suggesting a host-inducible defense mechanism that targets fungal effector secretion or production.

Host-induced gene silencing in barley powdery mildew reveals a class of ribonuclease-like effectors

Obligate biotrophic pathogens of plants must circumvent or counteract defenses to guarantee accommodation inside the host. To do so, they secrete a variety of effectors that regulate host immunity and facilitate the establishment of pathogen feeding structures called haustoria. The barley powdery mildew fungus Blumeria graminis f. sp. hordei produces a large number of proteins predicted to be secreted from haustoria. Fifty of these Blumeria effector candidates (BEC) were screened by host-induced gene silencing (HIGS), and eight were identified that contribute to infection. One shows similarity to β-1,3 glucosyltransferases, one to metallo-proteases, and two to microbial secreted ribonucleases; the remainder have no similarity to proteins of known function. Transcript abundance of all eight BEC increases dramatically in the early stages of infection and establishment of haustoria, consistent with a role in that process. Complementation analysis using silencing-insensitive synthetic cDNAs demonstrated that the ribonuclease-like BEC 1011 and 1054 are bona fide effectors that function within the plant cell. BEC1011 specifically interferes with pathogen-induced host cell death. Both are part of a gene superfamily unique to the powdery mildew fungi. Structural modeling was consistent, with BEC1054 adopting a ribonuclease-like fold, a scaffold not previously associated with effector function.

Communication between filamentous pathogens and plants at the biotrophic interface

Fungi and oomycetes that colonize living plant tissue form extensive interfaces with plant cells in which the cytoplasm of the microorganism is closely aligned with the host cytoplasm for an extended distance. In all cases, specialized biotrophic hyphae function to hijack host cellular processes across an interfacial zone consisting of a hyphal plasma membrane, a specialized interfacial matrix, and a plant-derived membrane. The interface is the site of active secretion by both players. This cross talk at the interface determines the winner in adversarial relationships and establishes the partnership in mutualistic relationships. Fungi and oomycetes secrete many specialized effector proteins for controlling the host, and they can stimulate remarkable cellular reorganization even in distant plant cells. Breakthroughs in live-cell imaging of fungal and oomycete encounter sites, including live-cell imaging of pathogens secreting fluorescently labeled effector proteins, have led to recent progress in understanding communication across the interface.

SNP design from 454 sequencing of Podosphaera plantaginis transcriptome reveals a genetically diverse pathogen metapopulation with high levels of mixed-genotype infection

Background

Molecular tools may greatly improve our understanding of pathogen evolution and epidemiology but technical constraints have hindered the development of genetic resources for parasites compared to free-living organisms. This study aims at developing molecular tools for Podosphaera plantaginis, an obligate fungal pathogen of Plantago lanceolata. This interaction has been intensively studied in the Åland archipelago of Finland with epidemiological data collected from over 4,000 host populations annually since year 2001.

Principal Findings

A cDNA library of a pooled sample of fungal conidia was sequenced on the 454 GS-FLX platform. Over 549,411 reads were obtained and annotated into 45,245 contigs. Annotation data was acquired for 65.2% of the assembled sequences. The transcriptome assembly was screened for SNP loci, as well as for functionally important genes (mating-type genes and potential effector proteins). A genotyping assay of 27 SNP loci was designed and tested on 380 infected leaf samples from 80 populations within the Åland archipelago. With this panel we identified 85 multilocus genotypes (MLG) with uneven frequencies across the pathogen metapopulation. Approximately half of the sampled populations contain polymorphism. Our genotyping protocol revealed mixed-genotype infection within a single host leaf to be common. Mixed infection has been proposed as one of the main drivers of pathogen evolution, and hence may be an important process in this pathosystem.

Significance

The developed SNP panel offers exciting research perspectives for future studies in this well-characterized pathosystem. Also, the transcriptome provides an invaluable novel genomic resource for powdery mildews, which cause significant yield losses on commercially important crops annually. Furthermore, the features that render genetic studies in this system a challenge are shared with the majority of obligate parasitic species, and hence our results provide methodological insights from SNP calling to field sampling protocols for a wide range of biological systems.

Structure and evolution of barley powdery mildew effector candidates

Background

Protein effectors of pathogenicity are instrumental in modulating host immunity and disease resistance. The powdery mildew pathogen of grasses Blumeria graminis causes one of the most important diseases of cereal crops. B. graminis is an obligate biotrophic pathogen and as such has an absolute requirement to suppress or avoid host immunity if it is to survive and cause disease.

Results

Here we characterise a superfamily predicted to be the full complement of Candidates for Secreted Effector Proteins (CSEPs) in the fungal barley powdery mildew parasite B. graminis f.sp. hordei. The 491 genes encoding these proteins constitute over 7% of this pathogen’s annotated genes and most were grouped into 72 families of up to 59 members. They were predominantly expressed in the intracellular feeding structures called haustoria, and proteins specifically associated with the haustoria were identified by large-scale mass spectrometry-based proteomics. There are two major types of effector families: one comprises shorter proteins (100–150 amino acids), with a high relative expression level in the haustoria and evidence of extensive diversifying selection between paralogs; the second type consists of longer proteins (300–400 amino acids), with lower levels of differential expression and evidence of purifying selection between paralogs. An analysis of the predicted protein structures underscores their overall similarity to known fungal effectors, but also highlights unexpected structural affinities to ribonucleases throughout the entire effector super-family. Candidate effector genes belonging to the same family are loosely clustered in the genome and are associated with repetitive DNA derived from retro-transposons.

Conclusions

We employed the full complement of genomic, transcriptomic and proteomic analyses as well as structural prediction methods to identify and characterize the members of the CSEPs superfamily in B. graminis f.sp. hordei. Based on relative intron position and the distribution of CSEPs with a ribonuclease-like domain in the phylogenetic tree we hypothesize that the associated genes originated from an ancestral gene, encoding a secreted ribonuclease, duplicated successively by repetitive DNA-driven processes and diversified during the evolution of the grass and cereal powdery mildew lineage.

Interaction of barley powdery mildew effector candidate CSEP0055 with the defence protein PR17c

A large number of effector candidates have been identified recently in powdery mildew fungi. However, their roles and how they perform their functions remain unresolved. In this study, we made use of host-induced gene silencing and confirmed that the secreted barley powdery mildew effector candidate, CSEP0055, contributes to the aggressiveness of the fungus. This result suggests that CSEP0055 is involved in the suppression of plant defence. A yeast two-hybrid screen indicated that CSEP0055 interacts with members of the barley pathogenesis-related protein families, PR1 and PR17. Interaction with PR17c was confirmed by bimolecular fluorescence complementation analyses. Down-regulation and over-expression of PR17c in epidermal cells of barley confirmed that this protein is important for penetration resistance against the powdery mildew fungus. In line with this, PR17c was found to be apoplastic, localizing to the papillae formed in response to this fungus. The CSEP0055 transcript did not start to accumulate until 24 h after inoculation. This suggests that this gene is expressed too late to influence primary penetration events, but rather sustains the fungus at sites of secondary penetration, where PR17c appears to be able to accumulate.

Genomic tillage and the harvest of fungal phytopathogens

Genome sequencing has been carried out on a small selection of major fungal ascomycete pathogens. These studies show that simple models whereby pathogens evolved from phylogenetically related saprobes by the acquisition or modification of a small number of key genes cannot be sustained.The genomes show that pathogens cannot be divided into three clearly delineated classes (biotrophs, hemibiotrophs and necrotrophs) but rather into a complex matrix of categories each with subtly different properties. It is clear that the evolution of pathogenicity is ancient, rapid and ongoing. Fungal pathogens have undergone substantial genomic rearrangements that can be appropriately described as 'genomic tillage'. Genomic tillage underpins the evolution and expression of large families of genes - known as effectors - that manipulate and exploit metabolic and defence processes of plants so as to allow the proliferation of pathogens.

A genetic mechanism for emergence of races in Fusarium oxysporum f. sp. lycopersici: inactivation of avirulence gene AVR1 by transposon insertion

Compatible/incompatible interactions between the tomato wilt fungus Fusarium oxysporum f. sp. lycopersici (FOL) and tomato Solanum lycopersicum are controlled by three avirulence genes (AVR1–3) in FOL and the corresponding resistance genes (I–I3) in tomato. The three known races (1, 2 and 3) of FOL carry AVR genes in different combinations. The current model to explain the proposed order of mutations in AVR genes is: i) FOL race 2 emerged from race 1 by losing the AVR1 and thus avoiding host resistance mediated by I (the resistance gene corresponding to AVR1), and ii) race 3 emerged when race 2 sustained a point mutation in AVR2, allowing it to evade I2-mediated resistance of the host. Here, an alternative mechanism of mutation of AVR genes was determined by analyses of a race 3 isolate, KoChi-1, that we recovered from a Japanese tomato field in 2008. Although KoChi-1 is race 3, it has an AVR1 gene that is truncated by the transposon Hormin, which belongs to the hAT family. This provides evidence that mobile genetic elements may be one of the driving forces underlying race evolution. KoChi-1 transformants carrying a wild type AVR1 gene from race 1 lost pathogenicity to cultivars carrying I, showing that the truncated KoChi-1 avr1 is not functional. These results imply that KoChi-1 is a new race 3 biotype and propose an additional path for emergence of FOL races: Race 2 emerged from race 1 by transposon-insertion into AVR1, not by deletion of the AVR1 locus; then a point mutation in race 2 AVR2 resulted in emergence of race 3.

Phylogeography and virulence structure of the powdery mildew population on its 'new' host triticale

BACKGROUND

Powdery mildew, caused by the obligate biotrophic fungus Blumeria graminis, is a major problem in cereal production as it can reduce quality and yield. B. graminis has evolved eight distinct formae speciales (f.sp.) which display strict host specialization. In the last decade, powdery mildew has emerged on triticale, the artificial intergeneric hybrid between wheat and rye. This emergence is probably triggered by a host range expansion of the wheat powdery mildew B. graminis f.sp. tritici. To gain more precise information about the evolutionary processes that led to this host range expansion, we pursued a combined pathological and genetic approach.

RESULTS

B. graminis isolates were sampled from triticale, wheat and rye from different breeding regions in Europe. Pathogenicity tests showed that isolates collected from triticale are highly pathogenic on most of the tested triticale cultivars. Moreover, these isolates were also able to infect several wheat cultivars (their previous hosts), although a lower aggressiveness was observed compared to isolates collected from wheat. Phylogenetic analysis of nuclear gene regions identified two statistically significant clades, which to a certain extent correlated with pathogenicity. No differences in virulence profiles were found among the sampled regions, but the distribution of genetic variation demonstrated to be geography dependent. A multilocus haplotype network showed that haplotypes pathogenic on triticale are distributed at different sites in the network, but always clustered at or near the tips of the network.

CONCLUSIONS

This study reveals a genetic structure in B. graminis with population differentiation according to geography and host specificity. In addition, evidence is brought forward demonstrating that the host range expansion of wheat isolates to the new host triticale occurred recently and multiple times at different locations in Europe.

Opprimo ergo sum--evasion and suppression in the root endophytic fungus Piriformospora indica

The genetically tractable endophytic fungus Piriformospora indica is able to colonize the root cortex of a great variety of different plant species with beneficial effects to its hosts, and it represents a suitable model system to study symbiotic interactions. Recent cytological studies in barley and Arabidopsis showed that, upon penetration of the root, P. indica establishes a biotrophic interaction during which fungal cells are encased by the host plasma membrane. Large-scale transcriptional analyses of fungal and plant responses have shown that perturbance of plant hormone homeostasis and secretion of fungal lectins and other small proteins (effectors) may be involved in the evasion and suppression of host defenses at these early colonization steps. At later stages, P. indica is found more often in moribund host cells where it secretes a large variety of hydrolytic enzymes that degrade proteins. This strategy of colonizing plants is reminiscent of that of hemibiotrophic fungi, although a defined shift to necrotrophy with massive host cell death is missing. Instead, the association with the plant root leads to beneficial effects for the host such as growth promotion, increased resistance to root as well as leaf pathogens, and increased tolerance to abiotic stresses. This review describes current advances in understanding the components of the P. indica endophytic lifestyle from molecular and genomic analyses.

The Top 10 fungal pathogens in molecular plant pathology

The aim of this review was to survey all fungal pathologists with an association with the journal Molecular Plant Pathology and ask them to nominate which fungal pathogens they would place in a 'Top 10' based on scientific/economic importance. The survey generated 495 votes from the international community, and resulted in the generation of a Top 10 fungal plant pathogen list for Molecular Plant Pathology. The Top 10 list includes, in rank order, (1) Magnaporthe oryzae; (2) Botrytis cinerea; (3) Puccinia spp.; (4) Fusarium graminearum; (5) Fusarium oxysporum; (6) Blumeria graminis; (7) Mycosphaerella graminicola; (8) Colletotrichum spp.; (9) Ustilago maydis; (10) Melampsora lini, with honourable mentions for fungi just missing out on the Top 10, including Phakopsora pachyrhizi and Rhizoctonia solani. This article presents a short resumé of each fungus in the Top 10 list and its importance, with the intent of initiating discussion and debate amongst the plant mycology community, as well as laying down a bench-mark. It will be interesting to see in future years how perceptions change and what fungi will comprise any future Top 10.

The Top 10 fungal pathogens in molecular plant pathology

The aim of this review was to survey all fungal pathologists with an association with the journal Molecular Plant Pathology and ask them to nominate which fungal pathogens they would place in a 'Top 10' based on scientific/economic importance. The survey generated 495 votes from the international community, and resulted in the generation of a Top 10 fungal plant pathogen list for Molecular Plant Pathology. The Top 10 list includes, in rank order, (1) Magnaporthe oryzae; (2) Botrytis cinerea; (3) Puccinia spp.; (4) Fusarium graminearum; (5) Fusarium oxysporum; (6) Blumeria graminis; (7) Mycosphaerella graminicola; (8) Colletotrichum spp.; (9) Ustilago maydis; (10) Melampsora lini, with honourable mentions for fungi just missing out on the Top 10, including Phakopsora pachyrhizi and Rhizoctonia solani. This article presents a short resumé of each fungus in the Top 10 list and its importance, with the intent of initiating discussion and debate amongst the plant mycology community, as well as laying down a bench-mark. It will be interesting to see in future years how perceptions change and what fungi will comprise any future Top 10.

Tomato immune receptor Ve1 recognizes effector of multiple fungal pathogens uncovered by genome and RNA sequencing

Fungal plant pathogens secrete effector molecules to establish disease on their hosts, and plants in turn use immune receptors to try to intercept these effectors. The tomato immune receptor Ve1 governs resistance to race 1 strains of the soil-borne vascular wilt fungi Verticillium dahliae and Verticillium albo-atrum, but the corresponding Verticillium effector remained unknown thus far. By high-throughput population genome sequencing, a single 50-Kb sequence stretch was identified that only occurs in race 1 strains, and subsequent transcriptome sequencing of Verticillium-infected Nicotiana benthamiana plants revealed only a single highly expressed ORF in this region, designated Ave1 (for Avirulence on Ve1 tomato). Functional analyses confirmed that Ave1 activates Ve1-mediated resistance and demonstrated that Ave1 markedly contributes to fungal virulence, not only on tomato but also on Arabidopsis. Interestingly, Ave1 is homologous to a widespread family of plant natriuretic peptides. Besides plants, homologous proteins were only found in the bacterial plant pathogen Xanthomonas axonopodis and the plant pathogenic fungi Colletotrichum higginsianum, Cercospora beticola, and Fusarium oxysporum f. sp. lycopersici. The distribution of Ave1 homologs, coincident with the presence of Ave1 within a flexible genomic region, strongly suggests that Verticillium acquired Ave1 from plants through horizontal gene transfer. Remarkably, by transient expression we show that also the Ave1 homologs from F. oxysporum and C. beticola can activate Ve1-mediated resistance. In line with this observation, Ve1 was found to mediate resistance toward F. oxysporum in tomato, showing that this immune receptor is involved in resistance against multiple fungal pathogens.

The genomics of obligate (and nonobligate) biotrophs

Biotrophy is a pervasive trait that evolved independently in plant pathogenic fungi and oomycetes. Comparative genomics of the first sequenced biotrophic pathogens highlight remarkable convergences, including gene losses in the metabolism of inorganic nitrogen, inorganic sulfur, and thiamine, and genes encoding carbohydrate active enzymes and secondary metabolism enzymes. Some biotrophs, but not all, display marked increases in overall genome size because of a proliferation of retrotransposons. I argue here that the release of constraints on transposon activity is driven by the advantages conferred by the genetic variability that results from transposition, in particular by the creation and diversification of broad palettes of effector genes. Increases in genome size and gene losses are the consequences of this trade-off. Genes that are not necessary for growth on a plant disappeared, but we still do not know what lost functions make some of these pathogens obligate.

Grapevine powdery mildew (Erysiphe necator): a fascinating system for the study of the biology, ecology and epidemiology of an obligate biotroph

Few plant pathogens have had a more profound effect on the evolution of disease management than Erysiphe necator, which causes grapevine powdery mildew. When the pathogen first spread from North America to England in 1845, and onwards to France in 1847, 'germ theory' was neither understood among the general populace nor even generally accepted within the scientific community. Louis Pasteur had only recently reported the microbial nature of fermentation, and it would be another 30 years before Robert Koch would publish his proofs of the microbial nature of certain animal diseases. However, within 6 years after the arrival of the pathogen, nearly 6 million grape growers in France were routinely applying sulphur to suppress powdery mildew on nearly 2.5 million hectares of vineyards (Campbell, 2006). The pathogen has remained a focus for disease management efforts ever since. Because of the worldwide importance of the crop and its susceptibility to the disease, and because conventional management with modern, organic fungicides has been compromised on several occasions since 1980 by the evolution of fungicide resistance, there has also been a renewed effort worldwide to explore the pathogen's biology and ecology, its genetics and molecular interactions with host plants, and to refine current and suggest new management strategies. These latter aspects are the subject of our review.

TAXONOMY

The most widely accepted classification follows. Family Erysiphaceae, Erysiphe necator Schw. [syn. Uncinula necator (Schw.) Burr., E. tuckeri Berk., U. americana Howe and U. spiralis Berk. & Curt; anamorph Oidium tuckeri Berk.]. Erysiphe necator var. ampelopsidis was found on Parthenocissus spp. in North America according to Braun (1987), although later studies revealed isolates whose host range spanned genera, making the application of this taxon somewhat imprecise (Gadoury and Pearson, 1991). The classification of the genera before 1980 was based on features of the mature ascocarp: (i) numbers of asci; and (ii) morphology of the appendages, in particular the appendage tips. The foregoing has been supplanted by phylogeny inferred from the internal transcribed spacer (ITS) of ribosomal DNA sequences (Saenz and Taylor, 1999), which correlates with conidial ontogeny and morphology (Braun et al., 2002).

HOST RANGE

The pathogen is obligately parasitic on genera within the Vitaceae, including Vitis, Cissus, Parthenocissus and Ampelopsis (Pearson and Gadoury, 1992). The most economically important host is grapevine (Vitis), particularly the European grape, V. vinifera, which is highly susceptible to powdery mildew. Disease symptoms and signs: In the strictest sense, macroscopically visible mildew colonies are signs of the pathogen rather than symptoms resulting from its infection, but, for convenience, we describe the symptoms and signs together as the collective appearance of colonized host tissues. All green tissues of the host may be infected. Ascospore colonies are most commonly found on the lower surface of the first-formed leaves near the bark of the vine, and may be accompanied by a similarly shaped chlorotic spot on the upper surface. Young colonies appear whitish and those that have not yet sporulated show a metallic sheen. They are roughly circular, ranging in size from a few millimetres to a centimetre or more in diameter, and can occur singly or in groups that coalesce to cover much of the leaf. Senescent colonies are greyish, and may bear cleistothecia in various stages of development. Dead epidermal cells often subtend the colonized area, as natural mortality in the mildew colony, the use of fungicides, mycoparasites or resistance responses in the leaf result in the deaths of segments of the mildew colony and infected epidermal cells. Severely affected leaves usually senesce, develop necrotic blotches and fall prematurely. Infection of stems initially produces symptoms similar to those on leaves, but colonies on shoots are eventually killed as periderm forms, producing a dark, web-like scar on the cane (Gadoury et al., 2011). Inflorescences and berries are most susceptible when young, and can become completely coated with whitish mildew. The growth of the berry epidermal tissue stops when severely infected, which may result in splitting as young fruit expand. Berries in a transitional stage between susceptible and resistant (generally between 3 and 4 weeks after anthesis) develop diffuse, nonsporulating mildew colonies only visible under magnification. Diffuse colonies die as berries continue to mature, leaving behind a network of necrotic epidermal cells (Gadoury et al., 2007). Survival over winter as mycelium in buds results in a distinctive foliar symptom. Shoots arising from these buds may be heavily coated with fungal growth, stark white in colour and stand out like white flags in the vine, resulting in the term 'flag shoots'. More commonly, colonization of a flag shoot is less extensive, and infection of a single leaf, or of leaves on one side of the shoot only, is observed (Gadoury et al., 2011).

Endophytic life strategies decoded by genome and transcriptome analyses of the mutualistic root symbiont Piriformospora indica

Recent sequencing projects have provided deep insight into fungal lifestyle-associated genomic adaptations. Here we report on the 25 Mb genome of the mutualistic root symbiont Piriformospora indica (Sebacinales, Basidiomycota) and provide a global characterization of fungal transcriptional responses associated with the colonization of living and dead barley roots. Extensive comparative analysis of the P. indica genome with other Basidiomycota and Ascomycota fungi that have diverse lifestyle strategies identified features typically associated with both, biotrophism and saprotrophism. The tightly controlled expression of the lifestyle-associated gene sets during the onset of the symbiosis, revealed by microarray analysis, argues for a biphasic root colonization strategy of P. indica. This is supported by a cytological study that shows an early biotrophic growth followed by a cell death-associated phase. About 10% of the fungal genes induced during the biotrophic colonization encoded putative small secreted proteins (SSP), including several lectin-like proteins and members of a P. indica-specific gene family (DELD) with a conserved novel seven-amino acids motif at the C-terminus. Similar to effectors found in other filamentous organisms, the occurrence of the DELDs correlated with the presence of transposable elements in gene-poor repeat-rich regions of the genome. This is the first in depth genomic study describing a mutualistic symbiont with a biphasic lifestyle. Our findings provide a significant advance in understanding development of biotrophic plant symbionts and suggest a series of incremental shifts along the continuum from saprotrophy towards biotrophy in the evolution of mycorrhizal association from decomposer fungi.

Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations

Fungi are of primary ecological, biotechnological and economic importance. Many fundamental biological processes that are shared by animals and fungi are studied in fungi due to their experimental tractability. Many fungi are pathogens or mutualists and are model systems to analyse effector genes and their mechanisms of diversification. In this study, we report the genome sequence of the phytopathogenic ascomycete Leptosphaeria maculans and characterize its repertoire of protein effectors. The L. maculans genome has an unusual bipartite structure with alternating distinct guanine and cytosine-equilibrated and adenine and thymine (AT)-rich blocks of homogenous nucleotide composition. The AT-rich blocks comprise one-third of the genome and contain effector genes and families of transposable elements, both of which are affected by repeat-induced point mutation, a fungal-specific genome defence mechanism. This genomic environment for effectors promotes rapid sequence diversification and underpins the evolutionary potential of the fungus to adapt rapidly to novel host-derived constraints.

Plant-parasite coevolution: bridging the gap between genetics and ecology

We review current ideas about coevolution of plants and parasites, particularly processes that generate genetic diversity. Frequencies of host resistance and parasite virulence alleles that interact in gene-for-gene (GFG) relationships coevolve in the familiar boom-and-bust cycle, in which resistance is selected when virulence is rare, and virulence is selected when resistance is common. The cycle can result in stable polymorphism when diverse ecological and epidemiological factors cause negative direct frequency-dependent selection (ndFDS) on host resistance, parasite virulence, or both, such that the benefit of a trait to fitness declines as its frequency increases. Polymorphism can also be stabilized by overdominance, when heterozygous hosts have greater resistance than homozygotes to diverse pathogens. Genetic diversity can also persist in the form of statistical polymorphism, sustained by random processes acting on gene frequencies and population size. Stable polymorphism allows alleles to be long-lived and genetic variation to be detectable in natural populations. In agriculture, many of the factors promoting stability in host-parasite interactions have been lost, leading to arms races of host defenses and parasite effectors.

Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism

Powdery mildews are phytopathogens whose growth and reproduction are entirely dependent on living plant cells. The molecular basis of this life-style, obligate biotrophy, remains unknown. We present the genome analysis of barley powdery mildew, Blumeria graminis f.sp. hordei (Blumeria), as well as a comparison with the analysis of two powdery mildews pathogenic on dicotyledonous plants. These genomes display massive retrotransposon proliferation, genome-size expansion, and gene losses. The missing genes encode enzymes of primary and secondary metabolism, carbohydrate-active enzymes, and transporters, probably reflecting their redundancy in an exclusively biotrophic life-style. Among the 248 candidate effectors of pathogenesis identified in the Blumeria genome, very few (less than 10) define a core set conserved in all three mildews, suggesting that most effectors represent species-specific adaptations.

HIGS: host-induced gene silencing in the obligate biotrophic fungal pathogen Blumeria graminis

Powdery mildew fungi are obligate biotrophic pathogens that only grow on living hosts and cause damage in thousands of plant species. Despite their agronomical importance, little direct functional evidence for genes of pathogenicity and virulence is currently available because mutagenesis and transformation protocols are lacking. Here, we show that the accumulation in barley (Hordeum vulgare) and wheat (Triticum aestivum) of double-stranded or antisense RNA targeting fungal transcripts affects the development of the powdery mildew fungus Blumeria graminis. Proof of concept for host-induced gene silencing was obtained by silencing the effector gene Avra10, which resulted in reduced fungal development in the absence, but not in the presence, of the matching resistance gene Mla10. The fungus could be rescued from the silencing of Avra10 by the transient expression of a synthetic gene that was resistant to RNA interference (RNAi) due to silent point mutations. The results suggest traffic of RNA molecules from host plants into B. graminis and may lead to an RNAi-based crop protection strategy against fungal pathogens.

Powdery mildew fungal effector candidates share N-terminal Y/F/WxC-motif

Background

Powdery mildew and rust fungi are widespread, serious pathogens that depend on developing haustoria in the living plant cells. Haustoria are separated from the host cytoplasm by a plant cell-derived extrahaustorial membrane. They secrete effector proteins, some of which are subsequently transferred across this membrane to the plant cell to suppress defense.

Results

In a cDNA library from barley epidermis containing powdery mildew haustoria, two-thirds of the sequenced ESTs were fungal and represented ~3,000 genes. Many of the most highly expressed genes encoded small proteins with N-terminal signal peptides. While these proteins are novel and poorly related, they do share a three-amino acid motif, which we named "Y/F/WxC", in the N-terminal of the mature proteins. The first amino acid of this motif is aromatic: tyrosine, phenylalanine or tryptophan, and the last is always cysteine. In total, we identified 107 such proteins, for which the ESTs represent 19% of the fungal clones in our library, suggesting fundamental roles in haustoria function. While overall sequence similarity between the powdery mildew Y/F/WxC-proteins is low, they do have a highly similar exon-intron structure, suggesting they have a common origin. Interestingly, searches of public fungal genome and EST databases revealed that haustoria-producing rust fungi also encode large numbers of novel, short proteins with signal peptides and the Y/F/WxC-motif. No significant numbers of such proteins were identified from genome and EST sequences from either fungi which do not produce haustoria or from haustoria-producing Oomycetes.

Conclusion

In total, we identified 107, 178 and 57 such Y/F/WxC-proteins from the barley powdery mildew, the wheat stem rust and the wheat leaf rust fungi, respectively. All together, our findings suggest the Y/F/WxC-proteins to be a new class of effectors from haustoria-producing pathogenic fungi.

Unique evolutionary pattern of numbers of gramineous NBS-LRR genes

Nucleotide binding site (NBS)-leucine-rich repeat (LRR) genes belong to the largest class of disease-resistance gene super groups in plants, and their intra- or interspecies nucleotide variations have been studied extensively to understand their evolution and function. However, little is known about the evolutionary patterns of their copy numbers in related species. Here, 129, 245, 239 and 508 NBSs were identified in maize, sorghum, brachypodium and rice, respectively, suggesting considerable variations of these genes. Based on phylogenetic relationships from a total of 496 ancestral branches of grass NBS families, three gene number variation patterns were categorized: conserved, sharing two or more species, and species-specific. Notably, the species-specific NBS branches are dominant (71.6%), while there is only a small percentage (3.83%) of conserved families. In contrast, the conserved families are dominant in 51 randomly selected house-keeping genes (96.1%). The opposite patterns between NBS and the other gene groups suggest that natural selection is responsible for the drastic number variation of NBS genes. The rapid expansion and/or contraction may be a fundamentally important strategy for a species to adapt to the quickly changing species-specific pathogen spectrum. In addition, the small proportion of conserved NBSs suggests that the loss of NBSs may be a general tendency in grass species.

Unique evolutionary pattern of numbers of gramineous NBS-LRR genes

Nucleotide binding site (NBS)-leucine-rich repeat (LRR) genes belong to the largest class of disease-resistance gene super groups in plants, and their intra- or interspecies nucleotide variations have been studied extensively to understand their evolution and function. However, little is known about the evolutionary patterns of their copy numbers in related species. Here, 129, 245, 239 and 508 NBSs were identified in maize, sorghum, brachypodium and rice, respectively, suggesting considerable variations of these genes. Based on phylogenetic relationships from a total of 496 ancestral branches of grass NBS families, three gene number variation patterns were categorized: conserved, sharing two or more species, and species-specific. Notably, the species-specific NBS branches are dominant (71.6%), while there is only a small percentage (3.83%) of conserved families. In contrast, the conserved families are dominant in 51 randomly selected house-keeping genes (96.1%). The opposite patterns between NBS and the other gene groups suggest that natural selection is responsible for the drastic number variation of NBS genes. The rapid expansion and/or contraction may be a fundamentally important strategy for a species to adapt to the quickly changing species-specific pathogen spectrum. In addition, the small proportion of conserved NBSs suggests that the loss of NBSs may be a general tendency in grass species.