Genome-Wide Host-Pathogen Interaction Unveiled by Transcriptomic Response of Diamondback Moth to Fungal Infection

Time-series RNA sequencing reveals that immulectin-1 (IML-1) is a potential target for biological control of fall armyworm, Spodoptera frugiperda

Insects mount robust immune responses to reduce the efficacy of biopesticides. Identification of critical immune targets for biocontrol strategies is thus essential. However, immune targets in the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae), a major pest of corn, remain largely unexplored. Here, we investigated the immune response of S. frugiperda to the entomopathogen Beauveria bassiana (Hypocreales: Cordycipitaceae). Time-series RNA sequencing revealed 863, 717, and 165 differentially expressed genes (DEGs) at 24, 48, and 60 h post-infection, respectively. GO and KEGG enrichment analyses revealed that DEGs were significantly enriched in immune-related terms and pathways. Among the DEGs, the gene encoding immulectin SfIML-1, a C-type lectin with dual carbohydrate recognition domains, was significantly upregulated during the early and intermediate stages of infection. Genome-wide analysis identified nine IML genes, with SfIML-1 as the only one highly responsive to B. bassiana. Recombinant SfIML-1 binds pathogen-associated molecular patterns and pathogens, enhancing cellular responses. RNAi of SfIML-1 decreased survival after fungal infection, while injection of recombinant SfIML-1 increased survival. Suppression of the NF-κB and JAK-STAT pathways diminished SfIML-1 expression in both fat body and Sf9 cells. These results identify SfIML-1 as a key pattern recognition receptor and a promising target for developing novel biocontrol strategies.

Genotype-environment associations reveal genes potentially linked to avian malaria infection in populations of an endemic island bird

Patterns of pathogen prevalence are, at least partially, the result of coevolutionary host-pathogen interactions. Thus, exploring the distribution of host genetic variation in relation to infection by a pathogen within and across populations can provide important insights into mechanisms of host defence and adaptation. Here, we use a landscape genomics approach (Bayenv) in conjunction with genome-wide data (ddRADseq) to test for associations between avian malaria (Plasmodium) prevalence and host genetic variation across 13 populations of the island endemic Berthelot's pipit (Anthus berthelotii). Considerable and consistent spatial heterogeneity in malaria prevalence was observed among populations over a period of 15 years. The prevalence of malaria infection was also strongly positively correlated with pox (Avipoxvirus) prevalence. Multiple host loci showed significant associations with malaria prevalence after controlling for genome-wide neutral genetic structure. These sites were located near to or within genes linked to metabolism, stress response, transcriptional regulation, complement activity and the inflammatory response, many previously implicated in vertebrate responses to malarial infection. Our findings identify diverse genes - not just limited to the immune system - that may be involved in host protection against malaria and suggest that spatially variable pathogen pressure may be an important evolutionary driver of genetic divergence among wild animal populations, such as Berthelot's pipit. Furthermore, our data indicate that spatio-temporal variation in multiple different pathogens (e.g. malaria and pox in this case) may have to be studied together to develop a more holistic understanding of host pathogen-mediated evolution.

Transcriptomic landscape of the interaction between the entomopathogenic fungus Beauveria bassiana and its tolerant host Tribolium castaneum revealed by dual RNA-seq

Entomopathogenic fungi such as Beauveria bassiana are the only insect pathogens able to start the infection process by penetrating through the host cuticle. However, some insects try to avoid fungal infection by embedding their cuticle with antifungal compounds. This is the case of the red flour beetle Tribolium castaneum, which generates economical loss of great significance in stored product environments worldwide. In this study, T. castaneum adults were fed during different time periods (from 3 to 72 h) on B. bassiana conidia-covered corn kernels. The progression of fungal infection was monitored using the dual RNA-seq technique to reconstruct the temporal transcriptomic profile and to perform gene enrichment analyses in both interacting organisms. After mapping the total reads with the B. bassiana genome, 904 genes were identified during this process. The more expressed fungal genes were related to carbon catabolite repression, cation binding, peptidase inhibition, redox processes, and stress response. Several immune-related genes from Toll, IMD, and JNK pathways, as well as genes related to chitin modification, were found to be differentially expressed in fungus-exposed T. castaneum. This study represents the first dual transcriptomic approach to help understand the interaction between the entomopathogenic fungus B. bassiana and its tolerant host T. castaneum.

Dual RNA Sequencing of Beauveria bassiana-Infected Spodoptera frugiperda Reveals a Fungal Protease with Entomopathogenic and Antiphytopathogenic Activities

Insect pests and phytopathogens significantly impact crop yield and quality. The fall armyworm (FAW) Spodoptera frugiperda and the phytopathogen Fusarium graminearum cause substantial economic losses in crops like barley and wheat. However, the entomopathogen Beauveria bassiana shows limited efficacy against FAW, and its antiphytopathogenic activities against F. graminearum remain unclear. Here, dual RNA sequencing was performed to identify differentially expressed genes in B. bassiana-infected FAW larvae. We found that the BbAorsin gene was significantly upregulated at 36 and 48 h post-infection. BbAorsin encodes a serine-carboxyl protease and is mainly expressed in blastospores and hyphae. Overexpression of BbAorsin in B. bassiana ARSEF2860 enhanced virulence against Galleria mellonella and FAW larvae and inhibited F. graminearum growth. The recombinant BbAorsin protein induced apoptosis and necrosis in FAW hemocytes and inhibited F. graminearum spore germination. These findings shed light on transcriptomic mechanisms governing insect-pathogen interactions, which could aid in developing dual-functional entomopathogens and anti-phytopathogens.

Transcriptomic Analysis Reveals the Impact of the Biopesticide Metarhizium anisopliae on the Immune System of Major Workers in Solenopsis invicta

The red imported fire ant (Solenopsis invicta Buren, 1972) is a globally significant invasive species, causing extensive agricultural, human health, and biodiversity damage amounting to billions of dollars worldwide. The pathogenic fungus Metarhizium anisopliae (Metchnikoff) Sorokin (1883), widely distributed in natural environments, has been used to control S. invicta populations. However, the interaction between M. anisopliae and the immune system of the social insect S. invicta remains poorly understood. In this study, we employed RNA-seq to investigate the effects of M. anisopliae on the immune systems of S. invicta at different time points (0, 6, 24, and 48 h). A total of 1313 differentially expressed genes (DEGs) were identified and classified into 12 expression profiles using short time-series expression miner (STEM) for analysis. Weighted gene co-expression network analysis (WGCNA) was employed to partition all genes into 21 gene modules. Upon analyzing the statistically significant WGCNA model and conducting Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis on the modules, we identified key immune pathways, including the Toll and Imd signaling pathways, lysosomes, autophagy, and phagosomes, which may collectively contribute to S. invicta defense against M. anisopliae infection. Subsequently, we conducted a comprehensive scan of all differentially expressed genes and identified 33 immune-related genes, encompassing various aspects such as recognition, signal transduction, and effector gene expression. Furthermore, by integrating the significant gene modules derived from the WGCNA analysis, we constructed illustrative pathway diagrams depicting the Toll and Imd signaling pathways. Overall, our research findings demonstrated that M. anisopliae suppressed the immune response of S. invicta during the early stages while stimulating its immune response at later stages, making it a potential biopesticide for controlling S. invicta populations. These discoveries lay the foundation for further understanding the immune mechanisms of S. invicta and the molecular mechanisms underlying its response to M. anisopliae.

Knockdown of CYP9A9 increases the susceptibility to lufenuron, methoxyfenozide and a mixture of both in Spodoptera exigua

Lufenuron (LUF) and Methoxyfenozide (MET) as Insect Growth Regulators (IGRs) contribute to the current control of the catastrophic crop pest, Spodoptera exigua (Lepidoptera, Noctuidae). Yet S. exigua has evolved resistance to LUF and MET, which is possibly mediated by cytochrome P450 monooxygenases (P450s), particularly from the CYP3 clade family, as it plays a key role in the detoxification of insecticides. However, a mixture of LUF and MET (MML) (optimal ratio: 6:4) remains highly insecticidal. Here, we analysed the response of S. exigua to sublethal concentrations of LUF, MET, and MML via transcriptomics. Twelve differentially expressed genes (DEGs) encoding CYP3 clade members were observed in transcriptomes and CYP9A9 was significantly upregulated after treatment with LUF, MET, and MML. Further, CYP9A9 was most highly expressed in the midgut of L4 S. exigua larvae. RNAi-mediated knockdown of CYP9A9 reduced the activity of CYP450 and increased the susceptibility of S. exigua larvae to LUF, MET, and MML. Thus, CYP9A9 plays a key role in the detoxification of LUF, MET, and MML in S. exigua. These findings provide new insights into insecticidal actions of IGRs, which can be applied to the establishment of novel pest management strategies.

Characterization and transcriptomic analysis of a native fungal pathogen against the rice pest Nilaparvata lugens

The brown planthopper (BPH), Nilaparvata lugens, is one of the most destructive pests of rice. Given the threats posed by insecticide resistance to its control, eco-friendly strategies based on microbial pathogens emerged as a promising biocontrol alternative. In the present study, we isolated a native fungal pathogen against BPH from infected BPH cadavers and preliminarily identified as a strain of Aspergillus fumigatus based on morphological and molecular methods. Laboratory bioassay revealed that this fungal strain was highly virulent to BPH both at nymphal and adult stages, with the median lethal times (LT₅₀) of 7.5 and 5.8 days under high conidial concentration of 1 × 10⁹ conidia mL^-1. A genome-wide view of gene expressions in BPH against fungal attack was analyzed by transcriptomic sequencing and consequently a large number of differentially expressed genes that mainly involved in host immune defense and cell detoxification were found. RNAi-mediated knockdown of an upregulated gene encoding a serine protease (NlSPN) could cause a significant decrease in BPH survival. Combination of dsRNA injection and fungal infection showed an additive effect on BPH mortality, which provided clues to develop new pest management strategies against BPH.

Transcriptomic analysis of entomopathogenic fungus Beauveria bassiana infected by a hypervirulent polymycovirus BbPmV-4

Polymycoviridae is a recently established family of mycoviruses. Beauveria bassiana polymycovirus 4 (BbPmV-4) was previously reported. However, the effect of the virus on host fungus B. bassiana was not clarified. Here, a comparison between virus-free and virus-infected isogenic lines of B. bassiana revealed that BbPmV-4 infection of B. bassiana changes morphology and could lead to decreases in conidiation and increases in virulence against Ostrinia furnacalis larvae. The differential expression of genes between virus-free and virus-infected strains was compared by RNA-Seq and was consistent with the phenotype of B. bassiana. The enhanced pathogenicity may be related to the significant up-regulation of genes encoding mitogen activated protein kinase, cytochrome P450, and polyketide synthase. The results enable studies of the mechanism of interaction between BbPmV-4 and B. bassiana.

Three Chitin Deacetylase Family Members of Beauveria bassiana Modulate Asexual Reproduction and Virulence of Fungi by Mediating Chitin Metabolism and Affect Fungal Parasitism and Saprophytic Life

As an important chitin-modifying enzyme, chitin deacetylase (CDA) has been characterized in many fungi, but its function in the entomopathogenic fungus Beauveria bassiana remains unclear. Three CDAs with conserved domains of the carbohydrate esterase 4 (CE-4) family were identified in B. bassiana. Disruption of CDA1 resulted in growth restriction of the fungus on medium with chitin as a carbon source or without a carbon source. Deletion of CDA1 and CDA2 led to defects in fungal conidial formation and conidial vitality compared with those of the wild type (WT), and the conidial yield decreased by 25.81% to 47.68%. Inactivation of three CDA genes resulted in a decrease of 20.23% to 27% in the blastospore yield. ΔCDA1 and ΔCDA3 showed 29.33% and 23.34% reductions in cuticular infection virulence, respectively. However, the CDA family may not contribute to hemocoel infection virulence. Additionally, the sporulation of the insect carcass showed that the three gene deletion mutants were 68.45%, 63.84%, and 56.65% less than WT. Penetration experiments with cicada wings and enzyme activity assays were used to further explore the effect of the fungus on chitin metabolism after gene deletion. Although the three gene deletion mutants penetrated the cicada wings successfully and continued to grow on the underlying medium, their colony sizes were reduced by 29.12% to 47.76%. The CDA enzyme activity of ΔCDA1 and ΔCDA3 decreased by 84.76% and 83.04%, respectively. These data showed that members of the CDA family play a different role in fungal growth, conidial quality, and virulence. IMPORTANCE In this study, we report the roles of CDA family in entomopathogenic fungus B. bassiana. Our results indicated that CDA modulates asexual development and regulates fungal virulence by altering chitin deacetylation and metabolic capacity. CDA affected the biological control potential and life history of B. bassiana by affecting its parasitic and saprophytic life. These findings provide novel insights into the roles of multiple CDA paralogues existing in fungal biocontrol agents.

Time-series gene expression patterns and their characteristics of Beauveria bassiana in the process of infecting pest insects

Beauveria bassiana has been widely used as an important biological control fungus for agricultural and forest pests, and clarifying the interaction mechanism between B. bassiana and its host will help to better exert the efficacy of the mycoinsecticide. Here, we proposed a novel pattern analysis (PA) method for analyzing time-series data and applied it to a transcriptomic data set of B. bassiana infecting Galleria mellonella. We screened out 14 patterns including 868 genes, which had some characteristics that were not inferior to differentially expressed genes (DEGs). Compared with the previous analysis of this data set, we had three novel discoveries during B. bassiana infection, including overall downregulation of gene expression, the more critical first 24 h, and enrichment of regulatory functions of downregulated genes. Our new PA method promises to be an important complement to DEGs analysis for time-series transcriptomic data, and our findings enrich our knowledge of molecular mechanisms of fungal-host interactions.

Fungal insecticidal activity elevated by non-risky markerless overexpression of an endogenous cysteine-free protein gene in Beauveria bassiana

BACKGROUND

Fungal insecticides are notorious for slow kill action, an intrinsic trait that can be improved by the genetic engineering of an exogenous or endogenous virulence factor. However, transgenic insecticides expressing exogenous toxin and herbicide-resistant marker genes may cause unexpected ecological risks and are hardly permitted for field release due to strict regulatory hurdles. It is necessary to improve biotechnology that can speed up fungal insect-killing action and exclude ecological risk source.

RESULTS

A markerless transformation system of Beauveria bassiana, a main source of wide-spectrum fungal insecticides, was reconstructed based on the fungal uridine auxotrophy (Δura3). The system was applied for overexpression of the small cysteine-free protein (120 amino acids) gene cfp previously characterized as a regulator of the fungal virulence and gene expression. Three cfp-overexpressed strains showed much faster kill action to Galleria mellonella larvae than the parental wild-type via normal cuticle infection but no change in vegetative growth and aerial condition. The faster kill action was achieved due to not only significant increases in conidial adherence to insect cuticle and total activity of secreted cuticle-degrading Pr1 proteases and of antioxidant enzymes crucial for collapse of insect immune defense but acceleration of hemocoel localization, proliferation in vivo and host death from mummification.

CONCLUSION

The markerless system is free of any foreign DNA fragment as a source of ecologic risk and provides a novel biotechnological approach to enhancing fungal insecticidal activity with non-risky endogenous genes and striding over the regulatory hurdles. © 2022 Society of Chemical Industry.

A high-quality de novo genome assembly based on nanopore sequencing of a wild-caught coconut rhinoceros beetle (Oryctes rhinoceros)

BACKGROUND

An optimal starting point for relating genome function to organismal biology is a high-quality nuclear genome assembly, and long-read sequencing is revolutionizing the production of this genomic resource in insects. Despite this, nuclear genome assemblies have been under-represented for agricultural insect pests, particularly from the order Coleoptera. Here we present a de novo genome assembly and structural annotation for the coconut rhinoceros beetle, Oryctes rhinoceros (Coleoptera: Scarabaeidae), based on Oxford Nanopore Technologies (ONT) long-read data generated from a wild-caught female, as well as the assembly process that also led to the recovery of the complete circular genome assemblies of the beetle's mitochondrial genome and that of the biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). As an invasive pest of palm trees, O. rhinoceros is undergoing an expansion in its range across the Pacific Islands, requiring new approaches to management that may include strategies facilitated by genome assembly and annotation.

RESULTS

High-quality DNA isolated from an adult female was used to create four ONT libraries that were sequenced using four MinION flow cells, producing a total of 27.2 Gb of high-quality long-read sequences. We employed an iterative assembly process and polishing with one lane of high-accuracy Illumina reads, obtaining a final size of the assembly of 377.36 Mb that had high contiguity (fragment N50 length = 12 Mb) and accuracy, as evidenced by the exceptionally high completeness of the benchmarked set of conserved single-copy orthologous genes (BUSCO completeness = 99.1%). These quality metrics place our assembly ahead of the published Coleopteran genomes, including that of an insect model, the red flour beetle (Tribolium castaneum). The structural annotation of the nuclear genome assembly contained a highly-accurate set of 16,371 protein-coding genes, with only 2.8% missing BUSCOs, and the expected number of non-coding RNAs. The number and structure of paralogous genes in a gene family like Sigma GST is lower than in another scarab beetle (Onthophagus taurus), but higher than in the red flour beetle (Tribolium castaneum), which suggests expansion of this GST class in Scarabaeidae. The quality of our gene models was also confirmed with the correct placement of O. rhinoceros among other members of the rhinoceros beetles (subfamily Dynastinae) in a phylogeny based on the sequences of 95 protein-coding genes in 373 beetle species from all major lineages of Coleoptera. Finally, we provide a list of 30 candidate dsRNA targets whose orthologs have been experimentally validated as highly effective targets for RNAi-based control of several beetles.

CONCLUSIONS

The genomic resources produced in this study form a foundation for further functional genetic research and management programs that may inform the control and surveillance of O. rhinoceros populations, and we demonstrate the efficacy of de novo genome assembly using long-read ONT data from a single field-caught insect.

Three Small Cysteine-Free Proteins (CFP1–3) Are Required for Insect-Pathogenic Lifestyle of Metarhizium robertsii

Unique CFP (cysteine-free protein; 120 aa) has been identified as an extraordinary virulence factor in Beauveria bassiana (Cordycipitaceae), a main source of wide-spectrum fungal insecticides. Its homologs exclusively exist in wide-spectrum insect pathogens of Hypocreales, suggesting their importance for a fungal insect-pathogenic lifestyle. In this study, all three CFP homologs (CFP1–3, 128–145 aa) were proven essential virulence factors in Metarhizium robertsii (Clavicipitaceae). Despite limited effects on asexual cycles in vitro, knockout mutants of cfp1,cfp2 and cfp3 were severely compromised in their capability for normal cuticle infection, in which most tested Galleria mellonella larvae survived. The blocked cuticle infection concurred with reduced secretion of extracellular enzymes, including Pr1 proteases required cuticle penetration. Cuticle-bypassing infection by intrahemocoel injection of ~250 conidia per larva resulted in a greater reduction in virulence in the mutant of cfp1 (82%) than of cfp2 (21%) or cfp3 (25%) versus the parental wild-type. Transcriptomic analysis revealed dysregulation of 604 genes (up/down ratio: 251:353) in the Δcfp1 mutant. Many of them were involved in virulence-related cellular processes and events aside from 154 functionally unknown genes (up/down ratio: 56:98). These results reinforce the essential roles of small CFP homologs in hypocrealean fungal adaptation to insect-pathogenic lifestyle and their exploitability for the genetic improvement of fungal insecticidal activity.

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

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

Virulence of the Bio-Control Fungus Purpureocillium lilacinum Against Myzus persicae (Hemiptera: Aphididae) and Spodoptera frugiperda (Lepidoptera: Noctuidae)

Eco-friendly entomopathogenic fungi are widely used to control agricultural insect pests. Purpureocillium lilacinum (Thom.) Luangsa-ard et al. (Hypocreales: Ophiocordycipitaceae) is a nematophagous fungus used for the bio-control of destructive root-knot nematodes. However, its insecticidal activities against agricultural insect pests haven't been widely studied. In this study, P. lilacinum PL-1 was isolated from soil (Hefei, China) and identified by molecular and morphological analyses. The growth rate, spore production, proteinase, and chitinase activities of the isolate were analyzed. Virulence tests against green peach aphid, Myzus persicae (Sulzer) (Hemiptera: Aphididae) and fall armyworm (FAW), Spodoptera frugiperda (Smith) (Lepidoptera: Noctuidae) were performed. The median lethal concentration (LC50) and median lethal time (LT50) against aphids (via immersion) and LT50 against FAW (via injection) were determined. FAW eggs immersed in aqueous conidia suspension were infected after 60 h. Differentially expressed genes (DEGs) in the infection of FAW larvae by P. lilacinum were analyzed by quantitative reverse transcription PCR. The significantly upregulated DEGs include FAW immune genes (antimicrobial peptides, C-type lectins, lysozymes, prophenoloxidase, and peptidoglycan recognition proteins) and fungal pathogenic genes (ligase, chitinase, and hydrophobin). Our data demonstrate that P. lilacinum can be used as an entomopathogenic fungus against agricultural insect pests.

Comprehensive insights into host-pathogen interaction between brown planthopper and a fungal entomopathogen by dual RNA sequencing

BACKGROUND

The brown planthopper (BPH) is one of the most destructive pests of rice, causing tremendous yield and economic losses every year. The fungal entomopathogen Metarhizium anisopliae was previously proved to have great potential for BPH biocontrol. Genome-wide insight into the insect-fungus interaction is crucial for genetic improvement of M. anisopliae to enhance its virulence to BPH but still has been poorly explored.

RESULTS

Using dual RNA-seq approach, we present here a global view of host and fungal gene expressions in BPH adults during the fungal infection. The results revealed that BPH could initiate strong defense responses against the fungal attack by upregulating the expressions of a large number of genes, including genes involved in cuticle formation, immune response, cell detoxification and biomacromolecule metabolism. Correspondingly, the fungal entomopathogen could induce a series of genes to infect and modulate BPH, including genes involved in fungal penetration, invasive growth, stress resistance and virulence. Three host defense-related genes (NlPCE4, NlPOD1 and NlCYP4DE1) were chosen for further function analysis. RNAi-mediated knockdown of NlPCE4 caused a significant decrease in BPH survival, but no obvious effects on the survival rates were detected by the suppression of NlPOD1 and NlCYP4DE1. Combination of dsRNA injection and fungal infection could significantly enhance the BPH-killing speed, as synergistic mortalities were observed in co-treatments of RNAi and M. anisopliae infection.

CONCLUSION

Our study provides a comprehensive insight into molecular mechanisms of host-pathogen interaction between BPH and M. anisopliae and contributes to future development of new efficient biocontrol strategies for BPH biocontrol.

Host-pathogen interaction between Asian citrus psyllid and entomopathogenic fungus (Cordyceps fumosorosea) is regulated by modulations in gene expression, enzymatic activity and HLB-bacterial population of the host

The host-pathogen interaction has been explored by several investigations, but the impact of fungal pathogens against insect resistance is still ambiguous. Therefore, we assessed the enzymatic activity and defense-related gene expression of Asian citrus psyllid (ACP) nymphal and adult populations on Huanglongbing-diseased citrus plants under the attack of Cordyceps fumosorosea. Overall, five enzymes viz. superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), glutathione S-transferase (GST), carboxylesterase (CarE), and four genes, namely SOD, 16S, CYP4C68, CYP4BD1, were selected for respective observations from ACP populations. Enzymatic activity of four enzymes (SOD, POD, GST, CarE) was significantly decreased after 5-days post-treatment (dpt) and 3-dpt fungal exposure in fungal treated ACP adult and nymphal populations, respectively, whereas the activity of CAT was boosted substantially post-treatment time schedule. Besides, we recorded drastic fluctuations in the expression of CYP4 genes among fungal treated ACP populations. After 24 hours post-treatment (hpt), expression of both CYP4 genes was boosted in fungal treated populations than controlled populations (adult and nymph). After 3-dpt, however, the expression of CYP4 genes was declined in the given populations. Likewise, fungal attack deteriorated the resistance of adult and nymphal of ACP population, as SOD expression was down-regulated in fungal-treated adult and nymphs after 5-dpt and 3-dpt exposure, respectively. Moreover, bacterial expression via the 16S gene was significantly increased in fungal-treated adult and nymphal ACP populations with increasing post-treatment time. Overall, our data illustrate that the fungal application disrupted the insect defense system. The expression of these genes and enzymes suppress the immune function of adult and nymphal ACP populations. As it is reported first time that the applications of C. fumosorosea against ACP reduce insect resistance by interfering with the CYP4 and SOD system. Therefore, we propose new strategies to discover the role of certain toxic compounds from fungus, which can reduce insect resistance, focusing on resistance-related genes and enzymes.

The Toxins of Beauveria bassiana and the Strategies to Improve Their Virulence to Insects

The long-term and excessive usage of pesticides is an enormous burden on the environment, which also increases pest resistance. To overcome this problem, research and application of entomopathogenic fungi, which are both environmentally friendly and cause lower resistance, have gained great momentum. Entomopathogenic fungi have a wide range of prospects. Apart from Bacillus thuringiensis, Beauveria bassiana is the most studied biopesticide. After invading insect hosts, B. bassiana produces a variety of toxins, which are secondary metabolites such as beauvericin, bassianin, bassianolide, beauverolides, tenellin, oosporein, and oxalic acid. These toxins help B. bassiana to parasitize and kill the hosts. This review unequivocally considers beauveria toxins highly promising and summarizes their attack mechanism(s) on the host insect immune system. Genetic engineering strategies to improve toxin principles, genes, or virulent molecules of B. bassiana have also been discussed. Lastly, we discuss the future perspective of Beauveria toxin research, including newly discovered toxins.

A Small Cysteine-Free Protein Acts as a Novel Regulator of Fungal Insect-Pathogenic Lifecycle and Genomic Expression

Small secreted proteins (SSPs), particularly cysteine-rich proteins secreted during fungal infection, comprise virulence effectors in plant-pathogenic fungi but remain unknown in insect-pathogenic fungi. We report here that only a small cysteine-free protein (CFP) is indispensable for insect pathogenicity of Beauveria bassiana among 10 studied SSPs (99 to 274 amino acids [aa]), including seven hypothetical proteins containing 0 to 12 Cys residues. CFP (120 aa) features an N-terminal signal peptide (residues 1 to 17), a nuclear localization signal motif (residues 24 to 57), and no predictable domain. Its homologs exist exclusively in insect-pathogenic Cordycipitaceae and Clavicipitaceae. Fluorescence-tagged CFP fusion protein was localized in the nucleus but extracellularly undetectable, suggesting an inability for CFP to be secreted out. Disruption of cfp resulted in abolished pathogenicity via normal cuticle infection, attenuated virulence via hemocoel injection, compromised conidiation capacity versus little growth defect, impaired conidial coat, blocked secretion of cuticle-degrading enzymes, impeded proliferation in vivo, disturbed cell cycle, reduced stress tolerance, and 1,818 dysregulated genes (genomic 17.54%). Hundreds of those genes correlated with phenotypic changes observed in the disruption mutant. Intriguingly, nearly 40% of those dysregulated genes encode hypothetical or unknown proteins, and another 13% encode transcription factors and enzymes or proteins collectively involved in genome-wide gene regulation. However, purified CFP showed no DNA-binding activity in an electrophoretic mobility shift assay. These findings unveil that CFP is a novel regulator of fungal insect-pathogenic life cycle and genomic expression and that cysteine richness is dispensable for distinguishing virulence effectors from putative SSPs in B. bassiana IMPORTANCE Small cysteine-rich proteins secreted during plant-pathogenic fungal infection comprise virulence effectors. Our study confirms that only a cysteine-free protein (CFP) is determinant to insect-pathogenic fungal virulence among 10 small putatively secreted proteins containing 0 to 12 Cys residues. Disruption of cfp abolished insect pathogenicity and caused not only a series of compromised cellular events associated with host infection and disease development but also dysregulation of 1,818 genes, although no DNA-binding activity was detected in purified CFP samples. Nearly 13% of those genes encode transcription factors and enzymes or proteins collectively involved in transcriptional regulation. Altogether, CFP serves as a novel regulator of the fungal insect-pathogenic life cycle and genomic expression. Cysteine richness is dispensable for distinguishing virulence effectors from the fungal SSPs.

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.

Immunotranscriptome analysis of Plutella xylostella reveals differences in innate immune responses to low- and high-virulence Beauveria bassiana strain challenges

BACKGROUND

Entomopathogenic fungi have developed multiple strategies to overcome the immune defenses of their target insects, whereas insect pests have devised various defense mechanisms to combat fungal infection. However, differences in the molecular mechanisms of the innate immune defense strategies of insects upon infection with different fungal strains from the same species have not been reported.

RESULTS

Two Beauveria bassiana strains were obtained that significantly varied in their pathogenicity but were comparable in terms of growth, conidial yield, and cuticle penetration. To investigate the molecular mechanisms underlying the immune response of Plutella xylostella infected with these two strains, RNA-Seq was performed 48 h after infection. A total of 1027 differentially expressed genes (DEGs) were identified, and more than 200 DEGs were enriched in Kyoto Encyclopedia of Genes and Genome (KEGG) pathways involved in disease response, revealing differences in the immune response of P. xylostella to different B. bassiana infections at 48 h. Twenty-eight of the DEGs were related to innate immune functions, such as pathogen recognition, immune system activation and antimicrobial reactions. RNA interference (RNAi)-mediated gene silencing assays showed that PxApoLIII and PxCSP played critical roles in the P. xylostella immune response. PxApoLIII was expressed at higher levels during infection with the high-virulence strain, whereas PxCSP showed the opposite expression pattern during infection with the low-virulence strain, indicating that PxApoLIII and PxCSP might participate in P. xylostella innate immune defense against high- and low-virulence B. bassiana strains.

CONCLUSION

The present findings demonstrate that strains of a single species of pathogenic fungi that differ in virulence can induce the expression of different genes in P. xylostella. These results advance our knowledge of the molecular mechanisms underlying fungi-pest interactions.

Nuclear Ssr4 Is Required for the In Vitro and In Vivo Asexual Cycles and Global Gene Activity of Beauveria bassiana

Ssr4 is known to serve as a cosubunit of chromatin-remodeling SWI/SNF and RSC complexes in yeasts but has not been functionally characterized in fungi. This study unveils for the first time the pleiotropic effects caused by deletion of ssr4 and its role in mediating global gene expression in a fungal insect pathogen. Our findings confirm an essential role of Ssr4 in hydrophobin biosynthesis and assembly required for growth, differentiation, and development of aerial hyphae for conidiation and conidial adhesion to insect surface and its essentiality for insect pathogenicity and virulence-related cellular events. Importantly, Ssr4 can regulate nearly one-fourth of all genes in the fungal genome in direct and indirect manners, including dozens involved in gene activity and hundreds involved in metabolism and/or transport of carbohydrates, amino acids, lipids, and/or inorganic ions. These findings highlight a significance of Ssr4 for filamentous fungal lifestyle.

Ssr4 serves as a cosubunit of chromatin-remodeling SWI/SNF and RSC complexes in yeasts but remains functionally uncharacterized due to its essentiality for yeast viability. Here, we report pleiotropic effects of the deletion of the ssr4 ortholog nonessential for cell viability in Beauveria bassiana, an asexual insect mycopathogen. The deletion of ssr4 resulted in severe growth defects on different carbon/nitrogen sources, increased hyphal hydrophilicity, blocked hyphal differentiation, and 98% reduced conidiation capacity compared to a wild-type standard. The limited Δssr4 conidia featured an impaired coat with disordered or obscure hydrophobin rodlet bundles, decreased hydrophobicity, increased size, and lost insect pathogenicity via normal cuticle infection and 90% of virulence via intrahemocoel injection. The expression of genes required for hydrophobin biosynthesis and assembly of the rodlet layer was drastically repressed in more hydrophilic Δssr4 cells. Transcriptomic analysis revealed 2,517 genes differentially expressed in the Δssr4 mutant, including 1,505 downregulated genes and 1,012 upregulated genes. The proteins encoded by hundreds of repressed genes were involved in metabolism and/or transport of carbohydrates, amino acids, and lipids, inorganic ion transport and energy production or conversion, including dozens involved in DNA replication, transcription, translation, and posttranslational modifications. However, purified Ssr4 samples showed no DNA-binding activity, implying that the role of Ssr4 in genome-wide gene regulation could rely upon its acting as a cosubunit of the two complexes. These findings provide the first insight into an essential role of Ssr4 in the asexual cycle in vitro and in vivo of B. bassiana and highlights its importance for the filamentous fungal lifestyle.

IMPORTANCE Ssr4 is known to serve as a cosubunit of chromatin-remodeling SWI/SNF and RSC complexes in yeasts but has not been functionally characterized in fungi. This study unveils for the first time the pleiotropic effects caused by deletion of ssr4 and its role in mediating global gene expression in a fungal insect pathogen. Our findings confirm an essential role of Ssr4 in hydrophobin biosynthesis and assembly required for growth, differentiation, and development of aerial hyphae for conidiation and conidial adhesion to insect surface and its essentiality for insect pathogenicity and virulence-related cellular events. Importantly, Ssr4 can regulate nearly one-fourth of all genes in the fungal genome in direct and indirect manners, including dozens involved in gene activity and hundreds involved in metabolism and/or transport of carbohydrates, amino acids, lipids, and/or inorganic ions. These findings highlight a significance of Ssr4 for filamentous fungal lifestyle.

Spatial and temporal transcriptomic analyses reveal locust initiation of immune responses to Metarhizium acridum at the pre-penetration stage

Insect hemocyte and fat body tissues play critical functional roles in insect immunity. Little, however, is known concerning the dynamic responses of these tissues to fungal infection. Here, we report on a time course of locust hemocyte and fat body transcriptomic responses to infection by the acridid specific fungal pathogen, Metarhizium acridum. Fat body responses were more pronounced at all infection stages as compared to hemocytes. Immune and other related genes were induced far earlier than previously considered including at pre-penetration stages. Differential expression in hemocyte and fat body tissues persisted throughout the course of infection up until host death. Our data indicate selective pressure on the host to recognize the infection as early as possible in order to limit its spread. Overall, fat body and hemocyte tissues launch a robust multi-tiered response to combat the fungal pathogen, with our data providing potential host targets for exploitation in pest control.

In vivo transcriptomic analysis of Beauveria bassiana reveals differences in infection strategies in Galleria mellonella and Plutella xylostella

BACKGROUND

Insect pests have evolved various defense mechanisms to combat fungal infection, and fungi have developed multiple strategies to overcome the immune defense responses of insects. However, transcriptomic analysis of fungal strategies for infecting different pests has not been reported.

RESULTS

Transcriptomic profiling of Beauveria bassiana was performed at 12, 24 and 48 h after infecting Galleria mellonella and Plutella xylostella, and 540, 847 and 932 differentially expressed genes were detected, respectively. Functional categorization showed that most of these genes are involved in the ribosome, nitrogen metabolism and oxidative phosphorylation pathways. Thirty-one differentially expressed virulence genes (including genes involved in adhesion, degradation, host colonization and killing, and secondary metabolism) were found, suggesting that different molecular mechanisms were used by the fungus during the infection of different pests, which was further confirmed by disrupting creA and fkh2. Virulence assay results showed that ΔcreA and Δfkh2 strains of B. bassiana had distinct fold changes in their 50% lethal time (LT₅₀ ) values (compared with the control stains) during infection of G. mellonella (ΔcreA: 1.38-fold > Δfkh2: 1.18-fold) and P. xylostella (ΔcreA: 1.44-fold < Δfkh2: 2.25-fold). creA was expressed at higher levels during the infection of G. mellonella compared with P. xylostella, whereas fkh2 showed the opposite expression pattern, demonstrating that creA and Fkh2 have different roles in B. bassiana during the infection of G. mellonella and P. xylostella.

CONCLUSION

These findings demonstrate that B. bassiana regulates different genes to infect different insects, advancing knowledge of the molecular mechanisms of Beauveria-pest interactions. © 2018 Society of Chemical Industry.

Rei1-like protein regulates nutritional metabolism and transport required for the asexual cycle in vitro and in vivo of a fungal insect pathogen

Rei1 is a cytoplasm-specific pre-60S subunit export factor that functions exclusively in cold-sensitive yeast growth but remains unexplored in filamentous fungi. Here, we report that Rei1-like BbRei1 is localized in both cytoplasm and nucleus and acts as a vital regulator in Beauveria bassiana. Deletion of BbRei1 resulted in delayed conidial germination, abnormally polarized germlings, severe growth defects on various carbon/nitrogen sources and reduced conidiation capacity as well as low temperature-sensitive growth. In ΔBbrei1, greatly attenuated virulence correlated with reduced activities of enzymes secreted for cuticular penetration and blocked formation of hyphal bodies in vivo essential for facilitation of host mummification. Revealed by transcriptomic analysis, 560 and 840 genes were significantly up- and down-regulated in ΔBbrei1 versus wild-type respectively, representing 13.5% of the fungal genome. Many repressed genes were involved in metabolism and transport of carbohydrates and amino acids. However, electrophoretic mobility shift assays presented no interactions of purified BbRei1 with 14 promoter DNA fragments. Conclusively, BbRei1 plays a pivotal role in gene expression and metabolism of nutrients and energy essential for the asexual cycle in vitro and in vivo of B. bassiana and functions much beyond the role for the yeast Rei1 in cold-sensitive cell growth.

Genes involved in Beauveria bassiana infection to Galleria mellonella

The ascomycete fungus Beauveria bassiana is a natural pathogen of hundreds of insect species and is commercially produced as an environmentally friendly mycoinsecticide. Many genes involved in fungal insecticide infection have been identified but few have been further explored. In this study, we constructed three transcriptomes of B. bassiana at 24, 48 and 72 h post infection of insect pests (BbI) or control (BbC). There were 3148, 3613 and 4922 genes differentially expressed at 24, 48 and 72 h post BbI/BbC infection, respectively. A large number of genes and pathways involved in infection were identified. To further analyze those genes, expression patterns across different infection stages (0, 12, 24, 36, 48, 60, 72 and 84 h) were studied using quantitative RT-PCR. This analysis showed that the infection-related genes could be divided into four patterns: highly expressed throughout the whole infection process (thioredoxin 1); highly expressed during early stages of infection but lowly expressed after the insect death (adhesin protein Mad1); lowly expressed during early infection but highly expressed after insect death (cation transporter, OpS13); or lowly expressed across the entire infection process (catalase protein). The data provide novel insights into the insect-pathogen interaction and help to uncover the molecular mechanisms involved in fungal infection of insect pests.

The Entomopathogenic Fungi Isaria fumosorosea Plays a Vital Role in Suppressing the Immune System of Plutella xylostella: RNA-Seq and DGE Analysis of Immunity-Related Genes

Most, if not all, entomopathogenic fungi have been used as alternative control agents to decrease the insect resistance and harmful effects of the insecticides on the environment. Among them, Isaria fumosorosea has also shown great potential to control different insect pests. In the present study, we explored the immune response of P. xylostella to the infection of I. fumosorosea at different time points by using RNA-Sequencing and differential gene expression technology at the genomic level. To gain insight into the host-pathogen interaction at the genomic level, five libraries of P. xylostella larvae at 12, 18, 24, and 36 h post-infection and a control were constructed. In total, 161 immunity-related genes were identified and grouped into four categories; immune recognition families, toll and Imd pathway, melanization, and antimicrobial peptides (AMPs). The results of differentially expressed immunity-related genes depicted that 15, 13, 53, and 14 up-regulated and 38, 51, 56, and 49 were down-regulated in P. xylostella at 12, 18, 24, and 36 h post-treatment, respectively. RNA-Seq results of immunity-related genes revealed that the expression of AMPs was reduced after treatment with I. fumosorosea. To validate RNA-Seq results by RT-qPCR, 22 immunity-related genes were randomly selected. In conclusion, our results demonstrate that I. fumosorosea has the potential to suppress the immune response of P. xylostella and can become a potential biopesticide for controlling P. xylostella.

Vital role for cyclophilin B (CypB) in asexual development, dimorphic transition and virulence of Beauveria bassiana

Cyclophilin B (CypB) was previously revealed as one of many putative secretory proteins in the transcriptome of Beauveria bassiana infection to a lepidopteran pest. Here we show a main localization of CypB in hyphal cell walls and septa and its essential role in the in vitro and in vivo asexual cycles of the fungal insect pathogen. Deletion of cypB reduced colony growth by 16-42% on two rich media and 30 scant media with different carbon or nitrogen sources. The deletion mutant suffered a delayed conidiation on a standard medium and a final 47% reduction in conidial yield, accompanied with drastic transcript depression of several key genes required for conidiation and conidial maturation. The mutant conidia required 10h longer to germinate 50% at optimal 25°C than wild-type conidia. Intriguingly, cultivation of the mutant conidia in a trehalose-peptone broth mimic to insect hemolymph resulted in 83% reduction in blastospore yield but only slight decrease in biomass level, indicating severe defects in transition of hyphae to blastospores. LT₅₀ for the deletion mutant against Galleria mellonella larvae through normal cuticle infection was prolonged to 7.4d from a wild-type estimate of 4.7d. During colony growth, additionally, the deletion mutant displayed hypersensitivity to Congo red, menadione, H₂O₂ and heat shock but increased tolerance to cyclosporine A and rapamycin. All of changes were restored by targeted gene complementation. Altogether, CypB takes part in sustaining normal growth, aerial conidiation, conidial germination, dimorphic transition, stress tolerance and pathogenicity in B. bassiana.

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

BACKGROUND

Entomopathogenic fungi such as Beauveria bassiana are promising biological agents for control of malaria mosquitoes. Indeed, infection with B. bassiana reduces the lifespan of mosquitoes in the laboratory and in the field. Natural isolates of B. bassiana show up to 10-fold differences in virulence between the most and the least virulent isolate. In this study, we sequenced the genomes of five isolates representing the extremes of low/high virulence and three RNA libraries, and applied a genome comparison approach to uncover genetic mechanisms underpinning virulence.

RESULTS

A high-quality, near-complete genome assembly was achieved for the highly virulent isolate Bb8028, which was compared to the assemblies of the four other isolates. Whole genome analysis showed a high level of genetic diversity between the five isolates (2.85-16.8 SNPs/kb), which grouped into two distinct phylogenetic clusters. Mating type gene analysis revealed the presence of either the MAT1-1-1 or the MAT1-2-1 gene. Moreover, a putative new MAT gene (MAT1-2-8) was detected in the MAT1-2 locus. Comparative genome analysis revealed that Bb8028 contains 163 genes exclusive for this isolate. These unique genes have a tendency to cluster in the genome and to be often located near the telomeres. Among the genes unique to Bb8028 are a Non-Ribosomal Peptide Synthetase (NRPS) secondary metabolite gene cluster, a polyketide synthase (PKS) gene, and five genes with homology to bacterial toxins. A survey of candidate virulence genes for B. bassiana is presented.

CONCLUSIONS

Our results indicate several genes and molecular processes that may underpin virulence towards mosquitoes. Thus, the genome sequences of five isolates of B. bassiana provide a better understanding of the natural variation in virulence and will offer a major resource for future research on this important biological control agent.

Perplexing Metabolomes in Fungal-Insect Trophic Interactions: A Terra Incognita of Mycobiocontrol Mechanisms

The trophic interactions of entomopathogenic fungi in different ecological niches viz., soil, plants, or insect themselves are effectively regulated by their maneuvered metabolomes and the plethora of metabotypes. In this article, we discuss a holistic framework of co-evolutionary metabolomes and metabotypes to model the interactions of biocontrol fungi especially with mycosed insects. Conventionally, the studies involving fungal biocontrol mechanisms are reported in the context of much aggrandized fungal entomotoxins while the adaptive response mechanisms of host insects are relatively overlooked. The present review asserts that the selective pressure exerted among the competing or interacting species drives alterations in their overall metabolomes which ultimately implicates in corresponding metabotypes. Quintessentially, metabolomics offers a most generic and tractable model to assess the fungal-insect antagonism in terms of interaction biomarkers, biosynthetic pathway plasticity, and their co-evolutionary defense. The fungi chiefly rely on a battery of entomotoxins viz., secondary metabolites falling in the categories of NRP's (non-ribosomal peptides), PK's (polyketides), lysine derive alkaloids, and terpenoids. On the contrary, insects overcome mycosis through employing different layers of immunity manifested as altered metabotypes (phenoloxidase activity) and overall metabolomes viz., carbohydrates, lipids, fatty acids, amino acids, and eicosanoids. Here, we discuss the recent findings within conventional premise of fungal entomotoxicity and the evolution of truculent immune response among host insect. The metabolomic frameworks for fungal-insect interaction can potentially transmogrify our current comprehensions of biocontrol mechanisms to develop the hypervirulent biocontrol strains with least environmental concerns. Moreover, the interaction metabolomics (interactome) in complementation with other -omics cascades could further be applied to address the fundamental bottlenecks of adaptive co-evolution among biological species.