LincRNA_XR209691.3 could promote Bombyx mori nucleopolyhedrovirus replication by interacting with BmHSP70

Long non-coding RNAs (lncRNAs), a class of transcripts exceeding 200 nucleotides and lacking protein coding potential, have been proven to play important roles in viral infection and host immunity. Bombyx mori nucleopolyhedrovirus (BmNPV) is an important pathogen, which causes the silkworm disease and leads to a huge challenge to the sericultural industry. At present, research on the roles of insect lncRNAs in host-virus interaction are relatively few. In this study, we explored the function of lincRNA_XR209691.3 that was significantly up-regulated in the silkworm fat body upon BmNPV infection. Firstly, the subcellular localization experiment confirmed that lincRNA_XR209691.3 was present in both the nucleus and cytoplasm. Enhancing the expression of lincRNA_XR209691.3 in BmN cells could promote the proliferation of BmNPV, while inhibition of lincRNA_XR209691.3 by RNA interference suppresses the proliferation of BmNPV. Combining RNA pull-down and mass spectrometry, we identified the host and BmNPV proteins that could interact with lincRNA_XR209691.3. Next, by using truncation experiment and RNA immunoprecipitation (RIP) assay, it was found that lincRNA_XR209691.3 could bind to the Actin domain of BmHSP70. Subsequently, overexpression of lncRNA_XR209691.3 in BmN cells promoted the expression of BmHSP70, while knockdown of BmHsp70 suppressed the replication of BmNPV. Based on the above results, it is speculated that lincRNA_XR209691.3 could promote the proliferation of BmNPV through interaction with BmHSP70, possibly by improving the stability of BmHSP70 and thereby enhancing the expression of BmHSP70. Our results shed light on the lncRNA function in insect-pathogen interactions and provide a new clue to elucidate the molecular mechanism of BmNPV infection.

Recessive embryonic lethal mutations uncovered in heterozygous condition in silkworm semiconsomic strains

In this study, we found two embryonic lethal mutations, t04 lethal (l-t04) and m04 lethal (l-m04), in semiconsomic strains T04 and M04, respectively. In these semiconsomic strains, the entire diploid genome, except for one chromosome 4 of the wild silkworm Bombyx mandarina, is substituted with chromosomes of the domesticated silkworm B. mori, and l-t04 and l-m04 mutations are located on B. mandarina-derived chromosome 4. To clarify the cause of the lethalities and the genes responsible for these mutations, positional cloning and CRISPR/Cas9 mediated knockout screening were performed. Finally, genetic complementation tests identified the mutations responsible for the l-t04 and l-m04 as the Bombyx homolog of imaginal discs arrested (Bmida) and TATA box binding protein-associated factor 5 (BmTaf5), respectively. Lethal stages of each knockout mutant indicated the importance of these genes in B. mori late embryogenesis. The lethal mutations responsible for l-t04 and l-m04 were not found in parental strains or wild B. mandarina collected from 39 distinct locations in Japan, indicating that both mutations were independently introduced during or after the development of the semiconsomic strains. We conclude that the recessive embryonic lethality in the T04 and M04 strains is due to deleterious mutations produced in B. mandarina-derived chromosome 4.

The tissue-specific expression of silkworm cuticle protein gene ASSCP2 is mediated by the Sox-2 transcription factor

The silk gland of silkworm is a unique organ in which silk proteins are synthesized, secreted, and transformed into fibers. The anterior silk gland (ASG) is located at the end of the silk gland, and is thought to be involved in silk fibrosis. In our previous study, a cuticle protein, ASSCP2, was identified. This protein is specifically and highly expressed in the ASG. In this work, the transcriptional regulation mechanism of ASSCP2 gene was studied by a transgenic route. The ASSCP2 promoter was analyzed, truncated sequentially, and used to initiate the expression of EGFP gene in silkworm larvae. After egg injection, seven transgenic silkworm lines were isolated. Molecular analysis revealed that the green fluorescent signal could not be detected when the promoter was truncated to -257 bp, suggesting that the -357 to -257 sequence is the key region responsible for the transcriptional regulation of the ASSCP2 gene. Furthermore, an ASG specific transcription factor Sox-2 was identified. EMSA assays showed that Sox-2 binds with the -357 to -257 sequence, and thus regulates the tissue-specific expression of ASSCP2. This study on the transcriptional regulation of ASSCP2 gene provides theoretical and experimental basis for further studies of the regulatory mechanism of tissue-specific genes.

Characterization of a glutamate-cysteine ligase in Bombyx mori

Glutamate-cysteine ligase (GCL) is a crucial enzyme involved in the synthesis of glutathione (GSH). Despite various studies on glutathione transferase, and its essential role in detoxification and resistance to oxidative stress, GSH synthesis has not been described in Bombyx mori (silkworms) to date. Silkworms form part of the lepidopterans that are considered as a model of agricultural pests. This study aimed to understand the GSH synthesis by GCL in silkworms, which may help in developing insecticides to tackle agricultural pests. Based on the amino acid sequence and phylogenetic tree, the B. mori GCL belongs to group 2, and is designated bmGCL. Recombinant bmGCL was overexpressed and purified to ensure homogeneity. Biochemical studies revealed that bmGCL uses ATP and Mg²⁺ to ligate glutamate and cysteine. High expression levels of bmgcl mRNA and GSH were observed in the silkworm fat body after exposure to insecticides and UV-B irradiation. Moreover, we found an increase in bmgcl mRNA and GSH content during pupation in the silkworm fat body. In this study, we characterized the B. mori GCL and analyzed its biochemical properties. These observations indicate that bmGCL might play an important role in the resistance to oxidative stress in the silkworms.

The effects of glyphosate exposure on gene transcription and immune function of the silkworm, Bombyx mori

Glyphosate is a widely used herbicide and crop desiccant. However, whether its extensive use has any effect on the species diversity of nontarget organisms is still unclear. In this study, we used the silkworm, Bombyx mori, as the research subject, and performed RNA sequencing to analyze the transcriptional profile of silkworm midgut after exposure to glyphosate at 2975.20 mg/L (a concentration commonly used at mulberry fields). A total of 125 significantly differentially expressed genes (DEGs) were detected in the midgut of glyphosate-exposed silkworm (q < 0.05), of which 53 were upregulated and 72 were downregulated. Gene ontology enrichment analysis showed that the DEGs were mainly enriched in biological process, cellular component, and molecular function. Kyoto encyclopedia of genes and genomes analysis showed that the differential genes were mainly related to oxidative stress, nutrient metabolism, and immune defense pathways, including oxidative stress-related Cat and Jafrac1, nutrient metabolism-related Fatp and Scpx, and immune-related CYP6AN2, UGT40B4, CTL11, serpin-2, and so forth. Experimental verification showed that glyphosate exposure led to a 4.35-fold increase in the mortality of silkworm after Beauveria bassiana infection, which might be caused by the decreased PO (phenoloxidase) activity and impaired immunity. These results provide evidence for the potential effects of residue glyphosate on the physiological functions of silkworm, and also provide a reference for the biosafety evaluation of glyphosate.

The piRNA cluster torimochi is an expanding transposon in cultured silkworm cells

PIWI proteins and PIWI-interacting RNAs (piRNAs) play a central role in repressing transposable elements in animal germ cells. It is thought that piRNAs are mainly produced from discrete genomic loci named piRNA clusters, which often contain many "dead" transposon remnants from past invasions and have heterochromatic features. In the genome of silkworm ovary-derived cultured cells called BmN4, a well-established model for piRNA research, torimochi was previously annotated as a unique and specialized genomic region that can capture transgenes and produce new piRNAs bearing a trans-silencing activity. However, the sequence identity of torimochi has remained elusive. Here, we carefully characterized torimochi by utilizing the updated silkworm genome sequence and the long-read sequencer MinION. We found that torimochi is in fact a full-length gypsy-like LTR retrotransposon, which is exceptionally active and has massively expanded its copy number in BmN4 cells. Many copies of torimochi in BmN4 cells have features of open chromatin and the ability to produce piRNAs. Therefore, torimochi may represent a young, growing piRNA cluster, which is still "alive" and active in transposition yet capable of trapping other transposable elements to produce de novo piRNAs.

Dimmed gene knockout shortens larval growth and reduces silk yield in the silkworm, Bombyx mori

The bHLH domain transcription factor, Bombyx mori-derived dimmed (Bmdimm), is directly regulated by the JH-BmMet/BmSRC-BmKr-h1 pathway and plays a key role in regulating the expression of FibH, which codes the main component of silk protein. However, the other roles of Bmdimm in silk protein synthesis remain unclear. Here, we established a Bmdimm knockout (KO) line containing a 7-bp deletion via CRISPR/Cas9 system, which led to the absence of the bHLH domain. The expression level of silk protein genes and silk yield decreased significantly in the Bmdimm KO line. Moreover, knocking out Bmdimm led to shortened larval stages and significant weight loss in larvae and adults. Bmdimm was found to be highly expressed in the silk gland, but it was also expressed in the fat body. The expression level of Bmkr-h1 in the fat body was significantly downregulated in the Bmdimm KO line. Exogenous JHA treatment upregulated Bmkr-h1 and rescued the phenotype of larval growth in the Bmdimm KO line. In conclusion, knocking out Bmdimm led to a shortened larval stage via the inhibition of Bmkr-h1 expression, then reduced silk yield. These findings help to elucidate the regulatory mechanism of fibroin synthesis and larval development in silkworms.

lncR26319/miR-2834/EndophilinA axis regulates oogenesis of the silkworm, Bombyx mori

Oocyte maturation is critical for insect reproduction. Vitellogenesis, the timely production and uptake of vitellogenin (Vg), is crucial for female fecundity. Vg is synthesized in fat body and absorbed by the oocytes through endocytosis during insect oogenesis. In the silkworm, Bombyx mori, we discovered that a nucleus-enriched long-noncoding RNA (lncRNA) lncR26319 regulates Endophilin A (EndoA) - a member of the endophilin family of endocytic proteins - through competitive binding to miR-2834. The lncR26319-miR-2834-EndoA axis was required for Vg endocytosis in the silkworm; loss of EndoA or overexpression of miR-2834 significantly reduced egg numbers in virgin moths. In addition, accumulation of miR-2834 resulted in pupal and adult deformation and reduced fecundity in females. The expression of Vg, 30-kDa (30K) protein, and egg-specific protein (Esp) decreased after knockdown of EndoA or overexpression of miR-2834, while knockdown of miR-2834 had an opposite effect on the expression of Vg, 30K protein gene, and Esp. These results suggest that the lncR26319-miR-2834-EndoA axis contributes to the endocytic activity in the Vg uptake and leads to the normal progression of oogenesis in the silkworm. Thus, miR-2834 and EndoA are crucial for female reproduction and could be potential targets for new pest management strategies in lepidopterans.

Genome-wide identification and comparative analysis of lipocalin families in Lepidoptera with an emphasis on Bombyx mori

Lipocalins exhibit functional diversity, including roles in retinol transport, invertebrate cryptic coloration, and stress response. However, genome-wide identification and characterization of lipocalin in the insect lineage have not been thoroughly explored. Here, we found that a lineage-specific expansion of the lipocalin genes in Lepidoptera occurred in large part due to tandem duplication events and several lipocalin genes involving insect coloration were expanded more via tandem duplication in butterflies. A comparative analysis of conserved motifs showed both conservation and divergence of lepidopteran lipocalin family protein structures during evolution. We observe dynamic changes in tissue expression preference of paralogs in Bombyx mori, suggesting differential contribution of paralogs to specific organ functions during evolution. Subcellular localization experiments revealed that lipocalins localize to the cytoplasm, nuclear membrane, or nucleus in BmN cells. Moreover, several lipocalin genes exhibited divergent responses to abiotic and biotic stresses, and 1 lipocalin gene was upregulated by 300 fold in B. mori. These results suggest that lipocalins act as signaling components in defense responses by mediating crosstalk between abiotic and biotic stress responses. This study deepens our understanding of the comprehensive characteristics of lipocalins in insects.

Temporal transcriptome reveals that circadian clock is involved in the dynamic regulation of immune response to bacterial infection in Bombyx mori

The circadian clock plays a critical role in the regulation of host immune defense. However, the mechanistic basis for this regulation is largely unknown. Herein, the core clock gene cryptochrome1 (cry1) knockout line in Bombyx mori, an invertebrate animal model, was constructed to obtain the silkworm with dysfunctional molecular clock, and the dynamic regulation of the circadian clock on the immune responsiveness within 24 h of Staphylococcus aureus infection was analyzed. We found that deletion of cry1 decreased viability of silkworms and significantly reduced resistance of larvae to S. aureus. Time series RNA-seq analysis identified thousands of rhythmically expressed genes, including immune response genes, in the larval immune tissue, fat bodies. Uninfected cry1 knockout silkworms exhibited expression patterns of rhythmically expressed genes similar to wild-type (WT) silkworms infected with S. aureus. However, cry1 knockout silkworms exhibited a seriously weakened response to S. aureus infection. The immune response peaked at 6 and 24 h after infection, during which "transcription storms" occurred, and the expression levels of the immune response genes, PGRP and antimicrobial peptides (AMPs), were significantly upregulated in WT. In contrast, cry1 knockout did not effectively activate Toll, Imd, or NF-κB signaling pathways during the immune adjustment period from 12 to 18 h after infection, resulting in failure to initiate the immune responsiveness peak at 24 h after infection. This may be related to inhibited silkworm fat body energy metabolism. These results demonstrated the dynamic regulation of circadian clock on silkworm immune response to bacterial infection and provided important insights into host antimicrobial defense mechanisms.

The LIM Domain Protein BmFHL2 Inhibits Egg Production in Female Silkworm, Bombyx mori

The female Bombyx mori accumulates a large amount of egg proteins, mainly Vg and 30K, during egg formation to provide nutrition for embryo development. The synthesis and transport of Vg have been extensively studied, particularly the regulation of Vg transcription induced by 20E; however, the mechanism of 30K protein synthesis is poorly studied. As a model organism of the order Lepidoptera, B. mori has high reproduction potential. In the present study, we found that the FHL2 homologous gene (BmFhl2) in B. mori is involved in inhibiting female egg formation by influencing the synthesis of 30K protein. Interference of BmFhl2 expression in silkworm females increased 30K protein synthesis, accelerated ovarian development, and significantly increased the number of eggs produced and laid; however, the 20E pathway was inhibited. The transcription levels of Vg and 30Kc19 were significantly downregulated following BmFhl2 overexpression in the silkworm ovarian cell line BmN. The Co-IP assay showed that the potential binding protein of BmFHL2 included three types of 30K proteins (30Kc12, 30Kc19, and 30Kc21). These results indicate that BmFHL2 participates in egg formation by affecting 30K protein in female B. mori.

The Mechanisms of Silkworm Resistance to the Baculovirus and Antiviral Breeding

Silkworm (Bombyx mori) is not only an economic insect but also a model organism for life science research. Bombyx mori nucleopolyhedrovirus (BmNPV) disease is a major infectious disease in the world's sericulture industry. The cocoon loss caused by this disease accounts for more than 60% of the total loss caused by all silkworm diseases. To date, there has been no effective solution for preventing and treating this disease. The most effective measure is to breed disease-resistant varieties. The quickest way to breed disease-resistant varieties is to apply genetic modification. However, this requires that we obtain disease resistance genes and know the mechanism of disease resistance. Since the discovery of disease-resistant resources in 1989, scholars in the sericulture industry around the world have been inspired to search for resistance genes. In the past two decades, with the help of multi-omics technologies, screening of resistance genes, gene localization, protein modification, virus-host interactions, etc., researchers have found some candidate genes that have been proposed to function at the cellular or individual level. Several disease-resistant varieties have been obtained and used in production through hybrid breeding, RNA interference, and genetic modification. This article summarizes and reviews the discovery of and research advances related to silkworm resistance to BmNPV. It is anticipated that the review will inspire scientific researchers to continue searching for disease resistance genes, clarify the molecular mechanism of silkworm disease resistance, and promote disease-resistant silkworm breeding.

Expression of tyrosine phosphatases in relation to PTTH-stimulated ecdysteroidogenesis in prothoracic glands of the silkworm, Bombyx mori

Protein tyrosine phosphorylation is a reversible, dynamic process regulated by the activities of tyrosine kinases and tyrosine phosphatases. Although the involvement of tyrosine kinases in the prothoracicotropic hormone (PTTH)-stimulated ecdysteroidogenesis in insect prothoracic glands (PGs) has been documented, few studies have been conducted on the involvement of protein tyrosine phosphatases (PTPs) in PTTH-stimulated ecdysteroidogenesis. In the present study, we investigated the correlation between PTPs and PTTH-stimulated ecdysteroidogenesis in Bombyx mori PGs. Our results showed that the basal PTP enzymatic activities exhibited development-specific changes during the last larval instar and pupation stage, with high activities being detected during the later stages of the last larval instar. PTP enzymatic activity was stimulated by PTTH treatment both in vitro and in vivo. Pretreatment with phenylarsine oxide (PAO) and benzylphosphonic acid (BPA), two chemical inhibitors of tyrosine phosphatase, reduced PTTH-stimulated enzymatic activity. Determination of ecdysteroid secretion showed that treatment with PAO and BPA did not affect basal ecdysteroid secretion, but greatly inhibited PTTH-stimulated ecdysteroid secretion, indicating that PTTH-stimulated PTP activity is indeed involved in ecdysteroid secretion. PTTH-stimulated phosphorylation of the extracellular signal-regulated kinase (ERK) and 4E-binding protein (4E-BP) was partially inhibited by pretreatment with either PAO or BPA, indicating the potential link between PTPs and phosphorylation of ERK and 4E-BP. In addition, we also found that in vitro treatment with 20-hydroxyecdysone did not affect PTP enzymatic activity. We further investigated the expressions of two important PTPs (PTP 1B (PTP1B) and the phosphatase and tension homologue (PTEN)) in Bombyx PGs. Our immunoblotting analysis showed that B. mori PGs contained the proteins of PTP1B and PTEN, with PTP1B protein undergoing development-specific changes. Protein levels of PTP1B and PTEN were not affected by PTTH treatment. The gene expression levels of PTP1B and PTEN showed development-specific changes. From these results, we suggest that PTTH-regulated PTP signaling may crosstalk with ERK and target of rapamycin (TOR) signaling pathways and is a necessary component for stimulation of ecdysteroid secretion.

Dipteran endoparasitoid Exorista bombycis utilizes antihemocyte components against host defense of silkworm Bombyx mori

Dipteran endoparasitoids avoid host immune response; however, antidefense components from the Dipterans are unknown. Infestation of commercial silkworm Bombyx mori Linnaeus (Lepidoptera: Bombycidae) by endoparasitoid Exorista bombycis Louis (Diptera: Tachinidae) induced immune reactions, cytotoxicity, granulation, degranulation, and augmented release of cytotoxic marker enzyme lactate dehydrogenase (LDH), and degranulation-mediator enzyme β-hexosaminidase in hemocytes. In this study, by reverse phase high-performance liquid chromatography, fractions of E. bombycis larval tissue protein with antihemocytic activity are separated. From the fraction, peptides of hemocyte aggregation inhibitor protein (HAIP) and pyridoxamine phosphate oxidase (PNPO) are identified by mass spectrometry. Interacting partners of HAIP and PNPO are retrieved that further enhance the virulence of the parasitoid. PNPO and HAIP genes showed a four- to seven fold increase in expression in the integument of the parasitoid larva. Together, the dipteran endoparasitoid E. bombycis exploit antihemocyte activity to inhibit host defense reactions in addition to defense evasion contemplated.

The inhibition of ecdysone signal pathway was the key of pyriproxyfen poisoning for silkworm, Bombyx mori

Pyriproxyfen is a juvenile hormone-like pesticide. Once intake occurs, it leads to a series of poisoning characters consequences in silkworm, Bombyx mori (ID: 7091, Lepidoptera), such as non- cocooning, non-pupation, production of low-active eggs, and extended stages. However, the poisoning mechanism is still unclear. Here, silkworms were fed mulberry leaves soaked with different pyriproxyfen concentrations, and the heads were dissected for transcriptome analysis, while the hemolymph was used for determinations of ecdysone and juvenile hormone titers. As a result, after conjoint analysis of 3 feeding groups and a control group, 555 differentially expressed genes (DEGs) were obtained, which were mainly involved in hormone metabolism, glycometabolism and protein metabolism. Meanwhile, 119 genes were significantly correlated with the pyriproxyfen concentrations, and they were mainly involved in drug metabolism and glycometabolism. The ecdysone titers in several feeding groups were significantly lower than those of the control group, while juvenile hormone was not detected in all groups, including the control and feeding groups. Correspondingly, due to activation of the juvenile hormone signaling pathway by pyriproxyfen, key genes in the ecdysone synthesis pathway were downregulated, and a large number of downstream genes were up- or downregulated. In addition, nearly all genes in the detoxification pathway were upregulated. These results suggested that, affected by the juvenile hormone signaling pathway, ecdysone titers decreased and further affected a series of downstream processes, and this was the key reason for pyriproxyfen poisoning in silkworm, B. mori, which could lay a foundation for the study of pyriproxyfen resistance in silkworm.

Genome-Wide Characterization and Comparative Genomic Analysis of the Serpin Gene Family in Microsporidian Nosema bombycis

Microsporidia are ubiquitous in the environment, infecting almost all invertebrates, vertebrates, and some protists. The microsporidian Nosema bombycis causes silkworms pébrine disease and leads to huge economic losses. Parasite secreted proteins play vital roles in pathogen-host interactions. Serine protease inhibitors (serpins), belonging to the largest and most broadly distributed protease inhibitor superfamily, are also found in Microsporidia. In this study, we characterized 19 serpins (NbSPNs) in N. bombycis; eight of them were predicted with signal peptides. All NbSPN proteins contain a typical conserved serpin (PF00079) domain. The comparative genomic analysis revealed that microsporidia serpins were only found in the genus Nosema. In addition to N. bombycis, a total of 34 serpins were identified in another six species of Nosema including N. antheraeae (11), N. granulosis (8), Nosema sp. YNPr (3), Nosema sp. PM-1 (3), N. apis (4), and N. ceranae (5). Serpin gene duplications in tandem obviously occurred in Nosema antheranae. Notably, the NbSPNs were phylogenetically clustered with serpins from the Chordopoxvirinae, the subfamily of Poxvirus. All 19 NbSPN transcripts were detected in the infected midgut and fat body, while 19 NbSPN genes except for NbSPN12 were found in the transcriptome of the infected silkworm embryonic cell line BmE-SWU1. Our work paves the way for further study of serpin function in microsporidia.

Construction of scFv Antibodies against the Outer Loops of the Microsporidium Nosema bombycis ATP/ADP-Transporters and Selection of the Fragment Efficiently Inhibiting Parasite Growth

Traditional sanitation practices remain the main strategy for controlling Bombyx mori infections caused by microsporidia Nosema bombycis. This actualizes the development of new approaches to increase the silkworm resistance to this parasite. Here, we constructed a mouse scFv library against the outer loops of N. bombycis ATP/ADP carriers and selected nine scFv fragments to the transporter, highly expressed in the early stages of the parasite intracellular growth. Expression of selected scFv genes in Sf9 cells, their infection with different ratios of microsporidia spores per insect cell, qPCR analysis of N. bombycis PTP2 and Spodoptera frugiperda COXI transcripts in 100 infected cultures made it possible to select the scFv fragment most effectively inhibiting the parasite growth. Western blot analysis of 42 infected cultures with Abs against the parasite β-tubulin confirmed its inhibitory efficiency. Since the VL part of this scFv fragment was identified as a human IgG domain retained from the pSEX81 phagemid during library construction, its VH sequence should be a key antigen-recognizing determinant. Along with the further selection of new recombinant Abs, this suggests the searching for its natural mouse VL domain or "camelization" of the VH fragment by introducing cysteine and hydrophilic residues, as well as the randomization of its CDRs.

Allele-specific knockouts reveal a role for apontic-like in the evolutionary loss of larval melanin pigmentation in the domesticated silkworm, Bombyx mori

The domesticated silkworm, Bombyx mori, and its wild progenitor, B. mandarina, are extensively studied as a model case of the evolutionary process of domestication. A conspicuous difference between these species is the dramatic reduction in melanin pigmentation in both larval and adult B. mori. Here we evaluate the efficiency of CRISPR/Cas9-targeted knockouts of pigment-related genes as a tool to understand their potential contributions to domestication-associated melanin pigmentation loss in B. mori. To demonstrate the efficacy of targeted knockouts in B. mandarina, we generated a homozygous CRISPR/Cas9-targeted knockout of yellow-y. In yellow-y knockout mutants, black body colour became lighter throughout the larval, pupal and adult stages, confirming a role for this gene in melanin pigment formation. Further, we performed allele-specific CRISPR/Cas9-targeted knockouts of the pigment-related transcription factor, apontic-like (apt-like) in B. mori × B. mandarina F1 hybrid individuals which exhibit B. mandarina-like larval pigmentation. Knockout of the B. mandarina allele of apt-like in F1 embryos results in white patches on the dorsal integument of larvae, whereas corresponding knockouts of the B. mori allele consistently exhibit normal F1 larval pigmentation. These results demonstrate a contribution of apt-like to the evolution of reduced melanin pigmentation in B. mori. Together, our results demonstrate the feasibility of CRISPR/Cas9-targeted knockouts as a tool for understanding the genetic basis of traits associated with B. mori domestication.

Bmmp influences wing morphology by regulating anterior-posterior and proximal-distal axes development

Insect wings are subject to strong selective pressure, resulting in the evolution of remarkably diverse wing morphologies that largely determine flight capacity. However, the genetic basis and regulatory mechanisms underlying wing size and shape development are not well understood. The silkworm Bombyx mori micropterous (mp) mutant exhibits shortened wing length and enlarged vein spacings, albeit without changes in total wing area. Thus, the mp mutant comprises a valuable genetic resource for studying wing development. In this study, we used molecular mapping to identify the gene responsible for the mp phenotype and designated it Bmmp. Phenotype-causing mutations were identified as indels and single nucleotide polymorphisms in noncoding regions. These mutations resulted in decreased Bmmp messenger RNA levels and changes in transcript isoform composition. Bmmp null mutants were generated by clustered regularly interspaced short palindromic repeats (CRISPR) / CRISPR-associated protein 9 and exhibited changed wing shape, similar to mp mutants, and significantly smaller total wing area. By examining the expression of genes critical to wing development in wildtype and Bmmp null mutants, we found that Bmmp exerts its function by coordinately modulating anterior-posterior and proximal-distal axes development. We also studied a Drosophila mp mutant and found that Bmmp is functionally conserved in Drosophila. The Drosophila mp mutant strain exhibits curly wings of reduced size and a complete loss of flight capacity. Our results increase our understanding of the mechanisms underpinning insect wing development and reveal potential targets for pest control.

Establishment and application of a silkworm CRISPR/Cas9 tool for conditionally manipulating gene disruption in the epidermis

Spatial or temporal specific gene knockout system is a valuable tool for studying the molecular mechanisms underlying developmental processes. The integument is essential for insect fitness and survival, but tools for dissecting function of genes in this tissue are lacking. In this study, we firstly identified an epidermis specifically expressed gene of the domesticated silkworm, BmCPG25, by comparative transcriptomic analysis. Furthermore, a transgenic silkworm expressing the RNA dependent CRISPR-Cas9 protein driven by the regulatory region of the BmCPG25 was established. Immunochemistry analysis showed the endonuclease was specifically expressed in the nuclear of epidermal cells. We also validated the efficiency of this system by disrupting the function of an epidermis specifically expressed cuticular protein gene (Cpr21) and a ubiquitously expressed cuticular gene (Cph18), respectively. In summary, we successfully constructed a conditional knockout toolkit to manipulate the gene editing in epidermal cells, which provides a valuable approach to study the molecular mechanism of integument development.

BmHen1 is essential for eupyrene sperm development in Bombyx mori but PIWI proteins are not

In lepidopteran insects, sperm dimorphism is a remarkable feature, in which males exhibit two different types of sperms. Both sperm morphs are essential for fertilization: Eupyrene sperm carry DNA and fertilize eggs, whereas apyrene sperm, which do not have nuclei, are necessary for transport of eupyrene sperm into eggs. In this study, we showed that the gene BmHen1, which encodes a methyltransferase that modifies piRNAs, is necessary for eupyrene sperm development in the lepidopteran model insect, Bombyx mori. Loss-of-function mutants of BmHen1 of both sexes were sterile. BmHen1 female mutants laid fewer eggs than wild-type females, and the eggs laid had morphological defects. Immunofluorescence analysis of BmHen1 male mutants revealed that nuclei formation in the eupyrene sperm was defective, whereas apyrene sperm were normal. In mice, worms, and flies, the components in piRNA biogenesis pathway play an important role in gonad development; therefore, we constructed mutations in genes encoding core elements in the piRNA biogenesis pathway, Siwi, and BmAgo3. To our surprise, no obvious phenotypes were observed in the male reproduction system in the Siwi and BmAgo3 mutants, which demonstrated that sperm development in B. mori does not depend on piRNAs. As the sperm development phenotype in BmHen1 mutants mimics the phenotype of the BmPnldc1 mutants, we then performed RNA sequencing analysis of sperm bundles from both mutants. We found that the defects in eupyrene sperm resulted from dysregulation of the expression of genes involved in energy metabolism. Taken together, our findings demonstrate the crucial functions of BmHen1 in the development of eupyrene sperm and provide evidence that spermatogenesis in B. mori is PIWI-independent. Our results suggest potential targets for lepidopteran pest control and broaden our knowledge of the reproduction in this order of insects.

Construction of gene co-expression networks in cultured silkworm cells and identification of previously uncharacterized lepidopteran-specific genes required for chromosome dynamics

Advances in sequencing technology and bioinformatics have accelerated gene discovery and homology-based functional annotation in many species, and numerous targeted gene studies have greatly expanded the understanding of gene functions. Nevertheless, there are still many genes that lack homology with genes in other evolutionary lineages and are left as genes with unknown functions. We constructed a gene co-expression network from the Bombyx mori ovary-derived cell line, BmN4, and attempted to infer the biological roles of uncharacterized genes based on the correlation between the function-known and unknown genes. Within this network, we focused on the co-expression modules involved in chromosome architecture, dynamics, and integrity, and selected the uncharacterized genes for subsequent RNAi-based phenotypic screening. This approach enabled the identification of 5 genes whose knockdown led to abnormalities in chromosome dynamics and spindle morphology in mitosis. One of them was a recently characterized gene, BmCenp-T, which plays a central role in building the kinetochore protein complex on the silkworm holocentric chromosomes. In this study, we suggest a method for constructing the gene co-expression network and selecting candidate genes for small-scale RNAi screening. This approach is complementary to homology-based annotation and may be useful for the analysis of lineage-specific uncharacterized genes such as orphan genes.

Toxic effects of triphenyltin on the silkworm Bombyx mori as a lepidopterous insect model

Triphenyltin (TPT) is a widely used reagent in various industries and agriculture, but is also known to accumulate in natural ecosystems and animal tissues. Hence, the aim of this study was to comprehensively assess the toxicity of TPT in the silkworm Bombyx mori as a model insect. The results showed that TPT exposure for the entire 5th instar larval stage significantly reduced the weight of silkworm pupa and inhibited development of the silkworm midgut. Following exposure to 2 μg/kg of TPT for 4 days, differentially expressed genes in midgut were associated with enriched pathways involved in the metabolism of carbohydrates, lipids, and amino acids, as determined by RNA sequencing. Furthermore, the metabolic profiles of the intestinal content of silkworms exposed to 2 μg/kg of TPT for 4 days were markedly altered and differential metabolites produced by metabolism of carbohydrates, lipids, and amino acids were enriched as determined by non-targeted GC-MS/MS metabolomics. This study provides novel insights into the mechanisms underlying the toxicity of TPT and emphasizes the risks posed by such pollutants released into the environment.

Accelerated degradation of cellulose in silkworm excrement by the interaction of housefly larvae and cellulose-degrading bacteria

The environmental pollution caused by silkworm (Bombyx mori) excrement is prominent, and rich in refractory cellulose is the bottleneck restricting the efficient recycling of silkworm excrement. This study was performed to investigate the effects of housefly larvae vermicomposting on the biodegradation of cellulose in silkworm excrement. After six days, a 58.90% reduction of cellulose content in treatment groups was observed, which was significantly higher than 11.5% of the control groups without housefly larvae. Three cellulose-degrading bacterial strains were isolated from silkworm excrement, which were identified as Bacillus licheniformis, Bacillus amyloliquefaciens, and Bacillus subtilis based on 16S rRNA gene sequence analysis. These three bacterial stains had a high cellulose degradation index (HC value ranged to between 1.86 and 5.97 and FPase ranged from 5.07 U/mL to 7.31 U/mL). It was found that housefly larvae increased the abundance of cellulose-degrading bacterial genus (Bacillus and Pseudomonas) by regulating the external environmental conditions (temperature and pH). Carbohydrate metabolism was the bacterial communities' primary function during vermicomposting based on the PICRUSt. The results of Tax4Fun indicated that the abundance of endo-β-1,4-glucanase and exo-β-1,4-glucanase increased rapidly and maintained at a higher level in silkworm excrement due to the addition of housefly larvae, which contributed to the accelerated degradation of cellulose in silkworm excrement. The finding of this investigation showed that housefly larvae can significantly accelerate the degradation of cellulose in silkworm excrement by increasing the abundance of cellulose-degrading bacterial genera and cellulase.

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

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