Identification and Characterization of 30 K Protein Genes Found in Bombyx mori (Lepidoptera: Bombycidae) Transcriptome

Evaluating bio-physicochemical properties of raw powder prepared from whole larvae containing liquid silk of the domestic silkworm

The domestic silkworm, Bombyx mori, has been widely used in silk production for centuries. It is also used as a bioreactor by the textile and pharmaceutical industries to mass produce recombinant bioactive proteins containing silk-based materials. Furthermore, silkworms are well-known as a source of food and have also been orally administered to prevent and treat several human disorders. In this study, we aimed to investigate the inherent bio-physicochemical properties of edible silkworms to accurately evaluate their clinical and nutritional potential. We prepared raw powder from whole larvae of silkworm. The yield rate of the powder derived from dried larvae was almost 100% (98.1-99.1% in replicates). As "percentage yield" translates to "Budomari" in Japanese, this raw powder was named "B100rw." We further prepared B100dn that was denatured through autoclaving. Thereafter, we examined whether B100rw sustained the original bio-physicochemical properties by comparing it with B100dn. There was no significant difference in nutritional content between B100rw and B100dn. B100rw contained proteins derived from silkworm larvae and mulberry leaves, whereas the proteins of B100dn were mostly degraded. On measuring the enzymatic activity of both powders using trehalase as an indicator enzyme, B100rw was found to maintain trehalase activity. B100rw also maintained a random coil conformation, similar to that of liquid silk. This suggested that B100rw sustained the unique bio-physicochemical properties of living larvae. These findings may facilitate the development of novel food products or orally administered vaccines.

BmNPV Bm60 is a key target gene used by a resistant strain of Bombyx mori to inhibit BmNPV proliferation

Bombyx mori nucleopolyhedrovirus (BmNPV) is a pathogen that causes significant losses to the silkworm industry. Numerous antiviral genes and proteins have been identified by studying silkworm resistance to BmNPV. However, the molecular mechanism of silkworm resistance to BmNPV is unclear. We analyzed the differences between the susceptible strain 871 and a near-isogenic resistant strain 871C. The survival of strain 871C was significantly greater than that of 871 after oral and subcutaneous exposure to BmNPV. Strain 871C exhibited a nearly 10,000-fold higher LD50 for BmNPV compared to 871. BmNPV proliferation was significantly inhibited in all tested tissues of strain 871C using HE strain and fluorescence analysis. Strain 871C exhibited cellular resistance to BmNPV rather than peritrophic membrane or serum resistance. Strain 871C suppressed the expression of the viral early gene Bm60. This led to the inhibition of BmNPV DNA replication and late structural gene transcription based on the cascade regulation of baculovirus gene expression. Bm60 could also interact with the viral DNA binding protein and alkaline nuclease, as well as host proteins Methylcrotonoyl-CoA carboxylase subunit alpha, mucin-2-like protein, and 30 K-8. Overexpression of 30 K-8 significantly inhibited BmNPV proliferation. These results increase understanding of the molecular mechanism behind silkworm resistance to BmNPV and suggest targets for the breeding of resistant silkworm strains and the controlling pest of Lepidoptera.

Influence of temperature variation on gene expression and cocoon production in Bombyx mori Linnaeus, 1758 (Lepidoptera: Bombycidae)

Silkworms (Bombyx mori) are lepidopterans of economic importance for global silk production. However, factors that directly affect the yield and quality of silkworm cocoon production, such as diseases and temperature fluctuations, cause great economic losses. Knowing how they respond to rearing temperature during the most critical stage of their life cycle (i.e., fifth instar) could provide information on their adaptation and improve silk production. In the current work, we analyzed transcriptional data from two groups of B. mori that were reared at 26 °C and 34 °C throughout the fifth instar. The silkworms and cocoons were weighed. In total, 3115 transcripts were differentially expressed (DE; including 1696 down-regulated and 1419 up-regulated) among the 29,157 sequences found by transcriptome assembly. We emphasize the genes associated with immunological response, transcription factors, silk biosynthesis, and heat shock proteins, among the DE transcripts in response to the temperature conditions. Silkworms reared at 34 °C presented a reduced mean body weight (-0.944 g in comparison to the 26 °C group), which had a direct impact on the weight of cocoons formed and the silk conversion rate. These changes were statistically significant when compared to silkworms reared at 26 °C. Mortality rates (6 and 9 %, at 26 °C and 34 °C, respectively) were similar to those obtained in breeding fields. The findings provide information on the biological processes involved in the temperature response mechanism of silkworms, as well as information that may be used in future climatization processes at rearing facilities and in breeding for improved thermotolerance.

LP30K protein manifested in hemocytes of Bombyx mori larva on Nosema bombycis infection and showed functional evolution based on glucose- binding domain

Infection by microsporidian Nosema bombycis induced appearance of exclusive protein conjugate of 190 kDa in hemocytes of silkworm Bombyx mori L (Lepidoptera: Bombycidae). Mass spectrometry of the band showed peptides of low molecular weight 30 kDa lipoprotein (LP30K). Six accessions of LP30K identified from the hemocytes comprised 30 K lipoprotein 1, 30 K protein 1, 2, 6, 7 and 11. Two uncharacterised proteins (UCP) identified from the hemocytes showed 100% similarity with LP30K sequence, altogether showed abundance after the infection. The LP30K accessions H9J4F6 (Q00802), E5EVW2 and the UCP accessions D4QGC0 and D4QGB9 showed presence of glucose binding protein I domain "ADSDVPNDILEEQLYNSIVVADYDSAVEK" that binds with fungal glucans to inhibit infection. However glucose binding protein II domain "TLAPRTDDVLAEQLYMSVVIGEYETAIAK" is absent in LP30K accessions from hemocytes showed loss of DNA sequences encoding the domain. The accessions H9J4F5, H9B440, A7LIK7 and H9B444 showed 92% identity with B. mori LP30K protein (NP_001095198.2) however the glucose binding domain I is absent in these accessions suggesting isoform- specific restricted fungal defense activity. Phylogeny tree of the LP30K homologues showed four groups including microvitellogenin and 30 kDa proteins showing functional diversity endorsed with evolutionary diversity. LP30K accessions with glucose binding domain diverged from that without glucose binding domain exemplify co-evolution for domain- dependent functional roles like storage and immune reactions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03685-x.

Purification, Identification and Neuroprotective Effects of Proteins from Bombyx batryticatus in Glu-Stimulated PC12 Cells

Bombyx batryticatus (BB) is one of the most commonly used Traditional Chinese Medicines (TCMs) in the treatment of convulsions and epilepsy. The antiepileptic effects of total proteins from BB (BBPs) have been proven in our previous research. In this study, BBPs were further purified, the neuroprotective effects were evaluated in Glu-stimulated PC12 cells, and the structure was identified by Matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF-MS). Six subfractions (PF-1 to PF-6) were obtained by DEAE-52 Sepharose FF ion-exchange chromatography. It was found that PF-1, PF-2, and PF-3, with similar protein compositions, possessed neuroprotective effects in Glu-stimulated PC12 cells by significantly increasing the GABA level, and decreasing the levels of IL-1β and TNF-α. The most active fraction (PF-2) was further separated by Sephadex G-75 gel filtration chromatography, and an effective protein component named PF-2-2 was obtained. Fluorescein isothiocyanate-labeled PF-2-2 (FITC-PF-2-2) was prepared, and the binding of FITC-F-2-2 to the PC12 cells was directly observed with a confocal microscope. PF-2-2 was found to first bind to the surface of PC12 cells and then internalize into the cells. The main band of PF-2-2 was then analyzed by MALDI-TOF/TOF-MS and searched in the MASCOT database; finally a protein named Low molecular mass 30 kDa lipoprotein 21G1 was identified. In conclusion, PF-2-2 and purified proteins isolated from BBPs have potential application prospects in the treatment of epilepsy.

Global Profiling of Genes Expressed in the Silk Glands of Philippine-Reared Mulberry Silkworms (Bombyx mori)

RNA sequencing was used to assemble transcriptome data for Philippine-reared silkworm and compare gene expression profiles of strains reared in high- and low-temperature environments. RNA was isolated from the silk glands of fifth instar larvae and mRNA-enriched libraries were sequenced using Illumina NextSeq 500. Transcriptome reads were assembled using reference-based and de novo assemblers, and assemblies were evaluated using different metrics for transcriptome quality, including the read mapping rate, E90N50, RSEM-eval, and the presence of single-copy orthologs. All transcriptome assemblies were able to reconstruct >40,000 transcripts. Differential expression analysis found 476 differentially expressed genes (DEGs; 222 upregulated, 254 downregulated) in strains reared in different temperatures. Among the top DEGs were myrosinase, heat shock proteins, serine protease inhibitors, dehydrogenases, and regulators of the juvenile hormone. Validation of some of the top DEGs using qPCR supported the findings of the in silico analysis. GO term enrichment analysis reveals an overrepresentation of GO terms related to nucleotide metabolism and biosynthesis, lipid and carbohydrate metabolic processes, regulation of transcription, nucleotide binding, protein binding, metal binding, catalytic activity, oxidoreductase activity, and hydrolase activity. The data provided here will serve as a resource for improving local strains and increasing silk production of Philippine-reared B. mori strains.

Expression and Role of Vitellogenin Genes in Ovarian Development of Zeugodacus cucurbitae

Vitellogenin (Vg) genes encode the major egg yolk protein precursor in arthropods. In this study, four Vgs were identified in Zeugodacus cucurbitae (Coquillett). Sequence analysis showed that four ZcVgs had the conserved Vg domain. Phylogenetic analysis indicated that four ZcVgs were homologous to the Vgs of Tephritidae insects. The temporal and spatial expression patterns of ZcVgs were analyzed by quantitative real-time polymerase chain reaction (RT-qPCR), and the four ZcVgs showed high expression levels in female adults, especially in the fat body. The expression of ZcVg1 and ZcVg3 was down-regulated by a low dosage (0.5 μg) of 20-hydroxyecdysone (20E), and ZcVg2, ZcVg3, and ZcVg4 were up-regulated by a high dosage (1.0 and 2.0 μg) of 20E. The expression of ZcVg1 and ZcVg2 was up-regulated by 5 μg of juvenile hormone (JH), while all of the ZcVgs were down-regulated by a low and high dosage of JH. Expression of ZcVgs was down-regulated after 24 h of starvation and recovered to normal after nutritional supplementation. After micro-injection of the gene-specific double-stranded RNA, the ZcVgs’ expression was significantly suppressed, and ovarian development was delayed in Z. cucurbitae females. The results indicate that RNA interference of reproduction-related genes is a potential pest control method that works by manipulating female fertility.

Allergenic characterization of Bomb m 4, a 30 kDa lipoprotein Bombyx mori lipoprotein 6 from silkworm pupa

BACKGROUND

Silkworm pupa (SWP) food anaphylaxis has been described frequently in Asian countries. However, false-positive reactions by skin pricks and serum IgE (sIgE) tests to the extract complicate diagnosis, requiring identification of clinically relevant major allergens.

OBJECTIVES

In this study, we characterized a novel SWP allergen, Bomb m 4, a 30 kDa lipoprotein, and evaluated its diagnostic sensitivity.

METHODS

Bomb m 4 was identified by a proteomic analysis. This recombinant (r)Bomb m 4 was overexpressed in Escherichia coli, and the IgE reactivity by ELISA was compared with other reported allergenic proteins: Bomb m 1 (arginine kinase), 27 kDa glycoprotein, Bomb m 3 (tropomyosin) using the serum samples from 17 SWP allergy patients and 11 asymptomatic sensitized subjects.

RESULTS

rBomb m 4 specific IgE was recognized by all 17 SWP allergy patients. The 27 kDa glycoprotein and Bomb m 1 sIgE were found in 35.3% and 0%, respectively in the SWP allergic patients. ELISA sIgE reactivity increased significantly, when 4 M urea was added in serum samples. However, only 16% inhibition of sIgE reactivity to the whole SWP extract was exhibited by rBomb m 4, whereas more than 93% of self-inhibition of rBomb m 4 sIgE was obtained, possibly due to the low abundance of Bomb m 4 in the extract. Three linear epitopes (81-95, 191-205, and 224-238 residues) of rBomb m 4 were identified. These epitopes are shown to be released by pepsin digestion. ROC analysis showed the highest diagnostic value of Bomb m 4 followed by Bomb m 1, 27 kDa glycoprotein, and Bomb m 3.

CONCLUSION

Bomb m 4 is the major allergen of SWP allergy patients. It has cryptic epitopes which are exposed to IgE antibodies with digestive enzymes. This recombinant Bomb m 4 allergen permits exact diagnosis of SWP allergy.

The acetylation modification regulates the stability of Bm30K-15 protein and its mechanism in silkworm, Bombyx mori

The 30 K proteins are the major silkworm hemolymph proteins and are involved in a variety of physiological processes, such as nutrient and energy storage, embryogenesis, immune response, and inhibition of apoptosis. The Bm30K-15 protein is one of the 30 K proteins and is abundant in the hemolymph of fifth instar silkworm larva. We previously found that the Bm30K-15 protein can be acetylated. In the present study, we found that acetylation can improve the protein stability of Bm30K-15. Further exploration confirmed that the increase in protein stability by acetylation was caused by competition between acetylation and ubiquitination. In summary, these findings aim to provide insight into the effect of acetylation modification on the protein level and stability of the Bm30K-15 and the possible molecular mechanism of its existence in silkworm, Bombyx mori.

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

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

Genome sequence and organization of the Mythimna (formerly Pseudaletia) unipuncta granulovirus Hawaiian strain

Purified occlusion bodies (OBs) of Mythimna (formerly Pseudaletia) unipuncta (the true armyworm) granulovirus Hawaiian strain (MyunGV-A) were observed, showing typical GV morphological characteristics under scanning and transmission electron microscopy (EM). The genome of MyunGV-A was completely sequenced and analysed. The genome is 176,677 bp in size, with a G+C content of 39.79%. It contains 183 open reading frames (ORFs) encoding 50 or more amino acids with minimal overlap. Comparison of MyunGV-A with TnGV, XcGV, and HearGV genomes revealed extensive sequence similarity and collinearity, and the four genomes contain the same nine homologous regions (hrs) with conserved structures and locations. Three unique genes, 12 baculovirus repeated ORF (bro), 2 helicase, and 3 enhancin genes, were identified. In particular, two repeated genes (ORF39 and 49) are present in the genome, in reverse and complementarily orientations. Twenty-four OB proteins were identified from the putative protein database of MyunGV-A. In addition, MyunGV-A belongs to the Betabaculovirus group and is most closely related to TnGV (99% amino acid identity) according to a phylogenetic tree based on the combined amino acid sequences of 38 core gene contents.

The stability and antiapoptotic activity of Bm30K-3 can be improved by lysine acetylation in the silkworm, Bombyx mori

Acetylation is an important, highly conserved, and reversible post-translational modification of proteins. Previously, we showed by nano-HPLC/MS/MS that many nutrient storage proteins in the silkworm are acetylated. Among these proteins, most of the known 30K proteins were shown to be acetylated, including 23 acetylated 30K proteins containing 49 acetylated sites (Kac), indicating the importance of the acetylation of 30K proteins in silkworm. In this study, Bm30K-3, a 30K protein containing three Kac sites, was further assessed in functional studies of its acetylation. Increasing the level of Bm30K-3 acetylation by adding the deacetylase inhibitor trichostatin A (TSA) increased the levels of this protein and further inhibited cellular apoptosis induced by H₂ O₂ . In contrast, decreasing the level of acetylation by adding the acetylase inhibitor C646 could reduce the level of Bm30K-3 and increase H₂ O₂ -induced apoptosis. Subsequently, BmN cells were treated with CHX and MG132, and increasing the acetylation level using TSA was shown to inhibit protein degradation and improve the stability of Bm30K-3. Furthermore, the acetylation of Bm30K-3 could compete with its ability to be ubiquitinated, suggesting that acetylation could inhibit the ubiquitin-mediated proteasome degradation pathway, improving the stability and accumulation of proteins in cells. These results further indicate that acetylation might regulate nutrition storage and utilization in Bombyx mori, which requires further study.

Identification of Genes that Control Silk Yield by RNA Sequencing Analysis of Silkworm (Bombyx mori) Strains of Variable Silk Yield

Silk is an important natural fiber of high economic value, and thus genetic study of the silkworm is a major area of research. Transcriptome analysis can provide guidance for genetic studies of silk yield traits. In this study, we performed a transcriptome comparison using multiple silkworms with different silk yields. A total of 22 common differentially expressed genes (DEGs) were identified in multiple strains and were mainly involved in metabolic pathways. Among these, seven significant common DEGs were verified by quantitative reverse transcription polymerase chain reaction, and the results coincided with the findings generated by RNA sequencing. Association analysis showed that BGIBMGA003330 and BGIBMGA005780 are significantly associated with cocoon shell weight and encode uridine nucleosidase and small heat shock protein, respectively. Functional annotation of these genes suggest that these play a role in silkworm silk gland development or silk protein synthesis. In addition, we performed principal component analysis (PCA) in combination with wild silkworm analysis, which indicates that modern breeding has a stronger selection effect on silk yield traits than domestication, and imply that silkworm breeding induces aggregation of genes related to silk yield.

Bombyx mori and Aedes aegypti form multi-functional immune complexes that integrate pattern recognition, melanization, coagulants, and hemocyte recruitment

The innate immune system of insects responds to wounding and pathogens by mobilizing multiple pathways that provide both systemic and localized protection. Key localized responses in hemolymph include melanization, coagulation, and hemocyte encapsulation, which synergistically seal wounds and envelop and destroy pathogens. To be effective, these pathways require a targeted deposition of their components to provide protection without compromising the host. Extensive research has identified a large number of the effectors that comprise these responses, but questions remain regarding their post-translational processing, function, and targeting. Here, we used mass spectrometry to demonstrate the integration of pathogen recognition proteins, coagulants, and melanization components into stable, high-mass, multi-functional Immune Complexes (ICs) in Bombyx mori and Aedes aegypti. Essential proteins common to both include phenoloxidases, apolipophorins, serine protease homologs, and a serine protease that promotes hemocyte recruitment through cytokine activation. Pattern recognition proteins included C-type Lectins in B. mori, while A. aegypti contained a protein homologous to Plasmodium-resistant LRIM1 from Anopheles gambiae. We also found that the B. mori IC is stabilized by extensive transglutaminase-catalyzed cross-linking of multiple components. The melanization inhibitor Egf1.0, from the parasitoid wasp Microplitis demolitor, blocked inclusion of specific components into the IC and also inhibited transglutaminase activity. Our results show how coagulants, melanization components, and hemocytes can be recruited to a wound surface or pathogen, provide insight into the mechanism by which a parasitoid evades this immune response, and suggest that insects as diverse as Lepidoptera and Diptera utilize similar defensive mechanisms.