Coevolution and life cycle specialization of plant cell wall degrading enzymes in a hemibiotrophic pathogen

The Molecular Dialogue Between Zymoseptoria tritici and Wheat

Zymoseptoria tritici is a highly damaging pathogen that causes high wheat yield losses in temperate climates. Z. tritici emerged during the domestication of wheat in the Fertile Crescent and has been extensively used as a model system for population genetic and genomic studies. New genetic tools and resources have provided a better understanding of the molecular components involved in the wheat-Z. tritici interaction, which is highlighted by the cloning of three wheat resistance genes and four Z. tritici avirulence genes. Despite the considerable progress made in the last few years, the mechanisms that mediate Z. tritici colonization remain largely unknown. In this review, we summarize the latest advances in understanding the molecular components mediating wheat-Z. tritici interactions, and we discuss future research lines to close current knowledge gaps. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

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.

Endogenous cell wall degrading enzyme LytD is important for the biocontrol activity of Bacillus subtilis

Autolysins are endogenous cell wall degrading enzymes (CWDEs) in bacteria that remodel the peptidoglycan layer of its own cell wall. In the Bacillus subtilis genome, at least 35 autolysin genes have been identified. However, the study of their roles in bacterial physiology has been hampered by their complexity and functional redundancy. B. subtilis GLB191 is an effective biocontrol strain against grape downy mildew disease, the biocontrol effect of which results from both direct effect against the pathogen and stimulation of the plant defense. In this study, we show that the autolysin N-acetylglucosaminidase LytD, a major autolysin of vegetative growth in B. subtilis, plays an important role in its biocontrol activity against grape downy mildew. Disruption of lytD resulted in reduced suppression of the pathogen Plasmopara viticola and stimulation of the plant defense. LytD is also shown to affect the biofilm formation and colonization of B. subtilis on grape leaves. This is the first report that demonstrates the role of an endogenous CWDE in suppressing plant disease infection of a biological control microorganism. These findings not only expand our knowledge on the biological function of autolysins but also provide a new target to promote the biocontrol activity of B. subtilis.

Sas3-mediated histone acetylation regulates effector gene activation in a fungal plant pathogen

IMPORTANCE

Pathogen infections require the production of effectors that enable host colonization. Effectors have diverse functions and are only expressed at certain stages of the infection cycle. Thus, effector genes are tightly regulated by several mechanisms, including chromatin remodeling. Here, we investigate the role of histone acetylation in effector gene activation in the fungal wheat pathogen Zymoseptoria tritici. We demonstrate that lysine acetyltransferases (KATs) are essential for the spatiotemporal regulation of effector genes. We show that the KAT Sas3 is involved in leaf symptom development and pycnidia formation. Importantly, our results indicate that Sas3 controls histone acetylation of effector loci and is a regulator of effector gene activation during stomatal penetration. Overall, our work demonstrates the key role of histone acetylation in regulating gene expression associated with plant infection.

A facultative ectomycorrhizal association is triggered by organic nitrogen

Increasing nitrogen (N) deposition often tends to negatively impact the functions of belowground ectomycorrhizal networks, although the exact molecular mechanisms underlying this trait are still unclear. Here, we assess how the root-associated fungus Clitopilus hobsonii establishes an ectomycorrhiza-like association with its host tree Populus tomentosa and how this interaction is favored by organic N over mineral N. The establishment of a functional symbiosis in the presence of organic N promotes plant growth and the transfer of ¹⁵N from the fungus to above ground plant tissues. Genomic traits and in planta transcriptional signatures suggest that C. hobsonii may have a dual lifestyle with saprotrophic and mutualistic traits. For example, several genes involved in the digestion of cellulose and hemicellulose are highly expressed during the interaction, whereas the expression of multiple copies of pectin-digesting genes is tightly controlled. Conversely, the nutritional mutualism is dampened in the presence of ammonium (NH₄⁺) or nitrate (NO₃^-). Increasing levels of NH₄⁺ led to a higher expression of pectin-digesting genes and a continuous increase in hydrogen peroxide production in roots, whereas the presence of NO₃^- resulted in toxin production. In summary, our results suggest that C. hobsonii is a facultative ectomycorrhizal fungus. Access to various forms of N acts as an on/off switch for mutualism caused by large-scale fungal physiological remodeling. Furthermore, the abundance of pectin-degrading enzymes with distinct expression patterns during functional divergence after exposure to NH₄⁺ or organic N is likely to be central to the transition from parasitism to mutualism.

Automated Quantitative Measurement of Yellow Halos Suggests Activity of Necrotrophic Effectors in Septoria tritici Blotch

Many necrotrophic plant pathogens utilize host-selective toxins or necrotrophic effectors during the infection process. We hypothesized that the chlorotic yellow halos frequently observed around necrotic lesions caused by the wheat pathogen Zymoseptoria tritici could result from the activity of necrotrophic effectors interacting with the products of toxin sensitivity genes. As an initial step toward testing this hypothesis, we developed an automated image analysis (AIA) workflow that could quantify the degree of yellow halo formation occurring in wheat leaves naturally infected by a highly diverse pathogen population under field conditions. This AIA based on statistical learning was applied to more than 10,000 naturally infected leaves collected from 335 wheat cultivars grown in a replicated field experiment. We estimated a high heritability (h² = 0.71) for the degree of yellow halo formation, suggesting that this quantitative trait has a significant genetic component. Using genome-wide association mapping, we identified six chromosome segments significantly associated with the yellow halo phenotype. Most of these segments contained candidate genes associated with targets of necrotrophic effectors in other necrotrophic pathogens. Our findings conform with the hypothesis that toxin sensitivity genes could account for a significant fraction of the observed variation in quantitative resistance to Septoria tritici blotch. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

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.

Genome of the destructive oomycete Phytophthora cinnamomi provides insights into its pathogenicity and adaptive potential

BACKGROUND

Phytophthora cinnamomi is an oomycete pathogen of global relevance. It is considered as one of the most invasive species, which has caused irreversible damage to natural ecosystems and horticultural crops. There is currently a lack of a high-quality reference genome for this species despite several attempts that have been made towards sequencing its genome. The lack of a good quality genome sequence has been a setback for various genetic and genomic research to be done on this species. As a consequence, little is known regarding its genome characteristics and how these contribute to its pathogenicity and invasiveness.

RESULTS

In this work we generated a high-quality genome sequence and annotation for P. cinnamomi using a combination of Oxford Nanopore and Illumina sequencing technologies. The annotation was done using RNA-Seq data as supporting gene evidence. The final assembly consisted of 133 scaffolds, with an estimated genome size of 109.7 Mb, N50 of 1.18 Mb, and BUSCO completeness score of 97.5%. Genome partitioning analysis revealed that P. cinnamomi has a two-speed genome characteristic, similar to that of other oomycetes and fungal plant pathogens. In planta gene expression analysis revealed up-regulation of pathogenicity-related genes, suggesting their important roles during infection and host degradation.

CONCLUSION

This study has provided a high-quality reference genome and annotation for P. cinnamomi. This is among the best assembled genomes for any Phytophthora species assembled to date and thus resulted in improved identification and characterization of pathogenicity-related genes, some of which were undetected in previous versions of genome assemblies. Phytophthora cinnamomi harbours a large number of effector genes which are located in the gene-poor regions of the genome. This unique genomic partitioning provides P. cinnamomi with a high level of adaptability and could contribute to its success as a highly invasive species. Finally, the genome sequence, its annotation and the pathogenicity effectors identified in this study will serve as an important resource that will enable future studies to better understand and mitigate the impact of this important pathogen.

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.

Nitric oxide signalling in plant interactions with pathogenic fungi and oomycetes

Nitric oxide (NO) and reactive nitrogen species have emerged as crucial signalling and regulatory molecules across all organisms. In plants, fungi, and fungi-like oomycetes, NO is involved in the regulation of multiple processes during their growth, development, reproduction, responses to the external environment, and biotic interactions. It has become evident that NO is produced and used as a signalling and defence cue by both partners in multiple forms of plant interactions with their microbial counterparts, ranging from symbiotic to pathogenic modes. This review summarizes current knowledge on the role of NO in plant-pathogen interactions, focused on biotrophic, necrotrophic, and hemibiotrophic fungi and oomycetes. Actual advances and gaps in the identification of NO sources and fate in plant and pathogen cells are discussed. We review the decisive role of time- and site-specific NO production in germination, oriented growth, and active penetration by filamentous pathogens of the host tissues, as well in pathogen recognition, and defence activation in plants. Distinct functions of NO in diverse interactions of host plants with fungal and oomycete pathogens of different lifestyles are highlighted, where NO in interplay with reactive oxygen species governs successful plant colonization, cell death, and establishment of resistance.

Recognition and defence of plant-infecting fungal pathogens

Attempted infections of plants with fungi result in diverse outcomes ranging from symptom-less resistance to severe disease and even death of infected plants. The deleterious effect on crop yield have led to intense focus on the cellular and molecular mechanisms that explain the difference between resistance and susceptibility. This research has uncovered plant resistance or susceptibility genes that explain either dominant or recessive inheritance of plant resistance with many of them coding for receptors that recognize pathogen invasion. Approaches based on cell biology and phytochemistry have contributed to identifying factors that halt an invading fungal pathogen from further invasion into or between plant cells. Plant chemical defence compounds, antifungal proteins and structural reinforcement of cell walls appear to slow down fungal growth or even prevent fungal penetration in resistant plants. Additionally, the hypersensitive response, in which a few cells undergo a strong local immune reaction, including programmed cell death at the site of infection, stops in particular biotrophic fungi from spreading into surrounding tissue. In this review, we give a general overview of plant recognition and defence of fungal parasites tracing back to the early 20th century with a special focus on Triticeae and on the progress that was made in the last 30 years.

Conditional promoters to investigate gene function during wheat infection by Zymoseptoria tritici

The fungus Zymoseptoria tritici causes Septoria tritici leaf blotch, which poses a serious threat to temperate-grown wheat. Recently, we described a raft of molecular tools to study the biology of this fungus in vitro. Amongst these are 5 conditional promoters (Pnar1, Pex1A, Picl1, Pgal7, PlaraB), which allow controlled over-expression or repression of target genes in cells grown in liquid culture. However, their use in the host-pathogen interaction in planta was not tested. Here, we investigate the behaviour of these promoters by quantitative live cell imaging of green-fluorescent protein-expressing cells during 6 stages of the plant infection process. We show that Pnar1 and Picl1 are repressed in planta and demonstrate their suitability for studying essential gene expression and function in plant colonisation. The promoters Pgal7 and Pex1A are not fully-repressed in planta, but are induced during pycnidiation. This indicates the presence of inducing galactose or xylose and/or arabinose, released from the plant cell wall by the activity of fungal hydrolases. In contrast, the PlaraB promoter, which normally controls expression of an α-l-arabinofuranosidase B, is strongly induced inside the leaf. This suggests that the fungus is exposed to L-arabinose in the mesophyll apoplast. Taken together, this study establishes 2 repressible promoters (Pnar1 and Picl1) and three inducible promoters (Pgal7, Pex1A, PlaraB) for molecular studies in planta. Moreover, we provide circumstantial evidence for plant cell wall degradation during the biotrophic phase of Z. tritici infection.

Asynchronous development of Zymoseptoria tritici infection in wheat

•

Zymoseptoria tritici passes 6 morphologically defined stages during infection.

•

Surface-located spores and hyphae are found for up to 17/18 days.

•

Entry through stomata occurs from 1 to 13 days post infection.

•

Mesophyll apoplast colonisation continuously increases during infection.

•

Up to 5 stages co-exist in infected leaves at a given time.

The fungus Zymoseptoria tritici causes Septoria tritici blotch of wheat. Pathogenicity begins with spore germination, followed by stomata invasion by hyphae, mesophyll colonization and fruiting body formation. It was previously found that entry into the plant via stomata occurs in a non-synchronized way over several days, while later developmental steps, such as early and late fruiting body formation, were reported to follow each other in time. This suggests synchronization of the pathogen population in planta prior to sporulation. Here, we image a fluorescent Z. tritici IPO323-derived strain during infection. We describe 6 morphologically distinct developmental stages, and determine their abundance in infected leaves, with time post inoculation. This demonstrates that 3-5 stages co-exist in infected tissues at any given time. Thus, later stages of pathogen development also occur asynchronously amongst the population of infecting cells. This merits consideration when interpreting transcriptomics or proteomics data gathered from infected plants.

Seed-Derived Microbial Colonization of Wild Emmer and Domesticated Bread Wheat (Triticum dicoccoides and T. aestivum) Seedlings Shows Pronounced Differences in Overall Diversity and Composition

The composition of the plant microbiota may be altered by ecological and evolutionary changes in the host population. Seed-associated microbiota, expected to be largely vertically transferred, have the potential to coadapt with their host over generations. Strong directional selection and changes in the genetic composition of plants during domestication and cultivation may have impacted the assembly and transmission of seed-associated microbiota. Nonetheless, the effect of plant speciation and domestication on the composition of these microbes is poorly understood. Here, we have investigated the composition of bacteria and fungi associated with the wild emmer wheat (Triticum dicoccoides) and domesticated bread wheat (Triticum aestivum). We show that vertically transmitted bacteria, but not fungi, of domesticated bread wheat species T. aestivum are less diverse and more inconsistent among individual plants compared to those of the wild emmer wheat species T. dicoccoides. We propagated wheat seeds under sterile conditions to characterize the colonization of seedlings by seed-associated microbes. Hereby, we show markedly different community compositions and diversities of leaf and root colonizers of the domesticated bread wheat compared to the wild emmer wheat. By propagating the wild emmer wheat and domesticated bread wheat in two different soils, we furthermore reveal a small effect of plant genotype on microbiota assembly. Our results suggest that domestication and prolonged breeding have impacted the vertically transferred bacteria, but only to a lesser extent have affected the soil-derived microbiota of bread wheat.IMPORTANCE Genetic and physiological changes associated with plant domestication have been studied for many crop species. Still little is known about the impact of domestication on the plant-associated microbiota. In this study, we analyze the seed-associated and soil-derived bacterial and fungal microbiota of domesticated bread wheat and wild emmer wheat. We show a significant difference in the seed-associated, but not soil-derived, bacterial communities of the wheat species. Interestingly, we find less pronounced effects on the fungal communities. Overall, this study provides novel insight into the diversity of vertically transmitted microbiota of wheat and thereby contributes to our understanding of wheat as a "metaorganism." Insight into the wheat microbiota is of fundamental importance for the development of improved crops.

Presence of ice-nucleating Pseudomonas on wheat leaves promotes Septoria tritici blotch disease (Zymoseptoria tritici) via a mutually beneficial interaction

Zymoseptoria tritici causes Septoria tritici blotch (STB) of wheat, an economically important disease causing yield losses of up to 10% despite the use of fungicides and resistant cultivars. Z. tritici infection is symptomless for around 10 days, during which time the fungus grows randomly across the leaf surface prior to entry through stomata. Wounded leaves show faster, more extensive STB, suggesting that wounds facilitate fungal entry. Wheat leaves also host epiphytic bacteria; these include ice-nucleating (INA+) bacteria, which induce frost damage at warmer temperatures than it otherwise occurs. Here, STB is shown to be more rapid and severe when wheat is exposed to both INA+ bacteria and sub-zero temperatures. This suggests that ice-nucleation-induced wounding of the wheat leaf provides additional openings for fungal entry. INA+ bacterial populations are shown to benefit from the presence of Z. tritici, indicating that this microbial interaction is mutualistic. Finally, control of INA+ bacteria is shown to reduce STB.

Chromatin Dynamics Contribute to the Spatiotemporal Expression Pattern of Virulence Genes in a Fungal Plant Pathogen

Dynamic changes in transcription profiles are key for the success of pathogens in colonizing their hosts. In many pathogens, genes associated with virulence, such as effector genes, are located in regions of the genome that are rich in transposable elements and heterochromatin. The contribution of chromatin modifications to gene expression in pathogens remains largely unknown. Using a combination of a reporter gene-based approach and chromatin immunoprecipitation, we show that the heterochromatic environment of effector genes in the fungal plant pathogen Zymoseptoria tritici is a key regulator of their specific spatiotemporal expression patterns. Enrichment in trimethylated lysine 27 of histone H3 dictates the repression of effector genes in the absence of the host. Chromatin decondensation during host colonization, featuring a reduction in this repressive modification, indicates a major role for epigenetics in effector gene induction. Our results illustrate that chromatin modifications triggered during host colonization determine the specific expression profile of effector genes at the cellular level and, hence, provide new insights into the regulation of virulence in fungal plant pathogens.IMPORTANCE Fungal plant pathogens possess a large repertoire of genes encoding putative effectors, which are crucial for infection. Many of these genes are expressed at low levels in the absence of the host but are strongly induced at specific stages of the infection. The mechanisms underlying this transcriptional reprogramming remain largely unknown. We investigated the role of the genomic environment and associated chromatin modifications of effector genes in controlling their expression pattern in the fungal wheat pathogen Zymoseptoria tritici Depending on their genomic location, effector genes are epigenetically repressed in the absence of the host and during the initial stages of infection. Derepression of effector genes occurs mainly during and after penetration of plant leaves and is associated with changes in histone modifications. Our work demonstrates the role of chromatin in shaping the expression of virulence components and, thereby, the interaction between fungal pathogens and their plant hosts.

The conifer root rot pathogens Heterobasidion irregulare and Heterobasidion occidentale employ different strategies to infect Norway spruce

Heterobasidion irregulare and H. occidentale are two closely related conifer root rot pathogens in the H. annosum sensu lato (s.l.) species complex. The two species H. irregulare and H. occidentale have different host preference with pine and non-pine tree species favored, respectively. The comparison of transcriptomes of H. irregulare and H. occidentale growing in Norway spruce bark, a susceptible host non-native to North America, showed large differences in gene expression. Heterobasidion irregulare induced more genes involved in detoxification of host compounds and in production of secondary metabolites, while the transcriptome induced in H. occidentale was more oriented towards carbohydrate degradation. Along with their separated evolutionary history, the difference might be driven by their host preferences as indicated by the differentially expressed genes enriched in particular Gene Ontology terms.

Diversity, structure, and synteny of the cutinase gene of Colletotrichum species

Colletotrichum species complexes are among the top 10 economically important fungal plant pathogens worldwide because they can infect climacteric and nonclimacteric fruit at the pre and/or postharvest stages. C. truncatum is the major pathogen responsible for anthracnose of green and red bell pepper fruit worldwide. C. brevisporum was recently reported to be a minor pathogen of red bell pepper fruit in Trinidad, but has recently been reported as pathogenic to other host species in other countries. The ability of these phytopathogens to produce and secrete cutinase is required for dismantling the cuticle of the host plant and, therefore, crucial to the necrotrophic phase of their infection strategy. In vitro bioassays using different lipid substrates confirmed the ability of C. truncatum and C. brevisporum isolates from green and red bell peppers to secrete cutinase. The diversity, structure and organization and synteny of the cutinase gene were determined among different Colletotrichum species. Cluster analysis indicated a low level of nucleotide variation among C. truncatum sequences. Nucleotide sequences of C. brevisporum were more related to C. truncatum cutinase nucleotide sequences than to C. gloeosporioides. Cluster patterns coincided with haplotype and there was evidence of significant positive selection with no recombination signatures. The structure of the cutinase gene included two exons with one intervening intron and, therefore, one splice variant. Although amino acid sequences were highly conserved among C. truncatum isolates, diversity "hot spots" were revealed when the 66-amino acid coding region of 200 fungal species was compared. Twenty cutinase orthologues were detected among different fungal species, whose common ancestor is Pezizomycotina and it is purported that these orthologues arose through a single gene duplication event prior to speciation. The cutinase domain was retained both in structure and arrangement among 34 different Colletotrichum species. The order of aligned genomic blocks between species and the arrangement of flanking protein domains were also conserved and shared for those domains immediately located at the N- and C-terminus of the cutinase domain. Among these were an RNA recognition motif, translation elongation factor, signal peptide, pentatricopeptide repeat, and Hsp70 family of chaperone proteins, all of which support the expression of the cutinase gene. The findings of this study are important to understanding the evolution of the cutinase gene in C. truncatum as a key component of the biotrophic-necrotrophic switch which may be useful in developing gene-targeting strategies to decrease the pathogenic potential of Colletotrichum species.

Does the Latent Period of Leaf Fungal Pathogens Reflect Their Trophic Type? A Meta-Analysis of Biotrophs, Hemibiotrophs, and Necrotrophs

We performed a meta-analysis to search for a relation between the trophic type and latent period of fungal pathogens. The pathogen incubation period and the level of resistance of the hosts were also investigated. This ecological knowledge would help us to more efficiently regulate crop epidemics for different types of pathogens. We gathered latent period data from 103 studies dealing with 51 fungal pathogens of the three major trophic types (25 biotrophs, 15 hemibiotrophs, and 11 necrotrophs), representing 2,542 mean latent periods. We show that these three trophic types display significantly different latent periods. Necrotrophs exhibited the shortest latent periods (<100 degree-days [DD]), biotrophs had intermediate ones (between 100 and 200 DD), and hemibiotrophs had the longest latent periods (>200 DD). We argue that this relation between trophic type and latent period points to two opposing host exploitation strategies: necrotrophs mount a rapid destructive attack on the host tissue, whereas biotrophs and hemibiotrophs avoid or delay the damaging phase. We query the definition of hemibiotrophic pathogens and discuss whether the length of the latent period is determined by the physiological limits inherent to each trophic type or by the adaptation of pathogens of different trophic types to the contrasting conditions experienced in their interaction with the host.

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.

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.

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.

Distinct Trajectories of Massive Recent Gene Gains and Losses in Populations of a Microbial Eukaryotic Pathogen

Differences in gene content are a significant source of variability within species and have an impact on phenotypic traits. However, little is known about the mechanisms responsible for the most recent gene gains and losses. We screened the genomes of 123 worldwide isolates of the major pathogen of wheat Zymoseptoria tritici for robust evidence of gene copy number variation. Based on orthology relationships in three closely related fungi, we identified 599 gene gains and 1,024 gene losses that have not yet reached fixation within the focal species. Our analyses of gene gains and losses segregating in populations showed that gene copy number variation arose preferentially in subtelomeres and in proximity to transposable elements. Recently lost genes were enriched in virulence factors and secondary metabolite gene clusters. In contrast, recently gained genes encoded mostly secreted protein lacking a conserved domain. We analyzed the frequency spectrum at loci segregating a gene presence–absence polymorphism in four worldwide populations. Recent gene losses showed a significant excess in low-frequency variants compared with genome-wide single nucleotide polymorphism, which is indicative of strong negative selection against gene losses. Recent gene gains were either under weak negative selection or neutral. We found evidence for strong divergent selection among populations at individual loci segregating a gene presence–absence polymorphism. Hence, gene gains and losses likely contributed to local adaptation. Our study shows that microbial eukaryotes harbor extensive copy number variation within populations and that functional differences among recently gained and lost genes led to distinct evolutionary trajectories.

Inheritance of phenotypic traits in the progeny of a Ceratocystis interspecific cross

Ceratocystis fimbriata is a fungal plant pathogen that causes black rot on Ipomoea batatas. Based on inoculation studies on numerous tree species, the pathogen is known to be host specific. The closely related species, Ceratocystis manginecans, causes severe wilt on a broad range of tree hosts, including Mangifera indica, Acacia mangium and other leguminous tree species. The genetic factors underlying the pathogenicity and host specificity of Ceratocystis species have rarely been investigated. In this study, an F₁ population of 70 recombinant progeny from a cross between C. fimbriata and C. manginecans was generated and the inheritance of various phenotypic traits was investigated. Results showed that colony colour, growth rate, asexual spore production and aggressiveness to I. batatas and A. mangium are all quantitative traits with high levels of heritability. However, conidia production and aggressiveness appeared to be regulated by a small number of genes. No correlation could be found between aggressiveness and other phenotypic traits, suggesting that these are inherited independently. This is the first study to consider genetic inheritance of pathogenicity and host specificity in Ceratocystis species and the results will contribute, in future, to the identification of quantitative trait loci and candidate genes associated with the traits investigated.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Evolutionary analyses of the avirulence effector AvrStb6 in global populations of Zymoseptoria tritici identify candidate amino acids involved in recognition

We analysed the population genetic diversity of AvrStb6, the first avirulence gene cloned from the wheat pathogen Zymoseptoria tritici, using 142 Z. tritici strains sampled from four wheat fields growing on three continents. Although AvrStb6 was located in a recombination hotspot, it was found in every strain, with 71 polymorphic sites that produced 41 distinct DNA haplotypes encoding 30 AvrStb6 protein isoforms. An AvrStb6 homologue was found in the closest known relative, Z. pseudotritici, but not in three other closely related Zymoseptoria species, indicating that this gene has emerged in Zymoseptoria quite recently. Two AvrStb6 homologues with nucleotide similarities greater than 70% were identified on chromosome 10 in all Z. tritici isolates, suggesting that AvrStb6 belongs to a multigene family of candidate effectors that has expanded recently through gene duplication. The AvrStb6 sequences exhibited strong evidence for non-neutral evolution, including a large number of non-synonymous mutations, with significant positive diversifying selection operating on nine of the 82 codons. It appears that balancing selection is operating across the entire gene in natural field populations. There was also evidence for co-evolving codons within the gene that may reflect compensatory mutations associated with the evasion of recognition by Stb6. Intragenic recombination also appears to have affected the diversity of AvrStb6.

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.

Evolutionary Analysis of Pectin Lyases of the Genus Colletotrichum

Pectin lyases (PNLs) are important enzymes that are involved in plant cell wall degradation during the infection process. Colletotrichum is a diverse genus of fungi, which allows the study of the evolution of PNLs and their possible role in pathogen-host interactions and lifestyle adaptations. The phylogenetic reconstruction of PNLs from Colletotrichum and analysis of selection pressures showed the formation of protein lineages by groups of species with different selection pressures and specific patterns. The analysis of positive selection at individual sites using different methods allowed for the identification of three codons with evidence of positive selection in the oligosaccharide-binding region and two codons on the antiparallel sheet, which may influence the interaction with the substrate. Seven codons on the surface of the protein, mainly in the peripheral helices of the PNLs, could have an important function in evasion of plant defenses, as has been proposed in other enzymes. According to our results, it is possible that events of genetic duplication occurred in ancestral lines, followed by episodes of genetic diversification and gene loss, probably influenced by differences in the composition of the host cell wall. Additionally, different patterns of evolution in Colletotrichum appear to be molded by a strong purifying selection and positive selection episodes that forged the observed evolutionary patterns, possibly influenced by host interaction or substrate specificity. This work represents a starting point for the study of sites that may be important for evasion of plant defenses and biotechnological purposes.

A role for random, humidity-dependent epiphytic growth prior to invasion of wheat by Zymoseptoria tritici

Zymoseptoria tritici causes Septoria leaf blotch of wheat. The prevailing paradigm of the Z. tritici-wheat interaction assumes fungal ingress through stomata within 24-48h, followed by days of symptomless infection. This is extrapolated from studies testing the mode of fungal ingress under optimal infection conditions. Here, we explicitly assess the timing of entry, using GFP-tagged Z. tritici. We show that early entry is comparatively rare, and extended epiphytic growth possible. We test the hypotheses that our data diverge from earlier studies due to: i. random ingress of Z. tritici into the leaf, with some early entry events; ii. previous reliance upon fungal stains, combined with poor attachment of Z. tritici to the leaf, leading to increased likelihood of observing internal versus external growth, compared to using GFP; iii. use of exceptionally high humidity to promote entry in previous studies. We combine computer simulation of leaf-surface growth with thousands of in planta observations to demonstrate that while spores germinate rapidly on the leaf, over 95% of fungi remain epiphytic, growing randomly over the leaf for ten days or more. We show that epiphytic fungi are easily detached from leaves by rinsing and that humidity promotes epiphytic growth, increasing infection rates. Together, these results explain why epiphytic growth has been dismissed and early ingress assumed. The prolonged epiphytic phase should inform studies of pathogenicity and virulence mutants, disease control strategies, and interpretation of the observed low in planta growth, metabolic quiescence and evasion of plant defences by Zymoseptoria during symptomless infection.

A fungal wheat pathogen evolved host specialization by extensive chromosomal rearrangements

Fungal pathogens can rapidly evolve virulence towards resistant crops in agricultural ecosystems. Gains in virulence are often mediated by the mutation or deletion of a gene encoding a protein recognized by the plant immune system. However, the loci and the mechanisms of genome evolution enabling rapid virulence evolution are poorly understood. We performed genome-wide association mapping on a global collection of 106 strains of Zymoseptoria tritici, the most damaging pathogen of wheat in Europe, to identify polymorphisms linked to virulence on two wheat varieties. We found 25 distinct genomic loci associated with reproductive success of the pathogen. However, no locus was shared between the host genotypes, suggesting host specialization. The main locus associated with virulence encoded a highly expressed, small secreted protein. Population genomic analyses showed that the gain in virulence was explained by a segregating gene deletion polymorphism. The deletion was likely adaptive by preventing detection of the encoded protein. Comparative genomics of closely related species showed that the locus emerged de novo since speciation. A large cluster of transposable elements in direct proximity to the locus generated extensive rearrangements leading to multiple independent gene losses. Our study demonstrates that rapid turnover in the chromosomal structure of a pathogen can drive host specialization.

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.

Comparative genomics to explore phylogenetic relationship, cryptic sexual potential and host specificity of Rhynchosporium species on grasses

BACKGROUND

The Rhynchosporium species complex consists of hemibiotrophic fungal pathogens specialized to different sweet grass species including the cereal crops barley and rye. A sexual stage has not been described, but several lines of evidence suggest the occurrence of sexual reproduction. Therefore, a comparative genomics approach was carried out to disclose the evolutionary relationship of the species and to identify genes demonstrating the potential for a sexual cycle. Furthermore, due to the evolutionary very young age of the five species currently known, this genus appears to be well-suited to address the question at the molecular level of how pathogenic fungi adapt to their hosts.

RESULTS

The genomes of the different Rhynchosporium species were sequenced, assembled and annotated using ab initio gene predictors trained on several fungal genomes as well as on Rhynchosporium expressed sequence tags. Structures of the rDNA regions and genome-wide single nucleotide polymorphisms provided a hypothesis for intra-genus evolution. Homology screening detected core meiotic genes along with most genes crucial for sexual recombination in ascomycete fungi. In addition, a large number of cell wall-degrading enzymes that is characteristic for hemibiotrophic and necrotrophic fungi infecting monocotyledonous hosts were found. Furthermore, the Rhynchosporium genomes carry a repertoire of genes coding for polyketide synthases and non-ribosomal peptide synthetases. Several of these genes are missing from the genome of the closest sequenced relative, the poplar pathogen Marssonina brunnea, and are possibly involved in adaptation to the grass hosts. Most importantly, six species-specific genes coding for protein effectors were identified in R. commune. Their deletion yielded mutants that grew more vigorously in planta than the wild type.

CONCLUSION

Both cryptic sexuality and secondary metabolites may have contributed to host adaptation. Most importantly, however, the growth-retarding activity of the species-specific effectors suggests that host adaptation of R. commune aims at extending the biotrophic stage at the expense of the necrotrophic stage of pathogenesis. Like other apoplastic fungi Rhynchosporium colonizes the intercellular matrix of host leaves relatively slowly without causing symptoms, reminiscent of the development of endophytic fungi. Rhynchosporium may therefore become an object for studying the mutualism-parasitism transition.

Population Genomics of Fungal and Oomycete Pathogens

We are entering a new era in plant pathology in which whole-genome sequences of many individuals of a pathogen species are becoming readily available. Population genomics aims to discover genetic mechanisms underlying phenotypes associated with adaptive traits such as pathogenicity, virulence, fungicide resistance, and host specialization, as genome sequences or large numbers of single nucleotide polymorphisms become readily available from multiple individuals of the same species. This emerging field encompasses detailed genetic analyses of natural populations, comparative genomic analyses of closely related species, identification of genes under selection, and linkage analyses involving association studies in natural populations or segregating populations resulting from crosses. The era of pathogen population genomics will provide new opportunities and challenges, requiring new computational and analytical tools. This review focuses on conceptual and methodological issues as well as the approaches to answering questions in population genomics. The major steps start with defining relevant biological and evolutionary questions, followed by sampling, genotyping, and phenotyping, and ending in analytical methods and interpretations. We provide examples of recent applications of population genomics to fungal and oomycete plant pathogens.

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.

Life cycle specialization of filamentous pathogens - colonization and reproduction in plant tissues

Filamentous plant pathogens explore host tissues to obtain nutrients for growth and reproduction. Diverse strategies for tissue invasion, defense manipulation, and colonization of inter and intra-cellular spaces have evolved. Most research has focused on effector molecules, which are secreted to manipulate plant immunity and facilitate infection. Effector genes are often found to evolve rapidly in response to the antagonistic host-pathogen co-evolution but other traits are also subject to adaptive evolution during specialization to the anatomy, biochemistry and ecology of different plant hosts. Although not directly related to virulence, these traits are important components of specialization but little is known about them. We present and discuss specific life cycle traits that facilitate exploration of plant tissues and underline the importance of increasing our insight into the biology of plant pathogens.

QTL mapping of temperature sensitivity reveals candidate genes for thermal adaptation and growth morphology in the plant pathogenic fungus Zymoseptoria tritici

Different thermal environments impose strong, differential selection on populations, leading to local adaptation, but the genetic basis of thermal adaptation is poorly understood. We used quantitative trait locus (QTL) mapping in the fungal wheat pathogen Zymoseptoria tritici to study the genetic architecture of thermal adaptation and identify candidate genes. Four wild-type strains originating from the same thermal environment were crossed to generate two mapping populations with 263 (cross 1) and 261 (cross 2) progeny. Restriction site-associated DNA sequencing was used to genotype 9745 (cross 1) and 7333 (cross 2) single-nucleotide polymorphism markers segregating within the mapping population. Temperature sensitivity was assessed using digital image analysis of colonies growing at two different temperatures. We identified four QTLs for temperature sensitivity, with unique QTLs found in each cross. One QTL had a logarithm of odds score >11 and contained only six candidate genes, including PBS2, encoding a mitogen-activated protein kinase kinase associated with low temperature tolerance in Saccharomyces cerevisiae. This and other QTLs showed evidence for pleiotropy among growth rate, melanization and growth morphology, suggesting that many traits can be correlated with thermal adaptation in fungi. Higher temperatures were highly correlated with a shift to filamentous growth among the progeny in both crosses. We show that thermal adaptation has a complex genetic architecture, with natural populations of Z. tritici harboring significant genetic variation for this trait. We conclude that Z. tritici populations have the potential to adapt rapidly to climate change and expand into new climatic zones.

Proteome catalog of Zymoseptoria tritici captured during pathogenesis in wheat

Zymoseptoria tritici is an economically important pathogen of wheat. However, the molecular basis of pathogenicity on wheat is still poorly understood. Here, we present a global survey of the proteins secreted by this fungus in the apoplast of resistant (cv. Shafir) and susceptible (cv. Obelisk) wheat cultivars after inoculation with reference Z. tritici strain IPO323. The fungal proteins present in apoplastic fluids were analyzed by gel electrophoresis and by data-independent acquisition liquid chromatography/mass spectrometry (LC/MS(E)) combined with data-dependent acquisition LC-MS/MS. Subsequent mapping mass spectrometry-derived peptide sequence data against the genome sequence of strain IPO323 identified 665 peptides in the MS(E) and 93 in the LC-MS/MS mode that matched to 85 proteins. The identified fungal proteins, including cell-wall degrading enzymes and proteases, might function in pathogenicity, but the functions of many remain unknown. Most fungal proteins accumulated in cv. Obelisk at the onset of necrotrophy. This inventory provides an excellent basis for future detailed studies on the role of these genes and their encoded proteins during pathogenesis in wheat.

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.

Genome-wide gene expression dynamics of the fungal pathogen Dothistroma septosporum throughout its infection cycle of the gymnosperm host Pinus radiata

We present genome-wide gene expression patterns as a time series through the infection cycle of the fungal pine needle blight pathogen, Dothistroma septosporum, as it invades its gymnosperm host, Pinus radiata. We determined the molecular changes at three stages of the disease cycle: epiphytic/biotrophic (early), initial necrosis (mid) and mature sporulating lesion (late). Over 1.7 billion combined plant and fungal reads were sequenced to obtain 3.2 million fungal-specific reads, which comprised as little as 0.1% of the sample reads early in infection. This enriched dataset shows that the initial biotrophic stage is characterized by the up-regulation of genes encoding fungal cell wall-modifying enzymes and signalling proteins. Later necrotrophic stages show the up-regulation of genes for secondary metabolism, putative effectors, oxidoreductases, transporters and starch degradation. This in-depth through-time transcriptomic study provides our first snapshot of the gene expression dynamics that characterize infection by this fungal pathogen in its gymnosperm host.

The Impact of Recombination Hotspots on Genome Evolution of a Fungal Plant Pathogen

Recombination has an impact on genome evolution by maintaining chromosomal integrity, affecting the efficacy of selection, and increasing genetic variability in populations. Recombination rates are a key determinant of the coevolutionary dynamics between hosts and their pathogens. Historic recombination events created devastating new pathogens, but the impact of ongoing recombination in sexual pathogens is poorly understood. Many fungal pathogens of plants undergo regular sexual cycles, and sex is considered to be a major factor contributing to virulence. We generated a recombination map at kilobase-scale resolution for the haploid plant pathogenic fungus Zymoseptoria tritici. To account for intraspecific variation in recombination rates, we constructed genetic maps from two independent crosses. We localized a total of 10,287 crossover events in 441 progeny and found that recombination rates were highly heterogeneous within and among chromosomes. Recombination rates on large chromosomes were inversely correlated with chromosome length. Short accessory chromosomes often lacked evidence for crossovers between parental chromosomes. Recombination was concentrated in narrow hotspots that were preferentially located close to telomeres. Hotspots were only partially conserved between the two crosses, suggesting that hotspots are short-lived and may vary according to genomic background. Genes located in hotspot regions were enriched in genes encoding secreted proteins. Population resequencing showed that chromosomal regions with high recombination rates were strongly correlated with regions of low linkage disequilibrium. Hence, genes in pathogen recombination hotspots are likely to evolve faster in natural populations and may represent a greater threat to the host.

RNA-Seq Analysis of the Expression of Genes Encoding Cell Wall Degrading Enzymes during Infection of Lupin (Lupinus angustifolius) by Phytophthora parasitica

RNA-Seq analysis has shown that over 60% (12,962) of the predicted transcripts in the Phytophthora parasitica genome are expressed during the first 60 h of lupin root infection. The infection transcriptomes included 278 of the 431 genes encoding P. parasitica cell wall degrading enzymes. The transcriptome data provide strong evidence of global transcriptional cascades of genes whose encoded proteins target the main categories of plant cell wall components. A major cohort of pectinases is predominantly expressed early but as infection progresses, the transcriptome becomes increasingly dominated by transcripts encoding cellulases, hemicellulases, β-1,3-glucanases and glycoproteins. The most highly expressed P. parasitica carbohydrate active enzyme gene contains two CBM1 cellulose binding modules and no catalytic domains. The top 200 differentially expressed genes include β-1,4-glucosidases, β-1,4-glucanases, β-1,4-galactanases, a β-1,3-glucanase, an α-1,4-polygalacturonase, a pectin deacetylase and a pectin methylesterase. Detailed analysis of gene expression profiles provides clues as to the order in which linkages within the complex carbohydrates may come under attack. The gene expression profiles suggest that (i) demethylation of pectic homogalacturonan occurs before its deacetylation; (ii) cleavage of the backbone of pectic rhamnogalacturonan I precedes digestion of its side chains; (iii) early attack on cellulose microfibrils by non-catalytic cellulose-binding proteins and enzymes with auxiliary activities may facilitate subsequent attack by glycosyl hydrolases and enzymes containing CBM1 cellulose-binding modules; (iv) terminal hemicellulose backbone residues are targeted after extensive internal backbone cleavage has occurred; and (v) the carbohydrate chains on glycoproteins are degraded late in infection. A notable feature of the P. parasitica infection transcriptome is the high level of transcription of genes encoding enzymes that degrade β-1,3-glucanases during middle and late stages of infection. The results suggest that high levels of β-1,3-glucanases may effectively degrade callose as it is produced by the plant during the defence response.

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

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

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

Genome-Wide Analysis in Three Fusarium Pathogens Identifies Rapidly Evolving Chromosomes and Genes Associated with Pathogenicity

Pathogens and hosts are in an ongoing arms race and genes involved in host–pathogen interactions are likely to undergo diversifying selection. Fusarium plant pathogens have evolved diverse infection strategies, but how they interact with their hosts in the biotrophic infection stage remains puzzling. To address this, we analyzed the genomes of three Fusarium plant pathogens for genes that are under diversifying selection. We found a two-speed genome structure both on the chromosome and gene group level. Diversifying selection acts strongly on the dispensable chromosomes in Fusarium oxysporum f. sp. lycopersici and on distinct core chromosome regions in Fusarium graminearum, all of which have associations with virulence. Members of two gene groups evolve rapidly, namely those that encode proteins with an N-terminal [SG]-P-C-[KR]-P sequence motif and proteins that are conserved predominantly in pathogens. Specifically, 29 F. graminearum genes are rapidly evolving, in planta induced and encode secreted proteins, strongly pointing toward effector function. In summary, diversifying selection in Fusarium is strongly reflected as genomic footprints and can be used to predict a small gene set likely to be involved in host–pathogen interactions for experimental verification.

Genomic analysis of a ginger pathogen Bacillus pumilus providing the understanding to the pathogenesis and the novel control strategy

Bacillus pumilus has been widely identified as a pathogen of plant and human, while the genetic information is rarely available for pathogenic B. pumilus strains. B. pumilus GR8 is a pathogen that causes ginger rhizome rot disease by invading ginger rhizome parenchymatous tissues, growing in the extracellular space, and producing plant cell wall-degrading enzymes to destroy ginger cells. In this study, the genome of GR8 was sequenced and characterized. This genome was the third completely sequenced genome of the B. pumilus species, and it exhibited high similarity to the genome of the B. pumilus strain B6033. The genome of GR8 was 3.67 Mb in length and encoded 3,713 putative ORFs. Among these predicted proteins, numerous plant cell wall-degrading enzymes and several proteins associated with invading and adapting to the environment in the extracellular space of the ginger rhizome parenchymatous tissue were found. The GR8 genome contained only one restriction-modification system and no CRISPR/Cas system. The lack of phage-resistant system suggested that phages might be potential agents for the control of GR8. The genomic analysis of GR8 provided the understanding to the pathogenesis and the phage-control strategy of pathogenic B. pumilus strains.

QTL mapping of fungicide sensitivity reveals novel genes and pleiotropy with melanization in the pathogen Zymoseptoria tritici

A major problem associated with the intensification of agriculture is the emergence of fungicide resistance. Azoles are ergosterol biosynthesis inhibitors that have been widely used in agriculture and medicine since the 1970s, leading to emergence of increasingly resistant fungal populations. The known genetic mechanisms underlying lower azole sensitivity include mutations affecting the CYP51 gene that encodes the target protein, but in many cases azole resistance is a more complex trait with an unknown genetic basis. We used quantitative trait locus (QTL) mapping to identify genes affecting azole sensitivity in two crosses of Zymoseptoria tritici, the most damaging wheat pathogen in Europe. Restriction site associated DNA sequencing (RADseq) was used to genotype 263 (cross 1) and 261 (cross 2) progeny at ∼ 8500 single nucleotide polymorphisms (SNP) and construct two dense linkage maps. Azole sensitivity was assessed using high-throughput digital image analysis of colonies growing on Petri dishes with or without the fungicide propiconazole. We identified three QTLs for azole sensitivity, including two that contained novel fungicide sensitivity genes. One of these two QTLs contained only 16 candidate genes, among which four most likely candidates were identified. The third QTL contained ERG6, encoding another protein involved in ergosterol biosynthesis. Known genes in QTLs affecting colony growth included CYP51 and PKS1, a gene affecting melanization in Z. tritici. PKS1 showed compelling evidence for pleiotropy, with a rare segregating allele that increased melanization while decreasing growth rate and propiconazole sensitivity. This study resolved the genetic architecture of an important agricultural trait and led to identification of novel genes that are likely to affect azole sensitivity in Z. tritici. It also provided insight into fitness costs associated with lowered azole sensitivity and suggests a novel fungicide mixture strategy.

Transcriptome and metabolite profiling of the infection cycle of Zymoseptoria tritici on wheat reveals a biphasic interaction with plant immunity involving differential pathogen chromosomal contributions and a variation on the hemibiotrophic lifestyle definition

The hemibiotrophic fungus Zymoseptoria tritici causes Septoria tritici blotch disease of wheat (Triticum aestivum). Pathogen reproduction on wheat occurs without cell penetration, suggesting that dynamic and intimate intercellular communication occurs between fungus and plant throughout the disease cycle. We used deep RNA sequencing and metabolomics to investigate the physiology of plant and pathogen throughout an asexual reproductive cycle of Z. tritici on wheat leaves. Over 3,000 pathogen genes, more than 7,000 wheat genes, and more than 300 metabolites were differentially regulated. Intriguingly, individual fungal chromosomes contributed unequally to the overall gene expression changes. Early transcriptional down-regulation of putative host defense genes was detected in inoculated leaves. There was little evidence for fungal nutrient acquisition from the plant throughout symptomless colonization by Z. tritici, which may instead be utilizing lipid and fatty acid stores for growth. However, the fungus then subsequently manipulated specific plant carbohydrates, including fructan metabolites, during the switch to necrotrophic growth and reproduction. This switch coincided with increased expression of jasmonic acid biosynthesis genes and large-scale activation of other plant defense responses. Fungal genes encoding putative secondary metabolite clusters and secreted effector proteins were identified with distinct infection phase-specific expression patterns, although functional analysis suggested that many have overlapping/redundant functions in virulence. The pathogenic lifestyle of Z. tritici on wheat revealed through this study, involving initial defense suppression by a slow-growing extracellular and nutritionally limited pathogen followed by defense (hyper) activation during reproduction, reveals a subtle modification of the conceptual definition of hemibiotrophic plant infection.

Extracellular peptidases of the cereal pathogen Fusarium graminearum

The plant pathogenic fungus Fusarium graminearum (Fgr) creates economic and health risks in cereals agriculture. Fgr causes head blight (or scab) of wheat and stalk rot of corn, reducing yield, degrading grain quality, and polluting downstream food products with mycotoxins. Fungal plant pathogens must secrete proteases to access nutrition and to breakdown the structural protein component of the plant cell wall. Research into the proteolytic activity of Fgr is hindered by the complex nature of the suite of proteases secreted. We used a systems biology approach comprising genome analysis, transcriptomics and label-free quantitative proteomics to characterize the peptidases deployed by Fgr during growth. A combined analysis of published microarray transcriptome datasets revealed seven transcriptional groupings of peptidases based on in vitro growth, in planta growth, and sporulation behaviors. A high resolution mass spectrometry-based proteomics analysis defined the extracellular proteases secreted by F. graminearum. A meta-classification based on sequence characters and transcriptional/translational activity in planta and in vitro provides a platform to develop control strategies that target Fgr peptidases.