Comparative genomics of Beauveria bassiana: uncovering signatures of virulence against mosquitoes

Novel genomic features in entomopathogenic fungus Beauveria bassiana ILB308: accessory genomic regions and putative virulence genes involved in the infection process of soybean pest Piezodorus guildinii

BACKGROUND

Biological control methods involving entomopathogenic fungi like Beauveria bassiana have been shown to be a valuable approach in integrated pest management as an environmentally friendly alternative to control pests and pathogens. Identifying genetic determinants of pathogenicity in B. bassiana is instrumental for enhancing its virulence against insects like the resistant soybean pest Piezodorus guildinii. This study focused on comparative genomics of different B. bassiana strains and gene expression analyses to identify virulence genes in the hypervirulent strain ILB308, especially in response to infection of P. guildinii and growth on hydrocarbon HC15, a known virulence enhancer.

RESULTS

Strain ILB308 showed the highest number of virulence-related features, such as candidate virulence proteins, effectors, small secreted proteins and biosynthetic gene clusters. ILB308 also had a high percentage of unique DNA sequences, including six accessory scaffolds. Gene expression analysis at 4 days post inoculation revealed upregulation of known virulence factors, including Tudor domain proteins, LysM motif-containing proteins, subtilisin-like proteases and novel genes encoding secreted effectors and heat-labile enterotoxins. Growth on HC15 led to the upregulation of genes associated with oxidoreductase activity related to cuticular alkane degradation and fermentation metabolism/antioxidant responses in the hemolymph. The low number of known B. bassiana virulence genes points to novel or unknown mechanisms acting on the interaction between P. guildinii and strain ILB308.

CONCLUSION

The presence of accessory genomic regions and unique virulence genes in ILB308 may contribute to its higher virulence. These genes could be considered as potential targets for enhancing fungal virulence through genetic manipulation. © 2025 Society of Chemical Industry.

Transcriptomic Analysis Reveals Molecular Mechanisms Underpinning Mycovirus-Mediated Hypervirulence in Beauveria bassiana Infecting Tenebrio molitor

Mycoviral infection can either be asymptomatic or have marked effects on fungal hosts, influencing them either positively or negatively. To fully understand the effects of mycovirus infection on the fungal host, transcriptomic profiling of four Beauveria bassiana isolates, including EABb 92/11-Dm that harbors mycoviruses, was performed 48 h following infection of Tenebrio molitor via topical application or injection. Genes that participate in carbohydrate assimilation and transportation, and those essential for fungal survival and oxidative stress tolerance, calcium uptake, and iron uptake, were found to be overexpressed in the virus-infected isolate during the mid-infection stage. Mycotoxin genes encoding bassianolide and oosporein were switched off in all isolates. However, beauvericin, a mycotoxin capable of inducing oxidative stress at the molecular level, was expressed in all four isolates, indicating an important contribution to virulence against T. molitor. These observations suggest that detoxification of immune-related (oxidative) defenses and nutrient scouting, as mediated by these genes, occurs in mid-infection during the internal growth phase. Consequently, we observe a symbiotic relationship between mycovirus and fungus that does not afflict the host; on the contrary, it enhances the expression of key genes leading to a mycovirus-mediated hypervirulence effect.

Computational identification and characterization of chitinase 1 and chitinase 2 from neotropical isolates of Beauveria bassiana

Background

The fungus Beauveria bassiana is widely used for agronomical applications, mainly in biological control. B. bassiana uses chitinase enzymes to degrade chitin, a major chemical component found in insect exoskeletons and fungal cell walls. However, until recently, genomic information on neotropical isolates, as well as their metabolic and biotechnological potential, has been limited.

Methods

Eight complete B. bassiana genomes of Neotropical origin and three references were studied to identify chitinase genes and its corresponding proteins, which were curated and characterized using manual curation and computational tools. We conducted a computational study to highlight functional differences and similarities for chitinase proteins in these Neotropical isolates.

Results

Eleven chitinase 1 genes were identified, categorized as chitinase 1.1 and chitinase 1.2. Five chitinase 2 genes were identified but presented a higher sequence conservation across all sequences. Interestingly, physicochemical parameters were more similar between chitinase 1.1 and chitinase 2 than between chitinase 1.1 and 1.2.

Conclusion

Chitinases 1 and 2 demonstrated variations, especially within chitinase 1, which presented a potential paralog. These differences were observed in their physical parameters. Additionally, CHIT2 completely lacks a signal peptide. This implies that CHIT1 might be associated with infection processes, while CHIT2 could be involved in morphogenesis and cellular growth. Therefore, our work highlights the importance of computational studies on local isolates, providing valuable resources for further experimental validation. Intrinsic changes within local species can significantly impact our understanding of complex pathogen-host interactions and offer practical applications, such as biological control.

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 accumulation of modular serine protease mediated by a novel circRNA sponging miRNA increases Aedes aegypti immunity to fungus

BACKGROUND

Mosquitoes transmit many infectious diseases that affect human health. The fungus Beauveria bassiana is a biological pesticide that is pathogenic to mosquitoes but harmless to the environment.

RESULTS

We found a microRNA (miRNA) that can modulate the antifungal immunity of Aedes aegypti by inhibiting its cognate serine protease. Fungal infection can induce the expression of modular serine protease (ModSP), and ModSP knockdown mosquitoes were more sensitive to B. bassiana infection. The novel miRNA-novel-53 is linked to antifungal immune response and was greatly diminished in infected mosquitoes. The miRNA-novel-53 could bind to the coding sequences of ModSP and impede its expression. Double fluorescence in situ hybridization (FISH) showed that this inhibition occurred in the cytoplasm. The amount of miRNA-novel-53 increased after miRNA agomir injection. This resulted in a significant decrease in ModSP transcript and a significant increase in mortality after fungal infection. An opposite effect was produced after antagomir injection. The miRNA-novel-53 was also knocked out using CRISPR-Cas9, which increased mosquito resistance to the fungus B. bassiana. Moreover, mosquito novel-circ-930 can affect ModSP mRNA by interacting with miRNA-novel-53 during transfection with siRNA or overexpression plasmid.

CONCLUSIONS

Novel-circ-930 affects the expression level of ModSP by a novel-circ-930/miRNA-novel-53/ModSP mechanism to modulate antifungal immunity, revealing new information on innate immunity in insects.

A bacterial-like Pictet-Spenglerase drives the evolution of fungi to produce β-carboline glycosides together with separate genes

Diverse β-carboline (βC) alkaloids are produced by microbes, plants, and animals with myriad bioactivities and drug potentials. However, the biosynthetic mechanism of βCs remains largely elusive, especially regarding the hydroxyl and glucosyl modifications of βCs. Here, we report the presence of the bacterial-like Pictet-Spenglerase gene Fcs1 in the entomopathogenic Beauveria fungi that can catalyze the biosynthesis of the βC skeleton. The overexpression of Fcs1 in Beauveria bassiana led to the identification of six βC methyl glycosides, termed bassicarbosides (BCSs) A-F. We verified that the cytochrome P450 (CYP) genes adjacent to Fcs1 cannot oxidize βCs. Alternatively, the separated CYP684B2 family gene Fcs2 was identified to catalyze βC hydroxylation together with its cofactor gene Fcs3. The functional homologue of Fcs2 is only present in the Fcs1-containing fungi and highly similar to the Fcs1-connected yet nonfunctional CYP. Both evolved quicker than those from fungi without Fcs1 homologues. Finally, the paired methyl/glucosyl transferase genes were verified to mediate the production of BCSs from hydroxy-βCs. All these functionally verified genes are located on different chromosomes of Beauveria, which is in contrast to the typical content-clustered feature of fungal biosynthetic gene clusters (BGCs). We also found that the production of BCSs selectively contributed to fungal infection of different insect species. Our findings shed light on the biosynthetic mechanism of βC glycosides, including the identification of a βC hydroxylase. The results of this study also propose an evolving process of fungal BGC formation following the horizontal transfer of a bacterial gene to fungi.

Genomic virulence features of Beauveria bassiana as a biocontrol agent for the mountain pine beetle population

BACKGROUND

The mountain pine beetle, Dendroctonus ponderosae, is an irruptive bark beetle that causes extensive mortality to many pine species within the forests of western North America. Driven by climate change and wildfire suppression, a recent mountain pine beetle (MPB) outbreak has spread across more than 18 million hectares, including areas to the east of the Rocky Mountains that comprise populations and species of pines not previously affected. Despite its impacts, there are few tactics available to control MPB populations. Beauveria bassiana is an entomopathogenic fungus used as a biological agent in agriculture and forestry and has potential as a management tactic for the mountain pine beetle population. This work investigates the phenotypic and genomic variation between B. bassiana strains to identify optimal strains against a specific insect.

RESULTS

Using comparative genome and transcriptome analyses of eight B. bassiana isolates, we have identified the genetic basis of virulence, which includes oosporein production. Genes unique to the more virulent strains included functions in biosynthesis of mycotoxins, membrane transporters, and transcription factors. Significant differential expression of genes related to virulence, transmembrane transport, and stress response was identified between the different strains, as well as up to nine-fold upregulation of genes involved in the biosynthesis of oosporein. Differential correlation analysis revealed transcription factors that may be involved in regulating oosporein production.

CONCLUSION

This study provides a foundation for the selection and/or engineering of the most effective strain of B. bassiana for the biological control of mountain pine beetle and other insect pests populations.

Male-killer symbiont screening reveals novel associations in Adalia ladybirds

While male-killing bacteria are known to infect across arthropods, ladybird beetles represent a hotspot for these symbioses. In some host species, there are multiple different symbionts that vary in presence and frequency between populations. To further our understanding of spatial and frequency variation, we tested for the presence of three male-killing bacteria: Wolbachia , Rickettsia and Spiroplasma , in two Adalia ladybird species from a previously unexplored UK population. The two-spot ladybird, A. bipunctata, is known to harbour all three male-killers, and we identified Spiroplasma infection in the Merseyside population for the first time. However, in contrast to previous studies on two-spot ladybirds from continental Europe, evidence from egg-hatch rates indicates the Spiroplasma strain present in the Merseyside population does not cause embryonic male-killing. In the related ten-spot ladybird, A. decempunctata, there is only one previous record of a male-killing symbiont, a Rickettsia , which we did not detect in the Merseyside sample. However, PCR assays indicated the presence of a Spiroplasma in a single A. decempunctata specimen. Marker sequence indicated that this Spiroplasma was divergent from that found in sympatric A. bipunctata. Genome sequencing of the Spiroplasma -infected A. decempunctata additionally revealed the presence of cobionts in the form of a Centistes parasitoid wasp and the parasitic fungi Beauveria. Further study of A. decempunctata from this population is needed to resolve whether it is the ladybird or wasp cobiont that harbours Spiroplasma , and to establish the phenotype of this strain. These data indicate first that microbial symbiont phenotype should not be assumed from past studies conducted in different locations, and second that cobiont presence may confound screening studies aimed to detect the frequency of a symbiont in field collected material from a focal host species.

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.

Genomic Characterization and Functional Description of Beauveria bassiana Isolates from Latin America

Beauveria bassiana is an entomopathogenic fungus used in agriculture as a biological controller worldwide. Despite being a well-studied organism, there are no genomic studies of B. bassiana isolates from Central American and Caribbean countries. This work characterized the functional potential of eight Neotropical isolates and provided an overview of their genomic characteristics, targeting genes associated with pathogenicity, the production of secondary metabolites, and the identification of CAZYmes as tools for future biotechnological applications. In addition, a comparison between these isolates and reference genomes was performed. Differences were observed according to geographical location and the lineages of the B. bassiana complex to which each isolate belonged.

Comparative mitogenome analysis reveals mitochondrial genome characteristics in eight strains of Beauveria

Despite the significant progress that has been made in the genome sequencing of Beauveria species, mitochondrial genome (mitogenome) used to examine genetic diversity within fungal populations. Complete mitogenomes of Beauveria species can be easily sequenced and assembled using various sequencing techniques. However, since mitogenome annotations are mainly derived from similar species comparison and software prediction, and are not supported by RNA-seq transcripts data, it leads to problems with the accuracy of mitochondrial annotations and the inability to understand RNA processing. In this study, we assembled and annotated the mitogenome of eight Beauveria strains using Illumina DNA and RNA sequencing data. The circular mitogenome of eight Beauveria strains ranged from 26,850 bp (B. caledonica strain ATCC 64970) to 35,999 bp (B. brongniartii strain GYU-BMZ03), with the intronic insertions accounting for most of the size variation, thus contributing to a total mitochondrial genome (mitogenome) size of 7.01% and 28.95%, respectively. Intron number variations were not directly related to the evolutionary relationship distance. Besides ribosomal protein S3 (rps3), most introns are lost too quickly and lack the stability of protein-coding genes. The short RNA-seq reads from next-generation sequencing can improve the mitochondrial annotation accuracy and help study polycistronic transcripts and RNA processing. The transcription initiation sites may be located in the control region. Most introns do not serve as taxonomic markers and also lack open reading frames (ORFs). We assumed that the poly A tail was added to the polycistronic transcript before splicing and one polycistronic transcript (trnM ₍₁₎-trnL ₍₁₎-trnA-trnF-trnK-trnL ₍₂₎-trnQ-trnH-trnM ₍₂₎-nad2-nad3-atp9-cox2-trnR ₍₁₎-nad4L-nad5-cob-trnC-cox1-trnR ₍₂₎-nad1-nad4-atp8-atp6-rns-trnY-trnD-trnS-trnN-cox3-trnG-nad6-trnV-trnI-trnS-trnW-trnP-rnl(rps3)-trnT-trnE-trnM ₍₃₎) was first processed from the mitogenome and was subsequently processed into smaller mono-, di-, or tricistronic RNAs.

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.

Isolation and identification of two Beauveria bassiana strains from silkworm, Bombyx mori

Silkworm diseases caused by fungi infection occur frequently in sericulture and brought huge economic loss to sericulture. However, on the other hand, some fungi such as Beauveria bassiana, as an important entomological fungus, play an important role in biological control of insect pests. Here, two fungal pathogens causing yellow muscardine were isolated from the silkworm and named as SZY1 and SZY2. These two strains showed almost the same conidial morphology which were smooth, near-spherical, spherical, or ovoid and 2.7 ± 0.6 µm × 2.5 ± 0.9 µm in size, and the hyphal growth rate was also similar. However, the conidia production of SZY2 was almost twice as many as that of SZY1. The complete ribosomal RNA gene was sequenced and analyzed. As a result, the gene sequences of internal transcript space 1 (ITS1)-5.8S rRNA-internal transcript space 2 (ITS2) of SZY1 and SZY2 were identical in sequence and size, and for 18S rRNA, 28S rRNA, and intergenic spacer (IGS), the gene identity of SZY1 to SZY2 was 99%, 99%, and 98%, respectively. Results of phylogenetic analysis based on either ITS1-5.8S rRNA-ITS2 or 18S rRNA showed that both SZY1 and SZY2 were closely related to Beauveria bassiana. These results revealed that the pathogens of yellow muscardine SZY1 and SZY2 were identified as two different strains of Beauveria bassiana, which could provide diagnostic evidence for silkworm muscardine and was helpful for the research and development of novel Bombyx batryticatus and fungal biological insecticide.

Draft Genome Sequences of Eight Isolates of Beauveria bassiana of Neotropical Origin

Beauveria bassiana, a well-known entomopathogenic fungus, has a worldwide distribution; however, genomes of isolates from the Neotropical region are scarce. Here, we report the draft genome sequences of eight B. bassiana isolates from Costa Rica, Puerto Rico, and Honduras.

The Entomopathogenic Fungus Beauveria bassiana Shows Its Toxic Side within Insects: Expression of Genes Encoding Secondary Metabolites during Pathogenesis

Entomopathogenic fungi are extensively used for the control of insect pests worldwide. Among them, Beauveria bassiana (Ascomycota: Hypocreales) produce a plethora of toxic secondary metabolites that either facilitate fungal invasion or act as immunosuppressive compounds. These toxins have different chemical natures, such as nonribosomal peptides and polyketides. Even though their precise role is poorly understood, they are usually linked to virulence. These fungal secondary metabolites are produced by the expression of gene clusters encoding the various proteins needed for their biosynthesis. Each cluster includes synthetases for nonribosomal peptides (NRPS), polyketides (PKS), or hybrid NRPS–PKS genes. The aim of this review is to summarize the information available from transcriptomics and quantitative PCR studies related to the expression of B. bassiana NRPS and PKS genes inside different insects as the infection progresses; as for the host immune response, to help understand the mechanisms that these toxins trigger as virulence factors, antimicrobials, or immunosuppressives within the context of a fungus–insect interaction.

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.

Passage and the evolution of virulence in invertebrate pathogens: Fundamental and applied perspectives

Understanding the ecological and genetic factors that determine the evolution of virulence has broad value for invertebrate pathology. In addition to helping us understand the fundamental biology of our study organisms this body of theory has important applications in microbial biocontrol. Experimental tests of virulence theory are often carried out in invertebrate models and yet theory rarely informs applied passage experiments that aim to increase or maintain virulence. This review summarizes recent progress in this field with a focus on work most relevant to biological control: the virulence of invertebrate pathogens that are 'obligate killers' and which require cadavers for the production of infectious propagules. We discuss recent theory and fundamental and applied experimental evolution with bacteria, fungi, baculoviruses and nematodes. While passage experiments using baculoviruses have a long history of producing isolates with increased virulence, studies with other pathogens have not been so successful. Recent passage experiments that have applied evolution of virulence frameworks based on cooperation (kin selection) have produced novel methods and promising mutants with increased killing power. Evolution of virulence theory can provide plausible explanations for the varied results of passage experiments as well as a predictive framework for improving artificial selection.

Lipidomic response of the entomopathogenic fungus Beauveria bassiana to pyrethroids

Pyrethroids are chemical insecticides that are widely used to control pests. Entomopathogenic fungi are considered environmentally safe alternatives to these compounds. Pyrethroids and entomopathogenic fungi not only co-exist in the environment but can also be applied together in pest control. They are often found in contact with each other, and thus, it seems important to understand their interactions at the cellular level. In this study, we analyzed whether pyrethroids could influence the phospholipid profile of Beauveria bassiana and whether membrane changes are one of the mechanisms by which these fungi adapt to unfavorable environmental conditions. The results of our study revealed that pyrethroids changed the phospholipid profile and increased the cell membrane permeability of B. bassiana, which enabled them to enter and accumulate within the fungal cells, resulting in oxidative stress. Pyrethroids influenced the amount of neutral lipids, caused a decrease in sodium content, and also temporarily lowered the level of the secondary metabolite oosporein in the studied fungi. These findings indicate that the effect of pyrethroids on entomopathogenic fungi may be more complex than originally thought and that lipidomic studies can aid in fully understanding the influence of these chemicals on the mentioned group of fungi.

The integrative taxonomy of Beauveria asiatica and B. bassiana species complexes with whole-genome sequencing, morphometric and chemical analyses

Fungi are rich in complexes of cryptic species that need a combination of different approaches to be delimited, including genomic information. Beauveria (Cordycipitaceae, Hypocreales) is a well-known genus of entomopathogenic fungi, used as a biocontrol agent. In this study we present a polyphasic taxonomy regarding two widely distributed complexes of Beauveria: B. asiatica and B. bassiana s.lat. Some of the genetic groups as previously detected within both taxa were either confirmed or fused using population genomics. High levels of divergence were found between two clades in B. asiatica and among three clades in B. bassiana, supporting their subdivision as distinct species. Morphological examination focusing on the width and the length of phialides and conidia showed no difference among the clades within B. bassiana while conidial length was significantly different among clades within B. asiatica. The secondary metabolite profiles obtained by liquid chromatography-mass spectrometry (LC-MS) allowed a distinction between B. asiatica and B. bassiana, but not between the clades therein. Based on these genomic, morphological, chemical data, we proposed a clade of B. asiatica as a new species, named B. thailandica, and two clades of B. bassiana to respectively represent B. namnaoensis and B. neobassiana spp. nov. Such closely related but divergent species with different host ranges have potential to elucidate the evolution of host specificity, with potential biocontrol application. Citation: Kobmoo N, Arnamnart N, Pootakham W, et al. 2021. The integrative taxonomy of Beauveria asiatica and B. bassiana species complexes with whole-genome sequencing, morphometric and chemical analyses. Persoonia 47: 136-150. https://doi.org/10.3767/persoonia.2021.47.04.

Distribution and Diversity of Beauveria in Boreal Forests of Northern European Russia

The distribution and genetic diversity of 91 of Beauveria isolates collected during a long-term survey in boreal forests of northern European Russia was studied. Based on morphological and sequence analysis of TEF and Bloc loci, three Beauveria spp. were identified: B. pseudobassiana, B. bassiana, and B. caledonica, with abundance of 81, 11, and 8%, respectively. Through multilocus sequencing, four haplotypes of B. bassiana and two haplotypes of B. caledonica were detected. Twelve haplotypes of B. pseudobassiana with non-random distribution were identified. Two haplotypes of B. pseudobassiana were the most abundant and widespread occurring across the whole study area, whereas others tended to be more specific to either the north or south of the study area, indicating the presence of different subpopulations. For further analysis of these putative subpopulations, southern and northern areas were separated along the boundary of the Köppen-Geiger climate zones (dfb and dfc), and the genetic structure was examined by analysis of molecular variance and spatial autocorrelation. Molecular evidence of intraspecific recombination of B. pseudobassiana and B. bassiana across northern European Russia area was indicated.

Empirical Support for the Pattern of Competitive Exclusion between Insect Parasitic Fungi

Fungal entomopathogens are largely facultative parasites and play an important role in controlling the density of insect populations in nature. A few species of these fungi have been used for biocontrol of insect pests. The pattern of the entomopathogen competition for insect individuals is still elusive. Here, we report the empirical competition for hosts or niches between the inter- and intra-species of the entomopathogens Metarhizium robertsii and Beauveria bassiana. It was found that the synergistic effect of coinfection on virulence increase was not evident, and the insects were largely killed and mycosed by M. robertsii independent of its initial co-inoculation dosage and infection order. For example, >90% dead insects were mycosed by M. robertsii even after immersion in a spore suspension with a mixture ratio of 9:1 for B. bassiana versus M. robertsii. The results thus support the pattern of competitive exclusion between insect pathogenic fungi that occurred from outside to inside the insect hosts. Even being inferior to compete for insects, B. bassiana could outcompete M. robertsii during co-culturing in liquid medium. It was also found that the one-sided mycosis of insects occurred during coinfection with different genotypic strains of either fungi. However, parasexual recombination was evident to take place between the compatible strains after coinfection. The data of this study can help explain the phenomena of the exclusive mycosis of insect individuals, but co-occurrence of entomopathogens in the fields, and suggest that the synergistic effect is questionable regarding the mixed use of fungal parasites for insect pest control.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

Influence of genetic diversity of seventeen Beauveria bassiana isolates from different hosts on virulence by comparative genomics

Background

Beauveria bassiana (B. bassiana) is a famous entomopathogenic fungus that could parasitize on hundreds of insect species, which are being used as an environmentally friendly mycoinsecticide. Nevertheless, the possible effect of genetic diversity of these B. bassiana isolates from different hosts on virulence has not been explored before. In order to explore that issue, we compared the genome sequences among seventeen B. bassiana isolates from 17 different insects using whole genome re-sequencing, with B. bassiana strain ARSEF 2860 as the reference genome.

Results

There were a total of 10,098 missense mutated genes, 720 positively selected genes were identified in 17 strains of B. bassiana. Among these, two genes with high frequency mutations encode the toxin-producing non-ribosomal peptide synthase (NRPS) protein. Seven genes undergoing positive selection were enriched in the two-component signaling pathway that is known to regulate the fungal toxicity. In addition, the domain changes of three positively selected genes are also directly related to the virulence plasticity. Besides, the functional categorization of mutated genes showed that most of them involved in the biological functions of toxic proteins involved in.

Conclusions

Based on our data, our results indicate that several mutated genes and positively selected genes may underpin virulence of B. bassiana towards hosts during infection process, which provide an insight into the potential effects of natural variation on the virulence of B. bassiana, which will be useful in screening out potential virulence factors in B. bassiana.

Genome Sequence of Beauveria bassiana Strain ATCC 74040, a Widely Employed Insect Pathogen

The broad-spectrum insecticidal activity of Beauveria bassiana strain ATCC 74040 is well documented. The whole-genome sequence of this strain is reported here, revealing a plethora of genes related to its insecticidal potential and providing new insights on the mechanism of action.

The broad-spectrum insecticidal activity of Beauveria bassiana strain ATCC 74040 is well documented. The whole-genome sequence of this strain is reported here, revealing a plethora of genes related to its insecticidal potential and providing new insights on the mechanism of action.

Population genomics and evolution of a fungal pathogen after releasing exotic strains to control insect pests for 20 years

Entomopathogenic fungi are one of the key regulators of insect populations in nature. Some species such as Beauveria bassiana with a wide host range have been developed as promising alternatives to chemical insecticides for the biocontrol of insect pests. However, the long-term persistence of the released strains, the effect on non-target hosts and local fungal populations remains elusive, but they are considerable concerns with respect to environmental safety. Here we report the temporal features of the Beauveria population genomics and evolution over 20 years after releasing exotic strains to control pine caterpillar pests. We found that the isolates within the biocontrol site were mostly of clonal origins. The released strains could persist in the environment for a long time but with low recovery rates. Similar to the reoccurrence of host jumping by local isolates, the infection of non-target insects by the released strains was evident to endemically occur in association with host seasonality. No obvious dilution effect on local population structure was evident by the releases. However, the population was largely replaced by genetically divergent isolates once per decade but evolved with a pattern of balancing selection and towards expansion through adaptation, non-random outcrossing and isolate migration. This study not only unveils the real-time features of entomopathogenic fungal population genomics and evolution but also provides added values to alleviate the concerns of environmental safety regarding the biocontrol application of mycoinsecticides.

Two genes in a pathogenicity gene cluster encoding secreted proteins are required for appressorial penetration and infection of the maize anthracnose fungus Colletotrichum graminicola

To avoid pathogen-associated molecular pattern recognition, the hemibiotrophic maize pathogen Colletotrichum graminicola secretes proteins mediating the establishment of biotrophy. Targeted deletion of 26 individual candidate genes and seven gene clusters comprising 32 genes of C. graminicola identified a pathogenicity cluster (CLU5) of five co-linear genes, all of which, with the exception of CLU5b, encode secreted proteins. Targeted deletion of all genes of CLU5 revealed that CLU5a and CLU5d are required for full appressorial penetration competence, with virulence deficiencies independent of the host genotype and organ inoculated. Cytorrhysis experiments and microscopy showed that Δclu5a mutants form pressurized appressoria, but they are hampered in forming penetration pores and fail to differentiate a penetration peg. Whereas Δclu5d mutants elicited WT-like papillae, albeit at increased frequencies, papillae induced by Δclu5a mutants were much smaller than those elicited by the WT. Synteny of CLU5 is not only conserved in Colletotrichum spp. but also in additional species of Sordariomycetes including insect pathogens and saprophytes suggesting importance of CLU5 for fungal biology. Since CLU5a and CLU5d also occur in non-pathogenic fungi and since they are expressed prior to plant invasion and even in vegetative hyphae, the encoded proteins probably do not act primarily as effectors.

Hiding in plain sight: New virus genomes discovered via a systematic analysis of fungal public transcriptomes

The distribution and diversity of RNA viruses in fungi is incompletely understood due to the often cryptic nature of mycoviral infections and the focused study of primarily pathogenic and/or economically important fungi. As most viruses that are known to infect fungi possess either single-stranded or double-stranded RNA genomes, transcriptomic data provides the opportunity to query for viruses in diverse fungal samples without any a priori knowledge of virus infection. Here we describe a systematic survey of all transcriptomic datasets from fungi belonging to the subphylum Pezizomycotina. Using a simple but effective computational pipeline that uses reads discarded during normal RNA-seq analyses, followed by identification of a viral RNA-dependent RNA polymerase (RdRP) motif in de novo assembled contigs, 59 viruses from 44 different fungi were identified. Among the viruses identified, 88% were determined to be new species and 68% are, to our knowledge, the first virus described from the fungal species. Comprehensive analyses of both nucleotide and inferred protein sequences characterize the phylogenetic relationships between these viruses and the known set of mycoviral sequences and support the classification of up to four new families and two new genera. Thus the results provide a deeper understanding of the scope of mycoviral diversity while also increasing the distribution of fungal hosts. Further, this study demonstrates the suitability of analyzing RNA-seq data to facilitate rapid discovery of new viruses.

Pangenome analyses of the wheat pathogen Zymoseptoria tritici reveal the structural basis of a highly plastic eukaryotic genome

Background

Structural variation contributes substantially to polymorphism within species. Chromosomal rearrangements that impact genes can lead to functional variation among individuals and influence the expression of phenotypic traits. Genomes of fungal pathogens show substantial chromosomal polymorphism that can drive virulence evolution on host plants. Assessing the adaptive significance of structural variation is challenging, because most studies rely on inferences based on a single reference genome sequence.

Results

We constructed and analyzed the pangenome of Zymoseptoria tritici, a major pathogen of wheat that evolved host specialization by chromosomal rearrangements and gene deletions. We used single-molecule real-time sequencing and high-density genetic maps to assemble multiple genomes. We annotated the gene space based on transcriptomics data that covered the infection life cycle of each strain. Based on a total of five telomere-to-telomere genomes, we constructed a pangenome for the species and identified a core set of 9149 genes. However, an additional 6600 genes were exclusive to a subset of the isolates. The substantial accessory genome encoded on average fewer expressed genes but a larger fraction of the candidate effector genes that may interact with the host during infection. We expanded our analyses of the pangenome to a worldwide collection of 123 isolates of the same species. We confirmed that accessory genes were indeed more likely to show deletion polymorphisms and loss-of-function mutations compared to core genes.

Conclusions

The pangenome construction of a highly polymorphic eukaryotic pathogen showed that a single reference genome significantly underestimates the gene space of a species. The substantial accessory genome provides a cradle for adaptive evolution.

Electronic supplementary material

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

Culture degeneration in conidia of Beauveria bassiana and virulence determinants by proteomics

The quality of Beauveria bassiana conidia directly affects the virulence against insects. In this study, continuous subculturing of B. bassiana on both rice grains and potato dextrose agar (PDA) resulted in 55 and 49 % conidial yield reduction after 12 passages and 68 and 60 % virulence reduction after 20 and 12 passages at four d post-inoculation, respectively. The passage through Tenebrio molitor and Spodoptera exigua restored the virulence of rice and PDA subcultures, respectively. To explore the molecular mechanisms underlying the conidial quality and the decline of virulence after multiple subculturing, we investigated the conidial proteomic changes. Successive subculturing markedly increased the protein levels in oxidative stress response, autophagy, amino acid homeostasis, and apoptosis, but decreased the protein levels in DNA repair, ribosome biogenesis, energy metabolism, and virulence. The nitro blue tetrazolium assay verified that the late subculture's colony and conidia had a higher oxidative stress level than the early subculture. A 2A-type protein phosphatase and a Pleckstrin homology domain protein Slm1, effector proteins of the target of rapamycin (TOR) complex 1 and 2, respectively, were dramatically increased in the late subculture. These results suggest that TOR signalling might be associated with ageing in B. bassiana late subculture, in turn affecting its physiological characteristics and virulence.

Biodegradation of gentamicin by bacterial consortia AMQD4 in synthetic medium and raw gentamicin sewage

Gentamicin, a broad spectrum antibiotic of the aminoglycoside class, is widely used for disease prevention of human beings as well as animals. Nowadays the environmental issue caused by the disposal of wastes containing gentamicin attracts increasing attention. In this study, a gentamicin degrading bacterial consortia named AMQD4, including Providencia vermicola, Brevundimonas diminuta, Alcaligenes sp. and Acinetobacter, was isolated from biosolids produced during gentamicin production for the removal of gentamicin in the environment. The component and structure of gentamicin have a great influence on its degradation and gentamicin C1a and gentamicin C2a were more prone to being degraded. AMQD4 could maintain relatively high gentamicin removal efficiency under a wide range of pH, especially in an alkaline condition. In addition, AMQD4 could remove 56.8% and 47.7% of gentamicin in unsterilized and sterilized sewage in a lab-scale experiment, respectively. And among the isolates in AMQD4, Brevundimonas diminuta BZC3 performed the highest gentamicin degradation about 50%. It was speculated that aac3iia was the gentamicin degradation gene and the main degradation product was 3'-acetylgentamicin. Our results suggest that AMQD4 and Brevundimonas diminuta BZC3 could be important candidates to the list of superior microbes for bioremediation of antibiotic pollution.

Omics data reveal the unusual asexual-fruiting nature and secondary metabolic potentials of the medicinal fungus Cordyceps cicadae

Background

Ascomycete Cordyceps species have been using as valued traditional Chinese medicines. Particularly, the fruiting bodies of Cordyceps cicadae (syn. Isaria cicadae) have long been utilized for the treatment of chronic kidney disease. However, the genetics and bioactive chemicals in this fungus have been largely unexplored.

Results

In this study, we performed comprehensive omics analyses of C. cicadae, and found that, in contrast to other Cordyceps fungi, C. cicadae produces asexual fruiting bodies with the production of conidial spores instead of the meiotic ascospores. Genome sequencing and comparative genomic analysis indicate that the protein families encoded by C. cicadae are typical of entomopathogenic fungi, including the expansion of proteases and chitinases for targeting insect hosts. Interestingly, we found that the MAT1-2 mating-type locus of the sequenced strain contains an abnormally truncated MAT1-1-1 gene. Gene deletions revealed that asexual fruiting of C. cicadae is independent of the MAT locus control. RNA-seq transcriptome data also indicate that, compared to growth in a liquid culture, the putative genes involved in mating and meiosis processes were not up-regulated during fungal fruiting, further supporting asexual reproduction in this fungus. The genome of C. cicadae encodes an array of conservative and divergent gene clusters for secondary metabolisms. Based on our analysis, the production of known carcinogenic metabolites by this fungus could be potentially precluded. However, the confirmed production of oosporein raises health concerns about the frequent consumption of fungal fruiting bodies.

Conclusions

The results of this study expand our knowledge of fungal genetics that asexual fruiting can occur independent of the MAT locus control. The obtained genomic and metabolomic data will benefit future investigations of this fungus for medicinal uses.

Electronic supplementary material

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

Managing the plant microbiome for biocontrol fungi: examples from Hypocreales

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

Experimental evolution to increase the efficacy of the entomopathogenic fungus Beauveria bassiana against malaria mosquitoes: Effects on mycelial growth and virulence

Entomopathogenic fungi such as Beauveria bassiana are currently considered as a potential control agent for malaria mosquitoes. The success of such strategies depends among others on the efficacy of the fungus to kill its hosts. As B. bassiana can use various resources for growth and reproduction, increasing the dependency on mosquitoes as a nutritional source may be instrumental for reaching this goal. Passage of entomopathogenic fungi through an insect host has been shown to increase its virulence. We evaluated the virulence, fungal outgrowth, mycelial growth rate, and sporulation rate of two B. bassiana isolates (Bb1520 and Bb8028) that underwent 10 consecutive selection cycles through malaria mosquitoes (Anopheles coluzzii) using an experimental evolution approach. This cycling resulted in an altered capacity of evolved B. Bassiana lineages to grow on different substrates while maintaining the ability to kill insects. Notably, however, there were no significant changes in virulence or speed of outgrowth when comparing the evolved lineages against their unevolved ancestors. These results suggest that fungal growth and sporulation evolved through successive and exclusive use of an insect host as a nutritional resource. We discuss the results in light of biocontrol and provide suggestions to increase fungal virulence.