Transgenic modifications of silkworms as a means to obtain therapeutic biomolecules and protein fibers with exceptional properties

Transgenic modification of Bombyx mori silkworms is a benign approach for the production of silk fibers with extraordinary properties and also to generate therapeutic proteins and other biomolecules for various applications. Silk fibers with fluorescence lasting more than a year, natural protein fibers with strength and toughness exceeding that of spider silk, proteins and therapeutic biomolecules with exceptional properties have been developed using transgenic technology. The transgenic modifications have been done primarily by modifying the silk sericin and fibroin genes and also the silk producing glands. Although the genetic modifications were typically performed using the sericin 1 and other genes, newer techniques such as CRISPR/Cas9 have enabled successful modifications of both the fibroin H-chain and L-chain. Such modifications have led to the production of therapeutic proteins and other biomolecules in reasonable quantities at affordable costs for tissue engineering and other medical applications. Transgenically modified silkworms also have distinct and long-lasting fluorescence useful for bioimaging applications. This review presents an overview of the transgenic techniques for modifications of B. mori silkworms and the properties obtained due to such modifications with particular focus on production of growth factors, fluorescent proteins, and high performance protein fibers.

Serpin-1a and serpin-6 regulate the Toll pathway immune homeostasis by synergistically inhibiting the Spätzle-processing enzyme CLIP2 in silkworm, Bombyx mori

The Toll receptor signaling pathway is an important innate immune response of insects to pathogen infection; its extracellular signal transduction involves serine protease cascade activation. However, excessive or constitutive activation of the Toll pathway can be detrimental. Hence, the balance between activation and inhibition of the extracellular protease cascade must be tightly regulated to achieve favorable outcomes. Previous studies have shown that serpins-serine protease inhibitors-negatively regulate insect innate immunity by inhibiting extracellular protease cascade signaling. Although the roles of serpins in insect innate immunity are well described, the physiological mechanisms underlying their synergistic effects remain poorly understand. Here, we characterize the molecular mechanism by which serpin-1a and serpin-6 synergistically maintain immune homeostasis of the silkworm Toll pathway under physiological and pathological conditions. Through in vitro biochemical assays and in vivo bioassays, we demonstrate that clip-domain serine protease 2 (CLIP2), as the Toll cascade-activating terminal protease, is responsible for processing proSpätzle1 to induce the expression of antimicrobial peptides. Further biochemical and genetic analyses indicate that constitutively expressed serpin-1a and inducible serpin-6 synergistically target CLIP2 to maintain homeostasis of the silkworm Toll pathway under physiological and pathological conditions. Taken together, this study provides new insights into the precise regulation of Toll cascade activation signals in insect innate immune responses and highlights the importance and complexity of insect immune homeostasis regulation.

Uric acid metabolism promotes apoptosis against Bombyx mori nucleopolyhedrovirus in silkworm, Bombyx mori

The white epidermis of silkworms is due to the accumulation of uric acid crystals. Abnormal silkworm uric acid metabolism decreases uric acid production, leading to a transparent or translucent phenotype. The oily silkworm op50 is a mutant strain with a highly transparent epidermis derived from the p50 strain. It shows more susceptibility to Bombyx mori nucleopolyhedrovirus (BmNPV) infection than the wild type; however, the underlying mechanism is unknown. This study analysed the changes in 34 metabolites in p50 and op50 at different times following BmNPV infection based on comparative metabolomics. The differential metabolites were mainly clustered in six metabolic pathways. Of these, the uric acid pathway was identified as critical for resistance in silkworms, as feeding with inosine significantly enhanced larval resistance compared to other metabolites and modulated other metabolic pathways. Additionally, the increased level of resistance to BmNPV in inosine-fed silkworms was associated with the regulation of apoptosis, which is mediated by the reactive oxygen species produced during uric acid synthesis. Furthermore, feeding the industrial strain Jingsong (JS) with inosine significantly increased the level of larval resistance to BmNPV, indicating its potential application in controlling the virus in sericulture. These results lay the foundation for clarifying the resistance mechanism of silkworms to BmNPV and provide new strategies and methods for the biological control of pests.

Expressions of sugar transporter genes during Bombyx mori embryonic development

Sugar transporters (Sts) play important roles in controlling carbohydrate transport and are responsible for mediating the movement of sugars into cells. Few studies have been conducted on expressions of Sts during insect embryonic development. In the present study, we investigated temporal expressions of St genes during the embryonic diapause process in Bombyx mori. We found that in HCl-treated developing eggs, high gene expressions of trehalose transporter 1 (Tret1) were detected during middle and later embryonic development. St4 and St3 gene expressions gradually increased during the early stages, reached a small peak on Day 3, and large peaks were again detected on Day 7. However, in diapause eggs, expression levels of the Tret1, St4, and St3 genes all remained at low levels. Differential temporal changes in expressions of the Tret1, St4, and St3 genes found between diapause and HCl-treated eggs were further confirmed using nondiapause eggs. Our results showed that nondiapause eggs exhibited similar changing patterns as those of HCl-treated eggs, thus clearly indicating potential correlations between expressions of these genes and embryonic development. In addition, high gene expressions of Tret1 were also detected when dechorionated eggs were incubated in the medium. The addition of LY294002 (a specific phosphatidylinositol 3-kinase [PI3K] inhibitor) and U0126 (a mitogen-activated protein kinase/extracellular signal-regulated kinase [ERK] kinase [MEK] inhibitor) partially inhibited Tret1 gene expression in dechorionated eggs, but did not affect either ecdysteroid-phosphate phosphatase gene expression or ecdysteroid biosynthesis, clearly indicating that both PI3K and ERK are involved in increased gene expression of Tret1 that was independent of ecdysteroid levels. To our knowledge, this is the first comprehensive report to demonstrate the transcriptional regulation of St genes during embryonic development, thus providing useful information for a clearer understanding of insect egg diapause mechanisms.

The acetyltransferase BmCBP changes the acetylation modification of BmSP3 and affects its protein expression in silkworm, Bombyx mori

BACKGROUND

Protein acetylation is an important post-translational modification (PTM) that widely exists in organisms. As a reversible PTM, acetylation modification can regulate the function of proteins with high efficiency. In the previous study, the acetylation sites of silkworm proteins were identified on a large scale by nano-HPLC/MS/MS (nanoscale high performance liquid chromatography-tandem secondary mass spectrometry), and a total of 11 acetylation sites were discovered on Bombyx mori nutrient-storage protein SP3 (BmSP3). The purpose of this study was to investigate the effect of acetylation level on BmSP3.

METHODS AND RESULTS

In this study, the acetylation of BmSP3 was further verified by immunoprecipitation (IP) and Western blotting. Then, it was confirmed that acetylation could up-regulate the expression of BmSP3 by improving its protein stability in BmN cells. Co-IP and RNAi experiments showed acetyltransferase BmCBP could bind to BmSP3 and catalyze its acetylation modification, then regulate the expression of BmSP3. Furthermore, the knock-down of BmCBP could improve the ubiquitination level of BmSP3. Both acetylation and ubiquitination occur on the side chain of lysine residues, therefore, we speculated that the acetylation of BmSP3 catalyzed by BmCBP could competitively inhibit its ubiquitination modification and improve its protein stability by inhibiting ubiquitin-mediated proteasome degradation pathway, and thereby increase the expression and intracellular accumulation.

CONCLUSIONS

BmCBP catalyzes the acetylation of BmSP3 and may improve the stability of BmSP3 by competitive ubiquitination. This conclusion provides a new functional basis for the extensive involvement of acetylation in the regulation of nutrient storage and utilization in silkworm, Bombyx mori.

Purification, Immunological Identification, and Characterization of the Novel Silkworm Pupae Allergen Bombyx mori Lipoprotein 3 (Bomb m 6)

Allergic reactions caused by silkworm pupae greatly limit their utilization, and studies suggest that silkworm pupae proteins of 25-30 kDa may be the principal allergens. To further understand these allergens, we attempted to purify a protein of about 30 kDa by ammonium sulfate salting, pH-graded precipitation, and ion-exchange chromatography. The protein was identified by mass spectrometry and characterized by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), western blot, enzyme-linked immunosorbent assays, circular dichroism, and fluorescence spectroscopy analyses. We identified the purified protein as Bombyx mori lipoprotein 3 (Bmlp3), which has high IgE reactivity and is a novel uncharacterized allergen that we named Bomb m 6 according to the WHO/IUIS Allergen Nomenclature Sub-Committee. This allergen is stable against heat, acids, bases, and digestion. In conclusion, we successfully purified and characterized a novel silkworm pupa allergen, which may inform the diagnosis and treatment of silkworm pupa allergies.

Low concentration of indoxacarb interferes with the growth and development of silkworm by damaging the structure of midgut cells

As an important economic insect, Bombyx mori plays an essential role in the development of the agricultural economy. Indoxacarb, a novel sodium channel blocker insecticide, has been widely used for the control of various pests in agriculture and forestry, and its environmental pollution caused by flight control operations has seriously affected the safe production of sericulture in recent years. However, the lethal toxicity and adverse effects of indoxacarb on silkworm remain largely unknown. In this study, the toxicity of indoxacarb on the 5th instar larvae of silkworm was determined, with an LC50 (72 h) of 2.07 mg/L. Short-term exposure (24 h) to a low concentration of indoxacarb (1/2 LC50) showed significantly reduced body weight and survival rate of silkworm larvae. In addition, indoxacarb also led to decreased cocoon weight and cocoon shell weight, but had no significant effects on pupation, adult eclosion, and oviposition. Histopathological and ultrastructural analysis indicated that indoxacarb could severely damage the structure of the midgut epithelial cells, and lead to physiological impairment of the midgut. A total of 3883 differentially expressed genes (DEGs) were identified by midgut transcriptome sequencing and functionally annotated using GO and KEGG. Furthermore, the transcription level and enzyme activity of the detoxification related genes were determined, and our results suggested that esterases (ESTs) might play a major role in metabolism of indoxacarb in the midgut of B. mori. Future studies to examine the detoxification or biotransformation function of candidate genes will greatly enhance our understanding of indoxacarb metabolism in B. mori. The results of this study provide a theoretical basis for elucidating the mechanism of toxic effects of indoxacarb on silkworm by interfering with the normal physiological functions of the midgut.

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.

Molecular diagnosis of 5 silkworm strains endemic to South Korea using single-nucleotide polymorphisms selected from whole-genome sequences

The domesticated silkworm, Bombyx mori Linnaeus (Lepidoptera: Bombycidae), often poses a challenge in strain identification due to similarities in morphology and genetic background. In South Korea, around 40 silkworm strains are classified as premium, including 5 endemic tri-molting strains: Goryeosammyeon, Sammyeonhonghoeback, Hansammyeon, Sun7ho, and Sandongsammyeon. These strains have potential for breeding programs in response to emerging industry demands, necessitating a reliable strain identification method. In this study, we established a molecular diagnosis approach for these 5 strains. We selected 2-4 single-nucleotide polymorphisms (SNPs) for each strain from whole-genome sequences of 39 strains, encompassing 37 previously studied and 2 newly added. These SNPs were utilized to construct decision trees for each endemic strain identification. The SNPs can be used to distinguish each target strain from the 38 nontarget strains by the tetra-primer amplification refractory mutation system-polymerase chain reaction, with the exception of HMS which needs the addition of PCR-restriction fragment length polymorphism method at the final step. This decision tree-based method using genomic SNPs, coupled with the 2 typing methods, produced consistent and accurate results, providing 100% accuracy. Additionally, the significant number of remaining SNPs identified in this study could be valuable for future diagnosis of the other strains.

Afidopyropen suppresses silkworm growth and vitality by affecting carbohydrate metabolism and immune function

Afidopyropen has strong insecticidal toxicity to sucking pests by silencing the vanilloid-type transient receptor potential (TRPV) channels. However, the toxicity of afidopyropen to the Lepidoptera model insect silkworm remain unknown. In this study, the LC50 of afidopyropen to the silkworm at 72 h exposure was 256.82 mg/L. This indicates that afidopyropen is moderately toxic to the silkworm. Long-term exposure to concentrations of 100 mg/L, or less, of afidopyropen, significantly reduced silkworm growth, vitality, silk protein synthesis, and fecundity. A total of 220 differentially expressed genes (DEGs) were detected by transcriptome sequencing, among which 166 were downregulated and 54 were upregulated. Gene Ontology (GO) enrichment analysis showed that the DEGs were enriched in the immune system, immune response and carbohydrate metabolism. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that DEGs were primarily concentrated in carbohydrate metabolism and biosynthesis of neomycin, kanamycin and gentamicin. Genes related to carbohydrate metabolism and immune system pathways in silkworm were detected by quantitative real-time PCR. The results showed that the genes related to carbohydrate metabolism, silk protein synthesis, and immune response were significantly downregulated. These genes included BCL-6 corepressor-like protein 1 (BCORL1), hexokinase type 2 (HEXO2), phosphoserine aminotransferase 1 (PSAT1), relish (Rel), peptidoglycan recognition protein 2 (PGRP2) and 27 kda glycoprotein precursor (P27K). The data demonstrated the toxic effects of afidopyropen against the silkworm and its regulation of genes responsible for immune function and abscissa carbohydrate metabolism.

Masculinizer gene controls male sex determination in the codling moth, Cydia pomonella

The molecular mechanisms of sex determination in moths and butterflies (Lepidoptera) with female heterogamety (WZ/ZZ) are poorly understood, except in the silkworm Bombyx mori. However, the Masculinizer (Masc) gene that controls male development and dosage compensation in B. mori, appears to be conserved in Lepidoptera, as its masculinizing function was recently confirmed in several moth species. In this work, we investigated the role of the Masc gene in sex determination of the codling moth Cydia pomonella (Tortricidae), a globally important pest of pome fruits and walnuts. The gene structure of the C. pomonella Masc ortholog, CpMasc, is similar to B. mori Masc. However, unlike B. mori, we identified 14 splice variants of CpMasc in the available transcriptomes. Subsequent screening for sex specificity and genetic variation using publicly available data and RT-PCR revealed three male-specific splice variants. Then qPCR analysis of these variants revealed sex-biased expression showing a peak only in early male embryos. Knockdown of CpMasc by RNAi during early embryogenesis resulted in a shift from male-to female-specific splicing of the C. pomonella doublesex (Cpdsx) gene, its downstream effector, in ZZ embryos, leading to a strongly female-biased sex ratio. These data clearly demonstrate that CpMasc functions as a masculinizing gene in the sex-determining cascade of C. pomonella. Our study also showed that CpMasc transcripts are provided maternally, as they were detected in unfertilized eggs after oviposition and in mature eggs dissected from virgin females. This finding is unique, as maternal provision of mRNA has rarely been studied in Lepidoptera.

Non-gonadal somatic piRNA pathways ensure sexual differentiation, larval growth, and wing development in silkworms

PIWI-interacting RNAs (piRNAs) guide PIWI proteins to target transposons in germline cells, thereby suppressing transposon activity to preserve genome integrity in metazoans' gonadal tissues. Piwi, one of three Drosophila PIWI proteins, is expressed in the nucleus and suppresses transposon activity by forming heterochromatin in an RNA cleavage-independent manner. Recently, Piwi was reported to control cell metabolism in Drosophila fat body, providing an example of piRNAs acting in non-gonadal somatic tissues. However, mutant flies of the other two PIWI proteins, Aubergine (Aub) and Argonaute3 (Ago3), show no apparent phenotype except for infertility, blurring the importance of the piRNA pathway in non-gonadal somatic tissues. The silkworm, Bombyx mori, possesses two PIWI proteins, Siwi (Aub homolog) and BmAgo3 (Ago3 homolog), whereas B. mori does not have a Piwi homolog. Siwi and BmAgo3 are mainly expressed in gonadal tissues and play a role in repressing transposon activity by cleaving transposon RNA in the cytoplasm. Here, we generated Siwi and BmAgo3 loss-of-function mutants of B. mori and found that they both showed delayed larval growth and failed to become adult moths. They also exhibited defects in wing development and sexual differentiation. Transcriptome analysis revealed that loss of somatic piRNA biogenesis pathways results in abnormal expression of not only transposons but also host genes, presumably causing severe growth defects. Our results highlight the roles of non-gonadal somatic piRNAs in B. mori development.

Characterization of glutathione-specific gamma glutamyl cyclotransferase (ChaC) in Bombyx mori

Glutathione (GSH) contributes to redox maintenance and detoxification of various xenobiotic and endogenous substances. γ-glutamyl cyclotransferase (ChaC) is involved in GSH degradation. However, the molecular mechanism underlying GSH degradation in silkworms (Bombyx mori) remains unknown. Silkworms are lepidopteran insects that are considered to be an agricultural pest model. We aimed to examine the metabolic mechanism underlying GSH degradation mediated by B. mori ChaC and successfully identified a novel ChaC gene in silkworms (herein, bmChaC). The amino acid sequence and phylogenetic tree revealed that bmChaC was closely related to mammalian ChaC2. We overexpressed recombinant bmChaC in Escherichia coli, and the purified bmChaC showed specific activity toward GSH. Additionally, we examined the degradation of GSH to 5-oxoproline and cysteinyl glycine via liquid chromatography-tandem mass spectrometry. Quantitative real-time polymerase chain reaction revealed that bmChaC mRNA expression was observed in various tissues. Our results suggest that bmChaC participates in tissue protection via GSH homeostasis. This study provides new insights into the activities of ChaC and the underlying molecular mechanisms that can aid the development of insecticides to control agricultural pests.

Effects of Chromium Exposure on the Gene Expression of the Midgut in Silkworms, Bombyx mori

Chromium is a severe heavy metal pollutant with significant environmental risks. The effects of Chromium on the digestion of Bombyx mori (silkworms) are of particular importance due to their ecological and economic significance. Herein, RNA sequencing was conducted on nine midgut samples from silkworms exposed to control, 12 g/kg and 24 g/kg Chromium chemical diets. Comparative transcriptomics revealed that under moderate Chromium exposure, there was a significant increase in up-regulated genes (1268 up-regulated to 857 down-regulated), indicating a stimulation response. At higher stress levels, a weakened survival response was observed, with a decrease in up-regulated genes and an increase in down-regulated genes (374 up-regulated to 399 down-regulated). A notable shift in cellular responses under medium chromium exposure was exposed, signifying the activation of crucial metabolic and transport systems and an elevation in cellular stress and toxicity mechanisms. The observation of up-regulated gene expression within xenobiotic metabolism pathways suggests a heightened defense against Chromium-induced oxidative stress, which was primarily through the involvement of antioxidant enzymes. Conversely, high-dose Chromium exposure down-regulates the folate biosynthesis pathway, indicating biological toxicity. Two novel genes responsive to pressure were identified, which could facilitate future stress adaptation understanding. The findings provide insights into the molecular mechanisms underlying silkworms' digestion response to Chromium exposure and could inform its biological toxicity.

Homeostatic Regulation of the Duox-ROS Defense System: Revelations Based on the Diversity of Gut Bacteria in Silkworms (Bombyx mori)

The Duox-ROS defense system plays an important role in insect intestinal immunity. To investigate the role of intestinal microbiota in Duox-ROS regulation herein, 16S rRNA sequencing technology was utilized to compare the characteristics of bacterial populations in the midgut of silkworm after different time-periods of treatment with three feeding methods: 1-4 instars artificial diet (AD), 1-4 instars mulberry leaf (ML) and 1-3 instars artificial diet + 4 instar mulberry leaf (TM). The results revealed simple intestinal microbiota in the AD group whilst microbiota were abundant and variable in the ML and TM silkworms. By analyzing the relationship among intestinal pH, reactive oxygen species (ROS) content and microorganism composition, it was identified that an acidic intestinal environment inhibited the growth of intestinal microbiota of silkworms, observed concurrently with low ROS content and a high activity of antioxidant enzymes (SOD, TPX, CAT). Gene expression associated with the Duox-ROS defense system was detected using RT-qPCR and identified to be low in the AD group and significantly higher in the TM group of silkworms. This study provides a new reference for the future improvement of the artificial diet feeding of silkworm and a systematic indicator for the further study of the relationship between changes in the intestinal environment and intestinal microbiota balance caused by dietary alterations.

Effects of artificial diet rearing during all instars on silk secretion and gene transcription in Bombyx mori (Lepidoptera: Bombycidae)

Silkworms (Bombyx mori) reared on artificial diets during all instars have the advantages of simplicity and efficiency, year-round production, and reduced risk of poisoning. However, low silk yield remains a challenge, limiting its industrial application. To address this issue, the spinning behavior, nutrient absorption, and transcriptomics of silkworms were investigated. Compared with silkworms reared on mulberry leaves during all instars, those fed with artificial diets showed significantly lower cocoon weight, cocoon shell weight, cocoon shell rate, and silk gland tissue somatic index at the end of the fifth instar (P < 0.01). The spinning duration and crawling distance of silkworms reared on artificial diets were also significantly lower than those reared on mulberry leaves (P < 0.01). Regarding nutrient absorption, the dietary efficiency indexes of silkworms fed with artificial diets were significantly lower than those fed with mulberry leaves, except for the efficiency conversion of digesta to cocoon (P < 0.01). Further RNA-Seq analysis revealed 386 differentially transcribed genes between the 2 groups, with 242 upregulated and 144 downregulated genes. GO enrichment analysis showed that differential transcriptional genes were mainly enriched in organic acid metabolism, oxidation-reduction, and drug catabolism. KEGG enrichment analysis showed that differential transcriptional genes were mainly enriched in genetic information processing and metabolism pathways. Our findings provide new insights into the silk secretion and can serve as a reference for future research and application of silkworms fed with artificial diets.

Overexpression of BmJHBPd2 Repressed Silk Synthesis by Inhibiting the JH/Kr-h1 Signaling Pathway in Bombyx mori

The efficient production of silkworm silk is crucial to the silk industry. Silk protein synthesis is regulated by the juvenile hormone (JH) and 20-Hydroxyecdysone (20E). Therefore, the genetic regulation of silk production is a priority. JH binding protein (JHBP) transports JH from the hemolymph to target organs and cells and protects it. In a previous study, we identified 41 genes containing a JHBP domain in the Bombyx mori genome. Only one JHBP gene, BmJHBPd2, is highly expressed in the posterior silk gland (PSG), and its function remains unknown. In the present study, we investigated the expression levels of BmJHBPd2 and the major silk protein genes in the high-silk-producing practical strain 872 (S872) and the low-silk-producing local strain Dazao. We found that BmJHBPd2 was more highly expressed in S872 than in the Dazao strain, which is consistent with the expression pattern of fibroin genes. A subcellular localization assay indicated that BmJHBPd2 is located in the cytoplasm. In vitro hormone induction experiments showed that BmJHBPd2 was upregulated by juvenile hormone analogue (JHA) treatment. BmKr-h1 upregulation was significantly inhibited by the overexpression of BmJHBPd2 (BmJHBPd2OE) at the cell level when induced by JHA. However, overexpression of BmJHBPd2 in the PSG by transgenic methods led to the inhibition of silk fibroin gene expression, resulting in a reduction in silk yield. Further investigation showed that in the transgenic BmJHBPd2OE silkworm, the key transcription factor of the JH signaling pathway, Krüppel homolog 1 (Kr-h1), was inhibited, and 20E signaling pathway genes, such as broad complex (Brc), E74A, and ultraspiracle protein (USP), were upregulated. Our results indicate that BmJHBPd2 plays an important role in the JH signaling pathway and is important for silk protein synthesis. Furthermore, our findings help to elucidate the mechanisms by which JH regulates silk protein synthesis.

fruitless is sex-differentially spliced and is important for the courtship behavior and development of silkmoth Bombyx mori

Sexual dimorphisms of the brain play essential roles in successful reproduction. Silkmoth Bombyx mori exhibits extensive sexual differences in sexual behavior, as well as their morphology. Although the neural circuits that transmit information about sex pheromone in the male brain are extensively analyzed, the molecular mechanisms that regulate their development are still elusive. In the present study, we focused on the silkmoth ortholog of fruitless (fru) as a candidate gene that regulates sexual dimorphisms of the brain. fru transcripts were expressed from multiple promoters in various tissues, and brain-specific transcripts were sex-specifically spliced, in a manner similar to Drosophila. Interestingly, fru was highly expressed in the adult female brain and the male larval testis. Analysis of CRISPR/Cas9-mediated fru knockout strains revealed that fru plays important roles in survival during late larval and pupal stages, testis development, and adult sexual behavior. fru mutant males exhibited highly reduced levels of courtship and low copulation rate, indicating that fru plays significant roles in the sexual behavior of silkmoths, although it is not absolutely necessary for copulation. In the fru mutant males, sexually dimorphic pattern of the odorant receptor expression was impaired, possibly causing the defects in courtship behavior. These results provide important clues to elucidate the development of sexual dimorphisms of silkmoth brains, as well as the evolution of fruitless gene in insects.

Diagnosis of Pebrine Disease in Silkworm Using Molecular Methods

Since pebrine disease, as the most important and dangerous disease in silkworms, spreads horizontally through the spores and vertically through the eggs, combating the disease and eliminating it completely from livestock production has been associated with numerous problems. This project aimed to identify the molecular cause of pebrine disease in silkworms using a sensitive, specific, and accurate method. To this purpose, a 136 bp fragment was selected based on the Nosema bombycis partial SSU rDNA sequence, and a pair of primers was designed. Afterward, using the conventional polymerase chain reaction (PCR) method, the target fragment was amplified and sequenced. After that, to determine the detection sensitivity, using the Real-Time PCR method, 5-fold serial dilutions of N. bombycis DNA were prepared, and the last dilution that produced a fluorescent signal was considered the minimum detection limit. All tests were performed in duplicates. Based on the results of the sensitivity test, the standard curve including Ct values ​​and DNA concentration was used for analysis. Moreover, 80 unknown samples examined by light microscope were evaluated using conventional PCR and Real-Time PCR. Both PCR results showed no amplification for the negative control samples. The findings demonstrated that the lowest detection limit for N. bombycis was less than 6 pg of DNA, while, this amount was 8 ng for conventional PCR. Out of 80 samples examined, 55, 60, and 62 samples were positive for light microscope, conventional PCR, and Real-Time PCR methods, respectively. The findings suggested that the Real-Time PCR method had a higher ability to detect the causative agent of pebrine disease than the conventional PCR method, and both methods were superior to light microscopy. Therefore, due to the fewer steps and higher accuracy of Real-Time PCR, it can be introduced as a suitable method for diagnosing pebrine disease.

Disruption of a BTB-ZF transcription factor causes female sterility and melanization in the larval body of the silkworm, Bombyx mori

The dilute black (bd) of the silkworm Bombyx mori is a recessive mutant that produces a grayish-black color in the larval integument, instead of the characteristic white color found in wild-type larvae. In addition, eggs produced by bd females are sterile due to a deficiency in the micropylar apparatus. We identified candidate genes responsible for the bd phenotype using publicly available RNA-seq data. One of these candidate genes was homologous to the maternal gene required for meiosis (mamo) of Drosophila melanogaster, which encodes a broad-complex, tramtrack, and bric-à-brac-zinc finger (BTB-ZF) transcription factor essential for female fertility. In three independent bd strains, the expression of the B. mori mamo (Bmmamo) was downregulated in the larval integument. Using a CRISPR/Cas9-mediated knockout strategy, we found that Bmmamo knockout mutants exhibit a grayish-black color in the larval integument and female infertility. Moreover, larvae obtained from the complementation cross between bd/+ mutants and heterozygous knockouts for the Bmmamo also exhibited a grayish-black color, indicating that Bmmamo is responsible for the bd phenotype. Gene expression analysis using Bmmamo knockout mutants suggested that the BmMamo protein suppresses the expression of melanin synthesis genes. Previous comparative genome analysis revealed that the Bmmamo was selected during silkworm domestication, and we found that Bmmamo expression in the larval integument is higher in B. mori than in the wild silkworm B. mandarina, suggesting that the Bmmamo is involved in domestication-associated pigmentation changes of the silkworm.

Diffusion and enrichment of high-risk antibiotic resistance genes (ARGs) via the transmission chain (mulberry leave, guts and feces of silkworm, and soil) in an ecological restoration area of manganese mining, China: Role of heavy metals

This study investigated the diffusion and enrichment of antibiotic resistance genes (ARGs) and pathogens via the transmission chain (mulberry leaves - silkworm guts - silkworm feces - soil) near a manganese mine restoration area (RA) and control area (CA, away from RA). Horizontal gene transfer (HGT) of ARGs was testified by an IncP a-type broad host range plasmid RP4 harboring ARGs (tetA) and conjugative genes (e.g., korB, trbA, and trbB) as an indicator. Compared to leaves, the abundances of ARGs and pathogens in feces after silkworms ingested leaves from RA increased by 10.8% and 52.3%, respectively, whereas their abundance in feces from CA dropped by 17.1% and 97.7%, respectively. The predominant ARG types in feces involved the resistances to β-lactam, quinolone, multidrug, peptide, and rifamycin. Therein, several high-risk ARGs (e.g., qnrB, oqxA, and rpoB) carried by pathogens were more enriched in feces. However, HGT mediated by plasmid RP4 in this transmission chain was not a main factor to promote the enrichment of ARGs due to the harsh survival environment of silkworm guts for the plasmid RP4 host E. coli. Notably, Zn, Mn, and As in feces and guts promoted the enrichment of qnrB and oqxA. Worriedly, the abundance of qnrB and oqxA in soil increased by over 4-fold after feces from RA were added into soil for 30 days regardless of feces with or without E. coli RP4. Overall, ARGs and pathogens could diffuse and enrich in environment via the sericulture transmission chain developed at RA, especially some high-risk ARGs carried by pathogens. Thus, greater attentions should be paid to dispel such high-risk ARGs to support benign development of sericulture industry in the safe utilization of some RAs.

Juvenile hormone suppresses the FoxO-takeout axis to shorten longevity in male silkworm

Juvenile hormone (JH) plays a crucial endocrine regulatory role in insect metamorphosis, reproduction, and longevity in multiple organisms, such as flies, honeybees, and migratory monarch butterflies. However, the molecular mechanism of JH affecting longevity remains largely unknown. In this study, we showed that JH III and its analog methoprene shortened the survival days significantly in the adulthood of male silkworm. At the same time, the allatostatin, a neuropeptide that inhibits the secretion of JH by the corpora allata, could extend the survival days dramatically after adult eclosion in male silkmoth. Interestingly, a central pro-longevity FoxO transcription factor was reduced upon JH stimulation in silkworm individuals and BmN-SWU1 cells. Furthermore, the analysis of the upstream sequence of the FoxO gene identified a JH response element which suggested that FoxO might be regulated as a target of JH. Surprisingly, we identified a Bmtakeout (BmTO) gene that encodes a JH-binding protein and contains a FoxO response element. As expected, FoxO overexpression and knockdown up- and down-regulated the expression of BmTO respectively, indicating that BmTO functions as a FoxO target. BmTO overexpression could release the inhibitory effect of JH on the BmFoxO gene by reducing JH bioavailability to block its signal transduction. Collectively, these results may provide insights into the mechanism of the JH-FoxO-TO axis in aging research and pest control.

Peptidoglycan Recognition Protein S5 of Bombyx mori Facilitates the Proliferation of Bombyx mori Cypovirus 1

Bombyx mori cypovirus 1 (BmCPV1), a primary pathogen of the silkworm, is a typical dsRNA virus belonging to the Reoviridae family. In this study, a total of 2520 differentially expressed genes (DEGs) were identified by RNA-seq analysis of the silkworm midgut after BmCPV1 infection and Gene Ontology (GO) functional annotation showed that the DEGs predominantly functioned in binding (molecular function), cell (cellular component), and cellular processes (biological process). Additionally, the Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation revealed that the DEGs were mainly distributed in global and overview metabolism maps, translation, and signal transduction. Among the identified DEGs, BmPGRP-S5 belongs to the peptidoglycan recognition protein (PGRP) family. Previous studies have revealed that PGRPs were involved in the interactions between silkworm and BmCPV1. Here, we explored the effect of BmPGRP-S5 on BmCPV1 replication and demonstrated that BmPGRP-S5 promotes the proliferation of BmCPV1 in BmN cells through overexpression or knockdown experiments. Knocking down of BmPGRP-S5 in silkworm larvae similarly promoted the proliferation of BmCPV1. Through experimental validation, we therefore determined that BmPGRP-S5 acts as a proviral host factor for BmCPV1 infection. This study clarifies the proliferation mechanism of BmCPV1 and provides new insights into the functional role of BmPGRP-S5.

Highly accurate genome assembly of an improved high-yielding silkworm strain, Nichi01

The silkworm (Bombyx mori) is an important lepidopteran model insect and an industrial domestic animal traditionally used for silk production. Here, we report the genome assembly of an improved Japanese strain Nichi01, in which the cocoon yield is comparable to that of commercial silkworm strains. The integration of PacBio Sequel II long-read and ddRAD-seq-based high-density genetic linkage map achieved the highest quality genome assembly of silkworms to date; 22 of the 28 pseudomolecules contained telomeric repeats at both ends, and only four gaps were present in the assembly. A total of 452 Mbp of the assembly with an N50 of 16.614 Mbp covered 99.3% of the complete orthologs of the lepidopteran core genes. Although the genome sequence of Nichi01 and that of the previously reported low-yielding tropical strain p50T assured their accuracy in most regions, we corrected several regions, misassembled in p50T, in our assembly. A total of 18,397 proteins were predicted using over 95 Gb of mRNA-seq derived from 10 different organs, covering 96.9% of the complete orthologs of the lepidopteran core genes. The final assembly and annotation files are available in KAIKObase (https://kaikobase.dna.affrc.go.jp/index.html) along with a genome browser and BLAST searching service, which would facilitate further studies and the breeding of silkworms and other insects.

Spindle-E is essential for gametogenesis in the silkworm, Bombyx mori

As a defense mechanism against transposable elements, the PIWI-interacting RNA (piRNA) pathway maintains genomic integrity and ensures proper gametogenesis in gonads. Numerous factors are orchestrated to ensure normal operation of the piRNA pathway. Spindle-E (Spn-E) gene was one of the first genes shown to participate in the piRNA pathway. In this study, we performed functional analysis of Spn-E in the model lepidopteran insect, Bombyx mori. Unlike the germline-specific expression pattern observed in Drosophila and mouse, BmSpn-E was ubiquitously expressed in all tissues tested, and it was highly expressed in gonads. Immunofluorescent staining showed that BmSpn-E was localized in both germ cells and somatic cells in ovary and was expressed in spermatocytes in testis. We used a binary transgenic CRISPR/Cas9 system to construct BmSpn-E mutants. Loss of BmSpn-E expression caused derepression of transposons in gonads. We also found that mutant gonads were much smaller than wild-type gonads and that the number of germ cells was considerably lower in mutant gonads. Quantitative real-time PCR analysis and TUNEL staining revealed that apoptosis was greatly enhanced in mutant gonads. Further, we found that the BmSpn-E mutation impacted gonadal development and gametogenesis at the early larval stage. In summary, our data provided the first evidence that BmSpn-E plays vital roles in gonadal development and gametogenesis in B. mori.