RNA-Seq Analyses for Two Silkworm Strains Reveals Insight into Their Susceptibility and Resistance to Beauveria bassiana Infection

Antifungal and antifeedant terpenoids from Paraphaeosphaeria sp. cultured by extract of host Ginkgo biloba

Five undescribed terpenoids including a polyketide-terpenoid hybrid paraphaone, and four eremophilane sesquiterpenoids, paraphaterpenes A-D, as well as two known compounds were isolated from the endophytic fungus Paraphaeosphaeria sp. cultured by extract of host Ginkgo biloba L. The structures were established by spectroscopic analyses, and the single-crystal X-ray diffraction. The antifungal activity of Paraphaeosphaeria sp. cultured by extract of G. biloba against the plant pathogen Alternaria alternata was significant and higher than that of PDB medium. Tested compounds indicated antifeedant activities against silkworms with feeding deterrence index at 10-70%. Paraphaone, paraphaterpenes A, C, D and alternariol methyl ether showed significant antifungal activities against the phytopathogens A. alternata, Aspergillus fumigatus, and entomopathogen Beauveria bassiana with MICs ≤4 μg/mL. And the preliminary structure-activity relationship of eremophilane sesquiterpenoids was exhibited. The culture of Paraphaeosphaeria sp. by host G. biloba medium afforded agricultural antibiotics.

Differential bumble bee gene expression associated with pathogen infection and pollen diet

BACKGROUND

Diet and parasitism can have powerful effects on host gene expression. However, how specific dietary components affect host gene expression that could feed back to affect parasitism is relatively unexplored in many wild species. Recently, it was discovered that consumption of sunflower (Helianthus annuus) pollen reduced severity of gut protozoan pathogen Crithidia bombi infection in Bombus impatiens bumble bees. Despite the dramatic and consistent medicinal effect of sunflower pollen, very little is known about the mechanism(s) underlying this effect. However, sunflower pollen extract increases rather than suppresses C. bombi growth in vitro, suggesting that sunflower pollen reduces C. bombi infection indirectly via changes in the host. Here, we analyzed whole transcriptomes of B. impatiens workers to characterize the physiological response to sunflower pollen consumption and C. bombi infection to isolate the mechanisms underlying the medicinal effect. B. impatiens workers were inoculated with either C. bombi cells (infected) or a sham control (un-infected) and fed either sunflower or wildflower pollen ad libitum. Whole abdominal gene expression profiles were then sequenced with Illumina NextSeq 500 technology.

RESULTS

Among infected bees, sunflower pollen upregulated immune transcripts, including the anti-microbial peptide hymenoptaecin, Toll receptors and serine proteases. In both infected and un-infected bees, sunflower pollen upregulated putative detoxification transcripts and transcripts associated with the repair and maintenance of gut epithelial cells. Among wildflower-fed bees, infected bees downregulated immune transcripts associated with phagocytosis and the phenoloxidase cascade.

CONCLUSIONS

Taken together, these results indicate dissimilar immune responses between sunflower- and wildflower-fed bumble bees infected with C. bombi, a response to physical damage to gut epithelial cells caused by sunflower pollen, and a strong detoxification response to sunflower pollen consumption. Identifying host responses that drive the medicinal effect of sunflower pollen in infected bumble bees may broaden our understanding of plant-pollinator interactions and provide opportunities for effective management of bee pathogens.

Changes of the gene expression in silkworm larvae and Cordyceps militaris at late stages of the pathogenesis

Cordyceps militaris is an entomopathogenic fungus that forms its fruiting body. The gene expression change in C. militaris and silkworm larvae were analyzed using RNA-seq to investigate the relationship of C. militaris with the host, silkworm larvae before the death by mycosis. At 144 h after the injection of C. militaris conidia, genes encoding proteases, protease inhibitors, and cuticle proteins in the fat body of silkworm larvae were upregulated, but genes encoding lipoproteins and other proteins in hemolymph were downregulated. On the other hand, at 168 h after the injection of C. militaris conidia, genes encoding amino acid and oligopeptide transporters and permeases in C. militaris were upregulated, suggesting that C. militaris may use peptides and amino acids in silkworm larvae as a nutrient to grow in vivo. Additionally, one gene cluster composed of genes putatively involved in the degradation of phenolic substrates was also upregulated. The addition of 4,5-dichlorocatechol, an inhibitor of catechol 1,2-dioxygenase, inhibited the in vivo growth of C. militaris, Beauveria bassiana and Metarhizium anisopliae. These results also suggest that the expression of the gene cluster may be crucial for the in vivo growth of C. militaris and entomopathogenic fungi. This study will clarify how C. militaris grows in insect hosts by avoiding host's immune systems.

Gene Expression Analysis of Adapted Insect Cells during Influenza VLP Production Using RNA-Sequencing

Adaptive laboratory evolution has been used to improve production of influenza hemagglutinin (HA)-displaying virus-like particles (VLPs) in insect cells. However, little is known about the underlying biological mechanisms promoting higher HA-VLP expression in such adapted cell lines. In this article, we present a study of gene expression patterns associated with high-producer insect High Five cells adapted to neutral pH, in comparison to non-adapted cells, during expression of influenza HA-VLPs. RNA-seq shows a decrease in the amount of reads mapping to host cell genomes along infection, and an increase in those mapping to baculovirus and transgenes. A total of 1742 host cell genes were found differentially expressed between adapted and non-adapted cells throughout infection, 474 of those being either up- or down-regulated at both time points evaluated (12 and 24 h post-infection). Interestingly, while host cell genes were found up- and down-regulated in an approximately 1:1 ratio, all differentially expressed baculovirus genes were found to be down-regulated in infected adapted cells. Pathway analysis of differentially expressed genes revealed enrichment of ribosome biosynthesis and carbohydrate, amino acid, and lipid metabolism. In addition, oxidative phosphorylation and protein folding, sorting and degradation pathways were also found to be overrepresented. These findings contribute to our knowledge of biological mechanisms of insect cells during baculovirus-mediated transient expression and will assist the identification of potential engineering targets to increase recombinant protein production in the future.

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.

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.

Beauveria bassiana ERL836 and JEF-007 with similar virulence show different gene expression when interacting with cuticles of western flower thrips, Frankniella occidentalis

Background

Insect-killing fungal species, Beauveria bassiana, is as an environment-friendly pest management tool, and many isolates are on the track of industrialization. However, some of B. bassiana isolates show similar morphology and virulence against insect pests, and so it is hard to differentiate them. Herein we used two patented isolates, ERL836 and JEF-007, and investigated their virulence against western flower thrips, Frankliniella occidentalis, and further analyzed genome structures and transcriptional responses when interacting with cuticles of thrips to see possible differences on the initial step of fungal infection.

Results

The two isolates showed no significant differences in fungal growth, conidial production, and virulence against thrips, and they were structurally similar in genome. But, in transcription level, ERL836 appeared to infect thrips easily, while JEF-007 appeared to have more difficulty. In the GO analysis of ERL836 DEGs (differentially expressed genes), the number of up-regulated genes was much larger than that of down-regulated genes, when compared to JEF-007 DEGs (more genes down-regulated). Interestingly, in the enrichment analysis using shared DEGs between two infecting isolates, plasma membrane-mediated transporter activity and fatty acid degradation pathway including cytochrome P450 were more active in infecting ERL836.

Conclusion

The two B. bassiana isolates had similar morphology and virulence as well as genome structure, but in transcription level they differently interacted with the cuticle of western flower thrips. This comparative approach using shared DEG analysis could be easily applied to characterize the difference of the two B. bassiana isolates, JEF-007 and ERL836.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07253-y.

Beauveria bassiana ribotoxin inhibits insect immunity responses to facilitate infection via host translational blockage

Entomopathogenic fungi are promising bio-pesticides. To facilitate infection, fungi recruit multiple virulence factors and deploy different molecular strategies to evade host immunity. Fungal ribotoxins are extracellular secreted ribonucleases (RNases) with ribotoxic cytotoxicity and insecticidal activity. However, it remains unclear whether they have further biological functions. Here we show that the entomopathogenic fungus Beauveria bassiana ribotoxin (Rib) contributes to fungal virulence by inhibiting insect host immunity. Gene deletion of Rib (ΔRib) resulted in attenuated fungal virulence during infection. Pathogenesis analysis demonstrated that Rib mainly inhibits insect immunity through modulating the reactive oxygen species (ROS) response, suppressing antimicrobial peptides (AMPs) production and retarding hyphae penetration from insect cuticles. To further confirm this immunosuppressive function, recombinant ribotoxin (rRib) protein was purified and co-injected with living or heat-killed bacteria, bacteria-derived peptidoglycan (PGN) and lipopolysaccharide (LPS) separately, which also significantly inhibited the AMPs production in Drosophila fat bodies. Furthermore, co-injection of rRib with Escherichia coli or Staphylococcus aureus significantly enhanced bacterial pathogenicity and facilitated infection. In addition, rRib injection resulted in a global inhibition of protein expression in different tissues of Drosophila adults. This work identified B. bassiana ribotoxin as a key virulence factor that inhibits insect immunity.

Quality Formation Mechanism of Stiff Silkworm, Bombyx batryticatus Using UPLC-Q-TOF-MS-Based Metabolomics

Bombyx batryticatus is a well-known animal in traditional Chinese medicine. The aim of the research was to reveal the quality formation mechanism of B. batryticatus and to screen out the characteristic component used for the quality control. The anticonvulsant effects of B. batryticatus with a stiff time of one, five, and nine days (D1, D5 and D9, respectively) and healthy silkworm of the same developmental stage (SW) were determined by animal experiment. The dynamic changes in chemical composition were analyzed using UPLC-Q-TOF-MS-based metabolomics. D5 and D9 B. batryticatus exhibited significant anticonvulsant effects (p < 0.05 and p < 0.01, respectively). Accordingly, principal component analysis (PCA) and partial least squares discrimination analysis (PLS-DA) indicated that the chemical composition of D5 and D9 B. batryticatus changed significantly. The different metabolites mainly consisted of primary metabolites such as lipids and amino acids and secondary metabolites such as flavonoids, beauvericin, and glycolipids. Interestingly, the relative abundance of quercetin-7-O-β-d-4-O-methylglucoside, the characteristic component of B. batryticatus, increased with stiff time and was promised to be used as an index component of quality control. The results expand our understanding of the quality formation mechanism of B. batryticatus. In addition, it highlights the potential of UPLC-Q-TOF-MS-based metabolomics for the quality control purpose of TCMs.

Identification of cordycepin biosynthesis-related genes through de novo transcriptome assembly and analysis in Cordyceps cicadae

Cordyceps cicadae (Chanhua) is a parasitic fungus that grows on Cicada flammata larvae and is used to relieve exhaustion and treat numerous diseases, in part through its active constituent, cordycepin. We used de novo Illumina HiSeq 4000 sequencing to obtain transcriptomes of C. cicadae mycelium, fruiting body, and sclerotium, and identify differentially expressed genes. In the mycelium versus sclerotium libraries, 1576 upregulated and 2300 downregulated genes were identified. In the mycelium versus fruiting body and fruiting body versus sclerotium body libraries, 1604 and 1474 upregulated and 1365 and 1320 downregulated genes, respectively, were identified. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses identified 19 genes differentially expressed in mycelium versus fruiting body as related to the purine pathway, along with 28 and 16 genes differentially expressed in the mycelium versus sclerotium and fruiting body versus sclerotium groups, respectively. Gene expression of six key enzymes was validated by quantitative polymerase chain reaction. Specifically, 5'-nucleotidase (c62060g1) and adenosine deaminase (c35629g1) in purine nucleotide metabolism, which are involved in cordycepin biosynthesis, were significantly upregulated in the sclerotium group. These findings improved our understanding of genes involved in the biosynthesis of cordycepin and other characteristic secondary metabolites in C. cicadae.

Comparative transcriptome analysis reveals significant metabolic alterations in eri-silkworm (Samia cynthia ricini) haemolymph in response to 1-deoxynojirimycin

Samia cynthia ricini (Lepidoptera: Saturniidae) is an important commercial silk-producing insect; however, in contrast to the silkworm, mulberry leaves are toxic to this insect because the leaves contain the component 1-deoxynojirimycin (DNJ). A transcriptomic analysis of eri-silkworm haemolymph was conducted to examine the genes related to different metabolic pathways and to elucidate the molecular mechanism underlying eri-silkworm haemolymph responses to DNJ. Eight hundred sixty-five differentially expressed genes (DEGs) were identified, among which 577 DEGs were up-regulated and 288 DEGs were down-regulated in the 2% DNJ group compared to control (ddH₂O) after 12h. Based on the results of the functional analysis, these DEGs were associated with ribosomes, glycolysis, N-glycan biosynthesis, and oxidative phosphorylation. In particular, according to the KEGG analysis, 138 DEGs were involved in energy metabolism, glycometabolism and lipid metabolism, and the changes in the expression of nine DEGs were confirmed by reverse transcription quantitative PCR (RT-qPCR). Thus, DNJ induced significant metabolic alterations in eri-silkworm haemolymph. These results will lay the foundation for research into the toxic effects of DNJ on eri-silkworm as a model and provide a reference for the exploitation of new drugs in humans.