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

The new CFEM protein CgCsa required for Fe3+ homeostasis regulates the growth, development, and pathogenicity of Colletotrichum gloeosporioides

Secreted common fungal extracellular membrane (CFEM) domain proteins have been implicated in multiple biological functions in fungi. However, it is still largely unknown whether the ferric iron (Fe3+), as an important trace element, was involved with the biological function of CFEM proteins. In this study, a new CFEM protein CgCsa, with high expression levels at the early inoculation stage on peppers by Colletotrichum gloeosporioides was investigated. Deletion of the targeted gene CgCsa revealed multiple biological roles in hyphal growth restriction, highly reduced conidial yield, delayed conidial germination, abnormal appressorium with elongated bud tubes, and significantly reduced virulence of C. gloeosporioides. Moreover, in CgCsa mutants, the expression levels of four cell wall synthesis-related genes were downregulated, and cell membrane permeability and electrical conductivity were increased. Compared to the wild-type, the CgCsa mutants downregulated expressions of iron transport-related genes, in addition, its three-dimensional structure was capable binding with iron. Increase in the Fe3+ concentration in the culture medium partially recovered the functions of ΔCgCsa mutant. This is probably the first report to show the association between CgCsa and iron homeostasis in C. gloeosporioides. The results suggest an alternative pathway for controlling plant fungal diseases by deplete their trace elements.

Quantitative pathogenicity and host adaptation in a fungal plant pathogen revealed by whole-genome sequencing

Knowledge of genetic determinism and evolutionary dynamics mediating host-pathogen interactions is essential to manage fungal plant diseases. Studies on the genetic architecture of fungal pathogenicity often focus on large-effect effector genes triggering strong, qualitative resistance. It is not clear how this translates to predominately quantitative interactions. Here, we use the Zymoseptoria tritici-wheat model to elucidate the genetic architecture of quantitative pathogenicity and mechanisms mediating host adaptation. With a multi-host genome-wide association study, we identify 19 high-confidence candidate genes associated with quantitative pathogenicity. Analysis of genetic diversity reveals that sequence polymorphism is the main evolutionary process mediating differences in quantitative pathogenicity, a process that is likely facilitated by genetic recombination and transposable element dynamics. Finally, we use functional approaches to confirm the role of an effector-like gene and a methyltransferase in phenotypic variation. This study highlights the complex genetic architecture of quantitative pathogenicity, extensive diversifying selection and plausible mechanisms facilitating pathogen adaptation.

Epiphytic proliferation of Zymoseptoria tritici isolates on resistant wheat leaves

The wheat pathogen Zymoseptoria tritici is capable of a long period of pre-invasive epiphytic growth. Studies have shown that virulent isolates vary in the extent, duration and growth form of this epiphytic growth, and the fungus has been observed to undergo behaviours such as asexual reproduction by budding and vegetative fusion of hyphae on the leaf surface. This epiphytic colonisation has been investigated very little during interactions in which an isolate of Z. tritici is unable to colonise the apoplast, as occurs during avirulence. However, avirulent isolates have been seen to undergo sexual crosses in the absense of leaf penetration, and it is widely accepted that the main point of distinction between virulent and avirulent isolates occurs at the point of attempted leaf penetration or attempted apoplastic growth, which fails in the avirulent case. In this work, we describe extensive epiphytic growth in three isolates which are unable or have very limited ability to invade the leaf, and show that growth form is as variable as for fully virulent isolates. We demonstrate that during certain interactions, Z. tritici isolates rarely invade the leaf and form pycnidia, but induce necrosis. These isolates are able to achieve higher epiphytic biomass than fully virulent isolates during asymptomatic growth, and may undergo very extensive asexual reproduction on the leaf surface. These findings have implications for open questions such as whether and how Z. tritici obtains nutrients on the leaf surface and the nature of its interaction with wheat defences.

H4K20me3 is important for Ash1-mediated H3K36me3 and transcriptional silencing in facultative heterochromatin in a fungal pathogen

Facultative heterochromatin controls development and differentiation in many eukaryotes. In metazoans, plants, and many filamentous fungi, facultative heterochromatin is characterized by transcriptional repression and enrichment with nucleosomes that are trimethylated at histone H3 lysine 27 (H3K27me3). While loss of H3K27me3 results in derepression of transcriptional gene silencing in many species, additional up- and downstream layers of regulation are necessary to mediate control of transcription in chromosome regions enriched with H3K27me3. Here, we investigated the effects of one histone mark on histone H4, namely H4K20me3, in the fungus Zymoseptoria tritici, a globally important pathogen of wheat. Deletion of kmt5, the gene encoding the sole methyltransferase responsible for H4K20 methylation, resulted in global derepression of transcription, especially in regions of facultative heterochromatin. Derepression in the absence of H4K20me3 not only affected known genes but also a large number of novel, previously undetected transcripts generated from regions of facultative heterochromatin on accessory chromosomes. Transcriptional activation in kmt5 deletion strains was accompanied by a complete loss of Ash1-mediated H3K36me3 and chromatin reorganization affecting H3K27me3 and H3K4me2 distribution in regions of facultative heterochromatin. Strains with H4K20L, M or Q mutations in the single histone H4 gene of Z. tritici recapitulated these chromatin changes, suggesting that H4K20me3 is important for Ash1-mediated H3K36me3. The ∆kmt5 mutants we obtained were more sensitive to genotoxic stressors than wild type and both, ∆kmt5 and ∆ash1, showed greatly increased rates of accessory chromosome loss. Taken together, our results provide insights into an unsuspected mechanism involved in the assembly and maintenance of facultative heterochromatin.

Quantitative and qualitative plant-pathogen interactions call upon similar pathogenicity genes with a spectrum of effects

Septoria leaf blotch is a foliar wheat disease controlled by a combination of plant genetic resistances and fungicides use. R-gene-based qualitative resistance durability is limited due to gene-for-gene interactions with fungal avirulence (Avr) genes. Quantitative resistance is considered more durable but the mechanisms involved are not well documented. We hypothesize that genes involved in quantitative and qualitative plant-pathogen interactions are similar. A bi-parental population of Zymoseptoria tritici was inoculated on wheat cultivar 'Renan' and a linkage analysis performed to map QTL. Three pathogenicity QTL, Qzt-I05-1, Qzt-I05-6 and Qzt-I07-13, were mapped on chromosomes 1, 6 and 13 in Z. tritici, and a candidate pathogenicity gene on chromosome 6 was selected based on its effector-like characteristics. The candidate gene was cloned by Agrobacterium tumefaciens-mediated transformation, and a pathology test assessed the effect of the mutant strains on 'Renan'. This gene was demonstrated to be involved in quantitative pathogenicity. By cloning a newly annotated quantitative-effect gene in Z. tritici that is effector-like, we demonstrated that genes underlying pathogenicity QTL can be similar to Avr genes. This opens up the previously probed possibility that 'gene-for-gene' underlies not only qualitative but also quantitative plant-pathogen interactions in this pathosystem.

Natural variation in Avr3D1 from Zymoseptoria sp. contributes to quantitative gene-for-gene resistance and to host specificity

Successful host colonization by plant pathogens requires the circumvention of host defense responses, frequently through sequence modifications in secreted pathogen proteins known as avirulence factors (Avrs). Although Avr sequences are often polymorphic, the contribution of these polymorphisms to virulence diversity in natural pathogen populations remains largely unexplored. We used molecular genetic tools to determine how natural sequence polymorphisms of the avirulence factor Avr3D1 in the wheat pathogen Zymoseptoria tritici contributed to adaptive changes in virulence. We showed that there is a continuous distribution in the magnitude of resistance triggered by different Avr3D1 isoforms and demonstrated that natural variation in an Avr gene can lead to a quantitative resistance phenotype. We further showed that homologues of Avr3D1 in two nonpathogenic sister species of Z. tritici are recognized by some wheat cultivars, suggesting that Avr-R gene-for-gene interactions can contribute to nonhost resistance. We suggest that the mechanisms underlying host range, qualitative resistance, and quantitative resistance are not exclusive.

Zymoseptoria tritici white-collar complex integrates light, temperature and plant cues to initiate dimorphism and pathogenesis

Transitioning from spores to hyphae is pivotal to host invasion by the plant pathogenic fungus Zymoseptoria tritici. This dimorphic switch can be initiated by high temperature in vitro (~27 °C); however, such a condition may induce cellular heat stress, questioning its relevance to field infections. Here, we study the regulation of the dimorphic switch by temperature and other factors. Climate data from wheat-growing areas indicate that the pathogen sporadically experiences high temperatures such as 27 °C during summer months. However, using a fluorescent dimorphic switch reporter (FDR1) in four wild-type strains, we show that dimorphic switching already initiates at 15-18 °C, and is enhanced by wheat leaf surface compounds. Transcriptomics reveals 1261 genes that are up- or down-regulated in hyphae of all strains. These pan-strain core dimorphism genes (PCDGs) encode known effectors, dimorphism and transcription factors, and light-responsive proteins (velvet factors, opsins, putative blue light receptors). An FDR1-based genetic screen reveals a crucial role for the white-collar complex (WCC) in dimorphism and virulence, mediated by control of PCDG expression. Thus, WCC integrates light with biotic and abiotic cues to orchestrate Z. tritici infection.

The Zymoseptoria tritici white collar-1 gene, ZtWco-1, is required for development and virulence on wheat

The fungus Zymoseptoria tritici causes Septoria Tritici Blotch (STB), which is one of the most devastating diseases of wheat in Europe. There are currently no fully durable methods of control against Z. tritici, so novel strategies are urgently required. One of the ways in which fungi are able to respond to their surrounding environment is through the use of photoreceptor proteins which detect light signals. Although previous evidence suggests that Z. tritici can detect light, no photoreceptor genes have been characterised in this pathogen. This study characterises ZtWco-1, a predicted photoreceptor gene in Z. tritici. The ZtWco-1 gene is a putative homolog to the blue light photoreceptor from Neurospora crassa, wc-1. Z. tritici mutants with deletions in ZtWco-1 have defects in hyphal branching, melanisation and virulence on wheat. In addition, we identify the putative circadian clock gene ZtFrq in Z. tritici. This study provides evidence for the genetic regulation of light detection in Z. tritici and it open avenues for future research into whether this pathogen has a circadian clock.

In Planta Gene Expression Analysis and Colonization of Fusarium oxysporum

In planta gene expression analysis and GFP-based confocal microscopy are two powerful techniques that may be coupled to assess the extent and dynamics of plant colonization by a fungal pathogen. Here we describe methods to prepare common bean plants for inoculation with a highly virulent strain of Fusarium oxysporum f. sp. phaseoli, quantify the extent of colonization by RT-qPCR, and visualize the colonized tissues by confocal microscopy.

Identification and Comparison of Colletotrichum Secreted Effector Candidates Reveal Two Independent Lineages Pathogenic to Soybean

Colletotrichum is one of the most important plant pathogenic genus of fungi due to its scientific and economic impact. A wide range of hosts can be infected by Colletotrichum spp., which causes losses in crops of major importance worldwide, such as soybean. Soybean anthracnose is mainly caused by C. truncatum, but other species have been identified at an increasing rate during the last decade, becoming one of the most important limiting factors to soybean production in several regions. To gain a better understanding of the evolutionary origin of soybean anthracnose, we compared the repertoire of effector candidates of four Colletotrichum species pathogenic to soybean and eight species not pathogenic. Our results show that the four species infecting soybean belong to two lineages and do not share any effector candidates. These results strongly suggest that two Colletotrichum lineages have acquired the capability to infect soybean independently. This study also provides, for each lineage, a set of candidate effectors encoding genes that may have important roles in pathogenicity towards soybean offering a new resource useful for further research on soybean anthracnose management.

Remote homology clustering identifies lowly conserved families of effector proteins in plant-pathogenic fungi

Plant diseases caused by fungal pathogens are typically initiated by molecular interactions between 'effector' molecules released by a pathogen and receptor molecules on or within the plant host cell. In many cases these effector-receptor interactions directly determine host resistance or susceptibility. The search for fungal effector proteins is a developing area in fungal-plant pathology, with more than 165 distinct confirmed fungal effector proteins in the public domain. For a small number of these, novel effectors can be rapidly discovered across multiple fungal species through the identification of known effector homologues. However, many have no detectable homology by standard sequence-based search methods. This study employs a novel comparison method (RemEff) that is capable of identifying protein families with greater sensitivity than traditional homology-inference methods, leveraging a growing pool of confirmed fungal effector data to enable the prediction of novel fungal effector candidates by protein family association. Resources relating to the RemEff method and data used in this study are available from https://figshare.com/projects/Effector_protein_remote_homology/87965.

Induction of Wheat Resistance to STB by the Endophytic Fungus serendipita Indica and pseudomonas Protegens

Background:

Septoria tritici blotch (STB) caused by fungus Zymoseptoria tritici, is one of the important wheat (Triticum aestivum L.) diseases difficult to control because of the lack of wheat resistant cultivars. The use of biological control agents is one possible way for triggering host plant resistance to biotic and abiotic stresses

Objective:

In this study, we examined the ability of Serendipita indica and Pseudomonas protegens CHA0-mCherry in inducing the local wheat cultivar Tajan resistance to STB.

Materials and Methods:

The interaction between biological control agents and the roots of wheat was evaluated. The experiment was conducted in a completely randomized design by three replicates. Spore suspension was supplied at concentrations of 10⁷ and 10⁹ for S. indica and bacteria isolate (CHA0-mCherry) respectively. Five treatments were applied including S. indica, CHA0-mCherry, S. indica and CHA0-mCherry co-inoculation, positive and negative control. Twenty-one days after inoculation, the interaction between biological agents and plant roots were evaluated through morphological traits and qPCR. The plant resistance, disease severity, and the correlation between resistance and disease severity were assessed. Pycnidial variation and agronomic traits were also evaluated.

Results:

Twenty-one days after inoculation, both biological agents clearly colonized all treated roots of all treatments except in control plants as demonstrated by qPCR analysis. Chlamydospores were observed in the S. indica-treated hosts with the CHA0-mCherry colonizing assessment showing 5×109 CFU g-1 in the root. The asexual phase of the fungal pathogen, pycnidial diameter, was reduced in S. indica treated plants more considerably than in the other treatments. There was a positive correlation between resistance and disease severity mean when calculated by Pearson’s correlation. There was a significant difference between the root length, fresh, and dry weight of root. Spore density was inversely correlated to resistance and disease severity, when compared with control, with CHA0-mCherry being the most effective in reducing the spore density. S. indica was the most effective in promoting root growth and stem biomass, when compared with control.

Conclusions:

Serendipita indica and Pseudomonas protegens CHA0-mCherry colonies showed a potential biological control activity and efficiently enhanced the plant resistance to Z. tritici in the treated wheat roots. The microbial biological control agents are very practical in crop protection against plant disease and can be very useful in sustainable agriculture.

Abbreviations:

PLSN: percentage of leave surface necrosis, DPI: day past inoculation, PLACL: percentage of leaf area covered by lesions, PPMLA: pycnidia per millimeter in leaf area.

Dissecting the Biology of the Fungal Wheat Pathogen Zymoseptoria tritici: A Laboratory Workflow

The fungus Zymoseptoria tritici is one of the most devastating pathogens of wheat. Aside from its importance as a disease-causing agent, this species has emerged as a powerful model system for evolutionary genetic studies of crop-infecting fungal pathogens. Z. tritici exhibits exceptionally high levels of genetic and phenotypic diversity as well as morphological plasticity, which can make experimental studies and comparability of results obtained in different laboratories, e.g., from infection assays, challenging. Therefore, standardized experimental methods are crucial for research on Z. tritici biology and the interaction of this fungus with its wheat host. Here, we describe a suite of well-tested and optimized protocols ranging from isolation of Z. tritici field specimens to analyses of virulence assays under controlled conditions. Several biological and technical aspects of working with Z. tritici under laboratory conditions are considered and carefully described in each protocol. © 2020 The Authors. Basic Protocol 1: Purification of Z. tritici field isolates from leaf material Basic Protocol 2: Molecular identification of Z. tritici isolates Support Protocol 1: Rapid extraction of Z. tritici genomic DNA Support Protocol 2: Extraction of high-quality Z. tritici genomic DNA Basic Protocol 3: In vitro culture and long-term storage of Z. tritici isolates Basic Protocol 4: Analysis of Z. tritici virulence in wheat Support Protocol 3: Preparation of Z. tritici inoculum.

The identification of a transposon affecting the asexual reproduction of the wheat pathogen Zymoseptoria tritici

Zymoseptoria tritici, the causal agent of Septoria tritici blotch, is a fungal wheat pathogen that causes significant global yield losses. Within Z. tritici populations, quantitative differences in virulence among different isolates are commonly observed; however, the genetic components that underpin these differences remain elusive. In this study, intraspecific comparative transcriptomic analysis was used to identify candidate genes that contribute to differences in virulence on the wheat cultivar WW2449. This led to the identification of a multicopy gene that was not expressed in the high-virulence isolate when compared to the medium- and low-virulence isolates. Further investigation suggested this gene resides in a 7.9-kb transposon. Subsequent long-read sequencing of the isolates used in the transcriptomic analysis confirmed that this gene did reside in an active Class II transposon, which is composed of four genes named REP9-1 to -4. Silencing and overexpression of REP9-1 in two distinct genetic backgrounds demonstrated that its expression alone reduces the number of pycnidia produced by Z. tritici during infection. The REP9-1 gene identified within a Class II transposon is the first discovery of a gene in a transposable element that influences the virulence of Z. tritici. This discovery adds further complexity to genetic loci that contribute to quantitative virulence in this important pathogen.

Ago1 Affects the Virulence of the Fungal Plant Pathogen Zymoseptoria tritici

In host-pathogen interactions RNA interference (RNAi) has emerged as a pivotal mechanism to modify both, the immune responses of the host as well as the pathogenicity and virulence of the pathogen. In addition, in some fungi RNAi is also known to affect chromosome biology via its effect on chromatin conformation. Previous studies reported no effect of the RNAi machinery on the virulence of the fungal plant pathogen Zymoseptoria tritici however the role of RNAi is still poorly understood in this species. Herein, we elucidate whether the RNAi machinery is conserved within the genus Zymoseptoria. Moreover, we conduct functional analyses of Argonaute and Dicer-like proteins and test if the RNAi machinery affects chromosome stability. We show that the RNAi machinery is conserved among closely related Zymoseptoria species while an exceptional pattern of allelic diversity was possibly caused by introgression. The deletion of Ago1 reduced the ability of the fungus to produce asexual propagules in planta in a quantitative matter. Chromosome stability of the accessory chromosome of Z. tritici was not prominently affected by the RNAi machinery. These results indicate, in contrast to previous finding, a role of the RNAi pathway during host infection, but not in the stability of accessory chromosomes in Z. tritici.

CgNPG1 as a Novel Pathogenic Gene of Colletotrichum gloeosporioides From Hevea brasiliensis in Mycelial Growth, Conidiation, and the Invasive Structures Development

The rubber tree (Hevea brasiliensis) is a tropical perennial crop for the primary source of natural rubber. Colletotrichum gloeosporioides from Hevea brasiliensis (C. gloeosporioides Hb) and Colletotrichum acutatum from Hevea brasiliensis (C. acutatum Hb) are the causal agents of rubber tree anthracnose and lead to serious loss of natural rubber production. Inoculation tests showed that C. gloeosporioides Hb possessed higher pathogenicity than C. acutatum Hb to the rubber tree. Genomic analysis revealed that an unknown gene, named CgNPG1 (a Novel Pathogenic Gene 1), was presented in the genome of C. gloeosporioides Hb but not identified in C. acutatum Hb. CgNPG1 was predicted to encode a small secretory protein without any conserved domain. To investigate the functions of CgNPG1 in C. gloeosporioides Hb and in C. acutatum Hb, the gene deletion and overexpression mutants were generated. The phenotype analysis showed that deletion of CgNPG1 led to changed conidia morphology, decreased mycelial growth, conidiation, conidia germination rate, appressorium formation rate, and pathogenicity of C. gloeosporioides Hb to the rubber tree. Meanwhile, heterogeneous expression of CgNPG1 in C. acutatum Hb significantly changed the conidia morphology and improved the mycelial growth rate, conidiation, conidia germination rate, appressorium formation rate, and the pathogenicity of C. acutatum Hb to the rubber tree. Consistently, CgNPG1 increased the expression level of CaCRZ1 and CaCMK1 in C. acutatum Hb. These data suggested that CgNPG1 contributed to mycelial growth, conidiation, the development of invasive structures, and the pathogenicity of Colletotrichum to the rubber tree, which might be related to the modulation of CaCRZ1 and mitogen-activated protein kinase CMK1. Our results provided new insight into CgNPG1 in regulating growth and pathogenicity of the Colletotrichum spp.

Recent loss of the Dim2 DNA methyltransferase decreases mutation rate in repeats and changes evolutionary trajectory in a fungal pathogen

DNA methylation is found throughout all domains of life, yet the extent and function of DNA methylation differ among eukaryotes. Strains of the plant pathogenic fungus Zymoseptoria tritici appeared to lack cytosine DNA methylation (5mC) because gene amplification followed by Repeat-Induced Point mutation (RIP) resulted in the inactivation of the dim2 DNA methyltransferase gene. 5mC is, however, present in closely related sister species. We demonstrate that inactivation of dim2 occurred recently as some Z. tritici isolates carry a functional dim2 gene. Moreover, we show that dim2 inactivation occurred by a different path than previously hypothesized. We mapped the genome-wide distribution of 5mC in strains with or without functional dim2 alleles. Presence of functional dim2 correlates with high levels of 5mC in transposable elements (TEs), suggesting a role in genome defense. We identified low levels of 5mC in strains carrying non-functional dim2 alleles, suggesting that 5mC is maintained over time, presumably by an active Dnmt5 DNA methyltransferase. Integration of a functional dim2 allele in strains with mutated dim2 restored normal 5mC levels, demonstrating de novo cytosine methylation activity of Dim2. To assess the importance of 5mC for genome evolution, we performed an evolution experiment, comparing genomes of strains with high levels of 5mC to genomes of strains lacking functional dim2. We found that presence of a functional dim2 allele alters nucleotide composition by promoting C to T transitions (C→T) specifically at CpA (CA) sites during mitosis, likely contributing to TE inactivation. Our results show that 5mC density at TEs is a polymorphic trait in Z. tritici populations that can impact genome evolution.

Whole-Genome Comparisons of Ergot Fungi Reveals the Divergence and Evolution of Species within the Genus Claviceps Are the Result of Varying Mechanisms Driving Genome Evolution and Host Range Expansion

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

Population Genomics of Fungal Plant Pathogens and the Analyses of Rapidly Evolving Genome Compartments

Genome sequencing of fungal pathogens have documented extensive variation in genome structure and composition between species and in many cases between individuals of the same species. This type of genomic variation can be adaptive for pathogens to rapidly evolve new virulence phenotypes. Analyses of genome-wide variation in fungal pathogen genomes rely on high quality assemblies and methods to detect and quantify structural variation. Population genomic studies in fungi have addressed the underlying mechanisms whereby structural variation can be rapidly generated. Transposable elements, high mutation and recombination rates as well as incorrect chromosome segregation during mitosis and meiosis contribute to extensive variation observed in many species. We here summarize key findings in the field of fungal pathogen genomics and we discuss methods to detect and characterize structural variants including an alignment-based pipeline to study variation in population genomic data.

Evolution and Adaptation of Forest and Crop Pathogens in the Anthropocene

Anthropocene marks the era when human activity is making a significant impact on earth, its ecological and biogeographical systems. The domestication and intensification of agricultural and forest production systems have had a large impact on plant and tree health. Some pathogens benefitted from these human activities and have evolved and adapted in response to the expansion of crop and forest systems, resulting in global outbreaks. Global pathogen genomics data including population genomics and high-quality reference assemblies are crucial for understanding the evolution and adaptation of pathogens. Crops and forest trees have remarkably different characteristics, such as reproductive time and the level of domestication. They also have different production systems for disease management with more intensive management in crops than forest trees. By comparing and contrasting results from pathogen population genomic studies done on widely different agricultural and forest production systems, we can improve our understanding of pathogen evolution and adaptation to different selection pressures. We find that in spite of these differences, similar processes such as hybridization, host jumps, selection, specialization, and clonal expansion are shaping the pathogen populations in both crops and forest trees. We propose some solutions to reduce these impacts and lower the probability of global pathogen outbreaks so that we can envision better management strategies to sustain global food production as well as ecosystem services.

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.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

A multi-omics approach to solving problems in plant disease ecology

The swift rise of omics-approaches allows for investigating microbial diversity and plant-microbe interactions across diverse ecological communities and spatio-temporal scales. The environment, however, is rapidly changing. The introduction of invasive species and the effects of climate change have particular impact on emerging plant diseases and managing current epidemics. It is critical, therefore, to take a holistic approach to understand how and why pathogenesis occurs in order to effectively manage for diseases given the synergies of changing environmental conditions. A multi-omics approach allows for a detailed picture of plant-microbial interactions and can ultimately allow us to build predictive models for how microbes and plants will respond to stress under environmental change. This article is designed as a primer for those interested in integrating -omic approaches into their plant disease research. We review -omics technologies salient to pathology including metabolomics, genomics, metagenomics, volatilomics, and spectranomics, and present cases where multi-omics have been successfully used for plant disease ecology. We then discuss additional limitations and pitfalls to be wary of prior to conducting an integrated research project as well as provide information about promising future directions.

Meiotic recombination in the offspring of Microbotryum hybrids and its impact on pathogenicity

BACKGROUND

Hybridization is a central mechanism in evolution, producing new species or introducing important genetic variation into existing species. In plant-pathogenic fungi, adaptation and specialization to exploit a host species are key determinants of evolutionary success. Here, we performed experimental crosses between the two pathogenic Microbotryum species, M. lychnidis-dioicae and M. silenes-acaulis that are specialized to different hosts. The resulting offspring were analyzed on phenotypic and genomic levels to describe genomic characteristics of hybrid offspring and genetic factors likely involved in host-specialization.

RESULTS

Genomic analyses of interspecific fungal hybrids revealed that individuals were most viable if the majority of loci were inherited from one species. Interestingly, species-specific loci were strictly controlled by the species' origin of the mating type locus. Moreover we detected signs of crossing over and chromosome duplications in the genomes of the analyzed hybrids. In Microbotryum, mitochondrial DNA was found to be uniparentally inherited from the a2 mating type. Genome comparison revealed that most gene families are shared and the majority of genes are conserved between the two species, indicating very similar biological features, including infection and pathogenicity processes. Moreover, we detected 211 candidate genes that were retained under host-driven selection of backcrossed lines. These genes and might therefore either play a crucial role in host specialization or be linked to genes that are essential for specialization.

CONCLUSION

The combination of genome analyses with experimental selection and hybridization is a promising way to investigate host-pathogen interactions. This study manifests genetic factors of host specialization that are required for successful biotrophic infection of the post-zygotic stage, but also demonstrates the strong influence of intra-genomic conflicts or instabilities on the viability of hybrids in the haploid host-independent stage.

Rapid evolution in plant-microbe interactions - an evolutionary genomics perspective

Access to greater genomic resolution through new sequencing technologies is transforming the field of plant pathology. As scientists embrace these new methods, some overarching patterns and observations come into focus. Evolutionary genomic studies are used to determine not only the origins of pathogen lineages and geographic patterns of genetic diversity, but also to discern how natural selection structures genetic variation across the genome. With greater and greater resolution, we can now pinpoint the targets of selection on a large scale. At multiple levels, crypsis and convergent evolution are evident. Host jumps and shifts may be more pervasive than once believed, and hybridization and horizontal gene transfer (HGT) likely play important roles in the emergence of genetic novelty.

A fungal pathogen induces systemic susceptibility and systemic shifts in wheat metabolome and microbiome composition

Yield losses caused by fungal pathogens represent a major threat to global food production. One of the most devastating fungal wheat pathogens is Zymoseptoria tritici. Despite the importance of this fungus, the underlying mechanisms of plant–pathogen interactions are poorly understood. Here we present a conceptual framework based on coinfection assays, comparative metabolomics, and microbiome profiling to study the interaction of Z. tritici in susceptible and resistant wheat. We demonstrate that Z. tritici suppresses the production of immune-related metabolites in a susceptible cultivar. Remarkably, this fungus-induced immune suppression spreads within the leaf and even to other leaves, a phenomenon that we term “systemic induced susceptibility”. Using a comparative metabolomics approach, we identify defense-related biosynthetic pathways that are suppressed and induced in susceptible and resistant cultivars, respectively. We show that these fungus-induced changes correlate with changes in the wheat leaf microbiome. Our findings suggest that immune suppression by this hemibiotrophic pathogen impacts specialized plant metabolism, alters its associated microbial communities, and renders wheat vulnerable to further infections.

The fungal plant pathogen Zymoseptoria tritici is a major threat to wheat yield. Here Seybold et al. show that Z. tritici can suppress immune responses not only in infected tissue but also on other leaves, a phenomenon termed “systemic induced susceptibility” that is correlated with systemic changes in metabolite accumulation.

Core endophyte communities of different citrus varieties from citrus growing regions in China

The native microbiomes of citrus trees play important roles in plant health, with good communication between the native microbiome and the host plant. Here, we report on the native endophytes in 24 citrus varieties in nine citrus growing regions in China; some of the trees were healthy and others had asymptomatic or symptomatic huanglongbing, which is caused by the pathogen Candidatus Liberibacter asiaticus (CLas). We used culture-dependent analysis and characterized the isolates by partial 16S rRNA gene sequencing. The endophytes were compared between different citrus varieties, regions, and disease states (healthy, asymptomatic, and symptomatic). The total number of endophytes isolated from most of the citrus varieties was 104-106 CFU/g of leaves, but it differed significantly by disease state, with the highest numbers in the healthy leaves and the lowest in the symptomatic leaves (p < 0.05). Among the citrus varieties, the Valencia variety had the maximum number of endophyte species (22). The most dominant endophytes were Bacillus subtilis, B. velezensis, Curtobacterium luteum, and Microbacterium testaceum. The higher frequency of B. subtilis in the healthy/asymptomatic plants compared to the symptomatic plants suggests that it has a role in huanglongbing resistance. Native endophyte communities in various citrus varieties could be used to improve citrus growth and combat CLas.

The transcription factor Zt107320 affects the dimorphic switch, growth and virulence of the fungal wheat pathogen Zymoseptoria tritici

Zymoseptoria tritici is a filamentous fungus causing Septoria tritici blotch in wheat. The pathogen has a narrow host range and infections of grasses other than susceptible wheat are blocked early after stomatal penetration. During these abortive infections, the fungus shows a markedly different gene expression pattern. However, the underlying mechanisms causing differential gene expression during host and non-host interactions are largely unknown, but likely include transcriptional regulators responsible for the onset of an infection programme in compatible hosts. MoCOD1, a member of the fungal Zn(II)₂ Cys₆ transcription factor family, has been shown to directly affect pathogenicity in the rice blast pathogen Magnaporthe oryzae. Here, we analyse the role of the putative transcription factor Zt107320, a homologue of MoCOD1, during infection of compatible and incompatible hosts by Z. tritici. We show for the first time that Zt107320 is differentially expressed in host versus non-host infections and that lower expression corresponds to an incompatible infection of non-hosts. Applying reverse genetics approaches, we further show that Zt107320 regulates the dimorphic switch as well as the growth rate of Z. tritici and affects fungal cell wall composition in vitro. Moreover, ∆Zt107320 mutants showed reduced virulence during compatible infections of wheat. We conclude that Zt107320 directly influences pathogen fitness and propose that Zt107320 is involved in the regulation of growth processes and pathogenicity during infection.

Taxonomically Restricted Wheat Genes Interact With Small Secreted Fungal Proteins and Enhance Resistance to Septoria Tritici Blotch Disease

Understanding the nuances of host/pathogen interactions are paramount if we wish to effectively control cereal diseases. In the case of the wheat/Zymoseptoria tritici interaction that leads to Septoria tritici blotch (STB) disease, a 10,000-year-old conflict has led to considerable armaments being developed on both sides which are not reflected in conventional model systems. Taxonomically restricted genes (TRGs) have evolved in wheat to better allow it to cope with stress caused by fungal pathogens, and Z. tritici has evolved specialized effectors which allow it to manipulate its' host. A microarray focused on the latent phase response of a resistant wheat cultivar (cv. Stigg) and susceptible wheat cultivar (cv. Gallant) to Z. tritici infection was mined for TRGs within the Poaceae. From this analysis, we identified two TRGs that were significantly upregulated in response to Z. tritici infection, Septoria-responsive TRG6 and 7 (TaSRTRG6 and TaSRTRG7). Virus induced silencing of these genes resulted in an increased susceptibility to STB disease in cvs. Gallant and Stigg, and significantly so in the latter (2.5-fold increase in STB disease). In silico and localization studies categorized TaSRTRG6 as a secreted protein and TaSRTRG7 as an intracellular protein. Yeast two-hybrid analysis and biofluorescent complementation studies demonstrated that both TaSRTRG6 and TaSRTRG7 can interact with small proteins secreted by Z. tritici (potential effector candidates). Thus we conclude that TRGs are an important part of the wheat-Z. tritici co-evolution story and potential candidates for modulating STB resistance.

The Current Epidemic of the Barley Pathogen Ramularia collo-cygni Derives from a Population Expansion and Shows Global Admixture

Ramularia leaf spot is becoming an ever-increasing problem in main barley-growing regions since the 1980s, causing up to 70% yield loss in extreme cases. Yet, the causal agent Ramularia collo-cygni remains poorly studied. The diversity of the pathogen in the field thus far remains unknown. Furthermore, it is unknown to what extent the pathogen has a sexual reproductive cycle. The teleomorph of R. collo-cygni has not been observed. To study the genetic diversity of R. collo-cygni and get more insights in its biology, we sequenced the genomes of 19 R. collo-cygni isolates from multiple geographic locations and diverse hosts. Nucleotide polymorphism analyses of all isolates shows that R. collo-cygni is genetically diverse worldwide, with little geographic or host specific differentiation. Next, we used two different methods to detect signals of recombination in our sample set. Both methods find putative recombination events, which indicate that sexual reproduction happens or has happened in the global R. collo-cygni population. Lastly, we used these data on recombination to perform historic population size analyses. These suggest that the effective population size of R. collo-cygni decreased during the domestication of barley and subsequently grew with the rise of agriculture. Our findings deepen our understanding of R. collo-cygni biology and can help us to understand the current epidemic. We discuss how our findings support possible global spread through seed transfer, and we highlight how recombination, clonal spreading, and lack of host specificity could amplify global epidemics of this increasingly important disease and suggest specific approaches to combat the pathogen.

Dynamic QTL for adult plant resistance to powdery mildew in common wheat (Triticum aestivum L.)

Agriculture will benefit from a rigorous characterization of genes for adult plant resistance (APR) since this gene class was recognized to provide more durable protection from plant diseases. The present study reports the identification of APR loci to powdery mildew in German winter wheat cultivars Cortez and Atlantis. Cortez was previously shown to carry all-stage resistance gene Pm3e. To avoid interference of Pm3e in APR studies, line 6037 that lacked Pm3e but showed field resistance from doubled-haploid (DH) population Atlantis/Cortez was used in two backcrosses to Atlantis for the establishment of DH population 6037/Atlantis//Atlantis. APR was assessed in the greenhouse 10, 15, and 20 days after inoculation (dai) from the 4-leaf stage onwards and combined with single-nucleotide polymorphism data in a genome-wide association study (GWAS) and a linkage map-based quantitative trait loci (QTL) analysis. In GWAS, two QTL were detected: one on chromosome 1BL 10 dai, the other on chromosome 2BL 20 dai. In conventional QTL analysis, both QTL were detected with all three disease ratings: the QTL on chromosome 1BL explained a maximum of 35.2% of the phenotypic variation 10 dai, whereas the QTL on chromosome 2BL explained a maximum of 43.5% of the phenotypic variation 20 dai. Compared with GWAS, linkage map-based QTL analysis allowed following the dynamics of QTL action. The two large-effect QTL for APR to powdery mildew with dynamic gene action can be useful for the enhancement of wheat germplasm.

What determines host specificity in hyperspecialized plant parasitic nematodes?

Background

In hyperspecialized parasites, the ability to grow on a particular host relies on specific virulence factors called effectors. These excreted proteins are involved in the molecular mechanisms of parasitism and distinguish virulent pathogens from non-virulent related species. The potato cyst nematodes (PCN) Globodera rostochiensis and G. pallida are major plant-parasitic nematodes developing on numerous solanaceous species including potato. Their close relatives, G. tabacum and G. mexicana are stimulated by potato root diffusate but unable to establish a feeding site on this plant host.

Results

RNA sequencing was used to characterize transcriptomic differences among these four Globodera species and to identify genes associated with host specificity. We identified seven transcripts that were unique to PCN species, including a protein involved in ubiquitination. We also found 545 genes that were differentially expressed between PCN and non-PCN species, including 78 genes coding for effector proteins, which represent more than a 6-fold enrichment compared to the whole transcriptome. Gene polymorphism analysis identified 359 homozygous non-synonymous variants showing a strong evidence for selection in PCN species.

Conclusions

Overall, we demonstrated that the determinant of host specificity resides in the regulation of essential effector gene expression that could be under the control of a single or of very few regulatory genes. Such genes are therefore promising targets for the development of novel and more sustainable resistances against potato cyst nematodes.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5853-4) contains supplementary material, which is available to authorized users.

Destabilization of chromosome structure by histone H3 lysine 27 methylation

Chromosome and genome stability are important for normal cell function as instability often correlates with disease and dysfunction of DNA repair mechanisms. Many organisms maintain supernumerary or accessory chromosomes that deviate from standard chromosomes. The pathogenic fungus Zymoseptoria tritici has as many as eight accessory chromosomes, which are highly unstable during meiosis and mitosis, transcriptionally repressed, show enrichment of repetitive elements, and enrichment with heterochromatic histone methylation marks, e.g., trimethylation of H3 lysine 9 or lysine 27 (H3K9me3, H3K27me3). To elucidate the role of heterochromatin on genome stability in Z. tritici, we deleted the genes encoding the methyltransferases responsible for H3K9me3 and H3K27me3, kmt1 and kmt6, respectively, and generated a double mutant. We combined experimental evolution and genomic analyses to determine the impact of these deletions on chromosome and genome stability, both in vitro and in planta. We used whole genome sequencing, ChIP-seq, and RNA-seq to compare changes in genome and chromatin structure, and differences in gene expression between mutant and wildtype strains. Analyses of genome and ChIP-seq data in H3K9me3-deficient strains revealed dramatic chromatin reorganization, where H3K27me3 is mostly relocalized into regions that are enriched with H3K9me3 in wild type. Many genome rearrangements and formation of new chromosomes were found in the absence of H3K9me3, accompanied by activation of transposable elements. In stark contrast, loss of H3K27me3 actually increased the stability of accessory chromosomes under normal growth conditions in vitro, even without large scale changes in gene activity. We conclude that H3K9me3 is important for the maintenance of genome stability because it disallows H3K27me3 in regions considered constitutive heterochromatin. In this system, H3K27me3 reduces the overall stability of accessory chromosomes, generating a "metastable" state for these quasi-essential regions of the genome.

What makes a specialized endophyte special?

Fungal plant symbionts can be highly specialized on a limited range of host genotypes and species. Understanding the genetic basis of this specialization, the mechanisms governing its establishment and the relationship between specialization and speciation is a major challenge for evolutionary biologists (Timms & Read, ). A deeper knowledge of evolutionary plant-microbe interactions could be exploited to improve agricultural management, by bringing fungal biodiversity and fungal biomass under greater and more durable human control. Previous studies on pathogens have shown that effectors, that is, small secreted proteins that modulate plant physiology to favour host colonization, play a key role in infection of novel hosts (e.g., Inoue et al., ) or in host specialization (e.g., Liao et al. ()). Like pathogens, endophytes also manipulate the physiology of their hosts and colonize novel hosts to which they specialize (Hardoim et al., ). These biological characteristics of endophytes raise the question of similarities in the protein arsenal contributing to the specialization of pathogens and endophytes. In this issue of Molecular Ecology, Schirrmann et al. () used a combination of divergence genome scans and tests for positive selection to investigate the genetic basis of specialization of two subspecies of the symbiont Epichloë typhina occurring on two different grass hosts. Their analyses suggest a key role of effectors as determinants of host specialization. This study paves the way towards the comparative analysis of the genomics of speciation among plant symbionts.

Highly flexible infection programs in a specialized wheat pathogen

Many filamentous plant pathogens exhibit high levels of genomic variability, yet the impact of this variation on host-pathogen interactions is largely unknown. We have addressed host specialization in the wheat pathogen Zymoseptoria tritici. Our study builds on comparative analyses of infection and gene expression phenotypes of three isolates and reveals the extent to which genomic variation translates into phenotypic variation. The isolates exhibit genetic and genomic variation but are similarly virulent. By combining confocal microscopy, disease monitoring, staining of ROS, and comparative transcriptome analyses, we conducted a detailed comparison of the infection processes of these isolates in a susceptible wheat cultivar. We characterized four core infection stages: establishment, biotrophic growth, lifestyle transition, and necrotrophic growth and asexual reproduction that are shared by the three isolates. However, we demonstrate differentiated temporal and spatial infection development and significant differences in the expression profiles of the three isolates during the infection stages. More than 20% of the genes were differentially expressed and these genes were located significantly closer to transposable elements, suggesting an impact of epigenetic regulation. Further, differentially expressed genes were enriched in effector candidates suggesting that isolate-specific strategies for manipulating host defenses are present in Z. tritici. We demonstrate that individuals of a host-specialized pathogen have highly differentiated infection programs characterized by flexible infection development and functional redundancy. This illustrates how high genetic diversity in pathogen populations results in highly differentiated infection phenotypes, which fact needs to be acknowledged to understand host-pathogen interactions and pathogen evolution.

A New Reference Genome Shows the One-Speed Genome Structure of the Barley Pathogen Ramularia collo-cygni

Ramularia leaf spot has recently emerged as a major threat to barley production world-wide, causing 25% yield loss in many barley growing regions. Here, we provide a new reference genome of the causal agent, the Dothideomycete Ramularia collo-cygni. The assembly of 32 Mb consists of 78 scaffolds. We used RNA-seq to identify 11,622 genes of which 1,303 and 282 are coding for predicted secreted proteins and putative effectors respectively.

The pathogen separated from its nearest sequenced relative, Zymoseptoria tritici ∼27 Ma. We calculated the divergence of the two species on protein level and see remarkably high synonymous and nonsynonymous divergence. Unlike in many other plant pathogens, the comparisons of transposable elements and gene distributions, show a very homogeneous genome for R. collo-cygni. We see no evidence for higher selective pressure on putative effectors or other secreted proteins and repetitive sequences are spread evenly across the scaffolds. These findings could be associated to the predominantly endophytic life-style of the pathogen. We hypothesize that R. collo-cygni only recently became pathogenic and that therefore its genome does not yet show the typical pathogen characteristics. Because of its high scaffold length and improved CDS annotations, our new reference sequence provides a valuable resource for the community for future comparative genomics and population genetics studies.

The Genome Biology of Effector Gene Evolution in Filamentous Plant Pathogens

Filamentous pathogens, including fungi and oomycetes, pose major threats to global food security. Crop pathogens cause damage by secreting effectors that manipulate the host to the pathogen's advantage. Genes encoding such effectors are among the most rapidly evolving genes in pathogen genomes. Here, we review how the major characteristics of the emergence, function, and regulation of effector genes are tightly linked to the genomic compartments where these genes are located in pathogen genomes. The presence of repetitive elements in these compartments is associated with elevated rates of point mutations and sequence rearrangements with a major impact on effector diversification. The expression of many effectors converges on an epigenetic control mediated by the presence of repetitive elements. Population genomics analyses showed that rapidly evolving pathogens show high rates of turnover at effector loci and display a mosaic in effector presence-absence polymorphism among strains. We conclude that effective pathogen containment strategies require a thorough understanding of the effector genome biology and the pathogen's potential for rapid adaptation.

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.

A fungal avirulence factor encoded in a highly plastic genomic region triggers partial resistance to septoria tritici blotch

Cultivar-strain specificity in the wheat-Zymoseptoria tritici pathosystem determines the infection outcome and is controlled by resistance genes on the host side, many of which have been identified. On the pathogen side, however, the molecular determinants of specificity remain largely unknown. We used genetic mapping, targeted gene disruption and allele swapping to characterise the recognition of the new avirulence factor Avr3D1. We then combined population genetic and comparative genomic analyses to characterise the evolutionary trajectory of Avr3D1. Avr3D1 is specifically recognised by wheat cultivars harbouring the Stb7 resistance gene, triggering a strong defence response without preventing pathogen infection and reproduction. Avr3D1 resides in a cluster of putative effector genes located in a genome region populated by independent transposable element insertions. The gene was present in all 132 investigated strains and is highly polymorphic, with 30 different protein variants identified. We demonstrated that specific amino acid substitutions in Avr3D1 led to evasion of recognition. These results demonstrate that quantitative resistance and gene-for-gene interactions are not mutually exclusive. Localising avirulence genes in highly plastic genomic regions probably facilitates accelerated evolution that enables escape from recognition by resistance proteins.

Genome-wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization

Host specialization is a key process in ecological divergence and speciation of plant-associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass-endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole-genome sequencing data, we combined F_ST and D_XY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant-fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.

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.

Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici

We conducted a comprehensive analysis of virulence in the fungal wheat pathogen Zymoseptoria tritici using quantitative trait locus (QTL) mapping. High-throughput phenotyping based on automated image analysis allowed the measurement of pathogen virulence on a scale and with a precision that was not previously possible. Across two mapping populations encompassing more than 520 progeny, 540 710 pycnidia were counted and their sizes and grey values were measured. A significant correlation was found between pycnidia size and both spore size and number. Precise measurements of percentage leaf area covered by lesions provided a quantitative measure of host damage. Combining these large and accurate phenotypic datasets with a dense panel of restriction site-associated DNA sequencing (RADseq) genetic markers enabled us to genetically dissect pathogen virulence into components related to host damage and those related to pathogen reproduction. We showed that different components of virulence can be under separate genetic control. Large- and small-effect QTLs were identified for all traits, with some QTLs specific to mapping populations, cultivars and traits and other QTLs shared among traits within the same mapping population. We associated the presence of four accessory chromosomes with small, but significant, increases in several virulence traits, providing the first evidence for a meaningful function associated with accessory chromosomes in this organism. A large-effect QTL involved in host specialization was identified on chromosome 7, leading to the identification of candidate genes having a large effect on virulence.

Forward Genetics Approach Reveals Host Genotype-Dependent Importance of Accessory Chromosomes in the Fungal Wheat Pathogen Zymoseptoria tritici

The fungal wheat pathogen Zymoseptoria tritici possesses a large complement of accessory chromosomes showing presence/absence polymorphism among isolates. These chromosomes encode hundreds of genes; however, their functional role and why the chromosomes have been maintained over long evolutionary times are so far not known. In this study, we addressed the functional relevance of eight accessory chromosomes in reference isolate IPO323. We induced chromosome losses by inhibiting the β-tubulin assembly during mitosis using carbendazim and generated several independent isogenic strains, each lacking one of the accessory chromosomes. We confirmed chromosome losses by electrophoretic karyotyping and whole-genome sequencing. To assess the importance of the individual chromosomes during host infection, we performed in planta assays comparing disease development results in wild-type and chromosome mutant strains. Loss of the accessory chromosomes 14, 16, 18, 19, and 21 resulted in increased virulence on wheat cultivar Runal but not on cultivars Obelisk, Titlis, and Riband. Moreover, some accessory chromosomes affected the switch from biotrophy to necrotrophy as strains lacking accessory chromosomes 14, 18, 19, and 21 showed a significantly earlier onset of necrosis than the wild type on the Runal cultivar. In general, we observed that the timing of the lifestyle switch affects the fitness of Z. tritici. Taking the results together, this study was the first to use a forward-genetics approach to demonstrate a cultivar-dependent functional relevance of the accessory chromosomes of Z. tritici during host infection.

Zymoseptoria tritici is a prominent fungal pathogen of wheat of worldwide distribution. This fungus shows a remarkable genome organization, with a large number of chromosomes that are present in only some isolates and therefore considered to be “accessory” chromosomes. To date, the function of these accessory chromosomes in Z. tritici has been unknown, although their maintenance in the species over evolutionary times suggests a functional relevance. Here we deleted whole accessory chromosomes to test the effect of these chromosomes on host specificity and virulence of the fungus. We show for the first time that some accessory chromosomes of Z. tritici affect the fitness of the fungus during host infection in a cultivar-dependent manner. These results show that the accessory chromosomes encode host-specific virulence determinants having a negative effect on fitness. Understanding the population dynamic of the accessory chromosomes and the molecular interaction of pathogen and plant traits is crucial to improve wheat-breeding strategies.