The genome sequence of the biocontrol fungus Metarhizium anisopliae and comparative genomics of Metarhizium species

Comparative genomics of Metarhizium brunneum strains V275 and ARSEF 4556: unraveling intraspecies diversity

Entomopathogenic fungi belonging to the Order Hypocreales are renowned for their ability to infect and kill insect hosts, while their endophytic mode of life and the beneficial rhizosphere effects on plant hosts have only been recently recognized. Understanding the molecular mechanisms underlying their different lifestyles could optimize their potential as both biocontrol and biofertilizer agents, as well as the wider appreciation of niche plasticity in fungal ecology. This study describes the comprehensive whole genome sequencing and analysis of one of the most effective entomopathogenic and endophytic EPF strains, Metarhizium brunneum V275 (commercially known as Lalguard Met52), achieved through Nanopore and Illumina reads. Comparative genomics for exploring intraspecies variability and analyses of key gene sets were conducted with a second effective EPF strain, M. brunneum ARSEF 4556. The search for strain- or species-specific genes was extended to M. brunneum strain ARSEF 3297 and other species of genus Metarhizium, to identify molecular mechanisms and putative key genome adaptations associated with mode of life differences. Genome size differed significantly, with M. brunneum V275 having the largest genome amongst M. brunneum strains sequenced to date. Genome analyses revealed an abundance of plant-degrading enzymes, plant colonization-associated genes, and intriguing intraspecies variations regarding their predicted secondary metabolic compounds and the number and localization of Transposable Elements. The potential significance of the differences found between closely related endophytic and entomopathogenic fungi, regarding plant growth-promoting and entomopathogenic abilities, are discussed, enhancing our understanding of their diverse functionalities and putative applications in agriculture and ecology.

Integrative taxonomy of Metarhizium anisopliae species complex, based on phylogenomics combined with morphometrics, metabolomics, and virulence data

Metarhizium anisopliae (Clavicipitaceae, Hypocreales) is a globally distributed entomopathogenic fungus, which has been largely studied and used in agriculture for its potent entomopathogenicity. Since its taxonomic establishment as a member of Metarhizium, many closely related taxa have been described with highly similar morphology (cryptic species). A holotype specimen of M. anisopliae is not extant, and the ex-neotype strain (CBS 130.71) does not form a monophyletic clade with other strains, up to now, recognized as M. anisopliae sensu stricto. In this study, we have conducted an integrative taxonomic treatment of M. anisopliae sensu lato by including the ex-neotype strain of M. anisopliae, other unknown strains from our collections identified as M. anisopliae s. lat., as well as other known species that have been previously delimited as closely related but distinct to M. anisopliae. By including whole-genome sequencing, morphometric analysis, LC-MS based metabolomics, and virulence assays, we have demonstrated that M. anisopliae s. str. should also include M. lepidiotae (synonym), and that M. anisopliae s. str. differentiates from the other species of the complex by its metabolome and less severe entomopathogenicity. New taxa, namely M. hybridum, M. neoanisopliae and M. parapingshaense spp. nov., are proposed. The novel taxa proposed here have strong phylogenomics support, corroborated by fine-scale differences in the length/width of conidia/phialides, while the metabolomics and virulence data still largely overlap. We have also demonstrated via population genomics data the existence of local clonal lineages, particularly the one corresponding to the persistence of a biocontrol candidate strain that has been used in the field application for three years. This study showcases the utility of combining various data sources for accurate delimitation of species within an important group of fungal biocontrol agents against pest insects.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

The Sordariomycetes: an expanding resource with Big Data for mining in evolutionary genomics and transcriptomics

Advances in genomics and transcriptomics accompanying the rapid accumulation of omics data have provided new tools that have transformed and expanded the traditional concepts of model fungi. Evolutionary genomics and transcriptomics have flourished with the use of classical and newer fungal models that facilitate the study of diverse topics encompassing fungal biology and development. Technological advances have also created the opportunity to obtain and mine large datasets. One such continuously growing dataset is that of the Sordariomycetes, which exhibit a richness of species, ecological diversity, economic importance, and a profound research history on amenable models. Currently, 3,574 species of this class have been sequenced, comprising nearly one-third of the available ascomycete genomes. Among these genomes, multiple representatives of the model genera Fusarium, Neurospora, and Trichoderma are present. In this review, we examine recently published studies and data on the Sordariomycetes that have contributed novel insights to the field of fungal evolution via integrative analyses of the genetic, pathogenic, and other biological characteristics of the fungi. Some of these studies applied ancestral state analysis of gene expression among divergent lineages to infer regulatory network models, identify key genetic elements in fungal sexual development, and investigate the regulation of conidial germination and secondary metabolism. Such multispecies investigations address challenges in the study of fungal evolutionary genomics derived from studies that are often based on limited model genomes and that primarily focus on the aspects of biology driven by knowledge drawn from a few model species. Rapidly accumulating information and expanding capabilities for systems biological analysis of Big Data are setting the stage for the expansion of the concept of model systems from unitary taxonomic species/genera to inclusive clusters of well-studied models that can facilitate both the in-depth study of specific lineages and also investigation of trait diversity across lineages. The Sordariomycetes class, in particular, offers abundant omics data and a large and active global research community. As such, the Sordariomycetes can form a core omics clade, providing a blueprint for the expansion of our knowledge of evolution at the genomic scale in the exciting era of Big Data and artificial intelligence, and serving as a reference for the future analysis of different taxonomic levels within the fungal kingdom.

Combatting insects mediated biotic stress through plant associated endophytic entomopathogenic fungi in horticultural crops

Horticultural production is a vital catalyst for economic growth, yet insect infestations reduce horticultural crop yield and quality. Pesticides and other pest control methods are used during planting to eliminate pests that cause direct and indirect losses. In such situations, endophytic entomo-pathogenic fungi (EEPF) can act as a potential tools for biological control. They protect plants by boosting growth, nutrition, morpho-physiology and salt or iron tolerance. Antixenosis, antibiosis and plant tolerance change insect performance and preferences. EEPF- plant colonisation slows herbivore development, food consumption, oviposition and larval survival. EEPF changes plant physio-chemical properties like volatile emission profile and secondary metabolite production to regulate insect pest defences. EEPF produces chitinases, laccases, amylases, and cellulases for plant defence. Recent studies focused on EEPF species' significance, isolation, identification and field application. Realizing their full potential is difficult due to insufficient mass production, storage stability and formulation. Genetic-molecular and bioinformatics can help to build EEPF-based biological control systems. Metagenomics helps study microbial EEPF taxonomy and function. Multi-omics and system biology can decode EEPF interactions with host plants and microorganisms. NGS (Next Generation Sequencing), comparative genomics, proteomics, transcriptomics, metabolomics, metatranscriptomics and microarrays are used to evaluate plant-EEPF relationships. IPM requires understanding the abiotic and biotic elements that influence plant-EEPF interaction and the physiological mechanisms of EEPF colonisation. Due to restricted research, there are hundreds of unexplored EEPFs, providing an urgent need to uncover and analyse them.

Comparative analysis of Penicillium genomes reveals the absence of a specific genetic basis for biocontrol in Penicillium rubens strain 212

Penicillium rubens strain 212 (PO212) is a filamentous fungus belonging to the division Ascomycete. PO212 acts as an effective biocontrol agent against several pathogens in a variety of horticultural crops including Fusarium oxysporum f.sp. lycopersici, causing vascular wilt disease in tomato plants. We assembled draft genomes of two P. rubens strains, the biocontrol agent PO212 and the soil isolate S27, which lacks biocontrol activity. We also performed comparative analyses of the genomic sequence of PO212 with that of the other P. rubens and P. chrysogenum strains. This is the first Penicillium strain with biocontrol activity whose genome has been sequenced and compared. PO212 genome size is 2,982 Mb, which is currently organized into 65 scaffolds and a total of 10,164 predicted Open Reading Frames (ORFs). Sequencing confirmed that PO212 belongs to P. rubens clade. The comparative analysis of the PO212 genome with the genomes of other P. rubens and Penicillium chrysogenum strains available in databases showed strong conservation among genomes, but a correlation was not found between these genomic data and the biocontrol phenotype displayed by PO212. Finally, the comparative analysis between PO212 and S27 genomes showed high sequence conservation and a low number of variations mainly located in ORF regions. These differences found in coding regions between PO212 and S27 genomes can explain neither the biocontrol activity of PO212 nor the absence of such activity in S27, opening a possible avenue toward transcriptomic and epigenetic studies that may shed light on this mechanism for fighting plant diseases caused by fungal pathogens. The genome sequences described in this study provide a useful novel resource for future research into the biology, ecology, and evolution of biological control agents.

High susceptibility of Tetranychus merganser (Acari: Tetranychidae), an emergent pest of the tropical crop Carica papaya, towards Metarhizium anisopliae s.l. and Beauveria bassiana strains

Background

The mite Tetranychus merganser is considered to be an emerging pest of various crops in tropical countries. It is one of the most detrimental pests in the papaya orchards of some regions of México. The current field control of Tetranychus spp. involves the extensive use of chemicals that have some degree of toxicity to humans or the environment and may cause selective resistance. The use of biological alternatives such as parasitoids and mite predators have limited effectiveness. In order to find effective but non-toxic alternatives for mite pest management, bio-products that are able to be mass produced and applied to large production areas have been sought, including the entomopathogen fungi. B. bassiana and M. anisopliae s.l. are the fungi most extensively used for the biological control of insect pests. Although they do not cause natural epizootic diseases in mites, there are reports that show that they infect T. urticae, and should be evaluated for use in the biological control of papaya's mite pests.

Methods

A T. merganser colony was established and the susceptibility of adult females to 30 entomopathogenic fungi strains was evaluated under laboratory conditions with an in vitro mass screening bioassay. Ten strains of Metarhizium anisopliae sensu lato (s.l.), eleven of Beauveria bassiana, nine of Lecanicillium sp. and one of Hirsutella thompsonii var. sinematosa were tested. The infectivity of adult females was evaluated calculating the percentage of mortality. To calculate the LC₅₀ and LT₅₀ of the most virulent strains, a bioassay was performed using serial concentrations (1×10⁴-1×10⁸ conidia/mL) for each strain. Strains showing ability to infect eggs laid were evaluated with a novel egg-infectivity bioassay. The internal transcribed spacer (ITS) region of the more lethal strains were sequenced.

Results

T. merganser and T. urticae were found in orchards of Carica papaya (Maradol variety and Tainung hybrid) in Campeche, México. All tested strains of M. anisopliae s.l. and B. bassiana were infectious to the adult female of T. merganser at a concentration of 1×10⁸ conidia/mL. Six strains of M. anisopliae (Ma002, Ma003, Ma004, Ma005, Ma014 and Ma034) caused 100% mortality, and one of B. bassiana (Bb016) caused 95% mortality. The most virulent was Ma034, with an LC₅₀ of 1.73×10⁶ conidia/mL followed by Ma005 and Ma003. Ma005 and Ma034 were the fastest strains to reach LT_50,achieving this in less than 3.7 days. Additionally, Ma034 and Ma014 strains were infectious to more than 70% of the eggs.

Conclusions

T. merganser and T. urticae are present in the papaya orchards of Campeche, México. The high susceptibility of T. merganser adult females and eggs toward several M. anisopliae s.l. or B. bassiana strains suggests that these fungi are a viable alternative to control this emergent pest. The most virulent strain, Ma034, was also infective to eggs, and is the most promising to be tested in the field.

Genomic and Experimental Analysis of the Insecticidal Factors Secreted by the Entomopathogenic Fungus Beauveria pseudobassiana RGM 2184

The entomopathogenic fungus Beauveria pseudobassiana strain RGM 2184 can reach a maximum efficacy of 80% against the quarantine pest Lobesia botrana in field assays. In this study, the RGM 2184 genome was sequenced, and genome mining analyses were performed to predict the factors involved in its insecticidal activity. Additionally, the metabolic profiling of the RMG 2184 culture's supernatants was analyzed by mass spectrometry, and the insecticidal activity from one of these extracts was evaluated in Galleria mellonella larvae. The genome analysis resulted in 114 genes encoding for extracellular enzymes, four biosynthetic gene clusters reported as producers of insecticidal and bactericidal factors (oosporein, beauvericin, desmethylbassianin, and beauveriolide), 20 toxins, and at least 40 undescribed potential biocontrol factors (polyketides and nonribosomal peptides). Comparative genomic analysis revealed that 65-95% of these genes are Beauveria genus-specific. Metabolic profiling of supernatant extracts from RGM 2184 cultures exhibited secondary metabolites such as beauveriolide, oosporein, inflatin C, and bassiatin. However, a number of detected metabolites still remain undescribed. The metabolite extract caused 79% mortality of Galleria mellonella larvae at 28 days. The results of this research lay the groundwork for the study of new insecticidal molecules.

Increasing Anthonomus grandis susceptibility to Metarhizium anisopliae through RNAi-induced AgraRelish knockdown: a perspective to combine biocontrol and biotechnology

BACKGROUND

The hemolymph and insect gut together have an essential role in the immune defense against microorganisms, including the production of antimicrobial peptides (AMP). AMPs are mainly induced by two specific signaling pathways, Toll and immune deficiency (IMD). Here, we characterize the expression profile of four genes from both pathways and describe the importance of AgraRelish in the immune defense of Anthonomus grandis against the entomopathogenic fungus Metarhizium anisopliae by RNA interference (RNAi).

RESULTS

To characterize the pathway that is activated early during the A. grandis-M. anisopliae interaction, we assessed the expression profiles of AgraMyD88 and AgraDorsal (Toll pathway), AgraIMD and AgraRelish (IMD pathway), and several AMP genes. Interestingly, we found that IMD pathway genes are upregulated early, and Toll pathway genes are upregulated just 3 days after inoculation (DAI). Furthermore, nine AMPs were upregulated 24 h after fungus inoculation, including attacins, cecropins, coleoptericins, and defensins. AgraRelish knockdown resulted in a reduction in median lethal time (LT₅₀ ) for M. anisopliae-treated insects of around 2 days compared to control treatments. In addition, AgraRelish remained knocked down at 3 DAI. Finally, we identified that AgraRelish knockdown increased fungal loads at 2 DAI compared to control treatments, possibly indicating a faster infection.

CONCLUSIONS

Our data indicate the influence of the IMD pathway on the antifungal response in A. grandis. Combining biocontrol and RNAi could significantly improve cotton boll weevil management. Hence, AgraRelish is a potential target for the development of biotechnological tools aimed at improving the efficacy of M. anisopliae against A. grandis.

Telomere length de novo assembly of all 7 chromosomes and mitogenome sequencing of the model entomopathogenic fungus, Metarhizium brunneum, by means of a novel assembly pipeline

Background

More accurate and complete reference genomes have improved understanding of gene function, biology, and evolutionary mechanisms. Hybrid genome assembly approaches leverage benefits of both long, relatively error-prone reads from third-generation sequencing technologies and short, accurate reads from second-generation sequencing technologies, to produce more accurate and contiguous de novo genome assemblies in comparison to using either technology independently. In this study, we present a novel hybrid assembly pipeline that allowed for both mitogenome de novo assembly and telomere length de novo assembly of all 7 chromosomes of the model entomopathogenic fungus, Metarhizium brunneum.

Results

The improved assembly allowed for better ab initio gene prediction and a more BUSCO complete proteome set has been generated in comparison to the eight current NCBI reference Metarhizium spp. genomes. Remarkably, we note that including the mitogenome in ab initio gene prediction training improved overall gene prediction. The assembly was further validated by comparing contig assembly agreement across various assemblers, assessing the assembly performance of each tool. Genomic synteny and orthologous protein clusters were compared between Metarhizium brunneum and three other Hypocreales species with complete genomes, identifying core proteins, and listing orthologous protein clusters shared uniquely between the two entomopathogenic fungal species, so as to further facilitate the understanding of molecular mechanisms underpinning fungal-insect pathogenesis.

Conclusions

The novel assembly pipeline may be used for other haploid fungal species, facilitating the need to produce high-quality reference fungal genomes, leading to better understanding of fungal genomic evolution, chromosome structuring and gene regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07390-y.

Virulence Factors of the Entomopathogenic Genus Metarhizium

The fungal genus Metarhizium has been used as an entomopathogen worldwide for approximately 140 years, and its mechanism of infection and its virulence factors have been studied. The present review is a compilation of virulence factors described in the literature to date and their participation in specific stages of the infection process.

An Overview of Genomics, Phylogenomics and Proteomics Approaches in Ascomycota

Fungi are among the most successful eukaryotes on Earth: they have evolved strategies to survive in the most diverse environments and stressful conditions and have been selected and exploited for multiple aims by humans. The characteristic features intrinsic of Fungi have required evolutionary changes and adaptations at deep molecular levels. Omics approaches, nowadays including genomics, metagenomics, phylogenomics, transcriptomics, metabolomics, and proteomics have enormously advanced the way to understand fungal diversity at diverse taxonomic levels, under changeable conditions and in still under-investigated environments. These approaches can be applied both on environmental communities and on individual organisms, either in nature or in axenic culture and have led the traditional morphology-based fungal systematic to increasingly implement molecular-based approaches. The advent of next-generation sequencing technologies was key to boost advances in fungal genomics and proteomics research. Much effort has also been directed towards the development of methodologies for optimal genomic DNA and protein extraction and separation. To date, the amount of proteomics investigations in Ascomycetes exceeds those carried out in any other fungal group. This is primarily due to the preponderance of their involvement in plant and animal diseases and multiple industrial applications, and therefore the need to understand the biological basis of the infectious process to develop mechanisms for biologic control, as well as to detect key proteins with roles in stress survival. Here we chose to present an overview as much comprehensive as possible of the major advances, mainly of the past decade, in the fields of genomics (including phylogenomics) and proteomics of Ascomycota, focusing particularly on those reporting on opportunistic pathogenic, extremophilic, polyextremotolerant and lichenized fungi. We also present a review of the mostly used genome sequencing technologies and methods for DNA sequence and protein analyses applied so far for fungi.

Several Metarhizium Species Produce Ergot Alkaloids in a Condition-Specific Manner

Genomic sequence data indicate that certain fungi in the genus Metarhizium have the capacity to produce lysergic acid-derived ergot alkaloids, but accumulation of ergot alkaloids in these fungi has not been demonstrated previously. We assayed several Metarhizium species grown under different conditions for accumulation of ergot alkaloids. Isolates of M. brunneum and M. anisopliae accumulated the lysergic acid amides lysergic acid α-hydroxyethyl amide, ergine, and ergonovine on sucrose-yeast extract agar but not on two other tested media. Isolates of six other Metarhizium species did not accumulate ergot alkaloids on sucrose-yeast extract agar. Conidia of M. brunneum lacked detectable ergot alkaloids, and mycelia of this fungus secreted over 80% of their ergot alkaloid yield into the culture medium. Isolates of M. brunneum, M. flavoviride, M. robertsii, M. acridum, and M. anisopliae produced high concentrations of ergot alkaloids in infected larvae of the model insect Galleria mellonella, but larvae infected with M. pingshaense, M. album, M. majus, and M. guizhouense lacked detectable ergot alkaloids. Alkaloid concentrations were significantly higher when insects were alive (as opposed to killed by freezing or gas) at the time of inoculation with M. brunneum Roots of corn and beans were inoculated with M. brunneum or M. flavoviride and global metabolomic analyses indicated that the inoculated roots were colonized, though no ergot alkaloids were detected. The data demonstrate that several Metarhizium species produce ergot alkaloids of the lysergic acid amide class and that production of ergot alkaloids is tightly regulated and associated with insect colonization.IMPORTANCE Our discovery of ergot alkaloids in fungi of the genus Metarhizium has agricultural and pharmaceutical implications. Ergot alkaloids produced by other fungi in the family Clavicipitaceae accumulate in forage grasses or grain crops; in this context they are considered toxins, though their presence also may deter or kill insect pests. Our data report ergot alkaloids in Metarhizium species and indicate a close association of ergot alkaloid accumulation with insect colonization. The lack of accumulation of alkaloids in spores of the fungi and in plants colonized by the fungi affirms the safety of using Metarhizium species as biocontrol agents. Ergot alkaloids produced by other fungi have been exploited to produce powerful pharmaceuticals. The class of ergot alkaloids discovered in Metarhizium species (lysergic acid amides) and their secretion into the growth medium make Metarhizium species a potential platform for future studies on ergot alkaloid synthesis and modification.

Meta-Analysis and Evaluation by Insect-Mediated Baiting Reveal Different Patterns of Hypocrealean Entomopathogenic Fungi in the Soils From Two Regions of China

A survey was carried out on forest soils and grassland soils from Hebei and Sichuan provinces using Tenebrio molitor larvae as a bait, and high-throughput DNA sequencing (HTS) of the fungal internal transcribed spacer-2 ribosomal DNA was used to monitor the natural distribution of three leading hypocrealean families of insect fungal pathogens (Clavicipitaceae, Cordycipitaceae, and Ophiocordycipitaceae). The occurrence of insect fungal pathogens in soil samples from 98 different sites was compared. The use of insect bait indicated that entomopathogenic fungi of the genus Metarhizium were predominant, followed by Beauveria and Isaria. Molecular characterization using the Mz_FG543 intergenic region revealed that the Metarhizium species pool was phylogenetically composed of three closely related species as follows; Metarhizium pingshaense (n = 74), Metarhizium robertsii (n = 51), and Metarhizium brunneum (n = 26), as well as one isolate which clustered with Metarhizium flavoviride. Nine potentially new phylogenetic species were delimited within the M. flavoviride species complex by sequencing of the 5′ elongation factor-1 alpha region locus. The Beauveria (n = 64) and Isaria (n = 5) isolates were characterized via sequence analyses of the Bloc region. An intergenic spacer phylogeny of the Beauveria isolate assemblage revealed the phylogenetic species within the Beauveria bassiana clade. Interestingly, the individuals of M. pingshaense (n = 18) and M. brunneum (n = 12) exhibited the presence of both mating types in Sichuan Province. Similarly, for the Beauveria isolates, reproductive mode assays demonstrated that all four B. bassiana subclades possessed bipolar outcrossing mating systems. Of these, 19 isolates contained two mating types, and the rest were fixed for single mating types, revealing opportunities for intra-lineage heterothallic mating. The HTS results showed a significantly higher occurrence of the Clavicipitaceae family and the Metarhizium genus in the soil samples. The Venn diagram showed Metarhizium anisopliae (senso lato), Isaria farinose, and B. bassiana as frequently abundant fungal pathogen operational taxonomic units (core) across sampling sites, while the baiting method showed that the genus of Isaria was isolated locally. The Mantel test verified that community dissimilarity increased significantly with geographical distance, suggesting that geographical coordinates are possible factors that influence the insect fungal pathogen community composition in the studied sites. This study is the first to highlight the usefulness of utilizing soil baiting and deep sequencing to investigate the population dynamics of entomopathogens in soil.

Genomic characterization of Parengyodontium americanum sp. nov

Modern genome analysis and phylogenomic methods have increased the number of fungal species, as well as enhanced appreciation of the degree of diversity within the fungal kingdom. In this context, we describe a new Parengyodontium species, P. americanum, which is phylogenetically related to the opportunistic human fungal pathogen P. album. Five unusual fungal isolates were recovered from five unique and confirmed coccidioidomycosis patients, and these isolates were subsequently submitted to detailed molecular and morphological identification procedures to determine identity. Molecular and morphological diagnostic analyses showed that the isolates belong to the Cordycipitaceae. Subsequently, three representative genomes were sequenced and annotated, and a new species, P. americanum, was identified. Using various genomic analyses, gene family expansions related to novel compounds and potential for ability to grow in diverse habitats are predicted. A general description of the genomic composition of this newly described species and comparison of genome content with Beauveria bassiana, Isaria fumosorosea and Cordyceps militaris shows a shared core genome of 6371 genes, and 148 genes that appear to be specific for P. americanum. This work provides the framework for future investigations of this interesting fungal species.

Host-dependent contributions of the Cfcdp1 protease gene to virulence in the entomopathogenic fungus Cordyceps fumosorosea

BACKGROUND

The genomes of broad host range insect pathogenic fungi, including Cordyceps fumosorosea, encode for a suite of secreted proteases implicated in targeting, penetration, and degradation of the host exoskeleton or cuticle. These cuticle-degrading proteases act as critical virulence factors, but their functions within the biological context, particularly in relation to host specificity, remain poorly characterized.

RESULT

A C. fumosorosea protease gene, Cfcdp1, was identified and a targeted gene-knockout strain constructed. Minor growth defects were observed for the ΔCfcdp1 strain when compared to the wild-type parent and complemented (ΔCfcdp1::Cfcdp1) strains, with delayed and decreased sporulation noted for the mutant. Decreased subtilisin-like protease activity was seen for the ΔCfcdp1 strain, although total secreted protease activity was similar between the mutant and wild-type strains. Insect bioassays using whitefly, Bemisia tabaci, and cabbageworm, Pieris rapae, showed decreased infectivity, i.e. 2.4-3.4-fold increase in lethal dose (LC₅₀ ) and an increased time to death (LT₅₀ ), for the ΔCfcdp1 strain. In contrast, insect bioassays using the diamondback moth, Plutella xylostella, or the brown planthopper, Nilaparvata lugens, showed increased infectivity, i.e. a 3-5-fold decrease in LC₅₀ , and a decreased LT₅₀ . Differential effects were also seen on the fecundity of B. tabaci infected by the different fungal strains.

CONCLUSION

These data reveal host-dependent effects of a protease implicated in cuticle degradation on C. fumosorosea virulence. The implications of these findings in suggesting context-dependent requirements of cuticle-degrading enzymes and their potentially differential roles in mediating virulence towards different hosts are discussed. © 2019 Society of Chemical Industry.

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 Sequencing of Cladobotryum protrusum Provides Insights into the Evolution and Pathogenic Mechanisms of the Cobweb Disease Pathogen on Cultivated Mushroom

Cladobotryum protrusum is one of the mycoparasites that cause cobweb disease on cultivated edible mushrooms. However, the molecular mechanisms of evolution and pathogenesis of C. protrusum on mushrooms are largely unknown. Here, we report a high-quality genome sequence of C. protrusum using the single-molecule, real-time sequencing platform of PacBio and perform a comparative analysis with closely related fungi in the family Hypocreaceae. The C. protrusum genome, the first complete genome to be sequenced in the genus Cladobotryum, is 39.09 Mb long, with an N50 of 4.97 Mb, encoding 11,003 proteins. The phylogenomic analysis confirmed its inclusion in Hypocreaceae, with its evolutionary divergence time estimated to be ~170.1 million years ago. The genome encodes a large and diverse set of genes involved in secreted peptidases, carbohydrate-active enzymes, cytochrome P450 enzymes, pathogen⁻host interactions, mycotoxins, and pigments. Moreover, C. protrusum harbors arrays of genes with the potential to produce bioactive secondary metabolites and stress response-related proteins that are significant for adaptation to hostile environments. Knowledge of the genome will foster a better understanding of the biology of C. protrusum and mycoparasitism in general, as well as help with the development of effective disease control strategies to minimize economic losses from cobweb disease in cultivated edible mushrooms.

Species of the Metarhizium anisopliae complex with diverse ecological niches display different susceptibilities to antifungal agents

Species of the Metarhizium anisopliae complex are globally ubiquitous soil-inhabiting and predominantly insect-pathogenic fungi. The Metarhizium genus contains species ranging from specialists, such as Metarhizium acridum that only infects acridids, to generalists, such as M. anisopliae, Metarhizium brunneum, and Metarhizium robertsii that infect a broad range of insects and can also colonize plant roots. There is little information available about the susceptibility of Metarhizium species to clinical and non-clinical antifungal agents. We determined the susceptibility of 16 isolates comprising four Metarhizium species with different ecological niches to seven clinical (amphotericin B, ciclopirox olamine, fluconazole, griseofulvin, itraconazole, tebinafine, and voriconazole) and one non-clinical (benomyl) antifungal agents. All isolates of the specialist M. acridum were clearly more susceptible to most antifungals than the isolates of the generalists M. anisopliae sensu lato, M. brunneum, and M. robertsii. All isolates of M. anisopliae, M. brunneum, and M. robertsii were resistant to fluconazole and some were also resistant to amphotericin B. The marked differences in susceptibility between the specialist M. acridum and the generalist Metarhizium species suggest that this characteristic is associated with their different ecological niches, and may assist in devising rational antifungal treatments for the rare cases of mycoses caused by Metarhizium species.

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.

Mode of Infection of Metarhizium spp. Fungus and Their Potential as Biological Control Agents

Chemical insecticides have been commonly used to control agricultural pests, termites, and biological vectors such as mosquitoes and ticks. However, the harmful impacts of toxic chemical insecticides on the environment, the development of resistance in pests and vectors towards chemical insecticides, and public concern have driven extensive research for alternatives, especially biological control agents such as fungus and bacteria. In this review, the mode of infection of Metarhizium fungus on both terrestrial and aquatic insect larvae and how these interactions have been widely employed will be outlined. The potential uses of Metarhizium anisopliae and Metarhizium acridum biological control agents and molecular approaches to increase their virulence will be discussed.

Agrobacterium-Mediated Co-transformation of Multiple Genes in Metarhizium robertsii

Fungi of the Metarhizium genus are a very versatile model for understanding pathogenicity in insects and their symbiotic relationship with plants. To establish a co-transformation system for the transformation of multiple M. robertsii genes using Agrobacterium tumefaciens, we evaluated whether the antibiotic nourseothricin has the same marker selection efficiency as phosphinothricin using separate vectors. Subsequently, in the two vectors containing the nourseothricin and phosphinothricin resistance cassettes were inserted eGFP and mCherry expression cassettes, respectively. These new vectors were then introduced independently into A. tumefaciens and used to transform M. robertsii either in independent events or in one single co-transformation event using an equimolar mixture of A. tumefaciens cultures. The number of transformants obtained by co-transformation was similar to that obtained by the individual transformation events. This method provides an additional strategy for the simultaneous insertion of multiple genes into M. robertsii.

Filamentous ascomycete genomes provide insights into Copia retrotransposon diversity in fungi

Background

The relative scarcity of Copia retrotransposons has been recently characterized in metazoans in comparison with the other superfamilies of LTR elements. Furthermore, Copia retrotransposons have often a particular dynamics that results in a highly predominant single clade of elements within a large host taxon, such as the GalEa-like retrotransposons in crustaceans. Taking advantage of the skyrocketing amount of genomic data available for fungi, we carried out the first large-scale comparative genomic analysis of the Copia clades in filamentous ascomycetes.

Results

Screening 30 completely sequenced genomes allowed us to identify more than 2500 Copia copies with conserved LTR, which are distributed in 138 families. Their characterization revealed that fungal Copia diversity is much broader than previously thought with at least 27 clades, 23 of which likely correspond to new ones. While the Copia copy number is low in most species, the two clades GalEa and FunCo1 are widely distributed and highly dominate Copia content as they both account for 80% of the detected sequences.

Conclusions

In Fungi, GalEa retrotransposons are restricted to Pezizomycotina in which they can make up an outstandingly high proportion of the genome (up to 10% in Cenococcum geophilum). At last, we revealed that fungal GalEa elements structurally differ from all other Copia elements with an absence of Primer Binding Site. These elements however harbor a Conserved Hairpin Site which is probably essential for their transposition.

Electronic supplementary material

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

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.

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.

Coevolution of parasitic fungi and insect hosts

Parasitic fungi and their insect hosts provide an intriguing model system for dissecting the complex co-evolutionary processes, which result in Red Queen dynamics. To explore the genetic basis behind host-parasite coevolution we chose two parasitic fungi (Beauveria bassiana and Metarhizium anisopliae, representing the most important entomopathogenic fungi used in the biological control of pest or vector insects) and two established insect model hosts (the greater wax moth Galleria mellonella and the red flour beetle Tribolium castaneum) for which sequenced genomes or comprehensive transcriptomes are available. Focusing on these model organisms, we review the knowledge about the interactions between fungal molecules operating as virulence factors and insect host-derived defense molecules mediating antifungal immunity. Particularly the study of the intimate interactions between fungal proteinases and corresponding host-derived proteinase inhibitors elucidated novel coevolutionary mechanisms such as functional shifts or diversification of involved effector molecules. Complementarily, we compared the outcome of coevolution experiments using the parasitic fungus B. bassiana and two different insect hosts which were initially either susceptible (Galleria mellonella) or resistant (Tribolium castaneum). Taking a snapshot of host-parasite coevolution, we show that parasitic fungi can overcome host barriers such as external antimicrobial secretions just as hosts can build new barriers, both within a relatively short time of coevolution.

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.

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.

Synthesis of silver nanoparticles from Metarhizium robertsii waste biomass extract after nonylphenol degradation, and their antimicrobial and catalytic potential

Waste biomass extracts of Metarhizium robertsii synthesize silver nanoparticles with high biocatalytic and antimicrobial activity.

Nitrogen fertilizer dose alters fungal communities in sugarcane soil and rhizosphere

Fungi play important roles as decomposers, plant symbionts and pathogens in soils. The structure of fungal communities in the rhizosphere is the result of complex interactions among selection factors that may favour beneficial or detrimental relationships. Using culture-independent fungal community profiling, we have investigated the effects of nitrogen fertilizer dosage on fungal communities in soil and rhizosphere of field-grown sugarcane. The results show that the concentration of nitrogen fertilizer strongly modifies the composition but not the taxon richness of fungal communities in soil and rhizosphere. Increased nitrogen fertilizer dosage has a potential negative impact on carbon cycling in soil and promotes fungal genera with known pathogenic traits, uncovering a negative effect of intensive fertilization.

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

Hirsutella thompsonii (Ht) is a fungal pathogen of acarines and the primary cause of epizootics among mites. The draft genomes of two isolates of Ht (MTCC 3556: Ht3, 34.6 Mb and MTCC 6686: Ht6, 34.7 Mb) are presented and compared with the genomes of Beauveria bassiana (Bb) ARSEF 2860 and Ophiocordyceps sinensis (Os) CO18. Comparative analysis of carbohydrate active enzymes, pathogen–host interaction genes, metabolism-associated genes, and genes involved in biosynthesis of secondary metabolites in the four genomes was carried out. Reduction in gene family sizes in Ht3 and Os as compared with Ht6 and Bb is observed. Analysis of the mating type genes in Ht reveals the presence of MAT idiomorphs which is suggestive of cryptic sexual traits in Ht. We further identify and classify putative chitinases that may function as virulence factors in fungal entomopathogens due to their role in degradation of arthropod cuticle.

Comparative genome analysis of entomopathogenic fungi reveals a complex set of secreted proteins

Background

Metarhizium anisopliae is an entomopathogenic fungus used in the biological control of some agricultural insect pests, and efforts are underway to use this fungus in the control of insect-borne human diseases. A large repertoire of proteins must be secreted by M. anisopliae to cope with the various available nutrients as this fungus switches through different lifestyles, i.e., from a saprophytic, to an infectious, to a plant endophytic stage. To further evaluate the predicted secretome of M. anisopliae, we employed genomic and transcriptomic analyses, coupled with phylogenomic analysis, focusing on the identification and characterization of secreted proteins.

Results

We determined the M. anisopliae E6 genome sequence and compared this sequence to other entomopathogenic fungi genomes. A robust pipeline was generated to evaluate the predicted secretomes of M. anisopliae and 15 other filamentous fungi, leading to the identification of a core of secreted proteins. Transcriptomic analysis using the tick Rhipicephalus microplus cuticle as an infection model during two periods of infection (48 and 144 h) allowed the identification of several differentially expressed genes. This analysis concluded that a large proportion of the predicted secretome coding genes contained altered transcript levels in the conditions analyzed in this study. In addition, some specific secreted proteins from Metarhizium have an evolutionary history similar to orthologs found in Beauveria/Cordyceps. This similarity suggests that a set of secreted proteins has evolved to participate in entomopathogenicity.

Conclusions

The data presented represents an important step to the characterization of the role of secreted proteins in the virulence and pathogenicity of M. anisopliae.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-822) contains supplementary material, which is available to authorized users.