Gene loss rather than gene gain is associated with a host jump from monocots to dicots in the Smut Fungus Melanopsichium pennsylvanicum

Genome comparisons reveal accessory genes crucial for the evolution of apple Glomerella leaf spot pathogenicity in Colletotrichum fungi

Apple Glomerella leaf spot (GLS) is an emerging fungal disease caused by Colletotrichum fructicola and other Colletotrichum species. These species are polyphyletic and it is currently unknown how these pathogens convergently evolved to infect apple. We generated chromosome-level genome assemblies of a GLS-adapted isolate and a non-adapted isolate in C. fructicola using long-read sequencing. Additionally, we resequenced 17 C. fructicola and C. aenigma isolates varying in GLS pathogenicity using short-read sequencing. Genome comparisons revealed a conserved bipartite genome architecture involving minichromosomes (accessory chromosomes) shared by C. fructicola and other closely related species within the C. gloeosporioides species complex. Moreover, two repeat-rich genomic regions (1.61 Mb in total) were specifically conserved among GLS-pathogenic isolates in C. fructicola and C. aenigma. Single-gene deletion of 10 accessory genes within the GLS-specific regions of C. fructicola identified three that were essential for GLS pathogenicity. These genes encoded a putative non-ribosomal peptide synthetase, a flavin-binding monooxygenase and a small protein with unknown function. These results highlight the crucial role accessory genes play in the evolution of Colletotrichum pathogenicity and imply the significance of an unidentified secondary metabolite in GLS pathogenesis.

Novel Secreted Effectors Conserved Among Smut Fungi Contribute to the Virulence of Ustilago maydis

Fungal pathogens deploy a set of molecules (proteins, specialized metabolites, and sRNAs), so-called effectors, to aid the infection process. In comparison to other plant pathogens, smut fungi have small genomes and secretomes of 20 Mb and around 500 proteins, respectively. Previous comparative genomic studies have shown that many secreted effector proteins without known domains, i.e., novel, are conserved only in the Ustilaginaceae family. By analyzing the secretomes of 11 species within Ustilaginaceae, we identified 53 core homologous groups commonly present in this lineage. By collecting existing mutants and generating additional ones, we gathered 44 Ustilago maydis strains lacking single core effectors as well as 9 strains containing multiple deletions of core effector gene families. Pathogenicity assays revealed that 20 of these 53 mutant strains were affected in virulence. Among the 33 mutants that had no obvious phenotypic changes, 13 carried additional, sequence-divergent, structurally similar paralogs. We report a virulence contribution of seven previously uncharacterized single core effectors and of one effector family. Our results help to prioritize effectors for understanding U. maydis virulence and provide genetic resources for further characterization. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Manipulation of Auxin Signaling by Smut Fungi during Plant Colonization

A common feature of many plant-colonizing organisms is the exploitation of plant signaling and developmental pathways to successfully establish and proliferate in their hosts. Auxins are central plant growth hormones, and their signaling is heavily interlinked with plant development and immunity responses. Smuts, as one of the largest groups in basidiomycetes, are biotrophic specialists that successfully manipulate their host plants and cause fascinating phenotypes in so far largely enigmatic ways. This review gives an overview of the growing understanding of how and why smut fungi target the central and conserved auxin growth signaling pathways in plants.

The Novel Effector Ue943 Is Essential for Host Plant Colonization by Ustilago esculenta

The smut fungus Ustilago esculenta obligately parasitizes Zizania latifolia and induces smut galls at the stem tips of host plants. Previous research identified a putative secreted protein, Ue943, which is required for the biotrophic phase of U. esculenta but not for the saprophytic phase. Here, we studied the role of Ue943 during the infection process. Conserved homologs of Ue943 were found in smut fungi. Ue943 can be secreted by U. esculenta and localized to the biotrophic interface between fungi and plants. It is required at the early stage of colonization. The Ue943 deletion mutant caused reactive oxygen species (ROS) production and callose deposition in the host plant at 1 and 5 days post inoculation, which led to failed colonization. The virulence deficiency was restored by overexpressing gene Ue943 or Ue943:GFP. Transcriptome analysis further showed a series of changes in plant hormones following ROS production when the host plant was exposed to ΔUe943. We hypothesize that Ue943 might be responsible for ROS suppression or avoidance of recognition by the plant immune system. The mechanism underlying Ue943 requires further study to provide more insights into the virulence of smut fungi.

Pan-Genomics Reveals a New Variation Pattern of Secreted Proteins in Pyricularia oryzae

(1) Background: Pyricularia oryzae, the causal agent of rice blast disease, is one of the major rice pathogens. The complex population structure of P. oryzae facilitates the rapid virulence variations, which make the blast disease a serious challenge for global food security. There is a large body of existing genomics research on P. oryzae, however the population structure at the pan-genome level is not clear, and the mechanism of genetic divergence and virulence variations of different sub-populations is also unknown. (2) Methods: Based on the genome data published in the NCBI, we constructed a pan-genome database of P. oryzae, which consisted of 156 strains (117 isolated from rice and 39 isolated from other hosts). (3) Results: The pan-genome contained a total of 24,100 genes (12,005 novel genes absent in the reference genome 70-15), including 16,911 (~70%) core genes (population frequency ≥95%) and 1378 (~5%) strain-specific genes (population frequency ≤5%). Gene presence-absence variation (PAV) based clustering analysis of the population structure of P. oryzae revealed four subgroups (three from rice and one from other hosts). Interestingly, the cloned avirulence genes and conventional secreted proteins (SPs, with signal peptides) were enriched in the high-frequency regions and significantly associated with transposable elements (TEs), while the unconventional SPs (without signal peptides) were enriched in the low-frequency regions and not associated significantly with TEs. This pan-genome will expand the breadth and depth of the rice blast fungus reference genome, and also serve as a new blueprint for scientists to further study the pathogenic mechanism and virulence variation of the rice blast fungus.

High Nucleotide Substitution Rates Associated with Retrotransposon Proliferation Drive Dynamic Secretome Evolution in Smut Pathogens

Transposable elements (TEs) play a pivotal role in shaping diversity in eukaryotic genomes. The covered smut pathogen on barley, Ustilago hordei, encountered a recent genome expansion. Using long reads, we assembled genomes of 6 U. hordei strains and 3 sister species, to study this genome expansion. We found that larger genome sizes can mainly be attributed to a higher genome fraction of long terminal repeat retrotransposons (LTR-RTs). In the studied smut genomes, LTR-RTs fractions are the largest in U. hordei and are positively correlated with the mating-type locus sizes, which is up to ~560 kb in U. hordei. Furthermore, LTR-RTs were found to be associated with higher nucleotide substitution levels, as these occur in specific genome regions of smut species with a recent LTR-RT proliferation. Moreover, genes in genome regions with higher nucleotide substitution levels generally reside closer to LTR-RTs than other genome regions. Genome regions with many nucleotide substitutions encountered an especially high fraction of CG substitutions, which is not observed for LTR-RT sequences. The high nucleotide substitution levels particularly accelerate the evolution of secretome genes, as their more accessory nature results in substitutions that often lead to amino acid alterations. IMPORTANCE Genomic alteration can be generated through various means, in which transposable elements (TEs) can play a pivotal role. Their mobility causes mutagenesis in itself and can disrupt the function of the sequences they insert into. They also impact genome evolution as their repetitive nature facilitates nonhomologous recombination. Furthermore, TEs have been linked to specific epigenetic genome organizations. We report a recent TE proliferation in the genome of the barley covered smut fungus, Ustilago hordei. This proliferation is associated with a distinct nucleotide substitution regime that has a higher rate and a higher fraction of CG substitutions. This different regime shapes the evolution of genes in subjected genome regions. We hypothesize that TEs may influence the error-rate of DNA polymerase in a hitherto unknown fashion.

Evolution of pathogenicity in obligate fungal pathogens and allied genera

Obligate fungal pathogens (ascomycetes and basidiomycetes) and oomycetes are known to cause diseases in cereal crop plants. They feed on living cells and most of them have learned to bypass the host immune machinery. This paper discusses some of the factors that are associated with pathogenicity drawing examples from ascomycetes, basidiomycetes and oomycetes, with respect to their manifestation in crop plants. The comparisons have revealed a striking similarity in the three groups suggesting convergent pathways that have arisen from three lineages independently leading to an obligate lifestyle. This review has been written with the intent, that new information on adaptation strategies of biotrophs, modifications in pathogenicity strategies and population dynamics will improve current strategies for breeding with stable resistance.

A Bacterial Type Three Secretion-Based Delivery System for Functional Characterization of Sporisorium scitamineum Plant Immune Suppressing Effector Proteins

The facultative biotrophic basidiomycete Sporisorium scitamineum causes smut disease in sugarcane. This study applied an assay to identify S. scitamineum candidate effectors (CEs) with plant immunity suppression activities by delivering them into Nicotiana benthamiana cells via the type-three secretion system of Pseudomonas fluorescens EtHAn. Six CEs were individually cloned into the pEDV6 vector and expressed by P. fluorescens EtHAn for translocation into the plant cells. Three CEs (g1052, g3890, and g5159) could suppress pattern-triggered immunity (PTI) responses with high reproducibility in different coinfiltration experiments with P. syringae pv. tomato DC3000. In addition, three CEs (g1052, g4549, and g5159) were also found to be AvrB-induced suppressors of effector-triggered immunity (ETI), demonstrating for the first time that S. scitamineum can defeat both PTI and ETI responses. A transcriptomic analysis at different stages of infection by the smut fungus of three sugarcane cultivars with contrasting responses to the pathogen revealed that suppressors g1052, g3890, g4549, and g5159 were induced at the early stage of infection. By contrast, the two CEs (g2666 and g6610) that did not exhibit suppression activities expressed only at the late stage of infection. Moreover, genomic structures of the CEs and searches for orthologs in other smut species suggested duplication events and further divergence in CEs evolution of S. scitamineum. Thus, the transient assay applied here demonstrated the potential of pEDV6 and P. fluorescens EtHAn as biological tools for identifying plant immune suppressors from S. scitamineum.

The Sporisorium reilianum Effector Vag2 Promotes Head Smut Disease via Suppression of Plant Defense Responses

Genome comparison between the maize pathogens Ustilago maydis and Sporisorium reilianum revealed a large diversity region (19-1) containing nearly 30 effector gene candidates, whose deletion severely hampers virulence of both fungi. Dissection of the S. reilianum gene cluster resulted in the identification of one major contributor to virulence, virulence-associated gene 2 (vag2; sr10050). Quantitative reverse-transcriptase polymerase chain reaction (qRT-PCR) experiments revealed high expression of vag2 during biotrophic growth of S. reilianum. Using the yeast secretion trap assay, we confirmed the existence of a functional signal peptide allowing protein secretion via the conventional secretory pathway. We identified the cytoplasmic maize chorismate mutase ZmCM2 by yeast two-hybrid screening as a possible interaction partner of Vag2. Interaction of the two proteins in planta was confirmed by bimolecular fluorescence complementation. qRT-PCR experiments revealed vag2-dependent downregulation of salicylic acid (SA)-induced genes, which correlated with higher SA levels in plant tissues colonized by Δvag2 deletion strains relative to S. reilianum wildtype strains. Metabolite analysis suggested rewiring of pathogen-induced SA biosynthesis by preferential conversion of the SA precursor chorismate into the aromatic amino acid precursor prephenate by ZmCM2 in the presence of Vag2. Possibly, the binding of Vag2 to ZmCM2 inhibits the back reaction of the ZmCM2-catalyzed interconversion of chorismate and prephenate, thus contributing to fungal virulence by lowering the plant SA-induced defenses.

Comparative Genomics and Gene Pool Analysis Reveal the Decrease of Genome Diversity and Gene Number in Rice Blast Fungi by Stable Adaption with Rice

Magnaporthe oryzae caused huge losses in rice and wheat production worldwide. Comparing to long-term co-evolution history with rice, wheat-infecting isolates were new-emerging. To reveal the genetic differences between rice and wheat blast on global genomic scale, 109 whole-genome sequences of M. oryzae from rice, wheat, and other hosts were reanalyzed in this study. We found that the rice lineage had gone through stronger selective sweep and fewer conserved genes than those of Triticum and Lolium lineages, which indicated that rice blast fungi adapted to rice by gene loss and rapid evolution of specific loci. Furthermore, 228 genes associated with host adaptation of M. oryzae were found by presence/absence variation (PAV) analyses. The functional annotation of these genes found that the fine turning of genes gain/loss involved with transport and transcription factor, thiol metabolism, and nucleotide metabolism respectively are major mechanisms for rice adaption. This result implies that genetic base of specific host plant may lead to gene gain/loss variation of pathogens, so as to enhance their adaptability to host. Further characterization of these specific loci and their roles in adaption and evaluation of the fungi may eventually lead to understanding of interaction mechanism and develop new strategies of the disease management.

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.

Cross-species analysis between the maize smut fungi Ustilago maydis and Sporisorium reilianum highlights the role of transcriptional change of effector orthologs for virulence and disease

The constitution and regulation of effector repertoires shape host-microbe interactions. Ustilago maydis and Sporisorium reilianum are two closely related smut fungi, which both infect maize but cause distinct disease symptoms. Understanding how effector orthologs are regulated in these two pathogens can therefore provide insights into the evolution of different infection strategies. We tracked the infection progress of U. maydis and S. reilianum in maize leaves and used two distinct infection stages for cross-species RNA-sequencing analyses. We identified 207 of 335 one-to-one effector orthologs as differentially regulated during host colonization, which might reflect the distinct disease development strategies. Using CRISPR-Cas9-mediated gene conversion, we identified two differentially expressed effector orthologs with conserved function between two pathogens. Thus, differential expression of functionally conserved genes might contribute to species-specific adaptation and symptom development. Interestingly, another differentially expressed orthogroup (UMAG_05318/Sr10075) showed divergent protein function, providing a possible case for neofunctionalization. Collectively, we demonstrated that the diversification of effector genes in related pathogens can be caused both by alteration on the transcriptional level and through functional diversification of the encoded effector proteins.

A 20-kb lineage-specific genomic region tames virulence in pathogenic amphidiploid Verticillium longisporum

Amphidiploid fungal Verticillium longisporum strains Vl43 and Vl32 colonize the plant host Brassica napus but differ in their ability to cause disease symptoms. These strains represent two V. longisporum lineages derived from different hybridization events of haploid parental Verticillium strains. Vl32 and Vl43 carry same-sex mating-type genes derived from both parental lineages. Vl32 and Vl43 similarly colonize and penetrate plant roots, but asymptomatic Vl32 proliferation in planta is lower than virulent Vl43. The highly conserved Vl43 and Vl32 genomes include less than 1% unique genes, and the karyotypes of 15 or 16 chromosomes display changed genetic synteny due to substantial genomic reshuffling. A 20 kb Vl43 lineage-specific (LS) region apparently originating from the Verticillium dahliae-related ancestor is specific for symptomatic Vl43 and encodes seven genes, including two putative transcription factors. Either partial or complete deletion of this LS region in Vl43 did not reduce virulence but led to induction of even more severe disease symptoms in rapeseed. This suggests that the LS insertion in the genome of symptomatic V. longisporum Vl43 mediates virulence-reducing functions, limits damage on the host plant, and therefore tames Vl43 from being even more virulent.

A Novel Core Effector Vp1 Promotes Fungal Colonization and Virulence of Ustilago maydis

The biotrophic fungus Ustilago maydis secretes a plethora of uncharacterized effector proteins and causes smut disease in maize. Among the effector genes that are up-regulated during the biotrophic growth in maize, we identified vp1 (virulence promoting 1), which has an expression that was up-regulated and maintained at a high level throughout the life cycle of the fungus. We characterized Vp1 by applying in silico analysis, reverse genetics, phenotypic assessment, microscopy, and protein localization and provided a fundamental understanding of the Vp1 protein in U. maydis. The reduction in fungal virulence and colonization in the vp1 mutant suggests the virulence-promoting function of Vp1. The deletion studies on the NLS (nuclear localization signal) sequence and the protein localization study revealed that the C-terminus of Vp1 is processed after secretion in plant apoplast and could localize to the plant nucleus. The Ustilago hordei ortholog UhVp1 lacks NLS localized in the plant cytoplasm, suggesting that the orthologs might have a distinct subcellular localization. Further complementation studies of the Vp1 orthologs in related smut fungi revealed that none of them could complement the virulence function of U. maydis Vp1, suggesting that UmVp1 could acquire a specialized function via sequence divergence.

The Pleiades are a cluster of fungal effectors that inhibit host defenses

Biotrophic plant pathogens secrete effector proteins to manipulate the host physiology. Effectors suppress defenses and induce an environment favorable to disease development. Sequence-based prediction of effector function is impeded by their rapid evolution rate. In the maize pathogen Ustilago maydis, effector-coding genes frequently organize in clusters. Here we describe the functional characterization of the pleiades, a cluster of ten effector genes, by analyzing the micro- and macroscopic phenotype of the cluster deletion and expressing these proteins in planta. Deletion of the pleiades leads to strongly impaired virulence and accumulation of reactive oxygen species (ROS) in infected tissue. Eight of the Pleiades suppress the production of ROS upon perception of pathogen associated molecular patterns (PAMPs). Although functionally redundant, the Pleiades target different host components. The paralogs Taygeta1 and Merope1 suppress ROS production in either the cytoplasm or nucleus, respectively. Merope1 targets and promotes the auto-ubiquitination activity of RFI2, a conserved family of E3 ligases that regulates the production of PAMP-triggered ROS burst in plants.

Obtaining deeper insights into microbiome diversity using a simple method to block host and nontargets in amplicon sequencing

Profiling diverse microbiomes is revolutionizing our understanding of biological mechanisms and ecologically relevant problems, including metaorganism (host + microbiome) assembly, functions and adaptation. Amplicon sequencing of multiple conserved, phylogenetically informative loci has therefore become an instrumental tool for many researchers. Investigations in many systems are hindered, however, since essential sequencing depth can be lost by amplification of nontarget DNA from hosts or overabundant microorganisms. Here, we introduce "blocking oligos", a low-cost and flexible method using standard oligonucleotides to block amplification of diverse nontargets and software to aid their design. We apply them primarily in leaves, where exceptional challenges with host amplification prevail. A. thaliana-specific blocking oligos applied in eight different target loci reduce undesirable host amplification by up to 90%. To expand applicability, we designed universal 16S and 18S rRNA gene plant blocking oligos for targets that are conserved in diverse plant species and demonstrate that they efficiently block five plant species from five orders spanning monocots and dicots (Bromus erectus, Plantago lanceolata, Lotus corniculatus, Amaranth sp., Arabidopsis thaliana). These can increase alpha diversity discovery without biasing beta diversity patterns and do not compromise microbial load information inherent to plant-derived 16S rRNA gene amplicon sequencing data. Finally, we designed and tested blocking oligos to avoid amplification of 18S rRNA genes of a sporulating oomycete pathogen, demonstrating their effectiveness in applications well beyond plants. Using these tools, we generated a survey of the A. thaliana leaf microbiome based on eight loci targeting bacterial, fungal, oomycete and other eukaryotic microorganisms and discuss complementarity of commonly used amplicon sequencing regions for describing leaf microbiota. This approach has potential to make questions in a variety of study systems more tractable by making amplicon sequencing more targeted, leading to deeper, systems-based insights into microbial discovery. For fast and easy design for blocking oligos for any nontarget DNA in other study systems, we developed a publicly available R package.

Population Genomics of the Maize Pathogen Ustilago maydis: Demographic History and Role of Virulence Clusters in Adaptation

The tight interaction between pathogens and their hosts results in reciprocal selective forces that impact the genetic diversity of the interacting species. The footprints of this selection differ between pathosystems because of distinct life-history traits, demographic histories, or genome architectures. Here, we studied the genome-wide patterns of genetic diversity of 22 isolates of the causative agent of the corn smut disease, Ustilago maydis, originating from five locations in Mexico, the presumed center of origin of this species. In this species, many genes encoding secreted effector proteins reside in so-called virulence clusters in the genome, an arrangement that is so far not found in other filamentous plant pathogens. Using a combination of population genomic statistical analyses, we assessed the geographical, historical, and genome-wide variation of genetic diversity in this fungal pathogen. We report evidence of two partially admixed subpopulations that are only loosely associated with geographic origin. Using the multiple sequentially Markov coalescent model, we inferred the demographic history of the two pathogen subpopulations over the last 0.5 Myr. We show that both populations experienced a recent strong bottleneck starting around 10,000 years ago, coinciding with the assumed time of maize domestication. Although the genome average genetic diversity is low compared with other fungal pathogens, we estimated that the rate of nonsynonymous adaptive substitutions is three times higher in genes located within virulence clusters compared with nonclustered genes, including nonclustered effector genes. These results highlight the role that these singular genomic regions play in the evolution of this pathogen.

Comparative genomics reveals dynamic genome evolution in host specialist ectomycorrhizal fungi

While there has been significant progress characterizing the 'symbiotic toolkit' of ectomycorrhizal (ECM) fungi, how host specificity may be encoded into ECM fungal genomes remains poorly understood. We conducted a comparative genomic analysis of ECM fungal host specialists and generalists, focusing on the specialist genus Suillus. Global analyses of genome dynamics across 46 species were assessed, along with targeted analyses of three classes of molecules previously identified as important determinants of host specificity: small secreted proteins (SSPs), secondary metabolites (SMs) and G-protein coupled receptors (GPCRs). Relative to other ECM fungi, including other host specialists, Suillus had highly dynamic genomes including numerous rapidly evolving gene families and many domain expansions and contractions. Targeted analyses supported a role for SMs but not SSPs or GPCRs in Suillus host specificity. Phylogenomic-based ancestral state reconstruction identified Larix as the ancestral host of Suillus, with multiple independent switches between white and red pine hosts. These results suggest that like other defining characteristics of the ECM lifestyle, host specificity is a dynamic process at the genome level. In the case of Suillus, both SMs and pathways involved in the deactivation of reactive oxygen species appear to be strongly associated with enhanced host specificity.

A cell surface-exposed protein complex with an essential virulence function in Ustilago maydis

Plant pathogenic fungi colonizing living plant tissue secrete a cocktail of effector proteins to suppress plant immunity and reprogramme host cells. Although many of these effectors function inside host cells, delivery systems used by pathogenic bacteria to translocate effectors into host cells have not been detected in fungi. Here, we show that five unrelated effectors and two membrane proteins from Ustilago maydis, a biotrophic fungus causing smut disease in corn, form a stable protein complex. All seven genes appear co-regulated and are only expressed during colonization. Single mutants arrest in the epidermal layer, fail to suppress host defence responses and fail to induce non-host resistance, two reactions that likely depend on translocated effectors. The complex is anchored in the fungal membrane, protrudes into host cells and likely contacts channel-forming plant plasma membrane proteins. Constitutive expression of all seven complex members resulted in a surface-exposed form in cultured U. maydis cells. As orthologues of the complex-forming proteins are conserved in smut fungi, the complex may become an interesting fungicide target.

What are the Top 10 Unanswered Questions in Molecular Plant-Microbe Interactions?

This article is part of the Top 10 Unanswered Questions in MPMI invited review series.The past few decades have seen major discoveries in the field of molecular plant-microbe interactions. As the result of technological and intellectual advances, we are now able to answer questions at a level of mechanistic detail that we could not have imagined possible 20 years ago. The MPMI Editorial Board felt it was time to take stock and reassess. What big questions remain unanswered? We knew that to identify the fundamental, overarching questions that drive our research, we needed to do this as a community. To reach a diverse audience of people with different backgrounds and perspectives, working in different areas of plant-microbe interactions, we queried the more than 1,400 participants at the 2019 International Congress on Molecular Plant-Microbe Interactions meeting in Glasgow. This group effort resulted in a list of ten, broad-reaching, fundamental questions that influence and inform our research. Here, we introduce these Top 10 unanswered questions, giving context and a brief description of the issues. Each of these questions will be the subject of a detailed review in the coming months. We hope that this process of reflecting on what is known and unknown and identifying the themes that underlie our research will provide a framework to use going forward, giving newcomers a sense of the mystery of the big questions and inspiring new avenues and novel insights.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

The Effector Repertoire of the Hop Downy Mildew Pathogen Pseudoperonospora humuli

Pseudoperonospora humuli is an obligate biotrophic oomycete that causes downy mildew (DM), one of the most destructive diseases of cultivated hop that can lead to 100% crop loss in susceptible cultivars. We used the published genome of P. humuli to predict the secretome and effectorome and analyze the transcriptome variation among diverse isolates and during infection of hop leaves. Mining the predicted coding genes of the sequenced isolate OR502AA of P. humuli revealed a secretome of 1,250 genes. We identified 296 RXLR and RXLR-like effector-encoding genes in the secretome. Among the predicted RXLRs, there were several WY-motif-containing effectors that lacked canonical RXLR domains. Transcriptome analysis of sporangia from 12 different isolates collected from various hop cultivars revealed 754 secreted proteins and 201 RXLR effectors that showed transcript evidence across all isolates with reads per kilobase million (RPKM) values > 0. RNA-seq analysis of OR502AA-infected hop leaf samples at different time points after infection revealed highly expressed effectors that may play a relevant role in pathogenicity. Quantitative RT-PCR analysis confirmed the differential expression of selected effectors. We identified a set of P. humuli core effectors that showed transcript evidence in all tested isolates and elevated expression during infection. These effectors are ideal candidates for functional analysis and effector-assisted breeding to develop DM resistant hop cultivars.

Antivirulence and avirulence genes in human pathogenic fungi

Opportunistic commensal and environmental fungi can cause superficial to systemic diseases in humans. But how did these pathogens adapt to infect us and how does host-pathogen co-evolution shape their virulence potential? During evolution toward pathogenicity, not only do microorganisms gain virulence genes, but they also tend to lose non-adaptive genes in the host niche. Additionally, virulence factors can become detrimental during infection when they trigger host recognition. The loss of non-adaptive genes as well as the loss of the virulence potential of genes by adaptations to the host has been investigated in pathogenic bacteria and phytopathogenic fungi, where they are known as antivirulence and avirulence genes, respectively. However, these concepts are nearly unknown in the field of pathogenic fungi of humans. We think that this unnecessarily limits our view of human-fungal interplay, and that much could be learned if we applied a similar framework to aspects of these interactions. In this review, we, therefore, define and adapt the concepts of antivirulence and avirulence genes for human pathogenic fungi. We provide examples for analogies to antivirulence genes of bacterial pathogens and to avirulence genes of phytopathogenic fungi. Introducing these terms to the field of pathogenic fungi of humans can help to better comprehend the emergence and evolution of fungal virulence and disease.

The functionally conserved effector Sta1 is a fungal cell wall protein required for virulence in Ustilago maydis

The biotrophic fungus Ustilago maydis causes the smut disease of maize. The interaction with its host and induction of characteristic tumors are governed largely by secreted effectors whose function is mostly unknown. To identify effectors with a prominent role in virulence, we used RNA sequencing and found that the gene sta1 is upregulated during early stages of infection. We characterized Sta1 by comparative genomics, reverse genetics, protein localization, stress assays, and microscopy. sta1 mutants show a dramatic reduction of virulence and show altered colonization of tissue neighboring the vascular bundles. Functional orthologues of Sta1 are found in related smut pathogens infecting monocot and dicot plants. Sta1 is secreted by budding cells but is attached to the cell wall of filamentous hyphae. Upon constitutive expression of Sta1, fungal filaments become susceptible to Congo red, β-glucanase, and chitinase, suggesting that Sta1 alters the structure of the fungal cell wall. Constitutive or delayed expression of sta1 during plant colonization negatively impacts on virulence. Our results suggest that Sta1 is a novel kind of effector, which needs to modify the hyphal cell wall to allow hyphae to be accommodated in tissue next to the vascular bundles.

Distinctive Gene Expression Profiles and Effectors Consistent With Host Specificity in Two Formae Speciales of Marssonina brunnea

The knowledge on the host specificity of a pathogen underlying an interaction is becoming an urgent necessity for global warming. In this study, the gene expression profiles and the roles of effectors in host specificity were integrally characterized with two formae speciales, multigermtubi and monogermtubi, of a hemibiotrophic pathogen Marssonina brunnea when they were infecting respective susceptible poplar hosts. With a functional genome comparison referring to a de novo transcriptome of M. brunnea and Pathogen-Host Interaction database functional annotations, the multigermtubi strain showed abundant and significant differentially expressed unigenes (DEGs) (more than 40%) in colonizing the initial invasion stage and in the necrotrophic stage. The monogermtubi strain induced less than 10% of DEGs in the initial invasion stage but which abruptly increased to more than 80% DEGs in the necrotrophic stage. Both strains induced the least DEGs in the biotrophic stage compared to the initial invasion and necrotrophic stages. The orthologs of the effector genes Ecp6, PemG1, XEG1, ACE1, and Mg3LysM were exclusively induced by one of the two formae speciales depending on the infection stages. Some unigenes homologous to carbohydrate lytic enzyme genes, especially pectate lyases, were notably induced with multigermtubi forma specialis infection but not expressed in the monogermtubi forma specialis at an earlier infection stage. The extraordinary differences in the functional genome level between the two formae speciales of M. brunnea could be fundamental to exploring their host specificity determinant and evolution. This study also firstly provided the fungal transcriptome of the monogermtubi forma specialis for M. brunnea.

The modular nature of protein evolution: domain rearrangement rates across eukaryotic life

BACKGROUND

Modularity is important for evolutionary innovation. The recombination of existing units to form larger complexes with new functionalities spares the need to create novel elements from scratch. In proteins, this principle can be observed at the level of protein domains, functional subunits which are regularly rearranged to acquire new functions.

RESULTS

In this study we analyse the mechanisms leading to new domain arrangements in five major eukaryotic clades (vertebrates, insects, fungi, monocots and eudicots) at unprecedented depth and breadth. This allows, for the first time, to directly compare rates of rearrangements between different clades and identify both lineage specific and general patterns of evolution in the context of domain rearrangements. We analyse arrangement changes along phylogenetic trees by reconstructing ancestral domain content in combination with feasible single step events, such as fusion or fission. Using this approach we explain up to 70% of all rearrangements by tracing them back to their precursors. We find that rates in general and the ratio between these rates for a given clade in particular, are highly consistent across all clades. In agreement with previous studies, fusions are the most frequent event leading to new domain arrangements. A lineage specific pattern in fungi reveals exceptionally high loss rates compared to other clades, supporting recent studies highlighting the importance of loss for evolutionary innovation. Furthermore, our methodology allows us to link domain emergences at specific nodes in the phylogenetic tree to important functional developments, such as the origin of hair in mammals.

CONCLUSIONS

Our results demonstrate that domain rearrangements are based on a canonical set of mutational events with rates which lie within a relatively narrow and consistent range. In addition, gained knowledge about these rates provides a basis for advanced domain-based methodologies for phylogenetics and homology analysis which complement current sequence-based methods.

Smut fungal strategies for the successful infection

The smut fungi include a large number of plant pathogens that establish obligate biotrophic relationships with their host. Throughout the whole life inside plant tissue, smut fungi keep plant cells alive and acquire nutrients via biotrophic interfaces. This mini-review mainly summarizes the interactions between smut fungi and their host plants during the infection process. Despite various strategies recruited by plants to defense invading pathogens, smut fungi successfully evolved an arsenal for colonization. Mating of two compatible haploids gives rise to parasitic mycelium, which can sense plant surface cues such as fatty acids and hydrophobic surface, and induce the formation of appressoria for surface penetration. Plants can recognize fungal invading and activate defense response, including callose and lignin deposition, programmed cell death, and SA signaling pathway. To suppress plant immunity and alter the metabolic pathway of host plants, a cocktail of effectors is secreted by smut fungi depending on the plant organ and cell type that is infected.

New insights into the evolution of host specificity of three Penicillium species and the pathogenicity of P. Italicum involving the infection of Valencia orange (Citrus sinensis)

Blue and green molds, the common phenotypes of post-harvest diseases in fruits, are mainly caused by Penicillium fungal species, including P. italicum, P. digitatum, and P. expansum. We sequenced and assembled the genome of a P. italicum strain, which contains 31,034,623 bp with 361 scaffolds and 627 contigs. The mechanisms underlying the evolution of host specificity among the analyzed Penicillium species were associated with the expansion of protein families, genome restructuring, horizontal gene transfer, and positive selection pressure. A dual-transcriptome analysis following the infection of Valencia orange (Citrus sinensis) by P. italicum resulted in the annotation of 9,307 P. italicum genes and 24,591 Valencia orange genes. The pathogenicity of P. italicum may be due to the activation of effectors, including 51 small secreted cysteine-rich proteins, 110 carbohydrate-active enzymes, and 12 G protein-coupled receptors. Additionally, 211 metabolites related to the interactions between P. italicum and Valencia orange were identified by gas chromatography-time of flight mass spectrography, three of which were further confirmed by ultra-high performance liquid chromatography triple quadrupole mass spectrometry. A metabolomics analysis indicated that P. italicum pathogenicity is associated with the sphingolipid and salicylic acid signaling pathways. Moreover, a correlation analysis between the metabolite contents and gene expression levels suggested that P. italicum induces carbohydrate metabolism in Valencia orange fruits as part of its infection strategy. This study provides useful information regarding the genomic determinants that drive the evolution of host specificity in Penicillium species and clarifies the host-plant specificity during the infection of Valencia orange by P. italicum.

IMPORTANCE

P. italicum GL_Gan1, a local strain in Guangzhou, China, was sequenced. Comparison of the genome of P. italicum GL_Gan1 with other pathogenic Penicillium species, P. digitatum and P. expansum, revealed that the expansion of protein families, genome restructuring, HGT, and positive selection pressure were related to the host range expansion of the analyzed Penicillium species. Moreover, gene gains or losses might be associated with the speciation of these Penicillium species. In addition, the molecular basis of host-plant specificity during the infection of Valencia orange (Citrus sinensis) by P. italicum was also elucidated by transcriptomic and metabolomics analysis. The data presented herein may be useful for further elucidating the molecular basis of the evolution of host specificity of Penicillium species and for illustrating the host-plant specificity during the infection of Valencia orange by P. italicum.

Two Is Better Than One: Studying Ustilago bromivora-Brachypodium Compatibility by Using a Hybrid Pathogen

Pathogenic fungi can have devastating effects on agriculture and health. One potential challenge in dealing with pathogens is the possibility of a host jump (i.e., when a pathogen infects a new host species). This can lead to the emergence of new diseases or complicate the management of existing threats. We studied host specificity by using a hybrid fungus formed by mating two closely related fungi: Ustilago bromivora, which normally infects Brachypodium spp., and U. hordei, which normally infects barley. Although U. hordei was unable to infect Brachypodium spp., the hybrid could. These hybrids also displayed the same mating-type bias that had been observed in U. bromivora and provide evidence of a dominant spore-killer-like system on the sex chromosome of U. bromivora. By analyzing the genomic composition of 109 hybrid strains, backcrossed with U. hordei over four generations, we identified three regions associated with infection on Brachypodium spp. and 75 potential virulence candidates. The most strongly associated region was located on chromosome 8, where seven genes encoding predicted secreted proteins were identified. The fact that we identified several regions relevant for pathogenicity on Brachypodium spp. but that none were essential suggests that host specificity, in the case of U. bromivora, is a multifactorial trait which can be achieved through different subsets of virulence factors.

An evolutionary framework for host shifts - jumping ships for survival

Host jumping is a process by which pathogens settle in new host groups. It is a cornerstone in the evolution of pathogens, as it leads to pathogen diversification. It is unsurprising that host jumping is observed in facultative pathogens, as they can reproduce even if they kill their hosts. However, host jumps were thought to be rare in obligate biotrophic pathogens, but molecular phylogenetics has revealed that the opposite is true. Here, I review some concepts and recent findings and present several hypotheses on the matter. In short, pathogens evolve and diversify via host jumps, followed by radiation, specialisation and speciation. Host jumps are facilitated by, for example, effector innovations, stress, compatible pathogens and physiological similarities. Host jumping, subsequent establishment, and speciation takes place rapidly - within centuries and millennia rather than over millions of years. If pathogens are unable to evolve into neutral or mutualistic interactions with their hosts, they will eventually be removed from the host population, despite balancing trade-offs. Thus, generally, plant pathogens only survive in the course of evolution if they jump hosts. This is also reflected by the diversity patterns observed in many genera of plant pathogens, where it leads to a mosaic pattern of host groups over time, in which the original host group becomes increasingly obscure.

Understanding Adaptation, Coevolution, Host Specialization, and Mating System in Castrating Anther-Smut Fungi by Combining Population and Comparative Genomics

Anther-smut fungi provide a powerful system to study host-pathogen specialization and coevolution, with hundreds of Microbotryum species specialized on diverse Caryophyllaceae plants, castrating their hosts through manipulation of the hosts' reproductive organs to facilitate disease transmission. Microbotryum fungi have exceptional genomic characteristics, including dimorphic mating-type chromosomes, that make this genus anexcellent model for studying the evolution of mating systems and their influence on population genetics structure and adaptive potential. Important insights into adaptation, coevolution, host specialization, and mating system evolution have been gained using anther-smut fungi, with new insights made possible by the recent advent of genomic approaches. We illustrate with Microbotryum case studies how using a combination of comparative genomics, population genomics, and transcriptomics approaches enables the integration of different evolutionary perspectives across different timescales. We also highlight current challenges and suggest future studies that will contribute to advancing our understanding of the mechanisms underlying adaptive processes in populations of fungal pathogens.

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.

A Dead Gene Walking: Convergent Degeneration of a Clade of MADS-Box Genes in Crucifers

Genes are "born," and eventually they "die." These processes shape the phenotypic evolution of organisms and are hence of great biological interest. If genes die in plants, they generally do so quite rapidly. Here, we describe the fate of GOA-like genes that evolve in a dramatically different manner. GOA-like genes belong to the subfamily of Bsister genes of MIKC-type MADS-box genes. Typical MIKC-type genes encode conserved transcription factors controlling plant development. We show that ABS-like genes, a clade of Bsister genes, are indeed highly conserved in crucifers (Brassicaceae) maintaining the ancestral function of Bsister genes in ovule and seed development. In contrast, their closest paralogs, the GOA-like genes, have been undergoing convergent gene death in Brassicaceae. Intriguingly, erosion of GOA-like genes occurred after millions of years of coexistence with ABS-like genes. We thus describe Delayed Convergent Asymmetric Degeneration, a so far neglected but possibly frequent pattern of duplicate gene evolution that does not fit classical scenarios. Delayed Convergent Asymmetric Degeneration of GOA-like genes may have been initiated by a reduction in the expression of an ancestral GOA-like gene in the stem group of Brassicaceae and driven by dosage subfunctionalization. Our findings have profound implications for gene annotations in genomics, interpreting patterns of gene evolution and using genes in phylogeny reconstructions of species.

Mating-type loci of Ustilago esculenta are essential for mating and development

Ustilago esculenta is closely related to the smut fungus Ustilago maydis and, in an endophytic-like life in the plant Zizania latifolia, only infects host stems and causes swollen stems to form edible galls called Jiaobai in China. In order to study its different modes of invasion and sites of symptom development from other smut fungi at the molecular level, we first characterized the a and b mating-type loci of U. esculenta. The a loci contained three a mating-type alleles, encoding two pheromones and one pheromone receptor per allele. The pheromone/receptor system controlled the conjugation formation, the initial step of mating, in which each pheromone was specific for recognition by only one mating partner. In addition, there are at least three b alleles identified in U. esculenta, encoding two subunits of heterodimeric homeodomain transcription factors bE and bW, responsible for hyphal growth and invasiveness. Hyphal formation, elongation and invasion after mating of two compatible partners occurred, only when a heterodimer complex was formed by the bE and bW proteins derived from different alleles. We also demonstrated that even with only one paired pheromone-pheromone receptor, the active b locus heterodimer triggered hyphal growth and infection.

A polyphasic method for the identification of aflatoxigenic Aspergilla from cashew nuts

The invasion of food by toxigenic fungi is a threat to public health. This study aimed at enumerating the microbial profile, detection of aflatoxin producing genes and quantification of the levels of aflatoxin contamination of cashew nuts meant for human consumption. A polyphasic method of analysis using newly formulated β-Cyclodextrin Neutral Red Desiccated coconut agar (β-CDNRDCA) and Yeast Extract Sucrose agar (YES) with Thin Layer Chromatography (TLC), Polymerase Chain Reaction (PCR) and High Performance Liquid Chromatographic (HPLC) method was adopted in determining the aflatoxigenic potential of the isolates, the presence of aflatoxin biosynthetic gene (aflM, aflD, aflR, aflJ omt-A) and estimation of the total aflatoxin content of the nuts. The fungal counts ranged from 2.0 to 2.4 log₁₀cfu/g and sixty-three fungal isolates belonging to 18 genera and 34 species were isolated. The Aspergillus spp. were the most frequently isolated (50.79%) while Trichoderma spp. (1.59%) were the least. and fluorescence production was enhanced on the newly formulated β-CDNRDCA by the aflatoxigenic species. The aflD gene was amplified in all the isolates while aflM, aflR and aflJ gene were each amplified in 77.77% of the isolates and omt-A gene in 70.37%. The aflatoxin content of the nuts ranged from 0.03 to 0.77 µg/kg and were below the 4 µg/kg EU recommended limit for total aflatoxins. The present work confirms that a single method of analysis may not be sufficient to screen for the presence of aflatoxins in foods, as with a combination of different methods.

The core effector Cce1 is required for early infection of maize by Ustilago maydis

The biotrophic pathogen Ustilago maydis, the causative agent of corn smut disease, infects one of the most important crops worldwide - Zea mays. To successfully colonize its host, U. maydis secretes proteins, known as effectors, that suppress plant defense responses and facilitate the establishment of biotrophy. In this work, we describe the U. maydis effector protein Cce1. Cce1 is essential for virulence and is upregulated during infection. Through microscopic analysis and in vitro assays, we show that Cce1 is secreted from hyphae during filamentous growth of the fungus. Strikingly, Δcce1 mutants are blocked at early stages of infection and induce callose deposition as a plant defense response. Cce1 is highly conserved among smut fungi and the Ustilago bromivora ortholog complemented the virulence defect of the SG200Δcce1 deletion strain. These data indicate that Cce1 is a core effector with apoplastic localization that is essential for U. maydis to infect its host.

Comparative whole-genome analysis reveals artificial selection effects on Ustilago esculenta genome

Ustilago esculenta, infects Zizania latifolia, and induced host stem swollen to be a popular vegetable called Jiaobai in China. It is the long-standing artificial selection that maximizes the occurrence of favourable Jiaobai, and thus maintaining the plant-fungi interaction and modulating the fungus evolving from plant pathogen to entophyte. In this study, whole genome of U. esculenta was sequenced and transcriptomes of the fungi and its host were analysed. The 20.2 Mb U. esculenta draft genome of 6,654 predicted genes including mating, primary metabolism, secreted proteins, shared a high similarity to related Smut fungi. But U. esculenta prefers RNA silencing not repeat-induced point in defence and has more introns per gene, indicating relatively slow evolution rate. The fungus also lacks some genes in amino acid biosynthesis pathway which were filled by up-regulated host genes and developed distinct amino acid response mechanism to balance the infection-resistance interaction. Besides, U. esculenta lost some surface sensors, important virulence factors and host range-related effectors to maintain the economic endophytic life. The elucidation of the U. esculenta genomic information as well as expression profiles can not only contribute to more comprehensive insights into the molecular mechanism underlying artificial selection but also into smut fungi-host interactions.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

Comparative Genomics of Smut Pathogens: Insights From Orphans and Positively Selected Genes Into Host Specialization

Host specialization is a key evolutionary process for the diversification and emergence of new pathogens. However, the molecular determinants of host range are poorly understood. Smut fungi are biotrophic pathogens that have distinct and narrow host ranges based on largely unknown genetic determinants. Hence, we aimed to expand comparative genomics analyses of smut fungi by including more species infecting different hosts and to define orphans and positively selected genes to gain further insights into the genetics basis of host specialization. We analyzed nine lineages of smut fungi isolated from eight crop and non-crop hosts: maize, barley, sugarcane, wheat, oats, Zizania latifolia (Manchurian rice), Echinochloa colona (a wild grass), and Persicaria sp. (a wild dicot plant). We assembled two new genomes: Ustilago hordei (strain Uhor01) isolated from oats and U. tritici (strain CBS 119.19) isolated from wheat. The smut genomes were of small sizes, ranging from 18.38 to 24.63 Mb. U. hordei species experienced genome expansions due to the proliferation of transposable elements and the amount of these elements varied among the two strains. Phylogenetic analysis confirmed that Ustilago is not a monophyletic genus and, furthermore, detected misclassification of the U. tritici specimen. The comparison between smut pathogens of crop and non-crop hosts did not reveal distinct signatures, suggesting that host domestication did not play a dominant role in shaping the evolution of smuts. We found that host specialization in smut fungi likely has a complex genetic basis: different functional categories were enriched in orphans and lineage-specific selected genes. The diversification and gain/loss of effector genes are probably the most important determinants of host specificity.

Comparative Methods for Molecular Determination of Host-Specificity Factors in Plant-Pathogenic Fungi

Many plant-pathogenic fungi are highly host-specific. In most cases, host-specific interactions evolved at the time of speciation of the respective host plants. However, host jumps have occurred quite frequently, and still today the greatest threat for the emergence of new fungal diseases is the acquisition of infection capability of a new host by an existing plant pathogen. Understanding the mechanisms underlying host-switching events requires knowledge of the factors determining host-specificity. In this review, we highlight molecular methods that use a comparative approach for the identification of host-specificity factors. These cover a wide range of experimental set-ups, such as characterization of the pathosystem, genotyping of host-specific strains, comparative genomics, transcriptomics and proteomics, as well as gene prediction and functional gene validation. The methods are described and evaluated in view of their success in the identification of host-specificity factors and the understanding of their functional mechanisms. In addition, potential methods for the future identification of host-specificity factors are discussed.

The birth and death of effectors in rapidly evolving filamentous pathogen genomes

Plant pathogenic fungi and oomycetes are major risks to food security due to their evolutionary success in overcoming plant defences. Pathogens produce effectors to interfere with host defences and metabolism. These effectors are often encoded in rapidly evolving compartments of the genome. We review how effector genes emerged and were lost in pathogen genomes drawing on the links between effector evolution and chromosomal rearrangements. Some new effectors entered pathogen genomes via horizontal transfer or introgression. However, new effector functions also arose through gene duplication or from previously non-coding sequences. The evolutionary success of an effector is tightly linked to its transcriptional regulation during host colonization. Some effectors converged on an epigenetic control of expression imposed by genomic defences against transposable elements. Transposable elements were also drivers of effector diversification and loss that led to mosaics in effector presence-absence variation. Such effector mosaics within species was the foundation for rapid pathogen adaptation.

Positively Selected Effector Genes and Their Contribution to Virulence in the Smut Fungus Sporisorium reilianum

Plants and fungi display a broad range of interactions in natural and agricultural ecosystems ranging from symbiosis to parasitism. These ecological interactions result in coevolution between genes belonging to different partners. A well-understood example is secreted fungal effector proteins and their host targets, which play an important role in pathogenic interactions. Biotrophic smut fungi (Basidiomycota) are well-suited to investigate the evolution of plant pathogens, because several reference genomes and genetic tools are available for these species. Here, we used the genomes of Sporisorium reilianum f. sp. zeae and S. reilianum f. sp. reilianum, two closely related formae speciales infecting maize and sorghum, respectively, together with the genomes of Ustilago hordei, Ustilago maydis, and Sporisorium scitamineum to identify and characterize genes displaying signatures of positive selection. We identified 154 gene families having undergone positive selection during species divergence in at least one lineage, among which 77% were identified in the two investigated formae speciales of S. reilianum. Remarkably, only 29% of positively selected genes encode predicted secreted proteins. We assessed the contribution to virulence of nine of these candidate effector genes in S. reilianum f. sp. zeae by deleting individual genes, including a homologue of the effector gene pit2 previously characterized in U. maydis. Only the pit2 deletion mutant was found to be strongly reduced in virulence. Additional experiments are required to understand the molecular mechanisms underlying the selection forces acting on the other candidate effector genes, as well as the large fraction of positively selected genes encoding predicted cytoplasmic proteins.

Widespread signatures of selection for secreted peptidases in a fungal plant pathogen

BACKGROUND

Fungal plant pathogens secrete a large arsenal of hydrolytic enzymes during the course of infection, including peptidases. Secreted peptidases have been extensively studied for their role as effectors. In this study, we combined transcriptomics, comparative genomics and evolutionary analyses to investigate all 39 secreted peptidases in the fungal wheat pathogen Zymoseptoria tritici and its close relatives Z. pseudotritici and Z. ardabiliae.

RESULTS

RNA-seq data revealed that a majority of the secreted peptidases displayed differential transcription during the course of Z. tritici infection, indicative of specialization for different stages in the life cycle. Evolutionary analyses detected widespread evidence of adaptive evolution acting on at least 28 of the peptidases. A few peptidases displayed lineage-specific rates of molecular evolution, suggesting altered selection pressure in Z. tritici following host specialization on domesticated wheat. The peptidases belonging to MEROPS families A1 and G1 emerged as a particularly interesting group that may play key roles in host-pathogen co-evolution, host adaptation and pathogenicity. Sister genes in the A1 and G1 families showed accelerated substitution rates after gene duplications.

CONCLUSIONS

These results suggest widespread evolution of secreted peptidases leading to novel gene functions, consistent with predicted models of "escape from adaptive conflict" and "neo-functionalization". Our analyses identified candidate genes worthy of functional analyses that may encode effector functions, for example by suppressing plant defenses during the biotrophic phase of infection.

Pangenome analyses of the wheat pathogen Zymoseptoria tritici reveal the structural basis of a highly plastic eukaryotic genome

Background

Structural variation contributes substantially to polymorphism within species. Chromosomal rearrangements that impact genes can lead to functional variation among individuals and influence the expression of phenotypic traits. Genomes of fungal pathogens show substantial chromosomal polymorphism that can drive virulence evolution on host plants. Assessing the adaptive significance of structural variation is challenging, because most studies rely on inferences based on a single reference genome sequence.

Results

We constructed and analyzed the pangenome of Zymoseptoria tritici, a major pathogen of wheat that evolved host specialization by chromosomal rearrangements and gene deletions. We used single-molecule real-time sequencing and high-density genetic maps to assemble multiple genomes. We annotated the gene space based on transcriptomics data that covered the infection life cycle of each strain. Based on a total of five telomere-to-telomere genomes, we constructed a pangenome for the species and identified a core set of 9149 genes. However, an additional 6600 genes were exclusive to a subset of the isolates. The substantial accessory genome encoded on average fewer expressed genes but a larger fraction of the candidate effector genes that may interact with the host during infection. We expanded our analyses of the pangenome to a worldwide collection of 123 isolates of the same species. We confirmed that accessory genes were indeed more likely to show deletion polymorphisms and loss-of-function mutations compared to core genes.

Conclusions

The pangenome construction of a highly polymorphic eukaryotic pathogen showed that a single reference genome significantly underestimates the gene space of a species. The substantial accessory genome provides a cradle for adaptive evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0457-4) contains supplementary material, which is available to authorized users.

Comparative genomics reveals cotton-specific virulence factors in flexible genomic regions in Verticillium dahliae and evidence of horizontal gene transfer from Fusarium

Verticillium dahliae isolates are most virulent on the host from which they were originally isolated. Mechanisms underlying these dominant host adaptations are currently unknown. We sequenced the genome of V. dahliae Vd991, which is highly virulent on its original host, cotton, and performed comparisons with the reference genomes of JR2 (from tomato) and VdLs.17 (from lettuce). Pathogenicity-related factor prediction, orthology and multigene family classification, transcriptome analyses, phylogenetic analyses, and pathogenicity experiments were performed. The Vd991 genome harbored several exclusive, lineage-specific (LS) genes within LS regions (LSRs). Deletion mutants of the seven genes within one LSR (G-LSR2) in Vd991 were less virulent only on cotton. Integration of G-LSR2 genes individually into JR2 and VdLs.17 resulted in significantly enhanced virulence on cotton but did not affect virulence on tomato or lettuce. Transcription levels of the seven LS genes in Vd991 were higher during the early stages of cotton infection, as compared with other hosts. Phylogenetic analyses suggested that G-LSR2 was acquired from Fusarium oxysporum f. sp. vasinfectum through horizontal gene transfer. Our results provide evidence that horizontal gene transfer from Fusarium to Vd991 contributed significantly to its adaptation to cotton and may represent a significant mechanism in the evolution of an asexual plant pathogen.

Lessons in Effector and NLR Biology of Plant-Microbe Systems

A diversity of plant-associated organisms secrete effectors-proteins and metabolites that modulate plant physiology to favor host infection and colonization. However, effectors can also activate plant immune receptors, notably nucleotide-binding domain and leucine-rich repeat region (NLR)-containing proteins, enabling plants to fight off invading organisms. This interplay between effectors, their host targets, and the matching immune receptors is shaped by intricate molecular mechanisms and exceptionally dynamic coevolution. In this article, we focus on three effectors, AVR-Pik, AVR-Pia, and AVR-Pii, from the rice blast fungus Magnaporthe oryzae (syn. Pyricularia oryzae), and their corresponding rice NLR immune receptors, Pik, Pia, and Pii, to highlight general concepts of plant-microbe interactions. We draw 12 lessons in effector and NLR biology that have emerged from studying these three little effectors and are broadly applicable to other plant-microbe systems.

Lessons in Effector and NLR Biology of Plant-Microbe Systems

A diversity of plant-associated organisms secrete effectors-proteins and metabolites that modulate plant physiology to favor host infection and colonization. However, effectors can also activate plant immune receptors, notably nucleotide-binding domain and leucine-rich repeat region (NLR)-containing proteins, enabling plants to fight off invading organisms. This interplay between effectors, their host targets, and the matching immune receptors is shaped by intricate molecular mechanisms and exceptionally dynamic coevolution. In this article, we focus on three effectors, AVR-Pik, AVR-Pia, and AVR-Pii, from the rice blast fungus Magnaporthe oryzae (syn. Pyricularia oryzae), and their corresponding rice NLR immune receptors, Pik, Pia, and Pii, to highlight general concepts of plant-microbe interactions. We draw 12 lessons in effector and NLR biology that have emerged from studying these three little effectors and are broadly applicable to other plant-microbe systems.