Novel genes differentially expressed between posterior and median silk gland identified by SAGE-aided transcriptome analysis

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.

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.

Effects of Osiris9a on Silk Properties in Bombyx mori Determined by Transgenic Overexpression

Osiris is an insect-specific gene family with multiple biological roles in development, phenotypic polymorphism, and protection. In the silkworm, we have previously identified twenty-five Osiris genes with high evolutionary conservation and remarkable synteny among several insects. Bombxy moriOsiris9a (BmOsi9a) is expressed only in the silk gland, particularly in the middle silk gland (MSG). However, the biological function of BmOsi9a is still unknown. In this study, we overexpressed BmOsi9a in the silk gland by germline transgene expression. BmOsi9a was overexpressed not only in the MSG but also in the posterior silk gland (PSG). Interestingly, BmOsi9a could be secreted into the lumen in the MSG but not in the PSG. In the silk fiber, overexpressed BmOsi9a interacted with Sericin1 in the MSG, as confirmed by a co-immunoprecipitation assay. The overexpression of BmOsi9a altered the secondary structure and crystallinity of the silk fiber, thereby changing the mechanical properties. These results provide insight into the mechanisms underlying silk proteins secretion and silk fiber formation.

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.

Gene expression analysis in the larval silk gland of the eri silkworm Samia ricini

Insects produce silk for a range of purposes. In the Lepidoptera, silk is utilized as a material for cocoon production and serves to protect larvae from adverse environmental conditions or predators. Species in the Saturniidae family produce an especially wide variety of cocoons, for example, large, golden colored cocoons and those with many small holes. Although gene expression in the silk gland of the domestic silkworm (Bombyx mori L.) has been extensively studied, considerably fewer investigations have focused on members of the saturniid family. Here, we established expression sequence tags from the silk gland of the eri silkworm (Samia ricini), a saturniid species, and used these to analyze gene expression. Although we identified the fibroin heavy chain gene in the established library, genes for other major silk proteins, such as fibroin light chain and fibrohexamerin, were absent. This finding is consistent with previous reports that these latter proteins are lacking in saturniid silk. Recently, a series of fibrohexamerin-like genes were identified in the Bombyx genome. We used this information to conduct a detailed analysis of the library established here. This analysis identified putative homologues of these genes. We also found several genes encoding small silk protein molecules that are also present in the silk of other Lepidoptera. Gene expression patterns were compared between eri and domestic silkworm, and both conserved and nonconserved expression patterns were identified for the tested genes. Such differential gene expression might be one of the major causes of the differences in silk properties between these species. We believe that our study can be of value as a basic catalogue for silk gland gene expression, which will yield to the further understanding of silk evolution.

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.

Systematic Identification and Characterization of Long Non-Coding RNAs in the Silkworm, Bombyx mori

Long noncoding RNAs (lncRNAs) are emerging as important regulators in various biological processes. However, to date, no systematic characterization of lncRNAs has been reported in the silkworm Bombyx mori. In the present study, we generated eighteen RNA-seq datasets with relatively high depth. Using an in-house designed lncRNA identification pipeline, 11,810 lncRNAs were identified for 5,556 loci. Among these lncRNAs, 474 transcripts were intronic lncRNAs (ilncRNAs), 6,250 transcripts were intergenic lncRNAs (lincRNAs), and 5,086 were natural antisense lncRNAs (lncNATs). Compared with protein-coding mRNAs, silkworm lncRNAs are shorter in terms of full length but longer in terms of exon and intron length. In addition, lncRNAs exhibit a lower level of sequence conservation, more repeat sequences overlapped and higher tissue specificity than protein-coding mRNAs in the silkworm. We found that 69 lncRNA transcripts from 33 gene loci may function as miRNA precursors, and 104 lncRNA transcripts from 72 gene loci may act as competing endogenous RNAs (ceRNAs). In total, 49.47% of all gene loci (2,749/5,556) for which lncRNAs were identified showed sex-biased expression. Co-expression network analysis resulted in 19 modules, 12 of which revealed relatively high tissue specificity. The highlighted darkgoldenrod module was specifically associated with middle and posterior silk glands, and the hub lncRNAs within this module were co-expressed with proteins involved in translation, translocation, and secretory processes, suggesting that these hub lncRNAs may function as regulators of the biosynthesis, translocation, and secretion of silk proteins. This study presents the first comprehensive genome-wide analysis of silkworm lncRNAs and provides an invaluable resource for genetic, evolutionary, and genomic studies of B. mori.

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

Silk proteins are synthesized in the middle and posterior silk glands of silkworms, then transit into the anterior of the silk gland, where the silk fibers are produced, stored and processed. The mechanism of formation and spinning of the silk fibers has not been fully elucidated, and transcriptome analyses specific to the anterior silk gland have not been reported. In the present study, we explored gene expression profiles in five regions of silk gland samples using the RNA-Seq method. As a result, there were 959,979,570 raw reads obtained, of which 583,068,172 reads were mapped to the silkworm genome. A total of 7419 genes were found to be expressed in terms of reads per kilobase of exon model per million mapped reads ≥ 5 in at least one sample. The gene numbers and expression levels of the expressed genes differed between these regions. The differentially expressed genes were analyzed, and 282 genes were detected as up-regulated in the anterior silk gland, compared with the other parts. Functions of these genes were addressed using the gene ontology and Kyoto Encyclopedia of Genes and Genomes databases, and seven key pathways were enriched. It suggested that the ion transportation, energy metabolism, protease inhibitors and cuticle proteins played essential roles in the process of silk formation and spinning in the anterior silk gland. In addition, 210 genes were found differently expressed between males and females, which should help to elucidate the mechanism of the quality difference in silk fibers from male and female silkworms.

Comparative analysis of the silk gland transcriptomes between the domestic and wild silkworms

Background

Bombyx mori was domesticated from the Chinese wild silkworm, Bombyx mandarina. Wild and domestic silkworms are good models in which to investigate genes related to silk protein synthesis that may be differentially expressed in silk glands, because their silk productions are very different. Here we used the mRNA deep sequencing (RNA-seq) approach to identify the differentially expressed genes (DEGs) in the transcriptomes of the median/posterior silk glands of two domestic and two wild silkworms.

Results

The results indicated that about 58% of the total genes were expressed (reads per kilo bases per million reads (RPKM) ≥ 1) in each silkworm. Comparisons of the domestic and wild silkworm transcriptomes revealed 32 DEGs, of which 16 were up-regulated in the domestic silkworms compared with in the wild silkworms, and the other 16 were up-regulated in the wild silkworms compared with in the domestic silkworms. Quantitative real-time polymerase chain reaction (qPCR) was performed for 15 randomly selected DEGs in domestic versus wild silkworms. The qPCR results were mostly consistent with the expression levels determined from the RNA-seq data. Based on a Gene Ontology (GO) enrichment analysis and manual annotation, five of the up-regulated DEGs in the wild silkworms were predicted to be involved in immune response, and seven of the up-regulated DEGs were related to the GO term “oxidoreductase activity”, which is associated with antioxidant systems. In the domestic silkworms, the up-regulated DEGs were related mainly to tissue development, secretion of proteins and metabolism.

Conclusions

The up-regulated DEGs in the two domestic silkworms may be involved mainly in the highly efficient biosynthesis and secretion of silk proteins, while the up-regulated DEGs in the two wild silkworms may play more important roles in tolerance to pathogens and environment adaptation. Our results provide a foundation for understanding the molecular mechanisms of the silk production difference between domestic and wild silkworms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1287-9) contains supplementary material, which is available to authorized users.