Developmental and transcriptomic features characterize defects of silk gland growth and silk production in silkworm naked pupa mutant

Effects of poly (ADP-ribose) polymerase 1 (PARP1) on silk proteins in the silkworm, Bombyx mori

Animal silk is economically important, while silk secretion is a complex and subtle mechanism regulated by many genes. We identified the poly (ADP-ribose) polymerase (PARP1) gene of the silkworm and successfully cloned its coding sequence (CDS) sequence. Using clustered regularly interspaced short palindromic repeat (CRISPR/Cas9) technology, we screened single guide RNA (sgRNA) with high knockout efficiency by cellular experiments and obtained PARP1 mutants by knocking out the PARP1 gene of the silkworm at the individual level. We found that the mutants mainly exhibited phenotypes such as smaller cocoon size and reduced cocoon shell rate than the wild type. We also detected the expression of silk protein genes in the mutant by quantitative real-time PCR (qPCR) and found that the expression of some silk protein genes was slightly down-regulated. Meanwhile, together with the results of transcriptomic analysis, we hypothesized that PARP1 may affect the synthesis of silk proteins, resulting in their failure to function properly. Our study may provide an important reference for future in-depth refinement of the molecular mechanism of silk protein expression in silk-producing animals, as well as a potential idea for future development of molecular breeding lines of silkworms to improve silk production.

Developmental and nuclear proteomic signatures characterize the temporal regulation of fibroin synthesis during the last molting-feeding transition of silkworm

Silkworm fibroins are natural proteinaceous macromolecules and provide core mechanical properties to silk fibers. The synthesis process of fibroins is posterior silk gland (PSG)-exclusive and appears active at the feeding stage and inactive at the molting stage. However, the molecular mechanisms controlling it remain elusive. Here, the silk gland's physiological and nuclear proteomic features were used to characterize changes in its structure and development from molting to feeding stages. The temporal expression profile and immunofluorescence analyses revealed a synchronous transcriptional on-off mode of fibroin genes. Next, the comparative nuclear proteome of the PSG during the last molting-feeding transition identified 798 differentially abundant proteins (DAPs), including 42 transcription factors and 15 epigenetic factors. Protein-protein interaction network analysis showed a "CTCF-FOX-HOX-SOX" association with activated expressions at the molting stage, suggesting a relatively complex and multifactorial regulation of the PSG at the molting stage. In addition, FAIRE-seq verification indicated "closed" and "open" conformations of fibroin gene promoters at the molting and feeding stages, respectively. Such proteome combined with chromatin accessibility analysis revealed the detailed signature of protein factors involved in the temporal regulation of fibroin synthesis and provided insights into silk gland development as well as silk production in silkworms.

Mutation in the Bombyx mori BmGMC2 gene impacts silk production and silk protein synthesis

Silk is a natural protein fiber that is predominantly comprised of fibroin and sericin. In addition, it contains seroins, protease inhibitors, enzymes, and other proteins. We found an ecdysone oxidase BmGMC2, notably, which is specifically and highly expressed only in the silk glands of silkworms (Bombyx mori L.). It is also one of the main components of non-cocoon silk, however, its precise function remains unclear. In this study, we examined the spatiotemporal expression pattern of this protein and obtained a homozygous mutant strain (K-GMC2) using the CRISPR-Cas9 system. Compared to the wild-type strain (WT), the silk production and main silk proteins significantly decreased in the larval stage, and the adhesive strength of native silk proteins decreased in the final instar. Proteomic data indicated the abundance of ribosomal proteins decreased significantly in K-GMC2, differentially expressed proteins (DEPs) were enriched in pathways related to neurodegenerative diseases and genetic information processing, indicating that knockout may lead to a certain degree of cell stress, affecting the synthesis of silk proteins. This study investigated the expression pattern and gene function of ecdysone oxidase BmGMC2 in silk and silk glands, laying the groundwork for understanding the role of enzymes in the production of silk fibers.

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

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

Transcriptomic analysis of the fat body of resistant and susceptible silkworm strains, Bombyx mori (Lepidoptera), after oral treatment with fenpropathrin

The widespread use of pyrethroid pesticides has brought serious economic losses in sericulture, but there is still no viable solution. The key to solving the problem is to improve silkworm resistance to pesticides, which depends on understanding the resistance mechanism of silkworms to pesticides. This study aimed to use transcriptomes to understand the underlying mechanism of silkworm resistance to fenpropathrin, which will provide a theoretical molecular reference for breeding pesticide-resistant silkworm varieties. In this study, the fat bodies of two strains with differential resistance after 12 h of fenpropathrin feeding were analyzed using RNA-Seq. After feeding fenpropathrin, 760 differentially expressed genes (DEGs) were obtained in the p50(r) strain and 671 DEGs in the 8y strain. The DEGs involved in resistance to fenpropathrin were further identified by comparing the two strains, including 207 upregulated DEGs in p50(r) and 175 downregulated DEGs in 8y. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that these fenpropathrin-related DEGs are mainly enriched in the metabolism and transporter pathways. Moreover, 28 DEGs involved in the metabolic pathway and 18 in the transporter pathway were identified. Furthermore, organic cation transporter protein 6 (BmOCT6), a transporter pathway member, was crucial in enhancing the tolerance of BmN cells to fenpropathrin. Finally, the knockdown of the expression of the homologs of BmOCT6 in Glyphodes pyloalis (G. pyloalis) significantly decreased the resistant level of larvae to fenpropathrin. The findings showed that the metabolism and transporter pathways are associated with resistance to fenpropathrin in silkworm, and OCT6 is an effective and potential target not only for silkworm breeding but also for pest biocontrol.

The silkworm (Bombyx mori) gut microbiota is involved in metabolic detoxification by glucosylation of plant toxins

Herbivores have evolved the ability to detoxify feed components through different mechanisms. The oligophagous silkworm feeds on Cudrania tricuspidata leaves (CTLs) instead of mulberry leaves for the purpose of producing special, high-quality silk. However, CTL-fed silkworms are found to have smaller bodies, slower growth and lower silk production than those fed mulberry leaves. Here, we show that the high content of prenylated isoflavones (PIFs) that occurred in CTLs is converted into glycosylated derivatives (GPIFs) in silkworm faeces through the silkworm gut microbiota, and this biotransformation is the key process in PIFs detoxification because GPIFs are found to be much less toxic, as revealed both in vitro and in vivo. Additionally, adding Bacillus subtilis as a probiotic to remodel the gut microbiota could beneficially promote silkworm growth and development. Consequently, this study provides meaningful guidance for silk production by improving the adaptability of CTL-fed silkworms.

Dynamic chromatin conformation and accessibility changes mediate the spatial-specific gene regulatory network in Bombyx mori

Silk gland genes of Bombyx mori can have strict spatial expression patterns, which impact their functions and silk quality; however, our understanding of their regulation mechanisms is currently insufficient. To address this, the middle silk gland (MSG) and posterior silk gland (PSG) of the silkworm were investigated. Gene ontology annotation showed that spatially specific expressed genes were involved in the formation of H3k9me and chromatin topology. Chromatin conformation data generated by Hi-C showed that the topologically associated domain boundaries around FibL and Sericin1 genes were significantly different between MSG and PSG. Changes in chromatin conformation led to changes in chromatin activity, which significantly affected the expression of nearby genes in silkworm. Chromatin accessibility regions of MSG and PSG were analyzed using FAIRE-seq, and 1006 transcription factor motifs were identified in open chromatin regions. Furthermore, the spatial-specific expression patterns of silk gland genes were mainly associated with homeobox-contained transcription factors, such as POU-M2, which was specifically bound and relatively highly expressed in the MSG. The regulatory network mediated by POU-M2 regulated most of the spatial-specific expressed genes in MSG, such as ADH1. These results can aid in improving silk performance, optimizing silkworm breeding, and improving the gene spatial regulatory model research for insects.

A Mechanical Performance Study of Dual Cured Thermoset Resin Systems 3D-Printed with Continuous Carbon Fiber Reinforcement

Additive manufacturing (AM) is one of the fastest-growing manufacturing technologies in modern times. One of the major challenges in the application of 3D-printed polymeric objects is expanding the applications to structural components, as they are often limited by their mechanical and thermal properties. To enhance the mechanical properties of 3D-printed thermoset polymer objects, reinforcing the polymer with continuous carbon fiber (CF) tow is an expanding direction of research and development. A 3D printer was constructed capable of printing with a continuous CF-reinforced dual curable thermoset resin system. Mechanical performance of the 3D-printed composites varied with the utilization of different resin chemistries. Three different commercially available violet light curable resins were mixed with a thermal initiator to improve curing by overcoming the shadowing effect of violet light by the CF. The resulting specimens’ compositions were analyzed, and then the specimens were mechanically characterized for comparison in tensile and flexural performance. The 3D-printed composites’ compositions were correlated to the printing parameters and resin characteristics. Slight enhancements in tensile and flexural properties from some commercially available resins over others appeared to be the result of better wet-out and adhesion.

Analysis of histomorphometric and proteome dynamics inside the silk gland lumen of Bombyx mori revealed the dynamic change of silk protein during the molt stage

Silkworms spin different silks at different growth stages for specific purposes. The silk spun before the end of each instar is stronger than that at the beginning of each instar and cocoon silk. However, the compositional changes in silk proteins during this process are unknown. Consequently, we performed histomorphological and proteomic analyses of the silk gland to characterize changes from the instar end to the next instar beginning. The silk glands were collected on day 3 of third- and fourth-instar larvae (III-3 and IV-3) and the beginning of fourth-instar larvae (IV-0). Proteomic analysis identified 2961 proteins from all silk glands. Silk proteins P25 and Ser5 were significantly more abundant in III-3 and IV-3 than in IV-0, and many cuticular proteins and protease inhibitors increased significantly in IV-0 compared with III-3 and IV-3. This shift may cause mechanical property differences between the instar end and beginning silk. Using section staining, qPCR, and western blotting, we found for the first time that silk proteins were degraded first and then resynthesized during the molting stage. Furthermore, we revealed that fibroinase mediated the changes of silk proteins during molting. Our results provide insights into the molecular mechanisms of silk proteins dynamic regulation during molting.

A molecular atlas reveals the tri-sectional spinning mechanism of spider dragline silk

The process of natural silk production in the spider major ampullate (Ma) gland endows dragline silk with extraordinary mechanical properties and the potential for biomimetic applications. However, the precise genetic roles of the Ma gland during this process remain unknown. Here, we performed a systematic molecular atlas of dragline silk production through a high-quality genome assembly for the golden orb-weaving spider Trichonephila clavata and a multiomics approach to defining the Ma gland tri-sectional architecture: Tail, Sac, and Duct. We uncovered a hierarchical biosynthesis of spidroins, organic acids, lipids, and chitin in the sectionalized Ma gland dedicated to fine silk constitution. The ordered secretion of spidroins was achieved by the synergetic regulation of epigenetic and ceRNA signatures for genomic group-distributed spidroin genes. Single-cellular and spatial RNA profiling identified ten cell types with partitioned functional division determining the tri-sectional organization of the Ma gland. Convergence analysis and genetic manipulation further validated that this tri-sectional architecture of the silk gland was analogous across Arthropoda and inextricably linked with silk formation. Collectively, our study provides multidimensional data that significantly expand the knowledge of spider dragline silk generation and ultimately benefit innovation in spider-inspired fibers.

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.

Comparative Silk Transcriptomics Illuminates Distinctive Impact of Artificial Selection in Silkworm Modern Breeding

Early domestication and the following improvement are two important processes in the cocoon silk evolution of silkworms. In contrast to early domestication, understanding of the improvement process is still fuzzy. By systematically comparing the larval silk gland transcriptomes of the wild, early domestic, and improved silkworms, we highlighted a novel landscape of transcriptome in the silk glands of improved ones. We first clarified that silk cocoon protein genes were up-regulated in modern breeding but not in early domestication. Furthermore, we found that differentially expressed genes (DEGs) between improved and early domestic silkworms (2711), as well as between improved and wild silkworms (2264), were obviously more than those between the early domestic and wild silkworms (158), with 1671 DEGs specific in the improved silkworm (IS-DEGs). Hierarchical clustering of all the DEGs consistently indicated that improved silkworms were significantly diverged from the early domestic and wild silkworms, suggesting that modern breeding might cause prompt and drastic dynamic changes of gene expression in the silk gland. We further paid attention to these 1671 IS-DEGs and were surprised to find that down-regulated genes were enriched in basic organonitrogen compound biosynthesis, RNA biosynthesis, and ribosome biogenesis processes, which are generally universally expressed, whereas those up-regulated genes were enriched in organonitrogen compound catabolic processes and functions involving in the dynamic regulation of protein post-translation of modification. We finally highlighted one candidate improvement gene among these up-regulated IS-DEGs, i.e., GDAP2, which may play roles in silk behavior and the overall robustness of the improved silkworm. The findings strongly suggest that modern breeding may facilitate effective control of the basic consumption of nitrogen and a stronger switch of nitrogen resources from other tissues to the silk glands, for an efficient supply for silk production, and implies the importance of brain behavior and robustness in silk yield improvement of modern breeding.

Comparative transcriptome and proteome reveal synergistic functions of differentially expressed genes and proteins implicated in an over-dominant silkworm heterosis of increased silk yield

We previously observed an over-dominant silkworm heterosis of increased yield in a cross of Bombyx mori nuclear polyhydrosis virus-resistant strain NB with a susceptible strain 306. In the present study, we found that heterosis also exists in crosses of NB with other susceptible strains, indicating it is a more general phenomenon. We performed comparative transcriptome and proteome and identified 1624 differentially expressed genes (DEGs) and 298 differentially expressed proteins (DEPs) in silk glands between parents and F1 hybrids, of which 24 DEGs/DEPs showed consistent expression at mRNA and protein levels revealed by Venn joint analysis. Their expressions are completely non-additive, mainly transgressive and under low-parent, suggesting recombination of parental genomes may be the major genetic mechanism for the heterosis. GO and KEGG analyses revealed that they may function in generally similar but distinctive aspects of metabolisms and processes with signal transduction and translation being most affected. Notably, they may not only up-regulate biosynthesis and transport of silk proteins but also down-regulate other unrelated processes, synergistically and globally remodelling the silk gland to increase yield and cause the heterosis. Our findings contribute insights into the understanding of silkworm heterosis and silk gland development and provide targets for transgenic manipulation to further increase the silk yield.

Combined analysis of silk synthesis and hemolymph amino acid metabolism reveal key roles for glycine in increasing silkworm silk yields

Rearing silkworms (Bombyx mori) using formula feed has revolutionized traditional mulberry feed strategies. However, low silk production efficiencies persist and have caused bottlenecks, hindering the industrial application of formula feed sericulture. Here, we investigated the effects of formula feed amino acid composition on silk yields. We showed that imbalanced amino acids reduced DNA proliferation, decreased Fib-H, Fib-L, and P25 gene expression, and caused mild autophagy in the posterior silk gland, reducing cocoon shell weight and ratio. When compared with mulberry leaves, Gly, Ala, Ser, and Tyr percentages of total amino acids in formula feed were decreased by 5.26%, while Glu and Arg percentages increased by 9.56%. These changes increased uric acid and several amino acids levels in the hemolymph of silkworms on formula feed. Further analyses showed that Gly and Thr (important synthetic Gly sources) increased silk yields, with Gly increasing amino acid conversion efficiencies to silk protein, and reducing urea levels in hemolymph. Also, Gly promoted endomitotic DNA synthesis in silk gland cells via phosphoinositide 3-kinase (PI3K)/Akt/target of rapamycin (TOR) signaling. In this study, we highlighted the important role of Gly in regulating silk yields in silkworms.

Mechanism of the growth and development of the posterior silk gland and silk secretion revealed by mutation of the fibroin light chain in silkworm

Silkworm, as a model organism, has very high economic value due to its silk secretion ability. Although a large number of studies have attempted to elucidate the mechanism of silk secretion, it remains unclear. In this study, the fibroin light chain (Fib-L) gene of silkworm was subjected to CRISPR/Cas9 editing, which yielded premature termination of translation at 135 aa. Compared with those of the wild type, the posterior silk glands (PSGs) of the homozygous mutants on the third day of the fifth instar showed obvious premature degeneration. Comparative transcriptome and proteomic analyses of the PSGs of wild-type individuals, heterozygous mutants and homozygous mutants were performed on the fourth day of the fifth instar. A GO enrichment analysis showed that the differentially expressed genes (DEGs) between homozygous mutants and wild-type individuals were enriched in cytoskeleton-related terms, and a KEGG enrichment analysis showed that the upregulated DEGs between homozygous mutants and wild-type individuals were enriched in the phagosome and apoptosis pathways. These results indicated that apoptosis was activated prematurely in the PSGs of homozygous mutants. Furthermore, autophagy and heat shock response were activated in the PSGs of homozygous mutants, as demonstrated by an analysis of the DEGs related to autophagy and heat shock. A comparative proteomic analysis further confirmed that autophagy, apoptosis and the heat shock response were activated in the PSGs of homozygous mutants, which led to premature degradation of the PSGs. These results provide insights for obtaining a more in-depth understanding of the mechanism of silk secretion in silkworms.

CRISPR-Mediated Endogenous Activation of Fibroin Heavy Chain Gene Triggers Cellular Stress Responses in Bombyx mori Embryonic Cells

Simple Summary

Based on a CRISPRa approach, activating endogenous fibroin heavy chain (FibH) gene expression in Bombyx mori embryonic (BmE) cells, which was driven by a combination of the dCas9-VPR (a tripartite activator, composed of VP64, p65, and Rta) and the sgRNA targeting to the promoter of FibH gene, was performed for investigating the biological roles of FibH in the development of silk gland cells. The activation of the endogenous FibH gene lead to up-regulation of cellular stress responses-related genes, which suggested a significant positive correlation between activated FibH gene expression and cellular stress responses. Moreover, the present findings might provide a potential model for studying the cellular stress responses caused by silk secretion disorder and lay a foundation for the understanding of silk gland development in silk-spinning insects.

Abstract

The silkworm Bombyx mori is an economically important insect, as it is the main producer of silk. Fibroin heavy chain (FibH) gene, encoding the core component of silk protein, is specifically and highly expressed in silk gland cells but not in the other cells. Although the silkworm FibH gene has been well studied in transcriptional regulation, its biological functions in the development of silk gland cells remain elusive. In this study, we constructed a CRISPRa system to activate the endogenous transcription of FibH in Bombyx mori embryonic (BmE) cells, and the mRNA expression of FibH was successfully activated. In addition, we found that FibH expression was increased to a maximum at 60 h after transient transfection of sgRNA/dCas9-VPR at a molar ratio of 9:1. The qRT-PCR analysis showed that the expression levels of cellular stress response-related genes were significantly up-regulated along with activated FibH gene. Moreover, the lyso-tracker red and monodansylcadaverine (MDC) staining assays revealed an apparent appearance of autophagy in FibH-activated BmE cells. Therefore, we conclude that the activation of FibH gene leads to up-regulation of cellular stress responses-related genes in BmE cells, which is essential for understanding silk gland development and the fibroin secretion process in B. mori.

Comparative mRNA and LncRNA Analysis of the Molecular Mechanisms Associated With Low Silk Production in Bombyx mori

Naked pupa sericin and Naked pupa are two mutant strains of Bombyx mori with extremely low or no fibroin production compared to the Qiufeng and Baiyu strains, both of which exhibit very high silk fibroin production. However, the molecular mechanisms by which long non-coding RNAs regulate fibroin synthesis need further study. In this study, we performed high-throughput RNA-seq to investigate lncRNA and mRNA expression profiles in the posterior silk gland of Qiufeng, Baiyu, Nd-s^D, and Nd silkworms at the third day of the 5th instar. Our efforts yielded 26,767 novel lncRNAs and 6,009 novel mRNAs, the expression levels of silk protein genes and silk gland transcription factors were decreased in Qiufeng vs. Nd-s^D and Qiufeng vs. Nd, while those of many genes related to autophagy, apoptosis, RNA degradation, ubiquitin-mediated proteolysis and heat shock proteins were increased. Moreover, the expression of a large number of genes responsible for protein synthesis and secretion was significantly decreased in Nd. GO and KEGG analysis results showed that nucleotide excision repair, mRNA surveillance pathways, amino acid degradation, protein digestion and absorption, ER-associated degradation and proteasome pathways were significantly enriched for the Qiufeng vs. Nd-s^D and Qiufeng vs. Nd comparisons. In conclusion, our findings contribute to the lncRNA and mRNA database of Bombyx mori, and the identified differentially expressed mRNAs and lncRNAs help to reveal the molecular mechanisms of low silk production in Nd-s^D and Nd, providing new insights for improvement of silk yield and elucidation of silk mechanical properties.

A transcriptome atlas of silkworm silk glands revealed by PacBio single-molecule long-read sequencing

The silk gland of the silkworm Bombyx mori is a specialized organ where silk proteins are efficiently synthesized under precise regulation that largely determines the properties of silk fibers. To understand the genes involved in the regulation of silk protein synthesis, considerable research has focused on the transcripts expressed in silk glands; however, the complete transcriptome profile of this organ has yet to be elucidated. Here, we report a full-length silk gland transcriptome obtained by PacBio single-molecule long-read sequencing technology. In total, 11,697 non-redundant transcripts were identified in mixed samples of silk glands dissected from larvae at five developmental stages. When compared with the published reference, the full-length transcripts optimized the structures of 3002 known genes, and a total of 9061 novel transcripts with an average length of 2171 bp were detected. Among these, 1403 (15.5%) novel transcripts were computationally revealed to be lncRNAs, 8135 (89.8%) novel transcripts were annotated to different protein and nucleotide databases, and 5655 (62.4%) novel transcripts were predicted to have complete ORFs. Furthermore, we found 1867 alternative splicing events, 2529 alternative polyadenylation events, 784 fusion events and 6596 SSRs. This study provides a comprehensive set of reference transcripts and greatly revises and expands the available silkworm transcript data. In addition, these data will be very useful for studying the regulatory mechanisms of silk protein synthesis.

SilkDB 3.0: visualizing and exploring multiple levels of data for silkworm

SilkDB is an open-accessibility database and powerful platform that provides comprehensive information on the silkworm (Bombyx mori) genome. Since SilkDB 2.0 was released 10 years ago, vast quantities of data about multiple aspects of the silkworm have been generated, including genome, transcriptome, Hi-C and pangenome. To visualize data at these different biological levels, we present SilkDB 3.0 (https://silkdb.bioinfotoolkits.net), a visual analytic tool for exploring silkworm data through an interactive user interface. The database contains a high-quality chromosome-level assembly of the silkworm genome, and its coding sequences and gene sets are more accurate than those in the previous version. SilkDB 3.0 provides a view of the information for each gene at the levels of sequence, protein structure, gene family, orthology, synteny, genome organization and gives access to gene expression information, genetic variation and genome interaction map. A set of visualization tools are available to display the abundant information in the above datasets. With an improved interactive user interface for the integration of large data sets, the updated SilkDB 3.0 database will be a valuable resource for the silkworm and insect research community.