Divergent and Convergent Evolution of Fungal Pathogenicity

Ophiostomatoid fungi associated with pine bark beetles and infested pines in south-eastern Australia, including Graphilbum ipis-grandicollis sp. nov

The ophiostomatoid fungi are an assemblage of ascomycetes which are arguably best-known for their associations with bark and ambrosia beetles (Curculonidae) and blue stain (sap stain) of many economically important tree species. These fungi are considered a significant threat to coniferous forests, which has resulted in numerous studies characterising the diversity of bark beetles and their ophiostomatoid associates globally. The diversity of ophiostomatoid fungi present in Australian pine plantations, however, remains largely undetermined. The aims of this study were therefore to reconsider the diversity of ophiostomatoid fungi associated with Pinus in Australia, and to establish the baseline of expected taxa found within these plantation ecosystems. To achieve this, we reviewed Australian plant pathogen reference collections, and analysed samples collected during forest health surveillance programs from the major pine growing regions in south-eastern Australia. In total, 135 ophiostomatoid isolates (15 from reference collections and 120 collected during the current study) were assessed using morphological identification and ITS screening which putatively distinguished 15 taxonomic groups. Whole genome sequencing (WGS) of representative isolates from each taxon was performed to obtain high-quality sequence data for multi-locus phylogenetic analysis. Our results revealed a greater than expected diversity, expanding the status of ophiostomatoid fungi associated with Pinus in Australia to include 14 species from six genera in the Ophiostomatales and a single species residing in the Microascales. While most of these were already known to science, our study includes seven first records for Australia and the description of one new species, Graphilbum ipis-grandicollis sp. nov.. This study also provides an early example of whole genome sequencing (WGS) approaches replacing traditional PCR-based methods for taxonomic surveys. This not only allowed for robust multi-locus sequence extraction during taxonomic assessment, but also permitted the rapid establishment of a curated genomic database for ophiostomatoid fungi which will continue to aid in the development of improved diagnostic resources and capabilities for Australian biosecurity.

Regulation of conidiation, polarity growth, and pathogenicity by MrSte12 transcription factor in entomopathogenic fungus, Metarhizium rileyi

Metarhizium rileyi, a well-known filamentous biocontrol fungus, is the main pathogen of numerous field pests, especially noctuid pests. To explore the potential factors involved in the fungal pathogenicity, MrSte12, an important and conserved functional transcription factor in mitogen-activated protein kinase pathway was carried out by functional analysis. Homologous recombination was used to disrupt the MrSte12 gene in M. rileyi. The deletant fungal strain exhibited malformed hyphae and impaired conidiogenesis, and conidia could not be collected from △MrSte12 in vitro towards SMAY medium. Although conidia could be collected again supplemented with KCl within SMAY medium, the conidial germination, growth and stress tolerance were much weaker compared with that in WT. Additionally, △MrSte12 showed a dramatic reduction in virulence in intra-hemolymph injections and no pathogenicity in topical inoculations against noctuid pests, which is due to the failure of appressorium formation. Moreover, the content of chitin and β-1, 3-glucan in cell wall significantly reduced in mutant conidia. These results indicate that the MrSte12 gene markedly contributes to invasive growth and conidiation, as well as the major pathogenicity in M. rileyi.

Mitochondrial genomes of the human pathogens Coccidioides immitis and Coccidioides posadasii

Fungal mitochondrial genomes encode genes involved in crucial cellular processes, such as oxidative phosphorylation and mitochondrial translation, and the molecule has been used as a molecular marker for population genetics studies. Coccidioides immitis and C. posadasii are endemic fungal pathogens that cause coccidioidomycosis in arid regions across both American continents. To date, approximately 150 Coccidioides isolates have been sequenced to infer patterns of variation in nuclear genomes. However, less attention has been given to the mitochondrial genomes of Coccidioides. In this report, we describe the assembly and annotation of mitochondrial reference genomes for two representative strains of C. posadasii and C. immitis, as well as assess population variation among 77 selected genomes. The sizes of the circular-mapping molecules are 68.2 Kb in C. immitis and 75.1 Kb in C. posadasii. We identify 14 mitochondrial protein-coding genes common to most fungal mitochondria, which are largely syntenic across different populations and species of Coccidioides. Both Coccidioides species are characterized by a large number of group I and II introns, harboring twice the number of elements as compared to closely related Onygenales. The introns contain complete or truncated ORFs with high similarity to homing endonucleases of the LAGLIDADG and GIY-YIG families. Phylogenetic comparisons of mitochondrial and nuclear genomes show extensive phylogenetic discordance suggesting that the evolution of the two types of genetic material is not identical. This work represents the first assessment of mitochondrial genomes among isolates of both species of Coccidioides, and provides a foundation for future functional work.

Comparative roles of three adhesin genes (adh1-3) in insect-pathogenic lifecycle of Beauveria bassiana

Adherence of conidia to insect integument is crucial for initiation of fungal infection through cuticular penetration and was previously reported to rely upon the Metarhizium-type adhesin Mad1 rather than Mad2, another adhesin crucial for conidial adherence of Metarhizium anisopliae to plant root surface. Mad1 and Mad2 have since been considered to function in fungal insect pathogenesis and plant root colonization respectively. Here, three adhesins were characterized in Beauveria bassiana, including Adh1/Mad1, Adh2/Mad2, and Adh3 known as filamentous hemagglutinin/adhesin and virulence factor in animal-pathogenic bacteria. Among those, only Adh2 was found to play a substantial role in sustaining the fungal virulence and some phenotypes associated with biological control potential. Disruption of adh2 resulted in decreased conidial adherence to insect wing cuticle, attenuated virulence via normal cuticle infection or cuticle-bypassing infection (injection), reduced blastospore production in an insect hemolymph-mimicking broth, largely reduced conidiation capacity, impaired conidial quality indicative of lowered viability, hydrophobicity, and UV resistance, but no growth defects on rich and scant media under normal or stressful culture conditions. The main phenotypic changes correlated well with repressed expression of developmental activator genes required for aerial conidiation and submerged blastospore production and of key hydrophobin genes essential for hydrophobin synthesis and assembly into rodlet bundles of conidial coat crucial for conidial adherence. In contrast, either adh1 or adh3 disruption caused insignificant changes in all phenotypes examined. These findings offer novel insight into a significance of Adh2, but a dispensability of Adh1 or Adh3, for insect-pathogenic lifecycle of B. bassiana. KEY POINTS: • Three adhesins (Adh1-3) of Beauvera bassiana are functionally characterized. • Adh2 plays a role in sustaining virulence and lifecycle-related cellular events. • Either Adh1 or Adh3 is dispensable for insect-pathogenic lifecycle of B. bassiana.

Predicting Lifestyle from Positive Selection Data and Genome Properties in Oomycetes

As evidenced in parasitism, host and niche shifts are a source of genomic and phenotypic diversification. Exemplary is a reduction in the core metabolism as parasites adapt to a particular host, while the accessory genome often maintains a high degree of diversification. However, selective pressures acting on the genome of organisms that have undergone recent lifestyle or host changes have not been fully investigated. Here, we developed a comparative genomics approach to study underlying adaptive trends in oomycetes, a eukaryotic phylum with a wide and diverse range of economically important plant and animal parasitic lifestyles. Our analysis reveals converging evolution on biological processes for oomycetes that have similar lifestyles. Moreover, we find that certain functions, in particular carbohydrate metabolism, transport, and signaling, are important for host and environmental adaptation in oomycetes. Given the high correlation between lifestyle and genome properties in our oomycete dataset, together with the known convergent evolution of fungal and oomycete genomes, we developed a model that predicts plant pathogenic lifestyles with high accuracy based on functional annotations. These insights into how selective pressures correlate with lifestyle may be crucial to better understand host/lifestyle shifts and their impact on the genome.

Comparative genomics reveals a core gene toolbox for lifestyle transitions in Hypocreales fungi

Fungi have evolved diverse lifestyles and adopted pivotal new roles in both natural ecosystems and human environments. However, the molecular mechanisms underlying their adaptation to new lifestyles are obscure. Here, we hypothesize that genes shared across all species with the same lifestyle, but absent in genera with alternative lifestyles, are crucial to that lifestyle. By analysing dozens of species within four genera in a fungal order, with each genus following a different lifestyle, we find that genus-specific genes are typically few in number. Notably, not all genus-specific genes appear to derive from de novo birth, with most instead reflecting recurrent loss across the fungi. Importantly, however, a subset of these genus-specific genes are shared by fungi with the same lifestyle in quite different evolutionary orders, thus supporting the view that some genus-specific genes are necessary for specific lifestyles. These lifestyle-specific genes are enriched for key functional classes and often exhibit specialized expression patterns. Genus-specific selection also contributes to lifestyle transitions, and is especially associated with intensity of pathogenesis. Our study, therefore, suggests that fungal adaptation to new lifestyles often requires just a small number of core genes, with gene turnover and positive selection playing complementary roles.

The G-protein coupled receptor GPRK contributes to fungal development and full virulence in Metarhizium robertsii

G-protein-coupled receptor K (GPRK), which is a class VI fungal G-protein-coupled receptor (GPCR), plays a critical role in plant immunity against pathogens by mediating the endocytic pathway, influencing metabolism in response to environmental signals, and regulating asexual reproduction and pathogenic development. However, the function of these proteins in entomopathogenic fungi has rarely been investigated. Accordingly, we characterized MrGPRK, a GPCR in the entomopathogenic fungus Metarhizium robertsii containing a C-terminal seven-transmembrane and a conserved regulator of G protein signaling domain, and found that it localized to endosomes. Mutant phenotype assays showed that a ΔMrGprk strain displayed increased defects in radial growth (~28%) and decreased conidial production (~80%) compared with a wild-type strain. Decreased conidiation rates coincided well with the repression of conidiation-related regulatory genes, including three key conidial transcription factors: brlA, abaA, and wetA. MrGprk deficiency impaired full virulence (both topical and injectable inoculations). Further analysis demonstrated that deleting fungal MrGprk decreased the rates of appressorium formation and suppressed the transcription of several genes contributing to appressorial turgor pressure, cuticle penetration, and pH regulation. Additionally, the ΔMrGprk strain showed higher cyclic (cAMP) levels, suggesting that this GPCR is critical for cAMP signal transduction. In summary, MrGPRK was found to contribute to vegetative growth, conidial production, and full virulence of M. robertsii. These findings are conducive to a better understanding of the roles of GPCRs in the development and pathogenicity of entomopathogenic fungi.

Research Progress on Phytopathogenic Fungi and Their Role as Biocontrol Agents

Phytopathogenic fungi decrease crop yield and quality and cause huge losses in agricultural production. To prevent the occurrence of crop diseases and insect pests, farmers have to use many synthetic chemical pesticides. The extensive use of these pesticides has resulted in a series of environmental and ecological problems, such as the increase in resistant weed populations, soil compaction, and water pollution, which seriously affect the sustainable development of agriculture. This review discusses the main advances in research on plant-pathogenic fungi in terms of their pathogenic factors such as cell wall-degrading enzymes, toxins, growth regulators, effector proteins, and fungal viruses, as well as their application as biocontrol agents for plant pests, diseases, and weeds. Finally, further studies on plant-pathogenic fungal resources with better biocontrol effects can help find new beneficial microbial resources that can control diseases.

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.

Effects of Cordycepin in Cordyceps militaris during Its Infection to Silkworm Larvae

Cordyceps militaris produces cordycepin, a secondary metabolite that exhibits numerous bioactive properties. However, cordycepin pharmacology in vivo is not yet understood. In this study, the roles of cordycepin in C. militaris during its infection were investigated. After the injection of conidia, C. militaris NBRC100741 killed silkworm larvae more rapidly than NBRC103752. At 96 and 120 h, Cmcns genes (Cmcns1-4), which are part of the cordycepin biosynthesis gene cluster, were expressed in fat bodies and cuticles. Thus, cordycepin may be produced in the infection of silkworm larvae. Further, cordycepin enhanced pathogenicity toward silkworm larvae of Metarhizium anisopliae and Beauveria bassiana, that are also entomopathogenic fungi and do not produce cordycepin. In addition, by RNA-seq analysis, the increased expression of the gene encoding a lipoprotein 30K-8 (Bmlp20, KWMTBOMO11934) and decreased expression of genes encoding cuticular proteins (KWMTBOMO13140, KWMTBOMO13167) and a serine protease inhibitor (serpin29, KWMTBOMO08927) were observed when cordycepin was injected into silkworm larvae. This result suggests that cordycepin may aid the in vivo growth of C. militaris in silkworm larvae by the influence of the expression of some genes in silkworm larvae.

MripacC regulates blastosphere budding and influences virulence of the pathogenic fungus Metarhizium rileyi

Fungal dimorphism is the ability of certain fungi to switch between two different cellular forms, yeast and mycelial forms, in response to external environmental factors. The pacC/Pal signal transduction pathway responds to neutral and alkaline environments and is also involved in the fungal dimorphic transition. In this study, we investigated the function of the pacC homolog, MripacC, which regulates the dimorphic transition and modulates virulence of the insect pathogenic fungus Metarhizium rileyi. MripacC expression was upregulated under alkaline condition, with increased number of yeast-like cells compared to the number of hyphae cells. A MripacC deletion mutant (ΔMripacC) was obtained by homologous replacement and exhibited decreased blastospore budding, with direct development of conidia into hyphae without entering the yeast-like stage when cultured on alkaline medium. Observation of host hemolymph morphology and analysis of samples to detect the main immune factors revealed a decreased ability of ΔMripacC to evade the host immune system. The results of insect bioassays showed that ΔMripacC had decreased virulence with extended median lethality time. Together, the results suggested that MripacC not only regulated adaptation to acidic and alkaline environments, but also influenced virulence by budding blastospores. This elucidation of the function of MripacC adds to our understanding of blastospore budding and virulence of this fungal pathogen.

Subtilisin-like Pr1 proteases marking the evolution of pathogenicity in a wide-spectrum insect-pathogenic fungus

Subtilisin-like Pr1 proteases of insect-pathogenic fungi are a large family of extracellular cuticle-degrading enzymes that presumably determine a capability of hyphal invasion into insect hemocoel through normal cuticle infection, but remain poorly understood although often considered as virulence factors for genetic improvement of fungal potential against pests. Here, we report that not all of 11 Pr1 family members necessarily function in Beauveria bassiana, an ancient wide-spectrum pathogen evolved insect pathogenicity ~200 million years ago. These Pr1 proteases are phylogenetically similar to or distinct from 11 homologues (Pr1A–K) early named in Metarhizium anisopliae complex, a young entomopathogen lineage undergoing molecular evolution toward Pr1 diversification, and hence renamed Pr1A1/A2, Pr1B1–B3, Pr1 C, Pr1F1–F4,4 and Pr1 G, respectively. Multiple analyses of all single gene-deleted and rescued mutants led to the recognition of five conserved members (Pr1C, Pr1G, Pr1A2, Pr1B1, and Pr1B2) contributing significantly to the fungal pathogenicity to insect. The conserved Pr1 proteases were proven to function only in cuticle degradation, individually contribute 19–29% to virulence, but play no role in post-infection cellular events critical for fungal killing action. Six other Pr1 proteases were not functional at all in either cuticle degradation during host infection or virulence-related cellular events post-infection. Therefore, only the five conserved proteases are collectively required for, and hence mark evolution of, insect pathogenicity in B. bassiana. These findings provide the first referable base for insight into the evolution of Pr1 family members in different lineages of fungal insect pathogens.

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.

Design of a diagnostic system based on molecular markers derived from the ascomycetes pan-genome analysis: The case of Fusarium dieback disease

A key factor to take actions against phytosanitary problems is the accurate and rapid detection of the causal agent. Here, we develop a molecular diagnostics system based on comparative genomics to easily identify fusariosis and specific pathogenic species as the Fusarium kuroshium, the symbiont of the ambrosia beetle Euwallaceae kuroshio Gomez and Hulcr which is responsible for Fusarium dieback disease in San Diego CA, USA. We performed a pan-genome analysis using sixty-three ascomycetes fungi species including phytopathogens and fungi associated with the ambrosia beetles. Pan-genome analysis revealed that 2,631 orthologue genes are only shared by Fusarium spp., and on average 3,941 (SD ± 1,418.6) are species-specific genes. These genes were used for PCR primer design and tested on DNA isolated from i) different strains of ascomycete species, ii) artificially infected avocado stems and iii) plant tissue of field-collected samples presumably infected. Our results let us propose a useful set of primers to either identify any species from Fusarium genus or, in a specific manner, species such as F. kuroshium, F. oxysporum, and F. graminearum. The results suggest that the molecular strategy employed in this study can be expanded to design primers against different types of pathogens responsible for provoking critical plant diseases.

OSIP1 is a self-assembling DUF3129 protein required to protect fungal cells from toxins and stressors

Secreted proteins are key players in fungal physiology and cell protection against external stressing agents and antifungals. Oak stress-induced protein 1 (OSIP1) is a fungal-specific protein with unknown function. By using Podospora anserina and Phanerochaete chrysosporium as models, we combined both in vivo functional approaches and biophysical characterization of OSIP1 recombinant protein. The P. anserina OSIP1^Δ mutant showed an increased sensitivity to the antifungal caspofungin compared to the wild type. This correlated with the production of a weakened extracellular exopolysaccharide/protein matrix (ECM). Since the recombinant OSIP1 from P. chrysosporium self-assembled as fibers and was capable of gelation, it is likely that OSIP1 is linked to ECM formation that acts as a physical barrier preventing drug toxicity. Moreover, compared to the wild type, the OSIP1^Δ mutant was more sensitive to oak extractives including chaotropic phenols and benzenes. It exhibited a strongly modified secretome pattern and an increased production of proteins associated to the cell-wall integrity signalling pathway, when grown on oak sawdust. This demonstrates that OSIP1 has also an important role in fungal resistance to extractive-induced stress.

A M35 family metalloprotease is required for fungal virulence against insects by inactivating host prophenoloxidases and beyond

A diverse family of metalloproteases (MPs) is distributed in eukaryotes. However, the functions of MPs are still understudied. We report that seven MPs belonging to the M35 family are encoded in the genome of the insect pathogenic fungus Metarhizium robertsii. By gene deletions and insect bioassays, we found that one of the M35-family MPs, i.e. MrM35-4, is required for fungal virulence against insect hosts. MrM35-4 is a secretable enzyme and shows a proteolytic activity implicated in facilitating fungal penetration of insect cuticles. After gene rescue and overexpression, insect bioassays indicated that MrM35-4 contributes to inhibiting insect cuticular and hemocyte melanization activities. Enzymatic cleavage assays revealed that the recombinant prophenoloxidases PPO1 and PPO2 of Drosophila melanogaster could be clipped by MrM35-4 in a manner differing from a serine protease that can activate PPO activities. In addition, it was found that MrM35-4 is involved in suppressing antifungal gene expression in insects. Consistent with the evident apoptogenic effect of MrM35-4 on host cells, we found that the PPO mutant flies differentially succumbed to the infections of the wild-type and mutant strains of M. robertsii. Thus, MrM35-4 plays a multifaceted role beyond targeting PPOs during fungus-insect interactions, which represents a previously unsuspected strategy employed by Metarhizium to outmaneuver insect immune defenses.

P-type Na+/K+ ATPases essential and nonessential for cellular homeostasis and insect pathogenicity of Beauveria bassiana

ENA1 and ENA2 are P-type IID/ENA Na+/K+-ATPases required for cellular homeostasis in yeasts but remain poorly understood in filamentous fungal insect pathogens. Here, we characterized seven genes encoding five ENA1/2 homologues (ENA1a-c and ENA2a/b) and two P-type IIC/NK Na+/K+-ATPases (NK1/2) in Beauveria bassiana, an insect-pathogenic fungus serving as a main source of fungal insecticides worldwide. Most of these genes were highly responsive to alkaline pH and Na+/K+ cues at transcription level. Cellular Na+, K+ and H+ homeostasis was disturbed only in the absence of ena1a or ena2b. The disturbed homeostasis featured acceleration of vacuolar acidification, elevation of cytosolic Na+/K+ level at pH 5.0 to 9.0, and stabilization of extracellular H+ level to initial pH 7.5 during a 5-day period of submerged incubation. Despite little defect in hyphal growth and asexual development, the Δena1a and Δena2b mutants were less tolerant to metal cations (Na+, K+, Li+, Zn2+, Mn2+ and Fe3+), cell wall perturbation, oxidation, non-cation hyperosmolarity and UVB irradiation, severely compromised in insect pathogenicity via normal cuticle infection, and attenuated in virulence via hemocoel injection. The deletion mutants of five other ENA and NK genes showed little change in vacuolar pH and all examined phenotypes. Therefore, only ENA1a and ENA2b evidently involved in both transmembrane and vacuolar activities are essential for cellular cation homeostasis, insect pathogenicity and multiple stress tolerance in B. bassiana. These findings provide a novel insight into ENA1a- and ENA2b-dependent vacuolar pH stability, cation-homeostatic process and fungal fitness to host insect and environment.

Comparative Genomic Data Provide New Insight on the Evolution of Pathogenicity in Sporothrix Species

Sporothrix species are commonly isolated from environmental and clinical samples. As common causes of zoonotic mycosis, Sporothrix species may result in localized or disseminated infections, posing considerable threat to animal and human health. However, the pathogenic profiles of different Sporothrix species varied, in virulence, geographic location and host ranges, which have yet to be explored. Analysing the genomes of Sporothrix species are useful for understanding their pathogenicity. In this study, we analyzed the whole genome of 12 Sporothrix species and six S. globosa isolates from different clinical samples in China. By combining comparative analyses with Kyoto Encyclopedia of Genes and Genomes (KEGG), Carbohydrate-Active enZymes (CAZy), antiSMASH, Pfam, and PHI annotations, Sporothrix species showed exuberant primary and secondary metabolism processes. The genome sizes of four main clinical species, i.e., S. brasiliensis, S. schenckii, S. globosa, and S. luriei were significantly smaller than other environmental and clinical Sporothrix species. The contracted genes included mostly CAZymes and peptidases genes that were usually associated with the decay of plants, as well as the genes that were associated with the loss of pathogenicity and the reduced virulence. Our results could, to some extent, explain a habitat shift of Sporothrix species from a saprobic life in plant materials to a pathogenic life in mammals and the increased pathogenicity during the evolution. Gene clusters of melanin and clavaric acid were identified in this study, which improved our understanding on their pathogenicity and possible antitumor effects. Moreover, our analyses revealed no significant genomic variations among different clinical isolates of S. globosa from different regions in China.

Secondary metabolites from hypocrealean entomopathogenic fungi: novel bioactive compounds

Covering: 2014 up to the third quarter of 2019 Entomopathogens constitute a unique, specialized trophic subgroup of fungi, most of whose members belong to the order Hypocreales (class Sordariomycetes, phylum Ascomycota). These Hypocrealean Entomopathogenic Fungi (HEF) produce a large variety of secondary metabolites (SMs) and their genomes rank highly for the number of predicted, unique SM biosynthetic gene clusters. SMs from HEF have diverse roles in insect pathogenicity as virulence factors by modulating various interactions between the producer fungus and its insect host. In addition, these SMs also defend the carcass of the prey against opportunistic microbial invaders, mediate intra- and interspecies communication, and mitigate abiotic and biotic stresses. Thus, these SMs contribute to the role of HEF as commercial biopesticides in the context of integrated pest management systems, and provide lead compounds for the development of chemical pesticides for crop protection. These bioactive SMs also underpin the widespread use of certain HEF as nutraceuticals and traditional remedies, and allowed the modern pharmaceutical industry to repurpose some of these molecules as life-saving human medications. Herein, we survey the structures and biological activities of SMs described from HEF, and summarize new information on the roles of these metabolites in fungal virulence.

MrHex1 is Required for Woronin Body Formation, Fungal Development and Virulence in Metarhizium robertsii

The Woronin body (WB) is a peroxisome-derived dense-core vesicle, a self-assembling hexagonal crystal of a single protein Hex1. This organelle is specific to the ascomycete fungi belonging to the Pezizomycotina subphylum by functioning in sealing septal pores in response to mycelium damage and the control of cell heterogeneity. We retrieved all available Hex1-domain containing proteins of different fungi from the GenBank database and found considerable length variations among 460 obtained Hex1 proteins. However, a highly conserved Hex1 domain containing 75 amino acid residues with a specific S/A-R/S-L consensus motif for targeting peroxisome is present at the carboxy-terminus of each protein. A homologous Hex1 gene, named MrHex1, was deleted in the entomopathogenic fungus Metarhizium robertsii. It was found that MrHex1 was responsible for WB formation in M. robertsii and involved in sealing septal pores to maintain cell integrity and heterogeneity. Different assays indicated that, relative to the wild-type (WT) strain, ∆Mrhex1 demonstrated a growth defect on a solid medium and substantial reductions of conidiation, appressorium formation and topical infectivity against insect hosts. However, there was no obvious virulence difference between WT and mutants during injection of insects. We also found that ∆MrHex1 could tolerate different stress conditions like the WT and the gene-rescued mutant of M. robertsii, which is in contrast to the reports of the stress-response defects of the Hex1 null mutants of other fungal species. In addition to revealing the phenotypic/functional alterations of the Hex1 deletion mutants between different pathotype fungi, the results of this study may benefit the understanding of the evolution and WB-control of fungal entomopathogenicity.

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.

GATA-type transcription factor MrNsdD regulates dimorphic transition, conidiation, virulence and microsclerotium formation in the entomopathogenic fungus Metarhizium rileyi

The GATA-type sexual development transcription factor NsdD has been implicated in virulence, secondary metabolism and asexual development in filamentous fungi. However, little is known about its function in the yeast-to-hypha transition and in microsclerotium formation. In the current study, the orthologous NsdD gene MrNsdD in the entomopathogenic fungus Metarhizium rileyi was characterized. Transcriptional analysis indicated that MrNsdD was involved in yeast-to-hypha transition, conidiation and microsclerotium formation. After targeted deletion of MrNsdD, dimorphic transition, conidiation, fungal virulence and microsclerotium formation were all impaired. Compared with the wild-type strain, the ΔMrNsdD mutants were hypersensitive to thermal stress. Furthermore, transcriptome sequencing analysis revealed that MrNsdD regulated a distinct signalling pathway in M. rileyi during the yeast-to-hypha transition or microsclerotium formation, but exhibited overlapping regulation of genes during the two distinct developmental stages. Taken together, characterization of the MrNsdD targets in this study will aid in the dissection of the molecular mechanisms of dimorphic transition and microsclerotium development.

Congruence Amidst Discordance between Sequence and Protein-Content Based Phylogenies of Fungi

Amid the genomic data explosion, phylogenomic analysis has resolved the tree of life of different organisms, including fungi. Genome-wide clustering has also been conducted based on gene content data that can lighten the issue of the unequal evolutionary rate of genes. In this study, using different fungal species as models, we performed phylogenomic and protein-content (PC)-based clustering analysis. The obtained sequence tree reflects the phylogenetic trajectory of examined fungal species. However, 15 PC-based trees constructed from the Pfam matrices of the whole genomes, four protein families, and ten subcellular locations largely failed to resolve the speciation relationship of cross-phylum fungal species. However, lifestyle and taxonomic associations were more or less evident between closely related fungal species from PC-based trees. Pairwise congruence tests indicated that a varied level of congruent or discordant relationships were observed between sequence- and PC-based trees, and among PC-based trees. It was intriguing to find that a few protein family and subcellular PC-based trees were more topologically similar to the phylogenomic tree than was the whole genome PC-based phylogeny. In particular, a most significant level of congruence was observed between sequence- and cell wall PC-based trees. Cophylogenetic analysis conducted in this study may benefit the prediction of the magnitude of evolutionary conservation, interactive associations, or networking between different family or subcellular proteins.

Identification of a key G-protein coupled receptor in mediating appressorium formation and fungal virulence against insects

Fungal G-protein coupled receptors (GPCRs) play essential roles in sensing environmental cues including host signals. The study of GPCR in mediating fungus-insect interactions is still limited. Here we report the evolution of GPCR genes encoded in the entomopathogenic Metarhizium species and found the expansion of Pth11-like GPCRs in the generalist species with a wide host range. By deletion of ten candidate genes MrGpr1-MrGpr10 selected from the six obtained subfamilies in the generalist M. robertsii, we found that each of them played a varied level of roles in mediating appressorium formation. In particular, deletion of MrGpr8 resulted in the failure of appressorium formation on different substrates and the loss of virulence during topical infection of insects but not during injection assays when compared with the wild-type (WT) strain. Further analysis revealed that disruption of MrGpr8 substantially impaired the nucleus translocation of the mitogen-activated protein kinase (MAPK) Mero-Fus3 but not the MAPK Mero-Slt2 during appressorium formation. We also found that the defect of AMrGpr8 could not be rescued with the addition of cyclic AMP for appressorium formation. Relative to the WT, differential expression of the selected genes have also been detected in AMrGpr8. The results of this study may benefit the understanding of fungus-interactions mediated by GPCRs.

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.

Genome-wide mapping using new AFLP markers to explore intraspecific variation among pathogenic Sporothrix species

Sporotrichosis is a chronic subcutaneous mycosis caused by Sporothrix species, of which the main aetiological agents are S. brasiliensis, S. schenckii, and S. globosa. Infection occurs after a traumatic inoculation of Sporothrix propagules in mammals' skin and can follow either a classic route through traumatic inoculation by plant debris (e.g., S. schenckii and S. globosa) or an alternative route through zoonotic transmission from animals (e.g., S. brasiliensis). Epizootics followed by a zoonotic route occur in Brazil, with Rio de Janeiro as the epicenter of a recent cat-transmitted epidemic. DNA-based markers are needed to explore the epidemiology of these Sporothrix expansions using molecular methods. This paper reports the use of amplified-fragment-length polymorphisms (AFLP) to assess the degree of intraspecific variability among Sporothrix species. We used whole-genome sequences from Sporothrix species to generate 2,304 virtual AFLP fingerprints. In silico screening highlighted 6 primer pair combinations to be tested in vitro. The protocol was used to genotype 27 medically relevant Sporothrix. Based on the overall scored AFLP markers (97-137 fragments), the values of polymorphism information content (PIC = 0.2552-0.3113), marker index (MI = 0.002-0.0039), effective multiplex ratio (E = 17.8519-35.2222), resolving power (Rp = 33.6296-63.1852), discriminating power (D = 0.9291-0.9662), expected heterozygosity (H = 0.3003-0.3857), and mean heterozygosity (Havp = 0.0001) demonstrated the utility of these primer combinations for discriminating Sporothrix. AFLP markers revealed cryptic diversity in species previously thought to be the most prevalent clonal type, such as S. brasiliensis, responsible for cat-transmitted sporotrichosis, and S. globosa responsible for large sapronosis outbreaks in Asia. Three combinations (#3 EcoRI-FAM-GA/MseI-TT, #5 EcoRI-FAM-GA/MseI-AG, and #6 EcoRI-FAM-TA/MseI-AA) provide the best diversity indices and lowest error rates. These methods make it easier to track routes of disease transmission during epizooties and zoonosis, and our DNA fingerprint assay can be further transferred between laboratories to give insights into the ecology and evolution of pathogenic Sporothrix species and to inform management and mitigation strategies to tackle the advance of sporotrichosis.

Bioactive Metabolites and Potential Mycotoxins Produced by Cordyceps Fungi: A Review of Safety

Ascomycete Cordyceps fungi such as C. militaris, C. cicadae, and C. guangdongensis have been mass produced on artificial media either as food supplements or health additives while the byproducts of culture substrates are largely used as animal feed. The safety concerns associated with the daily consumption of Cordyceps fungi or related products are still being debated. On the one hand, the known compounds from these fungi such as adenosine analogs cordycepin and pentostatin have demonstrated different beneficial or pharmaceutical activities but also dose-dependent cytotoxicities, neurological toxicities and or toxicological effects in humans and animals. On the other hand, the possibility of mycotoxin production by Cordyceps fungi has not been completely ruled out. In contrast to a few metabolites identified, an array of biosynthetic gene clusters (BGCs) are encoded in each genome of these fungi with the potential to produce a plethora of as yet unknown secondary metabolites. Conservation analysis of BGCs suggests that mycotoxin analogs of PR-toxin and trichothecenes might be produced by Cordyceps fungi. Future elucidation of the compounds produced by these functionally unknown BGCs, and in-depth assessments of metabolite bioactivity and chemical safety, will not only facilitate the safe use of Cordyceps fungi as human food or alternative medicine, but will also benefit the use of mass production byproducts as animal feed. To corroborate the long record of use as a traditional medicine, future efforts will also benefit the exploration of Cordyceps fungi for pharmaceutical purposes.

Transcriptome profiling reveals insertional mutagenesis suppressed the expression of candidate pathogenicity genes in honeybee fungal pathogen, Ascosphaera apis

Chalkbrood disease is caused by Ascosphaera apis which severely affects honeybee brood. Spore inoculation experiments shown pathogenicity varies among different strains and mutants, however, the molecular mechanism of pathogenicity is unclear. We sequenced, assembled and annotated the transcriptomes of wild type (SPE1) and three mutants (SPE2, SPE3 and SPE4) with reduced pathogenicity that were constructed in our previous study. Illumina sequencing generated a total of 394,910,604 clean reads and de novo Trinity-based assembled into 12,989 unigenes, among these, 9,598 genes were successfully annotated to known proteins in UniProt database. A total of 172, 3,996, and 650 genes were up-regulated and 4,403, 2,845, and 3,016 genes were down-regulated between SPE2-SPE1, SPE3-SPE1, and SPE4-SPE1, respectively. Overall, several genes with a potential role in fungal pathogenicity were detected down-regulated in mutants including 100 hydrolytic enzymes, 117 transcriptional factors, and 47 cell wall related genes. KEGG pathway enrichment analysis reveals 216 genes involved in nine pathways were down-regulated in mutants compared to wild type. The down-regulation of more pathways involved in pathogenicity in SPE2 and SPE4 than SPE3 supports their lower pathogenicity during in-vitro bioassay experiment. Expression of 12 down-regulated genes in mutants was validated by quantitative real time PCR. This study provides valuable information on transcriptome variation caused by mutation for further functional validation of candidate pathogenicity genes in A. apis.

Pathogenomics and Management of Fusarium Diseases in Plants

There is an urgency to supplant the heavy reliance on chemical control of Fusarium diseases in different economically important, staple food crops due to development of resistance in the pathogen population, the high cost of production to the risk-averse grower, and the concomitant environmental impacts. Pathogenomics has enabled (i) the creation of genetic inventories which identify those putative genes, regulators, and effectors that are associated with virulence, pathogenicity, and primary and secondary metabolism; (ii) comparison of such genes among related pathogens; (iii) identification of potential genetic targets for chemical control; and (iv) better characterization of the complex dynamics of host–microbe interactions that lead to disease. This type of genomic data serves to inform host-induced gene silencing (HIGS) technology for targeted disruption of transcription of select genes for the control of Fusarium diseases. This review discusses the various repositories and browser access points for comparison of genomic data, the strategies for identification and selection of pathogenicity- and virulence-associated genes and effectors in different Fusarium species, HIGS and successful Fusarium disease control trials with a consideration of loss of RNAi, off-target effects, and future challenges in applying HIGS for management of Fusarium diseases.

Analogous and Diverse Functions of APSES-Type Transcription Factors in the Morphogenesis of the Entomopathogenic Fungus Metarhizium rileyi

APSES-type transcription factors (TFs) have analogous and diverse functions in the regulation of fungal morphogenesis processes. However, little is known about these functions in microsclerotium formation. In this study, we characterized two orthologous APSES genes (MrStuA and MrXbp) in the entomopathogenic fungus Metarhizium rileyi Deletion of either MrStuA or MrXbp impaired dimorphic transition, conidiation, fungal virulence, and microsclerotium formation. Compared with the wild-type strain, ΔMrStuA and ΔMrXbp mutants were hypersensitive to thermal and oxidative stress. Furthermore, transcriptome sequencing analysis revealed that MrStuA and MrXbp independently regulate their own distinctive subsets of signaling pathways during dimorphic transition and microsclerotium formation, but they also show an overlapping regulation of genes during these two distinct morphogenesis processes. These results provide a global insight into vital roles of MrStuA and MrXbp in M. rileyi and aid in dissection of the interacting regulatory mechanisms of dimorphism transition and microsclerotium development.IMPORTANCE Transcription factors (TFs) are core components of the signaling pathway and play an important role in transcriptional regulation of gene expression during fungal morphogenesis processes. A prevailing theory suggests an interplay between different TFs regulating microsclerotial differentiation; however, the persisting issue remains that these interplay mechanisms are not clear. Here, we analyzed two members of the APSES-type TFs in Metarhizium rileyi using a gene deletion strategy and transcriptome analysis. Mutants were significantly impaired in microsclerotium formation and dimorphic transition. Transcriptome analysis provided evidence for interacting regulatory mechanisms by the two TFs in microsclerotium formation and dimorphic transition. Furthermore, we investigated their overlapping roles in mediating the expression of genes required for different fungal morphogenesis processes. Characterization of TFs in this study will aid in dissecting the interplay between regulatory mechanisms in fungal morphogenesis processes.

Seeking the Roles for Fungal Small-Secreted Proteins in Affecting Saprophytic Lifestyles

Small secreted proteins (SSPs) comprise 40–60% of the total fungal secretome and are present in fungi of all phylogenetic groups, representing the entire spectrum of lifestyles. They are characteristically shorter than 300 amino acids in length and have a signal peptide. The majority of SSPs are coded by orphan genes, which lack known domains or similarities to known protein sequences. Effectors are a group of SSPs that have been investigated extensively in fungi that interact with living hosts, either pathogens or mutualistic systems. They are involved in suppressing the host defense response and altering its physiology. Here, we aim to delineate some of the potential roles of SSPs in saprotrophic fungi, that have been bioinformatically predicted as effectors, and termed in this mini-review as “effector-like” proteins. The effector-like Ssp1 from the white-rot fungus Pleurotus ostreatus is presented as a case study, and its potential role in regulating the ligninolytic system, secondary metabolism, development, and fruiting body initiation are discussed. We propose that deciphering the nature of effector-like SSPs will contribute to our understanding of development and communication in saprophytic fungi, as well as help, to elucidate the origin, regulation, and mechanisms of fungal-host, fungal-fungal, and fungal-bacterial interactions.

Genome-wide analysis reveals molecular convergence underlying domestication in 7 bird and mammals

Background

In response to ecological niche of domestication, domesticated mammals and birds developed adaptively phenotypic homoplasy in behavior modifications like fearlessness, altered sociability, exploration and cognition, which partly or indirectly result in consequences for economic productivity. Such independent adaptations provide an excellent model to investigate molecular mechanisms and patterns of evolutionary convergence driven by artificial selection.

Results

First performing population genomic and brain transcriptional comparisons in 68 wild and domesticated chickens, we revealed evolutionary trajectories, genetic architectures and physiologic bases of adaptively behavioral alterations. To extensively decipher molecular convergence on behavioral changes thanks to domestication, we investigated selection signatures in hundreds of genomes and brain transcriptomes across chicken and 6 other domesticated mammals. Although no shared substitution was detected, a common enrichment of the adaptive mutations in regulatory sequences was observed, presenting significance to drive adaptations. Strong convergent pattern emerged at levels of gene, gene family, pathway and network. Genes implicated in neurotransmission, semaphorin, tectonic protein and modules regulating neuroplasticity were central focus of selection, supporting molecular repeatability of homoplastic behavior reshapes. Genes at nodal positions in trans-regulatory networks were preferably targeted. Consistent down-regulation of majority brain genes may be correlated with reduced brain size during domestication. Up-regulation of splicesome genes in chicken rather mammals highlights splicing as an efficient way to evolve since avian-specific genomic contraction of introns and intergenics. Genetic burden of domestication elicits a general hallmark. The commonly selected genes were relatively evolutionary conserved and associated with analogous neuropsychiatric disorders in human, revealing trade-off between adaption to life with human at the cost of neural changes affecting fitness in wild.

Conclusions

After a comprehensive investigation on genomic diversity and evolutionary trajectories in chickens, we revealed basis, pattern and evolutionary significance of molecular convergence in domesticated bird and mammals, highlighted the genetic basis of a compromise on utmost adaptation to the lives with human at the cost of high risk of neurophysiological changes affecting animals’ fitness in wild.

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.

Unveiling the function and regulation control of the DUF3129 family proteins in fungal infection of hosts

Many prokaryotic and eukaryotic proteins contain domains of unknown function (DUFs). A DUF3129 family of proteins is widely encoded in the genomes of fungal pathogens. A few studies in plant and insect pathogens indicated that the DUF3129 genes are required for fungal penetration of host cuticles with an unclear mechanism. We found that a varied number of DUF3129 proteins is present in different fungal species and the proteins are evolutionarily diverged from each other at the inter- and intra-specific levels. By using the insect pathogenic fungus Metarhizium robertsii as a model, we performed experiments and found that the seven DUF3129 proteins encoded by this fungus are localized to cellular lipid droplets (LDs). Individual deletion of these genes did not affect fungal formation of the infection structure appressoria and the accumulation of LDs in fungal conidia. When compared with the wild-type (WT) strain, insect bioassays revealed that the virulence of most null mutants were significantly impaired during topical infection but not during injection of insects. Carbon starvation and the subsequent Western blot analysis indicated that the LD-specific perilipin protein was completely degraded in the WT cells whereas varied levels of perilipin could be detected in the mutant cells, which signified that depletion of LD content was delayed in mutant cells, and DUF3129 proteins are therefore involved in LD degradation. We also provided biochemical evidence that these DUF3129 genes are transcriptionally regulated by a yeast Ste12-like transcription factor. The findings of this study not only unveil the function of DUF3129 proteins but also better understand the diverse mechanism of fungus-host interactions. This article is part of the theme issue 'Biotic signalling sheds light on smart pest management'.

Complete mitogenome of the entomopathogenic fungus Akanthomyces lecanii

In this study, the complete mitogenome of an entomopathogenic fungus Akanthomyces lecanii is assembled and annotated. This circular mitogenome is 24,643 bp in length and consists of 2 rRNA genes (rnl and rns), 26 tRNA genes and 14 standard protein-coding genes of the oxidative phosphorylation system. Only one intron (group IA) is identified, which invades rnl and carries an ORF coding for ribosomal protein S3. Phylogenetic analysis based on concatenated mitochondrial nucleotide sequences confirms A. lecanii in Cordycipitaceae, and A. lecanii clusters together with Akanthomyces muscarius.

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.

Optimizing cultivation of Cordyceps militaris for fast growth and cordycepin overproduction using rational design of synthetic media

Cordyceps militaris is an entomopathogenic fungus which is often used in Asia as a traditional medicine developed from age-old wisdom. Presently, cordycepin from C. militaris is a great interest in medicinal applications. However, cellular growth of C. militaris and the association with cordycepin production remain poorly understood. To explore the metabolism of C. militaris as potential cell factories in medical and biotechnology applications, this study developed a high-quality genome-scale metabolic model of C. militaris, iNR1329, based on its genomic content and physiological data. The model included a total of 1329 genes, 1821 biochemical reactions, and 1171 metabolites among 4 different cellular compartments. Its in silico growth simulation results agreed well with experimental data on different carbon sources. iNR1329 was further used for optimizing the growth and cordycepin overproduction using a novel approach, POPCORN, for rational design of synthetic media. In addition to the high-quality GEM iNR1329, the presented POPCORN approach was successfully used to rationally design an optimal synthetic medium with C:N ratio of 8:1 for enhancing 3.5-fold increase in cordycepin production. This study thus provides a novel insight into C. militaris physiology and highlights a potential GEM-driven method for synthetic media design and metabolic engineering application. The iNR1329 and the POPCORN approach are available at the GitHub repository: https://github.com/sysbiomics/Cordyceps_militaris-GEM.

MrMid2, encoding a cell wall stress sensor protein, is required for conidium production, stress tolerance, microsclerotium formation and virulence in the entomopathogenic fungus Metarhizium rileyi

Transmembrane proteins as sensors encoded by fungal genes activate specific intracellular signal pathways in response to stress cues to help the fungus survive in a changing environment. To better understand the role of the cell wall integrity (CWI) pathway in the entomopathogenic fungus Metarhizium rileyi, an ortholog encoding the transmembrane protein Mid2, MrMid2, was identified and characterized functionally. Transcriptional analysis indicated that MrMid2 was involved in dimorphic transition, conidiation, and microsclerotium formation. After a targeted deletion of MrMid2, all three traits were impaired. Compared with the wild-type strain, the △MrMid2 mutants were hypersensitive to thermal stress, and cell wall and oxidative stress. Insect bioassays revealed that △MrMid2 mutants had decreased virulence levels following topical (22.5%) and injection bioassays (38.7%). Furthermore, transcription analysis showed that other genes of the CWI pathway, with the exception of another major sensor protein encoding gene, MrWsc1, were down-regulated in △MrMid2 mutants. These results suggest that MrMid2 plays important roles in dimorphic transition, conidiation, the stress response, virulence, and microsclerotium development in M. rileyi.

Genetic and Pathogenic Variability of Mycogone perniciosa Isolates Causing Wet Bubble Disease on Agaricus bisporus in China

Wet bubble disease, caused by Mycogone perniciosa, is a major threat to Agaricus bisporus production in China. In order to understand the variability in genetic, pathogenicity, morphology, and symptom production of the fungus, 18 isolates of the pathogen were collected from diseased A. bisporus in different provinces in China. The isolates were characterized by a combination of morphological, cultural, molecular, and pathogenicity testing on different strains of A. bisporus and amplified fragment length polymorphism (AFLP) analysis. The 18 isolates were identified by Koch’s postulate and confirmed different pathogenic variability among them. The yellow to brown isolates were more virulent than the white isolates. AFLP markers clustered the isolates into two distinct groups based on their colony color, with a high level of polymorphism of Jaccard similarities range from 0.39% to 0.64%. However, there was no evidence of an association between the genetic diversity and the geographical origin of the isolates. Through knowledge of the genetic diversity, phenotypic virulence of M. perniciosa is a key factor for successful breeding of resistant strains of A. bisporus and developing of an integrated disease management strategy to manage wet bubble disease of A. bisporus.

Genome Sequence of Metarhizium rileyi, a Microbial Control Agent for Lepidoptera

Metarhizium rileyi (formerly known as Nomuraea rileyi) is a potential agent for microbial control of many insect pests from the order Lepidoptera, the damages of which can cause considerable loss of productivity in agriculture. We report the genome sequence and annotation of M. rileyi strain Cep018-CH2/ARSEF 7053.

Metarhizium rileyi (formerly known as Nomuraea rileyi) is a potential agent for microbial control of many insect pests from the order Lepidoptera, the damages of which can cause considerable loss of productivity in agriculture. We report the genome sequence and annotation of M. rileyi strain Cep018-CH2/ARSEF 7053.

Methylglucosylation of Phenolic Compounds by Fungal Glycosyltransferase-Methyltransferase Functional Modules

Glycosylation endows both natural and synthetic small molecules with modulated physicochemical and biological properties. Plant and bacterial glycosyltransferases capable of decorating various privileged scaffolds have been extensively studied, but those from kingdom Fungi still remain underexploited. Here, we use a combination of genome mining and heterologous expression techniques to identify four novel glycosyltransferase-methyltransferase (GT-MT) functional modules from Hypocreales fungi. These GT-MT modules display decent substrate promiscuity and regiospecificity, methylglucosylating a panel of natural products such as flavonoids, stilbenoids, anthraquinones, and benzenediol lactones. Native GT-MT modules can be split up and regrouped into hybrid modules with similar or even improved efficacy as compared with native pairs. Methylglucosylation of kaempferol considerably improves its insecticidal activity against the larvae of oriental armyworm Mythimna separata (Walker). Our work provides a set of efficient biocatalysts for the combinatorial biosynthesis of small molecule glycosides that may have significant importance to the pharmaceutical, agricultural, and food industries.

Fungal evolution: major ecological adaptations and evolutionary transitions

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

How Environmental Fungi Cause a Range of Clinical Outcomes in Susceptible Hosts

Environmental fungi are globally ubiquitous and human exposure is near universal. However, relatively few fungal species are capable of infecting humans, and among fungi, few exposure events lead to severe systemic infections. Systemic infections have mortality rates of up to 90%, cost the US healthcare system $7.2 billion annually, and are typically associated with immunocompromised patients. Despite this reputation, exposure to environmental fungi results in a range of outcomes, from asymptomatic latent infections to severe systemic infection. Here we discuss different exposure outcomes for five major fungal pathogens: Aspergillus, Blastomyces, Coccidioides, Cryptococcus, and Histoplasma species. These fungi include a mold, a budding yeast, and thermal dimorphic fungi. All of these species must adapt to dramatically changing environments over the course of disease. These dynamic environments include the human lung, which is the first exposure site for these organisms. Fungi must defend themselves against host immune cells while germinating and growing, which risks further exposing microbe-associated molecular patterns to the host. We discuss immune evasion strategies during early infection, from disruption of host immune cells to major changes in fungal cell morphology.

Genome sequence of Isaria javanica and comparative genome analysis insights into family S53 peptidase evolution in fungal entomopathogens

The fungus Isaria javanica is an important entomopathogen that parasitizes various insects and is effective for pest control. In this study, we sequenced and assembled the genomes (IJ1G and IJ2G) of two I. javanica strains isolated from different insects. The genomes were approximately 35 Mb in size with 11,441 and 11,143 protein-coding genes, respectively. Using a phylogenomic approach, we evaluated genome evolution across five entomopathogenic fungi in Cordycipitaceae. By comparative genome analysis, it was found that family S53 serine peptidases were expanded in Cordycipitaceae entomopathogens, particularly in I. javanica. Gene duplication events were identified based on phylogenetic relationships inferred from 82 S53 peptidases within six entomopathogenic fungal genomes. Moreover, we found that carbohydrate-active enzymes and proteinases were the largest secretory protein groups encoded in the I. javanica genome, especially chitinases (GH18), serine and aspartic peptidases (S53, S08, S10, A01). Pathogenesis-related genes and genes for bacterial-like toxins and secondary metabolites were also identified. By comparative transcriptome analysis, differentially expressed genes in response to insect nutrients (in vitro) were identified. Moreover, most S53 peptidases were detected to be significantly upregulated during the initial fungal infection process in insects (in vivo) by RT-qPCR. Our results provide new clues about understanding evolution of pathogenic proteases and may suggest that abundant S53 peptidases in the I. javanica genome may contribute to its effective parasitism on various insects.

Complete mitogenome of the entomopathogenic fungus Sporothrix insectorum RCEF 264 and comparative mitogenomics in Ophiostomatales

The fungal order Ophiostomatales contains numerous species important in medical fields, agriculture, and forestry, and several species have had available mitogenome information. The nuclear genome of the entomopathogenic fungus Sporothrix insectorum has been reported, while its mitogenome remains unknown. Herein, we firstly described the mitogenome of S. insectorum RCEF 264 and then compared Ophiostomatales mitogenomes from both interspecific and intraspecific perspectives. The mitogenome of S. insectorum RCEF 264 was 31,454 bp in length, containing typical fungal mitochondrial genes plus rnpB. Four group I introns interrupted rnl and cox1. Phylogenetic analyses confirmed the placement of S. insectorum RCEF 264 in Ophiostomatales. Comparison of mitogenomes among seven Ophiostomatales species revealed conserved gene contents and a high synteny, although there were also some differences among them. Their mitogenomes showed more than two-fold variations (26.6-65.1 kb) in size, with a total of 37 intron insertional loci from 11 genes (1-25 introns per species). The sole intron shared by all species was an rps3-encoding intron in rnl (mL2450), and this intron-based phylogeny was highly consistent with those constructed using mitochondrial/nuclear genes, suggesting convergent evolution of this intron with Ophiostomatales species. The dendrogram based on presence/absence patterns at all intron loci was quite different from those based on mitochondrial/nuclear genes. Comparison of mitogenomes among two to three intraspecific individuals in Ophiostoma novo-ulmi subsp. novo-ulmi and Sporothrix schenckii revealed mitogenome size variations due to single-nucleotide polymorphisms (SNPs) and indels but without fluctuation of intron numbers for each species. This study greatly enhanced our understanding of mitogenome evolution in Ophiostomatales.