Transcriptomic Analysis of the Anterior Silk Gland in the Domestic Silkworm (Bombyx mori) - Insight into the Mechanism of Silk Formation and Spinning

Ectopic expression of BmeryCA in Bombyx mori increases silk yield and mechanical properties by altering the pH of posterior silk gland

The silk glands are the specialized tissue where silk protein synthesis, secretion, and conformational transitions take place, with pH playing a critical role in both silk protein synthesis and fiber formation. In the present study, we have identified erythrocyte carbonic anhydrase (BmeryCA) belonging to the α-CA class in the silk gland, which is a Zn2+ dependent metalloenzyme capable of efficiently and reversibly catalyzing the hydrated reaction of CO2 to HCO3-, thus participating in the regulation of acid-base balance. Multiple sequence alignments revealed that the active site of BmeryCA was highly conserved. Tissue expression profiling showed that BmeryCA had relatively high expression levels in hemolymph and epidermis but is barely expressed in the posterior silk gland (PSG). By specifically overexpressing BmeryCA in the PSG, we generated transgenic silkworms. Ion-selective microelectrode (ISM) measurements demonstrated that specifically overexpression of BmeryCA in the PSG led to a shift in pH from weakly alkaline to slightly neutral conditions. Moreover, the resultant PSG-specific BmeryCA overexpression mutant strain displayed a significant increase in both silk yield and silk fiber mechanical properties. Our research provided new insights into enhancing silk yield and improving the mechanical properties of silk fibers.

Identification of genes associated with the silk gland size using multi-omics in silkworm (Bombyx mori)

Silk gland size in silkworms (Bombyx mori) affects silk output. However, the molecular mechanisms by which genes regulate silk gland size remain unclear. In this study, silk glands from three pure silkworm strains (A798, A306 and XH) with different silk gland weight phenotypes were compared using transcriptomics and proteomics to identify differentially expressed genes (DEGs) and proteins (DEPs). When comparing A798 to A306 and A798 to XH, 830 and 469 DEGs were up-regulated, respectively. These genes were related to the gene ontology terms, metabolic process, transport activity and biosynthesis process. In addition, 372 and 302 up-regulated differentially expressed proteins were detected in A798 to A306 and A798 to XH, respectively, related to the gene ontology terms, ribosome and protein export, ribosome and polypeptide biosynthesis processes. Moreover, combined transcriptomics, proteomics and weighted correlation network analyses showed that five genes (BGIBMGA002524, BGIBMGA002629, BGIBMGA005659, BGIBMGA005711 and BGIBMGA010889) were significantly associated with the silk gland weight. Reverse Transcription-quantitative real-time Polymerase Chain Reaction (RT-qPCR) and Enzyme linked immunosorbent assay (ELISA) were used to verify the mRNA and protein expression of five genes in the silk glands and tissues of 18 silkworm strains. The results showed that four genes have higher expression levels in heavier silk glands. These genes are associated with glycogen metabolism, fatty acid synthesis and branched chain amino acid metabolism, thus potentially promoting growth and silk protein synthesis. These findings provide valuable insights into the molecular mechanisms underlying the relationship between silk gland weight and silk yield in silkworms.

Silk Fibroin Bioink for 3D Printing in Tissue Regeneration: Controlled Release of MSC extracellular Vesicles

Sodium alginate (SA)-based hydrogels are often employed as bioink for three-dimensional (3D) scaffold bioprinting. They offer a suitable environment for cell proliferation and differentiation during tissue regeneration and also control the release of growth factors and mesenchymal stem cell secretome, which is useful for scaffold biointegration. However, such hydrogels show poor mechanical properties, fast-release kinetics, and low biological performance, hampering their successful clinical application. In this work, silk fibroin (SF), a protein with excellent biomechanical properties frequently used for controlled drug release, was blended with SA to obtain improved bioink and scaffold properties. Firstly, we produced a printable SA solution containing SF capable of the conformational change from Silk I (random coil) to Silk II (β-sheet): this transition is a fundamental condition to improve the scaffold's mechanical properties. Then, the SA-SF blends' printability and shape fidelity were demonstrated, and mechanical characterization of the printed hydrogels was performed: SF significantly increased compressive elastic modulus, while no influence on tensile response was detected. Finally, the release profile of Lyosecretome-a freeze-dried formulation of MSC-secretome containing extracellular vesicles (EV)-from scaffolds was determined: SF not only dramatically slowed the EV release rate, but also modified the kinetics and mechanism release with respect to the baseline of SA hydrogel. Overall, these results lay the foundation for the development of SA-SF bioinks with modulable mechanical and EV-release properties, and their application in 3D scaffold printing.

Molecular Characterization of the Functional Genes Associated with Silk Assembly, Transport, and Protection in the Silk Glands of Popular Multivoltine Breeds of Silkworm Bombyx mori. L

Bombyx mori is an agriculturally important insect used extensively for silk production. India, especially the eastern regions, is mostly dependent on the multivoltine breeds of silkworm Bombyx mori and their hybrids/crossbreeds. The multivoltine breeds are indigenous and superior in survival and hardiness but are relatively inferior in terms of qualitative traits, typically the silk quality. Therefore, it is highly relevant to understand the mechanism of silk production in the multivoltine breeds to decipher the reasons for the inferior quality of silk produced by the multivoltine breeds and thus gain leads to improve the quality of silk production in multivoltine breeds. With this background, study was carried to identify differential expression of the major genes associated with silk proteins in the silk gland region of the popular multivoltine breeds. Our results indicated that although fib-L, fib-H, Sericins, and P25 are the major genes associated with silk filament, a few other genes associated with silk assembly, transport, and protection in the silk glands are the ones that largely contribute towards efficient silk production. The differential expression of these genes had a major effect on the movement of silk proteins within the silk gland and the efficiency of silk production as well. The Pearson correlation revealed a positive correlation amongst the genes dealt with in this study, indicating that the concurrent increase in expression of both the types of genes in the silk glands, significantly improves the silk production.

The differences in cocoon and silk qualities among sex-related mulberry and silkworm feeding groups

Sexual dimorphism is seen in many dioecious plant and animal species, which may influence their trophic interactions. The differences in trophic interactions derived from sexual dimorphism in plants may influence herbivorous performance and population dynamics. Both silkworm (Bombyx mori L.) and mulberry (Morus alba L.) usually exhibit sexual dimorphism. However, few studies have been conducted on the effect of sex-related silkworm and mulberry pairings on cocoon and silk qualities, which are crucial in sericulture. Here, we compared the differences in cocoon and silk qualities under four feeding combinations (FS-FL: female silkworm fed with leaves from female mulberry trees; MS-FL: male silkworm fed with leaves from female mulberry trees; FS-ML: female silkworm fed with leaves from male mulberry trees; MS-ML: male silkworm fed with leaves from male mulberry trees). The results showed that silkworms exhibited male-biased herbivory with more male mulberry leaves digested. The FS-ML group had higher silk weight and silk ratio of fresh cocoons than the FS-FL group, and the MS-ML group had lower coarse points than the MS-FL group. Compared with groups FS-ML and FS-FL, both MS-FL and MS-ML had smaller cocoons with longer silk lengths and a higher silk ratio of the fresh cocoons. In addition, the Entropy Weight-TOPSIS method showed the cocoon quality rank as FS-ML > FS-FL > MS-FL > MS-ML, whereas silk quality rank was MS-ML > FS-FL > FS-ML > MS-FL. These results indicate that the quality of cocoon and silk is related to the interaction of silkworm and mulberry at the sex level. Furthermore, female silkworms fed with female and male tree leaves have a higher total yield in cocoon production, while male silkworms fed with male tree leaves produced higher silk quality.

A single-cell transcriptomic atlas characterizes the silk-producing organ in the silkworm

The silk gland of the domesticated silkworm Bombyx mori, is a remarkable organ that produces vast amounts of silk with exceptional properties. Little is known about which silk gland cells execute silk protein synthesis and its precise spatiotemporal control. Here, we use single-cell RNA sequencing to build a comprehensive cell atlas of the silkworm silk gland, consisting of 14,972 high-quality cells representing 10 distinct cell types, in three early developmental stages. We annotate all 10 cell types and determine their distributions in each region of the silk gland. Additionally, we decode the developmental trajectory and gene expression status of silk gland cells. Finally, we discover marker genes involved in the regulation of silk gland development and silk protein synthesis. Altogether, this work reveals the heterogeneity of silkworm silk gland cells and their gene expression dynamics, affording a deeper understanding of silk-producing organs at the single-cell level.

The molecular underpinning of silk-producing organs is not well characterized. Here the authors use single-cell RNA sequencing to build an atlas of the silkworm silk gland and reveal the heterogeneity of silk gland cells.

Chitin and cuticle proteins form the cuticular layer in the spinning duct of silkworm

Chitin is found in the exoskeleton and peritrophic matrix of arthropods, but recent studies have also identified chitin in the spinning duct of silk-spinning arthropods. Here, we report the presence and function of chitin and cuticle proteins ASSCP1 and ASSCP2 in the spinning duct of silkworm. We show that chitin and these proteins are co-located in the cuticular layer of the spinning duct. Ultrastructural analysis indicates that the cuticular layer has a multilayer structure by layered stacking of the chitin laminae. After knocking down ASSCP1 and ASSCP2, the fine structure of this layer was disrupted, which had negative impacts on the mechanical properties of silk. This work clarifies the function of chitin in the spinning duct of silkworm. Chitin and cuticle proteins are the main components of the cuticular layer, providing the shearing stress during silk fibrillogenesis and regulating the final mechanical properties of silk. STATEMENT OF SIGNIFICANCE: Recent studies have identified chitin in the spinning duct of silk-spinning arthropods. However, the role of chitin in this specific organ remains unclear. This study reports that chitin and cuticle proteins form the cuticular layer, a unique structure of the spinning duct of silkworm. This layer with a precise laminate structure gives the spinning duct flexible properties, provides shearing forces for silk fibrillogenesis, and contributes to silk final mechanical properties. Our work clarifies the component, ultrastructure, and biological significance of the silkworm cuticular layer, describes the specific process of silk fiber formation, and proposes new molecular targets (chitin and cuticle proteins) for the improvement of animal silks.

Co-Expression Network and Time-Course Expression Analyses to Identify Silk Protein Regulatory Factors in Bombyx mori

Simple Summary

Previous studies have reported how the silk production ability of Bombyx mori can be enhanced, but the mechanism that regulates silk protein genes remains unclear. We performed co-expression network analysis using networkz, an in-house program, which led to the identification of 91 transcription factors were co-expressed with silk protein genes. Of them, 13 transcripts were identified to be novel regulatory factors by time-course expression analysis during the fifth instar larvae stage. Their expression patterns were highly relevant to those of silk protein genes. Our results suggest that the two-step expression screening was robust and highly sensitive to screen relative genes, and a complex mechanism regulates silk protein production in B. mori. The novel candidates that were identified herein can serve as key genes to develop methods to enhance the silk protein production ability of B. mori.

Abstract

Bombyx mori is an important economic insect and an animal model in pharmacomedical research. Although its physiology has been studied for many years, the mechanism via which silk protein genes are regulated remains unclear. In this study, we performed two-step expression screening, namely co-expression network and time-course expression analyses to screen silk protein regulation factors. A co-expression network analysis using RNA-seq data that were obtained from various tissues, including the silk glands of B. mori, was performed to identify novel silk protein regulatory factors. Overall, 91 transcription factors, including some known ones, were found to be co-expressed with silk protein genes. Furthermore, time-course expression analysis during the fifth instar larvae stage revealed that the expression pattern of 13 novel transcription factors was highly relevant to that of silk protein genes and their known regulatory factor genes. In particular, the expression peak of several transcription factors (TFs) was detected before the expression of silk protein genes peak. These results indicated that a larger number of genes than expected may be involved in silk protein regulation in B. mori. Functional analyses of function-unknown transcription factors should enhance our understanding of this system.

Genome-wide survey and characterization of transcription factors in the silk gland of the silkworm, Bombyx mori

Transcription factors (TFs) are key proteins that modulate gene transcription and thereby lead to changes in the gene expression profile and the subsequent alteration of cellular functions. In the silk gland (SG) of silkworm Bombyx mori, an important silk-producing insect, TFs are of vital importance in the regulation of silk protein synthesis in this organ. However, which TFs exist and express in the SG remains largely unknown. Here, we report the large-scale identification of TFs in the SG based on available full-length transcript sequences and the most recent version of silkworm genome data. In total, 348 candidate TFs were identified by strict filtration and were classified into 56 TF families. Chromosomal distribution, motif composition, and phylogenetic relationship analyses revealed the typical characteristics of these TFs. In addition, the expression patterns of 348 TFs in various tissues of B. mori, especially the SG of fourth-molt (4LM) and day-3 and day-4 fifth-instar (5L3D and 5L4D) larvae, were investigated based on public RNA-seq data and gene microarray data, followed by spatiotemporal verification of TF expression levels by quantitative real-time PCR (qRT-PCR). This report describes the first comprehensive analysis of TFs in the B. mori SG. The results can serve as a baseline for further studies of the roles of TFs in the B. mori SG.

Transcriptomic analysis of the bagworm moth silk gland reveals a number of silk genes conserved within Lepidoptera

Lepidopteran insects produce cocoons with unique properties. The cocoons are made of silk produced in the larval tissue silk gland and our understanding of the silk genes is still very limited. Here, we investigated silk genes in the bagworm moth Eumeta variegata, a species that has recently been found to produce extraordinarily strong and tough silk. Using short-read transcriptomic analysis, we identified a partial sequence of the fibroin heavy chain gene and its product was found to have a C-terminal structure that is conserved within nonsaturniid species. This is in accordance with the presence of fibroin light chain/fibrohexamerin genes and it is suggested that the bagworm moth is producing silk composed of fibroin ternary complex. This indicates that the fibroin structure has been evolutionarily conserved longer than previously thought. Other than fibroins we identified candidates for sericin genes, expressed strongly in the middle region of the silk gland and encoding serine-rich proteins, and other silk genes, that are structurally conserved with other lepidopteran homologues. The bagworm moth is thus considered to be producing conventional lepidopteran type of silk. We further found a number of genes expressed in a specific region of the silk gland and some genes showed conserved expression with Bombyx mori counterparts. This is the first study allowing comprehensive silk gene identification and expression analysis in the lepidopteran Psychidae family and should contribute to the understanding of silk gene evolution as well as to the development of novel types of silk.

Reference Transcriptome Data in Silkworm Bombyx mori

Herein, we performed RNA-seq analysis of ten major tissues/subparts of silkworm larvae. The sequences were mapped onto the reference genome assembly and the reference transcriptome data were successfully constructed. The reference data provided a nearly complete sequence for sericin-1, a major silk gene with a complex structure. We also markedly improved the gene model for other genes. The transcriptomic expression was investigated in each tissue and a number of transcripts were identified that were exclusively expressed in tissues such as the testis. Transcripts strongly expressed in the midgut formed tight genomic clusters, suggesting that they originated from tandem gene duplication. Transcriptional factor genes expressed in specific tissues or the silk gland subparts were also identified. We successfully constructed reference transcriptome data in the silkworm and found that a number of transcripts showed unique expression profiles. These results will facilitate basic studies on the silkworm and accelerate its applications, which will contribute to further advances in lepidopteran and entomological research as well as the practical use of these insects.

Cytochrome P450 monooxygenase genes in the wild silkworm, Bombyx mandarina

Wild (Bombyx mandarina) and domestic silkworms (B. mori) are good models for investigating insect domestication, as 5000 years of artificial breeding and selection have resulted in significant differences between B. mandarina and B. mori. In this study, we improved the genome assemblies to the chromosome level and updated the protein-coding gene annotations for B. mandarina. Based on this updated genome, we identified 68 cytochrome P450 genes in B. mandarina. The cytochrome P450 repository in B. mandarina is smaller than in B. mori. Certain currently unknown key genes, rather than gene number, are critical for insecticide resistance in B. mandarina, which shows greater resistance to insecticides than B. mori. Based on the physical maps of B. mandarina, we located 66 cytochrome P450s on 18 different chromosomes, and 27 of the cytochrome P450 genes were concentrated into seven clusters. KEGG enrichment analysis of the P450 genes revealed the involvement of cytochrome P450 genes in hormone biosynthesis. Analyses of the silk gland transcriptome identified candidate cytochrome P450 genes (CYP306A) involved in ecdysteroidogenesis and insecticide metabolism in B. mandarina.

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.

Disruption of the Metal Ion Environment by EDTA for Silk Formation Affects the Mechanical Properties of Silkworm Silk

Silk fiber has become a research focus because of its comprehensive mechanical properties. Metal ions can influence the conformational transition of silk fibroin. Current research is mainly focused on the role of a single ion, rather than the whole metal ion environment. Here, we report the effects of the overall metal ion environment on the secondary structure and mechanical properties of silk fibers after direct injection and feeding of silkworms with EDTA. The metal composition of the hemolymph, silk gland, and silk fiber changed significantly post EDTA treatment. Synchrotron FTIR analysis indicated that the secondary structure of silk fiber after EDTA treatment changed dramatically; particularly, the β-sheets decreased and the β-turns increased. Post EDTA treatment, the silk fiber had significantly decreased strength, Young’s modulus, and toughness as compared with the control groups, while the strain exhibited no obvious change. These changes can be attributed to the change in the metal ion environment in the silk fibroin and sericin in the silk gland. Our investigation provides a new theoretical basis for the natural silk spinning process, and our findings could help develop a method to modify the mechanical properties of silk fiber using metal ions.

Deep Insight into the Transcriptome of the Single Silk Gland of Bombyx mori

The silk gland synthesizes and secretes a large amount of protein and stores liquid silk protein at an extremely high concentration. Interestingly, silk proteins and serine protease inhibitors are orderly arranged in the silk gland lumen and cocoon shells. Silk fiber formation and the spinning mechanism have not been fully elucidated. Therefore, we conducted a comparative transcriptome analysis of seven segments of the single silk gland to characterize internal changes in the silk gland during the 5th instar of mature larvae. In total, 3121 differentially expressed genes were identified in the seven segments. Genes highly expressed in the middle silk gland (MSG) were mainly involved in unsaturated fatty acid biosynthesis, fatty acid metabolism, apoptosis—fly, and lysosome pathways, whereas genes highly expressed in the posterior silk gland (PSG) were mainly involved in ribosome, proteasome, citrate cycle, and glycolysis/gluconeogenesis pathways. Thus, the MSG and PSG differ greatly in energy source use and function. Further, 773 gradually upregulated genes (from PSG to MSG) were involved in energy metabolism, silk protein synthesis, and secretion, suggesting that these genes play an important role in silk fiber formation. Our findings provide insights into the mechanism of silk protein synthesis and transport and silk fiber formation.

Comparative Proteome Analysis Reveals that Cuticular Proteins Analogous to Peritrophin-Motif Proteins are Involved in the Regeneration of Chitin Layer in the Silk Gland of Bombyx mori at the Molting Stage

The silk gland of silkworm produces silk proteins during larval development. Many studies have long focused on the silk gland of the fifth instar larvae, but few have investigated this gland at other larval stages. In the present study, the silk gland proteomes of the fourth instar and fourth molt are analyzed using liquid chromatography-tandem mass spectrometry. In total, 2654 proteins are identified from the silk gland. A high abundance of ribosomal proteins and RR-motif chitin-binding proteins is identified during day 2 of the fourth instar (IV-2) larval developmental stage, and the expression of cuticular proteins analogous to peritrophin (CPAP)-motif chitin-binding proteins is higher during the fourth molt (IV-M). In all, nine enzymes are found to be involved in the chitin regeneration pathway in the silk gland. Among them, two chitinase and two chitin deacetylases are identified as CPAP-motif proteins. Furthermore, the expression of CPAP3-G, the most abundant CPAP-motif cuticular protein in the silk gland during the IV-M stage, is investigated using western blot and immunofluorescence analyses; CPAP3-G shows a reverse changing trend with chitin in the silk gland. The findings of this study suggest that CPAP-motif chitin-binding proteins are involved in the degradation of the chitin layer in the silk gland. The data have been deposited to the ProteomeXchange with identifier PXD008677.

Comparative transcriptome analysis of silkworm, Bombyx mori colleterial gland suggests their functional role in mucous secretion

Colleterial glands (CG) present in the body of adult female of Bombyx mori, which can help adhere eggs on the surface of the host plants. Although this organ has been known for centuries, only morphology and its secretions have been studied. Their gene expression profiles and physiological roles remain largely unknown. Aided by high-throughput next generation sequencing (NGS), we reported the comparative transcriptome analysis of CG isolated from the H9 and the P50 strains of Bombyx mori. A total of 19,896,957 and 20,446,366 clean reads were obtained from CG of H9 and the P50 strains, respectively; then differential expression analysis was performed, and 1,509 differentially expressed genes (DEGs) were identified. Among them, 1,001 genes are up-regulated and 508 genes are down-regulated in P50 individuals compared with H9 individuals. The enrichment of GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) of DEGs confirmed that many DEGs were associated with "Amino acid transport and metabolism", "Nucleotide transport and metabolism", and "Inorganic ion transport and metabolism", 25 of the DEGs related to the "ECM-receptor interaction passway", "sphingolipid metabolism passway", and "amino sugar and nucleotide sugar metabolism passway" were potentially involved in the process of CG development and mucus secretion. According to these data, we hypothesized that CG play an important role in providing favorable physiological environment for the glue secretion formation. In addition, GO enrichment and differential expression analysis of the DEGs in the CG indicate that this gland may be involved in the transporting of small solutes such as sugars, ions, amino acids and nucleotide sugar to the CG. Our findings lay the foundation for further research on CG function.

Proteome profile of spinneret from the silkworm, Bombyx mori

The silkworm spinneret is an important tissue for silk fibrillogenesis and spinning. All biochemical processes during silk fibrillogenesis are correlated with silk properties. Understanding the role of spinneret in silk fibrillogenesis may help to reveal the mechanism of silk fibrillogenesis as well as improve silk quality for commercial purposes. Thus, we profiled the proteome of silkworm spinneret. A total of 1572 proteins and 232 differential abundance proteins were identified. Silk fibrillogenesis-related proteins, such as cuticle proteins, ion-transporting proteins, muscular proteins, and energy metabolic proteins, were abundant in spinneret. Metabolic pathway and GO enrichment analyses revealed that the identified proteins were involved in energy metabolism, chitin binding, and cuticle construction. Active energy metabolism may provide abundant energy for the muscle contraction as well as ion and water exchange. The chitin binding and cuticle construction process may provide sufficient shear forces for silk formation. Our data suggest that silkworm spinneret provides a suitable physiological and biochemical environment for silk fibrillogenesis. These proteins are potential targets for improving silk quality in the silk industry. Data are available via ProteomeXchange with identifier PXD004455.

Osiris9a is a major component of silk fiber in lepidopteran insects

In a previous high-throughput proteomics study, it was found that the silkworm cocoon contains hundreds of complex proteins, many of which have unknown functions, in addition to fibroins, sericins, and some protease inhibitors. Osiris was one of the proteins with no known function. In this study, we identified the Osiris gene family members and constructed a phylogenetic tree based on the sequences from different species. Our results indicate that the Osiris9 gene subfamily contains six members; it is specifically expressed in lepidopteran insects and has evolved by gene duplication. An Osiris gene family member from Bombyx mori was designated as BmOsiris9a (BmOsi9a) on the basis of its homology to Drosophila melanogaster Osiris9. The expression pattern of BmOsi9a showed that it was highly expressed only in the middle silk gland of silkworm larvae, similar to Sericin1 (Ser1). BmOsi9a was visualized as two bands in western blot analysis, implying that it probably undergoes post-translational modifications. Immunohistochemistry analysis revealed that BmOsi9a was synthesized and secreted into the lumen of the middle silk gland, and was localized in the sericin layer in the silk fiber. BmOsi9a was found in the silk fibers of not only three Bombycidae species, viz. B. mori, B. mandarina, and B. huttoni, but also in the fibers collected from Saturniidae species, including Antheraea assama, Antheraea mylitta, and Samia cynthia. Although the exact biological function of Osi9a in the silk fibers is unknown, our results are important because they demonstrate that Osi9a is a common structural component of silk fiber and is expressed widely among the silk-producing Bombycidae and Saturniidae insects. Our results should help in understanding the role of Osi9a in silk fibers.

Histomorphometric and transcriptomic features characterize silk glands' development during the molt to intermolt transition process in silkworm

The molt-intermolt cycle is an essential feature in holometabolous and hemimetabolous insects' development. In the silkworm, silk glands are under dramatic morphological and functional changes with fibroin genes' transcription being repeatedly turned off and on during the molt-intermolt cycles. However, the molecular mechanisms controlling it are still unknown. Here, silk gland's histomorphology and transcriptome analysis were used to characterize changes in its structure and gene expression patterns from molt to intermolt stages. By using section staining and transmission electron microscope, a renewable cell damage was detected in the silk gland at the molt stage, and an increased number of autophagosomes and lysosomes were found in silk gland cells' cytoplasm. Next, by using RNA sequencing, 54,578,413 reads were obtained, of which 85% were mapped to the silkworm reference genome. The expression level analysis of silk protein genes and silk gland transcription factors revealed that fibroin heavy chain, fibroin light chain, P25/fhx, sericin1, sericin3 and Dimm had consistent alteration trends in temporal expression. In addition, differentially expressed genes (DEGs) were identified, and most of the DEGs associated with ecdysone signal transduction, mRNA degradation, protein proteolysis, and autophagy were significantly down-regulated in the transition from molt to intermolt, suggesting that these pathways were activated for the silk gland renewal. These findings provide insights into the molecular mechanisms of silk gland development and silk protein genes transcriptional regulation during the molt to intermolt transition process.

Comparative transcriptome analysis of Bombyx mori spinnerets and Filippi's glands suggests their role in silk fiber formation

The spinneret is located at the end of silk glands in the silkworm, Bombyx mori. The Filippi's gland (FG), which communicates with the silk gland, is considered an accessory to the latter. Although these two organs have been known for centuries, only their morphology has been studied. Their gene expression profiles and physiological roles are still unknown. Aided by next-generation sequencing, we profiled more than 11000 transcripts from the spinneret and FG of silkworm larvae on day 3 of the fifth instar (L5D3) and wandering stage (W1) in this study. A total of 59 ion-transporting protein genes and 106 cuticle protein genes were identified in the spinneret. To analyze the dynamic changes in the expression of spinneret genes, differential expression analysis was performed, and 1452 differentially expressed genes (DEGs) were identified in spinneret tissue harvested on L5D3 and W1 of the silkworm larvae. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of the DEGs revealed that the spinneret had active ion-transporting, chitin binding, and energy metabolism processes at W1. Based on these data, we hypothesized that the role of the spinneret is to provide a favorable physiological environment for the silk fiber formation. Furthermore, differential expression analysis and GO enrichment of the DEGs in the FG suggested a possible role of this gland in transporting small solutes such as ions, sugars and amino acids to the silk gland. Our findings pave a way for further functional explanation of the spinneret and FG.