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

Entomopathogenic fungi in crops protection with an emphasis on bioactive metabolites and biological activities

Plant pathogens with their abundance are harmful and cause huge damage to different agricultural crops and economy of a country as well as lead towards the shortage of food for humans. For their management, the utilization of entomopathogenic fungi is an eco-friendly technique, sustainable to the environment, safe for humans and has promising effect over chemical-based pesticides. This process requires a biochemical mechanism, including the production of enzymes, toxins, and other metabolites that facilitate host infection and invasion. Essential enzymes such as chitinase, proteinase, and lipase play a direct role in breaking down the host cuticle, the primary barrier to EPF (Entomopathogenic Fungi) infection. Additionally, secondary metabolites such as destruxins in Metarhizium, beauvericin in Beauveria, hirsutellides in Hirsutella, isarolides in Isaria, cordyols in Cordyceps, and vertihemipterins in Verticillium, among others, act both directly and indirectly to disable the defense mechanisms of insect hosts, thereby accelerating the EPF infection process. The chemical composition of these secondary metabolites varies, ranging from simple non-peptide pigments such as oosporine to highly complex piperazine derivatives such as vertihemiptellides. The biocontrol efficacy of EPF is extensively studied, with numerous fungal strains commercially available on a large scale for managing arthropod pests. This review emphasizes the role of proteins and enzymes against crop pathogens, detailing their mode of action, and describing the metabolites from entomopathogenic fungi and their biological activities. In doing so, these findings contribute to establishing a symbiotic equilibrium between agricultural productivity and environmental conservation.

Clonal genomic population structure of Beauveria brongniartii and Beauveria pseudobassiana: Pathogens of the common European cockchafer (Melolontha melolontha L.)

Beauveria brongniartii is a fungal pathogen that infects the beetle Melolontha melolontha, a significant agricultural pest in Europe. While research has primarily focused on the use of B. brongniartii for controlling M. melolontha, the genomic structure of the B. brongniartii population remains unknown. This includes whether its structure is influenced by its interaction with M. melolontha, the timing of beetle-swarming flights, geographical factors, or reproductive mode. To address this, we analysed genome-wide SNPs to infer the population genomics of Beauveria spp., which were isolated from infected M. melolontha adults in an Alpine region. Surprisingly, only one-third of the isolates were identified as B. brongniartii, while two-thirds were distributed among cryptic taxa within B. pseudobassiana, a fungal species not previously recognized as a pathogen of M. melolontha. Given the prevalence of B. pseudobassiana, we conducted analyses on both species. We found no spatial or temporal genomic patterns within either species and no correlation with the population structure of M. melolontha, suggesting that the dispersal of the fungi is independent of the beetle. Both species exhibited clonal population structures, with B. brongniartii fixed for one mating type and B. pseudobassiana displaying both mating types. This implies that factors other than mating compatibility limit sexual reproduction. We conclude that the population genomic structure of Beauveria spp. is primarily influenced by predominant asexual reproduction and dispersal.

Fungal infection of insects: molecular insights and prospects

Entomopathogenic fungi (EPF) distribute in different fungal phyla with variable host ranges and play essential role in regulating insect populations by infecting hosts via cuticle penetration. The representative ascomycete EPF of Metarhizium and Beauveria species have been widely used in mechanistic investigations of fungus-insect interactions and as ecofriendly mycoinsecticides. Here, we review the function of diverse genes, pathways, and secondary metabolites associated with EPF stepwise infections. In particular, emerging evidence has shown that EPF have to outcompete insect ectomicrobiotas prior to penetrating cuticles, and subvert or evade host antifungal immunity by using effector-like proteins and chemicals like plant pathogens. Future prospects are discussed for a better understanding of fungal pathobiology, which will provide novel insights into microbe-animal interactions.

Fights on the surface prior to fungal invasion of insects

Entomopathogenic fungi (EPF) infect insects by landing on and penetrating cuticles. Emerging evidence has shown that, prior to the invasion of insects, fungal cells have to battle and overcome diverse challenges, including the host behavioral defenses, colonization resistance mediated by ectomicrobiotas, host recognition, and generation of enough penetration pressure. The ascomycete EPF such as Metarhizium and Beauveria can thus produce adhesive proteins and/or the exopolysaccharide mucilage to tightly glue fungal cells on cuticles. Producing antimicrobial peptides and chemical compounds can enable EPF to outcompete cuticular defensive microbes. The use of divergent membrane receptors, accumulation, and quick degradation of lipid droplets in conidial cells can help EPF recognize proper hosts and build up cellular turgor to breach cuticles for systematic invasion. Further investigations are still required to unveil the multifaceted and intricate relationships between EPF and insect hosts.

Fungal co-expression network analyses identify pathogen gene modules associated with host insect invasion

The broad host range fungal insect pathogen, Beauveria bassiana, has been commercialized as an alternative to chemical insecticides for pest control worldwide. B. bassiana represents a unique model system with which to examine host-pathogen interactions, and a wide range of genes and processes have been studied. However, significant aspects of virulence, particularly on the genomic scale, remain poorly studied. Here, we have combined available transcriptomes with three newly generated data sets for a combined total analysis of 76 deep-sequenced samples covering growth, development, stress responses, and infection during the life cycle of B. bassiana. Co-expression network analyses resulted in the identification of gene modules enriched during two critical stages of the infection process, namely (i) cuticle penetration and (ii) in vivo hyphal body (dimorphic transition) growth capable of avoiding innate and humoral immune defenses. These analyses identify unique signatures of metabolism, signaling, secondary metabolite production, host defense suppression, membrane reorganization, effector production, and secretion for each stage, including genetic regulators and epigenetic patterns. These data provide a comprehensive framework for understanding and probing fungal adaptations to its pathogenic life cycle and expand the candidate repertoire for continued dissection of the host-pathogen interaction. IMPORTANCE Insect fungal pathogens have evolved unique strategies for overcoming host structural and immunological defenses that span from the sclerotized cuticle to innate and humoral cellular responses. Two critical stages of the infection process involve (i) cuticle penetration and (ii) immune evasion within the insect hemocoel. A set of 76 global transcriptomic data for B. bassiana that include the cuticle penetration and hemocoel growth stages were analyzed for patterns (gene modules) of expression, yielding unique insights into these different life stages. These analyses integrate gene networks involved in fungal development, stress response and pathogenesis to further the systematic understanding of the global processes integral to the unique adaptation employed by fungal pathogens of insects.

Mycobiota of Potato-Cereal Soft Wraps and the Production Facility

The aim of this study was to investigate the mycobiota of potato-cereal soft wraps and in the processing area. Potato-cereal soft wraps are cooked, cold-stored and mashed potatoes that are kneaded together, usually with wheat flour, to form dough. In order to identify the main spoilage mould of this product, 150 visible mould colonies from mouldy wraps were identified. Five different mould species were isolated; Aspergillus niger, Penicillium brevicompactum, Penicillium commune, Penicllium corylophilum and Pencillium discolor. The dominating spoilage mould was Penicillium commune with 83.9% of the colonies. In order to study the mycobiota of the production area, 271 samples of air and surfaces were collected. In total, 647 mould colonies were isolated from air and surface samples. The mycobiota of air consisted of 27 different species within 9 different genera, and the mycobiota of surfaces consisted of 14 species within 4 different genera. Penicllium species were the dominating genera both in air and on surfaces, and Penicillium commune was the dominating species in the processing environment as well. Penicillium commune was found in the bakery and also in other production rooms. Spores from the flour and from soil on potatoes can disperse in the air as aerosols and may contaminate the wraps after baking when the product is cooled before packaging.

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.

The Entomopathogenic Fungus Beauveria bassiana Employs Autophagy as a Persistence and Recovery Mechanism during Conidial Dormancy

Many filamentous fungi develop a conidiation process as an essential mechanism for their dispersal and survival in natural ecosystems. However, the mechanisms underlying conidial persistence in environments are still not fully understood. Here, we report that autophagy is crucial for conidial lifespans (i.e., viability) and vitality (e.g., stress responses and virulence) in the filamentous mycopathogen Beauveria bassiana. Specifically, Atg11-mediated selective autophagy played an important, but not dominant, role in the total autophagic flux. Furthermore, the aspartyl aminopeptidase Ape4 was found to be involved in conidial vitality during dormancy. Notably, the vacuolar translocation of Ape4 was dependent on its physical interaction with autophagy-related protein 8 (Atg8) and associated with the autophagic role of Atg8, as determined through a truncation assay of a critical carboxyl-tripeptide. These observations revealed that autophagy acted as a subcellular mechanism for conidial recovery during dormancy in environments. In addition, a novel Atg8-dependent targeting route for vacuolar hydrolase was identified, which is essential for conidial exit from a long-term dormancy. These new insights improved our understanding of the roles of autophagy in the physiological ecology of filamentous fungi as well as the molecular mechanisms involved in selective autophagy. IMPORTANCE Conidial environmental persistence is essential for fungal dispersal in ecosystems while also serving as a determinant for the biocontrol efficacy of entomopathogenic fungi during integrated pest management. This study identified autophagy as a mechanism to safeguard conidial lifespans and vitality postmaturation. In this mechanism, the aspartyl aminopeptidase Ape4 translocates into vacuoles via its physical interaction with autophagy-related protein 8 (Atg8) and is involved in conidial vitality during survival. The study revealed that autophagy acted as a subcellular mechanism for maintaining conidial persistence during dormancy, while also documenting an Atg8-dependent targeting route for vacuolar hydrolase during conidial recovery from dormancy. Thus, these observations provided new insight into the roles of autophagy in the physiological ecology of filamentous fungi and documented novel molecular mechanisms involved in selective autophagy.

Suppression of the insect cuticular microbiomes by a fungal defensin to facilitate parasite infection

Insects can assemble defensive microbiomes on their body surfaces to defend against fungal parasitic infections. The strategies employed by fungal pathogens to combat host cuticular microbiotas remains unclear. Here, we report the identification and functional characterization of the defensin-like antimicrobial gene BbAMP1 encoded by the entomopathogenic fungus Beauveria bassiana. The mature peptide of BbAMP1 can coat fungal spores and can be secreted by the fungus to target and damage Gram-positive bacterial cells. Significant differences in insect survival were observed between the wild-type and BbAMP1 mutant strains during topical infection but not during injection assays that bypassed insect cuticles. Thus, BbAMP1 deletion considerably reduced fungal virulence while gene overexpression accelerated the fungal colonization of insects compared with the wild-type strain in natural infections. Topical infection of axenic Drosophila adults evidenced no difference in fly survivals between strains. However, the gnotobiotic infections with the addition of Gram-positive but not Gram-negative bacterial cells in fungal spore suspensions substantially increased the survival of the flies treated with ∆BbAMP1 compared to those infected by the wild-type and gene-overexpression strains. Bacterial colony counts and microbiome analysis confirmed that BbAMP1 could assist the fungus to manipulate insect surface bacterial loads. This study reveals that fungal defensin can suppress the host surface defensive microbiomes, which underscores the importance to extend the research scope of fungus-host interactions.

Promotion of Arabidopsis immune responses by a rhizosphere fungus via supply of pipecolic acid to plants and selective augment of phytoalexins

The ascomycete insect pathogenic fungi such as Metarhizium species have been demonstrated with the abilities to form the rhizosphere or endophytic relationships with different plants for nutrient exchanges. In this study, after the evident infeasibility of bacterial disease development in the boxed sterile soils, we established a hydroponic system for the gnotobiotic growth of Arabidopsis thaliana with the wild-type and transgenic strain of Metarhizium robertsii. The transgenic fungus could produce a high amount of pipecolic acid (PIP), a pivotal plant-immune-stimulating metabolite. Fungal inoculation experiments showed that M. robertsii could form a non-selective rhizosphere relationship with Arabidopsis. Similar to the PIP uptake by plants after exogenous application, PIP level increased in Col-0 and could be detected in the PIP-non-producing Arabidopsis mutant (ald1) after fungal inoculations, indicating that plants can absorb the PIP produced by fungi. The transgenic fungal strain had a better efficacy than the wild type to defend plants against the bacterial pathogen and aphid attacks. Contrary to ald1, fmo1 plants could not be boosted to resist bacterial infection after treatments. After fungal inoculations, the phytoalexins camalexin and aliphatic glucosinolate were selectively increased in Arabidopsis via both PIP-dependent and -independent ways. This study unveils the potential mechanism of the fungus-mediated beneficial promotion of plant immunity against biological stresses. The data also highlight the added values of M. robertsii to plants beyond the direct suppression of insect pest populations.

Towards Understanding the Function of Aegerolysins

Aegerolysins are remarkable proteins. They are distributed over the tree of life, being relatively widespread in bacteria and fungi, but also present in some insects, plants, protozoa, and viruses. Despite their abundance in cells of certain developmental stages and their presence in secretomes, only a few aegerolysins have been studied in detail. Their function, in particular, is intriguing. Here, we summarize previously published findings on the distribution, molecular interactions, and function of these versatile aegerolysins. They have very diverse protein sequences but a common fold. The machine learning approach of the AlphaFold2 algorithm, which incorporates physical and biological knowledge of protein structures and multisequence alignments, provides us new insights into the aegerolysins and their pore-forming partners, complemented by additional genomic support. We hypothesize that aegerolysins are involved in the mechanisms of competitive exclusion in the niche.

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.

Systematic contributions of CFEM domain-containing proteins to iron acquisition are essential for interspecies interaction of the filamentous pathogenic fungus Beauveria bassiana

Common in fungal extracellular membrane (CFEM) domain is unique in fungal proteins and some of which contribute to iron acquisition in yeast. However, their roles in iron acquisition remain largely unknown in filamentous fungi. In this study, 12 CFEM-containing proteins were bioinformatically identified in the filamentous entomopathogenic fungus Beauveria bassiana, and the roles of 11 genes were genetically characterized. Transmembrane helices were critical for their association with intracellular membranes, and their number varied among proteins. Eleven CFEM genes significantly contribute to vegetative growth under iron starvation and virulence. Notably, the virulence of most disruptants could be significantly weakened by a decrease in iron availability, in which the virulence of ΔBbcfem7 and 8 strains was partially recovered by exogenous hemin. ΔBbcfem7 and 8 mutants displayed defective competitiveness against the sister entomopathogenic fungus Beauveria brongniartii. All 11 disruptants displayed impaired growth in the antagonistic assay with the saprotrophic fungus Aspergillus niger, which could be repressed by exogenous ferric ions. These findings not only reveal the systematic contributions of CFEM proteins to acquire two forms of iron (i.e. heme and ferric ion) in the entire lifecycle of entomopathogenic fungi but also help to better understand the mechanisms of fungus-host and inter-fungus interactions.

Temporal changes in pathogen diversity in a perennial plant-pathogen-hyperparasite system

Hyperparasites can affect the evolution of pathosystems by influencing the stability of both pathogen and host populations. However, how pathogens of perennial hosts evolve in the presence of a hyperparasite has rarely been studied. Here, we investigated temporal changes in genetic diversity of the invasive chestnut blight pathogen Cryphonectria parasitica in the presence of its parasitic mycovirus Cryphonectria hypovirus 1 (CHV1). The virus reduces fungal virulence and represents an effective natural biocontrol agent against chestnut blight in Europe. We analysed genome‐wide diversity and CHV1 prevalence in C. parasitica populations in southern Switzerland that were sampled twice at an interval of about 30 years. Overall, we found that both pathogen population structure and CHV1 prevalence were retained over time. The results suggest that recent bottlenecks have influenced the structure of C. parasitica populations in southern Switzerland. Strong balancing selection signals were found at a single vegetative incompatibility (vic) locus, consistent with negative frequency‐dependent selection imposed by the vegetative incompatibility system. High levels of mating among related individuals (i.e., inbreeding) and genetic drift are probably at the origin of imbalanced allele ratios at vic loci and subsequently low vc type diversity. Virus infection rates were stable at ~30% over the study period and we found no significant impact of the virus on fungal population diversity. Consequently, the efficacy of CHV1‐mediated biocontrol was probably retained.

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.

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.

iTRAQ-based quantitative proteomic analysis of silkworm infected with Beauveria bassiana

Beauveria bassiana is a harmful pathogen to the economically important insect silkworm, always causes serious disease to the silkworm, which results in great losses to the sericulture industry. In order to explore the silkworm (Bombyx mori) response to B. bassiana infection, differential proteomes of the silkworm responsive to B. bassiana infection were identified with isobaric tags for relative and absolute quantitation (iTRAQ) at the different stage of the 3rd instar silkworm larvae. Among the 5040 proteins identified with confidence level of ≥95 %, total 937 proteins were differentially expressed, of which 488 proteins were up-regulated and 449 proteins were down-regulated. 23, 15, 250, 649 differentially expressed proteins (DEPs) were reliably quantified by iTRAQ analysis in the B. bassiana infected larvae at 18, 24, 36, 48 h post infection (hpi) respectively. Based on GO annotations, 6, 4, 128, 316 DEPs were involved in biological processes, 12, 5, 143, 376 DEPs were involved in molecular functions, and 6, 3, 108, 256 DEPs were involved in cell components at 18, 24, 36, 48 hpi respectively. KEGG pathway analysis displayed that 18, 12, 210, 548 DEPs separately participated in 63, 35, 201, 264 signal transduction pathways at different time of infection, and moreover a higher proportion of DEPs involved in metabolic pathways. The cluster analysis on the DEPs of different infection stages distinguished a co-regulated DEP, lysozyme precursor, which was up-regulated at both the mRNA level and the protein level, indicating that the lysozyme protein kept playing an important role in defending the silkworm against B. bassiana infection. This was the first report using an iTRAQ approach to analyze proteomes of the whole silkworm against B. bassiana infection, which contributes to better understanding the defense mechanisms of silkworm to B. bassiana infection and provides important experimental data for the identification of key factors involved in the interaction between the pathogenic fungus and its host.

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.

Production of Diverse Beauveriolide Analogs in Closely Related Fungi: a Rare Case of Fungal Chemodiversity

Fungal chemodiversity is well known in part due to the production of diverse analogous compounds by a single biosynthetic gene cluster (BGC). Usually, similar or the same metabolites are produced by closely related fungal species under a given condition, the foundation of fungal chemotaxonomy. Here, we report a rare case of the production of the cyclodepsipeptide beauveriolides (BVDs) in three insect-pathogenic fungi. We found that the more closely related fungi Beauveria bassiana and Beauveria brongniartii produced structurally distinct analogs of BVDs, whereas the less-close relatives B. brongniartii and Cordyceps militaris biosynthesized structurally similar congeners under the same growth condition. It was verified that a conserved BGC containing four genes is responsible for BVD biosynthesis in three fungi, including a polyketide synthase (PKS) for the production of 3-hydroxy fatty acids (FAs) with chain length variations. In contrast to BVD production patterns, phylogenetic analysis of the BGC enzymes or enzyme domains largely resulted in the congruence relationship with fungal speciation. Feeding assays demonstrated that an FA with a chain length of eight carbon atoms was preferentially utilized, whereas an FA with a chain longer than 10 carbon atoms could not be used as a substrate for BVD biosynthesis. Insect survival assays suggested that the contribution of BVDs to fungal virulence might be associated with the susceptibility of insect species. The results of this study enrich the knowledge of fungal secondary metabolic diversity that can question the reliability of fungal chemotaxonomy.IMPORTANCE Fungal chemotaxonomy is an approach to classify fungi based on the fungal production profile of metabolites, especially the secondary metabolites. We found an atypical example that could question the reliability of fungal chemical classifications in this study, i.e., the more closely related entomopathogenic species Beauveria bassiana and Beauveria brongniartii produced structurally different congeners of the cyclodepsipeptide beauveriolides, whereas the rather divergent species B. brongniartii and Cordyceps militaris biosynthesized similar analogs under the same growth condition. The conserved biosynthetic gene cluster (BGC) containing four genes present in each species is responsible for beauveriolide production. In contrast to the compound formation profiles, the phylogenies of biosynthetic enzymes or enzymatic domains show associations with fungal speciation. Dependent on the insect species, production of beauveriolides may contribute to fungal virulence against the susceptible insect hosts. The findings in this study augment the diversity of fungal secondary metabolisms.