Defining the predicted protein secretome of the fungal wheat leaf pathogen Mycosphaerella graminicola

Identification and analysis of the secretome of plant pathogenic fungi reveals lifestyle adaptation

The secretory proteome plays an important role in the pathogenesis of phytopathogenic fungi. However, the relationship between the large-scale secretome of phytopathogenic fungi and their lifestyle is not fully understood. In the present study, the secretomes of 150 plant pathogenic fungi were predicted and the characteristics associated with different lifestyles were investigated. In total, 94,974 secreted proteins (SPs) were predicted from these fungi. The number of the SPs ranged from 64 to 1,662. Among these fungi, hemibiotrophic fungi had the highest number (average of 970) and proportion (7.1%) of SPs. Functional annotation showed that hemibiotrophic and necrotroph fungi, differ from biotrophic and symbiotic fungi, contained much more carbohydrate enzymes, especially polysaccharide lyases and carbohydrate esterases. Furthermore, the core and lifestyle-specific SPs orthogroups were identified. Twenty-seven core orthogroups contained 16% of the total SPs and their motif function annotation was represented by serine carboxypeptidase, carboxylesterase and asparaginase. In contrast, 97 lifestyle-specific orthogroups contained only 1% of the total SPs, with diverse functions such as PAN_AP in hemibiotroph-specific and flavin monooxygenases in necrotroph-specific. Moreover, obligate biotrophic fungi had the largest number of effectors (average of 150), followed by hemibiotrophic fungi (average of 120). Among these effectors, 4,155 had known functional annotation and pectin lyase had the highest proportion in the functionally annotated effectors. In addition, 32 sets of RNA-Seq data on pathogen-host interactions were collected and the expression levels of SPs were higher than that of non-SPs, and the expression level of effector genes was higher in biotrophic and hemibiotrophic fungi than in necrotrophic fungi, while secretase genes were highly expressed in necrotrophic fungi. Finally, the secretory activity of five predicted SPs from Setosphearia turcica was experimentally verified. In conclusion, our results provide a foundation for the study of pathogen-host interaction and help us to understand the fungal lifestyle adaptation.

Inactivation of a candidate effector gene of Zymoseptoria tritici affects its sporulation

BACKGROUND

Wheat is one of the most important staple crops produced worldwide. Its susceptibility to plant diseases reduces its production significantly. One of the most important diseases of wheat is septoria tritici blotch, a devastating disease observed in fields with wet and temperate conditions. Z. tritici secretes effector proteins to influence the host's defense mechanisms, as is typical of plant pathogens. In this investigation, we evaluated the pathogenicity of some Zymoseptoria tritici effector candidate genes having a signal peptide for secretion with no known function.

METHODS AND RESULTS

Three genes named Mycgr3G104383, Mycgr3G104444 and Mycgr3G105826 were knocked out separately through homologous recombination, generating Z. tritici IPO323 mutants lacking the functional copy of the corresponding genes. While KO1 and KO3 mutants did not show any significant differences during phenotypic and virulence investigations, the KO2 mutant generated exclusively macropycnidiospores in artificial media, different from wild-type IPO323 which produce only micropycidiospores. The mycelial growth capability of KO2 was also severely attenuated in all of the investigated growth conditions. These changes were observed independent of growth media and growth temperatures, implying that changes were genetic and inherited through generations. Virulence of knockout mutants in wheat leaves was observed to be similar to the wild-type IPO323.

CONCLUSION

Understanding the biology of Z. tritici and its interactions with wheat will reveal new strategies to fight septoria tritici blotch, enabling breeding wheat cultivars resistant to a broader spectrum of Z. tritici strains. Furthermore, gene knockout via homologous recombination proved to be a powerful tool for discovering novel gene functions.

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

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

Secretome Analysis of Macrophomina phaseolina Identifies an Array of Putative Virulence Factors Responsible for Charcoal Rot Disease in Plants

Macrophomina phaseolina is a global devastating necrotrophic fungal pathogen. It causes charcoal rot disease in more than 500 host plants including major food crops, pulse crops, fiber crops, and oil crops. Despite having the whole-genome sequence of M. phaseolina, understanding the M. phaseolina genome-based plant–pathogen interactions is limited in the absence of direct experimental proof of secretion. Thus, it is essential to understand the host–microbe interaction and the disease pathogenesis, which can ensure global agricultural crop production and security. An in silico–predicted secretome of M. phaseolina is unable to represent the actual secretome. We could identify 117 proteins present in the secretome of M. phaseolina using liquid chromatography–electrospray ionization–tandem mass spectrometry. Data are available via ProteomeXchange with identifier PXD032749. An array of putative virulence factors of M. phaseolina were identified in the present study using solid-state culture. Similar virulence factors have been reported in other plant pathogenic fungi also. Among the secretory fungal proteins with positive economic impacts, lignocellulolytic enzymes are of prime importance. Further, we validated our results by detecting the cell wall–degrading enzymes xylanase, endoglucanase, and amylase in the secretome of M. phaseolina. The present study may provide a better understanding about the necrotrophic fungi M. phaseolina, which modulate the host plant defense barriers using secretory proteins.

Cross-Kingdom Gene Coexpression Analysis Using a Stemphylium botryosum-Lens ervoides System Revealed Plasticity of Intercommunication Between the Pathogen Secretome and the Host Immune Systems

Necrotrophic pathogens are responsible for significant declines in crop yield and quality worldwide. During the infection process, a pathogen releases a series of secretory proteins to counteract the plant immune system, and this interaction of pathogen and host molecules determines whether the pathogen will successfully invade the host plant tissues. In this study, we adopted co-transcriptomic approaches to analyze the Lens ervoides-Stemphylium botryosum system, with a focus on 1,216 fungal genes coding for secretory proteins and 8,810 disease-responsive genes of the host 48, 96, and 144 h postinoculation, captured in two F₉ recombinant inbred lines (RILs) displaying contrasting disease responses. By constructing in planta gene coexpression networks (GCNs) for S. botryosum, we found that the pathogen tended to co-upregulate genes regulating cell wall degradation enzymes, effectors, oxidoreductases, and peptidases to a much higher degree in the susceptible host LR-66-577 than in the resistant RIL LR-66-637, indicating that the promotion of these digestive enzymes and toxins increased S. botryosum virulence. Construction of cross-kingdom GCNs between pathogen and plant for the two RILs revealed that the co-upregulation of these fungal digestive enzymes and toxins simultaneously promoted a series of defense responses such as redox change, expression of membrane-related genes and serine/threonine kinase, and stress and disease responses in the susceptible RIL which was not observed in the resistant RIL, indicating that these activities exacerbated susceptibility to S. botryosum.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

In silico approach to predict pancreatic β-cells classically secreted proteins

Pancreatic β-cells, residents of the islets of Langerhans, are the unique insulin-producers in the body. Their physiology is a topic of intensive studies aiming to understand the biology of insulin production and its role in diabetes pathology. However, investigations about these cells' subset of secreted proteins, the secretome, are surprisingly scarce and a list describing islet/β-cell secretome upon glucose-stimulation is not yet available. In silico predictions of secretomes are an interesting approach that can be employed to forecast proteins likely to be secreted. In this context, using the rationale behind classical secretion of proteins through the secretory pathway, a Python tool capable of predicting classically secreted proteins was developed. This tool was applied to different available proteomic data (human and rodent islets, isolated β-cells, β-cell secretory granules, and β-cells supernatant), filtering them in order to selectively list only classically secreted proteins. The method presented here can retrieve, organize, search and filter proteomic lists using UniProtKB as a central database. It provides analysis by overlaying different sets of information, filtering out potential contaminants and clustering the identified proteins into functional groups. A range of 70-92% of the original proteomes analyzed was reduced generating predicted secretomes. Islet and β-cell signal peptide-containing proteins, and endoplasmic reticulum-resident proteins were identified and quantified. From the predicted secretomes, exemplary conservational patterns were inferred, as well as the signaling pathways enriched within them. Such a technique proves to be an effective approach to reduce the horizon of plausible targets for drug development or biomarkers identification.

Comparative analysis of extracellular proteomes reveals putative effectors of the boxwood blight pathogens, Calonectria henricotiae and C. pseudonaviculata

Calonectria henricotiae (Che) and C. pseudonaviculata (Cps) are destructive fungal pathogens causing boxwood blight, a persistent threat to horticultural production, landscape industries, established gardens, and native ecosystems. Although extracellular proteins including effectors produced by fungal pathogens are known to play a fundamental role in pathogenesis, the composition of Che and Cps extracellular proteins has not been examined. Using liquid chromatography-tandem mass spectrometry (LC-MS/MS) and bioinformatics prediction tools, 630 extracellular proteins and 251 cell membrane proteins of Che and Cps were identified in the classical secretion pathway in the present study. In the non-classical secretion pathway, 79 extracellular proteins were identified. The cohort of proteins belonged to 364 OrthoMCL clusters, with the majority (62%) present in both species, and a subset unique to Che (19%) and Cps (20%). These extracellular proteins were predicted to play important roles in cell structure, regulation, metabolism, and pathogenesis. A total of 124 proteins were identified as putative effectors. Many of them are orthologs of proteins with documented roles in suppressing host defense and facilitating infection processes in other pathosystems, such as SnodProt1-like proteins in the OrthoMCL cluster OG5_152723 and PhiA-like cell wall proteins in the cluster OG5_155754. This exploratory study provides a repository of secreted proteins and putative effectors that can provide insights into the virulence mechanisms of the boxwood blight pathogens.

A small secreted protein from Zymoseptoria tritici interacts with a wheat E3 ubiquitin ligase to promote disease

Septoria tritici blotch (STB), caused by the ascomycete fungus Zymoseptoria tritici, is a major threat to wheat production worldwide. The Z. tritici genome encodes many small secreted proteins (ZtSSPs) that are likely to play a key role in the successful colonization of host tissues. However, few of these ZtSSPs have been functionally characterized for their role during infection. In this study, we identified and characterized a small, conserved cysteine-rich secreted effector from Z. tritici which has homologues in other plant pathogens in the Dothideomycetes. ZtSSP2 was expressed throughout Z. tritici infection in wheat, with the highest levels observed early during infection. A yeast two-hybrid assay revealed an interaction between ZtSSP2 and wheat E3 ubiquitin ligase (TaE3UBQ) in yeast, and this was further confirmed in planta using bimolecular fluorescence complementation and co-immunoprecipitation. Down-regulation of this wheat E3 ligase using virus-induced gene silencing increased the susceptibility of wheat to STB. Together, these results suggest that TaE3UBQ is likely to play a role in plant immunity to defend against Z. tritici.

A small secreted protein from Zymoseptoria tritici interacts with a wheat E3 ubiquitin ligase to promote disease

Septoria tritici blotch (STB), caused by the ascomycete fungus Zymoseptoria tritici, is a major threat to wheat production worldwide. The Z. tritici genome encodes many small secreted proteins (ZtSSPs) that are likely to play a key role in the successful colonization of host tissues. However, few of these ZtSSPs have been functionally characterized for their role during infection. In this study, we identified and characterized a small, conserved cysteine-rich secreted effector from Z. tritici which has homologues in other plant pathogens in the Dothideomycetes. ZtSSP2 was expressed throughout Z. tritici infection in wheat, with the highest levels observed early during infection. A yeast two-hybrid assay revealed an interaction between ZtSSP2 and wheat E3 ubiquitin ligase (TaE3UBQ) in yeast, and this was further confirmed in planta using bimolecular fluorescence complementation and co-immunoprecipitation. Down-regulation of this wheat E3 ligase using virus-induced gene silencing increased the susceptibility of wheat to STB. Together, these results suggest that TaE3UBQ is likely to play a role in plant immunity to defend against Z. tritici.

Extracellular vesicles from the apoplastic fungal wheat pathogen Zymoseptoria tritici

BACKGROUND

The fungal pathogen Zymoseptoria tritici is a significant constraint to wheat production in temperate cropping regions around the world. Despite its agronomic impacts, the mechanisms allowing the pathogen to asymptomatically invade and grow in the apoplast of wheat leaves before causing extensive host cell death remain elusive. Given recent evidence of extracellular vesicles (EVs)-secreted, membrane-bound nanoparticles containing molecular cargo-being implicated in extracellular communication between plants and fungal pathogen, we have initiated an in vitro investigation of EVs from this apoplastic fungal wheat pathogen. We aimed to isolate EVs from Z. tritici broth cultures and examine their protein composition in relation to the soluble protein in the culture filtrate and to existing fungal EV proteomes.

RESULTS

Zymoseptoria tritici EVs were isolated from broth culture filtrates using differential ultracentrifugation (DUC) and examined with transmission electron microscopy (TEM) and nanoparticle tracking analysis (NTA). Z. tritici EVs were observed as a heterogeneous population of particles, with most between 50 and 250 nm. These particles were found in abundance in the culture filtrates of viable Z. tritici cultures, but not heat-killed cultures incubated for an equivalent time and of comparable biomass. Bottom-up proteomic analysis using LC-MS/MS, followed by stringent filtering revealed 240 Z. tritici EV proteins. These proteins were distinct from soluble proteins identified in Z. tritici culture filtrates, but were similar to proteins identified in EVs from other fungi, based on sequence similarity analyses. Notably, a putative marker protein recently identified in Candida albicans EVs was also consistently detected in Z. tritici EVs.

CONCLUSION

We have shown EVs can be isolated from the devastating fungal wheat pathogen Z. tritici and are similar to protein composition to previously characterised fungal EVs. EVs from human pathogenic fungi are implicated in virulence, but the role of EVs in the interaction of phytopathogenic fungi and their hosts is unknown. These in vitro analyses provide a basis for expanding investigations of Z. tritici EVs in planta, to examine their involvement in the infection process of this apoplastic wheat pathogen and more broadly, advance understanding of noncanonical secretion in filamentous plant pathogens.

Wheat Encodes Small, Secreted Proteins That Contribute to Resistance to Septoria Tritici Blotch

During plant-pathogen interactions, pathogens secrete many rapidly evolving, small secreted proteins (SSPs) that can modify plant defense and permit pathogens to colonize plant tissue. The fungal pathogen Zymoseptoria tritici is the causal agent of Septoria tritici blotch (STB), one of the most important foliar diseases of wheat, globally. Z. tritici is a strictly apoplastic pathogen that can secrete numerous proteins into the apoplast of wheat leaves to promote infection. We sought to determine if, during STB infection, wheat also secretes small proteins into the apoplast to mediate the recognition of pathogen proteins and/or induce defense responses. To explore this, we developed an SSP-discovery pipeline to identify small, secreted proteins from wheat genomic data. Using this pipeline, we identified 6,998 SSPs, representing 2.3% of all proteins encoded by the wheat genome. We then mined a microarray dataset, detailing a resistant and susceptible host response to STB, and identified 141 Z. tritici- responsive SSPs, representing 4.7% of all proteins encoded by Z. tritici - responsive genes. We demonstrate that a subset of these SSPs have a functional signal peptide and can interact with Z. tritici SSPs. Transiently silencing two of these wheat SSPs using virus-induced gene silencing (VIGS) shows an increase in susceptibility to STB, confirming their role in defense against Z. tritici.

EffHunter: A Tool for Prediction of Effector Protein Candidates in Fungal Proteomic Databases

Pathogens are able to deliver small-secreted, cysteine-rich proteins into plant cells to enable infection. The computational prediction of effector proteins remains one of the most challenging areas in the study of plant fungi interactions. At present, there are several bioinformatic programs that can help in the identification of these proteins; however, in most cases, these programs are managed independently. Here, we present EffHunter, an easy and fast bioinformatics tool for the identification of effectors. This predictor was used to identify putative effectors in 88 proteomes using characteristics such as size, cysteine residue content, secretion signal and transmembrane domains.

Prediction of the secretomes of endophytic and nonendophytic fungi reveals similarities in host plant infection and colonization strategies

Endophytic fungi are microorganisms that inhabit internal plant tissues without causing apparent damage. During the infection process, both endophytic and phytopathogenic fungi secrete proteins to resist or supplant the plant's defense mechanisms. This study analyzed the predicted secretomes of six species of endophytic fungi and compared them with predicted secretomes of eight fungal species with different lifestyles: saprophytic, necrotrophic, hemibiotrophic, and biotrophic. The sizes of the predicted secretomes varied from 260 to 1640 proteins, and the predicted secretomes have a wide diversity of CAZymes, proteases, and conserved domains. Regarding the CAZymes in the secretomes of the analyzed fungi, the most abundant CAZyme families were glycosyl hydrolase and serine proteases. Several predicted proteins have characteristics similar to those found in small, secreted proteins with effector characteristics (SSPEC). The most abundant conserved domains, besides those found in the SSPEC, have oxidation activities, indicating that these proteins can protect the fungus against oxidative stress, against domains with protease activity, which may be involved in the mechanisms of nutrition, or against lytic enzymes secreted by the host plant. This study demonstrates that secretomes of endophytic and nonendophytic fungi share an arsenal of proteins important in the process of infection and colonization of host plants.

Conservation and expansion of a necrosis-inducing small secreted protein family from host-variable phytopathogens of the Sclerotiniaceae

Fungal effector proteins facilitate host-plant colonization and have generally been characterized as small secreted proteins (SSPs). We classified and functionally tested SSPs from the secretomes of three closely related necrotrophic phytopathogens: Ciborinia camelliae, Botrytis cinerea, and Sclerotinia sclerotiorum. Alignment of predicted SSPs identified a large protein family that share greater than 41% amino acid identity and that have key characteristics of previously described microbe-associated molecular patterns (MAMPs). Strikingly, 73 of the 75 SSP family members were predicted within the secretome of the host-specialist C. camelliae with single-copy homologs identified in the secretomes of the host generalists S. sclerotiorum and B. cinerea. To explore the potential function of this family of SSPs, 10 of the 73 C. camelliae proteins, together with the single-copy homologs from S. sclerotiorum (SsSSP3) and B. cinerea (BcSSP2), were cloned and expressed as recombinant proteins. Infiltration of SsSSP3 and BcSSP2 into host tissue induced rapid necrosis. In contrast, only one of the 10 tested C. camelliae SSPs was able to induce a limited amount of necrosis. Analysis of chimeric proteins consisting of domains from both a necrosis-inducing and a non-necrosis-inducing SSP demonstrated that the C-terminus of the S. sclerotiorum SSP is essential for necrosis-inducing function. Deletion of the BcSSP2 homolog from B. cinerea did not affect growth or pathogenesis. Thus, this research uncovered a family of highly conserved SSPs present in diverse ascomycetes that exhibit contrasting necrosis-inducing functions.

Genomes and secretomes of Ascomycota fungi reveal diverse functions in plant biomass decomposition and pathogenesis

BACKGROUND

The dominant fungi in arid grasslands and shrublands are members of the Ascomycota phylum. Ascomycota fungi are important drivers in carbon and nitrogen cycling in arid ecosystems. These fungi play roles in soil stability, plant biomass decomposition, and endophytic interactions with plants. They may also form symbiotic associations with biocrust components or be latent saprotrophs or pathogens that live on plant tissues. However, their functional potential in arid soils, where organic matter, nutrients and water are very low or only periodically available, is poorly characterized.

RESULTS

Five Ascomycota fungi were isolated from different soil crust microhabitats and rhizosphere soils around the native bunchgrass Pleuraphis jamesii in an arid grassland near Moab, UT, USA. Putative genera were Coniochaeta, isolated from lichen biocrust, Embellisia from cyanobacteria biocrust, Chaetomium from below lichen biocrust, Phoma from a moss microhabitat, and Aspergillus from the soil. The fungi were grown in replicate cultures on different carbon sources (chitin, native bunchgrass or pine wood) relevant to plant biomass and soil carbon sources. Secretomes produced by the fungi on each substrate were characterized. Results demonstrate that these fungi likely interact with primary producers (biocrust or plants) by secreting a wide range of proteins that facilitate symbiotic associations. Each of the fungal isolates secreted enzymes that degrade plant biomass, small secreted effector proteins, and proteins involved in either beneficial plant interactions or virulence. Aspergillus and Phoma expressed more plant biomass degrading enzymes when grown in grass- and pine-containing cultures than in chitin. Coniochaeta and Embellisia expressed similar numbers of these enzymes under all conditions, while Chaetomium secreted more of these enzymes in grass-containing cultures.

CONCLUSIONS

This study of Ascomycota genomes and secretomes provides important insights about the lifestyles and the roles that Ascomycota fungi likely play in arid grassland, ecosystems. However, the exact nature of those interactions, whether any or all of the isolates are true endophytes, latent saprotrophs or opportunistic phytopathogens, will be the topic of future studies.

A novel starch-binding laccase from the wheat pathogen Zymoseptoria tritici highlights the functional diversity of ascomycete laccases

BACKGROUND

Laccases are multicopper oxidases, which are assigned into auxiliary activity family 1 (AA1) in the CAZy database. These enzymes, catalyzing the oxidation of phenolic and nonphenolic substrates coupled to reduction of O2 to H2O, are increasingly attractive as eco-friendly oxidation biocatalysts. Basidiomycota laccases are well characterized due to their potential in de-lignification of lignocellulose. By contrast, insight into the biochemical diversity of Ascomycota counterparts from saprophytes and plant pathogens is scarce.

RESULTS

Here, we report the properties of the laccase from the major wheat pathogen Zymoseptoria tritici (ZtrLac1A), distinguished from common plant fungal pathogens by an apoplastic infection strategy. We demonstrate that ZtrLac1A is appended to a functional starch-binding module and displays an activity signature disfavoring relatively apolar phenolic redox mediators as compared to the related biochemically characterized laccases. By contrast, the redox potential of ZtrLac1A (370 mV vs. SHE) is similar to ascomycetes counterparts. The atypical specificity is consistent with distinctive sequence substitutions and insertions in loops flanking the T1 site and the enzyme C-terminus compared to characterized laccases.

CONCLUSIONS

ZtrLac1A is the first reported modular laccase appended to a functional starch-specific carbohydrate binding module of family 20 (CBM20). The distinct specificity profile of ZtrLac1A correlates to structural differences in the active site region compared to previously described ascomycetes homologues. These differences are also highlighted by the clustering of the sequence of ZtrLac1A in a distinct clade populated predominantly by plant pathogens in the phylogenetic tree of AA1 laccases. The possible role of these laccases in vivo merits further investigations. These findings expand our toolbox of laccases for green oxidation and highlight the binding functionality of CBM-appended laccases as versatile immobilization tags.

Morphological changes in response to environmental stresses in the fungal plant pathogen Zymoseptoria tritici

During their life cycles, pathogens have to adapt to many biotic and abiotic environmental stresses to maximize their overall fitness. Morphological transitions are one of the least understood of the many strategies employed by fungal plant pathogens to adapt to constantly changing environments, even though different morphotypes may play important biological roles. Here, we first show that blastospores (the "yeast-like" form of the pathogen typically known only under laboratory conditions) can form from germinated pycnidiospores (asexual spores) on the surface of wheat leaves, suggesting that this morphotype can play an important role in the natural history of Z. tritici. Next, we characterized the morphological responses of this fungus to a series of environmental stresses to understand the effects of changing environments on fungal morphology and adaptation. All tested stresses induced morphological changes, but different responses were found among four strains. We discovered that Z. tritici forms chlamydospores and demonstrated that these structures are better able to survive extreme cold, heat and drought than other cell types. Finally, a transcriptomic analysis showed that morphogenesis and the expression of virulence factors are co-regulated in this pathogen. Our findings illustrate how changing environmental conditions can affect cellular morphology and lead to the formation of new morphotypes, with each morphotype having a potential impact on both pathogen survival and disease epidemiology.

Whole genome and transcriptome analysis reveal adaptive strategies and pathogenesis of Calonectria pseudoreteaudii to Eucalyptus

BACKGROUND

Leaf blight caused by Calonectria spp. is one of the most destructive diseases to affect Eucalyptus nurseries and plantations. These pathogens mainly attack Eucalyptus, a tree with a diversity of secondary metabolites employed as defense-related phytoalexins. To unravel the fungal adaptive mechanisms to various phytoalexins, we examined the genome of C. pseudoreteaudii, which is one of the most aggressive pathogens in southeast Asia.

RESULTS

A 63.7 Mb genome with 14,355 coding genes of C. pseudoreteaudii were assembled. Genomic comparisons identified 1785 species-specific gene families in C. pseudoreteaudii. Most of them were not annotated and those annotated genes were enriched in peptidase activity, pathogenesis, oxidoreductase activity, etc. RNA-seq showed that 4425 genes were differentially expressed on the eucalyptus(the resistant cultivar E. grandis×E.camaldulensis M1) tissue induced medium. The annotation of GO term and KEGG pathway indicated that some of the differential expression genes were involved in detoxification and transportation, such as genes encoding ABC transporters, degrading enzymes of aromatic compounds and so on.

CONCLUSIONS

Potential genomic determinants of phytoalexin detoxification were identified in C. pseudoreteaudii by comparison with 13 other fungi. This pathogen seems to employ membrane transporters and degradation enzymes to detoxify Eucalyptus phytoalexins. Remarkably, the Calonectria genome possesses a surprising number of secondary metabolism backbone enzyme genes involving toxin biosynthesis. It is also especially suited for cutin and lignin degradation. This indicates that toxin and cell wall degrading enzymes may act important roles in the establishment of Calonectria leaf blight. This study provides further understanding on the mechanism of pathogenesis in Calonectria.

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.

The role of reactive oxygen in the development of Ramularia leaf spot disease in barley seedlings

BACKGROUND AND AIMS

Ramularia collo-cygni is an ascomycete fungus that colonizes barley primarily as a benign endophyte, although this interaction can become pathogenic, causing the disease Ramularia leaf spot (RLS). Factors, particularly reactive oxygen species, that resulted in the transition of the fungus from endophyte to necrotrophic parasite and the development of disease symptoms were investigated.

METHODS

Disease development in artificially inoculated seedlings of barley varieties varying in partial resistance to RLS was related to exposure to abiotic stress prior to inoculation. Histochemical and molecular analysis determined the effect of R. collo-cygni colonization on accumulation of reactive oxygen species and antioxidant gene expression. Development of RLS on barley lines defective in antioxidant enzymes and with altered redox status or non-functional chloroplasts was compared with the accumulation of fungal biomass to determine how these factors affect disease symptom expression.

KEY RESULTS

Exposure to abiotic stress increased symptom development in all susceptible and most partially resistant barley varieties, in association with greater hydrogen peroxide (H2O2) levels in leaves. Decreased activity of the antioxidant enzymes superoxide dismutase and catalase in transgenic and mutant plants had no effect on the disease transition, whereas manipulation of H2O2 levels during asymptomatic growth of the fungus increased disease symptoms in most susceptible varieties but not in partially resistant plants. Barley mutants that undergo rapid loss of green leaf area when infected by R. collo-cygni or albino mutants with non-functional chloroplasts showed reduced development of RLS symptoms.

CONCLUSIONS

These results imply that in seedlings the pathogenic transition of the normally endophytic fungus R. collo-cygni does not result from senescence as such, but rather is promoted by factors that result in changes to host reactive oxygen species. Barley varieties vary in the extent to which these factors promote RLS disease.

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.

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.

A Puccinia striiformis f. sp. tritici secreted protein activates plant immunity at the cell surface

Pathogens secrete effector proteins to suppress host immunity, mediate nutrient uptake and subsequently enable parasitism. However, on non-adapted hosts, effectors can be detected as non-self by host immune receptors and activate non-host immunity. Nevertheless, the molecular mechanisms of effector triggered non-host resistance remain unknown. Here, we report that a small cysteine-rich protein PstSCR1 from the wheat rust pathogen Puccinia striiformis f. sp. tritici (Pst) activates immunity in the non-host solanaceous model plant Nicotiana benthamiana. PstSCR1 homologs were found to be conserved in Pst, and in its closest relatives, Puccinia graminis f. sp. tritici and Puccinia triticina. When PstSCR1 was expressed in N. benthamiana with its signal peptide, it provoked the plant immune system, whereas no stimulation was observed when it was expressed without its signal peptide. PstSCR1 expression in N. benthamiana significantly reduced infection capacity of the oomycete pathogens. Moreover, apoplast-targeted PstSCR1 triggered plant cell death in a dose dependent manner. However, in Brassinosteroid insensitive 1-Associated Kinase 1 (SERK3/BAK1) silenced N. benthamiana, cell death was remarkably decreased. Finally, purified PstSCR1 protein activated defence related gene expression in N. benthamiana. Our results show that a Pst-secreted protein, PstSCR1 can activate surface mediated immunity in non-adapted hosts and contribute to non-host resistance.

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.

The complete genome sequence of the phytopathogenic fungus Sclerotinia sclerotiorum reveals insights into the genome architecture of broad host range pathogens

Sclerotinia sclerotiorum is a phytopathogenic fungus with over 400 hosts including numerous economically important cultivated species. This contrasts many economically destructive pathogens that only exhibit a single or very few hosts. Many plant pathogens exhibit a “two-speed” genome. So described because their genomes contain alternating gene rich, repeat sparse and gene poor, repeat-rich regions. In fungi, the repeat-rich regions may be subjected to a process termed repeat-induced point mutation (RIP). Both repeat activity and RIP are thought to play a significant role in evolution of secreted virulence proteins, termed effectors. We present a complete genome sequence of S. sclerotiorum generated using Single Molecule Real-Time Sequencing technology with highly accurate annotations produced using an extensive RNA sequencing data set. We identified 70 effector candidates and have highlighted their in planta expression profiles. Furthermore, we characterized the genome architecture of S. sclerotiorum in comparison to plant pathogens that exhibit “two-speed” genomes. We show that there is a significant association between positions of secreted proteins and regions with a high RIP index in S. sclerotiorum but we did not detect a correlation between secreted protein proportion and GC content. Neither did we detect a negative correlation between CDS content and secreted protein proportion across the S. sclerotiorum genome. We conclude that S. sclerotiorum exhibits subtle signatures of enhanced mutation of secreted proteins in specific genomic compartments as a result of transposition and RIP activity. However, these signatures are not observable at the whole-genome scale.

Apoplastic recognition of multiple candidate effectors from the wheat pathogen Zymoseptoria tritici in the nonhost plant Nicotiana benthamiana

The fungus Zymoseptoria tritici is a strictly apoplastic, host-specific pathogen of wheat leaves and causal agent of septoria tritici blotch (STB) disease. All other plants are considered nonhosts, but the mechanism of nonhost resistance (NHR) to Z. tritici has not been addressed previously. We sought to develop Nicotiana benthamiana as a system to study NHR against Z. tritici. Fluorescence microscopy and quantitative reverse transcription polymerase chain reactions were used to establish the interaction between Z. tritici and N. benthamiana. Agrobacterium-mediated transient expression was used to screen putative Z. tritici effector genes for recognition in N. benthamiana, and virus-induced gene silencing (VIGS) was employed to determine the role of two receptor-like kinases (RLKs), NbBAK1 and NbSOBIR1, in Z. tritici effector recognition. Numerous Z. tritici putative effectors (14 of 63 tested) induced cell death or chlorosis in N. benthamiana. For most, phenotypes were light-dependent and required effector secretion to the leaf apoplastic space. Moreover, effector-induced host cell death was dependent on NbBAK1 and NbSOBIR1. Our results indicate widespread recognition of apoplastic effectors from a wheat-infecting fungal pathogen in a taxonomically distant nonhost plant species presumably by cell surface immune receptors. This suggests that apoplastic recognition of multiple nonadapted pathogen effectors may contribute to NHR.

Linear Correlation Analysis of Zymoseptoria tritici Aggressiveness with In Vitro Growth Rate

Zymoseptoria tritici is a globally distributed plant-pathogenic fungus causing Septoria tritici blotch of wheat. In this study, the in vitro growth rates and aggressiveness of 141 genetically distinct isolates sampled from four wheat fields on three continents were assessed to determine the association of these two ecological parameters. Aggressiveness was assessed on two spring wheat cultivars (‘Toronit’ and ‘Greina’) in a greenhouse using percentages of leaf area covered by lesions and pycnidia. We found a positive correlation between aggressiveness of pathogen strains on the two cultivars, consistent with a quantitative and host-nonspecific interaction in this pathosystem. We also found a positive correlation between aggressiveness and average growth rate at two temperatures, suggesting that in vitro pathogen growth rate may make a significant contribution to pathogen aggressiveness.

Schizophyllum commune has an extensive and functional alternative splicing repertoire

Recent genome-wide studies have demonstrated that fungi possess the machinery to alternatively splice pre-mRNA. However, there has not been a systematic categorization of the functional impact of alternative splicing in a fungus. We investigate alternative splicing and its functional consequences in the model mushroom forming fungus Schizophyllum commune. Alternative splicing was demonstrated for 2,285 out of 12,988 expressed genes, resulting in 20% additional transcripts. Intron retentions were the most common alternative splicing events, accounting for 33% of all splicing events, and 43% of the events in coding regions. On the other hand, exon skipping events were rare in coding regions (1%) but enriched in UTRs where they accounted for 57% of the events. Specific functional groups, including transcription factors, contained alternatively spliced genes. Alternatively spliced transcripts were regulated differently throughout development in 19% of the 2,285 alternatively spliced genes. Notably, 69% of alternatively spliced genes have predicted alternative functionality by loss or gain of functional domains, or by acquiring alternative subcellular locations. S. commune exhibits more alternative splicing than any other studied fungus. Taken together, alternative splicing increases the complexity of the S. commune proteome considerably and provides it with a rich repertoire of alternative functionality that is exploited dynamically.

Secretome analysis of virulent Pyrenophora teres f. teres isolates

Pyrenophora teres f. teres (Ptt) causes net form net blotch disease of barley, partially by producing necrosis-inducing proteins. The protein profiles of the culture filtrates of 28 virulent isolates were compared by a combination of 2DE and 1D-PAGE with 105 spots and 51 bands chosen for analysis by liquid chromatography electrospray ionization tandem mass spectrometry. A total of 259 individual proteins were identified with 63 of these proteins being common to the selected virulent isolates. Ptt secretes a broad spectrum of proteins including cell wall degrading enzymes; virulence factors and effectors; proteins associated with fungal pathogenesis and development; and proteins related to oxidation-reduction processes. Potential virulence factors and effectors identified included proteins with glucosidase activity, ricin B and concanavalin A-like lectins, glucanases, spherulin, cutinase, pectin lyase, leucine-rich repeat protein, and ceratoplatanin. Small proteins with unknown function but cysteine-rich, common to effectors, were also identified. Differences in the secretion profile of the Ptt isolates have also provided important insight into the different mechanisms contributing to virulence and the development of net form net blotch symptoms.

Global transcriptional analysis suggests Lasiodiplodia theobromae pathogenicity factors involved in modulation of grapevine defensive response

BACKGROUND

Lasiodiplodia theobromae is a fungus of the Botryosphaeriaceae that causes grapevine vascular disease, especially in regions with hot climates. Fungi in this group often remain latent within their host and become virulent under abiotic stress. Transcriptional regulation analysis of L. theobromae exposed to heat stress (HS) was first carried out in vitro in the presence of grapevine wood (GW) to identify potential pathogenicity genes that were later evaluated for in planta expression.

RESULTS

A total of 19,860 de novo assembled transcripts were obtained, forty-nine per cent of which showed homology to the Botryosphaeriaceae fungi, Neofusicoccum parvum or Macrophomina phaseolina. Three hundred ninety-nine have homology with genes involved in pathogenic processes and several belonged to expanded gene families in others fungal grapevine vascular pathogens. Gene expression analysis showed changes in fungal metabolism of phenolic compounds; where genes encoding for enzymes, with the ability to degrade salicylic acid (SA) and plant phenylpropanoid precursors, were up-regulated during in vitro HS response, in the presence of GW. These results suggest that the fungal L-tyrosine catabolism pathway could help the fungus to remove phenylpropanoid precursors thereby evading the host defense response. The in planta up-regulation of salicylate hydroxylase, intradiol ring cleavage dioxygenase and fumarylacetoacetase encoding genes, further supported this hypothesis. Those genes were even more up-regulated in HS-stressed plants, suggesting that fungus takes advantage of the increased phenylpropanoid precursors produced under stress. Pectate lyase was up-regulated while a putative amylase was down-regulated in planta, this could be associated with an intercellular growth strategy during the first stages of colonization.

CONCLUSIONS

L. theobromae transcriptome was established and validated. Its usefulness was demonstrated through the identification of genes expressed during the infection process. Our results support the hypothesis that heat stress facilitates fungal colonization, because of the fungus ability to use the phenylpropanoid precursors and SA, both compounds known to control host defense.

The genome of the emerging barley pathogen Ramularia collo-cygni

Background

Ramularia collo-cygni is a newly important, foliar fungal pathogen of barley that causes the disease Ramularia leaf spot. The fungus exhibits a prolonged endophytic growth stage before switching life habit to become an aggressive, necrotrophic pathogen that causes significant losses to green leaf area and hence grain yield and quality.

Results

The R. collo-cygni genome was sequenced using a combination of Illumina and Roche 454 technologies. The draft assembly of 30.3 Mb contained 11,617 predicted gene models. Our phylogenomic analysis confirmed the classification of this ascomycete fungus within the family Mycosphaerellaceae, order Capnodiales of the class Dothideomycetes. A predicted secretome comprising 1053 proteins included redox-related enzymes and carbohydrate-modifying enzymes and proteases. The relative paucity of plant cell wall degrading enzyme genes may be associated with the stealth pathogenesis characteristic of plant pathogens from the Mycosphaerellaceae. A large number of genes associated with secondary metabolite production, including homologs of toxin biosynthesis genes found in other Dothideomycete plant pathogens, were identified.

Conclusions

The genome sequence of R. collo-cygni provides a framework for understanding the genetic basis of pathogenesis in this important emerging pathogen. The reduced complement of carbohydrate-degrading enzyme genes is likely to reflect a strategy to avoid detection by host defences during its prolonged asymptomatic growth. Of particular interest will be the analysis of R. collo-cygni gene expression during interactions with the host barley, to understand what triggers this fungus to switch from being a benign endophyte to an aggressive necrotroph.

Electronic supplementary material

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

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.

Structure of a Berberine Bridge Enzyme-Like Enzyme with an Active Site Specific to the Plant Family Brassicaceae

Berberine bridge enzyme-like (BBE-like) proteins form a multigene family (pfam 08031), which is present in plants, fungi and bacteria. They adopt the vanillyl alcohol-oxidase fold and predominantly show bi-covalent tethering of the FAD cofactor to a cysteine and histidine residue, respectively. The Arabidopsis thaliana genome was recently shown to contain genes coding for 28 BBE-like proteins, while featuring four distinct active site compositions. We determined the structure of a member of the AtBBE-like protein family (termed AtBBE-like 28), which has an active site composition that has not been structurally and biochemically characterized thus far. The most salient and distinguishing features of the active site found in AtBBE-like 28 are a mono-covalent linkage of a histidine to the 8α-position of the flavin-isoalloxazine ring and the lack of a second covalent linkage to the 6-position, owing to the replacement of a cysteine with a histidine. In addition, the structure reveals the interaction of a glutamic acid (Glu426) with an aspartic acid (Asp369) at the active site, which appear to share a proton. This arrangement leads to the delocalization of a negative charge at the active site that may be exploited for catalysis. The structure also indicates a shift of the position of the isoalloxazine ring in comparison to other members of the BBE-like family. The dioxygen surrogate chloride was found near the C(4a) position of the isoalloxazine ring in the oxygen pocket, pointing to a rapid reoxidation of reduced enzyme by dioxygen. A T-DNA insertional mutant line for AtBBE-like 28 results in a phenotype, that is characterized by reduced biomass and lower salt stress tolerance. Multiple sequence analysis showed that the active site composition found in AtBBE-like 28 is only present in the Brassicaceae, suggesting that it plays a specific role in the metabolism of this plant family.

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.

Overview of genomic and bioinformatic resources for Zymoseptoria tritici

Zymoseptoria tritici (syn. Mycosphaerella graminicola, Septoria tritici) is a haploid fungus belonging to the class Dothideomycetes. It is the causal agent of septoria leaf blotch - one of the world's most significant diseases of wheat. Here we review the genomic and bioinformatic resources that have been generated for Z. tritici. These include the whole-genome reference assembly for isolate IPO323, genome resequencing of alternate isolates, mitochondrial genome sequences, transcriptome sequences and expression data, and annotations of gene structure and function. We also highlight important advances in our fundamental knowledge of genome evolution and its effects on adaptation and pathogenicity in Z. tritici that have been facilitated by these resources.

Dissecting the Molecular Interactions between Wheat and the Fungal Pathogen Zymoseptoria tritici

The Dothideomycete fungus Zymoseptoria tritici (previously known as Mycosphaerella graminicola and Septoria tritici) is the causative agent of Septoria tritici leaf blotch (STB) disease of wheat (Triticum aestivum L.). In Europe, STB is the most economically damaging disease of wheat, with an estimated ∼€1 billion per year in fungicide expenditure directed toward its control. Here, an overview of our current understanding of the molecular events that occur during Z. tritici infection of wheat leaves is presented. On the host side, this includes the contribution of (1) the pathogen-associated molecular pattern-triggered immunity (PTI) layer of plant defense, and (2) major Stb loci for resistance against Z. tritici. On the pathogen side of the interaction, we consolidate evidence from recent bioinformatic, transcriptomic and proteomic studies that begin to explain the contribution of Z. tritici effector proteins to the biphasic lifestyle of the fungus. This includes the discovery of chitin-binding proteins in the Z. tritici secretome, which contribute to evasion of immune surveillance by this pathogen, and the possible existence of 'necrotrophic' effectors from Z. tritici, which may actively stimulate host recognition in a manner similar to related necrotrophic fungal pathogens. We finish by speculating on how some of these recent fundamental discoveries might be harnessed to help improve resistance to STB in the world's second largest food crop.

Effector discovery in the fungal wheat pathogen Zymoseptoria tritici

Fungal plant pathogens, such as Zymoseptoria tritici (formerly known as Mycosphaerella graminicola), secrete repertoires of effectors to facilitate infection or trigger host defence mechanisms. The discovery and functional characterization of effectors provides valuable knowledge that can contribute to the design of new and effective disease management strategies. Here, we combined bioinformatics approaches with expression profiling during pathogenesis to identify candidate effectors of Z. tritici. In addition, a genetic approach was conducted to map quantitative trait loci (QTLs) carrying putative effectors, enabling the validation of both complementary strategies for effector discovery. In planta expression profiling revealed that candidate effectors were up-regulated in successive waves corresponding to consecutive stages of pathogenesis, contrary to candidates identified by QTL mapping that were, overall, expressed at low levels. Functional analyses of two top candidate effectors (SSP15 and SSP18) showed their dispensability for Z. tritici pathogenesis. These analyses reveal that generally adopted criteria, such as protein size, cysteine residues and expression during pathogenesis, may preclude an unbiased effector discovery. Indeed, genetic mapping of genomic regions involved in specificity render alternative effector candidates that do not match the aforementioned criteria, but should nevertheless be considered as promising new leads for effectors that are crucial for the Z. tritici-wheat pathosystem.

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.

Analysis of cytochrome b(5) reductase-mediated metabolism in the phytopathogenic fungus Zymoseptoria tritici reveals novel functionalities implicated in virulence

Septoria tritici blotch (STB) caused by the Ascomycete fungus Zymoseptoria tritici is one of the most economically damaging diseases of wheat worldwide. Z. tritici is currently a major target for agricultural fungicides, especially in temperate regions where it is most prevalent. Many fungicides target electron transfer enzymes because these are often important for cell function. Therefore characterisation of genes encoding such enzymes may be important for the development of novel disease intervention strategies. Microsomal cytochrome b5 reductases (CBRs) are an important family of electron transfer proteins which in eukaryotes are involved in the biosynthesis of fatty acids and complex lipids including sphingolipids and sterols. Unlike the model yeast Saccharomyces cerevisiae which possesses only one microsomal CBR, the fully sequenced genome of Z. tritici bears three possible microsomal CBRs. RNA sequencing analysis revealed that ZtCBR1 is the most highly expressed of these genes under all in vitro and in planta conditions tested, therefore ΔZtCBR1 mutant strains were generated through targeted gene disruption. These strains exhibited delayed disease symptoms on wheat leaves and severely limited asexual sporulation. ΔZtCBR1 strains also exhibited aberrant spore morphology and hyphal growth in vitro. These defects coincided with alterations in fatty acid, sphingolipid and sterol biosynthesis observed through GC-MS and HPLC analyses. Data is presented which suggests that Z. tritici may use ZtCBR1 as an additional electron donor for key steps in ergosterol biosynthesis, one of which is targeted by azole fungicides. Our study reports the first functional characterisation of CBR gene family members in a plant pathogenic filamentous fungus. This also represents the first direct observation of CBR functional ablation impacting upon fungal sterol biosynthesis.

Insight into the transcriptome of Arthrobotrys conoides using high throughput sequencing

Arthrobotrys conoides is a nematode-trapping fungus belonging to Orbiliales, Ascomycota group, and traps prey nematodes by means of adhesive network. Fungus has a potential to be used as a biocontrol agent against plant parasitic nematodes. In the present study, we characterized the transcriptome of A. conoides using high-throughput sequencing technology and characterized its virulence unigenes. Total 7,255 cDNA contigs with an average length of 425 bp were generated and 6184 (61.81%) transcripts were functionally annotated and characterized. Majority of unigenes were found analogous to the genes of plant pathogenic fungi. A total of 1749 transcripts were found to be orthologous with eukaryotic proteins of KOG database. Several carbohydrate active enzymes and peptidases were identified. We also analyzed classically and nonclassically secreted proteins and confirmed by BLASTP against fungal secretome database. A total of 916 contigs were analogous to 556 unique proteins of Pathogen Host Interaction (PHI) database. Further, we identified 91 unigenes homologous to the database of fungal virulence factor (DFVF). A total of 104 putative protein kinases coding transcripts were identified by BLASTP against KinBase database, which are major players in signaling pathways. This study provides a comprehensive look at the transcriptome of A. conoides and the identified unigenes might have a role in catching and killing prey nematodes by A. conoides.

An Interspecies Comparative Analysis of the Predicted Secretomes of the Necrotrophic Plant Pathogens Sclerotinia sclerotiorum and Botrytis cinerea

Phytopathogenic fungi form intimate associations with host plant species and cause disease. To be successful, fungal pathogens communicate with a susceptible host through the secretion of proteinaceous effectors, hydrolytic enzymes and metabolites. Sclerotinia sclerotiorum and Botrytis cinerea are economically important necrotrophic fungal pathogens that cause disease on numerous crop species. Here, a powerful bioinformatics pipeline was used to predict the refined S. sclerotiorum and B. cinerea secretomes, identifying 432 and 499 proteins respectively. Analyses focusing on S. sclerotiorum revealed that 16% of the secretome encoding genes resided in small, sequence heterogeneous, gene clusters that were distributed over 13 of the 16 predicted chromosomes. Functional analyses highlighted the importance of plant cell hydrolysis, oxidation-reduction processes and the redox state to the S. sclerotiorum and B. cinerea secretomes and potentially host infection. Only 8% of the predicted proteins were distinct between the two secretomes. In contrast to S. sclerotiorum, the B. cinerea secretome lacked CFEM- or LysM-containing proteins. The 115 fungal and oomycete genome comparison identified 30 proteins specific to S. sclerotiorum and B. cinerea, plus 11 proteins specific to S. sclerotiorum and 32 proteins specific to B. cinerea. Expressed sequence tag (EST) and proteomic analyses showed that 246 S. sclerotiorum secretome encoding genes had EST support, including 101 which were only expressed in vitro and 49 which were only expressed in planta, whilst 42 predicted proteins were experimentally proven to be secreted. These detailed in silico analyses of two important necrotrophic pathogens will permit informed choices to be made when candidate effector proteins are selected for function analyses in planta.

Previous bottlenecks and future solutions to dissecting the Zymoseptoria tritici-wheat host-pathogen interaction

Zymoseptoria tritici (previously Mycosphaerella graminicola, teleomorph, Septoria tritici, anamorph) causes Septoria tritici blotch, one of the most economically important diseases of wheat (Triticum aestivum). The host pathogenic interaction, as currently understood, is intriguing, and may distinguish Z. tritici from many of the current models for plant pathogenic fungi. Many important questions remain which require a deeper understanding including; the nature and biological significance of the characteristic long latent periods of symptomless plant infection; how/why the fungus then effectively transitions from this to cause disease and reproduce? Elements of this transition currently resemble a putative "hijack" on plant defence but how is Z. tritici able to do this without any form of plant cell penetration? This commentary provides a summary of the recent history of research into the host-pathogen interaction, whilst highlighting some of the challenges going forwards, which will be faced by improved technologies and a growing research community.

RNA-seq-Based Gene Annotation and Comparative Genomics of Four Fungal Grass Pathogens in the Genus Zymoseptoria Identify Novel Orphan Genes and Species-Specific Invasions of Transposable Elements

The fungal pathogen Zymoseptoria tritici (synonym Mycosphaerella graminicola) is a prominent pathogen of wheat. The reference genome of the isolate IPO323 is one of the best-assembled eukaryotic genomes and encodes more than 10,000 predicted genes. However, a large proportion of the previously annotated gene models are incomplete, with either no start or no stop codons. The availability of RNA-seq data allows better predictions of gene structure. We here used two different RNA-seq datasets, de novo transcriptome assemblies, homology-based comparisons, and trained ab initio gene callers to generate a new gene annotation of Z. tritici IPO323. The annotation pipeline was also applied to re-sequenced genomes of three closely related species of Z. tritici: Z. pseudotritici, Z. ardabiliae, and Z. brevis. Comparative analyses of the predicted gene models using the four Zymoseptoria species revealed sets of species-specific orphan genes enriched with putative pathogenicity-related genes encoding small secreted proteins that may play essential roles in virulence and host specificity. De novo repeat identification allowed us to show that few families of transposable elements are shared between Zymoseptoria species while we observe many species-specific invasions and expansions. The annotation data presented here provide a high-quality resource for future studies of Z. tritici and its sister species and provide detailed insight into gene and genome evolution of fungal plant pathogens.

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.

Comparative analysis of the in vitro and in planta secretomes from Mycosphaerella fijiensis isolates

Black Sigatoka, a devastating disease of bananas and plantains worldwide, is caused by the fungus Mycosphaerella fijiensis. Several banana cultivars such as 'Yangambi Km 5' and Calcutta IV, have been known to be resistant to the fungus, but the resistance has been broken in 'Yangambi Km 5' in Costa Rica. Since the resistance of this variety still persists in Mexico, the aim of this study was to compare the in vitro and in planta secretomes from two avirulent and virulent M. fijiensis isolates using proteomics and bioinformatics approaches. We aimed to identify differentially expressed proteins in fungal isolates that differ in pathogenicity and that might be responsible for breaking the resistance in 'Yangambi Km 5'. We were able to identify 90 protein spots in the secretomes of fungal isolates encoding 42 unique proteins and 35 differential spots between them. Proteins involved in carbohydrate transport and metabolism were more prevalent. Several proteases, pathogenicity-related, ROS detoxification and unknown proteins were also highly or specifically expressed by the virulent isolate in vitro or during in planta infection. An unknown protein representing a virulence factor candidate was also identified. These results demonstrated that the secretome reflects major differences between both M. fijiensis isolates.