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.

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.

[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.

Utilization of Diverse Molecules as Receptors by Cry Toxin and the Promiscuous Nature of Receptor-Binding Sites Which Accounts for the Diversity

By 2013, it had been shown that the genes cadherin-like receptor (Cad) and ATP-binding cassette transporter subfamily C2 (ABCC2) were responsible for insect resistance to several Cry1A toxins, acting as susceptibility-determining receptors, and many review articles have been published. Therefore, this review focuses on information about receptors and receptor-binding sites that have been revealed since 2014. Since 2014, studies have revealed that the receptors involved in determining susceptibility vary depending on the Cry toxin subfamily, and that binding affinity between Cry toxins and receptors plays a crucial role. Consequently, models have demonstrated that ABCC2, ABCC3, and Cad interact with Cry1Aa; ABCC2 and Cad with Cry1Ab and Cry1Ac; ABCC2 and ABCC3 with Cry1Fa; ABCB1 with Cry1Ba, Cry1Ia, Cry9Da, and Cry3Aa; and ABCA2 with Cry2Aa and Cry2Ba, primarily in the silkworm, Bombyx mori. Furthermore, since 2017, it has been suggested that the binding sites of BmCad and BmABCC2 on Cry1Aa toxin overlap in the loop region of domain II, indicating that Cry toxins use various molecules as receptors due to their ability to bind promiscuously in this region. Additionally, since 2017, several ABC transporters have been identified as low-efficiency receptors that poorly induce cell swelling in heterologously expressing cultured cells. In 2024, research suggested that multiple molecules from the ABC transporter subfamily, including ABCC1, ABCC2, ABCC3, ABCC4, ABCC10, and ABCC11, act as low-efficiency receptors for a single Cry toxin in the midgut of silkworm larvae. This observation led to the hypothesis that the presence of such low-efficiency receptors contributes to the evolution of Cry toxins towards the generation of highly functional receptors that determine the susceptibility of individual insects. Moreover, this evolutionary process is considered to offer valuable insights for the engineering of Cry toxins to overcome resistance and develop countermeasures against resistance.

Masculinizer gene controls sexual differentiation in Hyphantria cunea

The Masculinizer gene, Masc, encodes a lepidopteran-specific novel CCCH-type zinc finger protein, which controls sex determination and dosage compensation in Bombyx mori. Considering the potential application of it in pest control, it is necessary to investigate the function of Masc gene in Hyphantria cunea, a globally invasive forest pest. In the present study, we identified and functionally characterized the Masc gene, HcMasc, in H. cunea. Sequence analysis revealed that HcMASC contained the conserved CCCH-type zinc finger domain, nuclear localization signal, and male determining domain, in which the last was confirmed to be required for its masculinization in BmN cell line. However, expression data showed that unlike male-biased expression in B. mori, HcMasc gene expresses in main all developmental stages or tissues in both sexes. Clustered regularly interspaced palindromic repeats (CRISPR) / CRISPR-associated protein 9-based disruption of the common exons 1 and 3 of the HcMasc gene resulted in imbalanced sex ratio and abnormal external genitalia of both sexes. Our results suggest that the HcMasc gene is required for both male and female sexual differentiation and dosage compensation in H. cunea and provide a foundation for developing better strategies to control this pest.

Targeted DNA N6-methyladenine editing by dCas9 fused to METTL4 in the lepidopteran model insect Bombyx mori

We have established a novel CRISPR-dCas9-METTL4 epigenome editing tool that can methylate target regions to achieve site-specific DNA 6mA methylation in both hypermethylated and hypomethylated genes. Targeted methylation on genes by dCas9-METTL4 results in misexpression, allowing for the functional investigation of target genes of interest in silkworm.

Transcriptomic changes in the microsporidia proliferation and host responses in congenitally infected embryos and larvae

Congenital infection caused by vertical transmission of microsporidia N. bombycis can result in severe economic losses in the silkworm-rearing industry. Whole-transcriptome analyses have revealed non-coding RNAs and their regulatory networks in N. bombycis infected embryos and larvae. However, transcriptomic changes in the microsporidia proliferation and host responses in congenitally infected embryos and larvae remains unclear. Here, we simultaneously compared the transcriptomes of N. bombycis and its host B. mori embryos of 5-day and larvae of 1-, 5- and 10-day during congenital infection. For the transcriptome of N. bombycis, a comparison of parasite expression patterns between congenital-infected embryos and larva showed most genes related to parasite central carbon metabolism were down-regulated in larvae during infection, whereas the majority of genes involved in parasite proliferation and growth were up-regulated. Interestingly, a large number of distinct or shared differentially expressed genes (DEGs) were revealed by the Venn diagram and heat map, many of them were connected to infection related factors such as Ricin B lectin, spore wall protein, polar tube protein, and polysaccharide deacetylase. For the transcriptome of B. mori infected with N. bombycis, beyond numerous DEGs related to DNA replication and repair, mRNA surveillance pathway, RNA transport, protein biosynthesis, and proteolysis, with the progression of infection, a large number of DEGs related to immune and infection pathways, including phagocytosis, apoptosis, TNF, Toll-like receptor, NF-kappa B, Fc epsilon RI, and some diseases, were successively identified. In contrast, most genes associated with the insulin signaling pathway, 2-oxacarboxylic acid metabolism, amino acid biosynthesis, and lipid metabolisms were up-regulated in larvae compared to those in embryos. Furthermore, dozens of distinct and three shared DEGs that were involved in the epigenetic regulations, such as polycomb, histone-lysine-specific demethylases, and histone-lysine-N-methyltransferases, were identified via the Venn diagram and heat maps. Notably, many DEGs of host and parasite associated with lipid-related metabolisms were verified by RT-qPCR. Taken together, simultaneous transcriptomic analyses of both host and parasite genes lead to a better understanding of changes in the microsporidia proliferation and host responses in embryos and larvae in N. bombycis congenital infection.

Bmpali, Bmb1 and Bmcap are necessary for uric acid granule formation in Bombyx mori

Uric acid is the end-product of nitrogen metabolism of the silkworm and other lepidopterans. The accumulation of uric acid particles in the epidermis causes the larval silkworm to appear white and opaque. However, the mechanism of uric acid granule formation is still unclear. Silkworm epidermis color is linked to the genes responsible for uric acid particle formation. We first identified two genes in the Bombyx mori genome that encode subunits of the Bloc-1 (Biogenesis of Lysosome-related Organelles Complex-1) by homology to these genes in other eukaryotes, Bmpali and Bmb1. Mutation in these genes caused a transparent phenotype in the silkworm larvae, and the loss of BmBloc-1 subunit gene Bmcap resulted in the same phenotype. These three genes are highly conserved between human and silkworm. We discovered that Bmpali, Bmcap, and Bmb1 localize in the cytoplasm of BmN cells. Yeast two-hybrid assays demonstrated that the Bmpali physically interacts with both Bmcap and Bmb1. Investigating the roles of Bmpali, Bmb1, and Bmcap is essential for uric acid granule formation understanding in Bombyx mori. These mutants present a valuable silkworm model for studying the biogenesis of lysosome-related organelles (LROs).

Genetic and transcriptome analysis of the smb mutant in Bombyx mori

More than 600 mutations have been discovered in the history of silkworm domestication. It is important to study the formation mechanism of these mutations to further understand the life and development process of silkworms and agricultural pest control. The silkworm mutant smb was isolated from silkworm strain NCV, and transcriptome analysis was performed on the silkworm mutant. 796 differentially expressed genes (DEGs) were detected at 48 h of the second instar stage with 669 genes significantly upregulated and 127 genes significantly downregulated. During the GO enrichment analysis, it was found that the enrichment of biological processes was mainly concentrated in proteolysis, carbohydrate metabolism, aminoglycan metabolism, organic substance metabolism, protein metabolism and so on. Based on the analysis of KEGG pathways, it revealed that the pathways enriched in lysosomes, AMPK signaling, fatty acid metabolism, PPAR signaling, galactose metabolism, and protein digestion and absorption were the most significant. Through these most significantly enriched GO terms and KEGG pathways, DEGs consistent with the phenotypic characteristics of the smb mutant were identified, including small body size, slow development, and successive death after the fourth instar. These results provided experimental evidence for the potential formation mechanism of smb mutants.

Significance of the clock gene period in photoperiodism in larval development and production of diapause eggs in the silkworm Bombyx mori

Many insects living in seasonal environments sense seasonal changes from photoperiod and appropriately regulate their development and physiological activities. Genetic researches have indicated the importance of a circadian clock system in photoperiodic time-measurement for photoperiodic regulations. However, most previous studies have focused on the effects on a single photoperiodic phenotype, without elucidating whether the circadian clock is involved in the core photoperiodic mechanism or only in the production of one target phenotype, such as diapause. Here, we focused on two different phenotypes in a bivoltine Kosetsu strain of the silkworm Bombyx mori, namely, embryonic diapause and larval development, and examined their photoperiodic responses and relationship to the circadian clock gene period. Photoperiod during the larval stage clearly influenced the induction of embryonic diapause and duration of larval development in the Kosetsu strain; short-day exposure leaded to the production of diapause eggs and shortened the larval duration. Genetic knockout of period inhibited the short-day-induced embryonic diapause. Conversely, in the period-knockout silkworms, the larval duration was shortened, but the photoperiodic difference was maintained. In conclusion, our results indicate that the period gene is not causally involved in the photoperiodic response of larval development, while that is essential for the short-day-induced embryonic diapause.

High-throughput and genome-scale targeted mutagenesis using CRISPR in a nonmodel multicellular organism, Bombyx mori

Large-scale genetic mutant libraries are powerful approaches to interrogating genotype-phenotype correlations and identifying genes responsible for certain environmental stimuli, both of which are the central goal of life science study. We produced the first large-scale CRISPR-Cas9-induced library in a nonmodel multicellular organism, Bombyx mori We developed a piggyBac-delivered binary genome editing strategy, which can simultaneously meet the requirements of mixed microinjection, efficient multipurpose genetic operation, and preservation of growth-defect lines. We constructed a single-guide RNA (sgRNA) plasmid library containing 92,917 sgRNAs targeting promoters and exons of 14,645 protein-coding genes, established 1726 transgenic sgRNA lines following microinjection of 66,650 embryos, and generated 300 mutant lines with diverse phenotypic changes. Phenomic characterization of mutant lines identified a large set of genes responsible for visual phenotypic or economically valuable trait changes. Next, we performed pooled context-specific positive screens for tolerance to environmental pollutant cadmium exposure, and identified KWMTBOMO12902 as a strong candidate gene for breeding applications in sericulture industry. Collectively, our results provide a novel and versatile approach for functional B. mori genomics, as well as a powerful resource for identifying the potential of key candidate genes for improving various economic traits. This study also shows the effectiveness, practicality, and convenience of large-scale mutant libraries in other nonmodel organisms.

JH degradation pathway participates in hormonal regulation of larval development of Bombyx mori following λ-cyhalothrin exposure

λ-Cyhalothrin (λ-cyh), a widely utilized pyrethroid insecticide, poses serious threats to non-target organisms due to its persistence nature in the environment. Exposure to low concentrations of λ-cyh has been observed to result in prolonged larval development in Bombyx mori, leading to substantial financial losses in sericulture. The present study was undertaken to elucidate the underlying mechanisms for prolonged development caused by λ-cyh (LC10) exposure. The results showed that the JH Ⅲ titer was significantly increased at 24 h of λ-cyh exposure, and the JH interacting genes Methoprene-tolerant 2, Steroid Receptor Co-activator, Krüppel-homolog 1, and JH binding proteins were also up-regulated. Although the target of rapamycin (Tor) genes were induced by λ-cyh, the biosynthesis of JH in the corpora allata was not promoted. Notably, 13 JH degradation genes were found to be significantly down-regulated in the midgut of B. mori. The mRNA levels and enzyme activity assays indicated that λ-cyh had inhibitory effects on JH esterase, JH epoxide hydrolase, and JH diol kinase (JHDK). Furthermore, the suppression of JHDK (KWMTBOMO01580) was further confirmed by both western blot and immunohistochemistry. This study has offered a comprehensive perspective on the mechanisms underlying the prolonged development caused by insecticides, and our results also hold significant implications for the safe production of sericulture.

Super silkworm cocoons constructed by multi-silkworm larvae: Promising composites with dense structures and excellent mechanical properties

A normal silkworm cocoon (NSC) with a unique nonwoven structure is usually spun by a single silkworm larva. Notably, there is a special Bombyx mori genetic resource that many (three or more) mature larvae tend to collectively spin into one cocoon, which was named "multi-silkworm cocoon" ("MSC"). However, the MSCs display loose structure and poor mechanical properties which limits their further application. In this study, a series of hybrid silkworm cocoons (HMSCs) are obtained by hybridizing "MSC" with a selected commercial silkworm strain successfully. The morphology, microstructures, and mechanical properties of cocoons constructed by one to three silkworm larvae were characterized and compared. The results indicated that about 48.3 % of silkworm larvae could create double and triple cocoons in the F1 generation of the silkworm hybrid, displaying robust fiber networks and dense structures. The mechanical characteristics of the HMSCs, including the tensile, peeling, compression, and needle penetration resistance properties, exceeded those of MSCs, showing significant application potential for high-performance bio-composites. This study provides a practical approach for obtaining silkworm cocoons with controllable structures and mechanical properties to develop and fabricate natural composite and biomimetic materials.

Functional characterization of Bmcap in uric acid metabolism in the silkworm

After a millennium of domestication, numerous silkworm mutants have emerged that exhibit transparent epidermis, which is caused by abnormally low levels of uric acid. We identified the Bombyx mori gene Bmcap (BMSK0003832) as the homolog of cappuccino, a subunit of the biogenesis of lysosome-related organelles complex-1 (BLOC-1) that has been extensively characterized in human, mouse, and insect species, by analyzing the amino acid sequences of putative purine metabolism genes. Using the clustered regularly interspaced palindromic repeats (CRISPR) / CRISPR-associated protein 9 system, we disrupted Bmcap, resulting in decreased uric acid levels in the silkworm epidermis and a translucent skin phenotype. In the Bmcap mutant, the purine metabolism, nitrogen metabolism, pyrimidine metabolism, and membrane system were altered compared to the wild type. Biogenesis of lysosome-related organelle complex genes play a role in the pigmentation and biogenesis of lysosome-related organelles (LROs) in platelets, melanocytes, and megakaryocytes. LROs exhibit unique morphologies and functions in various tissues and cells. Investigation of the Bmcap mutant will enhance our understanding of the uric acid metabolic pathway in silkworms, and this mutant offers a valuable silkworm model for LRO studies.

An efficient and safe strategy for germ cell-specific automatic excision of foreign DNA in F1 hybrid transgenic silkworms

The safety of transgenic technology is a major obstacle in the popularization and use of transgenic silkworms and their products. In sericulture, only the first filial generation (F1 ) hybrid eggs produced by cross-breeding Japanese and Chinese original strains are usually used for the large-scale breeding of silkworms, but this may result in uncontrolled transgene dispersal during the popularization and application of the F1 hybrid transgenic eggs. To address this issue, we developed a safe and efficient strategy using the GAL4/Upstream activating sequence (UAS) system, the FLP/flippase recognition target (FRT) system, and the gonad-specific expression gene promoters (RSHP1p and Nanosp) for the germ cell-specific automatic excision of foreign DNA in the F1 hybrid transgenic silkworms. We established 2 types of activator strains, R1p::GAL4-Gr and Nsp::GAL4-Gr, containing the testis-specific GAL4 gene expression cassettes driven by RSHP1p or Nanosp, respectively, and 1 type of effector strain, UAS::FLP-Rg, containing the UAS-linked FLP gene expression cassette. The FLP recombinase-mediated sperm-specific complete excision of FRT-flanked target DNA in the F1 double-transgenic silkworms resulting from the hybridization of R1p::GAL4-Gr and UAS::FLP-Rg was 100%, whereas the complete excision efficiency resulting from the hybridization of Nsp::GAL4-Gr and UAS::FLP-Rg ranged from 13.73% to 80.3%. Additionally, we identified a gene, sw11114, that is expressed in both testis and ovary of Bombyx mori, and can be used to establish novel gonad-specific expression systems in transgenic silkworms. This strategy has the potential to fundamentally solve the safety issue in the production of F1 transgenic silkworm eggs and provides an important reference for the safety of transgenic technology in other insect species.

Cyclin B3 plays pleiotropic roles in female reproductive organogenesis and early embryogenesis in the silkworm, Bombyx mori

BACKGROUND

The reproductive system plays a crucial role in insect survival, reproduction and species specificity. Understanding the molecular mechanisms underlying reproductive organogenesis contributes to improving the efficiency of sterile insect technique marked by an eco-friendly pest management strategy. Lepidoptera is one of the largest orders of insects, most of which are major pests in agriculture and forestry. Our study aimed to screen the genes responsible for reproductive organogenesis and unravel the mechanism underlying female reproductive organ defects.

RESULTS

Morphological investigation of female reproductive organs showed a defective connection between oviductus geminus and oviductus communis on the second day of pupa (P2) in Speckled mutant silkworm. RNA_Seq identified a total of 18 049 transcripts that were expressed in the P2 female internal reproductive organs without ovary in Spc/+ compared to +Spc /+Spc . Differential expression analysis identified 312 up-regulated genes and 221 down-regulated genes in Spc/+. KEGG analysis identified 44 significantly enriched pathways. The results of qRT-PCR performed on 33 genes significantly matched the outcomes of the RNA_Seq. Dysfunction of Cyclin B3 resulted in a defective connection of the oviductus communis with the ovariole, dysfunction of oogenesis, and a petite body. Moreover, homozygous recessive lethality of Cyclin B3/Cyclin B3 occurred during early embryogenesis.

CONCLUSION

Our results suggest that Cyclin B3 is a pleiotropic functional gene that regulates early embryogenesis, oogenesis, development, and female reproductive organogenesis. These results showed that Cyclin B3 has significant effects on lepidopteran mortality, growth, and reproductive physiology, which might be considered a novel and potentially eco-friendly target for lepidopteran pest management. © 2023 Society of Chemical Industry.

Variation in Bombyx mori immune response against fungal pathogen Beauveria bassiana with variability in cell wall β-1,3-glucan

Entomopathogenic fungi are protected by a cell wall with dynamic structure for adapting to various environmental conditions. β-1,3-Glucan recognition proteins activate the innate immune system of insects by recognizing surface molecules of fungi. However, the associations between pathogenicity and the different components of entomopathogenic fungal cell walls remain unclear. Three Beauveria bassiana strains were selected that have significantly differing virulence against Bombyx mori. The molecular mechanisms underlying the immune response in B. mori were investigated using RNA sequencing, which revealed differences in the immune response to different B. bassiana strains at 12 h post-infection. Immunofluorescence assays revealed that β-1,3-glucan content had an opposite trend to that of fungal virulence. β-1,3-Glucan injection upregulated BmβGRP4 expression and significantly reduced the virulence of the high-virulence strain but not that of the medium-virulence or low-virulence strains. BmβGRP4 silencing in B. mori with RNA interference resulted in the opposite virulence pattern, indicating that the virulence of B. bassiana was affected by the cell walls' content of β-1,3-glucan, which could be recognized by BmβGRP4. Furthermore, interference with the gene CnA (calcineurin catalytic A subunit) involved in β-1,3-glucan synthesis eliminated differences in virulence between B. bassiana strains. These results indicate that strains of a single species of pathogenic fungi that have differing cell wall components are recognized differently by the innate immune system of B. mori.

Functional annotation of insecta transcriptomes: A cautionary tale from Lepidoptera

Functional annotation is a critical step in the analysis of genomic data, as it provides insight into the function of individual genes and the pathways in which they participate. Currently, there is no consensus on the best computational approach for assigning functional annotation. This study compares three functional annotation methods (BLAST, eggNOG-Mapper, and InterProScan) in their ability to assign Gene Ontology terms in two species of Insecta with differing levels of annotation, Bombyx mori and Manduca sexta. The methods were compared for their annotation coverage, number of term assignments, term agreement and non-overlapping terms. Here we show that there are large discrepancies in gene ontology term assignment among the three computational methods, which could lead to confounding interpretations of data and non-comparable results. This study provide insight into the strengths and weaknesses of each computational method and highlight the need for more standardized methods of functional annotation.

Expression of the Wolbachia male-killing factor Oscar impairs dosage compensation in lepidopteran embryos

Wolbachia are intracellular bacteria in insects that can manipulate the sexual development and reproduction by male killing or other methods. We have recently identified a Wolbachia protein named Oscar that acts as a male-killing factor for lepidopteran insects. Oscar interacts with the Masculinizer (Masc) protein, which is required for both masculinization and dosage compensation (DC) in lepidopteran insects. Embryonic expression of Oscar inhibits masculinization and causes male killing in two lepidopteran species, Ostrinia furnacalis and Bombyx mori. However, it remains unknown whether Oscar-induced male killing is caused by a failure of DC. Here, we performed a transcriptome analysis of Oscar complementary RNA-injected O. furnacalis and B. mori embryos, and found that Oscar primarily targets the Masc protein, resulting in male killing by interfering with DC in lepidopteran insects.

Exploring the apoptotic effects of sericin on HCT116 cells through comprehensive nanostring transcriptomics and proteomics analysis

Sericin, a silk protein from Bombyx mori (silkworms), has many applications, including cosmetics, anti-inflammation, and anti-cancer. Sericin complexes with nanoparticles have shown promise for breast cancer cell lines. Apoptosis, a programmed cell death mechanism, stops cancer cell growth. This study found that Sericin urea extract significantly affected HCT116 cell viability (IC50 = 42.00 ± 0.002 µg/mL) and caused apoptosis in over 80% of treated cells. S-FTIR analysis showed significant changes in Sericin-treated cells' macromolecule composition, particularly in the lipid and nucleic acid areas, indicating major cellular modifications. A transcriptomics study found upregulation of the apoptotic signaling genes FASLG, TNFSF10, CASP3, CASP7, CASP8, and CASP10. Early apoptotic proteins also showed that BAD, AKT, CASP9, p53, and CASP8 were significantly upregulated. A proteomics study illuminated Sericin-treated cells' altered protein patterns. Our results show that Sericin activated the extrinsic apoptosis pathway via the caspase cascade (CASP8/10 and CASP3/7) and the death receptor pathway, involving TNFSF10 or FASLG, in HCT116 cells. Upregulation of p53 increases CASP8, which activates CASP3 and causes HCT116 cell death. This multi-omics study illuminates the molecular mechanisms of Sericin-induced apoptosis, sheds light on its potential cancer treatment applications, and helps us understand the complex relationship between silk-derived proteins and cellular processes.

Structural Insights into the Substrate Binding of Farnesyl Diphosphate Synthase FPPS1 from Silkworm, Bombyx mori

Farnesyl diphosphate synthase (FPPS) is an important enzyme involved in the juvenile hormone (JH) biosynthesis pathway. Herein, we report the crystal structure of a type-I Lepidopteran FPPS from Bombyx mori (BmFPPS1) at 2.80 Å resolution. BmFPPS1 adopts an α-helix structure with a deep cavity at the center of the overall structure. Computational simulations combined with biochemical analysis allowed us to define the binding mode of BmFPPS1 to its substrates. Structural comparison revealed that BmFPPS1 adopts a structural pattern similar to that of type-II FPPS but exhibits a distinct substrate-binding site. These findings provide a structural basis for understanding substrate preferences and designing FPPS inhibitors. Furthermore, the expression profiles and RNA interference of BmFPPSs indicated that they play critical roles in the JH biosynthesis and larval-pupal metamorphosis. These findings enhance our understanding of the structural features of type-I Lepidopteran FPPS while providing direct evidence for the physiological role of BmFPPSs in silkworm development.

Modeling the 3-dimensional structure of the silkworm's spinning apparatus in silk production

The silk-spinning process of the silkworms transforms the liquid silk solution to a solid state under mild conditions, making it an attractive model for bioinspiration However, the precise mechanism behind silk expulsion remains largely unknown. Here we selected the silkworms as representative models to investigate the silk-spinning mechanism. We used serial block-face scanning electron microscopy (SBF-SEM) to reconstruct the three-dimensional structures of the spinnerets in silkworms at various stages and with different gene backgrounds. By comparing the musculature and duct deformation of these spinneret models during the spinning process, we were able to simulate the morphological changes of the spinneret. Based on the results, we proposed three essential factors for silkworm spinning: the pressure generated by the silk gland, the opening duct, and the pulling force generated by head movement. Understanding the silkworm spinning process provides insights into clarify the fluid-ejecting mechanism of a group of animals. Moreover, these findings are helpful to the development of biomimetic spinning device that mimics the push-and-pull dual-force system in silkworms. STATEMENT OF SIGNIFICANCE: The silkworms' spinning system produces fibers under mild conditions, making it an ideal candidate for bioinspiration. However, the mechanism of silk expulsion is unknown, and the three-dimensional structure of the spinneret is still uncertain. In this study, we reconstructed a detailed 3-dimensional model of the spinneret at near-nanometer resolution, and for the first time, we observed the changes that occur before and during the silk-spinning process. Our reconstructed models suggested that silkworms have the ability to control the spinning process by opening or closing the spinning duct. During the continuously spinning period, both the pressure generated by the silk gland and the pulling force resulting from head movement work in tandem to expel the silk solution. We believe that gaining a full understanding of the spinning process steps can advance our ability to spin synthetic fibers with properties comparable to those of native fibers by mimicking the natural spinning process.

SilkMeta: a comprehensive platform for sharing and exploiting pan-genomic and multi-omic silkworm data

The silkworm Bombyx mori is a domesticated insect that serves as an animal model for research and agriculture. The silkworm super-pan-genome dataset, which we published last year, is a unique resource for the study of global genomic diversity and phenotype-genotype association. Here we present SilkMeta (http://silkmeta.org.cn), a comprehensive database covering the available silkworm pan-genome and multi-omics data. The database contains 1082 short-read genomes, 546 long-read assembled genomes, 1168 transcriptomes, 294 phenotype characterizations (phenome), tens of millions of variations (variome), 7253 long non-coding RNAs (lncRNAs), 18 717 full length transcripts and a set of population statistics. We have compiled publications on functional genomics research and genetic stock deciphering (mutant map). A range of bioinformatics tools is also provided for data visualization and retrieval. The large batch of omics data and tools were integrated in twelve functional modules that provide useful strategies and data for comparative and functional genomics research. The interactive bioinformatics platform SilkMeta will benefit not only the silkworm but also the insect biology communities.

Rapid production of chimeric silkworm/spider silk with improved mechanical properties by infection of nonpermissive Bombyx mori with recombinant AcMNPV harboring native-size of spidroin genes

Spider silks with excellent mechanical properties attract more attention from scientists worldwide, and the dragline silk that serves as the framework of the spider's web is considered one of the strongest fibers. However, it is unfeasible for large-scale production of spider silk due to its highly territorial, cannibalistic, predatory, and solitary behavior. Herein, to alleviate some of these problems and explore aneasy way to produce spider fibers, we constructed recombinant baculovirus Autographa californica multiple nucleopolyhedrovirus (AcMNPV) simultaneously expressing Trichonephila clavipes native ampullate spidroin 2 (MaSp-G) and spidroin 1 (MaSp-C) driven by the promoters of silkworm fibroin genes, to infect the nonpermissive Bombyx mori larvae at the fifth instar. MaSp-G and MaSp-C were co-expressed in the posterior silk glands (PSGs) of infected silkworms and successfully secreted into the lumen of the silk gland for fibroin globule assembly. The integration of MaSp-G and MaSp-C into silkworm silk fibers significantly improved the mechanical properties of these chimeric silk fibers, especially the strength and extensibility, which may be caused by the increment of β-sheet in the chimeric silkworm/spider silk fiber. These results demonstrated that silkworms could be developed as the nonpermissive heterologous host for the mass production of chimeric silkworm/spider silk fibers via the recombinant baculovirus AcMNPV.

A major endogenous glycoside hydrolase mediating quercetin uptake in Bombyx mori

Quercetin is a common plant flavonoid which is involved in herbivore-plant interactions. Mulberry silkworms (domestic silkworm, Bombyx mori, and wild silkworm, Bombyx mandarina) take up quercetin from mulberry leaves and accumulate the metabolites in the cocoon, thereby improving its protective properties. Here we identified a glycoside hydrolase, named glycoside hydrolase family 1 group G 5 (GH1G5), which is expressed in the midgut and is involved in quercetin metabolism in the domestic silkworm. Our results suggest that this enzyme mediates quercetin uptake by deglycosylating the three primary quercetin glycosides present in mulberry leaf: rutin, quercetin-3-O-malonylglucoside, and quercetin-3-O-glucoside. Despite being located in an unstable genomic region that has undergone frequent structural changes in the evolution of Lepidoptera, GH1G5 has retained its hydrolytic activity, suggesting quercetin uptake has adaptive significance for mulberry silkworms. GH1G5 is also important in breeding: defective mutations which result in discoloration of the cocoon and increased silk yield are homozygously conserved in 27 of the 32 Japanese white-cocoon domestic silkworm strains and 12 of the 30 Chinese ones we investigated.

K64 acetylation of heat shock protein 90 suppresses nucleopolyhedrovirus replication in Bombyx mori

HSP90 is a highly conserved chaperone that facilitates the proliferation of many viruses, including silkworm (bombyx mori) nucleopolyhedrovirus (BmNPV), but the underlying regulatory mechanism was unclear. We found that suppression of HSP90 by 17-AAG, a HSP90-specific inhibitor, significantly reduced the expression of BmNPV capsid protein gp64 and viral genome replication, whereas overexpression of B. mori HSP90(BmHSP90) promoted BmNPV replication. Furthermore, in a recent study of the lysine acetylome of B. mori infected with BmNPV, we focused on the reduced viral proliferation due to changes of BmHSP90 lysine acetylation. Site-directed introduction of acetylated (K/Q) or deacetylated (K/R) mimic mutations into BmHSP90 revealed that lysine 64 (K64) acetylation activated the JAK/STAT pathway and reduced BmHSP90 ATPase activity, leading to diminished chaperone activity and ultimately inhibiting BmNPV proliferation. In this study, a single lysine 64 acetylation change of BmHSP90 was elucidated as a model of posttranslational modifications occurring in the wake of host-virus interactions, providing novel insights into potential antiviral strategies.

Single-nucleus sequencing of silkworm larval midgut reveals the immune escape strategy of BmNPV in the midgut during the late stage of infection

The midgut is an important barrier against microorganism invasion and proliferation, yet is the first tissue encountered when a baculovirus naturally invades the host. However, only limited knowledge is available how different midgut cell types contribute to the immune response and the clearance or promotion of viral infection. Here, single-nucleus RNA sequencing (snRNA seq) was employed to analyze the responses of various cell subpopulations in the silkworm larval midgut to B. mori nucleopolyhedrovirus (BmNPV) infection. We identified 22 distinct clusters representing enteroendocrine cells (EEs), enterocytes (ECs), intestinal stem cells (ISCs), Goblet cell-like and muscle cell types in the BmNPV-infected and uninfected silkworm larvae midgut at 72 h post infection. Further, our results revealed that the strategies for immune escape of BmNPV in the midgut at the late stage of infection include (1) inhibiting the response of antiviral pathways; (2) inhibiting the expression of antiviral host factors; (3) stimulating expression levels of genes promoting BmNPV replication. These findings suggest that the midgut, as the first line of defense against the invasion of the baculovirus, has dual characteristics of "resistance" and "tolerance". Our single-cell dataset reveals the diversity of silkworm larval midgut cells, and the transcriptome analysis provides insights into the interaction between host and virus infection at the single-cell level.

Proteomics Analysis to Explore the Resistance Genes of Silkworm to Bombyx mori Nuclear Polyhedrosis Virus

The resistance of silkworms to Bombyx mori nuclear polyhedrosis virus (BmNPV) is controlled by a major dominant gene and multiple modifying genes. Given the presence of modified genes, it is difficult to determine the main gene by positional cloning. In this study, the main anti-BmNPV gene of BmNPV-resistant silkworm variety N was introduced into the susceptible variety Su to breed the near-isogenic line SuN with BmNPV resistance. The infection process of BmNPV in the hemolymph of Su and SuN was analyzed using the cell analysis system TissueFAXS PLUS. According to the law of infection and proliferation, hemolymph was extracted every 6 h for two-dimensional electrophoresis (2-DE) analysis and quantitative real-time polymerase chain reaction (qRT-PCR). Seven DEPs were found in comparisons between Su and SuN by 2-DE analysis. Among them, acid phosphatase, storage protein, and phenoloxidase can prevent pathogen invasion, which may play a role against BmNPV. Polyamine oxidase plays an important role in energy metabolism, which may be indirectly involved in the process of resisting BmNPV. Most of the transcriptional expression profiles of the seven DEPs were consistent with the 2-DE results. This study can provide a reference for the identification of anti-BmNPV genes and the breeding of BmNPV-resistant silkworm varieties.

The Bmtret1 Gene Family and Its Potential Role in Response to BmNPV Stress in Bombyx mori

Trehalose is a non-reducing disaccharide and participates in physiological activities such as organ formation, energy metabolism, and stress resistance in insects. The Bmtret1 gene family is mainly involved in in the sugar metabolism of silkworm. In the present study, phylogenetic analysis divided 21 Bmtret1 orthologs into three clades. These genes are equally distributed on the nine chromosomes. The cis-elements in the promoter regions of Bmtret1s indicated the possible function of Bmtret1s in response to hormones and environmental stimulus. The qPCR analysis showed the significantly different expression levels of Bmtret1s in different tissues and organs, indicating possible functional divergence. In addition, most Bmtret1s showed disturbed expression levels in response to silkworm nuclear polyhedrosis virus (BmNPV) stresses. Our results provide a clue for further functional dissection of the Tret1s in Bombyx mori and implicate them as potential regulators of antiviral responses.

[Genome-wide identification of the BmAKR gene family in the silkworm (Bombyx mori) and their expression analysis in diapause eggs and nondiapause eggs]

The aldo-keto reductase super family (AKRs) has a wide range of substrate specificity. However, the systematic identification of insect AKR gene family members has not been reported. In this study, bioinformatics methods were used to predict the phylogenetic evolution, physical and chemical properties, chromosome location, conserved motifs, and gene structure of AKR family members in Bombyx mori (BmAKR). Transcriptome data or quantitative real time polymerase chain reaction (qRT-PCR) were used to analyze the expression level of BmAKR genes during different organizational periods and silkworm eggs in different developmental states. Moreover, Western blotting was used to detect the expression level of the BmAKR in silkworm eggs. The results showed that 11 BmAKR genes were identified. These genes were distributed on 4 chromosomes of the silkworm genome, all of which had the (α/β) 8-barrel motif, and their physical and chemical characteristics were relatively similar. Phylogenetic analysis showed that the BmAKR genes could be divided into 2 subgroups (AKR1 and AKR2). Transcriptome data analysis showed that the expression of BmAKR genes were quite different in different tissues and periods. Moreover, the expression analysis of BmAKR genes in silkworm eggs showed that some genes were expressed significantly higher in nondiapause eggs than in diapause eggs; but another gene, BmAKR1-1, was expressed significantly higher in diapause eggs than in nondiapause eggs. The detection of protein level found that the difference trend of BmAKR1-1 in diapause eggs and non-diapause eggs was consistent with the results of qRT-PCR. In conclusion, BmAKR1-1 was screened out as candidates through the identification and analysis of the BmAKR genes in silkworm, which may regulate silkworm egg development is worthy of further investigation.

[Characterization and immunofluorescence localization analysis of carboxypeptidase A in molt fluid of silkworm]

Molting is an important physiological phenomenon of many metamorphosis insects, during which the old and new epidermis are separated by enzymes present in the molting fluid. Various proteomic studies have discovered the presence of Bombyx mori carboxypeptidase A (Bm-CPA) in the molting fluid of silkworm, but its function remains unclear. In order to better understand the role of Bm-CPA in the molting process of silkworm, Bm-CPA was analyzed by bioinformatics analysis, real-time fluorescence quantitative PCR, antibody preparation, immunofluorescence staining, and expression in Pichia pastoris. The results showed that Bm-CPA had a conserved M14 zinc carboxypeptidase domain and glycosylation site. Its expression was regulated by ecdysone 20E, and large expression was observed in the epidermis of the upper cluster stage. Immunofluorescence staining showed that Bm-CPA was enriched in the epidermis during the molting stage, and the inhibitor of Bm-CPA led to the larval death due to the inability to molt. We also successfully obtained a large number of recombinant Bm-CPA proteins by Pichia pastoris expression in vitro. These results may facilitate further understanding the molting development process of silkworm.

[Analyzing the evolution of insect TMED gene and the expression pattern of silkworm TMED gene]

Transmembrane emp24 domain (TMED) gene is closely related to immune response, signal transduction, growth and disease development in mammals. However, only the Drosophila TMED gene has been reported on insects. We identified the TMED family genes of silkworm, Tribolium castaneum, tobacco moth and Italian bee from their genomes, and found that the TMED family gene composition patterns of one α-class, one β-class, one δ-class and several γ-classes arose in the common ancestor of pre-divergent Hymenoptera insects, while the composition of Drosophila TMED family members has evolved in a unique pattern. Insect TMED family γ-class genes have evolved rapidly, diverging into three separate subclasses, TMED6-like, TMED5-like and TMED3-like. The TMED5-like gene was lost in Hymenoptera, duplicated in the ancestors of Lepidoptera and duplicated in Drosophila. Insect TMED protein not only has typical structural characteristics of TMED, but also has obvious signal peptide. There are seven TMED genes in silkworm, distributed in six chromosomes. One of seven is single exon and others are multi-exons. The complete open reading frame (ORF) sequences of seven TMED genes of silkworm were cloned from larval tissues and registered in GenBank database. BmTMED1, BmTMED2 and BmTMED6 were expressed in all stages and tissues of the silkworm, and all genes were expressed in the 4th and 5th instar and silk gland of the silkworm. The present study revealed the composition pattern of TMED family members, their γ class differentiation and their evolutionary history, providing a basis for further studies on TMED genes in silkworm and other insects.

BmRRS1 Protein Inhibits the Proliferation of Baculovirus Autographa californica Nucleopolyhedrovirus in Silkworm, Bombyx mori

The study of functional genes involved in baculovirus infection is vital for its wide application in pest biocontrol. This study utilized the Autographa californica nucleopolyhedrovirus (AcMNPV) and silkworm as models to elucidate the role of BmRRS1, which has been found to exhibit notable differential expression between resistant and susceptible silkworm strains. The results showed that it was evolutionarily conserved in selected species. Among different tissues, it was expressed at the highest level in the gonads, followed by the hemolymph and silk glands; among the different developmental stages, it was the highest in the second instar, followed by the pupae and adults. Moreover, its vital role in suppressing AcMNPV infection was verified by the decreased expression of lef3 and vp39 protein after overexpression of BmRRS1 as well as by the increased expression of the viral gene lef3 and the viral protein vp39 after siRNA treatment against BmRRS1 expression in BmN cells. Additionally, the direct interaction between BmRRS1 and AcMNPV was detected by the GST pull-down assay. Finally, the homologue of BmRRS1 in Spodoptera frugiperda was found to be involved in larval resistance to AcMNPV. In a word, BmRRS1 plays a vital role in AcMNPV resistance in silkworms, and this might be related to the direct interaction with AcMNPV. The results of this study provide a potential target for protecting silkworm larvae from virus infection and controlling agricultural and forestry pests.

Cytotoxicity of AuCu-Cu2S Nanocomposites: Implications for Biological Evaluation of the Nanocomposite Effect on Bombyx mori Silkworms and Cell Lines

AuCu-Cu2S nanocomposites are unique materials with exceptional properties that have recently received a lot of interest. However, little is known about their potential toxicity in terrestrial organisms and their subsequent effects on the environment. Therefore, it is essential to develop effective methodologies for evaluating AuCu-Cu2S nanocomposites in biological systems. This study reports the biological evaluation of the AuCu-Cu2S nanocomposite from animal and cell entity levels. The Bombyx mori silkworm was used as a model organism to study the effects of different concentrations of AuCu-Cu2S on silkworm development. Transcriptome analysis was also carried out to examine the genetic modulation exerted by the treatment. Moreover, biocompatibility and cytotoxicity of AuCu-Cu2S were evaluated in human bronchial epithelial cells 16HBE, human lung adenocarcinoma, and the insect Spodoptera frugiperda cell sf9 cell lines. The results showed that although AuCu-Cu2S at ≤400 ppm can prolong the eating habit of silkworms and promote the weight of the cocoon layer, there was an increase in silkworm mortality and a decrease in moth formation at a concentration of ≥800 ppm. The genetic regulation by AuCu-Cu2S treatment showed varying effects in the silkworm, primarily related to functions such as transport and catabolism, metabolism of cofactors and vitamins, xenobiotic biodegradation, amino acid, and carbohydrate. 16HBE, PC-9, and sf9 treated with 300 ppm of AuCu-Cu2S showed viability percentages of 60, 20, and 90%, respectively. Thus, AuCu-Cu2S at low concentrations serves as a safe and biocompatible material for the sf9 cell lines but is lethal to 16HBE and PC-9. This research could aid in understanding the biological effects and biocompatibility of AuCu-Cu2S nanocomposites, particularly in the field of biochemistry; however, the mechanisms involved need further exploration.

An extended Tudor domain within Vreteno interconnects Gtsf1L and Ago3 for piRNA biogenesis in Bombyx mori

Piwi-interacting RNAs (piRNAs) direct PIWI proteins to transposons to silence them, thereby preserving genome integrity and fertility. The piRNA population can be expanded in the ping-pong amplification loop. Within this process, piRNA-associated PIWI proteins (piRISC) enter a membraneless organelle called nuage to cleave their target RNA, which is stimulated by Gtsf proteins. The resulting cleavage product gets loaded into an empty PIWI protein to form a new piRISC complex. However, for piRNA amplification to occur, the new RNA substrates, Gtsf-piRISC, and empty PIWI proteins have to be in physical proximity. In this study, we show that in silkworm cells, the Gtsf1 homolog BmGtsf1L binds to piRNA-loaded BmAgo3 and localizes to granules positive for BmAgo3 and BmVreteno. Biochemical assays further revealed that conserved residues within the unstructured tail of BmGtsf1L directly interact with BmVreteno. Using a combination of AlphaFold modeling, atomistic molecular dynamics simulations, and in vitro assays, we identified a novel binding interface on the BmVreteno-eTudor domain, which is required for BmGtsf1L binding. Our study reveals that a single eTudor domain within BmVreteno provides two binding interfaces and thereby interconnects piRNA-loaded BmAgo3 and BmGtsf1L.

Cold-stored mulberry leaves affect antioxidant system and silk proteins of silkworm (Bombyx mori) larva

BACKGROUND

Cold storage has been widely used to maintain the quality of vegetables, but whether eating cold-stored vegetables affects health remains unknown.

RESULTS

This study used silkworms as an animal model to evaluate the effects of nutrient changes in cold-stored mulberry leaves (CSML) on health. Compared with fresh mulberry leaves (FML), CSML contained lower vitamin C, soluble sugars and proteins, and higher H2 O2 , suggesting decreased antioxidant ability and nutrition. The CSML did not obviously affect larval survival rate, body weight or dry matter rate, cocoon shape, weight and size, or final rates of cluster and cocooning relative to the FML, suggesting CSML did not alter overall growth and development. However, the CSML increased the initial rates of cluster and cocooning and upregulated BmRpd3, suggesting CSML shortened larval lifespan and enhanced senescence. CSML upregulated BmNOX4, downregulated BmCAT, BmSOD and BmGSH-Px and increased H2 O2 in silkworms, suggesting CSML caused oxidative stress. CSML upregulated ecdysone biosynthesis and inactivation genes and elevated ecdysone concentration in silkworms, suggesting that CSML affected hormone homeostasis. CSML upregulated apoptosis-related genes, downregulated sericin and silk fibroin genes and decreased sericin content rate in silkworms, suggesting oxidative stress and protein deficiency.

CONCLUSION

Cold storage reduced nutrition and antioxidant capability of mulberry leaves. CSML did not influence growth and development of silkworm larva, but affected health by causing oxidative stress and reducing protein synthesis. The findings show that the ingredient changes in CSML had negative effects on health of silkworms. © 2023 Society of Chemical Industry.

Identification and characterization of extrachromosomal circular DNA in the silk gland of Bombyx mori

The silk gland cells of silkworm are special cells which only replicate DNA in the nucleus without cell division throughout the larval stage. The extrachromosomal circular DNAs (eccDNAs) have not yet been reported in the silk gland of silkworms. Herein, we have explored the characterization of eccDNAs in the posterior silk gland of silkworms. A total of 35 346 eccDNAs were identified with sizes ranging from 30 to 13 569 549 bp. Motif analysis revealed that dual direct repeats are flanking the 5' and 3' breaking points of eccDNA. The sequences exceeding 1 kb length in eccDNAs present palindromic sequence characteristics flanking the 5' and 3' breaking points of the eccDNA. These motifs might support possible models for eccDNA generation. Genomic annotation of the eccDNA population revealed that most eccDNAs (58.6%) were derived from intergenic regions, whereas full or partial genes were carried by 41.4% of eccDNAs. It was found that silk protein genes fib-H, fib-L, and P25, as well as the transcription factors SGF and sage, which play an important regulatory role in silk protein synthesis, could be carried by eccDNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses showed that the genes carried by eccDNAs were mainly associated with the development and metabolism-related signaling pathways. Moreover, it was found that eccDNAfib-L could promote the transcription of fib-L gene. Overall, the results of the present study not only provide a novel perspective on the mechanism of silk gland development and silk protein synthesis but also complement previously reported genome-scale eccDNA data supporting that eccDNAs are common in eukaryotes.

A secreted phospholipase A2 (BmsPLA2 ) regulates melanization of immunity through BmDDC in the silkworm Bombyx mori

Insect immune-associated phospholipase A2 (PLA2 ) is an important target of pathogen invasion. Melanization, an effective defense response, has significant correlations with other immune responses to coordinate immune attack against invaders. However, the effect of PLA2 on melanization has not yet been reported in insects or other arthropods. In this work, we cloned a PLA2 gene (BmsPLA2 ), and its protein had characteristic features of secreted PLA2 (sPLA2 ). After injection of bacteria, BmsPLA2 expression and sPLA2 activity in hemolymph significantly increased. BmsPLA2 fluorescence was transferred from the cytoplasm to the cell membranes of circulating hemocytes. These results indicated that BmsPLA2 was related to hemolymph immunity in silkworms. Interestingly, reducing BmsPLA2 by RNA interference decreased melanosis (melanistic hemocytes) levels in vivo and in vitro, while BmsPLA2 overexpression had the opposite effect. The larval survival and melanization rate in the hemocoel both slowed depending on the PLA2 inhibitor dosage. These results demonstrated that BmsPLA2 plays a role in melanization during the immune process of silkworms. Surprisingly, the level of BmDDC matched the degree of melanization in various observations. BmDDC expression showed a significant increase, with the peak occurring later than that of BmsPLA2 after injection of bacteria, implying that BmsPLA2 was activated prior to BmDDC. Moreover, the alteration of BmsPLA2 by RNA interference or overexpression led to altered BmDDC levels. These results suggested that BmsPLA2 regulates the melanization response in silkworms through BmDDC. Our study proposes a new regulatory mechanism of the melanization response and new directions for understanding the complex immune networks of insects.

ABC transporter subfamily B1 as a susceptibility determinant of Bombyx mori larvae to Cry1Ba, Cry1Ia and Cry9Da toxins

ATP binding cassette (ABC) transporters are a diverse family of transmembrane proteins. Specific subfamily members expressed in the lepidopteran midgut can act as susceptibility determinants for several insecticidal Bt Cry proteins. However, the susceptibility determinants to many Cry toxins still remain unclear. Therefore, we knocked out a series of ABC transporters that are highly expressed in the midgut of Bombyx mori larvae by transcription activator-like effector nuclease (TALEN)-mediated gene editing, and the lineages that became resistant to Cry toxins were searched by toxin overlay bioassay. As a result, the B. mori ABC transporter subfamily B1 (BmABCB1) knockout lineage showed 19.17-fold resistance to Cry1Ba, 876.2-fold resistance to Cry1Ia, and 29.1-fold resistance to Cry9Da, suggesting that BmABCB1 is the determinant of susceptibility to these toxins. BmABCC2 and BmABCC3 have been shown to be susceptibility determinants based on their function as receptors. Therefore, we next heterologously expressed these ABC transporters in HEK293T cells and performed a cell swelling assay to examine whether these molecules could exert receptor functions. As a result, BmABCB1-expressing cells showed swelling response to Cry1Ia and Cry9Da, and cells expressing PxABCB1, which is the Plutella xylostella ortholog of BmABCB1, showed swelling for Cry1Ba, suggesting that ABCB1 is a susceptibility determinant by functioning as a receptor to these toxins. Furthermore, in order to clarify how high binding affinity is based on receptor function, we performed surface plasmon resonance analysis and found that each KD of Cry1Ba, Cry1Ia, and Cry9Da to BmABCB1 were 7.69 × 10-8 M, 2.19 × 10-9 M, and 4.17 × 10-6 M respectively.

Inhibitory Effect on RT-PCR and Restriction Enzyme Activity by Ommochrome and Its Mechanism

To explore the physiological role and/or pharmacological effects of ommochrome, which is a natural organic pigment widely distributed in Protostomia, we attempted to investigate the influence of ommochrome on RT-PCR and activities of restriction enzymes. It was found that ommin, an ommochrome purified from the diapause eggs of Bombyx mori, inhibited the RT-PCR and restriction enzyme activities. The mechanism of these inhibitory reactions is assumed to be the direct binding of ommochrome to DNA rather than acting against the enzymes because, similarly to actinomycin D, there is a phenoxazine ring in the structure of ommin that is known to be intercalated to DNA. To reveal the ommin/DNA interaction, it was investigated by computational approaches such as molecular docking, molecular dynamics simulation, and free energy calculation. From the computational analyses, it was expected that ommin would bind to DNA with almost the same strength as actinomycin D and intercalate into DNA. This is the first report on the pharmacological effect of ommochrome and its inhibitory mechanism obtained from biochemical and computational analyses.

BmPxt1 mediated immune response by regulating PGE2 in silkworm, Bombyx mori

Prostaglandins (PGs) mediates the immune response of insects to multiple stimuli. Mammalian cyclooxygenase (COXs) is a key enzyme in the synthesis of PGs, and peroxinectin (Pxt) may have similar functions in some sequenced insect genomes. As a representative of Lepidoptera, the silkworm also contains PGs, but its synthetic pathway is not clear. We cloned a full-length cDNA encoding a Pxt, designated as BmPxt1, from silkworm. Sequence alignment analysis showed that the protein encoded by BmPxt1 has a conserved domain similar to Pxts, and its catalytic site is shared with the Pxt of Manduca sexta, which also produces PGs. The expression of BmPxt1 gene was the highest in the hemocytes and was induced by Nuclear Polyhedrosis Virus (NPV) challenge in the detected tissues. Moreover, we found that dsPxt1 treatment deficiency down-regulated BmPxt1 transcript levels and efficiently inhibiting hemocyte-spreading and nodule formation in silkworm. Hemocyte-spreading, nodule formation, phenoloxidase (PO) and AMP genes (attacin, defencin and moricin) were also inhibited by aspirin, a COX inhibitor. Treatment by PGE2 but not arachidonic acid (AA) rescued the immunosuppression; PGs concentrations was also inhibited by aspirin. PGE2, but not AA, treatment rescued the PGs concentrations. The COX inhibitor, aspirin, impaired the innate immune response including nodulation, encapsulation, and melanization in silkworm, while PGE2, but not arachidonic acid (AA), partially reversed these effects of aspirin. Recombinant BmsPxt1 significantly induced PO activation in larvae hemolymph, PGs concentrations and encapsulation of agarose beads. Injection of recombinant BmsPxt1 into larvae resulted in increased transcript levels of AMP genes. Our results confirmed that BmPxt1 was involved in the synthesis of PGs in the innate immune response of silkworm larvae, and provided new information for the role of BmsPxt1 secreted by silkworm in activating PO and antimicrobial peptides.

A novel chimeric RNA originating from BmCPV S4 and Bombyx mori HDAC11 transcripts regulates virus proliferation

Polymerases encoded by segmented negative-strand RNA viruses cleave 5'-m7G-capped host transcripts to prime viral mRNA synthesis ("cap-snatching") to generate chimeric RNA, and trans-splicing occurs between viral and cellular transcripts. Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), an RNA virus belonging to Reoviridae, is a major pathogen of silkworm (B. mori). The genome of BmCPV consists of 10 segmented double-stranded RNAs (S1-S10) from which viral RNAs encoding a protein are transcribed. In this study, chimeric silkworm-BmCPV RNAs, in which the sequence derived from the silkworm transcript could fuse with both the 5' end and the 3' end of viral RNA, were identified in the midgut of BmCPV-infected silkworms by RNA_seq and further confirmed by RT-PCR and Sanger sequencing. A novel chimeric RNA, HDAC11-S4 RNA 4, derived from silkworm histone deacetylase 11 (HDAC11) and the BmCPV S4 transcript encoding viral structural protein 4 (VP4), was selected for validation by in situ hybridization and Northern blotting. Interestingly, our results indicated that HDAC11-S4 RNA 4 was generated in a BmCPV RNA-dependent RNA polymerase (RdRp)-independent manner and could be translated into a truncated BmCPV VP4 with a silkworm HDAC11-derived N-terminal extension. Moreover, it was confirmed that HDAC11-S4 RNA 4 inhibited BmCPV proliferation, decreased the level of H3K9me3 and increased the level of H3K9ac. These results indicated that during infection with BmCPV, a novel mechanism, different from that described in previous reports, allows the genesis of chimeric silkworm-BmCPV RNAs with biological functions.

Advances in the extracellular signal-regulated kinase signaling pathway in silkworms, Bombyx mori (Lepidoptera)

Signaling pathways regulate the transmission of signals during organism growth and development, promoting the smooth and accurate completion of numerous physiological and biochemical reactions. The extracellular signal-regulated kinase (ERK) signaling pathway is an essential pathway involved in regulating various physiological processes, such as cell proliferation, differentiation, adhesion, migration, and more. This pathway also contributes to several important physiological processes in silkworms, including protein synthesis, reproduction, and immune defense against pathogens. Organizing related studies on the ERK signaling pathway in silkworms can provide a better understanding of its mechanism in Lepidopterans and develop a theoretical foundation for improving cocoon production and new strategies for pest biological control.

Determination of suitable reference genes for RT-qPCR normalisation in Bombyx mori (Lepidoptera: Bombycidae) infected by the parasitoid Exorista sorbillans (Diptera, Tachinidae)

The silkworm Bombyx mori (Lepidoptera: Bombycidae) is a lepidopteran model insect of great economic importance. The parasitoid Exorista sorbillans (Diptera, Tachinidae) is the major pest of B. mori and also a promising candidate for biological control. However, the molecular interactions between hosts and dipteran parasitoids have only partially been studied. Gene expression analysis by reverse-transcription quantitative real-time polymerase chain reaction (RT-qPCR) is indispensable to characterise their interactions. Accurate normalisation of RT-qPCR-based gene expression requires the use of reference genes that are constantly expressed irrespective of experimental conditions. In this study, the expression stability of 13 traditionally used reference genes was estimated by five statistical algorithms (ΔCt, geNorm, Normfinder, BestKeeper, and RefFinder) to determine the best reference genes for gene expression studies in different tissues of B. mori under E. sorbillans parasitism. Specifically, TATA-box-binding protein was the best reference gene in epidermis and testis, while elongation factor 1α was the most stable gene in prothoracic gland and midgut. Elongation factor 1γ, ribosomal protein L3, actin A1, ribosomal protein L40, glyceraldehyde-3-phosphate dehydrogenase and eukaryotic translation initiation factor 4A were the most suitable genes in head, silk gland, fat body, haemolymph, Malpighian tubule and ovary, respectively. Our study offers a set of suitable reference genes for gene expression normalisation in B. mori under the parasitic stress of E. sorbillans, which will benefit the in-depth exploration of host-dipteran parasitoid interactions, and also provide insights for further improvements of B. mori resistance against parasitoids and biocontrol efficacy of dipteran parasitoids.

Methyl-CpG binding domain (MBD)2/3 specifically recognizes and binds to the genomic mCpG site with a β-sheet in the MBD to affect embryonic development in Bombyx mori

Methyl-CpG (mCpG) binding domain (MBD) proteins especially bind with methylated DNA, and are involved in many important biological processes; however, the binding mechanism between insect MBD2/3 and mCpG remains unclear. In this study, we identified 2 isoforms of the MBD2/3 gene in Bombyx mori, MBD2/3-S and MBD2/3-L. Binding analysis of MBD2/3-L, MBD2/3-S, and 7 mutant MBD2/3-L proteins deficient in β1-β6 or α1 in the MBD showed that β2-β3-turns in the β-sheet of the MBD are necessary for the formation of the MBD2/3-mCpG complex; furthermore, other secondary structures, namely, β4-β6 and an α-helix, play a role in stabilizing the β-sheet structure to ensure that the MBD is able to bind mCpG. In addition, sequence alignment and binding analyses of different insect MBD2/3s indicated that insect MBD2/3s have an intact and conserved MBD that binds to the mCpG of target genes. Furthermore, MBD2/3 RNA interference results showed that MBD2/3-L plays a role in regulating B. mori embryonic development, similar to that of DNA methylation; however, MBD2/3-S without β4-β6 and α-helix does not alter embryonic development. These results suggest that MBD2/3-L recognizes and binds to mCpG through the intact β-sheet structure in its MBD, thus ensuring silkworm embryonic development.

Synthesis of 1,4-benzodioxan derivatives and the evaluation of their biological activity as a novel juvenile hormone signaling inhibitor

BACKGROUND

Juvenile hormone (JH) signaling inhibitors may be used as insect growth regulators because of their ability to control metamorphosis and reproduction in insects by regulating the action of JH.

RESULTS

We identified ethyl (E)-3-(4-{[7- (4-methoxycarbonylbenzyloxy)-1,4-benzodioxan-6-yl]methyl}phenyl)prop-2-enoate (EMBP) and observed its strong precocious metamorphosis-inducing activity against silkworm larvae. To further elucidate its mechanism of action, we investigated the effect of EMBP on the JH-mediated signaling pathway in vitro and in vivo. In a reporter assay using a Bombyx mori cell line, EMBP strongly suppressed the induction of reporter gene expression by Juvenile hormone I (JH I) in a concentration-dependent manner. A parallel rightward shift was observed in the dose-response curve of JH I after treatment with EMBP, indicating that EMBP competitively inhibited JH. Moreover, we monitored developmental changes in the JH-responsive gene, Krüppel homolog 1 (Kr-h1), and ecdysone-responsive gene, Broad-Complex (BRC), in EMBP-treated silkworm larvae. EMBP suppressed only the expression of Kr-h1 in third-instar larvae.

CONCLUSION

Our results demonstrated that EMBP specifically regulates the JH-mediated Kr-h1 signaling pathway. EMBP could be used as a lead compound in the development of new insect growth regulators. © 2023 Society of Chemical Industry.

Fine mapping of Green a, Ga, on chromosome 27 in Bombyx mori

Some strains of silkworms produce green cocoons of varying intensities. This results from quantitative and qualitative differences in flavonoid pigments, which are influenced by the environment and genetic background. We discovered that the appearance of a faint green cocoon is regulated by a gene (G27) located on chromosome 27. Through mating experiments, we found that G27 is identical to an essential flavonoid cocoon gene, Ga. This locus has not been previously described. Furthermore, we narrowed down the Ga region to 438 kbp using molecular markers. Within this region, several predicted genes for sugar transporters form a cluster structure, suggesting that Ga is among them.

Effects of non-lethal Cry1F toxin exposure on the growth, immune response, and intestinal microbiota of silkworm (Bombyx mori)

Bt (Bacillus thuringiensis) maize is expected to be commercial cultivated widely in China. When Bt maize is planted near mulberry trees, it renders silkworms (Bombyx mori) vulnerable, as they belong to the same class as the Lepidoptera insects targeted by Bt maize. Cry1F has been found to be highly toxic to silkworms, particularly in their early larval stages. In this study, we aimed to assess the effects of non-lethal Cry1F exposure on the growth, immune response, and intestinal microbiota in silkworms. The results showed that feeding silkworms with mulberry leaves soaked in 100 μg/mL Cry1F for 96 h had an impact on larval body weight acquisition, leading to a decrease in cocoon and pupae weight. Cry1F exposure disrupted the intestinal integrity of silkworms by affecting the columnar cells of the midgut. The activity of detoxification enzymes (CarE, AChE, and GST) as well as antioxidant enzymes (SOD, CAT, and POD) were also affected by Cry1F. After 96 h Cry1F exposure, the evenness of the bacterial community was disrupted, resulting in alterations in the structure of the intestinal microbiota. Additionally, Cry1F exposure affected the relative expression levels of the peritrophic membrane (PM) protein and the corresponding immune pathways genes of silkworms. Most of the immune-related gene expressions were inhibited after exposure to Cry1F toxin but increased with prolonged treatment. This study demonstrates that non-lethal Cry1F exposure can affect the growth, immune response, and intestinal microbiota of silkworm.