Dataset Information

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

ABSTRACT: 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.

SUBMITTER: Wu Y

PROVIDER: S-EPMC4714849 | biostudies-literature | 2016

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

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

Wu Yuqian Y Cheng Tingcai T Liu Chun C Liu Duolian D Zhang Quan Q Long Renwen R Zhao Ping P Xia Qingyou Q

PloS one 20160115 1

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 inte ...[more]

PMID: 26771876

Similar Datasets

Project description:BackgroundSmall non-coding RNAs (ncRNAs) are important regulators of gene expression in eukaryotes. Previously, only microRNAs (miRNAs) and piRNAs have been identified in the silkworm, Bombyx mori. Furthermore, only ncRNAs (50-500nt) of intermediate size have been systematically identified in the silkworm.ResultsHere, we performed a systematic identification and analysis of small RNAs (18-50nt) associated with the Bombyx mori argonaute2 (BmAgo2) protein. Using RIP-seq, we identified various types of small ncRNAs associated with BmAGO2. These ncRNAs showed a multimodal length distribution, with three peaks at ~20nt, ~27nt and ~33nt, which included tRNA-, transposable element (TE)-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. The tRNA-derived fragments (tRFs) were found at an extremely high abundance and accounted for 69.90% of the BmAgo2-associated small RNAs. Northern blotting confirmed that many tRFs were expressed or up-regulated only in the BmNPV-infected cells, implying that the tRFs play a prominent role by binding to BmAgo2 during BmNPV infection. Additional evidence suggested that there are potential cleavage sites on the D, anti-codon and TψC loops of the tRNAs. TE-derived small RNAs and piRNAs also accounted for a significant proportion of the BmAgo2-associated small RNAs, suggesting that BmAgo2 could be involved in the maintenance of genome stability by suppressing the activities of transposons guided by these small RNAs. Finally, Northern blotting was also used to confirm the Bombyx 5.8 s rRNA-derived small RNAs, demonstrating that various novel small RNAs exist in the silkworm.ConclusionsUsing an RIP-seq method in combination with Northern blotting, we identified various types of small RNAs associated with the BmAgo2 protein, including tRNA-, TE-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. Our findings provide new clues for future functional studies of the role of small RNAs in insect development and evolution.

Project description:The pathogenesis of Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) infection is unclear, although accumulating evidence indicates that circular RNAs (circRNAs), which act as competing endogenous RNAs or positive regulators, play important roles in regulating gene expression in eukaryotes and, thus, may play a role in BmCPV infections. To explore the expression and biological functions of circRNAs in the silkworm midgut following BmCPV infection, silkworm circRNA expression profiles of normal midgut tissue (control) and BmCPV-infected midgut tissue (test) were determined using high-through sequencing. A total of 9,753 and 7,475circRNAs were detected from the control and test samples, respectively. The two samples shared 6,085 circRNAs, while 646 and 737 circRNAs were expressed specifically in the control and test samples, respectively. A total of 3,638 circRNAs were shown to be differentially expressed, and 400 circRNAs were substantially differentially expressed with a fold-change ≥ 2.0 (p< 0.05 and a false discover rate < 0.05), of which 294 were up-regulated and 106 were down-regulated following infection. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were conducted to determine the principal functions of the substantially differentially regulated genes. circRNA-miRNA interaction networks were constructed based on a correlation analysis between the differentially expressed circRNAs and the nature of their microRNA (miRNA) binding sites. The network inferred that 13 miRNAs interacting with 193 circRNAs were among the 300 most abundant relationships. bmo-miR-3389-5p, bmo-miR-745-3p, and bmo-miR-3262 were related to 30, 34, and 34 circRNAs, respectively. circRNA_8115, circRNA_9444, circRNA_4553, circRNA_0827, and circRNA_6649 contained six, five, four, four, and four miRNA binding sites, respectively. We further found that alternative circularization of circRNAs is a common feature in silkworms and that the junction sites of many silkworm circRNAs are flanked by canonical GT/AG splicing signals. Our study is the first to show the circRNA response to virus infection. Thus, it provides a novel perspective on circRNA-miRNA interactions during BmCPV pathogenesis, and it lays the foundation for future research of the potential roles of circRNAs in BmCPV pathogenesis.

Project description:Increasing studies have shown that mature spermatozoa contain many transcripts including mRNAs and miRNAs. However, the expression profile of long non-coding RNAs (lncRNAs) in mammalian sperm has not been systematically investigated. Here, we used highly purified RNA to investigate lncRNA expression profiles in mouse mature sperm by stranded-specific RNA-seq. We identified 20,907 known and 4,088 novel lncRNAs transcripts, and the existence of intact lncRNAs was confirmed by RT-PCR and fluorescence in situ hybridization on two representative lncRNAs. Compared to round spermatids, 1,794 upregulated and 165 downregulated lncRNAs and 4,435 upregulated and 3,920 downregulated mRNAs were identified in sperm. Based on the "Cis and Trans" RNA-RNA interaction principle, we found 14,259 targeted coding genes of differently expressed lncRNAs. In terms of Gene ontology (GO) analysis, differentially expressed lncRNAs targeted genes mainly related to nucleic acid metabolic, protein modification, chromatin and histone modification, heterocycle compound metabolic, sperm function, spermatogenesis and other processes. In contrast, differentially expressed transcripts of mRNAs were highly enriched for protein metabolic process and RNA metabolic, spermatogenesis, sperm motility, cell cycle, chromatin organization, heterocycle and aromatic compound metabolic processes. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis showed that the differentially expressed lncRNAs were involved in RNA transport, mRNA surveillance pathway, PI3K-Akt signaling pathway, AMPK signaling pathway, protein processing in endoplasmic reticulum. Metabolic pathways, mRNA surveillance pathway, AMPK signaling pathway, cell cycle, RNA transport splicesome and endocytosis incorporated with the differentially expressed mRNA. Furthermore, many lncRNAs were specifically expressed in testis/sperm, and 880 lncRNAs were conserved between human and mouse. In summary, this study provides a preliminary database valuable for identifying lncRNAs critical in the late stage of spermatogenesis or important for sperm function regulation, fertilization and early embryo development.

Project description:BackgroundThe ATP-binding cassette (ABC) transporter superfamily is the largest transporter gene family responsible for transporting specific molecules across lipid membranes in all living organisms. In insects, ABC transporters not only have important functions in molecule transport, but also play roles in insecticide resistance, metabolism and development.ResultsFrom the genome of the silkworm, Bombyx mori, we have identified 51 putative ABC genes which are classified into eight subfamilies (A-H) by phylogenetic analysis. Gene duplication is very evident in the ABCC and ABCG subfamilies, whereas gene numbers and structures are well conserved in the ABCD, ABCE, ABCF, and ABCH subfamilies. Microarray analysis revealed that expression of 32 silkworm ABC genes can be detected in at least one tissue during different developmental stages, and the expression patterns of some of them were confirmed by quantitative real-time PCR. A large number of ABC genes were highly expressed in the testis compared to other tissues. One of the ABCG genes, BmABC002712, was exclusively and abundantly expressed in the Malpighian tubule implying that BmABC002712 plays a tissue-specific role. At least 5 ABCG genes, including BmABC005226, BmABC005203, BmABC005202, BmABC010555, and BmABC010557, were preferentially expressed in the midgut, showing similar developmental expression profiles to those of 20-hydroxyecdysone (20E)-response genes. 20E treatment induced the expression of these ABCG genes in the midgut and RNA interference-mediated knockdown of USP, a component of the 20E receptor, decreased their expression, indicating that these midgut-specific ABCG genes are 20E-responsive.ConclusionIn this study, a genome-wide analysis of the silkworm ABC transporters has been conducted. A comparison of ABC transporters from 5 insect species provides an overview of this vital gene superfamily in insects. Moreover, tissue- and stage-specific expression data of the silkworm ABCG genes lay a foundation for future analysis of their physiological function and hormonal regulation.

Dataset Information

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

Publications

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

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets