SECRETOOL: integrated secretome analysis tool for fungi

Intrinsically Disordered Kiwellin Protein-Like Effectors Target Plant Chloroplasts and are Extensively Present in Rust Fungi

The effector proteins produced by plant pathogens are one of the essential components of host-pathogen interaction. Despite being important, most of the effector proteins remain unexplored due to the diversity in their primary sequence generated by the high selection pressure of the host immune system. However to maintain the primary function in the infection process, these effectors may tend to maintain their native protein fold to perform the corresponding biological function. In the present study, unannotated candidate secretory effector proteins of sixteen major plant fungal pathogens were analyzed to find the conserved known protein folds using homology, ab initio, and Alpha Fold/Rosetta Fold protein dimensional (3D) structure approaches. Several unannotated candidate effector proteins were found to match various known conserved protein families potentially involved in host defense manipulation in different plant pathogens. Surprisingly a large number of plant Kiwellin proteins fold like secretory proteins (> 100) were found in studied rust fungal pathogens. Many of them were predicted as potential effector proteins. Furthermore, template independent modelling using Alpha Fold/Rosetta Fold analysis and structural comparison of these candidates also predicted them to match with plant Kiwellin proteins. We also found plant Kiwellin matching proteins outside rusts including several non-pathogenic fungi suggesting the broad function of these proteins. One of the highest confidently modeled Kiwellin matching candidates effectors, Pstr_13960 (97.8%), from the Indian P. striiformis race Yr9 was characterized using overexpression, localization, and deletion studies in Nicotiana benthamiana. The Pstr_13960 suppressed the BAX-induced cell death and localized in the chloroplast. Furthermore, the expression of the Kiwellin matching region (Pst_13960_kiwi) alone suppressed the BAX-induced cell death in N. benthamiana despite the change of location to the cytoplasm and nucleus, suggesting the novel function of the Kiwellin core fold in rust fungi. Molecular docking showed that Pstr_13960 can interact with plant Chorismate mutases (CMs) using three loops conserved in plant and rust Kiwellins. Further analysis of Pstr_13960 showed to contain Intrinsically disordered regions (IDRs) in place of the N-terminal β1/β2 region found in plant Kiwellins suggesting the evolution of rust Kiwellins-like effectors (KLEs). Overall, this study reports the presence of a Kiwellin protein-like fold containing a novel effector protein family in rust fungi depicting a classical example of the evolution of effectors at the structure level as Kiwellin effectors show very low significant similarity to plant Kiwellin at the sequence level.

Deciphering the molecular components of the quorum sensing system in the fungus Ophiostoma piceae

IMPORTANCE

This manuscript presents a comprehensive study on the molecular mechanisms triggered by the quorum sensing (QS) molecule farnesol in the biotechnologically relevant fungus Ophiostoma piceae. We present for the first time, using a multiomics approach, an in-depth analysis of the QS response in a saprotroph fungus, detailing the molecular components involved in the response and their possible mechanisms of action. We think that these results are particularly relevant in the knowledge of the functioning of the QS in eukaryotes, as well as for the exploitation of these mechanisms.

Population and genome-wide association studies of Sclerotinia sclerotiorum isolates collected from diverse host plants throughout the United States

Introduction

Sclerotinia sclerotiorum is a necrotrophic fungal pathogen causing disease and economic loss on numerous crop plants. This fungus has a broad host range and can infect over 400 plant species, including important oilseed crops such as soybean, canola, and sunflower. S. sclerotiorum isolates vary in aggressiveness of lesion formation on plant tissues. However, the genetic basis for this variation remains to be determined. The aims of this study were to evaluate a diverse collection of S. sclerotiorum isolates collected from numerous hosts and U.S. states for aggressiveness of stem lesion formation on sunflower, to evaluate the population characteristics, and to identify loci associated with isolate aggressiveness using genome-wide association mapping.

Methods

A total of 219 S. sclerotiorum isolates were evaluated for stem lesion formation on two sunflower inbred lines and genotyped using genotyping-by-sequencing. DNA markers were used to assess population differentiation across hosts, regions, and climatic conditions and to perform a genome-wide association study of isolate aggressiveness.

Results and discussion

We observed a broad range of aggressiveness for lesion formation on sunflower stems, and only a moderate correlation between aggressiveness on the two lines. Population genetic evaluations revealed differentiation between populations from warmer climate regions compared to cooler regions. Finally, a genome-wide association study of isolate aggressiveness identified three loci significantly associated with aggressiveness on sunflower. Functional characterization of candidate genes at these loci will likely improve our understanding of the virulence strategies used by this pathogen to cause disease on a wide array of agriculturally important host plants.

Whole genome sequencing and analysis of multiple isolates of Ceratocystis destructans, the causal agent of Ceratocystis canker of almond in California

Ceratocystis canker caused by Ceratocystis destructans is a severe disease of almond, reducing the longevity and productivity of infected trees. Once the disease has established in an individual tree, there is no cure, and management efforts are often limited to removing the infected area of cankers. In this study, we present the genome assemblies of five C. destructans isolates isolated from symptomatic almond trees. The genomes were assembled into a genome size of 27.2 ± 0.9 Mbp with an average of 6924 ± 135 protein-coding genes and an average GC content of 48.8 ± 0.02%. We concentrated our efforts on identifying putative virulence factors of canker pathogens. Analysis of the secreted carbohydrate-active enzymes showed that the genomes harbored 83.4 ± 1.8 secreted CAZymes. The secreted CAZymes covered all the known categories of CAZymes. AntiSMASH revealed that the genomes had at least 7 biosynthetic gene clusters, with one of the non-ribosomal peptide synthases encoding dimethylcoprogen, a conserved virulence determinant of plant pathogenic ascomycetes. From the predicted proteome, we also annotated cytochrome P450 monooxygenases, and transporters, these are well-established virulence determinants of canker pathogens. Moreover, we managed to identify 57.4 ± 2.1 putative effector proteins. Gene Ontology (GO) annotation was applied to compare gene content with two closely related species C. fimbriata, and C. albifundus. This study provides the first genome assemblies for C. destructans, expanding genomic resources for an important almond canker pathogen. The acquired knowledge provides a foundation for further advanced studies, such as molecular interactions with the host, which is critical for breeding for resistance.

Unraveling the Secrets of a Double-Life Fungus by Genomics: Ophiocordyceps australis CCMB661 Displays Molecular Machinery for Both Parasitic and Endophytic Lifestyles

Ophiocordyceps australis (Ascomycota, Hypocreales, Ophiocordycipitaceae) is a classic entomopathogenic fungus that parasitizes ants (Hymenoptera, Ponerinae, Ponerini). Nonetheless, according to our results, this fungal species also exhibits a complete set of genes coding for plant cell wall degrading Carbohydrate-Active enZymes (CAZymes), enabling a full endophytic stage and, consequently, its dual ability to both parasitize insects and live inside plant tissue. The main objective of our study was the sequencing and full characterization of the genome of the fungal strain of O. australis (CCMB661) and its predicted secretome. The assembled genome had a total length of 30.31 Mb, N50 of 92.624 bp, GC content of 46.36%, and 8,043 protein-coding genes, 175 of which encoded CAZymes. In addition, the primary genes encoding proteins and critical enzymes during the infection process and those responsible for the host-pathogen interaction have been identified, including proteases (Pr1, Pr4), aminopeptidases, chitinases (Cht2), adhesins, lectins, lipases, and behavioral manipulators, such as enterotoxins, Protein Tyrosine Phosphatases (PTPs), and Glycoside Hydrolases (GHs). Our findings indicate that the presence of genes coding for Mad2 and GHs in O. australis may facilitate the infection process in plants, suggesting interkingdom colonization. Furthermore, our study elucidated the pathogenicity mechanisms for this Ophiocordyceps species, which still is scarcely studied.

First genome-scale insights into the virulence of the snow mold causal fungus Microdochium nivale

Pink snow mold, caused by a phytopathogenic and psychrotolerant fungus, Microdochium nivale, is a severe disease of winter cereals and grasses that predominantly occurs under snow cover or shortly after its melt. Snow mold has significantly progressed during the past decade, often reaching epiphytotic levels in northern countries and resulting in dramatic yield losses. In addition, M. nivale gradually adapts to a warmer climate, spreading to less snowy territories and causing different types of plant diseases throughout the growing period. Despite its great economic importance, M. nivale is poorly investigated; its genome has not been sequenced and its crucial virulence determinants have not been identified or even predicted. In our study, we applied a hybrid assembly based on Oxford Nanopore and Illumina reads to obtain the first genome sequence of M. nivale. 11,973 genes (including 11,789 protein-encoding genes) have been revealed in the genome assembly. To better understand the genetic potential of M. nivale and to obtain a convenient reference for transcriptomic studies on this species, the identified genes were annotated and split into hierarchical three-level functional categories. A file with functionally classified M. nivale genes is presented in our study for general use. M. nivale gene products that best meet the criteria for virulence factors have been identified. The genetic potential to synthesize human-dangerous mycotoxins (fumonisin, ochratoxin B, aflatoxin, and gliotoxin) has been revealed for M. nivale. The transcriptome analysis combined with the assays for extracellular enzymatic activities (conventional virulence factors of many phytopathogens) was carried out to assess the effect of host plant (rye) metabolites on the M. nivale phenotype. In addition to disclosing plant-metabolite-upregulated M. nivale functional gene groups (including those related to host plant protein destruction and amino acid metabolism, xenobiotic detoxication (including phytoalexins benzoxazinoids), cellulose destruction (cellulose monooxygenases), iron transport, etc.), the performed analysis pointed to a crucial role of host plant lipid destruction and fungal lipid metabolism modulation in plant-M. nivale interactions.

Comparative Genomics of Three Aspergillus Strains Reveals Insights into Endophytic Lifestyle and Endophyte-Induced Plant Growth Promotion

Aspergillus includes both plant pathogenic and beneficial fungi. Although endophytes beneficial to plants have high potential for plant growth promotion and improving stress tolerance, studies on endophytic lifestyles and endophyte-plant interactions are still limited. Here, three endophytes belonging to Aspergillus, AS31, AS33, and AS42, were isolated. They could successfully colonize rice roots and significantly improved rice growth. The genomes of strains AS31, AS33, and AS42 were sequenced and compared with other Aspergillus species covering both pathogens and endophytes. The genomes of AS31, AS33, and AS42 were 36.8, 34.8, and 35.3 Mb, respectively. The endophytic genomes had more genes encoding carbohydrate-active enzymes (CAZymes) and small secreted proteins (SSPs) and secondary metabolism gene clusters involved in indole metabolism than the pathogens. In addition, these endophytes were able to improve Pi (phosphorus) accumulation and transport in rice by inducing the expression of Pi transport genes in rice. Specifically, inoculation with endophytes significantly increased Pi contents in roots at the early stage, while the Pi contents in inoculated shoots were significantly increased at the late stage. Our results not only provide important insights into endophyte-plant interactions but also provide strain and genome resources, paving the way for the agricultural application of Aspergillus endophytes.

An Excretory Protein of Echinococcus multilocularis Inhibits Complement Classical Pathway Activation

Introduction

Alveolar echinococcosis is a lethal zoonosis caused by Echinococcus multilocularis (E.m) larvae. The mechanism by which E.m evades host immune attacks and ensures long-term survival remains unexplained. The complement system is a cascade of sequentially activated complement proteins that results in opsonization-related phagocytosis or membrane lysis of invading organisms. Excretory/secretory proteins (ESPs) of parasites are the main antigens that induce the immune response and play important roles in the long-term survival.

Methods

We investigated the possibility that E.m inhibits complement activation through ESPs and examined the potential related mechanism. A haemolysis assay was used to determine if and how in vitro culture medium of E.m containing ESPs can inhibit complement activation. Potential ESPs were annotated using bioinformatics methods, and one ESP was subsequently expressed as a recombinant protein with a eukaryotic expression system. The ability of this protein to inhibit complement activation was also tested by haemolysis assay.

Results

These assays showed that in vitro culture medium of E.m inhibited activation of the complement classical pathway. EmuJ_000439500 encodes a protein containing seven Sushi domains, which was the only potential E.m-derived complement inhibitor (Em-CI, UniProt: A0A068Y4F2) annotated among the 653 ESPs. Recombinant Em-CI also displayed the ability to inhibit activation of the complement classical pathway.

Discussion

The discovery of Em-CI sheds light on the mechanism by which E.m escapes killing by the complement system and provides potential targets for immunotherapy for parasitic diseases.

Gene Mining for Conserved, Non-Annotated Proteins of Podosphaera xanthii Identifies Novel Target Candidates for Controlling Powdery Mildews by Spray-Induced Gene Silencing

The powdery mildew fungus Podosphaera xanthii is one of the most important limiting factors for cucurbit production worldwide. Despite the significant efforts made by breeding and chemical companies, effective control of this pathogen remains elusive to growers. In this work, we examined the suitability of RNAi technology called spray-induced gene silencing (SIGS) for controlling cucurbit powdery mildew. Using leaf disc and cotyledon infiltration assays, we tested the efficacy of dsRNA applications to induce gene silencing in P. xanthii. Furthermore, to identify new target candidate genes, we analyzed sixty conserved and non-annotated proteins (CNAPs) deduced from the P. xanthii transcriptome in silico. Six proteins presumably involved in essential functions, specifically respiration (CNAP8878, CNAP9066, CNAP10905 and CNAP30520), glycosylation (CNAP1048) and efflux transport (CNAP948), were identified. Functional analysis of these CNAP coding genes by dsRNA-induced gene silencing resulted in strong silencing phenotypes with large reductions in fungal growth and disease symptoms. Due to their important contributions to fungal development, the CNAP1048, CNAP10905 and CNAP30520 genes were selected as targets to conduct SIGS assays under plant growth chamber conditions. The spray application of these dsRNAs induced high levels of disease control, supporting that SIGS could be a sustainable approach to combat powdery mildew diseases.

Comparative genome analyses suggest a hemibiotrophic lifestyle and virulence differences for the beech bark disease fungal pathogens Neonectria faginata and Neonectria coccinea

Neonectria faginata and Neonectria coccinea are the causal agents of the insect-fungus disease complex known as beech bark disease (BBD), known to cause mortality in beech forest stands in North America and Europe. These fungal species have been the focus of extensive ecological and disease management studies, yet less progress has been made toward generating genomic resources for both micro- and macro-evolutionary studies. Here, we report a 42.1 and 42.7 mb highly contiguous genome assemblies of N. faginata and N. coccinea, respectively, obtained using Illumina technology. These species share similar gene number counts (12,941 and 12,991) and percentages of predicted genes with assigned functional categories (64 and 65%). Approximately 32% of the predicted proteomes of both species are homologous to proteins involved in pathogenicity, yet N. coccinea shows a higher number of predicted mitogen-activated protein kinase genes, virulence determinants possibly contributing to differences in disease severity between N. faginata and N. coccinea. A wide range of genes encoding for carbohydrate-active enzymes capable of degradation of complex plant polysaccharides and a small number of predicted secretory effector proteins, secondary metabolite biosynthesis clusters and cytochrome oxidase P450 genes were also found. This arsenal of enzymes and effectors correlates with, and reflects, the hemibiotrophic lifestyle of these two fungal pathogens. Phylogenomic analysis and timetree estimations indicated that the N. faginata and N. coccinea species divergence may have occurred at ∼4.1 million years ago. Differences were also observed in the annotated mitochondrial genomes as they were found to be 81.7 kb (N. faginata) and 43.2 kb (N. coccinea) in size. The mitochondrial DNA expansion observed in N. faginata is attributed to the invasion of introns into diverse intra- and intergenic locations. These first draft genomes of N. faginata and N. coccinea serve as valuable tools to increase our understanding of basic genetics, evolutionary mechanisms and molecular physiology of these two nectriaceous plant pathogenic species.

Secretome Analysis of Arabidopsis-Trichoderma atroviride Interaction Unveils New Roles for the Plant Glutamate:Glyoxylate Aminotransferase GGAT1 in Plant Growth Induced by the Fungus and Resistance against Botrytis cinerea

The establishment of plant-fungus mutualistic interaction requires bidirectional molecular crosstalk. Therefore, the analysis of the interacting organisms secretomes would help to understand how such relationships are established. Here, a gel-free shotgun proteomics approach was used to identify the secreted proteins of the plant Arabidopsis thaliana and the mutualistic fungus Trichoderma atroviride during their interaction. A total of 126 proteins of Arabidopsis and 1027 of T. atroviride were identified. Among them, 118 and 780 were differentially modulated, respectively. Bioinformatic analysis unveiled that both organisms' secretomes were enriched with enzymes. In T. atroviride, glycosidases, aspartic endopeptidases, and dehydrogenases increased in response to Arabidopsis. Additionally, amidases, protein-serine/threonine kinases, and hydro-lyases showed decreased levels. Furthermore, peroxidases, cysteine endopeptidases, and enzymes related to the catabolism of secondary metabolites increased in the plant secretome. In contrast, pathogenesis-related proteins and protease inhibitors decreased in response to the fungus. Notably, the glutamate:glyoxylate aminotransferase GGAT1 was secreted by Arabidopsis during its interaction with T. atroviride. Our study showed that GGAT1 is partially required for plant growth stimulation and on the induction of the plant systemic resistance by T. atroviride. Additionally, GGAT1 seems to participate in the negative regulation of the plant systemic resistance against B. cinerea through a mechanism involving H2O2 production.

A Multiomic Approach to Understand How Pleurotus eryngii Transforms Non-Woody Lignocellulosic Material

Pleurotus eryngii is a grassland-inhabiting fungus of biotechnological interest due to its ability to colonize non-woody lignocellulosic material. Genomic, transcriptomic, exoproteomic, and metabolomic analyses were combined to explain the enzymatic aspects underlaying wheat–straw transformation. Up-regulated and constitutive glycoside–hydrolases, polysaccharide–lyases, and carbohydrate–esterases active on polysaccharides, laccases active on lignin, and a surprisingly high amount of constitutive/inducible aryl–alcohol oxidases (AAOs) constituted the suite of extracellular enzymes at early fungal growth. Higher enzyme diversity and abundance characterized the longer-term growth, with an array of oxidoreductases involved in depolymerization of both cellulose and lignin, which were often up-regulated since initial growth. These oxidative enzymes included lytic polysaccharide monooxygenases (LPMOs) acting on crystalline polysaccharides, cellobiose dehydrogenase involved in LPMO activation, and ligninolytic peroxidases (mainly manganese-oxidizing peroxidases), together with highly abundant H₂O₂-producing AAOs. Interestingly, some of the most relevant enzymes acting on polysaccharides were appended to a cellulose-binding module. This is potentially related to the non-woody habitat of P. eryngii (in contrast to the wood habitat of many basidiomycetes). Additionally, insights into the intracellular catabolism of aromatic compounds, which is a neglected area of study in lignin degradation by basidiomycetes, were also provided. The multiomic approach reveals that although non-woody decay does not result in dramatic modifications, as revealed by detailed 2D-NMR and other analyses, it implies activation of the complete set of hydrolytic and oxidative enzymes characterizing lignocellulose-decaying basidiomycetes.

Interaction between the cell walls of microalgal host and fungal carbohydrate-activate enzymes is essential for the pathogenic parasitism process

Fungi can parasitize microalgae, exerting profound impacts on both the aquatic ecosystems and microalgal mass cultures. In this study, the unicellular green alga Haematococcus pluvialis and the blastocladialean fungus Paraphysoderma sedebokerense were used as a model system to address the mechanisms underlying the fungal parasitism on the algal host. High-throughput metabolic assay indicated that P. sedebokerense can utilize several carbon sources with a preference for mannose, glucose and their oligosaccharides, which was compatible with the profile of the host algal cell walls enriched with glucan and mannan. The results of dual transcriptomics analysis suggested that P. sedebokerense can upregulate a large number of putative carbohydrate-activate enzymes (CAZymes) encoding genes, including those coding for the endo-1,4-β-glucanase and endo-1,4-β-mannanase during the infection process. The cell walls of H. pluvialis can be decomposed by both P. sedebokerense and commercial CAZymes (e.g. cellulase and endo-1,4-β-mannanase) to produce mannooligomers, while several putative parasitism-related genes of P. sedebokerense can be in turn upregulated by mannooligomers. In addition, the parasitism can be blocked by interfering the selected CAZymes including glucanase, mannanase and lysozyme with the specific inhibitors, which provided a framework for screening suitable compounds for pathogen mitigation in algal mass culture.

Near-Chromosome-Level Genome Assembly of the Dark Septate Endophyte Laburnicola rhizohalophila: A Model for Investigating Root-Fungus Symbiosis

The novel DSE Laburnicola rhizohalophila (Pleosporales, Ascomycota) is frequently found in the halophytic seepweed (Suaeda salsa). In this article, we report a near-chromosome-level hybrid assembly of this fungus using a combination of short-read Illumina data to polish assemblies generated from long-read Nanopore data. The reference genome for L. rhizohalophila was assembled into 26 scaffolds with a total length of 64.0 Mb and a N50 length of 3.15 Mb. Of them, 17 scaffolds approached the length of intact chromosomes, and 5 had telomeres at one end only. A total of 10,891 gene models were predicted. Intriguingly, 27.5 Mb of repeat sequences that accounted for 42.97% of the genome was identified, and long terminal repeat retrotransposons were the most frequent known transposable elements, indicating that transposable element proliferation contributes to its increased genome size. BUSCO analyses using the Fungi_odb10 data set showed that 95.0% of genes were complete. In addition, 292 carbohydrate active enzymes, 33 secondary metabolite clusters, and 84 putative effectors were identified in silico. The resulting high-quality assembly and genome features are not only an important resource for further research on understanding the mechanism of root-fungi symbiotic interactions but will also contribute to comparative analyses of genome biology and evolution within Pleosporalean species.

Secretomic insight into the biomass hydrolysis potential of the phytopathogenic fungus Chrysoporthe cubensis

The phytopathogenic fungus Chrysoporthe cubensis has a great capacity to produce highly efficient enzymes for the hydrolysis of lignocellulosic biomass. The bioinfosecretome of C. cubensis was identified by computational predictions of secreted proteins combined with protein analysis using 1D-LC-MS/MS. The in silico secretome predicted 562 putative genes capable of encoding secreted proteins, including 273 CAZymes. Proteomics analysis confirmed the existence of 313 proteins, including 137 CAZymes classified as Glycosyl Hydrolases (GH), Polysaccharide Lyases (PL), Carbohydrate Esterases (CE) and Auxiliary Activities enzymes (AA), which indicates the presence of classical and oxidative cellulolytic mechanisms. The enzymes diversity in the extract shows fungal versatility to act in complex biomasses. This study provides an insight into the lignocellulose-degradation mechanisms by C. cubensis and allows the identification of the enzymes that are potentially useful in improving industrial process of bioconversion of lignocellulose. SIGNIFICANCE: Chrysoporthe cubensis is an important deadly canker pathogen of commercially cultivated Eucalyptus species. The effective depolymerisation of the recalcitrant plant cell wall performed by this fungus is closely related to its high potential of lignocellulolytic enzymes secretion. Since the degradation of biomass occurs in nature almost exclusively by enzyme secretion systems, it is reasonable to suggest that the identification of C. cubensis lignocellulolytic enzymes is relevant in contributing to new sustainable alternatives for industrial solutions. As far as we know, this work is the first accurate proteomic evaluation of the enzymes secreted by this species of fungus. The integration of the gel-based proteomic approach, the bioinformatic prediction of the secretome and the analyses of enzymatic activity are powerful tools in the evaluation of biotechnological potential of C. cubensis in producing carbohydrate-active enzymes. In addition, analysis of the C. cubensis secretome grown in wheat bran draws attention to this plant pathogen and its extracellular enzymatic machinery, especially regarding the identification of promising new enzymes for industrial applications. The results from this work allowed for explanation and reinforce previous research that revealed C. cubensis as a strong candidate to produce enzymes to hydrolyse sugarcane bagasse and similar substrates.

Differential Expression of Fungal Genes Determines the Lifestyle of Plectosphaerella Strains During Arabidopsis thaliana Colonization

The fungal genus Plectosphaerella comprises species and strains with different lifestyles on plants, such as P. cucumerina, which has served as model for the characterization of Arabidopsis thaliana basal and nonhost resistance to necrotrophic fungi. We have sequenced, annotated, and compared the genomes and transcriptomes of three Plectosphaerella strains with different lifestyles on A. thaliana, namely, PcBMM, a natural pathogen of wild-type plants (Col-0), Pc2127, a nonpathogenic strain on Col-0 but pathogenic on the immunocompromised cyp79B2 cyp79B3 mutant, and P0831, which was isolated from a natural population of A. thaliana and is shown here to be nonpathogenic and to grow epiphytically on Col-0 and cyp79B2 cyp79B3 plants. The genomes of these Plectosphaerella strains are very similar and do not differ in the number of genes with pathogenesis-related functions, with the exception of secreted carbohydrate-active enzymes (CAZymes), which are up to five times more abundant in the pathogenic strain PcBMM. Analysis of the fungal transcriptomes in inoculated Col-0 and cyp79B2 cyp79B3 plants at initial colonization stages confirm the key role of secreted CAZymes in the necrotrophic interaction, since PcBMM expresses more genes encoding secreted CAZymes than Pc2127 and P0831. We also show that P0831 epiphytic growth on A. thaliana involves the transcription of specific repertoires of fungal genes, which might be necessary for epiphytic growth adaptation. Overall, these results suggest that in-planta expression of specific sets of fungal genes at early stages of colonization determine the diverse lifestyles and pathogenicity of Plectosphaerella strains.

Secretomes of medically important fungi reflect morphological and phylogenetic diversity

Secretome represents a main target for understanding the mechanisms of fungal adaptation. In the present study, we focus on the secretomes of fungi associated with infections in humans and other mammals in order to explore relationships between the diverse morphological and phylogenetic groups. Almost all the mammalian pathogenic fungi analyzed have secretome sizes smaller than 1000 proteins and, secreted proteins comprise between 5% and 10% of the total proteome. As expected, the correlation pattern between the secretome size and the total proteome was similar to that described in previous secretome studies of fungi. With regard to the morphological groups, minimum secretome sizes of less than 250 secreted proteins and low values for the fraction of secreted proteins are shown in mammalian pathogenic fungi with reduced proteomes such as microsporidia, atypical fungi and some species of yeasts and yeast-like fungi (Malassezia). On the other hand, filamentous fungi have significantly more secreted proteins and the highest numbers are present in species of filamentous fungi that also are plant or insect pathogens (Fusarium verticilloides, Fusarium oxysporum and Basidiobolus meristosporus). With respect to phylogeny, there are also variations in secretome size across fungal subphyla: Microsporidia, Taphrinomycotina, Ustilagomycotina and Saccharomycotina contain small secretomes; whereas larger secretomes are found in Agaricomycotina, Pezizomycotina, Mucoromycotina and Entomophthoromycotina. Finally, principal component analysis (PCA) was conducted on the complete secretomes. The PCA results revealed that, in general, secretomes of fungi belonging to the same morphological group or subphyla cluster together. In conclusion, our results point out that in medically important fungi there is a relationship between the secretome and the morphological group or phylogenetic classification.

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.

Draft Genome Sequence of the Oomycete Globisporangium splendens Strain rgcb-1

Globisporangium splendens (formerly Pythium splendens) is an oomycete pathogen of many economically important vegetable crops. Here, we present the first draft genome of P. splendens, which comprises 197 scaffolds with a total length of 53.3 Mb and 17,350 predicted protein-coding genes.

Globisporangium splendens (formerly Pythium splendens) is an oomycete pathogen of many economically important vegetable crops. Here, we present the first draft genome of P. splendens, which comprises 197 scaffolds with a total length of 53.3 Mb and 17,350 predicted protein-coding genes.

Genomic Features of Cladobotryum dendroides, Which Causes Cobweb Disease in Edible Mushrooms, and Identification of Genes Related to Pathogenicity and Mycoparasitism

Cladobotryum dendroides, which causes cobweb disease in edible mushrooms, is one of the major fungal pathogens. Our previous studies focused on the genetic and morphological characterization of this fungus, as well as its pathogenicity and the identification of appropriate fungicides. However, little is known about the genome characters, pathogenic genes, and molecular pathogenic mechanisms of C. dendroides. Herein, we reported a high-quality de novo genomic sequence of C. dendroides and compared it with closely-related fungi. The assembled C. dendroides genome was 36.69 Mb, consisting of eight contigs, with an N50 of 4.76 Mb. This genome was similar in size to that of C. protrusum, and shared highly conserved syntenic blocks and a few inversions with C. protrusum. Phylogenetic analysis revealed that, within the Hypocreaceae, Cladobotryum was closer to Mycogone than to Trichoderma, which is consistent with phenotypic evidence. A significant number of the predicted expanded gene families were strongly associated with pathogenicity, virulence, and adaptation. Our findings will be instrumental for the understanding of fungi-fungi interactions, and for exploring efficient management strategies to control cobweb disease.

The haustorial transcriptome of the cucurbit pathogen Podosphaera xanthii reveals new insights into the biotrophy and pathogenesis of powdery mildew fungi

Background

Podosphaera xanthii is the main causal agent of powdery mildew disease in cucurbits and is responsible for important yield losses in these crops worldwide. Powdery mildew fungi are obligate biotrophs. In these parasites, biotrophy is determined by the presence of haustoria, which are specialized structures of parasitism developed by these fungi for the acquisition of nutrients and the delivery of effectors. Detailed molecular studies of powdery mildew haustoria are scarce due mainly to difficulties in their isolation. Therefore, their analysis is considered an important challenge for powdery mildew research. The aim of this work was to gain insights into powdery mildew biology by analysing the haustorial transcriptome of P. xanthii.

Results

Prior to RNA isolation and massive-scale mRNA sequencing, a flow cytometric approach was developed to isolate P. xanthii haustoria free of visible contaminants. Next, several commercial kits were used to isolate total RNA and to construct the cDNA and Illumina libraries that were finally sequenced by the Illumina NextSeq system. Using this approach, the maximum amount of information from low-quality RNA that could be obtained was used to accomplish the de novo assembly of the P. xanthii haustorial transcriptome. The subsequent analysis of this transcriptome and comparison with the epiphytic transcriptome allowed us to identify the importance of several biological processes for haustorial cells such as protection against reactive oxygen species, the acquisition of different nutrients and genetic regulation mediated by non-coding RNAs. In addition, we could also identify several secreted proteins expressed exclusively in haustoria such as cell adhesion proteins that have not been related to powdery mildew biology to date.

Conclusions

This work provides a novel approach to study the molecular aspects of powdery mildew haustoria. In addition, the results of this study have also allowed us to identify certain previously unknown processes and proteins involved in the biology of powdery mildews that could be essential for their biotrophy and pathogenesis.

Electronic supplementary material

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

Genome analysis of the ubiquitous boxwood pathogen Pseudonectria foliicola

Boxwood (Buxus spp.) are broad-leaved, evergreen landscape plants valued for their longevity and ornamental qualities. Volutella leaf and stem blight, caused by the ascomycete fungi Pseudonectria foliicola and P. buxi, is one of the major diseases affecting the health and ornamental qualities of boxwood. Although this disease is less severe than boxwood blight caused by Calonectria pseudonaviculata and C. henricotiae, its widespread occurrence and disfiguring symptoms have caused substantial economic losses to the ornamental industry. In this study, we sequenced the genome of P. foliicola isolate ATCC13545 using Illumina technology and compared it to other publicly available fungal pathogen genomes to better understand the biology of this organism. A de novo assembly estimated the genome size of P. foliicola at 28.7 Mb (425 contigs; N50 = 184,987 bp; avg. coverage 188×), with just 9,272 protein-coding genes. To our knowledge, P. foliicola has the smallest known genome within the Nectriaceae. Consistent with the small size of the genome, the secretome, CAzyme and secondary metabolite profiles of this fungus are reduced relative to two other surveyed Nectriaceae fungal genomes: Dactylonectria macrodidyma JAC15-245 and Fusarium graminearum Ph-1. Interestingly, a large cohort of genes associated with reduced virulence and loss of pathogenicity was identified from the P. foliicola dataset. These data are consistent with the latest observations by plant pathologists that P. buxi and most likely P. foliicola, are opportunistic, latent pathogens that prey upon weak and stressed boxwood plants.

Genome-wide analysis of Excretory/Secretory proteins in Trypanosoma brucei brucei: Insights into functional characteristics and identification of potential targets by immunoinformatics approach

Trypanosoma brucei brucei (T.b.brucei) is an extra-cellular parasite that causes Animal African Trypanosomiasis (AAT) disease in animals. Till day, this disease is more difficult to treat and control due to lack of efficient vaccines and early diagnosis of the parasite infection. T.b.brucei Excretory/Secretory (ES) proteins were involved in pathogenesis and key for understanding the host-parasite interactions. Functions of T.b.brucei's ES proteins were poorly investigated and experimental identification is expensive and time-consuming. Bioinformatics approaches are cost-effective by facilitating the experimental analysis of potential drug targets for parasitic diseases. Here we applied several bioinformatics tools to predict and functionalize the annotation of 1104 ES proteins and immunoinformatics approaches carried out to predict and evaluate the epitopes in T.b.brucei. Secretory information, functional annotations and potential epitopes of each ES proteins were available at http://tbb.insilico.in. This study provides functional information of T.b.brucei for experimental studies to identify potential targets for diagnosis and therapeutics development.

Comparative systems analysis of the secretome of the opportunistic pathogen Aspergillus fumigatus and other Aspergillus species

Aspergillus fumigatus and multiple other Aspergillus species cause a wide range of lung infections, collectively termed aspergillosis. Aspergilli are ubiquitous in environment with healthy immune systems routinely eliminating inhaled conidia, however, Aspergilli can become an opportunistic pathogen in immune-compromised patients. The aspergillosis mortality rate and emergence of drug-resistance reveals an urgent need to identify novel targets. Secreted and cell membrane proteins play a critical role in fungal-host interactions and pathogenesis. Using a computational pipeline integrating data from high-throughput experiments and bioinformatic predictions, we have identified secreted and cell membrane proteins in ten Aspergillus species known to cause aspergillosis. Small secreted and effector-like proteins similar to agents of fungal-plant pathogenesis were also identified within each secretome. A comparison with humans revealed that at least 70% of Aspergillus secretomes have no sequence similarity with the human proteome. An analysis of antigenic qualities of Aspergillus proteins revealed that the secretome is significantly more antigenic than cell membrane proteins or the complete proteome. Finally, overlaying an expression dataset, four A. fumigatus proteins upregulated during infection and with available structures, were found to be structurally similar to known drug target proteins in other organisms, and were able to dock in silico with the respective drug.

VerSeDa: vertebrate secretome database

Based on the current tools, de novo secretome (full set of proteins secreted by an organism) prediction is a time consuming bioinformatic task that requires a multifactorial analysis in order to obtain reliable in silico predictions. Hence, to accelerate this process and offer researchers a reliable repository where secretome information can be obtained for vertebrates and model organisms, we have developed VerSeDa (Vertebrate Secretome Database). This freely available database stores information about proteins that are predicted to be secreted through the classical and non-classical mechanisms, for the wide range of vertebrate species deposited at the NCBI, UCSC and ENSEMBL sites. To our knowledge, VerSeDa is the only state-of-the-art database designed to store secretome data from multiple vertebrate genomes, thus, saving an important amount of time spent in the prediction of protein features that can be retrieved from this repository directly.

Database URL: VerSeDa is freely available at http://genomics.cicbiogune.es/VerSeDa/index.php

Comparative genomics of Coniophora olivacea reveals different patterns of genome expansion in Boletales

Background

Coniophora olivacea is a basidiomycete fungus belonging to the order Boletales that produces brown-rot decay on dead wood of conifers. The Boletales order comprises a diverse group of species including saprotrophs and ectomycorrhizal fungi that show important differences in genome size.

Results

In this study we report the 39.07-megabase (Mb) draft genome assembly and annotation of C. olivacea. A total of 14,928 genes were annotated, including 470 putatively secreted proteins enriched in functions involved in lignocellulose degradation. Using similarity clustering and protein structure prediction we identified a new family of 10 putative lytic polysaccharide monooxygenase genes. This family is conserved in basidiomycota and lacks of previous functional annotation. Further analyses showed that C. olivacea has a low repetitive genome, with 2.91% of repeats and a restrained content of transposable elements (TEs). The annotation of TEs in four related Boletales yielded important differences in repeat content, ranging from 3.94 to 41.17% of the genome size. The distribution of insertion ages of LTR-retrotransposons showed that differential expansions of these repetitive elements have shaped the genome architecture of Boletales over the last 60 million years.

Conclusions

Coniophora olivacea has a small, compact genome that shows macrosynteny with Coniophora puteana. The functional annotation revealed the enzymatic signature of a canonical brown-rot. The annotation and comparative genomics of transposable elements uncovered their particular contraction in the Coniophora genera, highlighting their role in the differential genome expansions found in Boletales species.

Electronic supplementary material

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

Computational Prediction of Effector Proteins in Fungi: Opportunities and Challenges

Effector proteins are mostly secretory proteins that stimulate plant infection by manipulating the host response. Identifying fungal effector proteins and understanding their function is of great importance in efforts to curb losses to plant diseases. Recent advances in high-throughput sequencing technologies have facilitated the availability of several fungal genomes and 1000s of transcriptomes. As a result, the growing amount of genomic information has provided great opportunities to identify putative effector proteins in different fungal species. There is little consensus over the annotation and functionality of effector proteins, and mostly small secretory proteins are considered as effector proteins, a concept that tends to overestimate the number of proteins involved in a plant-pathogen interaction. With the characterization of Avr genes, criteria for computational prediction of effector proteins are becoming more efficient. There are 100s of tools available for the identification of conserved motifs, signature sequences and structural features in the proteins. Many pipelines and online servers, which combine several tools, are made available to perform genome-wide identification of effector proteins. In this review, available tools and pipelines, their strength and limitations for effective identification of fungal effector proteins are discussed. We also present an exhaustive list of classically secreted proteins along with their key conserved motifs found in 12 common plant pathogens (11 fungi and one oomycete) through an analytical pipeline.

Comparative Transcriptomic Analysis of Virulence Factors in Leptosphaeria maculans during Compatible and Incompatible Interactions with Canola

Leptosphaeria maculans is a hemibiotrophic fungus that causes blackleg of canola (Brassica napus), one of the most devastating diseases of this crop. In the present study, transcriptome profiling of L. maculans was performed in an effort to understand and define the pathogenicity genes that govern both the biotrophic and the necrotrophic phase of the fungus, as well as those that separate a compatible from an incompatible interaction. For this purpose, comparative RNA-seq analyses were performed on L. maculans isolate D5 at four different time points following inoculation on susceptible cultivar Topas-DH16516 or resistant introgression line Topas-Rlm2. Analysis of 1.6 billion Illumina reads readily identified differentially expressed genes that were over represented by candidate secretory effector proteins, CAZymes, and other pathogenicity genes. Comparisons between the compatible and incompatible interactions led to the identification of 28 effector proteins whose chronology and level of expression suggested a role in the establishment and maintenance of biotrophy with the plant. These included all known Avr genes of isolate D5 along with eight newly characterized effectors. In addition, another 15 effector proteins were found to be exclusively expressed during the necrotrophic phase of the fungus, which supports the concept that L. maculans has a separate and distinct arsenal contributing to each phase. As for CAZymes, they were often highly expressed at 3 dpi but with no difference in expression between the compatible and incompatible interactions, indicating that other factors were necessary to determine the outcome of the interaction. However, their significantly higher expression at 11 dpi in the compatible interaction confirmed that they contributed to the necrotrophic phase of the fungus. A notable exception was LysM genes whose high expression was singularly observed on the susceptible host at 7 dpi. In the case of TFs, their higher expression at 7 and 11 dpi on susceptible Topas support an important role in regulating the genes involved in the different pathogenic phases of L. maculans. In conclusion, comparison of the transcriptome of L. maculans during compatible and incompatible interactions has led to the identification of key pathogenicity genes that regulate not only the fate of the interaction but also lifestyle transitions of the fungus.

Genome-wide analysis of excretory/secretory proteins in Echinococcus multilocularis: insights into functional characteristics of the tapeworm secretome

Background

The cestode Echinococcus multilocularis is the causative agent of human alveolar echinococcosis (AE). However, this life-threatening disease is still difficult to treat and control, due to the lack of efficient drugs and vaccines. Excretory/secretory (ES) proteins are crucial for parasite survival and represent potential preferred targets for novel intervention strategies. However, the ES protein features in this parasite have been poorly investigated until now. The current study was carried out to identify and characterise a repertoire of ES proteins in E. multilocularis at the genome-wide level.

Methods

Here we predicted and functionally annotated the classical and non-classical ES proteins, and comprehensively compared the features and evolution of ES and non-ES proteins in E. multilocularis genome using an integration of bioinformatics tools. The intervention target and antigen potentials as well as the transcription information were also investigated.

Results

Computational analysis of the E. multilocularis proteins identified 673 putative ES proteins (6.4 %), of which 617 (91.68 %) could be supported by transcription analyses. The predicted ES proteins in E. multilocularis were mostly represented by molecular functions of protease inhibitors, proteases, glycoside hydrolases, immunoglobulin-like folds and growth factors. Analysis of the ratio between synonymous and non-synonymous substitution rates (dN/dS) revealed no significant difference of the evolution selection pressure on ES and non-ES protein coding genes. ES proteins showed higher antigenic density measured by AAR values as compared with the transmembrane proteins but had no significant difference of that feature with intracellular proteins. Additionally, 383 possible ideal drug targets were identified in ES proteins, of which four proteins have corresponding known drugs.

Conclusions

The present study is the first to identify a repertoire of predicted ES proteins at the genome-wide level in E. multilocularis. The comprehensive analysis provides some novel understanding of the parasite ES protein features and a valuable resource of potential targets for future experimental studies to develop new intervention tools against this parasite.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1282-7) contains supplementary material, which is available to authorized users.

IMA Genome-F 4: Draft genome sequences of Chrysoporthe austroafricana, Diplodia scrobiculata, Fusarium nygamai, Leptographium lundbergii, Limonomyces culmigenus, Stagonosporopsis tanaceti, and Thielaviopsis punctulata

The genomes of Chrysoporthe austroafricana, Diplodia scrobiculata, Fusarium nygami, Leptographium lundbergii, Limonomyces culmigenus, Stagonosporopsis tanaceti, and Thielaviopsis punctulata are presented in this genome announcement. These seven genomes are from endophytes, plant pathogens and economically important fungal species. The genome sizes range from 26.6 Mb in the case of Leptographium lundbergii to 44 Mb for Chrysoporthe austroafricana. The availability of these genome data will provide opportunities to resolve longstanding questions regarding the taxonomy of species in these genera, and may contribute to our understanding of the lifestyles through comparative studies with closely related organisms.

Destructuring plant biomass: focus on fungal and extremophilic cell wall hydrolases

The use of plant biomass as feedstock for biomaterial and biofuel production is relevant in the current bio-based economy scenario of valorizing renewable resources. Fungi, which degrade complex and recalcitrant plant polymers, secrete different enzymes that hydrolyze plant cell wall polysaccharides. The present review discusses the current research trends on fungal, as well as extremophilic cell wall hydrolases that can withstand extreme physico-chemical conditions required in efficient industrial processes. Secretomes of fungi from the phyla Ascomycota, Basidiomycota, Zygomycota and Neocallimastigomycota are presented along with metabolic cues (nutrient sensing, coordination of carbon and nitrogen metabolism) affecting their composition. We conclude the review by suggesting further research avenues focused on the one hand on a comprehensive analysis of the physiology and epigenetics underlying cell wall degrading enzyme production in fungi and on the other hand on the analysis of proteins with unknown function and metagenomics of extremophilic consortia. The current advances in consolidated bioprocessing, altered secretory pathways and creation of designer plants are also examined. Furthermore, recent developments in enhancing the activity, stability and reusability of enzymes based on synergistic, proximity and entropic effects, fusion enzymes, structure-guided recombination between homologous enzymes and magnetic enzymes are considered with a view to improving saccharification.

Secretome of the biocontrol agent metarhizium anisopliae induced by the cuticle of the cotton pest Dysdercus peruvianus reveals new insights into infection

Metarhizium anisopliae is an entomopathogenic fungus that has evolved specialized strategies to infect insect hosts. Here we analyzed secreted proteins related to Dysdercus peruvianus infection. Using shotgun proteomics, abundance changes in 71 proteins were identified after exposure to host cuticle. Among these proteins were classical fungal effectors secreted by pathogens to degrade physical barriers and alter host physiology. These include lipolytic enzymes, Pr1A, B, C, I, and J proteases, ROS-related proteins, oxidorreductases, and signaling proteins. Protein interaction networks were generated postulating interesting candidates for further studies, including Pr1C, based on possible functional interactions. On the basis of these results, we propose that M. anisopliae is degrading host components and actively secreting proteins to manage the physiology of the host. Interestingly, the secretion of these factors occurs in the absence of a host response. The findings presented here are an important step in understanding the host–pathogen interaction and developing more efficient biocontrol of D. peruvianus by M. anisopliae.