Comparative genomics and transcriptomics analyses reveal divergent lifestyle features of nematode endoparasitic fungus Hirsutella minnesotensis

Aspergillus fumigatus secondary metabolite pyripyropene is important for the dual biofilm formation with Pseudomonas aeruginosa

The human pathogenic fungus Aspergillus fumigatus establishes dual biofilm interactions in the lungs with the pathogenic bacterium Pseudomonas aeruginosa. Screening of 21 A. fumigatus null mutants revealed seven mutants (two G protein-coupled receptors, three mitogen-activated protein kinase receptors, a Gα protein, and one histidine kinase receptor) with reduced biofilm formation, specifically in the presence of P. aeruginosa. Transcriptional profiling and metabolomics analysis of secondary metabolites produced by one of these mutants, ΔgpaB (gpaB encodes a Gα protein), showed GpaB controls the production of several important metabolites for the dual biofilm interaction, including pyripyropene A, a potent inhibitor of mammalian acyl-CoA cholesterol acyltransferase. Deletion of pyr2, encoding a non-reducing polyketide synthase essential for pyripyropene biosynthesis, showed reduced A. fumigatus Δpyr2-P. aeruginosa biofilm growth, altered macrophage responses, and attenuated mouse virulence in a chemotherapeutic murine model. We identified pyripyropene as a novel player in the ecology and pathogenic interactions of this important human fungal pathogen.IMPORTANCEAspergillus fumigatus and Pseudomonas aeruginosa are two important human pathogens. Both organisms establish biofilm interactions in patients affected with chronic lung pulmonary infections, such as cystic fibrosis (CF) and chronic obstructive pulmonary disease. Colonization with A. fumigatus is associated with an increased risk of P. aeruginosa colonization in CF patients, and disease prognosis is poor when both pathogens are present. Here, we identified A. fumigatus genetic determinants important for the establishment of in vitro dual A. fumigatus-P. aeruginosa biofilm interactions. Among them, an A. fumigatus Gα protein GpaB is important for this interaction controlling the production of the secondary metabolite pyripyropene. We demonstrate that the lack of pyripyropene production decreases the dual biofilm interaction between the two species as well as the virulence of A. fumigatus in a chemotherapeutic murine model of aspergillosis. These results reveal a complete novel role for this secondary metabolite in the ecology and pathogenic interactions of this important human fungal pathogen.

Description of two new species of Ophiocordyceps: O.sinocampes and O.cystidiata (Ophiocordycipitaceae, Hypocreales) from typical karst landform forests in Guizhou, China

Karst habitats are hotspots of diversity and endemism. Their naturally fragmented distributions across broad geographic landscapes have led to a complex array of smaller evolutionary ecosystems. Comprehensive biodiversity assessments of karst habitats have revealed that these ecosystems contain a high level of endemism. During a survey of cordycipitoid fungi in the karst region of southwest Guizhou, China, we identified and proposed two new species, based on morphology and multi-locus (nrSSU, nrLSU, tef, rpb1, and rpb2) phylogenetic analyses. O.cystidiata sp. nov. is characterized by gray-white to yellow fertile part, verrucose phialides, and conidia enveloped in a cystic thickened mucus sheath, distinguishing it from other species. O.sinocampes sp. nov. is distinguished by long tapering phialides with inflated bases. Phylogenetic analyses using five loci reveal that O.cystidiata shares a close relationship with O.fenggangensis, O.musicaudata, O.alboperitheciata, and Hirsutellakuankuoshuiensis, while O.sinocampes is closely related to O.multiperitheciata and H.strigosa. Herein, we completed the descriptions, illustrations and molecular phylogeny of these two new species. The host diversity of O.sinocampes has also been documented within the orders Lepidoptera and Coccoidea. Our research further enriches the diversity of cordycipitoid species in the typical karst landform areas of Guizhou, China.

1,2,4-Oxadiazole-5-Carboxylic Acid Derivatives as Safe Seed Treatment Nematicides

Plant-parasitic nematodes pose a significant threat to crop production, impacting agricultural yields. In the search for new nematicides, a series of 1,2,4-oxadiazole-5-carboxylic acid derivatives containing amide or ester groups were designed and synthesized using an activity-based approach. Bioassay results showed that some compounds exhibited good nematicidal activity against Meloidogyne incognita, Aphelenchoides besseyi, and Bursaphelenchus xylophilus. Especially, compound f1 showed excellent nematicidal activity against A. besseyi with an LC50 value of 19.0 μg/mL at 48 h, outperforming tioxazafen (LC50 = 149 μg/mL) and fosthiazate (LC50 > 300 μg/mL). In semi-in vivo and pot experiments, compound f1 demonstrated superior nematicidal activity compared to fosthiazate against A. besseyi while also being safe for rice seeds. The mechanism of action revealed that compound f1 could inhibit the population and behavior of A. besseyi, damage the nematode cuticle, resulting in fluid leakage, affect the production of reactive oxygen species, lipofuscin, and lipids in nematodes, and inhibit acetylcholinesterase in A. besseyi. These findings suggest that compound f1 has the potential to be a safe and effective seed treatment nematicide.

Phylogeny and Systematics of the Genus Tolypocladium (Ophiocordycipitaceae, Hypocreales)

The taxonomy and phylogeny of the genus Tolypocladium are herein revised based on the most comprehensive dataset to date. Two species-level phylogenies of Tolypocladium were constructed: a single-gene phylogeny (ITS) of 35 accepted species and a multigene phylogeny (nrSSU, nrLSU, tef-1α, rpb1, and rpb2) of 27 accepted species. Three new species, Tolypocladium pseudoalbum sp. nov., Tolypocladium subparadoxum sp. nov., and Tolypocladium yunnanense sp. nov., are described in the present study. The genetic divergences of four markers (ITS, tef-1α, rpb1 and rpb2) among Tolypocladium species are also reported. The results indicated that species of Tolypocladium were best delimited by rpb1 sequence data, followed by the sequence data for the rpb2, tef-1α, and ITS provided regions. Finally, a key to the 48 accepted species of Tolypocladium worldwide is provided.

The fucose-specific lectin gene AOL_s00054g276 affects trap formation and nematocidal activity of the nematophagous fungus Arthrobotrys oligospora

Nematode-trapping fungi are natural enemies of nematodes in nature. Arthrobotrys oligospora, a typical nematode-trapping fungus with a clear genetic background, can capture and infect nematodes by forming adhesive three-dimensional networks. Lectins, a class of glycoproteins containing glycosyl-specific recognition domains, play an important role in biological recognition. However, the fucose-specific lectins have rarely been studied regarding the process of preying on nematodes. In this study, we characterized the biological role of the fucose-specific lectin encoding gene AOL_s00054g276 (g276) in A. oligospora. The gene g276 was first deleted based on homologous recombination, then the phenotype and nematocidal activity of the Δg276 mutant was evaluated. The results showed that the deletion of gene g276 delayed trap formation and weakened its nematocidal activity; however, mycelial growth, conidia production, conidial germination rates, and adaption to environmental stresses were not affected. Our results suggest that the fucose-specific lectin encoding gene g276 might be associated with the morphogenesis of this fungus, and its deletion resulted in a significantly low density of three-dimensional traps (P < 0.05) and a significantly low nematode-tapping efficiency (P < 0.001). These findings provide a basis for further elucidating the mechanism of A. oligospora preying on nematodes and lay a foundation for the development and utilization of fungal-derived lectins for nematode control in future.

Genome Size Versus Genome Assemblies: Are the Genomes Truly Expanded in Polyploid Fungal Symbionts?

Each day, as the amount of genomic data and bioinformatics resources grows, researchers are increasingly challenged with selecting the most appropriate approach to analyze their data. In addition, the opportunity to undertake comparative genomic analyses is growing rapidly. This is especially true for fungi due to their small genome sizes (i.e., mean 1C = 44.2 Mb). Given these opportunities and aiming to gain novel insights into the evolution of mutualisms, we focus on comparing the quality of whole genome assemblies for fungus-growing ants cultivars (Hymenoptera: Formicidae: Attini) and a free-living relative. Our analyses reveal that currently available methodologies and pipelines for analyzing whole-genome sequence data need refining. By using different genome assemblers, we show that the genome assembly size depends on what software is used. This, in turn, impacts gene number predictions, with higher gene numbers correlating positively with genome assembly size. Furthermore, the majority of fungal genome size data currently available are based on estimates derived from whole-genome assemblies generated from short-read genome data, rather than from the more accurate technique of flow cytometry. Here, we estimated the haploid genome sizes of three ant fungal symbionts by flow cytometry using the fungus Pleurotus ostreatus (Jacq.) P. Kumm. (1871) as a calibration standard. We found that published genome sizes based on genome assemblies are 2.5- to 3-fold larger than our estimates based on flow cytometry. We, therefore, recommend that flow cytometry is used to precalibrate genome assembly pipelines, to avoid incorrect estimates of genome sizes and ensure robust assemblies.

The Role of Melanin in the Biology and Ecology of Nematophagous Fungi

Melanin is a heteropolymer formed by the polymerization of phenolic and indolic compounds. It occurs in organisms across all biological kingdoms and has a range different of functions, thus indicating its important evolutionary role. The presence of melanin offers several protective advantages, including against ultraviolet radiation, traumatic damage, oxidative stress, extreme temperatures, and pressure. For many species of fungi, melanin also participates directly in the process of virulence and pathogenicity. These organisms can synthesize melanin in two main ways: using a substrate of endogenous origin, involving 1,8-dihydroxynaphthalene (DHN); alternatively, in an exogenous manner with the addition of L-3, 4-dihydroxyphenylalanine (L-DOPA or levodopa). As melanin is an amorphous and complex substance, its study requires expensive and inaccessible technologies and analyses are often difficult to perform with conventional biochemical techniques. As such, details about its chemical structure are not yet fully understood, particularly for nematophagous fungi that remain poorly studied. Thus, this review presents an overview of the different types of melanin, with an emphasis on fungi, and discusses the role of melanin in the biology and ecology of nematophagous fungi.

Current State and Future Directions of Genetics and Genomics of Endophytic Fungi for Bioprospecting Efforts

The bioprospecting of secondary metabolites from endophytic fungi received great attention in the 1990s and 2000s, when the controversy around taxol production from Taxus spp. endophytes was at its height. Since then, hundreds of reports have described the isolation and characterization of putative secondary metabolites from endophytic fungi. However, only very few studies also report the genetic basis for these phenotypic observations. With low sequencing cost and fast sample turnaround, genetics- and genomics-based approaches have risen to become comprehensive approaches to study natural products from a wide-range of organisms, especially to elucidate underlying biosynthetic pathways. However, in the field of fungal endophyte biology, elucidation of biosynthetic pathways is still a major challenge. As a relatively poorly investigated group of microorganisms, even in the light of recent efforts to sequence more fungal genomes, such as the 1000 Fungal Genomes Project at the Joint Genome Institute (JGI), the basis for bioprospecting of enzymes and pathways from endophytic fungi is still rather slim. In this review we want to discuss the current approaches and tools used to associate phenotype and genotype to elucidate biosynthetic pathways of secondary metabolites in endophytic fungi through the lens of bioprospecting. This review will point out the reported successes and shortcomings, and discuss future directions in sampling, and genetics and genomics of endophytic fungi. Identifying responsible biosynthetic genes for the numerous secondary metabolites isolated from endophytic fungi opens the opportunity to explore the genetic potential of producer strains to discover novel secondary metabolites and enhance secondary metabolite production by metabolic engineering resulting in novel and more affordable medicines and food additives.

Genome-Wide Identification and Analysis of Chitinase GH18 Gene Family in Mycogone perniciosa

Mycogone perniciosa causes wet bubble disease in Agaricus bisporus and various Agaricomycetes species. In a previous work, we identified 41 GH18 chitinase genes and other pathogenicity-related genes in the genome of M. perniciosa Hp10. Chitinases are enzymes that degrade chitin, and they have diverse functions in nutrition, morphogenesis, and pathogenesis. However, these important genes in M. perniciosa have not been fully characterized, and their functions remain unclear. Here, we performed a genome-wide analysis of M. perniciosa GH18 genes and analyzed the transcriptome profiles and GH18 expression patterns in M. perniciosa during the time course of infection in A. bisporus. Phylogenetic analysis of the 41 GH18 genes with those of 15 other species showed that the genes were clustered into three groups and eight subgroups based on their conserved domains. The GH18 genes clustered in the same group shared different gene structures but had the same protein motifs. All GH18 genes were localized in different organelles, were unevenly distributed on 11 contigs, and had orthologs in the other 13 species. Twelve duplication events were identified, and these had undergone both positive and purifying selection. The transcriptome analyses revealed that numerous genes, including transporters, cell wall degrading enzymes (CWDEs), cytochrome P450, pathogenicity-related genes, secondary metabolites, and transcription factors, were significantly upregulated at different stages of M. perniciosa Hp10 infection of A. bisporus. Twenty-three out of the 41 GH18 genes were differentially expressed. The expression patterns of the 23 GH18 genes were different and were significantly expressed from 3 days post-inoculation of M. perniciosa Hp10 in A. bisporus. Five differentially expressed GH18 genes were selected for RT-PCR and gene cloning to verify RNA-seq data accuracy. The results showed that those genes were successively expressed in different infection stages, consistent with the previous sequencing results. Our study provides a comprehensive analysis of pathogenicity-related and GH18 chitinase genes’ influence on M. perniciosa mycoparasitism of A. bisporus. Our findings may serve as a basis for further studies of M. perniciosa mycoparasitism, and the results have potential value for improving resistance in A. bisporus and developing efficient disease-management strategies to mitigate wet bubble disease.

Fungi-Nematode Interactions: Diversity, Ecology, and Biocontrol Prospects in Agriculture

Fungi and nematodes are among the most abundant organisms in soil habitats. They provide essential ecosystem services and play crucial roles for maintaining the stability of food-webs and for facilitating nutrient cycling. As two of the very abundant groups of organisms, fungi and nematodes interact with each other in multiple ways. Here in this review, we provide a broad framework of interactions between fungi and nematodes with an emphasis on those that impact crops and agriculture ecosystems. We describe the diversity and evolution of fungi that closely interact with nematodes, including food fungi for nematodes as well as fungi that feed on nematodes. Among the nematophagous fungi, those that produce specialized nematode-trapping devices are especially interesting, and a great deal is known about their diversity, evolution, and molecular mechanisms of interactions with nematodes. Some of the fungi and nematodes are significant pathogens and pests to crops. We summarize the ecological and molecular mechanisms identified so far that impact, either directly or indirectly, the interactions among phytopathogenic fungi, phytopathogenic nematodes, and crop plants. The potential applications of our understanding to controlling phytophagous nematodes and soilborne fungal pathogens in agricultural fields are discussed.

A New High-Quality Draft Genome Assembly of the Chinese Cordyceps Ophiocordyceps sinensis

Ophiocordyceps sinensis (Berk.) is an entomopathogenic fungus endemic to the Qinghai-Tibet Plateau. It parasitizes and mummifies the underground ghost moth larvae, then produces a fruiting body. The fungus-insect complex, called Chinese cordyceps or "DongChongXiaCao," is not only a valuable traditional Chinese medicine, but also a major source of income for numerous Himalayan residents. Here, taking advantage of rapid advances in single-molecule sequencing, we assembled a highly contiguous genome assembly of O. sinensis. The assembly of 23 contigs was ∼110.8 Mb with a N50 length of 18.2 Mb. We used RNA-seq and homologous protein sequences to identify 8,916 protein-coding genes in the IOZ07 assembly. Moreover, 63 secondary metabolite gene clusters were identified in the improved assembly. The improved assembly and genome features described in this study will further inform the evolutionary study and resource utilization of Chinese cordyceps.

Comparative Transcriptomics Reveals Features and Possible Mechanisms of Glucose-Mediated Soil Fungistasis Relief in Arthrobotrys oligospora

Soil-borne pest diseases result in large annual agricultural losses globally. Fungal bio-control agents are an alternative means of controlling pest diseases; however, soil fungistasis limits the effect of fungal agents. Nutrients can relieve soil fungistasis, but the mechanisms behind this process remain poorly understood. In this study, we determined and quantified the transcriptomes of Arthrobotrys oligospora, a nematode-trapping fungus, derived from samples of fresh conidia, germinated conidia, soil fungistatic conidia, and glucose-relieved conidia. The transcriptomes of fungistatic and glucose-relieved conidia were significantly different from those of the other two conidia samples. KEGG pathway analyses showed that those genes upregulated in fungistatic and glucose-relieved conidia were mainly involved in translation and substance metabolism, and the downregulated genes were mainly involved in MAPK pathway, autophagy, mitophagy, and endocytosis. As being different from the transcriptome of fungistatic conidia, upregulated genes in the transcriptome of glucose-relieved conidia are also related to replication and repair, spliceosome, oxidative phosphorylation, autophagy, and degradation pathway (lysosome, proteasome, and RNA degradation). And the upregulated genes resulted from comparison of glucose-relieved conidia and fungistatic conidia were enriched in metabolic pathways, cycle, DNA replication, and repair. The differentially splicing events in the transcriptome of glucose-relieved conidia are far more than that of other two transcriptomes, and genes regulated by differentially splicing were analyzed through KEGG pathway analysis. Furthermore, autophagy genes were proved to play important role in resisting soil fungistasis and glucose-mediated soil fungistasis relief. These data indicate that, in addition to being a carbon and energy source for conidia germination, glucose may also help to relieve soil fungistasis by activating many cellular processes, including autophagy, DNA replication and repair, RNA alternative splicing, and degradation pathways.

Analysis of the Indole Diterpene Gene Cluster for Biosynthesis of the Epoxy-Janthitrems in Epichloë Endophytes

Epoxy-janthitrems are a class of indole diterpenes with structural similarity to lolitrem B. Two taxa of asexual Epichloë endophytes have been reported to produce epoxy-janthitrems, LpTG-3 (Lolium perenne Taxonomic Group 3; e.g., NEA12) and LpTG-4 (e.g., E1). Epichloë epoxy-janthitrems are not well understood, the biosynthetic pathway and associated gene complement have not been described and while the literature suggests they are associated with superior protection against pasture insect pests and are tremorgenic in grazing mammals, these properties have not been confirmed using isolated and purified compounds. Whole genome sequence analysis was used to identify candidate genes for epoxy-janthitrem biosynthesis that are unique to epoxy-janthitrem producing strains of Epichloë. A gene, jtmD, was identified with homology to aromatic prenyl transferases involved in synthesis of indole diterpenes. The location of the epoxy-janthitrem biosynthesis gene cluster (JTM locus) was determined in the assembled nuclear genomes of NEA12 and E1. The JTM locus contains cluster 1 and cluster 2 of the lolitrem B biosynthesis gene cluster (LTM locus), as well as four genes jtmD, jtmO, jtm01, and jtm02 that are unique to Epichloë spp. that produce epoxy-janthitrems. Expression of each of the genes identified was confirmed using transcriptome analysis of perennial ryegrass-NEA12 and perennial ryegrass-E1 symbiota. Sequence analysis confirmed the genes are functionally similar to those involved in biosynthesis of related indole diterpene compounds. RNAi silencing of jtmD and in planta assessment in host-endophyte associations confirms the role of jtmD in epoxy-janthitrem production. Using LCMS/MS technologies, a biosynthetic pathway for the production of epoxy-janthitrems I-IV in Epichloë endophytes is proposed.

Fungal evolution: major ecological adaptations and evolutionary transitions

Fungi are a highly diverse group of heterotrophic eukaryotes characterized by the absence of phagotrophy and the presence of a chitinous cell wall. While unicellular fungi are far from rare, part of the evolutionary success of the group resides in their ability to grow indefinitely as a cylindrical multinucleated cell (hypha). Armed with these morphological traits and with an extremely high metabolical diversity, fungi have conquered numerous ecological niches and have shaped a whole world of interactions with other living organisms. Herein we survey the main evolutionary and ecological processes that have guided fungal diversity. We will first review the ecology and evolution of the zoosporic lineages and the process of terrestrialization, as one of the major evolutionary transitions in this kingdom. Several plausible scenarios have been proposed for fungal terrestralization and we here propose a new scenario, which considers icy environments as a transitory niche between water and emerged land. We then focus on exploring the main ecological relationships of Fungi with other organisms (other fungi, protozoans, animals and plants), as well as the origin of adaptations to certain specialized ecological niches within the group (lichens, black fungi and yeasts). Throughout this review we use an evolutionary and comparative-genomics perspective to understand fungal ecological diversity. Finally, we highlight the importance of genome-enabled inferences to envision plausible narratives and scenarios for important transitions.

Comparative evolutionary histories of fungal proteases reveal gene gains in the mycoparasitic and nematode-parasitic fungus Clonostachys rosea

Background

The ascomycete fungus Clonostachys rosea (order Hypocreales) can control several important plant diseases caused by plant pathogenic fungi and nematodes. Subtilisin-like serine proteases are considered to play an important role in pathogenesis in entomopathogenic, mycoparasitic, and nematophagous fungi used for biological control. In this study, we analysed the evolutionary histories of protease gene families, and investigated sequence divergence and regulation of serine protease genes in C. rosea.

Results

Proteases of selected hypocrealean fungal species were classified into families based on the MEROPS peptidase database. The highest number of protease genes (590) was found in Fusarium solani, followed by C. rosea with 576 genes. Analysis of gene family evolution identified non-random changes in gene copy numbers in the five serine protease gene families S1A, S8A, S9X, S12 and S33. Four families, S1A, S8A, S9X, and S33, displayed gene gains in C. rosea. A gene-tree / species-tree reconciliation analysis of the S8A family revealed that the gene copy number increase in C. rosea was primarily associated with the S08.054 (proteinase K) subgroup. In addition, regulatory and predicted structural differences, including twelve sites evolving under positive selection, among eighteen C. rosea S8A serine protease paralog genes were also observed. The C. rosea S8A serine protease gene prs6 was induced during interaction with the plant pathogenic species F. graminearum.

Conclusions

Non-random increases in S8A, S9X and S33 serine protease gene numbers in the mycoparasitic species C. rosea, Trichoderma atroviride and T. virens suggests an involvement in fungal-fungal interactions. Regulatory and predicted structural differences between C. rosea S8A paralogs indicate that functional diversification is driving the observed increase in gene copy numbers. The induction of prs6 expression in C. rosea during confrontation with F. graminearum suggests an involvement of the corresponding protease in fungal-fungal interactions. The results pinpoint the importance of serine proteases for ecological niche adaptation in C. rosea, including a potential role in the mycoparasitic attack on fungal prey.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1291-1) contains supplementary material, which is available to authorized users.

Introgression and gene family contraction drive the evolution of lifestyle and host shifts of hypocrealean fungi

Hypocrealean fungi (Ascomycota) are known for their diversity of lifestyles. Their vital influences on agricultural and natural ecosystems have resulted in a number of sequenced genomes, which provide essential data for genomic analysis. Totally, 45 hypocrealean fungal genomes constructed a phylogeny. The phylogeny showed that plant pathogens in Nectriaceae diverged earliest, followed by animal pathogens in Cordycipitaceae, Ophiocordycipitaceae and Clavicipitaceae with mycoparasites in Hypocreaceae. Insect/nematode pathogens and grass endophytes in Clavicipitaceae diverged at last. Gene families associated with host-derived nutrients are significantly contracted in diverged lineages compared with the ancestral species. Introgression was detected in certain lineages of hypocrealean fungi, and the main functions of the genes located in the introgressed regions are involved in host recognition, transcriptional regulation, stress response and cell growth regulation. These results indicate that contraction of gene families and introgression might be main mechanisms to drive lifestyle differentiation and evolution and host shift of hypocrealean fungi.

Genome and secretome analysis of Pochonia chlamydosporia provide new insight into egg-parasitic mechanisms

Pochonia chlamydosporia infects eggs and females of economically important plant-parasitic nematodes. The fungal isolates parasitizing different nematodes are genetically distinct. To understand their intraspecific genetic differentiation, parasitic mechanisms, and adaptive evolution, we assembled seven putative chromosomes of P. chlamydosporia strain 170 isolated from root-knot nematode eggs (~44 Mb, including 7.19% of transposable elements) and compared them with the genome of the strain 123 (~41 Mb) isolated from cereal cyst nematode. We focus on secretomes of the fungus, which play important roles in pathogenicity and fungus-host/environment interactions, and identified 1,750 secreted proteins, with a high proportion of carboxypeptidases, subtilisins, and chitinases. We analyzed the phylogenies of these genes and predicted new pathogenic molecules. By comparative transcriptome analysis, we found that secreted proteins involved in responses to nutrient stress are mainly comprised of proteases and glycoside hydrolases. Moreover, 32 secreted proteins undergoing positive selection and 71 duplicated gene pairs encoding secreted proteins are identified. Two duplicated pairs encoding secreted glycosyl hydrolases (GH30), which may be related to fungal endophytic process and lost in many insect-pathogenic fungi but exist in nematophagous fungi, are putatively acquired from bacteria by horizontal gene transfer. The results help understanding genetic origins and evolution of parasitism-related genes.

Comparative genomic analyses reveal the features for adaptation to nematodes in fungi

Nematophagous (NP) fungi are ecologically important components of the soil microbiome in natural ecosystems. Esteya vermicola (Ev) has been reported as a NP fungus with a poorly understood evolutionary history and mechanism of adaptation to parasitism. Furthermore, NP fungal genomic basis of lifestyle was still unclear. We sequenced and annotated the Ev genome (34.2 Mbp) and integrated genetic makeup and evolution of pathogenic genes to investigate NP fungi. The results revealed that NP fungi had some abundant pathogenic genes corresponding to their niche. A number of gene families involved in pathogenicity were expanded, and some pathogenic orthologous genes underwent positive selection. NP fungi with diverse morphological features exhibit similarities of evolutionary convergence in attacking nematodes, but their genetic makeup and microscopic mechanism are different. Endoparasitic NP fungi showed similarity in large number of transporters and secondary metabolite coding genes. Noteworthy, expanded families of transporters and endo-beta-glucanase implied great genetic potential of Ev in quickly perturbing nematode metabolism and parasitic behavior. These results facilitate our understanding of NP fungal genomic features for adaptation to nematodes and lay a solid theoretical foundation for further research and application.

Evaluation of Clonostachys rosea for Control of Plant-Parasitic Nematodes in Soil and in Roots of Carrot and Wheat

Biological control is a promising approach to reduce plant diseases caused by nematodes. We tested the effect of the fungus Clonostachys rosea strain IK726 inoculation on nematode community composition in a naturally nematode infested soil in a pot experiment, and the effect of C. rosea on plant health. The numbers of plant-parasitic nematode genera extracted from soil and plant roots decreased by 40 to 73% when C. rosea was applied, while genera of nonparasitic nematodes were not affected. Soil inoculation of C. rosea increased fresh shoot weight and shoot length of wheat plants by 20 and 24%, respectively, while only shoot dry weight increased by 48% in carrots. Light microscopy of in vitro C. rosea-nematode interactions did not reveal evidence of direct parasitism. However, culture filtrates of C. rosea growing in potato dextrose broth, malt extract broth and synthetic nutrient broth exhibited toxicity toward nematodes and immobilized 57, 62, and 100% of the nematodes, respectively, within 48 h. This study demonstrates that C. rosea can control plant-parasitic nematodes and thereby improve plant growth. The most likely mechanism responsible for the antagonism is antibiosis through production of nematicidal compounds, rather than direct parasitism.

Ant-infecting Ophiocordyceps genomes reveal a high diversity of potential behavioral manipulation genes and a possible major role for enterotoxins

Much can be gained from revealing the mechanisms fungal entomopathogens employ. Especially intriguing are fungal parasites that manipulate insect behavior because, presumably, they secrete a wealth of bioactive compounds. To gain more insight into their strategies, we compared the genomes of five ant-infecting Ophiocordyceps species from three species complexes. These species were collected across three continents, from five different ant species in which they induce different levels of manipulation. A considerable number of (small) secreted and pathogenicity-related proteins were only found in these ant-manipulating Ophiocordyceps species, and not in other ascomycetes. However, few of those proteins were conserved among them, suggesting that several different methods of behavior modification have evolved. This is further supported by a relatively fast evolution of previously reported candidate manipulation genes associated with biting behavior. Moreover, secondary metabolite clusters, activated during biting behavior, appeared conserved within a species complex, but not beyond. The independent co-evolution between these manipulating parasites and their respective hosts might thus have led to rather diverse strategies to alter behavior. Our data indicate that specialized, secreted enterotoxins may play a major role in one of these strategies.

Esteya Vermicola, a Nematophagous Fungus Attacking the Pine Wood Nematode, Harbors a Bacterial Endosymbiont Affiliated with Gammaproteobacteria

Symbioses have played pivotal roles in biological, ecological, and evolutionary diversification. Symbiotic bacteria affect the biology of hosts in a number of ways. Esteya vermicola, an endoparasitic nematophagous fungus, has high infectivity in the pine wood nematode (PWN), which causes devastating ecological damage and economic losses in Asia and Europe. An integration of molecular, phylogenetic, and morphological analyses revealed that surface-sterilized E. vermicola with septate hyphae from different geographic locations harbor bacterial endosymbionts. 16S rRNA gene sequences from four fungal strains all clustered in a well-supported monophyletic clade that was the most closely related to Pseudomonas stutzeri and affiliated with Gammaproteobacteria. The existence and intracellular location of endobacteria was revealed by fluorescent in situ hybridization (FISH). Our results showed that endobacteria were coccoid, vertically inherited, as yet uncultured, and essential symbionts. Ultrastructural observations indicated that young and old endobacteria differed in cell size, cell wall thickness, and the degree of reproduction. The results of the present study provide a fundamental understanding of the endobacteria inside E. vermicola and raise questions regarding the impact of endobacteria on the biology, ecology, and evolution of their fungal host.

The conidial mucilage, natural film coatings, is involved in environmental adaptability and pathogenicity of Hirsutella satumaensis Aoki

The Hirsutella genus is very special asexually-reproducing pathogens of insects by reduced sporulation, host specificity and spores covered by a thick mucilage layer. However, the ecological function of conidial mucilage remains elusive. In this study, the possible ecological role of conidial mucilage from the entomopathogenic fungus Hirsutella satumaensis was functionally investigated through tolerance, adherence and insect bioassays involving aerial conidia (AC) and mucilage-free conidia (MFC). Measurements of hydrophobicity using microbial adhesion to hydrocarbons (MATH) indicated that mucilage is main contributor to the surface hydrophobicity of AC. When subjected in tolerance assays to extreme temperatures, high chemical pressure, extended exposure to ultraviolet radiation and cold stress, AC produced more colonies, exhibited higher conidiation and germination percentages than those of MFC. In adhesion assays, MFC displayed an approximately 40% reduction in adherence to locust, dragonfly cuticle and onion epidermis when washed with 0.05% Tween 20. Similarly, Galleria mellonella and Plutella xylostella larvae infected with mucilage-producing AC experienced a relatively higher mortality rate. Our findings suggest that mucilage is critical to the ecological adaptability of H. satumaensis, where it plays positive roles on maintenance of spore surface hydrophobicity, enhancement of spore resistance to extreme environments and strengthening of spore adhesion and host pathogenicity.

Managing the plant microbiome for biocontrol fungi: examples from Hypocreales

Feeding an increasing global population requires continued improvements in agricultural efficiency and productivity. Meeting estimated future production levels requires the adoption of practices that increase output without environmental degradation associated with external inputs to supplement nutrition or control pests. Enriching the community of microbes associated with plants in agricultural systems for those providing ecosystem services such as pest control is one possible component towards achieving sustainable productivity increases. In this review we explore the current state of knowledge for Hypocreales fungi used in biological control. Advances in understanding the field ecology, diversity and genetic determinants of host range and virulence of hypocrealean fungi provide the means to improve their efficacy.

Mitochondrial genome of the nematode endoparasitic fungus Hirsutella vermicola reveals a high level of synteny in the family Ophiocordycipitaceae

Ophiocordycipitaceae is a diverse fungal family comprising multiple ecologically, economically, medicinally, and culturally important fungal species; however, only four species of the family have available mitochondrial genomes (mitogenomes). In this study, the complete mitogenome of the nematode endoparasitic fungus Hirsutella vermicola in Ophiocordycipitaceae was sequenced, and a comparative mitogenomic analysis of Ophiocordycipitaceae was performed. We found that the 53,793-bp circular mitogenome of H. vermicola, except for standard fungal mitochondrial genes, harbors seven introns acquired possibly through lateral transfer from other fungi and three free-standing open reading frames (ORFs) coding for hypothetical proteins. Phylogenetic analysis based on concatenated mitochondrial protein sequences confirmed its placement in Ophiocordycipitaceae. Comparison on five mitogenomes of Ophiocordycipitaceae revealed great variation on their sizes, from 35.2 kb in Tolypocladium ophioglossoides to 157.5 kb in Ophiocordyceps sinensis, mainly due to variable numbers of introns (from 7 to 54) as well as variable lengths of intergenic regions. The five mitogenomes, however, are highly syntenic to each other in terms of gene order, the presence of an intronic ORF encoding ribosomal protein S3 within rnl, and the nad2/nad3 joining pattern. Our study is the first report of the mitogenome of H. vermicola and has facilitated the understanding of mitogenome evolution of Ophiocordycipitaceae.

Characterization of a new toxin from the entomopathogenic fungus Metarhizium anisopliae: the ribotoxin anisoplin

Metarhizium anisopliae is an entomopathogenic fungus relevant in biotechnology with applications like malaria vector control. Studies of its virulence factors are therefore of great interest. Fungal ribotoxins are toxic ribonucleases with extraordinary efficiency against ribosomes and suggested as potential insecticides. Here we describe this ribotoxin characteristic activity in M. anisopliae cultures. Anisoplin has been obtained as a recombinant protein and further characterized. It is structurally similar to hirsutellin A, the ribotoxin from the entomopathogen Hirsutella thompsonii. Moreover, anisoplin shows the ribonucleolytic activity typical of ribotoxins and cytotoxicity against insect cells. How Metarhizium uses this toxin and possible applications are of interest.

Phylogenetic analyses reveal molecular signatures associated with functional divergence among Subtilisin like Serine Proteases are linked to lifestyle transitions in Hypocreales

Background

Subtilisin-like serine proteases or Subtilases in fungi are important for penetration and colonization of host. In Hypocreales, these proteins share several properties with other fungal, bacterial, plant and mammalian homologs. However, adoption of specific roles in entomopathogenesis may be governed by attainment of unique biochemical and structural features during the evolutionary course. Due to such functional shifts Subtilases coded by different family members of Hypocreales acquire distinct features according to respective hosts and lifestyle. We conducted phylogenetic and DIVERGE analyses and identified important protein residues that putatively assign functional specificity to Subtilases in fungal families/species under the order Hypocreales.

Results

A total of 161 Subtilases coded by 10 species from five different families under the fungal order Hypocreales was included in the analysis. Based on the presence of conserved domains, the Subtilase genes were divided into three subfamilies, Subtilisin (S08.005), Proteinase K (S08.054) and Serine-carboxyl peptidases (S53.001). These subfamilies were investigated for phylogenetic associations, protein residues under positive selection and functional divergence among paralogous clades. The observations were co-related with the life-styles of the fungal families/species. Phylogenetic and Divergence analyses of Subtilisin (S08.005) and Proteinase K (S08.054) families of proteins revealed that the paralogous clades were clear-cut representation of familial origin of the protein sequences. We observed divergence between the paralogous clades of plant-pathogenic fungi (Nectriaceae), insect-pathogenic fungi (Cordycipitaceae/Clavicipitaceae) and nematophagous fungi (Ophiocordycipitaceae). In addition, Subtilase genes from the nematode-parasitic fungus Purpureocillium lilacinum made a unique cluster which putatively indicated that the fungus might have developed distinctive mechanisms for nematode-pathogenesis. Our evolutionary genetics analysis revealed evidence of positive selection on the Subtilisin (S08.005) and Proteinase K (S08.054) protein sequences of the entomopathogenic and nematophagous species belonging to Cordycipitaceae, Clavicipitaceae and Ophiocordycipitaceae families of Hypocreales.

Conclusions

Our study provided new insights into the evolution of Subtilisin like serine proteases in Hypocreales, a fungal order largely consisting of biological control species. Subtilisin (S08.005) and Proteinase K (S08.054) proteins seemed to play important roles during life style modifications among different families and species of Hypocreales. Protein residues found significant in functional divergence analysis in the present study may provide support for protein engineering in future.

Electronic supplementary material

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

Population Genetics of Hirsutella rhossiliensis, a Dominant Parasite of Cyst Nematode Juveniles on a Continental Scale

Hirsutella rhossiliensis is a parasite of juvenile nematodes, effective against a diversity of plant-parasitic nematodes. Its global distribution on various nematode hosts and its genetic variation for several geographic regions have been reported, while the global population genetic structure and factors underlying patterns of genetic variation of H. rhossiliensis are unclear. In this study, 87 H. rhossiliensis strains from five nematode species (Globodera sp., Criconemella xenoplax, Rotylenchus robustus, Heterodera schachtii, and Heterodera glycines) in Europe, the United States, and China were investigated by multilocus sequence analyses. A total of 280 variable sites (frequency, 0.6%) at eight loci and six clustering in high accordance with geographic populations or host nematode-associated populations were identified. Although H. rhossiliensis is currently recognized as an asexual fungus, recombination events were frequently detected. In addition, significant genetic isolation by geography and nematode hosts was revealed. Overall, our analyses showed that recombination, geographic isolation, and nematode host adaptation have played significant roles in the evolutionary history of H. rhossiliensis IMPORTANCE: H. rhossiliensis has great potential for use as a biocontrol agent to control nematodes in a sustainable manner as an endoparasitic fungus. Therefore, this study has important implications for the use of H. rhossiliensis as a biocontrol agent and provides interesting insights into the biology of this species.

The transcription factor SKN7 regulates conidiation, thermotolerance, apoptotic-like cell death and parasitism in the nematode endoparasitic fungus Hirsutella minnesotensis

The transcription factor SKN7 is a highly conserved protein among fungi and was initially recognized as a response regulator that protects cells from oxidative stress and maintains cell wall integrity in yeast. Orthologs of SKN7 are extensively present in biocontrol agents of plant pathogens, but they had not been functionally characterized. Here, we identified and characterized the transcription factor SKN7 in the nematode endoparasitic fungus Hirsutella minnesotensis. Null mutant lacking HIM-SKN7 (HIM_03620), which was generated by a gene disruption strategy, demonstrated reduced conidiation, increased sensitivity to high temperature, hydrogen peroxide, mannitol and ethanol, and reduced fungal resistance to farnesol. However, over-expression mutant showed increased conidial production, thermotolerance and resistance to farnesol, suggesting that HIM-SKN7 regulates antiapoptotic-like cell death in H. minnesotensis. Moreover, the results showed that in null mutant, H. minnesotensis had decreased endoparasitic ability as compared to wild type and over-expression strain. During the infection process, the relative expression of the HIM-SKN7 gene was significantly induced in the wild type and over-expression strain. The results of the present study advance our understanding of the functions of the SKN7 gene in biocontrol agents, in particular, nematode endoparasitic fungi.

Biosynthesis of Antibiotic Leucinostatins in Bio-control Fungus Purpureocillium lilacinum and Their Inhibition on Phytophthora Revealed by Genome Mining

Purpureocillium lilacinum of Ophiocordycipitaceae is one of the most promising and commercialized agents for controlling plant parasitic nematodes, as well as other insects and plant pathogens. However, how the fungus functions at the molecular level remains unknown. Here, we sequenced two isolates (PLBJ-1 and PLFJ-1) of P. lilacinum from different places Beijing and Fujian. Genomic analysis showed high synteny of the two isolates, and the phylogenetic analysis indicated they were most related to the insect pathogen Tolypocladium inflatum. A comparison with other species revealed that this fungus was enriched in carbohydrate-active enzymes (CAZymes), proteases and pathogenesis related genes. Whole genome search revealed a rich repertoire of secondary metabolites (SMs) encoding genes. The non-ribosomal peptide synthetase LcsA, which is comprised of ten C-A-PCP modules, was identified as the core biosynthetic gene of lipopeptide leucinostatins, which was specific to P. lilacinum and T. ophioglossoides, as confirmed by phylogenetic analysis. Furthermore, gene expression level was analyzed when PLBJ-1 was grown in leucinostatin-inducing and non-inducing medium, and 20 genes involved in the biosynthesis of leucionostatins were identified. Disruption mutants allowed us to propose a putative biosynthetic pathway of leucinostatin A. Moreover, overexpression of the transcription factor lcsF increased the production (1.5-fold) of leucinostatins A and B compared to wild type. Bioassays explored a new bioactivity of leucinostatins and P. lilacinum: inhibiting the growth of Phytophthora infestans and P. capsici. These results contribute to our understanding of the biosynthetic mechanism of leucinostatins and may allow us to utilize P. lilacinum better as bio-control agent.

Comparative Genomic Analysis of Drechmeria coniospora Reveals Core and Specific Genetic Requirements for Fungal Endoparasitism of Nematodes

Drechmeria coniospora is an obligate fungal pathogen that infects nematodes via the adhesion of specialized spores to the host cuticle. D. coniospora is frequently found associated with Caenorhabditis elegans in environmental samples. It is used in the study of the nematode's response to fungal infection. Full understanding of this bi-partite interaction requires knowledge of the pathogen's genome, analysis of its gene expression program and a capacity for genetic engineering. The acquisition of all three is reported here. A phylogenetic analysis placed D. coniospora close to the truffle parasite Tolypocladium ophioglossoides, and Hirsutella minnesotensis, another nematophagous fungus. Ascomycete nematopathogenicity is polyphyletic; D. coniospora represents a branch that has not been molecularly characterized. A detailed in silico functional analysis, comparing D. coniospora to 11 fungal species, revealed genes and gene families potentially involved in virulence and showed it to be a highly specialized pathogen. A targeted comparison with nematophagous fungi highlighted D. coniospora-specific genes and a core set of genes associated with nematode parasitism. A comparative gene expression analysis of samples from fungal spores and mycelia, and infected C. elegans, gave a molecular view of the different stages of the D. coniospora lifecycle. Transformation of D. coniospora allowed targeted gene knock-out and the production of fungus that expresses fluorescent reporter genes. It also permitted the initial characterisation of a potential fungal counter-defensive strategy, involving interference with a host antimicrobial mechanism. This high-quality annotated genome for D. coniospora gives insights into the evolution and virulence of nematode-destroying fungi. Coupled with genetic transformation, it opens the way for molecular dissection of D. coniospora physiology, and will allow both sides of the interaction between D. coniospora and C. elegans, as well as the evolutionary arms race that exists between pathogen and host, to be studied.

Insights into Adaptations to a Near-Obligate Nematode Endoparasitic Lifestyle from the Finished Genome of Drechmeria coniospora

Nematophagous fungi employ three distinct predatory strategies: nematode trapping, parasitism of females and eggs, and endoparasitism. While endoparasites play key roles in controlling nematode populations in nature, their application for integrated pest management is hindered by the limited understanding of their biology. We present a comparative analysis of a high quality finished genome assembly of Drechmeria coniospora, a model endoparasitic nematophagous fungus, integrated with a transcriptomic study. Adaptation of D. coniospora to its almost completely obligate endoparasitic lifestyle led to the simplification of many orthologous gene families involved in the saprophytic trophic mode, while maintaining orthologs of most known fungal pathogen-host interaction proteins, stress response circuits and putative effectors of the small secreted protein type. The need to adhere to and penetrate the host cuticle led to a selective radiation of surface proteins and hydrolytic enzymes. Although the endoparasite has a simplified secondary metabolome, it produces a novel peptaibiotic family that shows antibacterial, antifungal and nematicidal activities. Our analyses emphasize the basic malleability of the D. coniospora genome: loss of genes advantageous for the saprophytic lifestyle; modulation of elements that its cohort species utilize for entomopathogenesis; and expansion of protein families necessary for the nematode endoparasitic lifestyle.

Genome Studies on Nematophagous and Entomogenous Fungi in China

The nematophagous and entomogenous fungi are natural enemies of nematodes and insects and have been utilized by humans to control agricultural and forestry pests. Some of these fungi have been or are being developed as biological control agents in China and worldwide. Several important nematophagous and entomogenous fungi, including nematode-trapping fungi (Arthrobotrys oligospora and Drechslerella stenobrocha), nematode endoparasite (Hirsutella minnesotensis), insect pathogens (Beauveria bassiana and Metarhizium spp.) and Chinese medicinal fungi (Ophiocordyceps sinensis and Cordyceps militaris), have been genome sequenced and extensively analyzed in China. The biology, evolution, and pharmaceutical application of these fungi and their interacting with host nematodes and insects revealed by genomes, comparing genomes coupled with transcriptomes are summarized and reviewed in this paper.

Development of a transformation system for Hirsutella spp. and visualization of the mode of nematode infection by GFP-labeled H. minnesotensis

Hirsutella rhossiliensis and H. minnesotensis are endoparasitic fungi of the second-stage juvenile (J2) of the soybean cyst nematode (Heterodera glycines) in nature. They also parasitize both H. glycines J2 and Caenorhabditis elegans on agar plates. Agrobacterium tumefaciens-mediated transformation conditions were established for these Hirsutella spp. The resulting transformants were similar to the corresponding wild-type strains. The infection processes of H. glycines J2 and C. elegans second larval stage (L2) by H. minnesotensis expressing ZsGreen were microscopically analyzed. Conidia of H. minnesotensis adhered to passing nematodes within 8 h post-inoculation (hpi), formed an infection peg between 8 and 12 hpi, and penetrated the nematode cuticle between 12 and 24 hpi for C. elegans L2 and between 12 and 32 hpi for H. glycines J2. Hyphal proliferation inside of the nematode coelom was observed at approximately 32 hpi for C. elegans L2 and at approximately 40 hpi for H. glycines J2. The fungus consumed the whole body and grew out to produce conidia at approximately 156 and 204 hpi for C. elegans L2 and H. glycines J2, respectively. The efficient transformation protocol and a better understanding of infection process provide a solid foundation for studying the molecular and cellular mechanisms underlying fungal parasitism of nematodes.

The complete mitochondrial genome of the nematode endoparasitic fungus Hirsutella minnesotensis

In this study, the complete mitochondrial genome of a nematode fungal pathogen, Hirsutella minnesotensis, is sequenced. The circular mitogenome of H. minnesotensis is 52,245 bp in length and consists of 2 rRNA, 25 tRNA and 14 standard protein-coding genes of the oxidative phosphorylation system as well as four additional free-standing ORFs encoding for an endonuclease or a hypothetical protein. Seven genes (rnl, cob, cox1, nad3, nad4, nad4L and nad5) are invaded by group I or unclassified introns, and these introns carry ORFs of ribosomal protein S3 and GIY-YIG/LAGLIDADG endonucleases or hypothetical proteins. Evidence of intron degeneration is significant for the nad4L intron and cox1-i1 due to unexpected stop codons and/or frame shifting. Phylogenetic analysis based on concatenated protein sequences confirms H. minnesotensis as a member of the fungal order Hypocreales. In this study, we present the complete mitogenome sequence of H. minnesotensis, which is the first complete mitogenome of the family Ophiocordycipitaceae.

A vision for the innovative study of fungal biology in China: Presidential address

I am proud to be elected as the sixth president of the Mycological Society of China, and highly pleased to have a chance to share my personal opinion here with my fellow mycologists and students regarding the innovative performance of fungal biology studies in China. A stepwise buildup of knowledge and sharp scientific vision is the prerequisite for innovative studies. Taken together with the most advanced techniques and elegant experimental designs, the scholars would have a better chance to acquire novel and conceptual results rather than the “me too” stories by focusing on the mechanisms related with fungal unique biology.