The Spatio-Temporal Expression Profiles of Silkworm Pseudogenes Provide Valuable Insights into Their Biological Roles.

Pseudogenes are sequences that have lost the ability to transcribe RNA molecules or encode truncated but possibly functional proteins. While they were once considered to be meaningless remnants of evolution, recent researches have shown that pseudogenes play important roles in various biological processes. However, the studies of pseudogenes in the silkworm, an important model organism, are limited and have focused on single or only a few specific genes. To fill these gaps, we present a systematic genome-wide studies of pseudogenes in the silkworm. We identified the pseudogenes in the silkworm using the silkworm genome assemblies, transcriptome, protein sequences from silkworm and its related species. Then we used transcriptome datasets from 832 RNA-seq analyses to construct spatio-temporal expression profiles for these pseudogenes. Additionally, we identified tissue-specifically expressed and differentially expressed pseudogenes to further understand their characteristics. Finally, the functional roles of pseudogenes as lncRNAs were systematically analyzed. We identified a total of 4410 pseudogenes, which were grouped into 4 groups, including duplications (DUPs), unitary pseudogenes (Unitary), processed pseudogenes (retropseudogenes, RETs), and fragments (FRAGs). The most of pseudogenes in the domestic silkworm were generated before the divergence of wild and domestic silkworm, however, the domestication may also involve in the accumulation of pseudogenes. These pseudogenes were clearly divided into 2 cluster, a highly expressed and a lowly expressed, and the posterior silk gland was the tissue with the most tissue-specific pseudogenes (199), implying these pseudogenes may be involved in the development and function of silkgland. We identified 3299 lncRNAs in these pseudogenes, and the target genes of these lncRNAs in silkworm pseudogenes were enriched in the egg formation and olfactory function. This study replenishes the genome annotations for silkworm, provide valuable insights into the biological roles of pseudogenes. It will also contribute to our understanding of the complex gene regulatory networks in the silkworm and will potentially have implications for other organisms as well.

Wan L ，Su S ，Liu J ，Zou B ，Jiang Y ，Jiao B ，Tang S ，Zhang Y ，Deng C ，Xiao W ... -  《Evolutionary Bioinformatics》

A Comparison of Co-expression Networks in Silk Gland Reveals the Causes of Silk Yield Increase During Silkworm Domestication.

Long-term domestication and selective breeding have increased the silk yield of the domestic silkworm (Bombyx mori) by several times the amount of the silk yield of its wild ancestor (Bombyx mandarina). However, little is known about the molecular mechanisms behind the increase in silk yield during domestication. Based on dynamic patterns of functional divergence in the silk gland between domestic and wild silkworms, we found that at early and intermediate stages of silk gland development, the up-regulated genes of the domestic silkworm were mainly involved in DNA integration, nucleic acid binding, and transporter activity, which are related to the division and growth of cells. This has led to the posterior silk gland (PSG) of the domestic silkworm having significantly more cells ("factories" of fibroin protein synthesis) than that of the wild silkworm. At the late stage of silk gland development, the up-regulated genes in the domestic silkworm was enriched in protein processing and ribosome pathways, suggesting protein synthesis efficiency is greatly improved during silkworm domestication. While there was an increase in fibroin protein synthesis, the production of sericin protein was simultaneously reduced in the silk gland of the domestic silkworm. This reflects that domestic and wild silkworms have been under different selection pressures. Importantly, we found that the network co-expressed with the silk-coding genes of the domestic silkworm was larger than that of the wild silkworm. Furthermore, many more genes co-expressed with silk-coding genes in the domestic silkworm were subjected to artificial selection than those in the wild silkworm. Our results revealed that the increase of silk yield during silkworm domestication is involved in improvement of a biological system which includes not only expansion of "factories" (cells of PSG) of protein synthesis, but also a high expression of silk-coding genes and silk production-related genes such as biological energy, transport, and ribosome pathway genes.

Zhou QZ ，Fu P ，Li SS ，Zhang CJ ，Yu QY ，Qiu CZ ，Zhang HB ，Zhang Z ... -  《Frontiers in Genetics》

Identification and comparison of long non-coding RNAs in the silk gland between domestic and wild silkworms.

Under long-term artificial selection, the domestic silkworm (Bombyx mori) has increased its silk yield tremendously in comparison with its wild progenitor, Bombyx mandarina. However, the molecular mechanism of silk yield increase is still unknown. Comparative analysis of long non-coding RNAs (lncRNAs) may provide some insights into understanding this phenotypic variation. In this study, using RNA sequencing technology data of silk gland in domestic and wild silkworms, we identified 599 lncRNAs in the silk gland of the silkworm. Compared with protein-coding genes, the silk gland lncRNA genes tend to have fewer exon numbers, shorter transcript length and lower GC-content. Moreover, we found that three lncRNA genes are significantly and differentially expressed between domestic and wild silkworms. The potential targets of two differentially expressed lncRNAs (DELs) (dw4sg_0040 and dw4sg_0483) and the expression-correlated genes with the two DELs are mainly enriched in the related processes of silk protein translation. This implies that these DELs may affect the phenotypic variation in silk yield between the domestic and wild silkworms through the post-transcriptional regulation of silk protein.

Zhou QZ ，Fang SM ，Zhang Q ，Yu QY ，Zhang Z ... -  《-》

Three vital RNA functions and interactions in the process of silk gland apoptosis in silkworm Bombyx mori.

The Silkworm Bombyx mori is an important insect in terms of economics and a model organism with a complete metamorphosis. The economic importance of silkworms is dependent on the functions of the silkgland, a specialized organ that synthesizes silk proteins. The silk gland undergoes massive degeneration during the larval to pupal stage, which involves in cell apoptosis. In this paper, high throughput sequencing was used to detect the expression of messenger RNA (mRNA), long noncoding RNA (lncRNA), and microRNA (miRNA) from silk glands of Day 3 in the fifth instar larvae (L5D3) and the spinning 36h (sp36h). We analyzed the Gene Ontology (GO) functions of target genes of the differentially expressed lncRNAs and miRNAs. We investigated the regulations of mRNA, lncRNA, and miRNA on silk gland apoptosis in L5D3 and sp36h. In total, 10,947 lncRNAs were detected in the silk gland and the index number TCONS-00021360 lncRNA may be involved in the process of apoptosis. In addition, 344 miRNAs targeted 285 mRNAs were related to the death process under GO entry. The results indicated that miRNAs play an important role in the molecular regulation of the silk gland apoptosis compared with that of lncRNAs. Finally, we screened 746 lncRNAs and 20 miRNAs that might interact with BmDredd, and drew an interaction network among them.

Chen RT ，Xiao Y ，Liu Z ，Li LL ，Lu Y ，Jiao P ，Miao YG ... -  《-》

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

Wu Y ，Cheng T ，Liu C ，Liu D ，Zhang Q ，Long R ，Zhao P ，Xia Q ... -  《PLoS One》