Allelic variation in the effector genes of the tomato pathogen Cladosporium fulvum reveals different modes of adaptive evolution

Sequential breakdown of the Cf-9 leaf mould resistance locus in tomato by Fulvia fulva

Leaf mould, caused by Fulvia fulva, is a devastating disease of tomato plants. In many commercial tomato cultivars, resistance to this disease is governed by the Cf-9 locus, which encodes five paralogous receptor-like proteins. Two of these proteins confer resistance: Cf-9C recognises the previously identified F. fulva effector Avr9 and provides resistance during all plant growth stages, while Cf-9B recognises the yet-unidentified F. fulva effector Avr9B and provides mature plant resistance only. In recent years, F. fulva strains have emerged that can overcome the Cf-9 locus, with Cf-9C circumvented through Avr9 deletion. To understand how Cf-9B is circumvented, we set out to identify Avr9B. Comparative genomics, transient expression assays and gene complementation experiments were used to identify Avr9B, while gene sequencing was used to assess Avr9B allelic variation across a world-wide strain collection. A strict correlation between Avr9 deletion and resistance-breaking mutations in Avr9B was observed in strains recently collected from Cf-9 cultivars, whereas Avr9 deletion but no mutations in Avr9B were observed in older strains. This research showcases how F. fulva has evolved to sequentially break down the Cf-9 locus and stresses the urgent need for commercial tomato cultivars that carry novel, stacked resistance genes active against this pathogen.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BACKGROUND

Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

RESULTS

Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.

CONCLUSIONS

Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

Beyond the genomes of Fulvia fulva (syn. Cladosporium fulvum) and Dothistroma septosporum: New insights into how these fungal pathogens interact with their host plants

Fulvia fulva and Dothistroma septosporum are closely related apoplastic pathogens with similar lifestyles but different hosts: F. fulva is a pathogen of tomato, whilst D. septosporum is a pathogen of pine trees. In 2012, the first genome sequences of these pathogens were published, with F. fulva and D. septosporum having highly fragmented and near-complete assemblies, respectively. Since then, significant advances have been made in unravelling their genome architectures. For instance, the genome of F. fulva has now been assembled into 14 chromosomes, 13 of which have synteny with the 14 chromosomes of D. septosporum, suggesting these pathogens are even more closely related than originally thought. Considerable advances have also been made in the identification and functional characterization of virulence factors (e.g., effector proteins and secondary metabolites) from these pathogens, thereby providing new insights into how they promote host colonization or activate plant defence responses. For example, it has now been established that effector proteins from both F. fulva and D. septosporum interact with cell-surface immune receptors and co-receptors to activate the plant immune system. Progress has also been made in understanding how F. fulva and D. septosporum have evolved with their host plants, whilst intensive research into pandemics of Dothistroma needle blight in the Northern Hemisphere has shed light on the origins, migration, and genetic diversity of the global D. septosporum population. In this review, we specifically summarize advances made in our understanding of the F. fulva-tomato and D. septosporum-pine pathosystems over the last 10 years.

Progress in pathogenesis research of Ustilago maydis, and the metabolites involved along with their biosynthesis

Ustilago maydis is a pathogenic fungus that causes corn smut. Because of its easy cultivation and genetic transformation, U. maydis has become an important model organism for plant-pathogenic basidiomycetes. U. maydis is able to infect maize by producing effectors and secreted proteins as well as surfactant-like metabolites. In addition, the production of melanin and iron carriers is also associated with its pathogenicity. Here, advances in our understanding of the pathogenicity of U. maydis, the metabolites involved in the pathogenic process, and the biosynthesis of these metabolites, are reviewed and discussed. This summary will provide new insights into the pathogenicity of U. maydis and the functions of associated metabolites, as well as new clues for deciphering the biosynthesis of metabolites.

A highly polymorphic effector protein promotes fungal virulence through suppression of plant-associated Actinobacteria

Plant pathogens secrete effector proteins to support host colonization through a wide range of molecular mechanisms, while plant immune systems evolved receptors to recognize effectors or their activities to mount immune responses to halt pathogens. Importantly, plants do not act as single organisms, but rather as holobionts that actively shape their microbiota as a determinant of health. The soil-borne fungal pathogen Verticillium dahliae was recently demonstrated to exploit the VdAve1 effector to manipulate the host microbiota to promote vascular wilt disease in the absence of the corresponding immune receptor Ve1. We identify a multiallelic V. dahliae gene displaying c. 65% sequence similarity to VdAve1, named VdAve1-like (VdAve1L), which shows extreme sequence variation, including alleles that encode dysfunctional proteins, indicative of selection pressure to overcome host recognition. We show that the orphan cell surface receptor Ve2, encoded at the Ve locus, does not recognize VdAve1L. Additionally, we demonstrate that the full-length variant VdAve1L2 possesses antimicrobial activity, like VdAve1, yet with a divergent activity spectrum, that is exploited by V. dahliae to mediate tomato colonization through the direct suppression of antagonistic Actinobacteria in the host microbiota. Our findings open up strategies for more targeted biocontrol against microbial plant pathogens.

Improving sustainable crop protection using population genetics concepts

Growing genetically resistant plants allows pathogen populations to be controlled and reduces the use of pesticides. However, pathogens can quickly overcome such resistance. In this context, how can we achieve sustainable crop protection? This crucial question has remained largely unanswered despite decades of intense debate and research effort. In this study, we used a bibliographic analysis to show that the research field of resistance durability has evolved into three subfields: (i) 'plant breeding' (generating new genetic material), (ii) 'molecular interactions' (exploring the molecular dialogue governing plant-pathogen interactions) and (iii) 'epidemiology and evolution' (explaining and forecasting of pathogen population dynamics resulting from selection pressure(s) exerted by resistant plants). We argue that this triple split of the field impedes integrated research progress and ultimately compromises the sustainable management of genetic resistance. After identifying a gap among the three subfields, we argue that the theoretical framework of population genetics could bridge this gap. Indeed, population genetics formally explains the evolution of all heritable traits, and allows genetic changes to be tracked along with variation in population dynamics. This provides an integrated view of pathogen adaptation, in particular via evolutionary-epidemiological feedbacks. In this Opinion Note, we detail examples illustrating how such a framework can better inform best practices for developing and managing genetically resistant cultivars.

Identification of Pathogenicity-Related Effector Proteins and the Role of Piwsc1 in the Virulence of Penicillium italicum on Citrus Fruits

Blue mold caused by Penicillium italicum is one of the two major postharvest diseases of citrus fruits. The interactions of pathogens with their hosts are complicated, and virulence factors that mediate pathogenicity have not yet been identified. In present study, a prediction pipeline approach based on bioinformatics and transcriptomic data is designed to determine the effector proteins of P. italicum. Three hundred and seventy-five secreted proteins of P. italicum were identified, many of which (29.07%) were enzymes for carbohydrate utilization. Twenty-nine candidates were further analyzed and the expression patterns of 12 randomly selected candidate effector genes were monitored during the early stages of growth on PDA and infection of Navel oranges for validation. Functional analysis of a cell wall integrity-related gene Piwsc1, a core candidate, was performed by gene knockout. The deletion of Piwsc1 resulted in reduced virulence on citrus fruits, as presented by an approximate 57% reduction in the diameter of lesions. In addition, the mycelial growth rate, spore germination rate, and sporulation of ΔPiwsc1 decreased. The findings provide us with new insights to understand the pathogenesis of P. italicum and develop an effective and sustainable control method for blue mold.

The Diversity of Passalora fulva Isolates Collected from Tomato Plants in U.S. High Tunnels

High tunnels extend the growing season of high value crops, including tomatoes, but the environmental conditions within high tunnels favor the spread of the tomato leaf mold pathogen, Passalora fulva (syn. Cladosporium fulvum). Tomato leaf mold results in defoliation, and if severe, losses in yield. Despite substantial research, little is known regarding the genetic structure and diversity of populations of P. fulva associated with high tunnel tomato production in the United States. From 2016 to 2019, a total of 50 P. fulva isolates were collected from tomato leaf samples in high tunnels in the Northeast and Minnesota. Other Cladosporium species were also isolated from the leaf surfaces. Koch's postulates were conducted to confirm that P. fulva was the cause of the disease symptoms observed. Race determination experiments revealed that the isolates belonged to either race 0 (six isolates) or race 2 (44 isolates). Polymorphisms were identified within four previously characterized effector genes: Avr2, Avr4, Avr4e, and Avr9. The largest number of polymorphisms were observed for Avr2. Both mating type genes, MAT1-1-1 and MAT1-2-1, were present in the isolate collection. For further insights into the pathogen diversity, the 50 isolates were genotyped at 7,514 single-nucleotide polymorphism loci using genotyping-by-sequencing. Differentiation by region but not by year was observed. Within the collection of 50 isolates, there were 18 distinct genotypes. Information regarding P. fulva population diversity will enable better management recommendations for growers, as high tunnel production of tomatoes expands.

A chromosome-scale genome assembly of the tomato pathogen Cladosporium fulvum reveals a compartmentalized genome architecture and the presence of a dispensable chromosome

Cladosporium fulvum is a fungal pathogen that causes leaf mould of tomato. The reference genome of this pathogen was released in 2012 but its high repetitive DNA content prevented a contiguous assembly and further prohibited the analysis of its genome architecture. In this study, we combined third generation sequencing technology with the Hi-C chromatin conformation capture technique, to produce a high-quality and near complete genome assembly and gene annotation of a Race 5 isolate of C. fulvum. The resulting genome assembly contained 67.17 Mb organized into 14 chromosomes (Chr1-to-Chr14), all of which were assembled telomere-to-telomere. The smallest of the chromosomes, Chr14, is only 460 kb in size and contains 25 genes that all encode hypothetical proteins. Notably, PCR assays revealed that Chr14 was absent in 19 out of 24 isolates of a world-wide collection of C. fulvum, indicating that Chr14 is dispensable. Thus, C. fulvum is currently the second species of Capnodiales shown to harbour dispensable chromosomes. The genome of C. fulvum Race 5 is 49.7 % repetitive and contains 14 690 predicted genes with an estimated completeness of 98.9%, currently one of the highest among the Capnodiales. Genome structure analysis revealed a compartmentalized architecture composed of gene-dense and repeat-poor regions interspersed with gene-sparse and repeat-rich regions. Nearly 39.2 % of the C. fulvum Race 5 genome is affected by Repeat-Induced Point (RIP) mutations and evidence of RIP leakage toward non-repetitive regions was observed in all chromosomes, indicating the RIP plays an important role in the evolution of this pathogen. Finally, 345 genes encoding candidate effectors were identified in C. fulvum Race 5, with a significant enrichment of their location in gene-sparse regions, in accordance with the 'two-speed genome' model of evolution. Overall, the new reference genome of C. fulvum presents several notable features and is a valuable resource for studies in plant pathogens.

A single amino acid polymorphism in a conserved effector of the multihost blast fungus pathogen expands host-target binding spectrum

Accelerated gene evolution is a hallmark of pathogen adaptation and specialization following host-jumps. However, the molecular processes associated with adaptive evolution between host-specific lineages of a multihost plant pathogen remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), host specialization on different grass hosts is generally associated with dynamic patterns of gain and loss of virulence effector genes that tend to define the distinct genetic lineages of this pathogen. Here, we unravelled the biochemical and structural basis of adaptive evolution of APikL2, an exceptionally conserved paralog of the well-studied rice-lineage specific effector AVR-Pik. Whereas AVR-Pik and other members of the six-gene AVR-Pik family show specific patterns of presence/absence polymorphisms between grass-specific lineages of M. oryzae, APikL2 stands out by being ubiquitously present in all blast fungus lineages from 13 different host species. Using biochemical, biophysical and structural biology methods, we show that a single aspartate to asparagine polymorphism expands the binding spectrum of APikL2 to host proteins of the heavy-metal associated (HMA) domain family. This mutation maps to one of the APikL2-HMA binding interfaces and contributes to an altered hydrogen-bonding network. By combining phylogenetic ancestral reconstruction with an analysis of the structural consequences of allelic diversification, we revealed a common mechanism of effector specialization in the AVR-Pik/APikL2 family that involves two major HMA-binding interfaces. Together, our findings provide a detailed molecular evolution and structural biology framework for diversification and adaptation of a fungal pathogen effector family following host-jumps.

Multiple Mechanisms Drive the Evolutionary Adaptation of Phytophthora infestans Effector Avr1 to Host Resistance

Effectors, a group of small proteins secreted by pathogens, play a central role in antagonistic interactions between plant hosts and pathogens. The evolution of effector genes threatens plant disease management and sustainable food production, but population genetic analyses to understand evolutionary mechanisms of effector genes are limited compared to molecular and functional studies. Here we investigated the evolution of the Avr1 effector gene from 111 Phytophthora infestans isolates collected from six areas covering three potato cropping regions in China using a population genetic approach. High genetic variation of the effector gene resulted from diverse mechanisms including base substitution, pre-termination, intragenic recombination and diversifying selection. Nearly 80% of the 111 sequences had a point mutation in the 512th nucleotide (T512G), which generated a pre-termination stop codon truncating 38 amino acids in the C-terminal, suggesting that the C-terminal may not be essential to ecological and biological functions of P. infestans. A significant correlation between the frequency of Avr1 sequences with the pre-termination and annual mean temperature in the collection sites suggests that thermal heterogeneity might be one of contributors to the diversifying selection, although biological and biochemical mechanisms of the likely thermal adaptation are not known currently. Our results highlight the risk of rapid adaptation of P. infestans and possibly other pathogens as well to host resistance, and the application of eco-evolutionary principles is necessary for sustainable disease management in agricultural ecosystems.

Genome-wide mapping of histone modifications during axenic growth in two species of Leptosphaeria maculans showing contrasting genomic organization

Leptosphaeria maculans 'brassicae' (Lmb) and Leptosphaeria maculans 'lepidii' (Lml) are closely related phytopathogenic species that exhibit a large macrosynteny but contrasting genome structure. Lmb has more than 30% of repeats clustered in large repeat-rich regions, while the Lml genome has only a small amount of evenly distributed repeats. Repeat-rich regions of Lmb are enriched in effector genes, expressed during plant infection. The distinct genome structures of Lmb and Lml provide an excellent model for comparing the organization of pathogenicity genes in relation to the chromatin landscape in two closely related phytopathogenic fungi. Here, we performed chromatin immunoprecipitation (ChIP) during axenic culture, targeting histone modifications typical for heterochromatin or euchromatin, combined with transcriptomic analysis to analyze the influence of chromatin organization on gene expression. In both species, we found that facultative heterochromatin is enriched with genes lacking functional annotation, including numerous effector and species-specific genes. Notably, orthologous genes located in H3K27me3 domains are enriched with effector genes. Compared to other fungal species, including Lml, Lmb is distinct in having large H3K9me3 domains associated with repeat-rich regions that contain numerous species-specific effector genes. Discovery of these two distinctive heterochromatin landscapes now raises questions about their involvement in the regulation of pathogenicity, the dynamics of these domains during plant infection and the selective advantage to the fungus to host effector genes in H3K9me3 or H3K27me3 domains.

Torres D, Usami T, Schouten HJ, Bai Y, Seidl MF, Thomma BPHJ. Comparative genomics reveals the in planta-secreted Verticillium dahliae Av2 effector protein recognized in tomato plants that carry the V2 resistance locus

Plant pathogens secrete effector molecules during host invasion to promote colonization. However, some of these effectors become recognized by host receptors to mount a defence response and establish immunity. Recently, a novel resistance was identified in wild tomato, mediated by the single dominant V2 locus, to control strains of the soil-borne vascular wilt fungus Verticillium dahliae that belong to race 2. With comparative genomics of race 2 strains and resistance-breaking race 3 strains, we identified the avirulence effector that activates V2 resistance, termed Av2. We identified 277 kb of race 2-specific sequence comprising only two genes encoding predicted secreted proteins that are expressed during tomato colonization. Subsequent functional analysis based on genetic complementation into race 3 isolates and targeted deletion from the race 1 isolate JR2 and race 2 isolate TO22 confirmed that one of the two candidates encodes the avirulence effector Av2 that is recognized in V2 tomato plants. Two Av2 allelic variants were identified that encode Av2 variants that differ by a single acid. Thus far, a role in virulence could not be demonstrated for either of the two variants.

Evolution of a guarded decoy protease and its receptor in solanaceous plants

Rcr3 is a secreted protease of tomato that is targeted by fungal effector Avr2, a secreted protease inhibitor of the fungal pathogen Cladosporium fulvum. The Avr2-Rcr3 complex is recognized by receptor-like protein Cf-2, triggering hypersensitive cell death (HR) and disease resistance. Avr2 also targets Rcr3 paralog Pip1, which is not required for Avr2 recognition but contributes to basal resistance. Thus, Rcr3 acts as a guarded decoy in this interaction, trapping the fungus into a recognition event. Here we show that Rcr3 evolved > 50 million years ago (Mya), whereas Cf-2 evolved <6Mya by co-opting the pre-existing Rcr3 in the Solanum genus. Ancient Rcr3 homologs present in tomato, potato, eggplants, pepper, petunia and tobacco can be inhibited by Avr2 with the exception of tobacco Rcr3. Four variant residues in Rcr3 promote Avr2 inhibition, but the Rcr3 that co-evolved with Cf-2 lacks three of these residues, indicating that the Rcr3 co-receptor is suboptimal for Avr2 binding. Pepper Rcr3 triggers HR with Cf-2 and Avr2 when engineered for enhanced inhibition by Avr2. Nicotiana benthamiana (Nb) is a natural null mutant carrying Rcr3 and Pip1 alleles with deleterious frame-shift mutations. Resurrected NbRcr3 and NbPip1 alleles were active proteases and further NbRcr3 engineering facilitated Avr2 inhibition, uncoupled from HR signalling. The evolution of a receptor co-opting a conserved pathogen target contrasts with other indirect pathogen recognition mechanisms.

Allelic Variation and Selection in Effector Genes of Phytophthora infestans (Mont.) de Bary

Phytophthora infestans is a devastating plant pathogen in several crops such as potato (Solanum tuberosum), tomato (Solanum lycopersicum) and Andean fruits such as tree tomato (Solanum betaceum), lulo (Solanum quitoense), uchuva (Physalis peruviana) and wild species in the genus Solanum sp. Despite intense research performed around the world, P. infestans populations from Colombia, South America, are poorly understood. Of particular importance is knowledge about pathogen effector proteins, which are responsible for virulence. The present work was performed with the objective to analyze gene sequences coding for effector proteins of P. infestans from isolates collected from different hosts and geographical regions. Several genetic parameters, phylogenetic analyses and neutrality tests for non-synonymous and synonymous substitutions were calculated. Non-synonymous substitutions were identified for all genes that exhibited polymorphisms at the DNA level. Significant negative selection values were found for two genes (PITG_08994 and PITG_12737) suggesting active coevolution with the corresponding host resistance proteins. Implications for pathogen virulence mechanisms and disease management are discussed.

Divergent Evolution of PcF/SCR74 Effectors in Oomycetes Is Associated with Distinct Recognition Patterns in Solanaceous Plants

Plants deploy cell surface receptors known as pattern-recognition receptors (PRRs) that recognize non-self molecules from pathogens and microbes to defend against invaders. PRRs typically recognize microbe-associated molecular patterns (MAMPs) that are usually widely conserved, some even across kingdoms. Here, we report an oomycete-specific family of small secreted cysteine-rich (SCR) proteins that displays divergent patterns of sequence variation in the Irish potato famine pathogen Phytophthora infestans A subclass that includes the conserved effector PcF from Phytophthora cactorum activates immunity in a wide range of plant species. In contrast, the more diverse SCR74 subclass is specific to P. infestans and tends to trigger immune responses only in a limited number of wild potato genotypes. The SCR74 response was recently mapped to a G-type lectin receptor kinase (G-LecRK) locus in the wild potato Solanum microdontum subsp. gigantophyllum. The G-LecRK locus displays a high diversity in Solanum host species compared to other solanaceous plants. We propose that the diversification of the SCR74 proteins in P. infestans is driven by a fast coevolutionary arms race with cell surface immune receptors in wild potato, which contrasts the presumed slower dynamics between conserved apoplastic effectors and PRRs. Understanding the molecular determinants of plant immune responses to these divergent molecular patterns in oomycetes is expected to contribute to deploying multiple layers of disease resistance in crop plants.IMPORTANCE Immune receptors at the plant cell surface can recognize invading microbes. The perceived microbial molecules are typically widely conserved and therefore the matching surface receptors can detect a broad spectrum of pathogens. Here we describe a family of Phytophthora small extracellular proteins that consists of conserved subfamilies that are widely recognized by solanaceous plants. Remarkably, one subclass of SCR74 proteins is highly diverse, restricted to the late blight pathogen Phytophthora infestans and is specifically detected in wild potato plants. The diversification of this subfamily exhibits signatures of a coevolutionary arms race with surface receptors in potato. Insights into the molecular interaction between these potato-specific receptors and the recognized Phytophthora proteins are expected to contribute to disease resistance breeding in potato.

Genes Encoding Recognition of the Cladosporium fulvum Effector Protein Ecp5 Are Encoded at Several Loci in the Tomato Genome

The molecular interactions between tomato and Cladosporium fulvum have been an important model for molecular plant pathology. Complex genetic loci on tomato chromosomes 1 and 6 harbor genes for resistance to Cladosporium fulvum, encoding receptor like-proteins that perceive distinct Cladosporium fulvum effectors and trigger plant defenses. Here, we report classical mapping strategies for loci in tomato accessions that respond to Cladosporium fulvum effector Ecp5, which is very sequence-monomorphic. We screened 139 wild tomato accessions for an Ecp5-induced hypersensitive response, and in five accessions, the Ecp5-induced hypersensitive response segregated as a monogenic trait, mapping to distinct loci in the tomato genome. We identified at least three loci on chromosomes 1, 7 and 12 that harbor distinct Cf-Ecp5 genes in four different accessions. Our mapping showed that the Cf-Ecp5 in Solanum pimpinellifolium G1.1161 is located at the Milky Way locus. The Cf-Ecp5 in Solanum pimpinellifolium LA0722 was mapped to the bottom arm of chromosome 7, while the Cf-Ecp5 genes in Solanum lycopersicum Ontario 7522 and Solanum pimpinellifolium LA2852 were mapped to the same locus on the top arm of chromosome 12. Bi-parental crosses between accessions carrying distinct Cf-Ecp5 genes revealed putative genetically unlinked suppressors of the Ecp5-induced hypersensitive response. Our mapping also showed that Cf-11 is located on chromosome 11, close to the Cf-3 locus. The Ecp5-induced hypersensitive response is widely distributed within tomato species and is variable in strength. This novel example of convergent evolution could be used for choosing different functional Cf-Ecp5 genes according to individual plant breeding needs.

Time-Course RNAseq Reveals Exserohilum turcicum Effectors and Pathogenicity Determinants

Exserohilum turcicum (sexual stage Setosphaeria turcica) is the hemibiotrophic causal agent of northern leaf blight of maize and sorghum. This study aimed to identify the genes involved in host colonization during the biotrophic and necrotrophic phases of infection. It also aimed to identify race-specific differences in gene expression. RNAseq of maize seedlings inoculated with a race 13N or 23N E. turcicum isolate was conducted before inoculation and at 2, 5, 7, and 13 days post-inoculation (dpi). Biological replicates were pooled per time point for each race and sequenced. A bioinformatics pipeline was used to identify candidate effectors, and expression was validated for selected candidates. Fungal biomass was positively correlated with the percentages of E. turcicum reads mapped, which were low at early time points (2-7 dpi) with a significant increase at 13 dpi, indicating a lifestyle switch from biotrophy to necrotrophy between 7 and 13 dpi. AVRHt1 is the putative E. turcicum effector recognized by the maize resistance gene Ht1. Consistent with this, AVRHt1 was expressed in planta by race 23N, but transcripts were absent in race 13N. In addition, specific transposable elements were expressed in 23N only. Genes encoding the virulence-associated peptidases leupeptin-inhibiting protein 1 and fungalysin were expressed in planta. Transcriptional profiles of genes involved in secondary metabolite synthesis or cell wall degradation revealed the importance of these genes during late stages of infection (13 dpi). A total of 346 expressed candidate effectors were identified, including Ecp6 and proteins similar to the secreted in xylem (SIX) effectors common to formae speciales of Fusarium oxysporum, SIX13 and SIX5. Expression profiling of Ecp6 and SIX13-like indicated a peak in expression at 5 and 7 dpi compared to 2 and 13 dpi. Sequencing of SIX13-like from diverse isolates of E. turcicum revealed host-specific polymorphisms that were mostly non-synonymous, resulting in two groups of SIX13-like proteins that corresponded to the maize or sorghum origin of each isolate. This study suggests putative mechanisms whereby E. turcicum causes disease. Identification of the candidate effector SIX13-like is consistent with the infection mode of E. turcicum through the xylem of susceptible hosts.

DsEcp2-1 is a polymorphic effector that restricts growth of Dothistroma septosporum in pine

The detrimental effect of fungal pathogens on forest trees is an increasingly important problem that has implications for the health of our planet. Despite this, the study of molecular plant-microbe interactions in forest trees is in its infancy, and very little is known about the roles of effector molecules from forest pathogens. Dothistroma septosporum causes a devastating needle blight disease of pines, and intriguingly, is closely related to Cladosporium fulvum, a tomato pathogen in which pioneering effector biology studies have been carried out. Here, we studied D. septosporum effectors that are shared with C. fulvum, by comparing gene sequences from global isolates of D. septosporum and assessing effector function in both host and non-host plants. Many of the effectors were predicted to be non-functional in D. septosporum due to their pseudogenization or low expression in planta, suggesting adaptation to lifestyle and host. Effector sequences were polymorphic among a global collection of D. septosporum isolates, but there was no evidence for positive selection. The DsEcp2-1 effector elicited cell death in the non-host plant Nicotiana tabacum, whilst D. septosporum DsEcp2-1 mutants showed increased colonization of pine needles. Together these results suggest that DsEcp2-1 might be recognized by an immune receptor in both angiosperm and gymnosperm plants. This work may lead to the identification of plant targets for DsEcp2-1 that will provide much needed information on the molecular basis of gymnosperm-pathogen interactions in forests, and may also lead to novel methods of disease control.

Chitin-Binding Protein of Verticillium nonalfalfae Disguises Fungus from Plant Chitinases and Suppresses Chitin-Triggered Host Immunity

During fungal infections, plant cells secrete chitinases, which digest chitin in the fungal cell walls. The recognition of released chitin oligomers via lysin motif (LysM)-containing immune host receptors results in the activation of defense signaling pathways. We report here that Verticillium nonalfalfae, a hemibiotrophic xylem-invading fungus, prevents these digestion and recognition processes by secreting a carbohydrate-binding motif 18 (CBM18)-chitin-binding protein, VnaChtBP, which is transcriptionally activated specifically during the parasitic life stages. VnaChtBP is encoded by the Vna8.213 gene, which is highly conserved within the species, suggesting high evolutionary stability and importance for the fungal lifestyle. In a pathogenicity assay, however, Vna8.213 knockout mutants exhibited wilting symptoms similar to the wild-type fungus, suggesting that Vna8.213 activity is functionally redundant during fungal infection of hop. In a binding assay, recombinant VnaChtBP bound chitin and chitin oligomers in vitro with submicromolar affinity and protected fungal hyphae from degradation by plant chitinases. Moreover, the chitin-triggered production of reactive oxygen species from hop suspension cells was abolished in the presence of VnaChtBP, indicating that VnaChtBP also acts as a suppressor of chitin-triggered immunity. Using a yeast-two-hybrid assay, circular dichroism, homology modeling, and molecular docking, we demonstrated that VnaChtBP forms dimers in the absence of ligands and that this interaction is stabilized by the binding of chitin hexamers with a similar preference in the two binding sites. Our data suggest that, in addition to chitin-binding LysM (CBM50) and Avr4 (CBM14) fungal effectors, structurally unrelated CBM18 effectors have convergently evolved to prevent hydrolysis of the fungal cell wall against plant chitinases and to interfere with chitin-triggered host immunity.

Molecular Interactions Between Smut Fungi and Their Host Plants

Smut fungi are a large group of biotrophic plant pathogens that infect mostly monocot species, including economically relevant cereal crops. For years, Ustilago maydis has stood out as the model system to study the genetics and cell biology of smut fungi as well as the pathogenic development of biotrophic plant pathogens. The identification and functional characterization of secreted effectors and their role in virulence have particularly been driven forward using the U. maydis-maize pathosystem. Today, advancing tools for additional smut fungi such as Ustilago hordei and Sporisorium reilianum, as well as an increasing number of available genome sequences, provide excellent opportunities to investigate in parallel the effector function and evolution associated with different lifestyles and host specificities. In addition, genome analyses revealed similarities in the genomic signature between pathogenic smuts and epiphytic Pseudozyma species. This review elaborates on how knowledge about fungal lifestyles, genome biology, and functional effector biology has helped in understanding the biology of this important group of fungal pathogens. We highlight the contribution of the U. maydis model system but also discuss the differences from other smut fungi, which raises the importance of comparative genomic and genetic analyses in future research.

Effector proteins of Rhizophagus proliferus: conserved protein domains may play a role in host-specific interaction with different plant species

Arbuscular mycorrhizal (AM) fungi show high promiscuity in terms of host. Effector proteins expressed by AM fungi are found important in establishing interaction with host. However, the mechanistic underlying host-specific interactions of the fungi remain unknown. The present study aimed (i) to identify effectors encoded by Rhizophagus proliferus and (ii) to understand molecular specificity encoded in effectors for interaction with specific plant species. The effectors predicted from the whole genome sequence were annotated by homology search in NCBI non-redundant protein, Interproscan, and pathogen-host interaction (PHI) databases. In total, 416 small secreted peptides (SSPs) were predicted, which were effector peptides with presence of nuclear localization signal, small cysteine-rich, and repeat-containing proteins domains. Similar to the functionally validated SP7 effectors in Rhizophagus irregularis, two proteins (RP8598 and RP23081) were identified in R. proliferus. To understand whether interaction between SP7 and the plant target protein, ERF19, is specific in nature, we examined protein-peptide interaction using in silico molecular docking. Pairwise interaction of RP8598 and RP23081 with the ethylene-responsive factors (ERF19) coded by five different plant species (Lotus japonicus, Solanum lycopersicum, Ocimum tenuiflorum, Medicago truncatula, Diospyros kaki) was investigated. Prediction of high-quality interaction of SP7 effector with ERF19 protein expressed only by specific plant species was observed in in silico molecular docking, which may reiterate the role of effectors in host specificity. The outcomes from our study indicated that sequence precision encoded in the effector peptides of AM fungi and immunomodulatory proteins of host may regulate host specificity in these fungi.

Dynamic Insertion of Pot3 in AvrPib Prevailing in a Field Rice Blast Population in the Philippines Led to the High Virulence Frequency Against the Resistance Gene Pib in Rice

The Magnaporthe oryzae avirulence gene AvrPib is required for the resistance mediated by its cognate resistance gene Pib, which has been intensively used in indica rice breeding programs in many Asian countries. However, the sequence diversity of AvrPib among geographically distinct M. oryzae populations was recently shown to be increasing. Here, we selected a field population consisting of 248 rice blast isolates collected from a disease hotspot in Philippine for the analysis of AvrPib haplotypes and their pathogenicity against Pib. We found that all of the isolates were virulent to Pib and each of them contained an insertion of Pot3 transposon in AvrPib. Moreover, Pot3 insertion was detected in different genomic positions, resulting in three different AvrPib haplotypes, designated avrPib-H1 to H3. We further conducted a genome-wide Pot2 fingerprinting analysis by repetitive element palindromic polymerase chain reaction (PCR) and identified seven different lineages out of 47 representative isolates. The isolates belonging to the same lineage often had the same AvrPib haplotype. In contrast, the isolates having the same AvrPib haplotypes did not always belong to the same lineages. Both mating types MAT1-1 and MAT1-2 were identified in the population in Bohol and the latter appeared dominant. On the host side, we found that 32 of 52 released rice varieties in the Philippines contained Pib diagnosed by PCR gene-specific primers and DNA sequencing of gene amplicons, suggesting that it was widely incorporated in different rice varieties. Our study highlights the genetic dynamics of rice blast population at both the AvrPib locus and the genome-wide levels, providing insight into the mechanisms of the mutations in AvrPib leading to the breakdown of Pib-mediated resistance in rice.

Specific Hypersensitive Response-Associated Recognition of New Apoplastic Effectors from Cladosporium fulvum in Wild Tomato

Tomato leaf mold disease is caused by the biotrophic fungus Cladosporium fulvum. During infection, C. fulvum produces extracellular small secreted protein (SSP) effectors that function to promote colonization of the leaf apoplast. Resistance to the disease is governed by Cf immune receptor genes that encode receptor-like proteins (RLPs). These RLPs recognize specific SSP effectors to initiate a hypersensitive response (HR) that renders the pathogen avirulent. C. fulvum strains capable of overcoming one or more of all cloned Cf genes have now emerged. To combat these strains, new Cf genes are required. An effectoromics approach was employed to identify wild tomato accessions carrying new Cf genes. Proteomics and transcriptome sequencing were first used to identify 70 apoplastic in planta-induced C. fulvum SSPs. Based on sequence homology, 61 of these SSPs were novel or lacked known functional domains. Seven, however, had predicted structural homology to antimicrobial proteins, suggesting a possible role in mediating antagonistic microbe-microbe interactions in planta. Wild tomato accessions were then screened for HR-associated recognition of 41 SSPs, using the Potato virus X-based transient expression system. Nine SSPs were recognized by one or more accessions, suggesting that these plants carry new Cf genes available for incorporation into cultivated tomato.

Distinct regions of the Phytophthora essential effector Avh238 determine its function in cell death activation and plant immunity suppression

Phytophthora pathogens secrete effectors to manipulate host innate immunity, thus facilitating infection. Among the RXLR effectors highly induced during Phytophthora sojae infection, Avh238 not only contributes to pathogen virulence but also triggers plant cell death. However, the detailed molecular basis of Avh238 functions remains largely unknown. We mapped the regions responsible for Avh238 functions in pathogen virulence and plant cell death induction using a strategy that combines investigation of natural variation and large-scale mutagenesis assays. The correlation between cellular localization and Avh238 functions was also evaluated. We found that the 79^th residue (histidine or leucine) of Avh238 determined its cell death-inducing activity, and that the 53 amino acids in its C-terminal region are responsible for promoting Phytophthora infection. Transient expression of Avh238 in Nicotiana benthamiana revealed that nuclear localization is essential for triggering cell death, while Avh238-mediated suppression of INF1-triggered cell death requires cytoplasmic localization. Our results demonstrate that a representative example of an essential Phytophthora RXLR effector can evolve to escape recognition by the host by mutating one nucleotide site, and can also retain plant immunosuppressive activity to enhance pathogen virulence in planta.

Comparative genomics of Fusarium oxysporum f. sp. melonis reveals the secreted protein recognized by the Fom-2 resistance gene in melon

Development of resistant crops is the most effective way to control plant diseases to safeguard food and feed production. Disease resistance is commonly based on resistance genes, which generally mediate the recognition of small proteins secreted by invading pathogens. These proteins secreted by pathogens are called 'avirulence' proteins. Their identification is important for being able to assess the usefulness and durability of resistance genes in agricultural settings. We have used genome sequencing of a set of strains of the melon wilt fungus Fusarium oxysporum f. sp. melonis (Fom), bioinformatics-based genome comparison and genetic transformation of the fungus to identify AVRFOM2, the gene that encodes the avirulence protein recognized by the melon Fom-2 gene. Both an unbiased and a candidate gene approach identified a single candidate for the AVRFOM2 gene. Genetic complementation of AVRFOM2 in three different race 2 isolates resulted in resistance of Fom-2-harbouring melon cultivars. AvrFom2 is a small, secreted protein with two cysteine residues and weak similarity to secreted proteins of other fungi. The identification of AVRFOM2 will not only be helpful to select melon cultivars to avoid melon Fusarium wilt, but also to monitor how quickly a Fom population can adapt to deployment of Fom-2-containing cultivars in the field.

A conserved proline residue in Dothideomycete Avr4 effector proteins is required to trigger a Cf-4-dependent hypersensitive response

CfAvr4, a chitin-binding effector protein produced by the Dothideomycete tomato pathogen Cladosporium fulvum, protects the cell wall of this fungus against hydrolysis by secreted host chitinases during infection. However, in the presence of the Cf-4 immune receptor of tomato, CfAvr4 triggers a hypersensitive response (HR), which renders the pathogen avirulent. Recently, several orthologues of CfAvr4 have been identified from phylogenetically closely related species of Dothideomycete fungi. Of these, DsAvr4 from Dothistroma septosporum also triggers a Cf-4-dependent HR, but CaAvr4 and CbAvr4 from Cercospora apii and Cercospora beticola, respectively, do not. All, however, bind chitin. To identify the region(s) and specific amino acid residue(s) of CfAvr4 and DsAvr4 required to trigger a Cf-4-dependent HR, chimeric and mutant proteins, in which specific protein regions or single amino acid residues, respectively, were exchanged between CfAvr4 and CaAvr4 or DsAvr4 and CbAvr4, were tested for their ability to trigger an HR in Nicotiana benthamiana plants transgenic for the Cf-4 immune receptor gene. Based on this approach, a single region common to CfAvr4 and DsAvr4 was determined to carry a conserved proline residue necessary for the elicitation of this HR. In support of this result, a Cf-4-dependent HR was triggered by mutant CaAvr4 and CbAvr4 proteins carrying an arginine-to-proline substitution at this position. This study provides the first step in deciphering how Avr4 orthologues from different Dothideomycete fungi trigger a Cf-4-dependent HR.

A Survey on Occurrence of Cladosporium fulvum Identifies Race 0 and Race 2 in Tomato-Growing Areas of Argentina

The presence of Cladosporium fulvum (syn. Passalora fulva), causal agent of tomato leaf mold, was confirmed in the two main greenhouse-production areas for tomato in Argentina. Using both morphological characters and internal transcribed spacer sequencing, we confirmed the presence of physiological races of this pathogen. A diagnostic multiplex polymerase chain reaction (PCR) was also developed, using primers derived from C. fulvum avirulence (Avr) genes. In all, 20 isolates of Cladosporium spp. were obtained as monospore cultures and 12 were identified as C. fulvum. By this method, we showed that, of these 12 isolates, 5 were race 0 (carrying functional Avr2, Avr4, Avr4E, and Avr9 genes) and 7 were race 2 (lacking the Avr2 gene). Race identity was confirmed by testing their virulence on a set of tomato differentials carrying different Cf resistance genes. All Avr genes could be amplified in single or multiplex PCR using DNA isolated from in vitro grown monospore cultures but only three Avr could be amplified when genomic DNA was isolated from C. fulvum-infected necrotic leaf tissue.

Rapidly Evolving Genes Are Key Players in Host Specialization and Virulence of the Fungal Wheat Pathogen Zymoseptoria tritici (Mycosphaerella graminicola)

The speciation of pathogens can be driven by divergent host specialization. Specialization to a new host is possible via the acquisition of advantageous mutations fixed by positive selection. Comparative genome analyses of closely related species allows for the identification of such key substitutions via inference of genome-wide signatures of positive selection. We previously used a comparative genomics framework to identify genes that have evolved under positive selection during speciation of the prominent wheat pathogen Zymoseptoria tritici (synonym Mycosphaerella graminicola). In this study, we conducted functional analyses of four genes exhibiting strong signatures of positive selection in Z. tritici. We deleted the four genes in Z. tritici and confirm a virulence-related role of three of the four genes ΔZt80707, ΔZt89160 and ΔZt103264. The two mutants ΔZt80707 and ΔZt103264 show a significant reduction in virulence during infection of wheat; the ΔZt89160 mutant causes a hypervirulent phenotype in wheat. Mutant phenotypes of ΔZt80707, ΔZt89160 and ΔZt103264 can be restored by insertion of the wild-type genes. However, the insertion of the Zt80707 and Zt89160 orthologs from Z. pseudotritici and Z. ardabiliae do not restore wild-type levels of virulence, suggesting that positively selected substitutions in Z. tritici may relate to divergent host specialization. Interestingly, the gene Zt80707 encodes also a secretion signal that targets the protein for cell secretion. This secretion signal is however only transcribed in Z. tritici, suggesting that Z. tritici-specific substitutions relate to a new function of the protein in the extracellular space of the wheat-Z. tritici interaction. Together, the results presented here highlight that Zt80707, Zt103264 and Zt89160 represent key genes involved in virulence and host-specific disease development of Z. tritici. Our findings illustrate that evolutionary predictions provide a powerful tool for the identification of novel traits crucial for host adaptation and pathogen evolution.

Zymoseptoria spp provides a unique model system to study the underlying genetics of host specialization of plant pathogens. Closely related Zymoseptoria species, including the prominent wheat pathogen Z. tritici, have recently specialized to distinct grass hosts. Positively selected substitutions have played a central role in the acquisition of new host specificities. We have identified a small set of genes showing signatures of positive selection. We demonstrate that three of these four candidate genes play an important role during host infection. Two mutants of Z. tritici were impaired in virulence; a third mutant showed a hypervirulent phenotype. New protein specificities not only include changes at the amino acid sequence level but also at the level of the protein structure. We conducted a gene replacement experiment to test if mutant phenotypes in Z. tritici could be complemented by the insertion of orthologous genes from the two closely related species Z. pseudotritici and Z. ardabiliae. For two genes, we confirm that the species-specific protein changes are essential for proper protein functioning in Z. tritici; key traits involved in the evolution of virulence and host specificity of this prominent pathogen can be characterized via a combination of evolutionary predictions and functional analyses.

The fungal core effector Pep1 is conserved across smuts of dicots and monocots

The secreted fungal effector Pep1 is essential for penetration of the host epidermis and establishment of biotrophy in the Ustilago maydis-maize pathosystem. Previously, Pep1 was found to be an inhibitor of apoplastic plant peroxidases, which suppresses the oxidative burst, a primary immune response of the host plant and enables fungal colonization. To investigate the conservation of Pep1 in other pathogens, genomes of related smut species were screened for pep1 orthologues. Pep1 proteins were produced in Escherichia coli for functional assays. The biological function of Pep1 was tested by heterologous expression in U. maydis and Hordeum vulgare. Pep1 orthologues revealed a remarkable degree of sequence conservation, indicating that this effector might play a fundamental role in virulence of biotrophic smut fungi. Pep1 function and its role in virulence are conserved in different pathogenic fungi, even across the monocot-dicot border of host plants. The findings described in this study classify Pep1 as a phylogenetically conserved fungal core effector. Furthermore, we documented the influence of Pep1 on the disease caused by Blumeria graminis f. sp. hordei which is a non-smut-related pathosystem.

Novel Mutations Detected in Avirulence Genes Overcoming Tomato Cf Resistance Genes in Isolates of a Japanese Population of Cladosporium fulvum

Leaf mold of tomato is caused by the biotrophic fungus Cladosporium fulvum which complies with the gene-for-gene system. The disease was first reported in Japan in the 1920s and has since been frequently observed. Initially only race 0 isolates were reported, but since the consecutive introduction of resistance genes Cf-2, Cf-4, Cf-5 and Cf-9 new races have evolved. Here we first determined the virulence spectrum of 133 C. fulvum isolates collected from 22 prefectures in Japan, and subsequently sequenced the avirulence (Avr) genes Avr2, Avr4, Avr4E, Avr5 and Avr9 to determine the molecular basis of overcoming Cf genes. Twelve races of C. fulvum with a different virulence spectrum were identified, of which races 9, 2.9, 4.9, 4.5.9 and 4.9.11 occur only in Japan. The Avr genes in many of these races contain unique mutations not observed in races identified elsewhere in the world including (i) frameshift mutations and (ii) transposon insertions in Avr2, (iii) point mutations in Avr4 and Avr4E, and (iv) deletions of Avr4E, Avr5 and Avr9. New races have developed by selection pressure imposed by consecutive introductions of Cf-2, Cf-4, Cf-5 and Cf-9 genes in commercially grown tomato cultivars. Our study shows that molecular variations to adapt to different Cf genes in an isolated C. fulvum population in Japan are novel but overall follow similar patterns as those observed in populations from other parts of the world. Implications for breeding of more durable C. fulvum resistant varieties are discussed.

Repeat-containing protein effectors of plant-associated organisms

Many plant-associated organisms, including microbes, nematodes, and insects, deliver effector proteins into the apoplast, vascular tissue, or cell cytoplasm of their prospective hosts. These effectors function to promote colonization, typically by altering host physiology or by modulating host immune responses. The same effectors however, can also trigger host immunity in the presence of cognate host immune receptor proteins, and thus prevent colonization. To circumvent effector-triggered immunity, or to further enhance host colonization, plant-associated organisms often rely on adaptive effector evolution. In recent years, it has become increasingly apparent that several effectors of plant-associated organisms are repeat-containing proteins (RCPs) that carry tandem or non-tandem arrays of an amino acid sequence or structural motif. In this review, we highlight the diverse roles that these repeat domains play in RCP effector function. We also draw attention to the potential role of these repeat domains in adaptive evolution with regards to RCP effector function and the evasion of effector-triggered immunity. The aim of this review is to increase the profile of RCP effectors from plant-associated organisms.

Transcriptome sequencing uncovers the Avr5 avirulence gene of the tomato leaf mold pathogen Cladosporium fulvum

The Cf-5 gene of tomato confers resistance to strains of the fungal pathogen Cladosporium fulvum carrying the avirulence gene Avr5. Although Cf-5 has been cloned, Avr5 has remained elusive. We report the cloning of Avr5 using a combined bioinformatic and transcriptome sequencing approach. RNA-Seq was performed on the sequenced race 0 strain (0WU; carrying Avr5), as well as a race 5 strain (IPO 1979; lacking a functional Avr5 gene) during infection of susceptible tomato. Forty-four in planta-induced C. fulvum candidate effector (CfCE) genes of 0WU were identified that putatively encode a secreted, small cysteine-rich protein. An expressed transcript sequence comparison between strains revealed two polymorphic CfCE genes in IPO 1979. One of these conferred avirulence to IPO 1979 on Cf-5 tomato following complementation with the corresponding 0WU allele, confirming identification of Avr5. Complementation also led to increased fungal biomass during infection of susceptible tomato, signifying a role for Avr5 in virulence. Seven of eight race 5 strains investigated escape Cf-5-mediated resistance through deletion of the Avr5 gene. Avr5 is heavily flanked by repetitive elements, suggesting that repeat instability, in combination with Cf-5-mediated selection pressure, has led to the emergence of race 5 strains deleted for the Avr5 gene.

Positive selection and intragenic recombination contribute to high allelic diversity in effector genes of Mycosphaerella fijiensis, causal agent of the black leaf streak disease of banana

Previously, we have determined the nonhost-mediated recognition of the MfAvr4 and MfEcp2 effector proteins from the banana pathogen Mycosphaerella fijiensis in tomato, by the cognate Cf-4 and Cf-Ecp2 resistance proteins, respectively. These two resistance proteins could thus mediate resistance against M. fijiensis if genetically transformed into banana (Musa spp.). However, disease resistance controlled by single dominant genes can be overcome by mutated effector alleles, whose products are not recognized by the cognate resistance proteins. Here, we surveyed the allelic variation within the MfAvr4, MfEcp2, MfEcp2-2 and MfEcp2-3 effector genes of M. fijiensis in a global population of the pathogen, and assayed its impact on recognition by the tomato Cf-4 and Cf-Ecp2 resistance proteins, respectively. We identified a large number of polymorphisms that could reflect a co-evolutionary arms race between host and pathogen. The analysis of nucleotide substitution patterns suggests that both positive selection and intragenic recombination have shaped the evolution of M. fijiensis effectors. Clear differences in allelic diversity were observed between strains originating from South-East Asia relative to strains from other banana-producing continents, consistent with the hypothesis that M. fijiensis originated in the Asian-Pacific region. Furthermore, transient co-expression of the MfAvr4 effector alleles and the tomato Cf-4 resistance gene, as well as of MfEcp2, MfEcp2-2 and MfEcp2-3 and the putative Cf-Ecp2 resistance gene, indicated that effector alleles able to overcome these resistance genes are already present in natural populations of the pathogen, thus questioning the durability of resistance that can be provided by these genes in the field.

Secretome analysis reveals effector candidates associated with broad host range necrotrophy in the fungal plant pathogen Sclerotinia sclerotiorum

BACKGROUND

The white mold fungus Sclerotinia sclerotiorum is a devastating necrotrophic plant pathogen with a remarkably broad host range. The interaction of necrotrophs with their hosts is more complex than initially thought, and still poorly understood.

RESULTS

We combined bioinformatics approaches to determine the repertoire of S. sclerotiorum effector candidates and conducted detailed sequence and expression analyses on selected candidates. We identified 486 S. sclerotiorum secreted protein genes expressed in planta, many of which have no predicted enzymatic activity and may be involved in the interaction between the fungus and its hosts. We focused on those showing (i) protein domains and motifs found in known fungal effectors, (ii) signatures of positive selection, (iii) recent gene duplication, or (iv) being S. sclerotiorum-specific. We identified 78 effector candidates based on these properties. We analyzed the expression pattern of 16 representative effector candidate genes on four host plants and revealed diverse expression patterns.

CONCLUSIONS

These results reveal diverse predicted functions and expression patterns in the repertoire of S. sclerotiorum effector candidates. They will facilitate the functional analysis of fungal pathogenicity determinants and should prove useful in the search for plant quantitative disease resistance components active against the white mold.

[Transposable elements reshaping genomes and favouring the evolutionary and adaptive potential of fungal phytopathogens]

Phytopathogenic fungi are a major threat for global food security and show an extreme plasticity in pathogenicity behaviours. They often have a high adaptive potential allowing them to rapidly counteract the control methods used by men in agrosystems. In this paper, we evaluate the link between genome plasticity and adaptive potential using genomics and comparative genomics approaches. Our model is the evolutionary series Leptosphaeria maculans-Leptosphaeria biglobosa, encompassing five distinct entities, whose conspecificity or heterospecificity status is unclear, and which all are pathogens of cruciferous plants. They however differ by their host range and pathogenicity. Compared to other species of the species complex, the species best adapted to oilseed rape, L. maculans "brassicae", causing important losses in the crop, has a genome that was submitted to a recent and massive burst of transposition by a few families of transposable elements (TEs). Whether the genome invasion contributed to speciation is still unclear to-date but there is a coincidence between this burst of TEs and divergence between two species. This TE burst contributed to diversification of effector proteins and thus to generation of novel pathogenic specificities. In addition, the location of effector genes within genome regions enriched in TEs has direct consequences on adaptation to plant resistance and favours a multiplicity of mutation events allowing "breakdown" of resistance. These data are substantiated by other examples in the literature showing that fungi tend to have a "two-speed" genome, in which a plastic compartment enriched in TE host genes is involved in pathogenicity and adaptation to host.

Pseudogenization in pathogenic fungi with different host plants and lifestyles might reflect their evolutionary past

Pseudogenes are genes with significant homology to functional genes, but contain disruptive mutations (DMs) leading to the production of non- or partially functional proteins. Little is known about pseudogenization in pathogenic fungi with different lifestyles. Here, we report the identification of DMs causing pseudogenes in the genomes of the fungal plant pathogens Botrytis cinerea, Cladosporium fulvum, Dothistroma septosporum, Mycosphaerella fijiensis, Verticillium dahliae and Zymoseptoria tritici. In these fungi, we identified 1740 gene models containing 2795 DMs obtained by an alignment-based gene prediction method. The contribution of sequencing errors to DMs was minimized by analyses of resequenced genomes to obtain a refined dataset of 924 gene models containing 1666 true DMs. The frequency of pseudogenes varied from 1% to 5% in the gene catalogues of these fungi, being the highest in the asexually reproducing fungus C. fulvum (4.9%), followed by D. septosporum (2.4%) and V. dahliae (2.1%). The majority of pseudogenes do not represent recent gene duplications, but members of multi-gene families and unitary genes. In general, there was no bias for pseudogenization of specific genes in the six fungi. Single exceptions were those encoding secreted proteins, including proteases, which appeared more frequently pseudogenized in C. fulvum than in D. septosporum. Most pseudogenes present in these two phylogenetically closely related fungi are not shared, suggesting that they are related to adaptation to a different host (tomato versus pine) and lifestyle (biotroph versus hemibiotroph).

When virulence originates from non-agricultural hosts: new insights into plant breeding

Monogenic plant resistance breakdown is a model for testing evolution in action in pathogens. As a rule, plant pathologists argue that virulence - the allele that allows pathogens to overcome resistance - is due to a new mutation at the avirulence locus within the native/endemic population that infects susceptible crops. In this article, we develop an alternative and neglected scenario where a given virulence pre-exists in a non-agricultural host and might be accidentally released or introduced on the matching resistant cultivar in the field. The main difference between the two scenarios is the divergence time expected between the avirulent and the virulent populations. As a consequence, population genetic approaches such as genome scans and Approximate Bayesian Computation methods allow explicit testing of the two scenarios by timing the divergence. This review then explores the fundamental implications of this alternative scenario for plant breeding, including the invasion of virulence or the evolution of more aggressive hybrids, and proposes concrete solutions to achieve a sustainable resistance.

The Phytophthora sojae Avr1d gene encodes an RxLR-dEER effector with presence and absence polymorphisms among pathogen strains

Soybean root and stem rot is caused by the oomycete pathogen Phytophthora sojae. The interaction between P. sojae and soybean fits the "gene-for-gene" hypothesis. Although more than 10 P. sojae avirulence (Avr) effectors have been genetically identified, nearly half of genetically defined avr genes have been cloned. In a previous bioinformatic and global transcriptional analysis, we identified a P. sojae RxLR effector, Avr1d, which was 125 amino acids in length. Mapping data demonstrated that Avr1d presence or absence in the genome was co-segregated with the Avr1d avirulence phenotype in F2 populations. Transient expression of the Avr1d gene using co-bombardment in soybean isogenic lines revealed that this gene triggered a hypersensitive response (HR) in the presence of Rps1d. Sequencing of Avr1d genes in different P. sojae strains revealed two Avr1d alleles. Although polymorphic, the two Avr1d alleles could trigger Rps1d-mediated HR. P. sojae strains carrying either of the alleles were avirulent on Rps1d soybean lines. Avr1d was upregulated during the germinating cyst and early infection stages. Furthermore, transient expression of Avr1d in Nicotiana benthamiana suppressed BAX-induced cell death and enhanced P. capsici infection. Avr1d also suppressed effector-triggered immunity induction by associating with Avr1b and Rps1b, suggestive of a role in suppressing plant immunity.

Genome structure and reproductive behaviour influence the evolutionary potential of a fungal phytopathogen

Modern agriculture favours the selection and spread of novel plant diseases. Furthermore, crop genetic resistance against pathogens is often rendered ineffective within a few years of its commercial deployment. Leptosphaeria maculans, the cause of phoma stem canker of oilseed rape, develops gene-for-gene interactions with its host plant, and has a high evolutionary potential to render ineffective novel sources of resistance in crops. Here, we established a four-year field experiment to monitor the evolution of populations confronted with the newly released Rlm7 resistance and to investigate the nature of the mutations responsible for virulence against Rlm7. A total of 2551 fungal isolates were collected from experimental crops of a Rlm7 cultivar or a cultivar without Rlm7. All isolates were phenotyped for virulence and a subset was genotyped with neutral genetic markers. Virulent isolates were investigated for molecular events at the AvrLm4-7 locus. Whilst virulent isolates were not found in neighbouring crops, their frequency had reached 36% in the experimental field after four years. An extreme diversity of independent molecular events leading to virulence was identified in populations, with large-scale Repeat Induced Point mutations or complete deletion of AvrLm4-7 being the most frequent. Our data suggest that increased mutability of fungal genes involved in the interactions with plants is directly related to their genomic environment and reproductive system. Thus, rapid allelic diversification of avirulence genes can be generated in L. maculans populations in a single field provided that large population sizes and sexual reproduction are favoured by agricultural practices.

Phytotoxic secondary metabolites and peptides produced by plant pathogenic Dothideomycete fungi

Many necrotrophic plant pathogenic fungi belonging to the class of Dothideomycetes produce phytotoxic metabolites and peptides that are usually required for pathogenicity. Phytotoxins that affect a broad range of plant species are known as non-host-specific toxins (non-HSTs), whereas HSTs affect only a particular plant species or more often genotypes of that species. For pathogens producing HSTs, pathogenicity and host specificity are largely defined by the ability to produce the toxin, while plant susceptibility is dependent on the presence of the toxin target. Non-HSTs are not the main determinants of pathogenicity but contribute to virulence of the producing pathogen. Dothideomycetes are remarkable for the production of toxins, particularly HSTs because they are the only fungal species known so far to produce them. The synthesis, regulation, and mechanisms of action of the most important HSTs and non-HSTs will be discussed. Studies on the mode of action of HSTs have highlighted the induction of programed cell death (PCD) as an important mechanism. We discuss HST-induced PCD and the plant hypersensitive response upon recognition of avirulence factors that share common pathways. In this respect, although nucleotide-binding-site-leucine-rich repeat types of resistance proteins mediate resistance against biotrophs, they can also contribute to susceptibility toward necrotrophs.

Estimating frequencies of virulent isolates in field populations of a plant pathogenic fungus, Leptosphaeria maculans, using high-throughput pyrosequencing

AIM

To develop a pyrosequencing assay to monitor the frequency of alleles of an avirulence gene, AvrLm4, in populations of sexual spores of Leptosphaeria maculans, a fungal pathogen of canola (Brassica napus).

METHODS AND RESULTS

The predominant mutation in AvrLm4 responsible for virulence to the corresponding resistance gene, Rlm4, is a single nucleotide polymorphism (SNP) at base 358. Pyrosequencing primers were designed to amplify a 90-bp region that included this SNP. The assay was developed and validated by analysing the frequency of AvrLm4 in isolate mixtures of different proportions. Furthermore, the frequency of avrLm4 (virulence allele) determined by pyrosequencing of populations of sexual spores was consistent with the frequency of avrLm4 determined by Sanger sequencing of the entire AvrLm4 gene from single isolates cultured from the same stubble.

CONCLUSION

This high-throughput assay can play an important role in predicting the risk of resistance breakdown in crops.

SIGNIFICANCE AND IMPACT OF THE STUDY

Similar assays can be applied to monitor frequencies of fungicide resistance in pathogens of crops and to assay diversity in microbial soil communities such as in soil samples from bat caves where white-nose syndrome has been detected.

In silico characterization and molecular evolutionary analysis of a novel superfamily of fungal effector proteins

Most fungal plant pathogens secrete effector proteins during pathogenesis to manipulate their host's defense and promote disease. These are so highly diverse in sequence and distribution, they are essentially considered as species-specific. However, we have recently shown the presence of homologous effectors in fungal species of the Dothideomycetes class. One such example is Ecp2, an effector originally described in the tomato pathogen Cladosporium fulvum but later detected in the plant pathogenic fungi Mycosphaerella fijiensis and Mycosphaerella graminicola as well. Here, using in silico sequence-similarity searches against a database of 135 fungal genomes and GenBank, we extend our queries for homologs of Ecp2 to the fungal kingdom and beyond, and further study their history of diversification. Our analyses show that Ecp2 homologs are members of an ancient and widely distributed superfamily of putative fungal effectors, which we term Hce2 for Homologs of C. fulvum Ecp2. Molecular evolutionary analyses show that the superfamily originated and diversified within the fungal kingdom, experiencing multiple lineage-specific expansions and losses that are consistent with the birth-and-death model of gene family evolution. Newly formed paralogs appear to be subject to diversification early after gene duplication events, whereas at later stages purifying selection acts to preserve diversity and the newly evolved putative functions. Some members of the Hce2 superfamily are fused to fungal Glycoside Hydrolase family 18 chitinases that show high similarity to the Zymocin killer toxin from the dairy yeast Kluyveromyces lactis, suggesting an analogous role in antagonistic interactions. The observed high rates of gene duplication and loss in the Hce2 superfamily, combined with diversification in both sequence and possibly functions within and between species, suggest that Hce2s are involved in adaptation to stresses and new ecological niches. Such findings address the need to rationalize effector biology and evolution beyond the perspective of solely host-microbe interactions.

Quantitative variation in effector activity of ToxA isoforms from Stagonospora nodorum and Pyrenophora tritici-repentis

ToxA is a proteinaceous necrotrophic effector produced by Stagonospora nodorum and Pyrenophora tritici-repentis. In this study, all eight mature isoforms of the ToxA protein were purified and compared. Circular dichroism spectra indicated that all isoforms were structurally intact and had indistinguishable secondary structural features. ToxA isoforms were infiltrated into wheat lines that carry the sensitivity gene Tsn1. It was observed that different wheat lines carrying identical Tsn1 alleles varied in sensitivity to ToxA. All ToxA isoforms induced necrosis when introduced into any Tsn1 wheat line but we observed quantitative variation in effector activity, with the least-active version found in isolates of P. tritici-repentis. Pathogen sporulation increased with higher doses of ToxA. The isoforms that induced the most rapid necrosis also induced the most sporulation, indicating that pathogen fitness is affected by differences in ToxA activity. We show that differences in toxin activity encoded by a single gene can contribute to the quantitative inheritance of necrotrophic virulence. Our findings support the hypothesis that the variation at ToxA results from selection that favors increased toxin activity.

Analysis of polymorphism and transcription of the effector gene Avr1b in Phytophthora sojae isolates from China virulent to Rps1b

The effector gene Avr1b-1 of Phytophthora sojae determines the efficacy of the resistance gene Rps1b in soybean. The sequences of the Avr1b-1 locus in 34 Chinese isolates of P. sojae were obtained and analysed by polymerase chain reaction (PCR) and inverse PCR. Four different alleles and a complete deletion mutation of the Avr1b-1 gene were identified. Molecular analysis of the deletion breakpoints in the Avr1b-1 locus revealed that an 8-kb DNA sequence containing Avr1b-1 was deleted and a 12.7-kb DNA sequence was inserted at the same locus. A truncated transposase gene was found and five transposable elements were predicted in the inserted sequence, suggesting that the deletion of Avr1b-1 might be attributed to transposon movement. The transcription of Avr1b-1 was analysed in virulent isolates containing four alleles of Avr1b-1 by real-time reverse transcription-PCR. In all virulent isolates, only those isolates containing the second allele transcripted Avr1b-1.

Identification of Hyaloperonospora arabidopsidis transcript sequences expressed during infection reveals isolate-specific effectors

Biotrophic plant pathogens secrete effector proteins that are important for infection of the host. The aim of this study was to identify effectors of the downy mildew pathogen Hyaloperonospora arabidopsidis (Hpa) that are expressed during infection of its natural host Arabidopsis thaliana. Infection-related transcripts were identified from Expressed Sequence Tags (ESTs) derived from leaves of the susceptible Arabidopsis Ws eds1-1 mutant inoculated with the highly virulent Hpa isolate Waco9. Assembly of 6364 ESTs yielded 3729 unigenes, of which 2164 were Hpa-derived. From the translated Hpa unigenes, 198 predicted secreted proteins were identified. Of these, 75 were found to be Hpa-specific and six isolate Waco9-specific. Among 42 putative effectors identified there were three Elicitin-like proteins, 16 Cysteine-rich proteins and 18 host-translocated RXLR effectors. Sequencing of alleles in different Hpa isolates revealed that five RXLR genes show signatures of diversifying selection. Thus, EST analysis of Hpa-infected Arabidopsis is proving to be a powerful method for identifying pathogen effector candidates expressed during infection. Delivery of the Waco9-specific protein RXLR29 in planta revealed that this effector can suppress PAMP-triggered immunity and enhance disease susceptibility. We propose that differences in host colonization can be conditioned by isolate-specific effectors.