The piRNA response to BmNPV infection in the silkworm fat body and midgut

Bombyx mori nucleopolyhedrovirus (BmNPV) is a DNA virus that causes huge losses to the silkworm industry but the piRNA responses during BmNPV infection in the silkworm remain uninvestigated. Here, silkworm piRNA profiles of uninfected and BmNPV-infected fat body and midgut were determined by high-through sequencing in the early stages of BmNPV infection. A total of 2675 and 3396 genome-derived piRNAs were identified from fat body and midgut, respectively. These genome-derived piRNAs mainly originated from unannotated instead of transposon regions in the silkworm genome. In total, 572 piRNAs were associated with 280 putative target genes in fat body and 805 piRNAs with 380 target genes in midgut. Compared to uninfected tissues, 322 and 129 piRNAs were significantly upregulated in BmNPV-infected fat body and midgut, respectively. In addition, 276 and 117 piRNAs were significantly downregulated. Moreover, differentially expressed (DE) piRNAs during BmNPV infection differed significantly between fat body and midgut. Putative DE piRNA-targeted genes were associated with "response to stimulus" and "environmental information processing" in fat body after infection with BmNPV, which may indicate an active piRNA response to BmNPV infection in fat body. This study may lay the foundation for future research of the potential roles of the piRNA pathway and specific piRNAs in BmNPV pathogenesis.

The digestive proteinase trypsin, alkaline A contributes to anti-BmNPV activity in silkworm (Bombyx mori)

Bombyx mori nucleopolyhedrovirus (BmNPV) is a serious pathogenic microorganism that causes tremendous loss to sericulture. Previous studies have found that some proteins of serine protease family in the digestive juice of B. mori larvae have anti-BmNPV activity. In our previous publication about proteome analysis of the digestive juice of B. mori larvae, the digestive enzyme trypsin, alkaline A (BmTA) was filtered as a differentially expressed protein possibly involved in BmNPV resistance. Here, the biological characteristics and anti-BmNPV functions of BmTA were comprehensively analysed. The cDNA sequence of BmTA had an ORF of 768 nucleotides encoding 255 amino acid residues. Domain architecture analysis showed that BmTA contained a signal peptide and a typical Tryp_SPc domain. Quantitative real-time PCR analysis showed that BmTA was highly expressed in the larval stages and specifically expressed in the midgut of B. mori larvae. The expression level of BmTA in BmNPV resistant strain A35 was higher than that in susceptible strain P50. After BmNPV infection, the expression of BmTA increased in both strains from 24 to 72 h. Virus amplification analysis showed that the relative levels of VP39 in B. mori larvae and BmN cells infected with the appropriate concentration of recombinant-BmTA-treated BmNPV were significantly lower than in the control groups. Moreover, overexpression of BmTA in BmN cells significantly inhibited the amplification of BmNPV. Taken together, the results of this study indicated that BmTA possessed anti-BmNPV activity in B. mori, which broadens the horizon for virus-resistant breeding of silkworms.

Genome-wide CRISPR-Cas9 screening in Bombyx mori reveals the toxicological mechanisms of environmental pollutants, fluoride and cadmium

Fluoride and cadmium, two typical environmental pollutants, have been extensively existed in the ecosystem and severely injured various organisms including humans. To explore the toxicological properties and the toxicological mechanism of fluoride and cadmium in silkworm, we perform a CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) -based functional genomic screen, which can directly measure the genetic requirement of genes in response to the pollutants. Our screen identifies 751 NaF-resistance genes, 753 NaF-sensitive genes, 757 CdCl₂-resistance genes, and 725 CdCl₂-sensitive genes. The top-ranked resistant genes are experimentally verified and the results show that their loss conferred resistance to fluoride or cadmium. Functional analysis of the resistant- and sensitive-genes demonstrates enrichment of multiple signaling pathways, among which the MAPK signaling pathway and DNA damage and repair are both required for fluoride- or cadmium-induced cell death, whereas the Toll and Imd signaling pathway and Autophagy are fluoride- or cadmium-specific. Moreover, we confirm that these pathways are truly involved in the toxicological mechanism in both cultured cells and individual tissues. Our results supply potential targets for rescuing the biohazards of fluoride and cadmium in silkworm, and reveal the feasible toxicological mechanism, which highlights the role of functional genomic screens in elucidating the toxicity mechanisms of environmental pollutants.

Proteomic Identification of Bombyx mori Organelles Using the Engineered Ascorbate Peroxidase APEX and Development of Silkworm Organelle Proteome Database (SilkOrganPDB)

Silkworm Bombyx mori is an economically important insect and a lepidopteran model. Organelle proteome is vital to understanding gene functions; however, it remains to be identified in silkworm. Here, using the engineered ascorbate peroxidase APEX, we constructed transgenic B. mori embryo cells (BmE) expressing APEX-NLS, COX4-APEX, APEX-Rev, and APEX-KDEL in nucleus, mitochondrial matrix (MM), cytosol, and endoplasmic reticulum (ER), and isolated the biotin-labeled proteins using streptavidin-affinity purification, respectively. The isolated proteins were determined using LC-MS/MS and annotated by searching B. mori genomes downloaded from GenBank, SilkBase, SilkDB 2.0, and SilkDB 3.0, resulting in 842, 495, 311, and 445 organelle proteins identified, respectively. We mapped the 296 MM proteins annotated in the GenBank data to mitochondrial protein databases of the fly, human, and mouse, and found that 140 (47%) proteins are homologous to 80 fly proteins, and 65 (22%) proteins match to 31 and 29 human and mouse proteins, respectively. Protein orthology was predicted in multiple insects using OrthoMCL, producing 460 families containing 839 proteins we identified. Out of 460 families, 363 were highly conserved and found in all insects, leaving only three proteins without orthology in other insects, indicating that the identified proteins are highly conserved and probably play important roles in insects. A gene ontology enrichment analysis by clusterProfiler revealed that the nucleus proteins significantly enriched in cellular component terms of nucleus and nucleolus, the MM proteins markedly enriched in molecular function terms of nucleotide binding, and the cytosol proteins mainly enriched in biological process terms of small molecule metabolism. To facilitate the usage and analysis of our data, we developed an open-access database, Silkworm Organelle Proteome Database (SilkOrganPDB), which provides multiple modules for searching, browsing, downloading, and analyzing these proteins, including BLAST, HMMER, Organelle Proteins, Protein Locations, Sequences, Gene Ontology, Homologs, and Phylogeny. In summary, our work revealed the protein composition of silkworm BmE organelles and provided a database resource helpful for understanding the functions and evolution of these proteins.

Global Metabolic Profiling of Baculovirus Infection in Silkworm Hemolymph Shows the Importance of Amino-Acid Metabolism

Viruses rely on host cell metabolism to provide the necessary energy and biosynthetic precursors for successful viral replication. Infection of the silkworm, Bombyx mori, by Bombyx mori nucleopolyhedrovirus (BmNPV), has been studied extensively in the past to unravel interactions between baculoviruses and their lepidopteran hosts. To understand the interaction between the host metabolic responses and BmNPV infection, we analyzed global metabolic changes associated with BmNPV infection in silkworm hemolymph. Our metabolic profiling data suggests that amino acid metabolism is strikingly altered during a time course of BmNPV infection. Amino acid consumption is increased during BmNPV infection at 24 h post infection (hpi), but their abundance recovered at 72 hpi. Central carbon metabolism, on the other hand, particularly glycolysis and glutaminolysis, did not show obvious changes during BmNPV infection. Pharmacologically inhibiting the glycolytic pathway and glutaminolysis also failed to reduce BmNPV replication, revealing that glycolysis and glutaminolysis are not essential during BmNPV infection. This study reveals a unique amino acid utilization process that is implemented during BmNPV infection. Our metabolomic analysis of BmNPV-infected silkworm provides insights as to how baculoviruses induce alterations in host metabolism during systemic infection.

Sage controls silk gland development by regulating Dfd in Bombyx mori

Silk gland is an organ that produces and secretes silk proteins. The development of the silk gland is essential for high silk production yield and silk quality. Although Sage reportedly plays a pivotal role in embryonic silk gland development, the mechanism underlying its action remains unclear. Our study aimed to determine the genes downstream of Sage through which it regulates the development of the silk gland. After chromatin immunoprecipitation and sequencing, Dfd was identified as a downstream target gene of Sage and it was confirmed that Sage could inhibit Dfd expression by competing with SGF1. When Dfd was knocked down through RNA interference (RNAi), the number of cells in the middle silk gland decreased, and the posterior silk gland was straightened. Simultaneously, the expression of Ser1 and silk fibroin genes was no longer strictly regional. These changes eventually led to an alteration in the composition of the Dfd RNAi cocoon. In conclusion, our research contributes to a deeper understanding of the development of silk glands.

Transcriptome analysis reveals changes in silkworm energy metabolism during Nosema bombycis infection

Energy metabolism is important for the proliferation of microsporidia in infected host cells, but there is limited information on the host response. The energy metabolism response of silkworm (Bombyx mori) to microsporidia may help manage Nosema bombycis infections. We analyzed differentially expressed genes in the B.mori midgut transcriptome at two significant time points of microsporidia infection. A total of 1448 genes were up-regulated, while 315 genes were down-regulated. A high proportion of genes were involved in the phosphatidylinositol signaling system, protein processing in the endoplasmic reticulum, and glycerolipid metabolism at 48 h post infection (h p.i.), and a large number of genes were involved in the TCA cycle and protein processing at 120 h p.i. These results showed that the early stages of microsporidia infection affected the basic metabolism and biosynthesis processes of the silkworm. Knockout of Bm_nscaf2860_46 (Bombyx mori isocitrate dehydrogenase, BmIDH) and Bm_nscaf3027_062 (Bombyx mori hexokinase, BmHXK) reduced the production of ATP and inhibited microsporidia proliferation. Host fatty acid degradation, glycerol metabolism, glycolysis pathway, and TCA cycle response to microsporidia infection were also analyzed, and their importance to microsporidia proliferation was verified. These results increase our understanding of the molecular mechanisms involved in N. bombycis infection and provide new insights for research on microsporidia control. IMPORTANCE: Nosema bombycis can be vertically transmitted in silkworm eggs. The traditional prevention and control strategies for microsporidia are difficult and time-consuming, and this is a problem in silkworm culture. Research has mainly focused on host gene functions related to microsporidia infection and host immune responses after microsporidia infection. Little is known about the metabolic changes occurring in the host after infection. Understanding the metabolic changes in the silkworm host could aid in the recognition of host genes important for microsporidia infection and growth. We analyzed host metabolic changes and the main participating pathways at two time points after microsporidia infection and screened the microsporidia-dependent host energy metabolism genes BmIDH and BmHXK. The results revealed genes that are important for the proliferation of Nosema bombycis. These results illustrate how microsporidia hijack the host genome for their growth and reproduction.

Low concentration acetamiprid-induced oxidative stress hinders the growth and development of silkworm posterior silk glands

Acetamiprid is a new type of nicotinic insecticide that is widely used in pest control. Its environmental residues may cause silkworm cocooning disorder. In this study, silkworms that received continuous feeding of low concentration acetamiprid (0.15 mg/L) showed significantly decreased silk gland index and cocooning rate. Gene expression profiling of posterior silk glands (PSGs) revealed that the differentially expressed genes were significantly enriched in oxidative stress-related signal pathways with significant up-regulation. The contents of both H₂O₂ and MDA were increased, along with significantly elevated SOD and CAT activities, all of which reached maximal values at 48 h when H₂O₂ and MDA's contents were 10.46 and 7.98 nmol/mgprot, respectively, and SOD and CAT activities were 5.51 U/mgprot and 33.48 U/gprot, respectively. The transcription levels of antioxidant enzyme-related genes SOD, Mn-SOD, CuZn-SOD, CAT, TPX and GPX were all up-regulated, indicating that exposure to low concentration acetamiprid led to antioxidant response in silkworm PSG. The key genes in the FoxO/CncC/Keap1 signaling pathway that regulates antioxidant enzyme activity, FoxO, CncC, Keap1, NQO1, HO-1 and sMaf were all up-regulated during the whole process of treatment, with maximal values being reached at 72 h with 2.91, 1.46, 1.82, 2.52, 2.32 and 4.01 times of increases, respectively. These results demonstrate that exposure to low concentration acetamiprid causes oxidative stress in silkworm PSG, which may be the cause of cocooning disorder in silkworm. Our study provides a reference for the safety evaluation of environmental residues of acetamiprid on non-target insects.

Genetic analysis and transcriptome analysis of the mini mutant of the silkworm, Bombyx mori

The expression levels of some intrinsic genes, protease activity, and regulation of signaling pathways were distinct during different growth and development stages in the silkworm, Bombyx mori. The silkworm mutant mini was discovered from the normal silkworm strain S8V, and the body-size of the mini mutant was smaller than the wild-type from the second-instar and the difference became more significant in the following stages. In this study, genetic analysis of mini mutant showed that mini mutant was controlled by a single recessive gene, manifested as homozygous lethal. Then, the transcriptome analysis of the mini mutant indicated that 2944 differentially expressed genes (DEGs) were identified from the silkworm in the 48 h of the second-instar, of which 1638 genes in the mini mutants were upregulated and 1306 genes were downregulated. These DEGs were mainly distributed in the biological process, cellular component, and molecular function. The functional annotation based on the KEGG database showed that these genes were mainly clustered in metabolic pathways, fatty acid metabolism pathways, ribosome biogenesis in eukaryotes, and so on. Further analysis indicated that some genes involved in the growth and metabolism including enzyme genes, juvenile hormone, and ecdysone exhibited different transcriptional levels. These results provided new experimental evidence regarding the mechanism of the underlying formation of mini mutants.

Estrogen-Related Receptor Influences the Hemolymph Glucose Content by Regulating Midgut Trehalase Gene Expression in the Last Instar Larvae of Bombyx mori

The expression of trehalase in the midgut of insects plays an important role in glucose supply to the hemolymph. Energy metabolism is usually regulated by the estrogen-related receptor (ERR). A decrease in ATP levels is caused by the ERR hindering glycolysis. However, the relationship between trehalose accumulation and ERR expression is still unclear. Here, we found that silkworm ERR (BmERR) is concentrated and BmERR expression is strongly correlated with trehalase in the midgut during the last instar silkworm larval stage. We cloned the promoter of the trehalase from Bombyx mori (BmTreh) and found that the ERR bound directly to the core response elements of the promoter. Cell level interference and the overexpression of ERR can reduce or enhance BmTreh transcription and promoter activity. Overexpressed transgenic BmERR can significantly increase the expression of BmTreh in the midgut of the last instar silkworm larvae, thereby hydrolyzing trehalose into glucose and releasing it into the hemolymph. Additionally, increased hemolymph glucose content reduces silkworm pupa weight but does not affect silk protein production from the silk gland. Our results suggest a novel function for BmERR through its involvement in BmTreh regulation and expand the understanding of ERR functions in insect trehalose metabolism.

Characterization of nuclear localization signal in Ostrinia furnacalis Masculinizer protein

Bombyx mori Masculinizer protein (BmMasc) is essential for both masculinization and dosage compensation in B. mori. We previously identified a bipartite nuclear localization signal (NLS) of BmMasc and two essential residues (lysine at 274 [K274] and arginine at 275 [R275]) implicated in its function. Sequence comparison showed the presence of putative NLSs in lepidopteran Masc proteins, but their functional properties and critical residues are unknown. Here we characterized a putative NLS of Ostrinia furnacalis Masc (OfMasc) using B. mori ovary-derived BmN-4 cell line. Deletion and alanine scanning mutagenesis revealed that a putative NLS is required for nuclear localization of OfMasc. However, mutations at both K227 and R228, which correspond to K274 and R275 of BmMasc, respectively, do not greatly abolish the NLS activity. Additional mutagenesis analysis revealed that triple mutations at K227, R228, and K240 almost completely inhibited OfMasc nuclear localization. These results suggest that lepidopteran Masc proteins possess a common functional NLS, but the critical residues for its activity are different. Moreover, we examined the masculinizing activity of OfMasc derivatives and found that nuclear localization is not required for the masculinizing activity of OfMasc. The results from our studies indicate that lepidopteran Masc proteins function in the cytoplasm to drive masculinizing cascade.

Scavenger receptor B8 improves survivability by mediating innate immunity in silkworm, Bombyx mori

Scavenger receptor class B (SR-B) is an extracellular transmembrane glycoprotein that plays a vital role in innate immunity. Although SR-Bs have been widely studied in vertebrates, their functions remained to elucidate in insects. Here, we identified and characterized a scavenger receptor class B member from the silkworm, Bombyx mori (designated as BmSCRB8). BmSCRB8 is broadly expressed in various immune tissues/organs, including fat body, gut, and hemocyte. Its expression is dramatically enhanced after challenge with different types of bacteria or pathogen-associated molecular patterns (PAMPs). The recombinant BmSCRB8 protein can detect different types of bacteria by directly binding to PAMPs and significantly improve the bacterial clearance in vivo. After knockdown of BmSCRB8, the pathogenic bacterial clearance was strongly impaired, and several AMP genes were down-regulated following E. coli challenge. Moreover, pathogenic bacteria's treatment following the depletion of BmSCRB8 remarkably decreased silkworm larvae's survival rate. Taken together, these results demonstrate that BmSCRB8 acts as a pattern recognition protein and plays an essential role in silkworm innate immunity by enhancing bacterial clearance and contributing to the production of AMPs in vivo.

Identification of long noncoding RNAs in silkworm larvae infected with Bombyx mori cypovirus

Bombyx mori cypovirus (BmCPV) is one of the most important pathogens causing severe disease to silkworm. Emerging evidence indicates that long noncoding RNAs (lncRNAs) play importantly regulatory roles in virus infection and host immune response. To better understand the interaction between silkworm, Bombyx mori and BmCPV, we performed a comparative transcriptome analysis on lncRNAs and mRNAs between the virus-infected and noninfected silkworm larvae midgut at two time points postinoculation. A total of 16,753 genes and 1845 candidate lncRNAs were identified, among which 356 messenger RNA (mRNAs) and 41 lncRNAs were differentially expressed (DE). Target gene prediction revealed that most of DEmRNAs (123) were coexpressed with 28 DElncRNAs, suggesting that the expression of mRNA is mainly affected through trans- regulation by BmCPV-induced lncRNAs, and a regulatory network of DElncRNAs and DEmRNAs was then constructed. According to the network, many genes involved in apoptosis, autophagy, and antiviral response, such as ATG3, PDCD6, IBP2, and MFB1, could be targeted by different DElncRNAs, implying the essential roles of these genes and lncRNAs in BmCPV infection. In all, our studies revealed for the first time the alteration of lncRNA expression in BmCPV-infected larvae and its potential influence on BmCPV replication, providing a new perspective for host-cypovirus interaction studies.

Bombyx mori nucleopolyhedrovirus downregulates transcription factor BmFoxO to elevate virus infection

Forkhead-box O (FoxO) is the primary transcriptional effector of the insulin-like signaling pathway that enhances gluconeogenesis through transcriptional activation of PEPCK and G6Pase in mammals. We have previously demonstrated the involvement of phosphoenolpyruvate carboxykinase (BmPEPCK-2) in antiviral immunity against the multiplication of Bombyx mori nuclearpolyhedrosisvirus (BmNPV) in silkworm. Therefore, we speculated that BmFoxO might suppress BmNPV by regulating the expression of PEPCK in silkworm. In the present study, we found that the expression of BmFoxO decreased after BmNPV infection in Bombyx mori; this finding was consistent with BmPEPCK-2 expression. In addition, the expression of BmFoxO was altered, and it was found that reduced expression of BmFoxO (dsBmFoxO) downregulated the expression of BmPEPCK-2 and increased the viral fluorescence and content in silkworm embryonic cell line BmE cells, and vice versa. BmFoxO could upregulate the expression of BmPEPCK-2 by binding to the BmPEPCK-2 promoter. Moreover, overexpression of BmFoxO significantly increased the expression of autophagy genes ATG6/7/8 after infection with BmNPV, consistent with BmPEPCK-2. These results indicate that BmNPV downregulates transcription factor BmFoxO to elevate virus infection, and BmFoxO overexpression upregulates BmPEPCK-2 expression and enhances silkworm antiviral resistance.

An update of KAIKObase, the silkworm genome database

KAIKObase was established in 2009 as the genome database of the domesticated silkworm Bombyx mori. It provides several gene sets and genetic maps as well as genome annotation obtained from the sequencing project of the International Silkworm Genome Consortium in 2008. KAIKObase has been used widely for silkworm and insect studies even though there are some erroneous predicted genes due to misassembly and gaps in the genome. In 2019, we released a new silkworm genome assembly, showing improvements in gap closure and covering more and longer gene models. Therefore, there is a need to include new genome and new gene models to KAIKObase. In this article, we present the updated contents of KAIKObase and the methods to generate, integrate and analyze the data sets. Database URL: https://kaikobase.dna.affrc.go.jp.

20E-mediated regulation of BmKr-h1 by BmKRP promotes oocyte maturation

BACKGROUND

Krüppel homolog 1 (Kr-h1) is a critical transcription factor for juvenile hormone (JH) signaling, known to play a key role in regulating metamorphosis and adult reproduction in insects. Kr-h1 can also be induced by molting hormone 20-hydroxyecdysone (20E), however, the underlying mechanism of 20E-induced Kr-h1 expression remains unclear. In the present study, we investigated the molecular mechanism of Kr-h1 induction by 20E in the reproductive system of a model lepidopteran insect, Bombyx mori.

RESULTS

Developmental and tissue-specific expression analysis revealed that BmKr-h1 was highly expressed in ovaries during the late pupal and adult stages and the expression was induced by 20E. RNA interference (RNAi)-mediated depletion of BmKr-h1 in female pupae severely repressed the transcription of vitellogenin receptor (VgR), resulting in the reduction in vitellogenin (Vg) deposition in oocytes. BmKr-h1 specifically bound the Kr-h1 binding site (KBS) between - 2818 and - 2805 nt upstream of BmVgR and enhanced the transcription of BmVgR. A 20E cis-regulatory element (CRE) was identified in the promoter of BmKr-h1 and functionally verified using luciferase reporter assay, EMSA and DNA-ChIP. Using pull-down assays, we identified a novel transcription factor B. mori Kr-h1 regulatory protein (BmKRP) that specifically bound the BmKr-h1 CRE and activated its transcription. CRISPR/Cas9-mediated knockout of BmKRP in female pupae suppressed the transcription of BmKr-h1 and BmVgR, resulting in arrested oogenesis.

CONCLUSION

We identified BmKRP as a new transcription factor mediating 20E regulation of B. mori oogenesis. Our data suggests that induction of BmKRP by 20E regulates BmKr-h1 expression, which in turn induces BmVgR expression to facilitate Vg uptake and oogenesis.

The Sex Determination Cascade in the Silkworm

In insects, sex determination pathways involve three levels of master regulators: primary signals, which determine the sex; executors, which control sex-specific differentiation of tissues and organs; and transducers, which link the primary signals to the executors. The primary signals differ widely among insect species. In Diptera alone, several unrelated primary sex determiners have been identified. However, the doublesex (dsx) gene is highly conserved as the executor component across multiple insect orders. The transducer level shows an intermediate level of conservation. In many, but not all examined insects, a key transducer role is performed by transformer (tra), which controls sex-specific splicing of dsx. In Lepidoptera, studies of sex determination have focused on the lepidopteran model species Bombyx mori (the silkworm). In B. mori, the primary signal of sex determination cascade starts from Fem, a female-specific PIWI-interacting RNA, and its targeting gene Masc, which is apparently specific to and conserved among Lepidoptera. Tra has not been found in Lepidoptera. Instead, the B. mori PSI protein binds directly to dsx pre-mRNA and regulates its alternative splicing to produce male- and female-specific transcripts. Despite this basic understanding of the molecular mechanisms underlying sex determination, the links among the primary signals, transducers and executors remain largely unknown in Lepidoptera. In this review, we focus on the latest findings regarding the functions and working mechanisms of genes involved in feminization and masculinization in Lepidoptera and discuss directions for future research of sex determination in the silkworm.

Small nucleolar RNA of silkworm can translocate from the nucleolus to the cytoplasm under abiotic stress

Small nucleolar RNAs (snoRNAs) are thought to be exclusively nuclear and guide nucleotide modifications of ribosomal RNAs. Recently, more and more evidence has suggested that the nucleolus is a stress sensor for changes in growth status and that snoRNAs may orchestrate the response to environmental stress through molecular interactions outside of the nucleus. We previously showed that a box C/D snoRNA Bm-15 had both nuclear and cytoplasmic location in BmN4 cell line of the silkworm, Bombyx mori. To further study the functional roles of Bm-15, changes in expression level and cellular location of Bm-15 were examined in BmN4 cells subjected to serum starvation and ultraviolet (UV) ray radiation. Results indicated that total RNA level of Bm-15 was unchanged after 24 h serum starvation, but exhibited 3-fold increases in the cytoplasm, and the nuclear-to-cytosolic distribution ratio was reduced from 5:1 to 2:1. Moreover, UV radiation also causes rapid decline in nuclear Bm-15 and progressive cytoplasmic accumulation with a percentage of 22% and 57% after 6 and 24 h UV radiation. UV treatment results in a dramatic decrease in Bm-15 nuclear-to-cytosolic ratio from 7:1 to 2:1 and 2:1 to 1:20 after 6 and 24 h UV radiation, respectively. We show here for the first time that box C/D snoRNAs can translocate from the nucleus to the cytoplasm under the abiotic stress of nutritional deficiency and UV radiation. The rapid translocation of snoRNAs from nucleus to cytoplasm may slow down the maturation of rRNAs and synthesis of ribosomes to enhance the stress resistance of cells.

The validation of the role of several genes related to Bombyx mori nucleopolyhedrovirus infection in vivo

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of primary silkworm pathogens and causes a serious damage of cocoon losses every year. Recent years, many works have been done to clarify the silkworm anti-BmNPV mechanism, and a significant progress has been made in screening and studying of genes and proteins related to BmNPV infection, but several of them lacked the proofs in vivo. In this study, to further validate the function of seven newly reported genes in vivo, including BmAtlatin-n, Bmferritin-heavy chain (BmFerHCH), Bmthymosin (BmTHY), Bmseroin1, Bmseroin2, Bmnuclear hormone receptors 96 (BmNHR96), and BmE3 ubiquitin-protein ligase SINA-like 10 (BmSINAL10), the response of them in the midgut, fat body, and hemolymph of differentially resistant strains (resistant strain YeA and susceptible strain YeB) at 48 h following BmNPV infection were analyzed. The results showed that the relative stable or upregulated expression level of BmAtlatin-n, BmTHY, Bmseroin1, and Bmseroin2 in YeA resistant strain following BmNPV infection further indicated their antiviral role in vivo, compared with susceptible YeB strain. Moreover, the significant downregulation of BmFerHCH, BmNHR96, and BmSINAL10 in both strains following BmNPV infection revealed their role in benefiting virus infection, as well as the upregulation of BmFerHCH in YeB midgut and BmSINAL10 in YeB hemolymph. These data could be used to complementary the proofs of the function of these genes in response to BmNPV infection.

Transcriptomic analysis reveals the regulation network of BmKrüppel homolog 1 in the oocyte development of Bombyx mori

Krüppel homolog 1 (Kr-h1), a zinc finger transcription factor, is involved in the metamorphosis and adult reproduction of insects. However, the role of Kr-h1 in reproduction of holometabolic insects remains to be elucidated. The regulation network of Kr-h1-associated genes in the reproduction in Bombyx mori was investigated in this study. The higher expression level of BmKr-h1 in the ovaries was detected during the late pupal stage and adults. RNA interference (RNAi)-mediated depletion of BmKr-h1 in the female at day 6 of pupae resulted in abnormal oocytes at 48 h post-double-stranded RNA treatment, which showed less yolk protein deposition and partially transparent chorion. RNA-seq and subsequent differentially expressed transcripts analysis showed that knockdown of BmKr-h1 caused a decrease in the expression of 2882 genes and an increase in the expression of 2565 genes in the oocytes at day 8 of pupae. Totally, 27 genes coding for transcription factors were down-regulated, while six genes coding for other transcription factors were up-regulated. BmKr-h1 bound to the Kr-h1 binding site of the transcription factors AP-1 (activating protein-1) and FOXG1 to increase their messenger RNA transcripts in the BmN cells, respectively. Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses of that positively co-expressed with AP-1 and FOXG1 transcripts showed mainly enrichment in the metabolic-related pathways, the nutrient absorption and the yolk protein absorption processes. These data suggested that BmKr-h1 might directly regulate the metabolic-related pathways, the nutrient absorption and the yolk protein absorption processes or probably through AP-1 and /or FOXG1 to regulate oocyte development.

Target antifungal peptides of immune signalling pathways in silkworm, Bombyx mori, against Beauveria bassiana

Antifungal innate immunity is an important defence used by insects against entomogenous fungi. However, the downstream target antifungal peptides of different immune signalling pathways are unknown. We found that the Toll, Janus kinase/signal transducer and activator of transcription (Jak/STAT) and Immunodeficiency (IMD) signalling pathways in the silkworm, Bombyx mori, can be activated by Beauveria bassiana. Inhibition of the Toll, IMD and Jak/STAT signalling pathways reduced the antifungal activities of silkworm haemolymph. We verified the target antifungal peptides of different immune signalling pathways. The expression patterns of five anti-fungal peptide genes in silkworm larvae and BmN cells were detected after blocking or over-expressing the immune signalling pathways. The Toll signalling pathways mediated the expression of Bmcecropin A, Bmattacin 1 and Bmgloverin 2; IMD signalling pathways mediated Bmenbocin 1, Bmgloverin 2 and Bmattacin 1; Jak/STAT signalling pathways mediated Bmstorage protein 30K-19G1 (Bmsp 1), Bmattacin 1 and Bmcecropin A. These data indicated that anti-microbial peptide genes in B. mori evolved through expansion and selection of existing genes to adapt to the challenge of invasive microorganisms such as fungi. This information provides insight into the antifungal immune responses in B. mori and aids understanding of insect immune regulation mechanisms.

The biological characters of Bmelav-like genes in the development of Bombyx mori

The ELAV/Hu family is a conserved multigene family of pan-neuronal RNA-binding protein involved in post-transcriptional regulation in metazoans. In Drosophila, three members of this family, ELAV, RBP9 and FNE, are involved in neuronal differentiation, gene expression regulation and so on. This family is less well characterized in Bombyx mori. Two orthologs BmELAV-like-1 (BmEL-1) and BmELAV-like-2 (BmEL-2) share 55%-71% and 47%-62% identity with that of in Drosophila and humans, respectively. Bmel-1 is ubiquitously expressed while Bmel-2 is expressed in the head and ovaries specifically. Proteins encoded by both genes are localized in nuclear and cytoplasm. The weight of body, cocoon, pupae and cocoon shell are differently affected in Bmel-1^- /-2^- mutants created using CRISPR/Cas9 technology. Mutations of both genes increase the expression of four silk protein genes, Fib-L, Fib-H, P25 and Ser-1. In addition, the oviposition ability of Bmel-2^- females is decreased. This study not only provides valuable insights into the functional roles of Bmelav-like genes in the growth, cocoon characters and regulation of silk protein genes expression, but also provides useful information for silkworm variety breeding.

Peptidoglycan recognition protein-S1 acts as a receptor to activate AMP expression through the IMD pathway in the silkworm Bombyx mori

Peptidoglycan recognition proteins (PGRPs) are the most important pattern recognition receptors (PRRs) in insects. PGRPs can recognize pathogenic microorganism peptidoglycans (PGs) and play an important role in innate immunity. Twelve PGRPs have been identified in silkworms. However, the specific roles played by these PGPRs in the silkworm innate immune system have not been elucidated to date. In this study, we systematically investigated the biological functions of BmPGRP-S1 in silkworms. We observed that BmPGRP-S1 was highly expressed in silkworm immune-related organs and was upregulated in response to bacterial challenges. Furthermore, we determined that BmPGRP-S1 can bind to bacteria or PGs and activate antimicrobial peptide (AMP) expression. Inhibition of the expression of BmPGRP-S1 by siRNA reduced AMP gene expression in silkworms. Further experiments demonstrated that BmPGRP-S1 is involved in IMD pathway activation to induce AMP expression. Taken together, these results demonstrate that BmPGRP-S1 serves as a receptor to activate AMP gene expression through the IMD pathway to address bacterial challenges.

Genetic bioengineering of overexpressed guanylate binding protein family BmAtlastin-n enhances silkworm resistance to Nosema bombycis

Microsporidia are obligate single-celled eukaryote parasites. Microsporidian infection can cause large economic losses to beneficial insects such as silkworms and honey bees. Identification of resistance biomacromolecules and breeding of transgenic lines resistant to the microsporidian Nosema bombycis are important for disease management. We previously used transcriptome analysis to identify a guanylate binding protein family BmAtlastin-n gene that was significantly upregulated after Nosema bombycis infection, and we determined that the molecule was highly expressed in resistance-related tissues such as the midgut, fat body and the epidermis. The transgenic silkworm line overexpressing BmAtlastin-n biomolecules had economic characters similar to those of non-transgenic lines. The transgenic OE-BmAtlastin-n lines had significantly improved survival after microspore infection. We used RT-PCR and H&E staining to show that the number of spores in the transgenic lines was significantly lower than in the control lines. In this study, we identified a BmAtlastin-n macromolecule with resistance to N. bombycis and developed a transgenic line. The results improved understanding of the GBP protein family and provided biomacromolecule material for the treatment and prevention of microsporidia.

Maternal GABAergic and GnRH/corazonin pathway modulates egg diapause phenotype of the silkworm Bombyx mori

Diapause represents a major developmental switch in insects and is a seasonal adaptation that evolved as a specific subtype of dormancy in most insect species to ensure survival under unfavorable environmental conditions and synchronize populations. However, the hierarchical relationship of the molecular mechanisms involved in the perception of environmental signals to integration in morphological, physiological, behavioral, and reproductive responses remains unclear. In the bivoltine strain of the silkworm Bombyx mori, embryonic diapause is induced transgenerationally as a maternal effect. Progeny diapause is determined by the environmental temperature during embryonic development of the mother. Here, we show that the hierarchical pathway consists of a γ-aminobutyric acid (GABA)ergic and corazonin signaling system modulating progeny diapause induction via diapause hormone release, which may be finely tuned by the temperature-dependent expression of plasma membrane GABA transporter. Furthermore, this signaling pathway possesses similar features to the gonadotropin-releasing hormone (GnRH) signaling system for seasonal reproductive plasticity in vertebrates.

Functional conservation and diversification of yellow-y in lepidopteran insects

The diverse colors and patterns found in Lepidoptera are important for success of these species. Similar to the wings of adult butterflies, lepidopteran larvae exhibit diverse color variations to adapt to their habitats. Compared with butterfly wings, however, less attention has been paid to larval body colorations and patterns. In the present study, we focus on the yellow-y gene, which participates in the melanin synthesis pathway. We conducted CRISPR/Cas9-mediated targeted mutagenesis of yellow-y in the tobacco cutworm Spodoptera litura. We analyzed the role of S. litura yellow-y in pigmentation by morphological observation and discovered that yellow-y is necessary for normal black pigmentation in S. litura. We also showed species- and tissue-specific requirements of yellow-y in pigmentation in comparison with those of Bombyx mori yellow-y mutants. Furthermore, we found that almost none of the yellow-y mutant embryos hatched unaided. We provide evidence that S. litura yellow-y has a novel important function in egg hatching, in addition to pigmentation. The present study will enable a greater understanding of the functions and diversification of the yellow-y gene in insects.

D-3-phosphoglycerate dehydrogenase from the silkworm Bombyx mori: Identification, functional characterization, and expression

D-3-phosphoglycerate dehydrogenase (PHGDH) is a key enzyme involved in the synthesis of l-serine. Despite the high serine content in silk proteins and the crucial role of PHGDH in serine biosynthesis, PHGDH has not been described in silkworms to date. Here, we identified PHGDH in the silkworm Bombyx mori and evaluated its biochemical properties. On the basis of the amino acid sequence and phylogenetic tree, this PHGDH has been categorized as a new type and designated as bmPHGDH. The recombinant bmPHGDH was overexpressed and purified to homogeneity. Kinetic studies revealed that PHGDH uses NADH as a coenzyme to reduce phosphohydroxypyruvate. High expression levels of bmphgdh messenger RNA (mRNA) were observed in the middle part of the silk gland and midgut in a standard strain of silkworm. Moreover, a sericin-deficient silkworm strain displayed reduced expression of bmphgdh mRNA. These findings indicate that bmPHGDH might play a crucial role in the provision of l-serine in the larva of B. mori.

Transcriptome analysis of Nosema assamensis infecting muga silkworms (Antheraea assamensis) reveals insights into candidate pathogenicity related genes and molecular pathways required for pathogenesis

Muga silkworms are often prone to many diseases since, these are non-domesticated and are reared outdoors. Microsporidia, an obligate intracellular pathogen with spore as its active form, causes pebrine disease in these silkworms. The study has attempted to categorise the transcript data of the Nosema obtained from the infected muga silkworm using gene ontology and KEGG pathway studies. A total of 2850 unigene sets were identified out of which 2739 unigenes were placed under biological, cellular as well as molecular function categories based on the gene ontology (GO) terms. 1620 out of these unigenes sets found their orthologous partner in the corresponding Nosema bombycis transcriptome. The unigenes were found to be enriched under organic substance metabolic process, organic cyclic compound binding and intracellular anatomical structure for biological process, molecular function and cellular components respectively. The KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis indicated majority of the enzymes were found to be mapped under purine and thiamine metabolic pathways, indicating an increase in the energy metabolism required to establish the infection in the silkworms. The putative virulence genes identified in this study are PTP2, PTP3, SWP12 and SWP26 which were found to be expressed in other Nosema species indigenous to India, indicating a probable conservation of these genes, which are primarily involved in establishing host pathogen interactions. The expression of these genes was in detectable levels in the infected silkworm samples. These genes may be validated further through bioassay in order understand their roles in establishing the infection and propagation of the spores. The identified virulence genes may be further targeted to develop diagnostic tools for identification of the pathogen at early stages of infection.

Molecular cloning, expression, and characterization of BmSOD3 in silkworm (Bombyx mori)

Superoxide dismutases (SODs) play an essential role in eliminating excess reactive oxygen species and maintaining the redox balance of the immune system. To study the function of BmSOD3 in silkworm, 543-bp full-length complementary DNA-encoding BmSOD3 was cloned from silkworm. The BmSOD3 amino acids were compared to their homologs, and several highly conserved regions were analyzed. We also carried out phylogenetic analyses of the SOD gene. Our results showed that the BmSOD3 gene belonged with the ecCu/Zn SOD gene. The BmSOD3 gene was transformed into the pET28a vector for functional expression in Escherichia coli. The sodium salt-polyacrylamide gel electrophoresis results showed that the molecular weight of recombinant BmSOD3 was about 22 kDa. The recombinant protein BmSOD3 was purified to detect its properties. After purification analyses, the enzyme activity showed Cu/Zn SOD activity, and the specific activity of the purified enzyme was 0.51 U/mg. The BmSOD3 transcripts showed tissue-specific expression in the midgut and malpighian tubule. The immune microarray data for BmSOD3 showed an expression signal that had a strong response to the induction of four pathogens (Bacillus bombyseptieus, Beauveria bassiana, E. coli, and nuclear polyhedrosis virus), particularly after infection for 24 h, which indicates that the BmSOD3 gene plays a key role in response to bacterial, fungal, and viral invasion. The fusion protein also showed antibacterial activity against E. coli in vitro. Thus, the fusion protein BmSOD3 exhibits antibacterial activity and may be used in production to combat diseases caused by bacteria in silkworm.

Transgenic expression of antimicrobial peptides from black soldier fly enhance resistance against entomopathogenic bacteria in the silkworm, Bombyx mori

Antimicrobial peptides (AMPs) are effective molecules produced by the innate immune system of most organisms to fend off invading microbes and regarded as promising alternatives to conventional antibiotics due to their potent antimicrobial activities. The larvae of black soldier fly (BSF), Hermetia illucens, inhabit microbe-rich environments and its insect genome encodes a broad repertoire of AMPs. In the present study, three AMPs encoded by BSF Hidefensin-1, Hidiptericin-1 and HiCG13551 were cloned, expressed and purified in a recombinant Escherichia coli expression system. In vitro, both Hidefensin-1 and Hidiptericin-1 inhibited the growth of Streptococcus pneumoniae and Escherichia coli, while HiCG13551 inhibited the growth of Staphylococcus aureus and E. coli. Transmission electron microscopy showed that Hidiptericin-1 inhibited bacterial growth through bacterial membrane lysis. We also constructed a transgenic silkworm line constitutively expressing an AMP cassette HiAMP4516 encoding all the three AMPs, and the silkworms showed an increased resistance to both gram-positive and gram-negative entomopathogenic bacteria. These results provide insights into the antibacterial activities of BSF AMPs both in vitro and in vivo and suggest a great potential of exploiting insect-derived AMPs in silkworm disease resistance breeding.

Structural insight into the substrate specificity of Bombyx mori β-fructofuranosidase belonging to the glycoside hydrolase family 32

Sucrose-hydrolyzing enzymes are largely divided into β-fructofuranosidase and sucrose α-glucosidase. The domestic silkworm Bombyx mori possesses both enzymes, BmSUC1 and BmSUH, belonging to the glycoside hydrolase family 32 (GH32) and GH13, respectively. BmSUC1 was presumed to be acquired by horizontal gene transfer from bacteria based on phylogenetic analysis and related to tolerance to sugar-mimic alkaloids contained in mulberry latex. Here we investigated the substrate specificity of recombinant BmSUC1 that can hydrolyze not only sucrose but also fructooligosaccharides and fructans, and revealed that the enzyme was competitively inhibited by 1,4-dideoxy-1,4-imino-D-arabinitol, one of the alkaloids. Moreover, the crystal structures of BmSUC1 in apo form and complex with sucrose were determined, and the active site pocket was shallow and suitable for shorter substrates but was related to more relaxed substrate specificity than the strict sucrose α-glucosidase BmSUH. Considering together with the distribution of BmSUC1-orthologous genes in many lepidopterans, our results suggest that BmSUC1 contributes to the digestion of fructooligosaccharides and fructans derived from feed plants.

Signal transducer and activator of transcription is involved in the expression regulation of ecdysteroid-induced insulin-like growth factor-like peptide in the pupal wing disc of silkworm, Bombyx mori

In insects, 20-hydroxyecdysone (20E) and insulin-like growth factor-like peptides (IGFLPs) regulate the development of imaginal discs. However, how IGFLPs are up-regulated to impact the development of the pupal wing disc is still unclear. In this study, we investigated the expression regulation of IGFLP in the pupal wing disc of silkworm, Bombyx mori. We confirmed that B. mori IGFLP (BmIGFLP) was mainly expressed in the pupal wing disc and the expression of BmIGFLP could be significantly induced by 20E. Bioinformatics analysis of BmIGFLP promoter sequence revealed three cis-regulation elements (CREs) of signal transducer and activator of transcription (STAT), which is a key component in the Janus-activated kinase / STAT pathway. Luciferase activity assays showed that two CREs enhanced the transcriptional activity of BmIGFLP. Electrophoretic mobility shift and chromatin immunoprecipitation assays demonstrated that BmSTAT proteins in the nuclear extracts of B. mori pupal wing discs and BmN cells could only bind to the STAT CRE3, indicating that STAT CRE3 activated by BmSTAT enhances BmIGFLP expression at pupal stages. Although 20E could not enhance the expression of BmSTAT, 20E enhanced the nucleus translocation of BmSTAT to bind with the STAT CRE3 in the BmIGFLP promoter. The increase of transcriptional activity of the STAT CRE3 by overexpression of BmSTAT and addition of 20E in BmN cells confirmed this result. Taken together, all data indicate that BmSTAT is one of the transcription factors activating 20E-induced BmIGFLP expression in the pupal wing disc.

The mechanism of damage to the posterior silk gland by trace amounts of acetamiprid in the silkworm, Bombyx mori

Acetamiprid is a new neonicotinoid insecticide widely used in the prevention and control of pests in agriculture. However, its residues in the environment affect the cocooning of the silkworm, Bombyx mori (B. mori), a non-target insect. To investigate the mechanism of damage, B. mori larvae were fed with trace amounts of acetamiprid (0.15 mg/L). At 96 h after exposure, the larvae showed signs of poisoning and decreased body weight, resulting in reduced survival and ratio of cocoon shell. At 48 h and 96 h after exposure, the residues in the posterior silk gland (PSG), which is responsible for synthesizing silk fibroin, were 0.72 μg/mg and 1.21 μg/mg, respectively, as measured by high performance liquid chromatography, indicating that acetamiprid can accumulate in the PSG. Moreover, pathological sections and transmission electron microscopy also demonstrate the damage of the PSG by acetamiprid. Digital gene expression (DGE) and KEGG pathway enrichment analysis revealed that genes related to metabolism, stress responses and inflammation were significantly up-regulated after exposure. Quantitative RT-PCR analysis showed that the transcript levels of FMBP-1 and FTZ-F1 (transcription factors for synthesizing silk protein) were up-regulated by 2.55-and 1.56-fold, respectively, and the transcript levels of fibroin heavy chain (Fib-H), fibroin light chain (Fib-L), P25, Bmsage and Bmdimm were down-regulated by 0.75-, 0.76-, 0.65-, 0.44- and 0.40-fold, respectively. The results indicate that accumulated acetamiprid causes damage to the PSG and leads to reduced expression of genes responsible for synthesizing silk fibroin. Our data provide reference for evaluating the safety of acetamiprid residues in the environment for non-target insects.

Induction of ER stress, antioxidant and detoxification response by sublethal doses of chlorantraniliprole in the silk gland of silkworm, Bombyx mori

Sublethal doses of chlorantraniliprole (CAP) disrupt spinning disorder in the silkworm Bombyx mori (B. mori) and cause reduced cocoon production. In the present study, we investigated the effects of trace amounts of CAP on morphology and gene expression of the B. mori silk gland, found the posterior silk gland cells were possessed of disintegrated Endoplasmic reticulum (ER), unevenly distributed chromatin after exposure to CAP (0.01 mg/L). Gene expression analysis revealed that IRE1 and ATF6 ER stress-signaling pathways were inhibited, the PERK/CncC pathway was activated. Digital gene expression (DGE) analysis showed that detoxification-related genes, antioxidant genes and genes involved in ER protein processing pathway were expressed differentially in CAP-treated silkworm larvae. Notably, the transcript levels of the detoxification-related genes (CYP4M5, CYP6AB4, GSTD3 and GSTS1) and the antioxidant genes (CAT, TPX and SOD) were significantly increased, and the expression of ER protein processing-related genes (Sec61β, Sec61γ, Sec23α and ERGIC-53) was significantly decreased after CAP exposure. The results showed that sublethal doses of CAP exposure caused ER stress, oxidative damage to the silk gland and the perturbation of protein processing in ER, thereby probably leading to abnormal growth of the silk gland and triggering the spinning failure in silkworm.

Expression profile analysis of circular RNAs in BmN cells (Bombyx mori) upon BmNPV infection

The disease caused by Bombyx mori nucleopolyhedrovirus (BmNPV) has always been difficult to control, resulting in tremendous economic losses in the sericulture industry. Although much has been learned about the impact of noncoding RNAs on pathogenesis, the role of circular RNA (circRNA) in insect immunity remains unclear. To explore circRNA regulation involved in BmNPV infection, we used transcriptome analysis of BmN cells with or without BmNPV infection to generate circRNA data set. A total of 444 novel circRNAs were identified in BmN cells, with 198 pervasively distributed both in the control group and BmNPV-infection group. The host genes were enriched inMAPK signaling pathway, dorso-ventral axis formation, and ECM-receptor interaction, which were required for the normal larval growth. A total of 75 circRNAs were differentially expressed (DE) on BmNPV infection. Six downregulated circRNAs were validated by Sanger sequencing and qRT-PCR. DEcircRNA-miRNA-DEmRNA network was constructed based on the six validated circRNAs. Pathway analysis indicated that the predicted target genes were mainly enriched in the metabolic pathway and immune-related signaling pathway. Our results may provide a basis for further studies on circRNA function in BmN cells challenged by BmNPV infection and offer an insight into the molecular mechanism on silkworm-virus interaction.

Identification and effect of Zf-AD-containing C2H2 zinc finger genes on BmNPV replication in the silkworm (Bombyx mori)

Zf-AD-containing C2H2 zinc -finger genes (ZAD) are uniquely present and have lineage-specific expansion in arthropods. Arthropods are also the hosts of Baculoviruses. We studied the possible relationship between the lineage-specific expansion of ZAD genes and arthropod-Baculovirus co-evolution. We used the silkworm (Bombyx mori) as a model. We identified 73 ZAD genes (BmZAD) in the silkworm. Sequence-based similarity analysis showed that nine clusters involving 28 BmZADs may have undergone species-specific expansion in the silkworm. Expression pattern analysis showed that the BmZADs were divided into five groups. Group I comprised 10 genes with high expression in multiple tissues, suggesting that BmZADs may play roles in the development of various tissues. We identified six BmZADs that could be induced by the Nucleopolyhedrovirus (BmNPV). Among them, BmZAD69 expression is capable of responding to BmNPV infection, and the ZAD domain is indispensable for the function of BmZAD69 in BmNPV replication. We also detected a 3 bp deletion at 1.7 kb upstream of BmZAD69, which may make it more sensitive to BmNPV infection, and thus elevate the BmNPV resistance in Qiufeng_N, a strain with strong virus resistance. These data suggest that BmZADs may be involved in BmNPV infection and that ZAD genes may play a role in arthropod-Baculovirus co-evolution.

A defect in purine nucleotide metabolism in the silkworm, Bombyx mori, causes a translucent larval integument and male infertility

p-oily (op) is a novel mutant of Bombyx mori exhibiting translucent larval integument and male infertility. Elucidation of the causative gene of the op mutant will help understand the genetic mechanism underlying larval integument coloration and male fertility. Using polymorphisms between B. mori and B. mandarina, the op locus was narrowed down to a 375-kb region. Using RNA-seq analysis, we found that op mutants have a frameshift mutation in the KWMTBOMO13770 gene located in the 375-kb region. A database search indicated that this gene is the human cytosolic 5'-nucleotidase II gene (cN-II) homolog in Bombyx, which mediates the conversion of inosine monophosphate (IMP) to inosine, a precursor of uric acid. CRISPR/Cas9-mediated knockout mutants of the Bm-cN-II gene showed translucent integuments, and there appeared translucent larvae in the crosses between knockout moths and +/op moths. Moreover, the translucent phenotype of, and decreased uric acid content in the larval integument caused by the mutations in the Bm-cN-II gene were rescued by oral administration of inosine. These results indicated that the Bm-cN-II gene is responsible for the op phenotype and that the molecular function of the Bm-cN-II gene is the conversion of IMP to inosine. We also discuss the genetic relationship between the Bm-cN-II gene and male fertility.

Construction and analysis of the protein-protein interaction network for the olfactory system of the silkworm Bombyx mori

Olfaction plays an essential role in feeding and information exchange in insects. Previous studies on the olfaction of silkworms have provided a wealth of information about genes and proteins, yet, most studies have only focused on a single gene or protein related to the insect's olfaction. The aim of the current study is to determine key proteins in the olfactory system of the silkworm, and further understand protein-protein interactions (PPIs) in the olfactory system of Lepidoptera. To achieve this goal, we integrated Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and network analyses. Furthermore, we selected 585 olfactory-related proteins and constructed a (PPI) network for the olfactory system of the silkworm. Network analysis led to the identification of several key proteins, including GSTz1, LOC733095, BGIBMGA002169-TA, BGIBMGA010939-TA, GSTs2, GSTd2, Or-2, and BGIBMGA013255-TA. A comprehensive evaluation of the proteins showed that glutathione S-transferases (GSTs) had the highest ranking. GSTs also had the highest enrichment levels in GO and KEGG. In conclusion, our analysis showed that key nodes in the biological network had a significant impact on the network, and the key proteins identified via network analysis could serve as new research targets to determine their functions in olfaction.

Role of Bm30kc6 gene in cell apoptosis and the silk gland degradation signaling pathway in Bombyx mori L

Apoptosis is a process of programmed cell death that is regulated by genes independently. The Bm30kc6 gene is a kind of small molecular lipoprotein about 30 kDa, expressed highly in the late stage of the silkworm hemolymph. Our study showed that overexpression of Bm30kc6 could decrease caspase-3 activation. Meanwhile, activation of caspase-3 increased when Bm30kc6 expression was disturbed by small interfering RNA (siRNA). Cell apoptosis was decreased when Bm30kc6 was overexpressed under UV treatment. The apoptosis rate induced by actinomycin D is similar to the trend by UV. It was inferred that Bm30kc6 has an inhibitory effect on the apoptosis of silkworm cells. The apoptosis-related genes, such as BmFadd, BmDredd, and BmDaxx were increased after overexpression of Bm30kc6 or decreased after interference of siRNA. It was speculated that there was an interactive relationship between Bm30kc6, BmDaxx, BmFadd, and BmDredd in the apoptosis signaling pathways. We investigated the transcription expression of the Bm30kc6 gene in different growth stages and tissues of the silkworm. The results showed that Bm30kc6 reached its peak in the hemolymph during the 6th to 7th days of the 5th instar, or in spinning post 24 h of the silk gland. In the silkworm BmN cells treated with caspase-3/7 inhibitor, the caspase-3 enzyme activity, and the expression levels of Bm30kc6, BmFadd, BmDredd, and BmDaxx were significantly reduced. The expression levels of Bm30kc6 increased sharply when silkworms were treated by molting hormone at Day 3 or 5 of the 5th instar. The results indicated that the expression of the Bm30kc6 gene was affected by the molting hormone and was likely to be its downstream target. In conclusion, the results suggest that the Bm30kc6 gene is involved in the regulation of the apoptotic signaling pathway and plays a role in the apoptotic process.

Effects of Quercetin on the Growth and Expression of Immune-Pathway-Related Genes in Silkworm (Lepidoptera: Bombycidae)

Quercetin is a flavonoid produced as a defense by plants. The effects of 1% quercetin on the growth and development of Bombyx mori were studied. The activities of the enzymes superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), carboxy-lesterase (CarEs), and glutathione S-transferase (GST) were all measured at 24, 48, 72, and 96 h after quercetin exposure. The results show that quercetin induces the activities of antioxidant and detoxification enzymes. With longer exposure times, enzyme activity first increased and then decreased. The relative expressions of AMP (defensin, CecA), the Toll pathway (cactus, Spatzle, and Rel), the IMD pathway (Imd, Fadd, and Dorsal), the JAK-STAT pathway (STAT, HOP, and Pi3k60), and the Melanization gene (DDC and PAH) were analyzed using quantitative polymerase chain reaction (qPCR). The results indicated that long-term exposure to quercetin could inhibit the expression of immune-related pathway genes in silkworms. This suggests that it can inhibit the activities of antioxidant and detoxifying enzymes, thus inhibiting the immune system and affecting the growth and development, resulting in an increase in the death rate in silkworm. This study provides the novel conclusion that quercetin accumulation inhibits the immune system of silkworm and increases its death rate, a result that may promote the development and utilization of better biopesticides that avoid environmental pollution.

Identification and Characterization of Aldehyde Oxidase 5 in the Pheromone Gland of the Silkworm (Lepidoptera: Bombycidae)

Aldehyde oxidases (AOXs) are a subfamily of cytosolic molybdo-flavoenzymes that play critical roles in the detoxification and degradation of chemicals. Active AOXs, such as AOX1 and AOX2, have been identified and functionally analyzed in insect antennae but are rarely reported in other tissues. This is the first study to isolate and characterize the cDNA that encodes aldehyde oxidase 5 (BmAOX5) in the pheromone gland (PG) of the silkworm, Bombyx mori. The size of BmAOX5 cDNA is 3,741 nucleotides and includes an open reading frame, which encodes a protein of 1,246 amino acid residues. The theoretical molecular weight and isoelectric point of BmAOX5 are approximately 138 kDa and 5.58, respectively. BmAOX5 shares a similar primary structure with BmAOX1 and BmAOX2, containing two [2Fe-2S] redox centers, a FAD-binding domain, and a molybdenum cofactor (MoCo)-binding domain. RT-PCR revealed BmAOX5 to be particularly highly expressed in the PG (including ovipositor) of the female silkworm moth, and the expression was further confirmed by in situ hybridization, AOX activity staining, and anti-BmAOX5 western blotting. Further, BmAOX5 was shown to metabolize aromatic aldehydes, such as benzaldehyde, salicylaldehyde, and vanillic aldehyde, and fatty aldehydes, such as heptaldehyde and propionaldehyde. The maximum reaction rate (Vmax) of benzaldehyde as substrate was 21 mU and Km was 1.745 mmol/liter. These results suggested that BmAOX5 in the PG could metabolize aldehydes in the cytoplasm for detoxification or participate in the degradation of aldehyde pheromone substances and odorant compounds to identify mating partners and locate suitable spawning sites.

Formation of the CenH3-Deficient Holocentromere in Lepidoptera Avoids Active Chromatin

Despite the essentiality for faithful chromosome segregation, centromere architectures are diverse among eukaryotes^1,2 and embody two main configurations: mono- and holocentromeres, referring, respectively, to localized or unrestricted distribution of centromeric activity. Of the two, some holocentromeres offer the curious condition of having arisen independently in multiple insects, most of which have lost the otherwise essential centromere-specifying factor CenH3³ (first described as CENP-A in humans).^4-7 The loss of CenH3 raises intuitive questions about how holocentromeres are organized and regulated in CenH3-lacking insects. Here, we report the first chromatin-level description of CenH3-deficient holocentromeres by leveraging recently identified centromere components^6,7 and genomics approaches to map and characterize the holocentromeres of the silk moth Bombyx mori, a representative lepidopteran insect lacking CenH3. This uncovered a robust correlation between the distribution of centromere sites and regions of low chromatin activity along B. mori chromosomes. Transcriptional perturbation experiments recapitulated the exclusion of B. mori centromeres from active chromatin. Based on reciprocal centromere occupancy patterns observed along differentially expressed orthologous genes of Lepidoptera, we further found that holocentromere formation in a manner that is recessive to chromatin dynamics is evolutionarily conserved. Our results help us discuss the plasticity of centromeres in the context of a role for the chromosome-wide chromatin landscape in conferring centromere identity rather than the presence of CenH3. Given the co-occurrence of CenH3 loss and holocentricity in insects,⁷ we further propose that the evolutionary establishment of holocentromeres in insects was facilitated through the loss of a CenH3-specified centromere.

Identification of Peritrophins and Antiviral Effect of Bm01504 against BmNPV in the Silkworm, Bombyx mori

The insect midgut secretes a semi-permeable, acellular peritrophic membrane (PM) that maintains intestinal structure, promotes digestion, and protects the midgut from food particles and pathogenic microorganisms. Peritrophin is an important PM protein (PMP) in the PM. Here, we identified 11 peritrophins with 1–16 chitin binding domains (CBDs) comprising 50–56 amino acid residues. Multiple CBDs in the same peritrophin clustered together, rather than by species. The CBD contained six highly conserved cysteine residues, with the key feature of amino acids between them being CX_11-15CX₅CX_9-14CX_11-12CX_6-7C. Peritrophins with 2 and 4 CBDs (Bm09641 and Bm01504, respectively), and with 1, 8, and 16 CBDs (Bm11851, Bm00185, and Bm01491, respectively) were mainly expressed in the anterior midgut, and throughout the midgut, respectively. Survival rates of transgenic silkworms with Bm01504 overexpression (Bm01504-OE) and knockout (Bm01504-KO) infected with B. morinucleopolyhedrovirus (BmNPV) were significantly higher and lower, whereas expression of the key viral gene, p10, were lower and higher, respectively, compared with wild type (WT). Therefore, Bm01504-OE and Bm01504-KO transgenic silkworms were more and less resistant, respectively, to BmNPV. Bm01504 plays important roles in resisting BmNPV invasion. We provide a new perspective for studying PM function, and reveal how the silkworm midgut resists invasive exogenous pathogenic microorganisms.

Semaphorin-1a-like gene plays an important role in the embryonic development of silkworm, Bombyx mori

Fuyin-lethal red egg (Fuyin-lre) is a red egg mutant discovered from the germplasm resource Fuyin of Bombyx mori. The embryo of Fuyin-lre stops developing at the late stage of gastrulation due to chromosome structural variation. In this work, precise mutation sites at both ends of the mutated region were determined, and two inserted sequences with lengths of 1232 bp and 1845 bp were obtained at both ends of the mutation region. Interestingly, a bmmar1 transposon was detected in the inserted 1845 bp sequence. Bmmar1 possesses features of the Tcl/mariner superfamily of transposable elements (TEs), which belongs to class II TEs that use a DNA-mediated "cut and paste" mechanism to transpose. This finding suggests that Fuyin-lre mutation might be related to the "cut and paste" action of bmmar1. The mutation resulted in the deletion of 9 genes in the mutation region, of which the red egg gene re (BMSK0002766) did not affect embryonic development of B. mori, and the BMSK0002765 gene was unexpressed during the early stage of embryonic development. The RNA interference results of the remaining 7 genes suggest that the semaphorin-1a-like gene (BMSK0002764) had a major contribution to the embryonic lethality of Fuyin-lre.

Enhanced silk yield in transgenic silkworm (Bombyx mori) via ectopic expression of BmGT1-L in the posterior silk gland

The silkworm is an economically important insect producing plentiful silk fibre in the silk gland. In this study, we reported a cross-talk between the fat body, silk gland and midgut through a glycine-serine biosynthetic pathway in the silkworm. Amino acid sequence and functional domains of glycine transporter gene BmGT1-L were mapped. Our results indicated that BmGT1-L was specifically expressed in the midgut microvilli and persistently expressed during the feeding stages. RNA interference of BmGT1-L activated glycine biosynthesis, and BmGT1-L overexpression facilitated serine biosynthesis in the BmN4-SID1 cell. In addition, silkworms after FibH gene knock-out or silk gland extirpation showed markedly decreased BmGT1-L transcripts in the midgut and disturbed glycine-serine biosynthesis as silk yield decreased. Finally, BmGT1-L ectopic expression in the posterior silk gland promoted glycine biosynthesis, and enhanced silk yield via increasing fibroin synthesis. These results suggested that cross-talk between tissues can be used for enhancing silk yield in the silkworm.

A novel system for separation and purification of egg pigments from the nondiapause red-egg mutant Re-nd in the silkworm, Bombyx mori

Re-nd, which was induced from the wild-type C108 by the chemical mutagen N-methane-N-methylnitrourea, is a nondiapause red-egg mutant of silkworm Bombyx mori. The special significance of the Re-nd mutant is that it is an independent dominant mutant. The aim of this study was to establish the type of pigment responsible for the red coloration in the Re-nd mutant eggs in silkworm. We compared the eggs of Re-nd mutants with those of the other B. mori egg color strains and confirmed that the Re-nd mutant is the only strain with red color and red pigment granules in nondiapause, showing this mutant belongs to the pigmentation in the serosa. We speculated that the red substance, which contributed to the bright red pigmentation for nondiapause eggs of the Re-nd mutant, could potentially be a novel pigment according to its solubility, optimum absorption peak, and oxidation-reduction reaction. Moreover, we have successfully constituted the system for enrichment, extraction, and purification of the red substance responsible for the Re-nd mutant, providing a new method for the separation and purification of other known and unknown pigments or substances.

Functional characterization of an aldose reductase (bmALD1) obtained from the silkworm Bombyx mori

We describe a new member of the aldo-keto reductase (AKR) superfamily in the silkworm Bombyx mori. On the basis of its amino acid sequence and phylogenetic tree, this AKR belongs to the AKR1B family and has been designated as bmALD1. In the current study, recombinant bmALD1 was overexpressed, purified to homogeneity and kinetically characterized. We discovered that bmALD1 uses NADPH as a coenzyme to reduce carbonyl compounds such as DL-glyceraldehyde, glucose and 2-nonenal. No NADH-dependent activity was detected. To the best of our knowledge, bmALD1 is only the third AKR characterized in silkworm which, given its substrate specificity, could play a major role in glucose metabolism and antioxidant reactions. Our data provide an increased understanding of insect AKR function.

Effect of pyriproxyfen exposure on cocooning and gene expression in the silk gland of Bombyx mori （Linnaeus, 1758）

Bombyx mori（Linnaeus, 1758） is an important economical insect, and the sericulture is a flourishing industry in many developing countries. Pyriproxyfen, a juvenile hormone pesticide, is often applied to cultivations widely in the world, and its exposure often resulted in silk yield reduction and non-cocooning. However, the effect of pyriproxyfen exposure on cocooning and gene expression level in the silk gland of B. mori has not been studied yet, and this study focused on the above issues. The result indicated that pyriproxyfen exposure can lead to silk gland injury, reduction of silk yield and cocooning rate. Furthermore, the expression levels of silk protein synthesis related genes were down regulated significantly. The same change trends were shown between PI3K/Akt and CncC/Keap1 pathway, which is the expressions of key genes can be elevated by pyriproxyfen exposure. In addition, the activity of detoxification enzymes (P450, GST and CarE) and the expression levels of detoxification genes were elevated after pyriproxyfen exposure, suggesting that detoxification enzymes may play an important role in detoxification of pyriproxyfen in silk gland. These results provided possible clues to the silk gland injury and gene transcriptional level changes in silkworm after pyriproxyfen exposure.

Cocoonase is indispensable for Lepidoptera insects breaking the sealed cocoon

Many insects spin cocoons to protect the pupae from unfavorable environments and predators. After emerging from the pupa, the moths must escape from the sealed cocoons. Previous works identified cocoonase as the active enzyme loosening the cocoon to form an escape-hatch. Here, using bioinformatics tools, we show that cocoonase is specific to Lepidoptera and that it probably existed before the occurrence of lepidopteran insects spinning cocoons. Despite differences in cocooning behavior, we further show that cocoonase evolved by purification selection in Lepidoptera and that the selection is more intense in lepidopteran insects spinning sealed cocoons. Experimentally, we applied gene editing techniques to the silkworm Bombyx mori, which spins a dense and sealed cocoon, as a model of lepidopteran insects spinning sealed cocoons. We knocked out cocoonase using the CRISPR/Cas9 system. The adults of homozygous knock-out mutants were completely formed and viable but stayed trapped and died naturally in the cocoon. This is the first experimental and phenotypic evidence that cocoonase is the determining factor for breaking the cocoon. This work led to a novel silkworm strain yielding permanently intact cocoons and provides a new strategy for controlling the pests that form cocoons.

Spinning and cocooning are the instincts of many insects, providing a shelter to the residing pupae to resist adverse factors. After the metamorphosis of pupa into adult, the adult must break open the cocoon to emerge, which is called decocooning. We have bioinformatically identified that cocoonase is specific to Lepidoptera, and demonstrated that it is the determining factor for breaking the sealed cocoon experimentally for the first time. This work led to a novel silkworm strain yielding permanently intact cocoons and provides a new strategy for controlling the pests that form cocoons and for breeding.