Comparative analysis of LC-MS/MS and real-time PCR assays for efficient detection of potential allergenic silkworm

The silkworm (Bombyx mori) has long been valued food and feed in East Asia for its abundant nutritional and medicinal attributes, conversely, it can elicit allergic responses in susceptible individuals. Therefore, the development of silkworm detection method is required to avert allergenic incidents. In this study, two methodologies, tandem mass spectrometry (LC-MS/MS) and real-time PCR, were developed to achieve effective silkworm detection. These methods exhibited exceptional sensitivity in identifying silkworm presence in processed foods. Furthermore, model cookies spiked with silkworm were used to validate the sensitivities of LC-MS/MS (0.0005%) and real-time PCR (0.001%). Overall, these techniques were useful for trace silkworm detection in food products; therefore, they may help prevent allergic reactions. To the best of our knowledge, this study represents the first comparison of LC-MS/MS and real-time PCR methods for silkworm detection, marking an important contribution to the field. Data are available from ProteomeXchange under identifier PXD042494.

Study on anti-BmNPV mechanism of branched-chain amino acid aminotransferases in silkworm

Bombyx mori nucleopolyhedrovirus (BmNPV) is the most important virus that threatens sericulture industry. At present, there is no effective treatment for BmNPV infection in silkworms, and lncRNA plays an important role in biological immune response and host-virus interaction, but there are relatively few studies in silkworms. In this study, the four midgut tissue samples of the resistance strain NB (NB) and susceptible strain 306 (306) and the NB and 306 continuously infected with BmNPV for 96 h are used for whole transcriptome sequencing to analyze the differences in the genetic background of NB and 306 and the differences after inoculation of BmNPV, and the significantly different mRNA, miRNA and lnRNA between NB and 306 after BmNPV inoculation were screened. By comparing NB and 306, 2651 significantly different mRNAs, 57 significantly different miRNAs and 198 significantly different lncRNAs were screened. By comparing NB and 306 after BmNPV inoculation, 2684 significantly different mRNAs, 39 significantly different miRNAs and 125 significantly different lncRNAs were screened. According to the significantly different mRNA, miRNA and lncRNA screened from NB and 306 and NB and 306 after virus inoculation, the mRNA-miRNA-lncRNA regulatory network was constructed before and after virus inoculation, and the BmBCAT-Bomo_chr7_8305-MSTRG.3236.2 regulatory axis was screened from them, and it was found that BmBCAT was not Bomo_chr7_8305 regulated in the genetic background, after viral infection, MSTRG.3236.2 competes for binding Bomo_chr7_8305 regulates BmBCAT. The whole transcriptome sequencing results were verified by qPCR and the time-series expression analysis was performed to prove the reliability of the regulatory network. The BmBCAT-Bomo_chr7_8305-MSTRG.3236.2 regulatory axis may play a potential role in the interaction between silkworms and BmNPV. These results provide new insights into the interaction mechanism between silkworms and BmNPV.

Physiological and transcriptomic responses of silkworms to graphene oxide exposure

The growing use of nanomaterials has sparked significant interest in assessing the insect toxicities of nanoparticles. The silkworm, as an economically important insect, serves as a promising model for studying how insects respond to harmful substances. Here, we conducted a comprehensive investigation on the impact of graphene oxide (GO) on silkworms using a combination of physiological and transcriptome analyses. GO can enter the midguts and posterior silk glands of silkworms. High GO concentrations (> 25 mg/L) significantly (P < 0.01) inhibited larval growth. Additionally, GO (> 5 mg/L) significantly reduced the cocooning rate, and GO (> 15 mg/L) hindered oviduct development and egg laying in silkworms. GO increased the reactive oxygen species content and regulated catalase activity, suggesting that it may affect insect growth by regulating reactive oxygen detoxification. The transcriptome data analysis showed that 35 metabolism-related genes and 20 ribosome biogenesis-related genes were differentially expressed in response to GO, and their expression levels were highly correlated. Finally, we propose that a Ribosome biogenesis-Metabolic signaling network is involved in responses to GO. The research provides a new perspective on the molecular responses of insects to GO.

Codon Optimization-based Whole-gene Scanning Identifies Hidden Nucleotides Essential for Bombyx mori Nucleopolyhedrovirus polyhedrin Hyperexpression

During the late stage of infection, alphabaculoviruses produce many occlusion bodies (OBs) in the nuclei of the insect host's cells through the hyperexpression of polyhedrin (POLH), a major OB component encoded by polh. The strong polh promoter has been used to develop a baculovirus expression vector system for recombinant protein expression in cultured insect cells and larvae. However, the relationship between POLH accumulation and the polh coding sequence remains largely unelucidated. This study aimed to assess the importance of polh codon usage and/or nucleotide sequences in POLH accumulation by generating a baculovirus Bombyx mori nucleopolyhedrovirus (BmNPV) expressing mutant polh (co-polh) optimized according to the codon preference of its host insect. Although the deduced amino acid sequence of CO-POLH was the same as that of wild-type POLH, POLH accumulation was significantly lower in cells infected with the co-polh mutant. This reduction was due to decreased polh mRNA levels rather than translational repression. Analysis of mutant viruses with chimeric polh revealed that a 30 base-pair (bp) 5' proximal polh coding region was necessary for maintaining high polh mRNA levels. Sequence comparison of wild-type polh and co-polh identified five nucleotide differences in this region, indicating that these nucleotides were critical for polh hyperexpression. Furthermore, luciferase reporter assays showed that the 30 bp 5' coding region was sufficient for maintaining the polh promoter-driven high level of polh mRNA. Thus, our whole-gene scanning by codon optimization identified important hidden nucleotides for polh hyperexpression in alphabaculoviruses.

Identification of long non-coding RNAs in response to microsporidia infection in Silkworm, Bombyx mori

Microsporidia Nosema bombycis (Nb) is a cellular parasite responsible for pébrine disease in silkworms, significantly impacting the sericulture industry. Long non-coding RNAs (lncRNAs), which are RNA fragments longer than 200 nucleotides, are pivotal in a range of cellular and physiological functions. However, the potential role of silkworm lncRNAs in response to Nb infection remains unknown. This study conducted transcriptome sequencing on both larvae and Nb-infected midguts of silkworms, identifying 1,440 lncRNAs across all examined midgut samples. Within the Nb-infected group, 42 differentially expressed lncRNAs (DElncRNAs) and 305 differentially expressed mRNAs (DEmRNAs) were detected. Functional annotation and pathway analysis showed that these DEmRNAs are mostly involved in metabolism, apoptosis, autophagy, and other key pathways. The co-expression network of DEmRNAs and DElncRNAs illustrates that 1 gene could be regulated by multiple lncRNAs and 1 lncRNA may target multiple genes, indicating that the regulation of lncRNA is intricate and networked. In addition, the DElncRNA-miRNA-mRNA network showed that some DElncRNAs may be involved in the immune response and metabolism through miRNA. Notably, the study observed an increase in lncRNA MSTRG857.1 following Nb infection, which may promote Nb proliferation. These findings offer insights into the complex interplay between insects and microsporidia.

Perilipin1 inhibits Nosema bombycis proliferation by promoting Domeless- and Hop-mediated JAK-STAT pathway activation in Bombyx mori

Lipid droplets (LDs) are dynamic organelles that participate in the regulation of lipid metabolism and cellular homeostasis inside of cells. LD-associated proteins, also known as perilipins (PLINs), are a family of proteins found on the surface of LDs that regulate lipid metabolism, immunity, and other functions. In silkworms, pébrine disease caused by infection by the microsporidian Nosema bombycis (Nb) is a severe threat to the sericultural industry. Although we found that Nb relies on lipids from silkworms to facilitate its proliferation, the relationship between PLINs and Nb proliferation remains unknown. Here, we found Nb infection caused the accumulation of LDs in the fat bodies of silkworm larvae. The characterized perilipin1 gene (plin1) promotes the accumulation of intracellular LDs and is involved in Nb proliferation. plin1 is similar to perilipin1 in humans and is conserved in all insects. The expression of plin1 was mostly enriched in the fat body rather than in other tissues. Knockdown of plin1 enhanced Nb proliferation, whereas overexpression of plin1 inhibited its proliferation. Furthermore, we confirmed that plin1 increased the expression of the Domeless and Hop in the JAK-STAT immune pathway and inhibited Nb proliferation. Taken together, our current findings demonstrate that plin1 inhibits Nb proliferation by promoting the JAK-STAT pathway through increased expression of Domeless and Hop. This study provides new insights into the complicated connections among microsporidia pathogens, LD surface proteins, and insect immunity.IMPORTANCELipid droplets (LDs) are lipid storage sites in cells and are present in almost all animals. Many studies have found that LDs may play a role in host resistance to pathogens and are closely related to innate immunity. The present study found that a surface protein of insect lipid droplets could not only regulate the morphological changes of lipid droplets but also inhibit the proliferation of a microsporidian pathogen Nosema bombycis (Nb) by activating the JAK-STAT signaling pathway. This is the first discovery of the relationship between microsporidian pathogen and insect lipid surface protein perilipin and insect immunity.

BmNPV p35 regulates apoptosis in Bombyx mori via a novel target of interaction with the BmVDAC2-BmRACK1 complex

Voltage-dependent anion channel 2 (VDAC2) is an important channel protein that plays a crucial role in the host response to viral infection. The receptor for activated C kinase 1 (RACK1) is also a key host factor involved in viral replication. Our previous research revealed that Bombyx mori VDAC2 (BmVDAC2) and B. mori RACK1 (BmRACK1) may interact with Bombyx mori nucleopolyhedrovirus (BmNPV), though the specific molecular mechanism remains unclear. In this study, the interaction between BmVDAC2 and BmRACK1 in the mitochondria was determined by various methods. We found that BmNPV p35 interacts directly with BmVDAC2 rather than BmRACK1. BmNPV infection significantly reduced the expression of BmVDAC2, and activated the mitochondrial apoptosis pathway. Overexpression of BmVDAC2 in BmN cells inhibited BmNPV-induced cytochrome c (cyto c) release, decrease in mitochondrial membrane potential as well as apoptosis. Additionally, the inhibition of cyto c release by BmVDAC2 requires the involvement of BmRACK1 and protein kinase C. Interestingly, overexpression of p35 inhibited cyto c release during mitochondrial apoptosis in a RACK1 and VDAC2-dependent manner. Even the mutant p35, which loses Caspase inhibitory activity, could still bind to VDAC2 and inhibit cyto c release. In summary, our results indicated that BmNPV p35 interacts with the VDAC2-RACK1 complex to regulate apoptosis by inhibiting cyto c release. These findings confirm the interaction between BmVDAC2 and BmRACK1, the interaction between p35 and the VDAC2-RACK1 complex, and a novel target that BmNPV p35 regulates apoptosis in Bombyx mori via interaction with the BmVDAC2-BmRACK1 complex. The result provide an initial exploration of the function of this interaction in the BmNPV-induced mitochondrial apoptosis pathway.

Comprehensive transcriptome sequencing of silkworm Midguts: Uncovering extensive isoform diversity and alternative splicing in BmNPV-Sensitive and BmNPV-resistant strains

The silkworm, Bombyx mori, stands out as one of the few economically valuable insects within the realm of model organisms. However, Bombyx mori nucleopolyhedrovirus (BmNPV) poses a significant threat, decreasing the quality and quantity of silkworm cocoons. Over the past few decades, a multitude of researchers has delved into the mechanisms that underlie silkworm resistance to BmNPV, employing diverse methodologies and approaching the problem from various angles. Despite this extensive research, the role of alternative splicing (AS) in the silkworm's response to BmNPV infection has been largely unexplored. This study leveraged both third-generation (Oxford Nanopore Technologies) and second-generation (Illumina) high-throughput sequencing technologies to meticulously identify and analyze AS patterns in the context of BmNPV response, utilizing two distinct silkworm strains-the susceptible strain 306 and the resistant strain NB. Consequently, we identified five crucial genes (Dsclp, LOC692903, LOC101743583, LOC101742498, LOC101743809) that are linked to the response to BmNPV infection through AS and differential expression. Additionally, a thorough comparative analysis was conducted on their diverse transcriptomic expression profiles, including alternative polyadenylation, simple sequence repeats, and transcription factors.

Studying the role of Bombyx mori molybdenum cofactor sulfurase in Bombyx mori nucleopolyhedrovirus infection

Molybdenum cofactor sulfurase (MoCoS) is a key gene involved in the uric acid metabolic pathway that activates xanthine dehydrogenase to synthesise uric acid. Uric acid is harmful to mammals but plays crucial roles in insects, one of which is the immune responses. However, the function of Bombyx mori MoCoS in response to BmNPV remains unclear. In this study, BmMoCoS was found to be relatively highly expressed in embryonic development, gonads and the Malpighian tubules. In addition, the expression levels of BmMoCoS were significantly upregulated in three silkworm strains with different levels of resistance after virus infection, suggesting a close link between them. Furthermore, RNAi and overexpression studies showed that BmMoCoS was involved in resistance to BmNPV infection, and its antivirus effects were found to be related to the regulation of uric acid metabolism, which was uncovered by inosine- and febuxostat-coupled RNAi and overexpression. Finally, the BmMoCoS-mediated uric acid pathway was preliminarily confirmed to be a potential target to protect silkworms from BmNPV infection. Overall, this study provides new evidence for elucidating the molecular mechanism of silkworms in response to BmNPV infection and new strategies for the prevention of viral infections in sericulture.

Bmo-miR-6498-5p suppresses Bombyx mori nucleopolyhedrovirus infection by down-regulating BmPLPP2 to modulate pyridoxal phosphate content in B. mori

The RNA interference pathway mediated by microRNAs (miRNAs) is one of the methods to defend against viruses in insects. Recent studies showed that miRNAs participate in viral infection by binding to target genes to regulate their expression. Here, we found that the Bombyx mori miRNA, miR-6498-5p was down-regulated, whereas its predicted target gene pyridoxal phosphate phosphatase PHOSPHO2 (BmPLPP2) was up-regulated upon Bombyx mori nucleopolyhedrovirus (BmNPV) infection. Both in vivo and in vitro experiments showed that miR-6498-5p targets BmPLPP2 and suppresses its expression. Furthermore, we found miR-6498-5p inhibits BmNPV genomic DNA (gDNA) replication, whereas BmPLPP2 promotes BmNPV gDNA replication. As a pyridoxal phosphate (PLP) phosphatase (PLPP), the overexpression of BmPLPP2 results in a reduction of PLP content, whereas the knockdown of BmPLPP2 leads to an increase in PLP content. In addition, exogenous PLP suppresses the replication of BmNPV gDNA; in contrast, the PLP inhibitor 4-deoxypyridoxine facilitates BmNPV gDNA replication. Taken together, we concluded that miR-6498-5p has a potential anti-BmNPV role by down-regulating BmPLPP2 to modulate PLP content, but BmNPV induces miR-6498-5p down-regulation to promote its proliferation. Our findings provide valuable insights into the role of host miRNA in B. mori-BmNPV interaction. Furthermore, the identification of the antiviral molecule PLP offers a novel perspective on strategies for preventing and managing viral infection in sericulture.

The Bombyx mori G protein β subunit 1 (BmGNβ1) gene inhibits BmNPV infection

G protein β subunit 1 (GNβ1) has several functions, including cell growth regulation, the control of second messenger levels, and ion channel switching. Previous transcriptome analyses in our laboratory have shown that BmGNβ1 transcription is reduced following infection with Bombyx mori nucleopolyhedrovirus (BmNPV), but it is unknown what role this gene may have in the host response to BmNPV infection. In this study, the BmGNβ1 gene was cloned using the RACE method. After BmNPV infection, BmGNβ1 was downregulated in Baiyu strains in tissues such as the hemolymph and midgut. Indirect immunofluorescence showed that BmGNβ1 was localized to the cytoplasm. We further constructed a BmGNβ1-pIZ/V5-His-mCherry overexpression plasmid and designed siRNA to evaluate the role of BmGNβ1 in host response to infection. The results showed that BmGNβ1 overexpression inhibited BmNPV proliferation, while knockdown of BmGNβ1 was correlated with increased BmNPV proliferation. The siRNA-mediated reduction of BmGNβ1 was correlated with an increase in BmNPV infection of BmN cells, increased BmNPV vp39 transcription, and reduced survival time of BmNPV-infected B. mori. Overexpression of BmGNβ1 in BmN cells was also correlated with apoptosis and a modification in transcript levels of genes involved in host response to BmNPV infection (PI3K, AKT, Bmp53, BmFOXO, Caspase-1, Bmp21, BmPKN and BmCREB), suggesting that BmGNβ1 may influence the apoptotic host response of infected B. mori through the PI3K-AKT pathway. This study provides potential targets and theoretical support for breeding BmNPV-resistant silkworm varieties.

An engineered Pichia pastoris platform for the biosynthesis of silk-based nanomaterials with therapeutic potential

Silk fibroin (SF) from the cocoon of silkworm has exceptional mechanical properties and biocompatibility and is used as a biomaterial in a variety of fields. Sustainable, affordable, and scalable manufacturing of SF would enable its large-scale use. We report for the first time the high-level secretory production of recombinant SF peptides in engineered Pichia pastoris cell factories and the processing thereof to nanomaterials. Two SF peptides (BmSPR3 and BmSPR4) were synthesized and secreted by P. pastoris using signal peptides and appropriate spacing between hydrophilic sequences. By strain engineering to reduce protein degradation, increase glycyl-tRNA supply, and improve protein secretion, we created the optimized P. pastoris chassis PPGSP-8 to produce BmSPR3 and BmSPR4. The SF fed-batch fermentation titers of the resulting two P. pastoris cell factories were 11.39 and 9.48 g/L, respectively. Protein self-assembly was inhibited by adding Tween 80 to the medium. Recombinant SF peptides were processed to nanoparticles (NPs) and nanofibrils. The physicochemical properties of nanoparticles R3NPs and R4NPs from the recombinant SFs synthesized in P. pastoris cell factories were similar or superior to those of RSFNPs (Regenerated Silk Fibroin NanoParticles) originating from commercially available SF. Our work will facilitate the production by microbial fermentation of functional SF for use as a biomaterial.

Transcriptome analysis indicates the mechanisms of BmNPV resistance in Bombyx mori midgut

Bombyx mori nucleopolyhedrovirus (BmNPV) caused serious economic losses in sericulture. Analyzing the molecular mechanism of silkworms (B. mori) resistance to BmNPV is of great significance for the prevention and control of silkworm virus diseases and the biological control of agricultural lepidopteran pests. In order to clarify the defense mechanisms of silkworms against BmNPV, we constructed a near isogenic line BC8 with high resistance to BmNPV through the highly BmNPV-resistant strain NB and the highly BmNPV-susceptible strain 306. In this study, RNA-Seq technique was used to analyze the transcriptome level differences in the midgut of BC8 and 306 following BmNPV infection. A total of 1350 DEGs were identified. Clustering analysis showed that these genes could be divided into 8 clusters with different expression patterns. Functional annotations based on GO and KEGG analysis indicated that they were involved in various metabolism pathways. Finally, 32 BmNPV defense responsive genes were screened. They were involved in metabolism, reactive oxygen species (ROS), signal transduction and immune response, and insect hormones. The further verification shows that HSP70 should participate in resistance responses of anti-BmNPV. These findings have paved the way in further functional characterization of candidate genes and subsequently can be used in breeding of BmNPV resistance dominant silkworms.

Effective expression and characterization of the receptor binding domains in SARS-CoV-2 Spike proteins from original strain and variants of concern using Bombyx mori nucleopolyhedrovirus in silkworm

A new coronavirus, known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is responsible for the global pandemic of COVID-19 in 2020. Through structural analysis, it was found that several amino acid residues in the human angiotensin-converting enzyme-2 (hACE2) receptor directly interact with those in the receptor binding domain (RBD) of the spike glycoprotein (S-protein). Various cell lines, including HEK293, HeLa cells, and the baculovirus expression vector system (BEVS) with the insect cell line Sf9, have been utilized to produce the RBD. In this study, we investigated the use of Bombyx mori nucleopolyhedrovirus (BmNPV) and BEVS. For efficient production of a highly pure recombinant RBD protein, we designed it with two tags (His tag and STREP tag) at the C-terminus and a solubilizing tag (SUMO) at the N-terminus. After expressing the protein using BmNPV and silkworm and purifying it with a HisTrap excel column, the eluted protein was digested with SUMO protease and further purified using a Strep-Tactin Superflow column. As a result, we obtained the RBD as a monomer with a yield of 2.6 mg/10 mL serum (equivalent to 30 silkworms). The RBD showed an affinity for the hACE2 receptor. Additionally, the RBDs from the Alpha, Beta, Gamma, Delta, and Omicron variants were expressed and purified using the same protocol. It was found that the RBD from the Alpha, Beta, Gamma, and Delta variants could be obtained with yields of 1.4-2.6 mg/10 mL serum and had an affinity to the hACE2 receptor.

A parasitoid regulates 20E synthesis and antibacterial activity of the host for development by inducing host nitric oxide production

Parasitoids are important components of the natural enemy guild in the biological control of insect pests. They depend on host resources to complete the development of a specific stage or whole life cycle and thus have evolved towards optimal host exploitation strategies. In the present study, we report a specific survival strategy of a fly parasitoid Exorista sorbillans (Diptera: Tachinidae), which is a potential biological control agent for agricultural pests and a pest in sericulture. We found that the expression levels of nitric oxide synthase (NOS) and nitric oxide (NO) production in host Bombyx mori (Lepidoptera: Bombycidae) were increased after E. sorbillans infection. Reducing NOS expression and NO production with an NOS inhibitor (NG-nitro-L-arginine methyl ester hydrochloride) in infected B. mori significantly impeded the growth of E. sorbillans larvae. Moreover, the biosynthesis of 20-hydroxyecdysone (20E) in infected hosts was elevated with increasing NO production, and inhibiting NOS expression lowered 20E biosynthesis. More importantly, induced NO synthesis was required to eliminate intracellular bacterial pathogens that presumably competed for shared host resources. Inhibiting NOS expression down-regulated the transcription of antimicrobial peptide genes and increased the number of bacteria in parasitized hosts. Collectively, this study revealed a new perspective on the role of NO in host-parasitoid interactions and a novel mechanism for parasitoid regulation of host physiology to support its development.

Identification and function of the Pax gene Bmgsb in the silk gland of Bombyx mori

Paired box (Pax) genes are highly conserved throughout evolution, and the Pax protein is an important transcription factor of embryonic development. The Pax gene Bmgsb is expressed in the silk glands of silkworm, but its biological functions remain unclear. This study aimed to investigate the expression pattern of Bmgsb in the silk gland and explore its functions using RNA interference (RNAi). Here, we identified eight Pax genes in Bombyx mori. Phylogenetic analysis showed that the B. mori Pax genes were highly homologous to the Pax genes in other insects and highly evolutionarily conserved. The tissue expression profile showed that Bmgsb was expressed in the anterior silk gland and anterior part of the middle silk gland (AMSG). RNAi of Bmgsb resulted in defective development of the AMSG, and the larvae were mostly unable to cocoon in the wandering stage. RNA-seq analysis showed that the fibroin genes fib-l, fib-h and p25, cellular heat shock response-related genes and phenol oxidase genes were considerably upregulated upon Bmgsb knockdown. Furthermore, quantitative reverse transcription-PCR results showed that the fibroin genes and ubiquitin proteolytic enzyme-related genes were significantly upregulated in the AMSG after Bmgsb knockdown. This study provides a foundation for future research on the biological functions of B. mori Pax genes. In addition, it demonstrates the important roles of Bmgsb in the transcriptional regulation of fibroin genes and silk gland development.

Contrasting gut bacteriomes unveiled between wild Antheraea assamensis Helfer (Lepidoptera: Saturniidae) and domesticated Bombyx mori L. (Lepidoptera: Bombycidae) silkworms

BACKGROUND

Insect gut microbiomes play a fundamental role in various aspects of insect physiology, including digestion, nutrient metabolism, detoxification, immunity, growth and development. The wild Muga silkworm, Antheraea assamensis Helfer holds significant economic importance, as it produces golden silk.

METHODS AND RESULTS

In the current investigation, we deciphered its intricate gut bacteriome through high-throughput 16S rRNA amplicon sequencing. Further, to understand bacterial community dynamics among silkworms raised under outdoor environmental conditions, we compared its gut bacteriomes with those of the domesticated mulberry silkworm, Bombyx mori L. Most abundant bacterial phyla identified in the gut of A. assamensis were Proteobacteria (78.1%), Bacteroidetes (8.0%) and Firmicutes (6.6%), whereas the most-abundant phyla in B. mori were Firmicutes (49-86%) and Actinobacteria (10-36%). Further, Gammaproteobacteria (57.1%), Alphaproteobacteria (10.47%) and Betaproteobacteria (8.28%) were the dominant bacterial classes found in the gut of A. assamensis. The predominant bacterial families in A. assamensis gut were Enterobacteriaceae (27.7%), Comamonadaceae (9.13%), Pseudomonadaceae (9.08%) Flavobacteriaceae (7.59%) Moraxellaceae (7.38%) Alteromonadaceae (6.8%) and Enterococcaceae (4.46%). In B. mori, the most-abundant bacterial families were Peptostreptococcaceae, Enterococcaceae, Lactobacillaceae and Bifidobacteriaceae, though all showed great variability among the samples. The core gut bacteriome of A. assamensis consisted of Pseudomonas, Acinetobacter, Variovorax, Myroides, Alteromonas, Enterobacter, Enterococcus, Sphingomonas, Brevundimonas, Oleispira, Comamonas, Oleibacter Vagococcus, Aminobacter, Marinobacter, Cupriavidus, Aeromonas, and Bacillus. Comparative gut bacteriome analysis revealed a more complex gut bacterial diversity in wild A. assamensis silkworms than in domesticated B. mori silkworms, which contained a relatively simple gut bacteriome as estimated by OTU richness. Predictive functional profiling of the gut bacteriome suggested that gut bacteria in A. assamensis were associated with a wide range of physiological, nutritional, and metabolic functions, including biodegradation of xenobiotics, lipid, amino acid, carbohydrate metabolism, and biosynthesis of secondary metabolites and amino acids.

CONCLUSIONS

These results showed great differences in the composition and diversity of gut bacteria between the two silkworm species. Both insect species harbored core bacterial taxa commonly found in insects, but the relative abundance and composition of these taxa varied markedly.

The mechanoreceptor Piezo is required for spermatogenesis in Bombyx mori

BACKGROUND

The animal sperm shows high diversity in morphology, components, and motility. In the lepidopteran model insect, the silkworm Bombyx mori, two types of sperm, including nucleate fertile eupyrene sperm and anucleate unfertile apyrene sperm, are generated. Apyrene sperm assists fertilization by facilitating the migration of eupyrene spermatozoa from the bursa copulatrix to the spermatheca. During spermatogenesis, eupyrene sperm bundles extrude the cytoplasm by peristaltic squeezing, while the nuclei of the apyrene sperm bundles are discarded with the same process, forming matured sperm.

RESULTS

In this study, we describe that a mechanoreceptor BmPiezo, the sole Piezo ortholog in B. mori, plays key roles in larval feeding behavior and, more importantly, is essential for eupyrene spermatogenesis and male fertility. CRISPR/Cas9-mediated loss of BmPiezo function decreases larval appetite and subsequent body size and weight. Immunofluorescence analyses reveal that BmPiezo is intensely localized in the inflatable point of eupyrene sperm bundle induced by peristaltic squeezing. BmPiezo is also enriched in the middle region of apyrene sperm bundle before peristaltic squeezing. Cytological analyses of dimorphic sperm reveal developmental arrest of eupyrene sperm bundles in BmPiezo mutants, while the apyrene spermatogenesis is not affected. RNA-seq analysis and q-RT-PCR analyses demonstrate that eupyrene spermatogenic arrest is associated with the dysregulation of the actin cytoskeleton. Moreover, we show that the deformed eupyrene sperm bundles fail to migrate from the testes, resulting in male infertility due to the absence of eupyrene sperm in the bursa copulatrix and spermatheca.

CONCLUSIONS

In conclusion, our studies thus uncover a new role for Piezo in regulating spermatogenesis and male fertility in insects.

Role of small nucleolar RNAs in alternative splicing of the doublesex gene in the silkworm, Bombyx mori

More and more evidence shows that small noncoding RNAs (ncRNAs) play diverse roles in development, stress response and other cellular processes, but functional study of intermediate-size ncRNAs is still rare. Here, the expression profile of 16 intermediate-size ncRNAs in ovary and testis of silkworm Bombyx mori were analyzed. Twelve ncRNAs, including 5 small nucleolar RNAs (snoRNAs) and 7 unclassified ncRNAs, accumulated more in the testis than in the ovary of silkworm, especially Bm-163, Bm-51 and Bm-68. Four ncRNAs (including three orphan snoRNAs and one unclassified ncRNA) had higher expression level in the ovary than in the testis, especially Bm-86. Overexpression of the testis-enriched snoRNA Bm-68 in the female led to the accumulation of male-specific isoform of doublesex (BmdsxM) and increased the expression ratio of BmdsxM: BmdsxF. While overexpression of ovary-enriched snoRNA Bm-86 in the male decreased the expression ratio of BmdsxM: BmdsxF, indicating the roles of the two snoRNAs played in the alternative splicing of Bmdsx of silkworm, which will provide new clues for the functional study of snoRNAs in insects.

The dual role of Spn-E in supporting heterotypic ping-pong piRNA amplification in silkworms

The PIWI-interacting RNA (piRNA) pathway plays a crucial role in silencing transposons in the germline. piRNA-guided target cleavage by PIWI proteins triggers the biogenesis of new piRNAs from the cleaved RNA fragments. This process, known as the ping-pong cycle, is mediated by the two PIWI proteins, Siwi and BmAgo3, in silkworms. However, the detailed molecular mechanism of the ping-pong cycle remains largely unclear. Here, we show that Spindle-E (Spn-E), a putative ATP-dependent RNA helicase, is essential for BmAgo3-dependent production of Siwi-bound piRNAs in the ping-pong cycle and that this function of Spn-E requires its ATPase activity. Moreover, Spn-E acts to suppress homotypic Siwi-Siwi ping-pong, but this function of Spn-E is independent of its ATPase activity. These results highlight the dual role of Spn-E in facilitating proper heterotypic ping-pong in silkworms.

Bombyx moriSuppressor of Hairless is involved in the regulation of the silkworm cell cycle and endoreplication of the silk glands

The Notch signaling pathway is important in cell cycle regulation and cell proliferation. The transcriptional repressor Suppressor of Hairless [Su(H)] is a molecular switch for downstream target genes of the Notch signaling pathway but the regulatory mechanism of the Su(H) gene in the cell cycle is unclear. We determined the function of the Notch signaling pathway and Bombyx mori Su(H) [BmSu(H)] in the regulation of the silkworm cell cycle. Inhibition of Notch signaling promoted the replication of DNA in silkworm gland cells and expression of the BmSu(H) gene was significantly reduced. Overexpression of the BmSu(H) gene inhibited DNA replication and cell proliferation of silkworm cells, whereas knockout of the BmSu(H) gene promoted DNA replication and cell proliferation. Knockout of the BmSu(H) in silkworms improved the efficiency of silk gland cell endoreplication and increased important economic traits. We demonstrated that BmSu(H) protein can directly bind to the promoters of BmCyclinA, BmCyclinE and BmCDK1 genes, inhibiting or promoting their transcription at the cell and individual level. This study identified molecular targets for genetic improvement of the silkworm and also provided insights into the regulatory mechanism of the cell cycle.

Biochemical characterization of Bombyx mori Dicer-2 that dices double-stranded RNAs into 20-nt small RNA

We detected enzymatic activity that generates 20-nucleotide (nt) RNA from double-stranded RNAs (dsRNAs) in crude extracts prepared from various silkworm (Bombyx mori) organs. The result using knocked-down cultured cells indicated that this dicing activity originated from B. mori Dicer-2 (BmDcr2). Biochemical analyses revealed that BmDcr2 preferentially cleaves 5'-phosphorylated dsRNAs at the 20-nt site-counted from the 5'-phosphorylated end-and required ATP and magnesium ions for the dicing reaction. This is the first report of the biochemical characterization of Dicer-2 in lepidopteran insects. This enzymatic property of BmDcr2 in vitro is consistent with the in vivo small interfering RNA profile in virus-infected silkworm cells.

[Identification and expression pattern analysis of a secretory phospholipase A2 gene in Bombyx mori]

Phospholipase A2 (PLA2) is widely distributed in animals, plants, and microorganisms, and it plays an important role in many physiological activities. In a previous study, we have identified a secretory PLA2 in Bombyx mori (BmsPLA2-1-1). In this study, we further identified four new sPLA2 genes (BmsPLA2-1-2, BmsPLA2-2, BmsPLA2-3, and BmsPLA2-4) in B. mori genome. All four genes exhibits the characteristic features of sPLA2, including the sPLA2 domain, metal binding sites, and highly conserved catalytic domain. This study completed the cloning, in vitro expression, and expression pattern analysis of the BmsPLA2-4 gene in B. mori. The full length of BmsPLA2-4 is 585 bp, and the recombinant protein obtained through prokaryotic expression has an estimated size of 25 kDa. qRT-PCR analysis revealed that the expression level of BmsPLA2-4 reached its peak on the first day of the fifth instar larval stage. Tissue expression profiling analysis showed that BmsPLA2-4 had the highest expression level in the midgut, followed by the epidermis and fat body. Western blotting analysis results were consistent with those of qRT-PCR. Furthermore, after infecting fifth instar 1-day-old larvae with Escherichia coli and Staphylococcus aureus, the expression level of the BmsPLA2-4 gene significantly increased in 24 h. The findings of this study provides a theoretical basis and valuable experimental data for future related research.

The BTB-ZF gene Bm-mamo regulates pigmentation in silkworm caterpillars

The color pattern of insects is one of the most diverse adaptive evolutionary phenotypes. However, the molecular regulation of this color pattern is not fully understood. In this study, we found that the transcription factor Bm-mamo is responsible for black dilute (bd) allele mutations in the silkworm. Bm-mamo belongs to the BTB zinc finger family and is orthologous to mamo in Drosophila melanogaster. This gene has a conserved function in gamete production in Drosophila and silkworms and has evolved a pleiotropic function in the regulation of color patterns in caterpillars. Using RNAi and clustered regularly interspaced short palindromic repeats (CRISPR) technology, we showed that Bm-mamo is a repressor of dark melanin patterns in the larval epidermis. Using in vitro binding assays and gene expression profiling in wild-type and mutant larvae, we also showed that Bm-mamo likely regulates the expression of related pigment synthesis and cuticular protein genes in a coordinated manner to mediate its role in color pattern formation. This mechanism is consistent with the dual role of this transcription factor in regulating both the structure and shape of the cuticle and the pigments that are embedded within it. This study provides new insight into the regulation of color patterns as well as into the construction of more complex epidermal features in some insects.

Harvesting Silk Fibers from Plodia interpunctella: Role of Environmental Rearing Conditions in Fiber Production and Properties

Silk fibers are produced by a wide variety of insects. The silkworm Bombyx mori (Bombyx) was domesticated because the physical properties of its silk fibers were amenable to the production of fine textiles. Subsequently, engineers have regenerated silk fibroin to form biomaterials. The monocular focus on Bombyx silk has underutilized the expanse of diverse silk proteins produced by more than 100,000 other arthropods. This vast array of silk fibers could be utilized for biomedical engineering challenges if sufficient rearing and purification processes are developed. Herein, we show that the moth, Plodia interpunctella (Plodia), represents an alternative silk source that is easily reared in highly regulated culture environments allowing for greater consistency in the silk produced. We controlled the temperature, resource availability (larvae/gram diet), and population density (larvae/mL) with the goal of increasing silk fiber production and improving homogeneity in Plodia silk proteins. We determined that higher temperatures accelerated insect growth and reduced life cycle length. Furthermore, we established initial protocols for the production of Plodia silk with optimal silk production occurring at 24 °C, with a resource availability of 10 larvae/gram and a population density of 0.72 larvae/mL. Population density was shown to be the most prominent driving force of Plodia silk mat formation among the three parameters assessed. Future work will need to link gene expression, protein production and purification, and resulting mechanical properties as a function of environmental cues to further transition Plodia silk into regenerated silk fibroin biomaterials.