Genome sequence, comparative analysis, and evolutionary insights into chitinases of entomopathogenic fungus Hirsutella thompsonii

Genomic and Transcriptome Analysis Reveals the Biosynthesis Network of Cordycepin in Cordyceps militaris

Cordycepin is the primary active compound of Cordyceps militaris. However, the definitive genetic mechanism governing cordycepin synthesis in fruiting body growth and development remains elusive, necessitating further investigation. This study consists of 64 C. militaris strains collected from northeast China. The high-yielding cordycepin strain CMS19 was selected for the analysis of cordycepin production and the genetic basis of cordycepin anabolism. First, the whole-genome sequencing of CMS19 yielded a final size of 30.96 Mb with 8 contigs and 9781 protein-coding genes. The genome component revealed the presence of four additional secondary metabolite gene clusters compared with other published genomes, suggesting the potential for the production of new natural products. The analyses of evolutionary and genetic differentiation revealed a close relationship between C. militaris and Beauveria bassiana. The population of strains distributed in northeast China exhibited the significant genetic variation. Finally, functional genes associated with cordycepin synthesis were identified using a combination of genomic and transcriptomic analyses. A large number of functional genes associated with energy and purine metabolism were significantly enriched, facilitating the reconstruction of a hypothetical cordycepin metabolic pathway. Therefore, our speculation of the cordycepin metabolism pathway involved 24 genes initiating from the glycolysis and pentose phosphate pathways, progressing through purine metabolism, and culminating in the core region of cordycepin synthesis. These findings could offer fundamental support for scientific utilizations of C. militaris germplasm resources and standardized cultivation for cordycepin production.

The transcriptome of the entomopathogenic fungus Culicinomyces clavisporus contains an ortholog of the insecticidal ribotoxin Hirsutellin

The entomopathogenic fungus Culicinomyces clavisporus is known to infect and kill mosquito larvae and therefore has been seen as a potential biological control agent against disease vector mosquitoes. Whereas most fungal entomopathogens infect hosts by penetrating the external cuticle, C. clavisporus initiates infection through ingestion (per os). This unique infection strategy suggests that the C. clavisporus genome may be mined for novel pathogenicity factors with potential for vector control. To this end, an Isoseq-based transcriptome analysis was initiated, and resulted in a total of 3,512,145 sequences, with an average length of 1,732 bp. Transcripts assembly and annotation suggested that the C. clavisporus transcriptome lacked the cuticle-degrading proteins that have been associated with other entomopathogenic fungi, supporting the per os pathogenicity process. Furthermore, mining of the sequence data unexpectedly revealed C. clavisporus transcripts homologous to the Hirsutellin toxin. Comparative sequence analyses indicated that the C. clavisporus Hirsutellin predicted protein has retained the canonical molecular features that have been associated with the ribotoxic and insecticidal properties of the original toxin isolated from Hirsutella thompsonii. The identification of an Hirsutellin ortholog in C. clavisporus was supported by phylogenetic analyses demonstrating that Culicinomyces and Hirsutella were closely related genera in the Ophiocordycipitaceae family. Validation of the mosquitocidal activity of this novel C. clavisporus protein has yet to be performed but may help position Hirsutellin orthologs as prime candidates for the development of alternative biocontrol approaches complementing the current toolbox of vector mosquito management strategies.

Chitinases as key virulence factors in microbial pathogens: Understanding their role and potential as therapeutic targets

Chitinases are crucial for the survival of bacterial and fungal pathogens both during host infection and outside the host in the environment. Chitinases facilitate adhesion onto host cells, act as virulence factors during infection, and provide protection from the host immune system, making them crucial factors in the survival of microbial pathogens. Understanding the mechanisms behind chitinase action is beneficial to design novel therapeutics to control microbial infections. This review explores the role of chitinases in the pathogenesis of bacterial, fungal, and viral infections. The mechanisms underlying the action of chitinases of bacterial, fungal, and viral pathogens in host cells are thoroughly reviewed. The evolutionary relationships between chitinases of various bacterial, fungal, and viral pathogens are discussed to determine their involvement in processes, such as adhesion and host immune system modulation. Gaining a better understanding of the distribution and activity of chitinases in these microbial pathogens can help elucidate their role in the invasion and infection of host cells.

Identification, Culture Characteristics and Whole-Genome Analysis of Pestalotiopsis neglecta Causing Black Spot Blight of Pinus sylvestris var. mongolica

Black spot needle blight is a serious conifer disease of Pinus sylvestris var. mongolica occurring in Northeast China, which is usually caused by the plant pathogenic fungus Pestalotiopsis neglecta. From the diseased pine needles collected in Honghuaerji, the P. neglecta strain YJ-3 was isolated and identified as the phytopathogen, and its culture characteristics were studied. Then, we generated a highly contiguous 48.36-Mbp genome assembly (N50 = 6.62 Mbp) of the P. neglecta strain YJ-3 by combining the PacBio RS II Single Molecule Real Time (SMRT) and Illumina HiSeq X Ten sequencing platforms. The results showed that a total of 13,667 protein-coding genes were predicted and annotated using multiple bioinformatics databases. The genome assembly and annotation resource reported here will be useful for the study of fungal infection mechanisms and pathogen-host interaction.

Genomic Determinants of Entomopathogenic Fungi and Their Involvement in Pathogenesis

Entomopathogenic fungi offer an effective and eco-friendly alternative to curb insect populations in biocontrol strategy. The evolutionary history of selected entomopathogenic fungi indicates their ancestral relationship with plant endophytes. During this host shifting, entomopathogenic fungi must have acquired multiple mechanisms, including a combination of various biomolecules that make them distinguishable from other fungi. In this review, we focus on understanding various biochemical and molecular mechanisms involved in entomopathogenesis. In particular, we attempt to explain the indispensable role of enlarged gene families of various virulent factors, viz. chitinases, proteases, lipases, specialized metabolites, and cytochrome P450, in entomopathogenesis. Our analysis suggests that entomopathogenic fungi recruit a different set of gene products during the progression of pathogenesis. Knowledge of these bio-molecular interactions between fungi and insect hosts will allow researchers to execute pointed efforts towards the development of improved entomopathogenic fungal strains.

Study of Pathogenesis Using Fluorescent Strain of Cordyceps farinosa Revealed Infection of Thitarodes armoricanus Larvae via Digestive Tract

Cordyceps farinosa is often utilized as a biocontrol agent because of its wide host range, strong lethality, and safety for mammals. Artificial rearing of Thitarodes armoricanus larvae is a prerequisite for the artificial cultivation of Chinese cordyceps, and C. farinosa is the most lethal pathogenic fungus during the rearing process. However, the infection process of C. farinosa is still unclear. In this study, we cloned the promoter of the C. farinosa glyceraldehyde 3-phosphate dehydrogenase gene, constructed the EGFP expression cassette, and integrated it into the C. farinosa genome via Agrobacterium transformation. We obtained a fluorescent strain for better observation of the infection process. Using two different inoculation methods of the fluorescent strain, we observed the traditional infection process through the body surface as well as through the digestive tract via feeding. Both infection modes can lead to larval death and mummification. Our findings demonstrated that during the artificial rearing of T. armoricanus, preventing C. farinosa pollution should be an important part of the disinfection of the rearing environment.

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.

Entomopathogenic fungi Beauveria bassiana and Metarhizium anisopliae play roles of maize (Zea mays) growth promoter

Beauveria bassiana and Metarhizium anisopliae are two of the most important and widely used entomopathogenic fungi (EPFs) to control insect pests. Recent studies have revealed their function in promoting plant growth after artificial inoculation. To better assess fungal colonization and growth-promoting effects of B. bassiana and M. anisopliae on crops, maize Zea mays seedlings were treated separately with 13 B. bassiana and 73 M. anisopliae as rhizosphere fungi in a hydroponic cultural system. Plant growth indexes, including plant height, root length, fresh weight, etc., were traced recorded for 35 days to prove the growth promoting efficiency of the EPFs inoculation. Fungal recovery rate (FRR) verified that both B. bassiana and M. anisopliae could endophytically colonize in maize tissues. The recovery rates of B. bassiana in stems and leaves were 100% on the 7th day, but dropped to 11.1% in the stems and 22.2% in the leaves on the 28th day. Meanwhile, B. bassiana was not detected in the roots until the 28th day, reaching a recovery rate of 33.3%. M. anisopliae strains were isolated from the plant roots, stems and leaves throughout the tracing period with high recovery rates. The systematical colonization of B. bassiana and M. anisopliae in different tissues were further corroborated by PCR amplification of fungus-specified DNA band, which showed a higher detection sensitivity of 100% positive reaction. Fungal density comparing to the initial value in the hydroponic solution, dropped to be well below 1% on the 21st day. Thus, the two selected entomopathogenic fungal strains successfully established endophytic colonization rather than rhizospheric colonization in maize, and significantly promoted its growth in a hydroponic cultural system. Entomopathogenic fungi have great application potential in eco-agricultural fields including biopesticides and biofertilizers.

Microbial polyketides and their roles in insect virulence: from genomics to biological functions

Covering: May 1966 up to January 2022Entomopathogenic microorganisms have potential for biological control of insect pests. Their main secondary metabolites include polyketides, nonribosomal peptides, and polyketide-nonribosomal peptide (PK-NRP) hybrids. Among these secondary metabolites, polyketides have mainly been studied for structural identification, pathway engineering, and for their contributions to medicine. However, little is known about the function of polyketides in insect virulence. This review focuses on the role of bacterial and fungal polyketides, as well as PK-NRP hybrids in insect infection and killing. We also discuss gene distribution and evolutional relationships among different microbial species. Further, the role of microbial polyketides and the hybrids in modulating insect-microbial symbiosis is also explored. Understanding the mechanisms of polyketides in insect pathogenesis, how compounds moderate the host-fungus interaction, and the distribution of PKS genes across different fungi and bacteria will facilitate the discovery and development of novel polyketide-derived bio-insecticides.

Phylogenomics of a new fungal phylum reveals multiple waves of reductive evolution across Holomycota

Compared to multicellular fungi and unicellular yeasts, unicellular fungi with free-living flagellated stages (zoospores) remain poorly known and their phylogenetic position is often unresolved. Recently, rRNA gene phylogenetic analyses of two atypical parasitic fungi with amoeboid zoospores and long kinetosomes, the sanchytrids Amoeboradix gromovi and Sanchytrium tribonematis, showed that they formed a monophyletic group without close affinity with known fungal clades. Here, we sequence single-cell genomes for both species to assess their phylogenetic position and evolution. Phylogenomic analyses using different protein datasets and a comprehensive taxon sampling result in an almost fully-resolved fungal tree, with Chytridiomycota as sister to all other fungi, and sanchytrids forming a well-supported, fast-evolving clade sister to Blastocladiomycota. Comparative genomic analyses across fungi and their allies (Holomycota) reveal an atypically reduced metabolic repertoire for sanchytrids. We infer three main independent flagellum losses from the distribution of over 60 flagellum-specific proteins across Holomycota. Based on sanchytrids' phylogenetic position and unique traits, we propose the designation of a novel phylum, Sanchytriomycota. In addition, our results indicate that most of the hyphal morphogenesis gene repertoire of multicellular fungi had already evolved in early holomycotan lineages.

In vitro transcriptomes analysis identifies some special genes involved in pathogenicity difference of the Beauveria bassiana against different insect hosts

Typical entomopathogenic filamentous fungi such as Beauveria bassiana infect susceptible hosts via penetration of insect cuticle. The pathogenicity of B. bassiana strain to diverse insect hosts is different. While the molecular mechanisms of B. bassiana adapt to different insects are not well clear. B. bassiana GXsk1011 is a hyper-virulent strain from silkworm, which was investigated on the metabolic responses to three cuticle extracts of Bombyx mori, Helicoverpa armigera and Clanis bilineata at 24 h by RNA-seq method. A total of 638 up- and 400 down-regulated differentially expressed genes (DEGs) were identified in B. bassiana grown on H. armigera compared with B. mori, and 910 up- and 401 down-regulated genes for C. bilineata compared with B. mori. Functional categorization showed that DEGs are mainly involved in metabolic processes, localization, catalytic activity and transporter activity. Analysis of 20 highest fold change genes in DEGs showed that when B. bassiana transferred to non-original hosts as H. armigera and C. bilineata, the adhesion (Mad1), protease (Pr2) and cell surface protein (BBA_09174), etc. were down-regulated. While the class III chitinase ChiA2 (BBA_05353, Bbchi-17), major allergen Asp f 2-like protein (BBA_05395, Bb-f2) and nonribosomal peptide synthase, etc. were up-regulated. The secretory lipase that responded to H. armigera and the phosphate permease responded to C. bilineata were also up-regulated in the Top 20 DEGs. These special expressed genes indicate when the B. bassiana transferred to non-original hosts (or called as non-natural hosts), the strain appeared the changes of metabolic response and infection strategies to adapt to new hosts, and implied the key actions of infected adaptation were to break the barrier of different cuticle chitin component and against the immune stress of hosts. This study provided an insight into the B. bassiana that with wide host ranges how to adapt to infect different insect hosts, which will help us to further understand the pathogenesis of B. bassiana infection.

Metagenomics survey unravels diversity of biogas microbiomes with potential to enhance productivity in Kenya

The obstacle to optimal utilization of biogas technology is poor understanding of biogas microbiomes diversities over a wide geographical coverage. We performed random shotgun sequencing on twelve environmental samples. Randomized complete block design was utilized to assign the twelve treatments to four blocks, within eastern and central regions of Kenya. We obtained 42 million paired-end reads that were annotated against sixteen reference databases using two ENVO ontologies, prior to β-diversity studies. We identified 37 phyla, 65 classes and 132 orders. Bacteria dominated and comprised 28 phyla, 42 classes and 92 orders, conveying substrate's versatility in the treatments. Though, Fungi and Archaea comprised 5 phyla, the Fungi were richer; suggesting the importance of hydrolysis and fermentation in biogas production. High β-diversity within the taxa was largely linked to communities' metabolic capabilities. Clostridiales and Bacteroidales, the most prevalent guilds, metabolize organic macromolecules. The identified Cytophagales, Alteromonadales, Flavobacteriales, Fusobacteriales, Deferribacterales, Elusimicrobiales, Chlamydiales, Synergistales to mention but few, also catabolize macromolecules into smaller substrates to conserve energy. Furthermore, δ-Proteobacteria, Gloeobacteria and Clostridia affiliates syntrophically regulate PH2 and reduce metal to provide reducing equivalents. Methanomicrobiales and other Methanomicrobia species were the most prevalence Archaea, converting formate, CO2(g), acetate and methylated substrates into CH4(g). Thermococci, Thermoplasmata and Thermoprotei were among the sulfur and other metal reducing Archaea that contributed to redox balancing and other metabolism within treatments. Eukaryotes, mainly fungi were the least abundant guild, comprising largely Ascomycota and Basidiomycota species. Chytridiomycetes, Blastocladiomycetes and Mortierellomycetes were among the rare species, suggesting their metabolic and substrates limitations. Generally, we observed that environmental and treatment perturbations influenced communities' abundance, β-diversity and reactor performance largely through stochastic effect. Understanding diversity of biogas microbiomes over wide environmental variables and its' productivity provided insights into better management strategies that ameliorate biochemical limitations to effective biogas production.

Metarhizium: jack of all trades, master of many

The genus Metarhizium and Pochonia chlamydosporia comprise a monophyletic clade of highly abundant globally distributed fungi that can transition between long-term beneficial associations with plants to transitory pathogenic associations with frequently encountered protozoans, nematodes or insects. Some very common 'specialist generalist' species are adapted to particular soil and plant ecologies, but can overpower a wide spectrum of insects with numerous enzymes and toxins that result from extensive gene duplications made possible by loss of meiosis and associated genome defence mechanisms. These species use parasexuality instead of sex to combine beneficial mutations from separate clonal individuals into one genome (Vicar of Bray dynamics). More weakly endophytic species which kill a narrow range of insects retain sexuality to facilitate host-pathogen coevolution (Red Queen dynamics). Metarhizium species can fit into numerous environments because they are very flexible at the genetic, physiological and ecological levels, providing tractable models to address how new mechanisms for econutritional heterogeneity, host switching and virulence are acquired and relate to diverse sexual life histories and speciation. Many new molecules and functions have been discovered that underpin Metarhizium associations, and have furthered our understanding of the crucial ecology of these fungi in multiple habitats.

Ecologically controlling insect and mite pests of tea plants with microbial pesticides: a review

Insect and mite pests are damaging stressors that are threatening the cultivation of tea plants, which result in enormous crop loss. Over the years, the effectiveness of synthetic pesticides has allowed for its prominent application as a control strategy. However, the adverse effects of synthetic pesticides in terms of pesticide residue, environmental contamination and insect pest resistance have necessitated the need for alternative strategies. Meanwhile, microbial pesticides have been applied to tackle the damaging activities of the insect and mite pests of tea plants, and their performances were scientifically adjudged appreciable and environmental friendly. Herein, entomopathogenic microbes that were effective against tea geometrid (Ectropis obliqua Prout), tea green leafhopper (Empoasca onukii Matsuda), paraguay tea ampul (Gyropsylla spegazziniana), tea mosquito bug (Helopeltis theivora Waterhouse) and red spider mite (Oligonychus coffea Nietner) have been reviewed. The current findings revealed that microbial pesticides were effective and showed promising performances against these pests. Overall, this review has provided the basic and integrative information on the integrated pest management (IPM) tool(s) that can be utilized towards successful control of the aforementioned insect and mite pests.

Identification and characterization of the glycoside hydrolase family 18 genes from the entomopathogenic fungus Isaria cicadae genome

Fungal chitinases play essential roles in chitin degradation, cell wall remodeling, chitin recycling, nutrition acquisition, autolysis, and virulence. In this study, 18 genes of the glycoside hydrolase 18 (GH18) family were identified in the Isaria cicadae genome. Seventeen of the genes belonged to chitinases and one was an endo-β-N-acetylglucosaminidase (ENGase). According to phylogenetic analysis, the 17 chitinases were designated as subgroups A (7 chitinases), B (7), and C (3). The exon-intron organizations of these genes were analyzed. The conserved regions DxxDxDxE and S/AxGG and the domains CBM1, CBM18, and CBM50 were detected in I. cicadae chitinases and ENGase. The results of analysis of expression patterns showed that genes ICchiA1, ICchiA6, ICchiB1, and ICchiB4 had high transcript levels in the different growth conditions or developmental stages. Subgroup A chitinase genes had higher transcript levels than the genes of all other chitinases. Subgroup B chitinase genes (except ICchiB7) presented higher transcript levels in chitin medium compared with other conditions. ICchiC2 and ICchiC3 were mainly transcribed in autolysis medium and in blastospores, respectively. Moreover, ICchiB1 presented higher transcript levels than genes of other chitinases. This work provides an overview of the GH18 chitinases and ENGase in I. cicadae and provides a context for the chitinolytic potential, functions, and biological controls of these enzymes of entomopathogenic fungi.

Genome Analysis of Hypomyces perniciosus, the Causal Agent of Wet Bubble Disease of Button Mushroom (Agaricus bisporus)

The mycoparasitic fungus Hypomyces perniciosus causes wet bubble disease of mushrooms, particularly Agaricus bisporus. The genome of a highly virulent strain of H. perniciosus HP10 was sequenced and compared to three other fungi from the order Hypocreales that cause disease on A. bisporus. H. perniciosus genome is ~44 Mb, encodes 10,077 genes and enriched with transposable elements up to 25.3%. Phylogenetic analysis revealed that H. perniciosus is closely related to Cladobotryum protrusum and diverged from their common ancestor ~156.7 million years ago. H. perniciosus has few secreted proteins compared to C. protrusum and Trichoderma virens, but significantly expanded protein families of transporters, protein kinases, CAZymes (GH 18), peptidases, cytochrome P450, and SMs that are essential for mycoparasitism and adaptation to harsh environments. This study provides insights into H. perniciosus evolution and pathogenesis and will contribute to the development of effective disease management strategies to control wet bubble disease.

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.

Chromosome level assembly and secondary metabolite potential of the parasitic fungus Cordyceps militaris

Background

Cordyceps militaris is an insect pathogenic fungus that is prized for its use in traditional medicine. This and other entomopathogenic fungi are understudied sources for the discovery of new bioactive molecules. In this study, PacBio SMRT long read sequencing technology was used to sequence the genome of C. militaris with a focus on the genetic potential for secondary metabolite production in the genome assembly of this fungus.

Results

This is first chromosome level assembly of a species in the Cordyceps genera. In this seven chromosome assembly of 33.6 Mba there were 9371 genes identified. Cordyceps militaris was determined to have the MAT 1-1-1 and MAT 1-1-2 mating type genes. Secondary metabolite analysis revealed the potential for at least 36 distinct metabolites from a variety of classes. Three of these gene clusters had homology with clusters producing desmethylbassianin, equisetin and emericellamide that had been studied in other fungi.

Conclusion

Our assembly and analysis has revealed that C. militaris has a wealth of gene clusters for secondary metabolite production distributed among seven chromosomes. The identification of these gene clusters will facilitate the future study and identification of the secondary metabolites produced by this entomopathogenic fungus.

Electronic supplementary material

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

The Insect Pathogens

Fungi are the most common disease-causing agents of insects; aside from playing a crucial role in natural ecosystems, insect-killing fungi are being used as alternatives to chemical insecticides and as resources for biotechnology and pharmaceuticals. Some common experimentally tractable genera, such as Metarhizium spp., exemplify genetic diversity and dispersal because they contain numerous intraspecific variants with distinct environmental and insect host ranges. The availability of tools for molecular genetics and multiple sequenced genomes has made these fungi ideal experimental models for answering basic questions on the genetic and genomic processes behind adaptive phenotypes. For example, comparative genomics of entomopathogenic fungi has shown they exhibit diverse reproductive modes that often determine rates and patterns of genome evolution and are linked as cause or effect with pathogenic strategies. Fungal-insect pathogens represent lifestyle adaptations that evolved numerous times, and there are significant differences in host range and pathogenic strategies between the major groups. However, typically, spores landing on the cuticle produce appressoria and infection pegs that breach the cuticle using mechanical pressure and cuticle-degrading enzymes. Once inside the insect body cavity, fungal pathogens face a potent and comprehensively studied immune defense by which the host attempts to eliminate or reduce an infection. The Fungal Kingdom stands alone in the range, extent, and complexity of their manipulation of arthropod behavior. In part, this is because most only sporulate on cadavers, so they must ensure the dying host positions itself to allow efficient transmission.

Draft genome of the fungus-growing termite pathogenic fungus Ophiocordyceps bispora (Ophiocordycipitaceae, Hypocreales, Ascomycota)

This article documents the public availability of genome sequence data and assembled contigs representing the partial draft genome of Ophiocordyceps bispora. As one of the few known pathogens of fungus-farming termites, a draft genome of O. bispora represents the opportunity to further the understanding of disease and resistance in these complex termite societies. With the ongoing attempts to resolve the taxonomy of the Hypocralaean family, more genetic data will also help to shed light on the phylogenetic relationship between sexual and asexual life stages. Next generation sequence data is available from the European Nucleotide Archive (ENA) under accession PRJEB13655; run numbers: ERR1368522, ERR1368523, and ERR1368524. Genome assembly available from ENA under accession numbers: FKNF01000001-FKNF01000302. Gene prediction available as protein fasta, nucleotide fasta and GFF file from Mendeley Data with accession doi:10.17632/r99fd6g3s4.2 (http://dx.doi.org/10.17632/r99fd6g3s4.2).

Insect Pathogenic Fungi: Genomics, Molecular Interactions, and Genetic Improvements

Entomopathogenic fungi play a pivotal role in the regulation of insect populations in nature, and representative species have been developed as promising environmentally friendly mycoinsecticides. Recent advances in the genome biology of insect pathogenic fungi have revealed genomic features associated with fungal adaptation to insect hosts and different host ranges, as well as the evolutionary relationships between insect and noninsect pathogens. By using species in the Beauveria and Metarhizium genera as models, molecular biology studies have revealed the genes that function in fungus-insect interactions and thereby contribute to fungal virulence. Taken together with efforts toward genetic improvement of fungal virulence and stress resistance, knowledge of entomopathogenic fungi will potentiate cost-effective applications of mycoinsecticides for pest control in the field. Relative to our advanced insights into the mechanisms of fungal pathogenesis in plants and humans, future studies will be necessary to unravel the gene-for-gene relationships in fungus-insect interactive models.

PHI-base: a new interface and further additions for the multi-species pathogen-host interactions database

The pathogen-host interactions database (PHI-base) is available at www.phi-base.org PHI-base contains expertly curated molecular and biological information on genes proven to affect the outcome of pathogen-host interactions reported in peer reviewed research articles. In addition, literature that indicates specific gene alterations that did not affect the disease interaction phenotype are curated to provide complete datasets for comparative purposes. Viruses are not included. Here we describe a revised PHI-base Version 4 data platform with improved search, filtering and extended data display functions. A PHIB-BLAST search function is provided and a link to PHI-Canto, a tool for authors to directly curate their own published data into PHI-base. The new release of PHI-base Version 4.2 (October 2016) has an increased data content containing information from 2219 manually curated references. The data provide information on 4460 genes from 264 pathogens tested on 176 hosts in 8046 interactions. Prokaryotic and eukaryotic pathogens are represented in almost equal numbers. Host species belong ∼70% to plants and 30% to other species of medical and/or environmental importance. Additional data types included into PHI-base 4 are the direct targets of pathogen effector proteins in experimental and natural host organisms. The curation problems encountered and the future directions of the PHI-base project are briefly discussed.

Glycoside hydrolases family 20 (GH20) represent putative virulence factors that are shared by animal pathogenic oomycetes, but are absent in phytopathogens

BACKGROUND

Although interest in animal pathogenic oomycetes is increasing, the molecular basis mediating oomycete-animal relationships remains virtually unknown. Crinkler (CRN) genes, which have been traditionally associated with the cytotoxic activity displayed by plant pathogenic oomycetes, were recently detected in transcriptome sequences from the entomopathogenic oomycete Lagenidium giganteum, suggesting that these genes may represent virulence factors conserved in both animal and plant pathogenic oomycetes. In order to further characterize the L. giganteum pathogenome, an on-going genomic survey was mined to reveal novel putative virulence factors, including canonical oomycete effectors Crinkler 13 (CRN13) orthologs. These novel sequences provided a basis to initiate gene expression analyses and determine if the proposed L. giganteum virulence factors are differentially expressed in the presence of mosquito larvae (Aedes aegypti).

RESULTS

Sequence analyses revealed that L. giganteum express CRN13 transcripts. The predicted proteins, like other L. giganteum CRNs, contained a conserved LYLA motif at the N terminal, but did not display signal peptides. In contrast, other potential virulence factors, such as Glycoside Hydrolases family 20 (hexosaminidase) and 37 (trehalase) proteins (GH20 and GH37), contained identifiable signal peptides. Genome mining demonstrated that GH20 genes are absent from phytopathogenic oomycete genomes, and that the L. giganteum GH20 sequence is the only reported peronosporalean GH20 gene. All other oomycete GH20 homologs were retrieved from animal pathogenic, saprolegnialean genomes. Furthermore, phylogenetic analyses demonstrated that saprolegnialean and peronosporalean GH20 protein sequences clustered in unrelated clades. The saprolegnialean GH20 sequences appeared as a strongly supported, monophyletic group nested within an arthropod-specific clade, suggesting that this gene was acquired via a lateral gene transfer event from an insect or crustacean genome. In contrast, the L. giganteum GH20 protein sequence appeared as a sister taxon to a plant-specific clade that included exochitinases with demonstrated insecticidal activities. Finally, gene expression analyses demonstrated that the L. giganteum GH20 gene expression level is significantly modulated in the presence of mosquito larvae. In agreement with the protein secretion predictions, CRN transcripts did not show any differential expression.

CONCLUSIONS

These results identified GH20 enzymes, and not CRNs, as potential pathogenicity factors shared by all animal pathogenic oomycetes.

Genome sequence and comparative analysis of clavicipitaceous insect-pathogenic fungus Aschersonia badia with Metarhizium spp

Background

Aschersonia badia [(Ab) Teleomorph: Hypocrella siamensis] is an entomopathogenic fungus that specifically infects scale insects and whiteflies. We present the whole genome sequence of Ab and its comparison with two clavicipitaceous fungi Metarhizium robertsii (MR: generalist entomopathogen) and M. acridum (MAC: acridid-specific entomopathogen) that exhibit variable host preferences. Here, through comparative analysis of pathogen-host interacting genes, carbohydrate active enzymes, secondary metabolite biosynthesis genes, and sexuality genes, we explore the proteins with possible virulence functions in clavicipitaceous fungi. Comprehensive overview of GH18 family chitinases has been provided to decipher the role of chitinases in claviceptaceous fungi that are either host specific or generalists.

Results

We report the 28.8 Mb draft genome of Ab and its comparative genome analysis with MR and MAC. The comparative analyses suggests expansion in pathogen-host interacting gene families and carbohydrate active enzyme families in MR, whilst their contraction in Ab and MAC genomes. The multi-modular NRPS gene (dtxS1) responsible for biosynthesis of the secondary metabolite destruxin in MR is not conserved in Ab, similar to the specialist pathogen MAC. An additional siderophore biosynthetic gene responsible for acquisition of iron was identified in MR. Further, the domain survey of chitinases suggest that the CBM50 (LysM) domains, which participate in chitin-binding functions, were not observed in MAC, but were present in Ab and MR. However, apparent differences in frequency of CBM50 domains associated with chitinases of Ab and MR was identified, where MR chitinases displayed a higher proportion of associated CBM50 domains than Ab chitinases.

Conclusions

This study suggests differences in distribution of dtxS1 and chitinases in specialists (Ab and MAC) and generalists (MR) fungi. Our analysis also suggests the presence of a siderophore biosynthetic gene in the MR genome which perhaps aids in enhanced virulence potential and host range. The variation in association of CBMs, being higher in generalists (MR) and lower in specialists (Ab and MAC) fungi may further be responsible for the differences in host affiliation.

Electronic supplementary material

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

Divergent and Convergent Evolution of Fungal Pathogenicity

Fungal pathogens of plants and animals have multifarious effects; they cause devastating damages to agricultures, lead to life-threatening diseases in humans, or induce beneficial effects by reducing insect pest populations. Many virulence factors have been determined in different fungal pathogens; however, the molecular determinants contributing to fungal host selection and adaptation are largely unknown. In this study, we sequenced the genomes of seven ascomycete insect pathogens and performed the genome-wide analyses of 33 species of filamentous ascomycete pathogenic fungi that infect insects (12 species), plants (12), and humans (9). Our results revealed that the genomes of plant pathogens encode more proteins and protein families than the insect and human pathogens. Unexpectedly, more common orthologous protein groups are shared between the insect and plant pathogens than between the two animal group pathogens. We also found that the pathogenicity of host-adapted fungi evolved multiple times, and that both divergent and convergent evolutions occurred during pathogen–host cospeciation thus resulting in protein families with similar features in each fungal group. However, the role of phylogenetic relatedness on the evolution of protein families and therefore pathotype formation could not be ruled out due to the effect of common ancestry. The evolutionary correlation analyses led to the identification of different protein families that correlated with alternate pathotypes. Particularly, the effector-like proteins identified in plant and animal pathogens were strongly linked to fungal host adaptation, suggesting the existence of similar gene-for-gene relationships in fungus–animal interactions that has not been established before. These results well advance our understanding of the evolution of fungal pathogenicity and the factors that contribute to fungal pathotype formation.

Advances in Genomics of Entomopathogenic Fungi

Fungi are the commonest pathogens of insects and crucial regulators of insect populations. The rapid advance of genome technologies has revolutionized our understanding of entomopathogenic fungi with multiple Metarhizium spp. sequenced, as well as Beauveria bassiana, Cordyceps militaris, and Ophiocordyceps sinensis among others. Phylogenomic analysis suggests that the ancestors of many of these fungi were plant endophytes or pathogens, with entomopathogenicity being an acquired characteristic. These fungi now occupy a wide range of habitats and hosts, and their genomes have provided a wealth of information on the evolution of virulence-related characteristics, as well as the protein families and genomic structure associated with ecological and econutritional heterogeneity, genome evolution, and host range diversification. In particular, their evolutionary transition from plant pathogens or endophytes to insect pathogens provides a novel perspective on how new functional mechanisms important for host switching and virulence are acquired. Importantly, genomic resources have helped make entomopathogenic fungi ideal model systems for answering basic questions in parasitology, entomology, and speciation. At the same time, identifying the selective forces that act upon entomopathogen fitness traits could underpin both the development of new mycoinsecticides and further our understanding of the natural roles of these fungi in nature. These roles frequently include mutualistic relationships with plants. Genomics has also facilitated the rapid identification of genes encoding biologically useful molecules, with implications for the development of pharmaceuticals and the use of these fungi as bioreactors.

The N-acetylglucosamine catabolic gene cluster in Trichoderma reesei is controlled by the Ndt80-like transcription factor RON1

Chitin is an important structural constituent of fungal cell walls composed of N-acetylglucosamine (GlcNAc) monosaccharides, but catabolism of GlcNAc has not been studied in filamentous fungi so far. In the yeast Candida albicans, the genes encoding the three enzymes responsible for stepwise conversion of GlcNAc to fructose-6-phosphate are clustered. In this work, we analysed GlcNAc catabolism in ascomycete filamentous fungi and found that the respective genes are also clustered in these fungi. In contrast to C. albicans, the cluster often contains a gene for an Ndt80-like transcription factor, which we named RON1 (regulator of N-acetylglucosamine catabolism 1). Further, a gene for a glycoside hydrolase 3 protein related to bacterial N-acetylglucosaminidases can be found in the GlcNAc gene cluster in filamentous fungi. Functional analysis in Trichoderma reesei showed that the transcription factor RON1 is a key activator of the GlcNAc gene cluster and essential for GlcNAc catabolism. Furthermore, we present an evolutionary analysis of Ndt80-like proteins in Ascomycota. All GlcNAc cluster genes, as well as the GlcNAc transporter gene ngt1, and an additional transcriptional regulator gene, csp2, encoding the homolog of Neurospora crassa CSP2/GRHL, were functionally characterised by gene expression analysis and phenotypic characterisation of knockout strains in T. reesei.