Expression profiling of the wheat pathogen Zymoseptoria tritici reveals genomic patterns of transcription and host-specific regulatory programs

Combined pangenomics and transcriptomics reveals core and redundant virulence processes in a rapidly evolving fungal plant pathogen

BACKGROUND

Studying genomic variation in rapidly evolving pathogens potentially enables identification of genes supporting their "core biology", being present, functional and expressed by all strains or "flexible biology", varying between strains. Genes supporting flexible biology may be considered to be "accessory", whilst the "core" gene set is likely to be important for common features of a pathogen species biology, including virulence on all host genotypes. The wheat-pathogenic fungus Zymoseptoria tritici represents one of the most rapidly evolving threats to global food security and was the focus of this study.

RESULTS

We constructed a pangenome of 18 European field isolates, with 12 also subjected to RNAseq transcription profiling during infection. Combining this data, we predicted a "core" gene set comprising 9807 sequences which were (1) present in all isolates, (2) lacking inactivating polymorphisms and (3) expressed by all isolates. A large accessory genome, consisting of 45% of the total genes, was also defined. We classified genetic and genomic polymorphism at both chromosomal and individual gene scales. Proteins required for essential functions including virulence had lower-than average sequence variability amongst core genes. Both core and accessory genomes encoded many small, secreted candidate effector proteins that likely interact with plant immunity. Viral vector-mediated transient in planta overexpression of 88 candidates failed to identify any which induced leaf necrosis characteristic of disease. However, functional complementation of a non-pathogenic deletion mutant lacking five core genes demonstrated that full virulence was restored by re-introduction of the single gene exhibiting least sequence polymorphism and highest expression.

CONCLUSIONS

These data support the combined use of pangenomics and transcriptomics for defining genes which represent core, and potentially exploitable, weaknesses in rapidly evolving pathogens.

Fungal plant pathogen "mutagenomics" reveals tagged and untagged mutations in Zymoseptoria tritici and identifies SSK2 as key morphogenesis and stress-responsive virulence factor

"Mutagenomics" is the combination of random mutagenesis, phenotypic screening, and whole-genome re-sequencing to uncover all tagged and untagged mutations linked with phenotypic changes in an organism. In this study, we performed a mutagenomics screen on the wheat pathogenic fungus Zymoseptoria tritici for altered morphogenetic switching and stress sensitivity phenotypes using Agrobacterium-mediated "random" T-DNA mutagenesis (ATMT). Biological screening identified four mutants which were strongly reduced in virulence on wheat. Whole genome re-sequencing defined the positions of the T-DNA insertion events and revealed several unlinked mutations potentially affecting gene functions. Remarkably, two independent reduced virulence mutant strains, with similarly altered stress sensitivities and aberrant hyphal growth phenotypes, were found to have a distinct loss of function mutations in the ZtSSK2 MAPKKK gene. One mutant strain had a direct T-DNA insertion affecting the predicted protein's N-terminus, while the other possessed an unlinked frameshift mutation towards the C-terminus. We used genetic complementation to restore both strains' wild-type (WT) function (virulence, morphogenesis, and stress response). We demonstrated that ZtSSK2 has a non-redundant function with ZtSTE11 in virulence through the biochemical activation of the stress-activated HOG1 MAPK pathway. Moreover, we present data suggesting that SSK2 has a unique role in activating this pathway in response to specific stresses. Finally, dual RNAseq-based transcriptome profiling of WT and SSK2 mutant strains revealed many HOG1-dependent transcriptional changes in the fungus during early infection and suggested that the host response does not discriminate between WT and mutant strains during this early phase. Together these data define new genes implicated in the virulence of the pathogen and emphasise the importance of a whole genome sequencing step in mutagenomic discovery pipelines.

Genome Organization and Copy-Number Variation Reveal Clues to Virulence Evolution in Coccidioides posadasii

The human fungal pathogen Coccidioides spp. causes valley fever, a treatment-refractory and sometimes deadly disease prevalent in arid regions of the western hemisphere. Fungal virulence in the mammalian host hinges on a switch between growth as hyphae and as large spherules containing infectious spores. How these virulence programs are encoded in the genome remains poorly understood. Drawing on Coccidioides genomic resources, we first discovered a new facet of genome organization in this system: spherule-gene islands, clusters of genes physically linked in the genome that exhibited specific mRNA induction in the spherule phase. Next, we surveyed copy-number variation genome-wide among strains of C. posadasii. Emerging from this catalog were spherule-gene islands with striking presence-absence differentiation between C. posadasii populations, a pattern expected from virulence factors subjected to different selective pressures across habitats. Finally, analyzing single-nucleotide differences across C. posadasii strains, we identified signatures of natural selection in spherule-expressed genes. Together, our data establish spherule-gene islands as candidate determinants of virulence and targets of selection in Coccidioides.

Heterochromatin in the fungal plant pathogen, Zymoseptoria tritici: Control of transposable elements, genome plasticity and virulence

Heterochromatin is a repressive chromatin state that plays key roles in the functional organisation of eukaryotic genomes. In fungal plant pathogens, effector genes that are required for host colonization tend to be associated with heterochromatic regions of the genome that are enriched with transposable elements. It has been proposed that the heterochromatin environment silences effector genes in the absence of host and dynamic chromatin remodelling facilitates their expression during infection. Here we discuss this model in the context of the key wheat pathogen, Zymoseptoria tritici. We cover progress in understanding the deposition and recognition of heterochromatic histone post translational modifications in Z. tritici and the role that heterochromatin plays in control of genome plasticity and virulence.

Strive or thrive: Trends in Phytophthora capsici gene expression in partially resistant pepper

Partial resistance in plants generally exerts a low selective pressure on pathogens, and thus ensuring their durability in agrosystems. However, little is known about the effect of partial resistance on the molecular mechanisms of pathogenicity, a knowledge that could advance plant breeding for sustainable plant health. Here we investigate the gene expression of Phytophthora capsici during infection of pepper (Capsicum annuum L.), where only partial genetic resistance is reported, using Illumina RNA-seq. Comparison of transcriptomes of P. capsici infecting susceptible and partially resistant peppers identified a small number of genes that redirected its own resources into lipid biosynthesis to subsist on partially resistant plants. The adapted and non-adapted isolates of P. capsici differed in expression of genes involved in nucleic acid synthesis and transporters. Transient ectopic expression of the RxLR effector genes CUST_2407 and CUST_16519 in pepper lines differing in resistance levels revealed specific host-isolate interactions that either triggered local necrotic lesions (hypersensitive response or HR) or elicited leave abscission (extreme resistance or ER), preventing the spread of the pathogen to healthy tissue. Although these effectors did not unequivocally explain the quantitative host resistance, our findings highlight the importance of plant genes limiting nutrient resources to select pepper cultivars with sustainable resistance to P. capsici.

Induction of distinct plant cell death programs by secreted proteins from the wheat pathogen Zymoseptoria tritici

Cell death processes in eukaryotes shape normal development and responses to the environment. For plant-microbe interactions, initiation of host cell death plays an important role in determining disease outcomes. Cell death pathways are frequently initiated following detection of pathogen-derived molecules which can lead to resistance or susceptibility to disease depending on pathogen lifestyle. We previously identified several small secreted proteins (SSPs) from the wheat-infecting fungus Zymoseptoria tritici that induce rapid cell death in Nicotiana benthamiana following Agrobacterium-mediated delivery and expression (agroinfiltration). Here we investigated whether the execution of host cells was mechanistically similar in response to different Z. tritici SSPs. Using RNA sequencing, we found that transient expression of four Z. tritici SSPs led to massive transcriptional reprogramming within 48 h of agroinfiltration. We observed that distinct host gene expression profiles were induced dependent on whether cell death occurs in a cell surface immune receptor-dependent or -independent manner. These gene expression profiles involved differential transcriptional networks mediated by WRKY, NAC and MYB transcription factors. In addition, differential expression of genes belonging to different classes of receptor-like proteins and receptor-like kinases was observed. These data suggest that different Z. tritici SSPs trigger differential transcriptional reprogramming in plant cells.

Transcriptional response to host chemical cues underpins the expansion of host range in a fungal plant pathogen lineage

The host range of parasites is an important factor in assessing the dynamics of disease epidemics. The evolution of pathogens to accommodate new hosts may lead to host range expansion, a process the molecular bases of which are largely enigmatic. The fungus Sclerotinia sclerotiorum has been reported to parasitize more than 400 plant species from diverse eudicot families while its close relative, S. trifoliorum, is restricted to plants from the Fabaceae family. We analyzed S. sclerotiorum global transcriptome reprogramming on hosts from six botanical families and reveal a flexible, host-specific transcriptional program. We generated a chromosome-level genome assembly for S. trifoliorum and found near-complete gene space conservation in two representative strains of broad and narrow host range Sclerotinia species. However, S. trifoliorum showed increased sensitivity to the Brassicaceae defense compound camalexin. Comparative analyses revealed a lack of transcriptional response to camalexin in the S. trifoliorum strain and suggest that regulatory variation in detoxification and effector genes at the population level may associate with the genetic accommodation of Brassicaceae in the Sclerotinia host range. Our work proposes transcriptional plasticity and the co-existence of signatures for generalist and polyspecialist adaptive strategies in the genome of a plant pathogen.

Two Novel Dimorphism-Related Virulence Factors of Zymoseptoria tritici Identified Using Agrobacterium-Mediated Insertional Mutagenesis

Diseases caused by dimorphic phytopathogenic and systemic dimorphic fungi have markedly increased in prevalence in the last decades, and understanding the morphogenic transition to the virulent state might yield novel means of controlling dimorphic fungi. The dimorphic fungus Z. tritici causes significant economic impact on wheat production, and yet the regulation of the dimorphic switch, a key first step in successful plant colonization, is still largely unexplored in this fungus. The fungus is amenable to suppression by fungicides at this switch point, and the identification of the factors controlling the dimorphic switch provides a potential source of novel targets to control Septoria tritici blotch (STB). Inhibition of the dimorphic switch can potentially prevent penetration and avoid any damage to the host plant. The aim of the current work was to unveil genetic determinants of the dimorphic transition in Z. tritici by using a forward genetics strategy. Using this approach, we unveiled two novel factors involved in the switch to the pathogenic state and used reverse genetics and complementation to confirm the role of the novel virulence factors and further gained insight into the role of these genes, using transcriptome analysis via RNA-Seq. The transcriptomes generated potentially contain key determinants of the dimorphic transition.

Fatty acid dioxygenase-cytochrome P450 fusion enzymes of filamentous fungal pathogens

Oxylipins designate oxygenated unsaturated C₁₈ fatty acids. Many filamentous fungi pathogens contain dioxygenases (DOX) in oxylipin biosynthesis with homology to human cyclooxygenases. They contain a DOX domain, which is often fused to a functional cytochrome P450 at the C-terminal end. A Tyr radical in the DOX domain initiates dioxygenation of linoleic acid by hydrogen abstraction with formation of 8-, 9-, or 10-hydroperoxy metabolites. The P450 domains can catalyze heterolytic cleavage of 8- and 10-hydroperoxides with oxidation of the heme thiolate iron for hydroxylation at C-5, C-7, C-9, or C-11 and for epoxidation of the 12Z double bond; thus displaying linoleate diol synthase (LDS) and epoxy alcohol synthase (EAS) activities. LSD activities are present in the rice blast pathogen Magnaporthe oryzae, Botrytis cinerea causing grey mold and the black scurf pathogen Rhizoctonia solani. 10R-DOX-EAS has been found in M. oryzae and Fusarium oxysporum. The P450 domains may also catalyze homolytic cleavage of 8- and 9-hydroperoxy fatty acids and dehydration to produce epoxides with an adjacent double bond, i.e., allene oxides, thus displaying 8- and 9-DOX-allene oxide synthases (AOS). F. oxysporum, F. graminearum, and R. solani express 9S-DOX-AOS and Zymoseptoria tritici 8S-and 9R-DOX-AOS. Homologues are present in endemic human-pathogenic fungi with extensive studies in Aspergillus fumigatus, A. flavus (also a plant pathogen) as well as the genetic model A. nidulans. 8R-and 10R-DOX appear to bind fatty acids "headfirst" in the active site, whereas 9S-DOX binds them "tail first" in analogy with cyclooxygenases. The biological relevance of 8R-DOX-5,8-LDS (also designated PpoA) was first discovered in relation to sporulation of A. nidulans and recently for development and programmed hyphal branching of A. fumigatus. Gene deletion DOX-AOS homologues in F. verticillioides, A. flavus, and A. nidulans alters, inter alia, mycotoxin production, sporulation, and gene expression.

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.

Dynamics of transposable elements in recently diverged fungal pathogens: lineage-specific transposable element content and efficiency of genome defenses

Transposable elements (TEs) impact genome plasticity, architecture, and evolution in fungal plant pathogens. The wide range of TE content observed in fungal genomes reflects diverse efficacy of host-genome defense mechanisms that can counter-balance TE expansion and spread. Closely related species can harbor drastically different TE repertoires. The evolution of fungal effectors, which are crucial determinants of pathogenicity, has been linked to the activity of TEs in pathogen genomes. Here, we describe how TEs have shaped genome evolution of the fungal wheat pathogen Zymoseptoria tritici and four closely related species. We compared de novo TE annotations and repeat-induced point mutation signatures in 26 genomes from the Zymoseptoria species-complex. Then, we assessed the relative insertion ages of TEs using a comparative genomics approach. Finally, we explored the impact of TE insertions on genome architecture and plasticity. The 26 genomes of Zymoseptoria species reflect different TE dynamics with a majority of recent insertions. TEs associate with accessory genome compartments, with chromosomal rearrangements, with gene presence/absence variation, and with effectors in all Zymoseptoria species. We find that the extent of RIP-like signatures varies among Z. tritici genomes compared to genomes of the sister species. The detection of a reduction of RIP-like signatures and TE recent insertions in Z. tritici reflects ongoing but still moderate TE mobility.

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.

WITHDRAWN: Fatty acid dioxygenase-cytochrome P450 fusion enzymes of the top 10 fungal pathogens in molecular plant pathology and human-pathogenic fungi

This article has been withdrawn at the request of the author(s) and/or editor. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at https://www.elsevier.com/about/our-business/policies/article-withdrawal

Evolution and Plasticity of the Transcriptome Under Temperature Fluctuations in the Fungal Plant Pathogen Zymoseptoria tritici

Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.

Genome compartmentalization predates species divergence in the plant pathogen genus Zymoseptoria

Background

Antagonistic co-evolution can drive rapid adaptation in pathogens and shape genome architecture. Comparative genome analyses of several fungal pathogens revealed highly variable genomes, for many species characterized by specific repeat-rich genome compartments with exceptionally high sequence variability. Dynamic genome structure may enable fast adaptation to host genetics. The wheat pathogen Zymoseptoria tritici with its highly variable genome, has emerged as a model organism to study genome evolution of plant pathogens. Here, we compared genomes of Z. tritici isolates and of sister species infecting wild grasses to address the evolution of genome composition and structure.

Results

Using long-read technology, we sequenced and assembled genomes of Z. ardabiliae, Z. brevis, Z. pseudotritici and Z. passerinii, together with two isolates of Z. tritici. We report a high extent of genome collinearity among Zymoseptoria species and high conservation of genomic, transcriptomic and epigenomic signatures of compartmentalization. We identify high gene content variability both within and between species. In addition, such variability is mainly limited to the accessory chromosomes and accessory compartments. Despite strong host specificity and non-overlapping host-range between species, predicted effectors are mainly shared among Zymoseptoria species, yet exhibiting a high level of presence-absence polymorphism within Z. tritici. Using in planta transcriptomic data from Z. tritici, we suggest different roles for the shared orthologs and for the accessory genes during infection of their hosts.

Conclusion

Despite previous reports of high genomic plasticity in Z. tritici, we describe here a high level of conservation in genomic, epigenomic and transcriptomic composition and structure across the genus Zymoseptoria. The compartmentalized genome allows the maintenance of a functional core genome co-occurring with a highly variable accessory genome.

Host-specialized transcriptome of plant-associated organisms

Living organisms respond to their immediate environment by modulating their genetic programme to perform adapted functions. Eukaryotic organisms that associate with plants (fungi, oomycetes, insects, …) alter their transcriptome in a host-specific manner. Recent comparative transcriptomic studies revealed that host-specialized transcriptomes consist of a limited set of genes. Such a set typically encodes proteins that modulate host structures and functions (predicted effectors and other secreted proteins), control nutrient assimilation (proteases, transporters), and maintain cellular homeostasis (oxidoreductases, detoxification enzymes). We conclude by discussing open mechanistic and evolutionary questions and integrated approaches to move beyond descriptive studies.

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.

Hiding in plain sight: New virus genomes discovered via a systematic analysis of fungal public transcriptomes

The distribution and diversity of RNA viruses in fungi is incompletely understood due to the often cryptic nature of mycoviral infections and the focused study of primarily pathogenic and/or economically important fungi. As most viruses that are known to infect fungi possess either single-stranded or double-stranded RNA genomes, transcriptomic data provides the opportunity to query for viruses in diverse fungal samples without any a priori knowledge of virus infection. Here we describe a systematic survey of all transcriptomic datasets from fungi belonging to the subphylum Pezizomycotina. Using a simple but effective computational pipeline that uses reads discarded during normal RNA-seq analyses, followed by identification of a viral RNA-dependent RNA polymerase (RdRP) motif in de novo assembled contigs, 59 viruses from 44 different fungi were identified. Among the viruses identified, 88% were determined to be new species and 68% are, to our knowledge, the first virus described from the fungal species. Comprehensive analyses of both nucleotide and inferred protein sequences characterize the phylogenetic relationships between these viruses and the known set of mycoviral sequences and support the classification of up to four new families and two new genera. Thus the results provide a deeper understanding of the scope of mycoviral diversity while also increasing the distribution of fungal hosts. Further, this study demonstrates the suitability of analyzing RNA-seq data to facilitate rapid discovery of new viruses.

The genomic determinants of adaptive evolution in a fungal pathogen

Unravelling the strength, frequency, and distribution of selective variants along the genome as well as the underlying factors shaping this distribution are fundamental goals of evolutionary biology. Antagonistic host-pathogen coevolution is thought to be a major driver of genome evolution between interacting species. While rapid evolution of pathogens has been documented in several model organisms, the genetic mechanisms of their adaptation are still poorly understood and debated, particularly the role of sexual reproduction. Here, we apply a population genomic approach to infer genome-wide patterns of selection among 13 isolates of Zymoseptoria tritici, a fungal pathogen characterized by extremely high genetic diversity, gene density, and recombination rates. We report that the genome of Z. tritici undergoes a high rate of adaptive substitutions, with 44% of nonsynonymous substitutions being adaptive on average. This fraction reaches 68% in so-called effector genes encoding determinants of pathogenicity, and the distribution of fitness effects differs in this class of genes as they undergo adaptive mutations with stronger positive fitness effects, but also more slightly deleterious mutations. Besides the globally high rate of adaptive substitutions, we report a negative relationship between pN/pS and the fine-scale recombination rate and a strong positive correlation between the rate of adaptive nonsynonymous substitutions (ωa) and recombination rate. This result suggests a pervasive role of both background selection and Hill-Robertson interference even in a species with an exceptionally high recombination rate (60 cM/Mb on average). While transposable elements (TEs) have been suggested to contribute to adaptation by creating compartments of fast-evolving genomic regions, we do not find a significant effect of TEs on the rate of adaptive mutations. Overall our study suggests that sexual recombination is a significant driver of genome evolution, even in rapidly evolving organisms subject to recurrent mutations with large positive effects.

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.

Dissection of Ramularia Leaf Spot Disease by Integrated Analysis of Barley and Ramularia collo-cygni Transcriptome Responses

Ramularia leaf spot disease (RLS), caused by the ascomycete fungus Ramularia collo-cygni, has emerged as a major economic disease of barley. No substantial resistance has been identified, so far, among barley genotypes and, based on the epidemiology of the disease, a quantitative genetic determinacy of RLS has been suggested. The relative contributions of barley and R. collo-cygni genetics to disease infection and epidemiology are practically unknown. Here, we present an integrated genome-wide analysis of host and pathogen transcriptome landscapes identified in a sensitive barley cultivar following infection by an aggressive R. collo-cygni isolate. We compared transcriptional responses in the infected and noninfected leaf samples in order to identify which molecular events are associated with RLS symptom development. We found a large proportion of R. collo-cygni genes to be expressed in planta and that many were also closely associated with the infection stage. The transition from surface to apoplastic colonization was associated with downregulation of cell wall-degrading genes and upregulation of nutrient uptake and resistance to oxidative stresses. Interestingly, the production of secondary metabolites was dynamically regulated within the fungus, indicating that R. collo-cygni produces a diverse panel of toxic compounds according to the infection stage. A defense response against R. collo-cygni was identified in barley at the early, asymptomatic infection and colonization stages. We found activation of ethylene signaling, jasmonic acid signaling, and phenylpropanoid and flavonoid pathways to be highly induced, indicative of a classical response to necrotrophic pathogens. Disease development was found to be associated with gene expression patterns similar to those found at the onset of leaf senescence, when nutrients, possibly, are used by the infecting fungus. These analyses, combining both barley and R. collo-cygni transcript profiles, demonstrate the activation of complex transcriptional programs in both organisms.

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.

Extraordinary Genome Instability and Widespread Chromosome Rearrangements During Vegetative Growth

The haploid genome of the pathogenic fungus Zymoseptoria tritici is contained on "core" and "accessory" chromosomes. While 13 core chromosomes are found in all strains, as many as eight accessory chromosomes show presence/absence variation and rearrangements among field isolates. The factors influencing these presence/absence polymorphisms are so far unknown. We investigated chromosome stability using experimental evolution, karyotyping, and genome sequencing. We report extremely high and variable rates of accessory chromosome loss during mitotic propagation in vitro and in planta Spontaneous chromosome loss was observed in 2 to >50% of cells during 4 weeks of incubation. Similar rates of chromosome loss in the closely related Zymoseptoria ardabiliae suggest that this extreme chromosome dynamic is a conserved phenomenon in the genus. Elevating the incubation temperature greatly increases instability of accessory and even core chromosomes, causing severe rearrangements involving telomere fusion and chromosome breakage. Chromosome losses do not affect the fitness of Zymoseptoria tritici in vitro, but some lead to increased virulence, suggesting an adaptive role of this extraordinary chromosome instability.

Pan-Parastagonospora Comparative Genome Analysis-Effector Prediction and Genome Evolution

We report a fungal pan-genome study involving Parastagonospora spp., including 21 isolates of the wheat (Triticum aestivum) pathogen Parastagonospora nodorum, 10 of the grass-infecting Parastagonospora avenae, and 2 of a closely related undefined sister species. We observed substantial variation in the distribution of polymorphisms across the pan-genome, including repeat-induced point mutations, diversifying selection and gene gains and losses. We also discovered chromosome-scale inter and intraspecific presence/absence variation of some sequences, suggesting the occurrence of one or more accessory chromosomes or regions that may play a role in host-pathogen interactions. The presence of known pathogenicity effector loci SnToxA, SnTox1, and SnTox3 varied substantially among isolates. Three P. nodorum isolates lacked functional versions for all three loci, whereas three P. avenae isolates carried one or both of the SnTox1 and SnTox3 genes, indicating previously unrecognized potential for discovering additional effectors in the P. nodorum-wheat pathosystem. We utilized the pan-genomic comparative analysis to improve the prediction of pathogenicity effector candidates, recovering the three confirmed effectors among our top-ranked candidates. We propose applying this pan-genomic approach to identify the effector repertoire involved in other host-microbe interactions involving necrotrophic pathogens in the Pezizomycotina.

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.

Accessories Make the Outfit: Accessory Chromosomes and Other Dispensable DNA Regions in Plant-Pathogenic Fungi

Fungal pathogen genomes can often be divided into core and accessory regions. Accessory regions ARs) may be comprised of either ARs (within core chromosomes (CCs) or wholly dispensable (accessory) chromosomes (ACs). Fungal ACs and ARs typically accumulate mutations and structural rearrangements more rapidly over time than CCs and many harbor genes relevant to host-pathogen interactions. These regions are of particular interest in plant pathology and include host-specific virulence factors and secondary metabolite synthesis gene clusters. This review outlines known ACs and ARs in fungal genomes, methods used for their detection, their common properties that differentiate them from the core genome, and what is currently known of their various roles in pathogenicity. Reports on the evolutionary processes generating and shaping AC and AR compartments are discussed, including repeat induced point mutation and breakage fusion bridge cycles. Previously ACs have been studied extensively within key genera, including Fusarium, Zymoseptoria, and Alternaria, but are growing in frequency of observation and perceived importance across a wider range of fungal species. Recent advances in sequencing technologies permit affordable genome assembly and resequencing of populations that will facilitate further discovery and routine screening of ACs.

Fine-Scale Recombination Maps of Fungal Plant Pathogens Reveal Dynamic Recombination Landscapes and Intragenic Hotspots

Meiotic recombination is an important driver of evolution. Variability in the intensity of recombination across chromosomes can affect sequence composition, nucleotide variation, and rates of adaptation. In many organisms, recombination events are concentrated within short segments termed recombination hotspots. The variation in recombination rate and positions of recombination hotspot can be studied using population genomics data and statistical methods. In this study, we conducted population genomics analyses to address the evolution of recombination in two closely related fungal plant pathogens: the prominent wheat pathogen Zymoseptoria tritici and a sister species infecting wild grasses Z. ardabiliae We specifically addressed whether recombination landscapes, including hotspot positions, are conserved in the two recently diverged species and if recombination contributes to rapid evolution of pathogenicity traits. We conducted a detailed simulation analysis to assess the performance of methods of recombination rate estimation based on patterns of linkage disequilibrium, in particular in the context of high nucleotide diversity. Our analyses reveal overall high recombination rates, a lack of suppressed recombination in centromeres, and significantly lower recombination rates on chromosomes that are known to be accessory. The comparison of the recombination landscapes of the two species reveals a strong correlation of recombination rate at the megabase scale, but little correlation at smaller scales. The recombination landscapes in both pathogen species are dominated by frequent recombination hotspots across the genome including coding regions, suggesting a strong impact of recombination on gene evolution. A significant but small fraction of these hotspots colocalize between the two species, suggesting that hotspot dynamics contribute to the overall pattern of fast evolving recombination in these species.

Stress and sexual reproduction affect the dynamics of the wheat pathogen effector AvrStb6 and strobilurin resistance

Host resistance and fungicide treatments are cornerstones of plant-disease control. Here, we show that these treatments allow sex and modulate parenthood in the fungal wheat pathogen Zymoseptoria tritici. We demonstrate that the Z. tritici-wheat interaction complies with the gene-for-gene model by identifying the effector AvrStb6, which is recognized by the wheat resistance protein Stb6. Recognition triggers host resistance, thus implying removal of avirulent strains from pathogen populations. However, Z. tritici crosses on wheat show that sex occurs even with an avirulent parent, and avirulence alleles are thereby retained in subsequent populations. Crossing fungicide-sensitive and fungicide-resistant isolates under fungicide pressure results in a rapid increase in resistance-allele frequency. Isolates under selection always act as male donors, and thus disease control modulates parenthood. Modeling these observations for agricultural and natural environments reveals extended durability of host resistance and rapid emergence of fungicide resistance. Therefore, fungal sex has major implications for disease control.

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.

Non-parent of Origin Expression of Numerous Effector Genes Indicates a Role of Gene Regulation in Host Adaption of the Hybrid Triticale Powdery Mildew Pathogen

Powdery mildew is an important disease of cereals. It is caused by one species, Blumeria graminis, which is divided into formae speciales each of which is highly specialized to one host. Recently, a new form capable of growing on triticale (B.g. triticale) has emerged through hybridization between wheat and rye mildews (B.g. tritici and B.g. secalis, respectively). In this work, we used RNA sequencing to study the molecular basis of host adaptation in B.g. triticale. We analyzed gene expression in three B.g. tritici isolates, two B.g. secalis isolates and two B.g. triticale isolates and identified a core set of putative effector genes that are highly expressed in all formae speciales. We also found that the genes differentially expressed between isolates of the same form as well as between different formae speciales were enriched in putative effectors. Their coding genes belong to several families including some which contain known members of mildew avirulence (Avr) and suppressor (Svr) genes. Based on these findings we propose that effectors play an important role in host adaptation that is mechanistically based on Avr-Resistance gene-Svr interactions. We also found that gene expression in the B.g. triticale hybrid is mostly conserved with the parent-of-origin, but some genes inherited from B.g. tritici showed a B.g. secalis-like expression. Finally, we identified 11 unambiguous cases of putative effector genes with hybrid-specific, non-parent of origin gene expression, and we propose that they are possible determinants of host specialization in triticale mildew. These data suggest that altered expression of multiple effector genes, in particular Avr and Svr related factors, might play a role in mildew host adaptation based on hybridization.

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.

A conserved fungal glycosyltransferase facilitates pathogenesis of plants by enabling hyphal growth on solid surfaces

Pathogenic fungi must extend filamentous hyphae across solid surfaces to cause diseases of plants. However, the full inventory of genes which support this is incomplete and many may be currently concealed due to their essentiality for the hyphal growth form. During a random T-DNA mutagenesis screen performed on the pleomorphic wheat (Triticum aestivum) pathogen Zymoseptoria tritici, we acquired a mutant unable to extend hyphae specifically when on solid surfaces. In contrast "yeast-like" growth, and all other growth forms, were unaffected. The inability to extend surface hyphae resulted in a complete loss of virulence on plants. The affected gene encoded a predicted type 2 glycosyltransferase (ZtGT2). Analysis of >800 genomes from taxonomically diverse fungi highlighted a generally widespread, but discontinuous, distribution of ZtGT2 orthologues, and a complete absence of any similar proteins in non-filamentous ascomycete yeasts. Deletion mutants of the ZtGT2 orthologue in the taxonomically un-related fungus Fusarium graminearum were also severely impaired in hyphal growth and non-pathogenic on wheat ears. ZtGT2 expression increased during filamentous growth and electron microscopy on deletion mutants (ΔZtGT2) suggested the protein functions to maintain the outermost surface of the fungal cell wall. Despite this, adhesion to leaf surfaces was unaffected in ΔZtGT2 mutants and global RNAseq-based gene expression profiling highlighted that surface-sensing and protein secretion was also largely unaffected. However, ΔZtGT2 mutants constitutively overexpressed several transmembrane and secreted proteins, including an important LysM-domain chitin-binding virulence effector, Zt3LysM. ZtGT2 likely functions in the synthesis of a currently unknown, potentially minor but widespread, extracellular or outer cell wall polysaccharide which plays a key role in facilitating many interactions between plants and fungi by enabling hyphal growth on solid matrices.

Adaptation to the Host Environment by Plant-Pathogenic Fungi

Many fungi can live both saprophytically and as endophyte or pathogen inside a living plant. In both environments, complex organic polymers are used as sources of nutrients. Propagation inside a living host also requires the ability to respond to immune responses of the host. We review current knowledge of how plant-pathogenic fungi do this. First, we look at how fungi change their global gene expression upon recognition of the host environment, leading to secretion of effectors, enzymes, and secondary metabolites; changes in metabolism; and defense against toxic compounds. Second, we look at what is known about the various cues that enable fungi to sense the presence of living plant cells. Finally, we review literature on transcription factors that participate in gene expression in planta or are suspected to be involved in that process because they are required for the ability to cause disease.

Karyotype Variability in Plant-Pathogenic Fungi

Recent advances in genetic and molecular technologies gradually paved the way for the transition from traditional fungal karyotyping to more comprehensive chromosome biology studies. Extensive chromosomal polymorphisms largely resulting from chromosomal rearrangements (CRs) are widely documented in fungal genomes. These extraordinary CRs in fungi generate substantial genome plasticity compared to other eukaryotic organisms. Here, we review the most recent findings on fungal CRs and their underlying mechanisms and discuss the functional consequences of CRs for adaptation, fungal evolution, host range, and pathogenicity of fungal plant pathogens in the context of chromosome biology. In addition to a complement of permanent chromosomes called core chromosomes, the genomes of many fungal pathogens comprise distinct unstable chromosomes called dispensable chromosomes (DCs) that also contribute to chromosome polymorphisms. Compared to the core chromosomes, the structural features of DCs usually differ for gene density, GC content, housekeeping genes, and recombination frequency. Despite their dispensability for normal growth and development, DCs have important biological roles with respect to pathogenicity in some fungi but not in others. Therefore, their evolutionary origin is also reviewed in relation to overall fungal physiology and pathogenicity.

A Matter of Scale and Dimensions: Chromatin of Chromosome Landmarks in the Fungi

Chromatin and chromosomes of fungi are highly diverse and dynamic, even within species. Much of what we know about histone modification enzymes, RNA interference, DNA methylation, and cell cycle control was first addressed in Saccharomyces cerevisiae, Schizosaccharomyces pombe, Aspergillus nidulans, and Neurospora crassa. Here, we examine the three landmark regions that are required for maintenance of stable chromosomes and their faithful inheritance, namely, origins of DNA replication, telomeres and centromeres. We summarize the state of recent chromatin research that explains what is required for normal function of these specialized chromosomal regions in different fungi, with an emphasis on the silencing mechanism associated with subtelomeric regions, initiated by sirtuin histone deacetylases and histone H3 lysine 27 (H3K27) methyltransferases. We explore mechanisms for the appearance of "accessory" or "conditionally dispensable" chromosomes and contrast what has been learned from studies on genome-wide chromosome conformation capture in S. cerevisiae, S. pombe, N. crassa, and Trichoderma reesei. While most of the current knowledge is based on work in a handful of genetically and biochemically tractable model organisms, we suggest where major knowledge gaps remain to be closed. Fungi will continue to serve as facile organisms to uncover the basic processes of life because they make excellent model organisms for genetics, biochemistry, cell biology, and evolutionary biology.

Transcriptome analysis of smut fungi reveals widespread intergenic transcription and conserved antisense transcript expression

Background

Biotrophic fungal plant pathogens cause billions of dollars in losses to North American crops annually. The model for functional investigation of these fungi is Ustilago maydis. Its 20.5 Mb annotated genome sequence has been an excellent resource for investigating biotrophic plant pathogenesis. Expressed-sequence tag libraries and microarray hybridizations have provided insight regarding the type of transcripts produced by U. maydis but these analyses were not comprehensive and there were insufficient data for transcriptome comparison to other smut fungi. To improve transcriptome annotation and enable comparative analyses, comprehensive strand-specific RNA-seq was performed on cell-types of three related smut species: U. maydis (common smut of corn), Ustilago hordei (covered smut of barley), and Sporisorium reilianum (head smut of corn).

Results

In total, >1 billion paired-end sequence reads were obtained from haploid cell, dikaryon and teliospore RNA of U. maydis, haploid cell RNA of U. hordei, and haploid and dikaryon cell RNA of S. reilianum. The sequences were assembled into transfrags using Trinity, and updated gene models were created using PASA and categorized with Cufflinks Cuffcompare. Representative genes that were predicted for the first time with these RNA-seq analyses and genes with novel annotation features were independently assessed by reverse transcriptase PCR. The analyses indicate hundreds more predicted proteins, relative to the previous genome annotation, could be produced by U. maydis from altered transcript forms, and that the number of non-coding RNAs produced, including transcribed intergenic sequences and natural antisense transcripts, approximately equals the number of mRNAs. This high representation of non-coding RNAs appears to be a conserved feature of the smut fungi regardless of whether they have RNA interference machinery. Approximately 50% of the identified NATs were conserved among the smut fungi.

Conclusions

Overall, these analyses revealed: 1) smut genomes encode a number of transcriptional units that is twice the number of annotated protein-coding genes, 2) a small number of intergenic transcripts may encode proteins with characteristics of fungal effectors, 3) the vast majority of intergenic and antisense transcripts do not contain ORFs, 4) a large proportion of the identified antisense transcripts were detected at orthologous loci among the smut fungi, and 5) there is an enrichment of functional categories among orthologous loci that suggests antisense RNAs could have a genome-wide, non-RNAi-mediated, influence on gene expression in smut fungi.

Electronic supplementary material

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

The novel fungal-specific gene FUG1 has a role in pathogenicity and fumonisin biosynthesis in Fusarium verticillioides

Fusarium verticillioides is a globally important pathogen of maize, capable of causing severe yield reductions and economic losses. In addition, F. verticillioides produces toxic secondary metabolites during kernel colonization that pose significant threats to human and animal health. Fusarium verticillioides and other plant-pathogenic fungi possess a large number of genes with no known or predicted function, some of which could encode novel virulence factors or antifungal targets. In this study, we identified and characterized the novel gene FUG1 (Fungal Unknown Gene 1) in F. verticillioides through functional genetics. Deletion of FUG1 impaired maize kernel colonization and fumonisin biosynthesis. In addition, deletion of FUG1 increased sensitivity to the antimicrobial compound 2-benzoxazolinone and to hydrogen peroxide, which indicates that FUG1 may play a role in mitigating stresses associated with host defence. Transcriptional profiling via RNA-sequencing (RNA-seq) identified numerous fungal genes that were differentially expressed in the kernel environment following the deletion of FUG1, including genes involved in secondary metabolism and mycelial development. Sequence analysis of the Fug1 protein provided evidence for nuclear localization, DNA binding and a domain of unknown function associated with previously characterized transcriptional regulators. This information, combined with the observed transcriptional reprogramming in the deletion mutant, suggests that FUG1 represents a novel class of fungal transcription factors or genes otherwise involved in signal transduction.

A small secreted protein in Zymoseptoria tritici is responsible for avirulence on wheat cultivars carrying the Stb6 resistance gene

Zymoseptoria tritici is the causal agent of Septoria tritici blotch, a major pathogen of wheat globally and the most damaging pathogen of wheat in Europe. A gene-for-gene (GFG) interaction between Z. tritici and wheat cultivars carrying the Stb6 resistance gene has been postulated for many years, but the genes have not been identified. We identified AvrStb6 by combining quantitative trait locus mapping in a cross between two Swiss strains with a genome-wide association study using a natural population of c. 100 strains from France. We functionally validated AvrStb6 using ectopic transformations. AvrStb6 encodes a small, cysteine-rich, secreted protein that produces an avirulence phenotype on wheat cultivars carrying the Stb6 resistance gene. We found 16 nonsynonymous single nucleotide polymorphisms among the tested strains, indicating that AvrStb6 is evolving very rapidly. AvrStb6 is located in a highly polymorphic subtelomeric region and is surrounded by transposable elements, which may facilitate its rapid evolution to overcome Stb6 resistance. AvrStb6 is the first avirulence gene to be functionally validated in Z. tritici, contributing to our understanding of avirulence in apoplastic pathogens and the mechanisms underlying GFG interactions between Z. tritici and wheat.

Comparative Transcriptome Analyses in Zymoseptoria tritici Reveal Significant Differences in Gene Expression Among Strains During Plant Infection

Zymoseptoria tritici is an ascomycete fungus that causes Septoria tritici blotch, a globally distributed foliar disease on wheat. Z. tritici populations are highly polymorphic and exhibit significant quantitative variation for virulence. Despite its importance, the genes responsible for quantitative virulence in this pathogen remain largely unknown. We investigated the expression profiles of four Z. tritici strains differing in virulence in an experiment conducted under uniform environmental conditions. Transcriptomes were compared at four different infection stages to characterize the regulation of gene families thought to be involved in virulence and to identify new virulence factors. The major components of the fungal infection transcriptome showed consistent expression profiles across strains. However, strain-specific regulation was observed for many genes, including some encoding putative virulence factors. We postulate that strain-specific regulation of virulence factors can determine the outcome of Z. tritici infections. We show that differences in gene expression may be major determinants of virulence variation among Z. tritici strains, adding to the already known contributions to virulence variation based on differences in gene sequence and gene presence/absence polymorphisms.

Codon optimization underpins generalist parasitism in fungi

The range of hosts that parasites can infect is a key determinant of the emergence and spread of disease. Yet, the impact of host range variation on the evolution of parasite genomes remains unknown. Here, we show that codon optimization underlies genome adaptation in broad host range parasites. We found that the longer proteins encoded by broad host range fungi likely increase natural selection on codon optimization in these species. Accordingly, codon optimization correlates with host range across the fungal kingdom. At the species level, biased patterns of synonymous substitutions underpin increased codon optimization in a generalist but not a specialist fungal pathogen. Virulence genes were consistently enriched in highly codon-optimized genes of generalist but not specialist species. We conclude that codon optimization is related to the capacity of parasites to colonize multiple hosts. Our results link genome evolution and translational regulation to the long-term persistence of generalist parasitism.

New fungal defensin-like peptides provide evidence for fold change of proteins in evolution

Defensins containing a consensus cystine framework, Cys_[1]…Cys_[2]X₃Cys_[3]…Cys_[4]… Cys_[5]X₁Cys_[6] (X, any amino acid except Cys; …, variable residue numbers), are extensively distributed in a variety of multicellular organisms (plants, fungi and invertebrates) and essentially involved in immunity as microbicidal agents. This framework is a prerequisite for forming the cysteine-stabilized α-helix and β-sheet (CSαβ) fold, in which the two invariant motifs, Cys_[2]X₃Cys_[3]/Cys_[5]X₁Cys_[6], are key determinants of fold formation. By using a computational genomics approach, we identified a large superfamily of fungal defensin-like peptides (fDLPs) in the phytopathogenic fungal genus – Zymoseptoria, which includes 132 structurally typical and 63 atypical members. These atypical fDLPs exhibit an altered cystine framework and accompanying fold change associated with their secondary structure elements and disulfide bridge patterns, as identified by protein structure modelling. Despite this, they definitely are homologous with the typical fDLPs in view of their precise gene structure conservation and identical precursor organization. Sequence and structural analyses combined with functional data suggest that most of Zymoseptoria fDLPs might have lost their antimicrobial activity. The present study provides a clear example of fold change in the evolution of proteins and is valuable in establishing remote homology among peptide superfamily members with different folds.

The HAP Complex Governs Fumonisin Biosynthesis and Maize Kernel Pathogenesis in Fusarium verticillioides

Contamination of maize ( Zea mays ) with fumonisins produced by the fungus Fusarium verticillioides is a global concern for food safety. Fumonisins are a group of polyketide-derived secondary metabolites linked to esophageal cancer and neural tube birth defects in humans and numerous toxicoses in livestock. Despite the importance of fumonisins in global maize production, the regulation of fumonisin biosynthesis during kernel pathogenesis is poorly understood. The HAP complex is a conserved, heterotrimeric transcriptional regulator that binds the consensus sequence CCAAT to modulate gene expression. Recently, functional characterization of the Hap3 subunit linked the HAP complex to the regulation of secondary metabolism and stalk rot pathogenesis in F. verticillioides . Here, we determine the involvement of HAP3 in fumonisin biosynthesis and kernel pathogenesis. Deletion of HAP3 suppressed fumonisin biosynthesis on both nonviable and live maize kernels and impaired pathogenesis in living kernels. Transcriptional profiling via RNA sequencing indicated that the HAP complex regulates at least 1,223 genes in F. verticillioides , representing nearly 10% of all predicted genes. Disruption of the HAP complex caused the misregulation of biosynthetic gene clusters underlying the production of secondary metabolites, including fusarins. Taken together, these results reveal that the HAP complex is a central regulator of fumonisin biosynthesis and kernel pathogenesis and works as both a positive and negative regulator of secondary metabolism in F. verticillioides .

Comparative transcriptomic analyses of Zymoseptoria tritici strains show complex lifestyle transitions and intraspecific variability in transcription profiles

Zymoseptoria tritici causes Septoria tritici blotch (STB) on wheat. The disease interaction is characterized by clearly defined temporal phases of infection, ultimately resulting in the death of host tissue. Zymoseptoria tritici is a highly polymorphic species with significant intraspecific variation in virulence profiles. We generated a deep transcriptomic sequencing dataset spanning the entire time course of an infection using a previously uncharacterized, highly virulent Z. tritici strain isolated from a Swiss wheat field. We found that seven clusters of gene transcription profiles explained the progression of the infection. The earliest highly up-regulated genes included chloroperoxidases, which may help the fungus cope with plant defences. The onset of necrotrophy was characterized by a concerted up-regulation of proteases, plant cell wall-degrading enzymes and lipases. Functions related to nutrition and growth characterized late necrotrophy and the transition to saprotrophic growth on dead plant tissue. We found that the peak up-regulation of genes essential for mating coincided with the necrotrophic phase. We performed an intraspecies comparative transcriptomics analysis using a comparable time course infection experiment of the genome reference isolate IPO323. Major components of the fungal infection transcriptome were conserved between the two strains. However, individual small, secreted proteins, proteases and cell wall-degrading enzymes showed strongly differentiated transcriptional profiles between isolates. Our analyses illustrate that successful STB infections involve complex transcriptomic remodelling to up-regulate distinct gene functions. Heterogeneity in transcriptomes among isolates may explain some of the considerable variation in virulence and host specialization found within the species.

Utilizing Gene Tree Variation to Identify Candidate Effector Genes in Zymoseptoria tritici

Zymoseptoria tritici is a host-specific, necrotrophic pathogen of wheat. Infection by Z. tritici is characterized by its extended latent period, which typically lasts 2 wks, and is followed by extensive host cell death, and rapid proliferation of fungal biomass. This work characterizes the level of genomic variation in 13 isolates, for which we have measured virulence on 11 wheat cultivars with differential resistance genes. Between the reference isolate, IPO323, and the 13 Australian isolates we identified over 800,000 single nucleotide polymorphisms, of which ∼10% had an effect on the coding regions of the genome. Furthermore, we identified over 1700 probable presence/absence polymorphisms in genes across the Australian isolates using de novo assembly. Finally, we developed a gene tree sorting method that quickly identifies groups of isolates within a single gene alignment whose sequence haplotypes correspond with virulence scores on a single wheat cultivar. Using this method, we have identified < 100 candidate effector genes whose gene sequence correlates with virulence toward a wheat cultivar carrying a major resistance gene.

Next-generation re-sequencing as a tool for rapid bioinformatic screening of presence and absence of genes and accessory chromosomes across isolates of Zymoseptoria tritici

The wheat pathogen Zymoseptoria tritici possesses a large number of accessory chromosomes that may be present or absent in its genome. The genome of the reference isolate IPO323 has been assembled to a very high standard and contains 21 full length chromosome sequences, 8 of which represent accessory chromosomes. The IPO323 reference, when combined with low-cost next-generation sequencing and bioinformatics, can be used as a powerful tool to assess the presence or absence of accessory chromosomes. We present an outline of a range of bioinformatics techniques that can be applied to the analysis of presence-absence variation among accessory chromosomes across 13 novel isolates of Z. tritici.

Cell biology of Zymoseptoria tritici: Pathogen cell organization and wheat infection

•

Cell biology of the infection begins to shed light on the host–pathogen interaction.

•

The cell biology of the fungal pathogen is highly understudied.

•

Intensified cell biology research promises new fungicide targets and will help mode-of-action studies.

Cell biological research in the wheat pathogen Zymoseptoria tritici (formerly Mycosphaerella graminicola) has led to a good understanding of the histology of the infection process. Expression profiling and bioinformatic approaches, combined with molecular studies on signaling pathways, effectors and potential necrosis factors provides first insight into the complex interplay between the host and the pathogen. Cell biological studies will help to further our understanding of the infection strategy of the fungus. The cellular organization and intracellular dynamics of the fungus itself is largely unexplored. Insight into essential cellular processes within the pathogen will expand our knowledge of the basic biology of Z. tritici, thereby providing putative new anti-fungal targets.