Experimental RNomics and genomic comparative analysis reveal a large group of species-specific small non-message RNAs in the silkworm Bombyx mori

Role of small nucleolar RNAs in alternative splicing of the doublesex gene in the silkworm, Bombyx mori

More and more evidence shows that small noncoding RNAs (ncRNAs) play diverse roles in development, stress response and other cellular processes, but functional study of intermediate-size ncRNAs is still rare. Here, the expression profile of 16 intermediate-size ncRNAs in ovary and testis of silkworm Bombyx mori were analyzed. Twelve ncRNAs, including 5 small nucleolar RNAs (snoRNAs) and 7 unclassified ncRNAs, accumulated more in the testis than in the ovary of silkworm, especially Bm-163, Bm-51 and Bm-68. Four ncRNAs (including three orphan snoRNAs and one unclassified ncRNA) had higher expression level in the ovary than in the testis, especially Bm-86. Overexpression of the testis-enriched snoRNA Bm-68 in the female led to the accumulation of male-specific isoform of doublesex (BmdsxM) and increased the expression ratio of BmdsxM: BmdsxF. While overexpression of ovary-enriched snoRNA Bm-86 in the male decreased the expression ratio of BmdsxM: BmdsxF, indicating the roles of the two snoRNAs played in the alternative splicing of Bmdsx of silkworm, which will provide new clues for the functional study of snoRNAs in insects.

Probing small ribosomal subunit RNA helix 45 acetylation across eukaryotic evolution

NAT10 is an essential enzyme that catalyzes N4-acetylcytidine (ac4C) in eukaryotic transfer RNA and 18S ribosomal RNA. Recent studies suggested that rRNA acetylation is dependent on SNORD13, a box C/D small nucleolar RNA predicted to base-pair with 18S rRNA via two antisense elements. However, the selectivity of SNORD13-dependent cytidine acetylation and its relationship to NAT10's essential function remain to be defined. Here, we demonstrate that SNORD13 is required for acetylation of a single cytidine of human and zebrafish 18S rRNA. In-depth characterization revealed that SNORD13-dependent ac4C is dispensable for human cell growth, ribosome biogenesis, translation and development. This loss of function analysis inspired a cross-evolutionary survey of the eukaryotic rRNA acetylation 'machinery' that led to the characterization of many novel metazoan SNORD13 genes. This includes an atypical SNORD13-like RNA in Drosophila melanogaster which guides ac4C to 18S rRNA helix 45 despite lacking one of the two rRNA antisense elements. Finally, we discover that Caenorhabditis elegans 18S rRNA is not acetylated despite the presence of an essential NAT10 homolog. Our findings shed light on the molecular mechanisms underlying SNORD13-mediated rRNA acetylation across eukaryotic evolution and raise new questions regarding the biological and evolutionary relevance of this highly conserved rRNA modification.

Long Noncoding RNA lncR17454 Regulates Metamorphosis of Silkworm Through let-7 miRNA Cluster

A number of long noncoding RNAs (lncRNAs) have been identified in silkworm, but little is known about their functions. Recent study showed that the let-7 miRNA cluster (contains let-7, miR-2795, and miR-100) was transcribed from the last exon of lncRNA lncR17454 in silkworm. To investigate the functional role of lncR17454, dsRNAs of lncR17454 were injected into the hemolymph of 1-d-old third-instar larvae of Bombyx mori, repression of lncR17454 led to molting arrestment during the larval-larval and larval-pupal transition of silkworm, which was consistent to the result as let-7 knockdown in other studies. The expression level of mature let-7, miR-100, and miR-2795 decreased 40%, 36%, and 40%, respectively, while the mRNA level of two predicted target genes of let-7, the Broad Complex isoform 2 (BR-C-Z2) and the BTB-Zinc finger transcription repression factor gene Abrupt (Ab), increased significantly after lncR17454 knockdown. In contrast, when adding the 20-Hydroxyecdysone (20E) to silkworm BmN4 cell lines, the expression level of lncR17454 and let-7 cluster all increased significantly, but the expression of Abrupt, the predicted target gene of let-7, was repressed. Dual-luciferase reporter assays confirmed Abrupt was the real target of let-7. Here we found that the lncRNA lncR17454 can play regulator roles in the metamorphosis of silkworm through let-7 miRNA cluster and the ecdysone signaling pathway, which will provide new clues for lepidopteran pest control.

Small nucleolar RNAs: continuing identification of novel members and increasing diversity of their molecular mechanisms of action

Identified five decades ago amongst the most abundant cellular RNAs, small nucleolar RNAs (snoRNAs) were initially described as serving as guides for the methylation and pseudouridylation of ribosomal RNA through direct base pairing. In recent years, however, increasingly powerful high-throughput genomic approaches and strategies have led to the discovery of many new members of the family and surprising diversity in snoRNA functionality and mechanisms of action. SnoRNAs are now known to target RNAs of many biotypes for a wider range of modifications, interact with diverse binding partners, compete with other binders for functional interactions, recruit diverse players to targets and affect protein function and accessibility through direct interaction. This mini-review presents the continuing characterization of the snoRNome through the identification of new snoRNA members and the discovery of their mechanisms of action, revealing a highly versatile noncoding family playing central regulatory roles and connecting the main cellular processes.

Small nucleolar RNA of silkworm can translocate from the nucleolus to the cytoplasm under abiotic stress

Small nucleolar RNAs (snoRNAs) are thought to be exclusively nuclear and guide nucleotide modifications of ribosomal RNAs. Recently, more and more evidence has suggested that the nucleolus is a stress sensor for changes in growth status and that snoRNAs may orchestrate the response to environmental stress through molecular interactions outside of the nucleus. We previously showed that a box C/D snoRNA Bm-15 had both nuclear and cytoplasmic location in BmN4 cell line of the silkworm, Bombyx mori. To further study the functional roles of Bm-15, changes in expression level and cellular location of Bm-15 were examined in BmN4 cells subjected to serum starvation and ultraviolet (UV) ray radiation. Results indicated that total RNA level of Bm-15 was unchanged after 24 h serum starvation, but exhibited 3-fold increases in the cytoplasm, and the nuclear-to-cytosolic distribution ratio was reduced from 5:1 to 2:1. Moreover, UV radiation also causes rapid decline in nuclear Bm-15 and progressive cytoplasmic accumulation with a percentage of 22% and 57% after 6 and 24 h UV radiation. UV treatment results in a dramatic decrease in Bm-15 nuclear-to-cytosolic ratio from 7:1 to 2:1 and 2:1 to 1:20 after 6 and 24 h UV radiation, respectively. We show here for the first time that box C/D snoRNAs can translocate from the nucleus to the cytoplasm under the abiotic stress of nutritional deficiency and UV radiation. The rapid translocation of snoRNAs from nucleus to cytoplasm may slow down the maturation of rRNAs and synthesis of ribosomes to enhance the stress resistance of cells.

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.

Small nucleolar RNA Sf-15 regulates proliferation and apoptosis of Spodoptera frugiperda Sf9 cells

BACKGROUND

Small nucleolar RNAs (snoRNAs) function in guiding 2'-O-methylation and pseudouridylation of ribosomal RNAs (rRNAs) and small nuclear RNAs (snRNAs). In recent years, more and more snoRNAs have been found to play novel roles in mRNA regulation, such as pre-mRNA splicing or RNA editing. In our previous study, we found a silkworm C/D box snoRNA Bm-15 can interact with Notch receptor gene in vitro. To further study the function of Bm-15, we cloned its homolog Sf-15 from Spodoptera frugiperda and investigate the function of Sf-15 in Sf9 cells.

RESULTS

We showed that knocking down of Sf-15 can inhibit the proliferation, then induce apoptosis of insect S. frugiperda Sf9 cells, but the results were reversed when Sf-15 was overexpressed. De novo sequencing of transcriptome of Sf9 cells showed that the expression of 21 apoptosis-related genes were increased upon Sf-15 repression. Further analysis showed that a Ca²⁺-induced cell death pathway gene Cn (PPP3C, the serine/threonine-protein phosphatase 2B catalytic subunit), was significantly increased upon Sf-15 depression but decreased when Sf-15 was overexpressed, which indicated that Cn might be a potential target of Sf-15.

CONCLUSIONS

We conclude that C/D box snoRNA Sf-15 can participate in apoptosis through regulating the expression of Ca²⁺-induced cell death pathway gene Cn in Sf9 cells. This is the first time that we found snoRNAs exhibiting dual functions in insect, which reveals a novel layer of ncRNA modulation in cell growth and death.

Identification and functional characterization of intermediate-size non-coding RNAs in maize

BACKGROUND

The majority of eukaryote genomes can be actively transcribed into non-coding RNAs (ncRNAs), which are functionally important in development and evolution. In the study of maize, an important crop for both humans and animals, aside from microRNAs and long non-coding RNAs, few studies have been conducted on intermediate-size ncRNAs.

RESULTS

We constructed a homogenized cDNA library of 50-500 nt RNAs in the maize inbred line Chang 7-2. Sequencing revealed 169 ncRNAs, which contained 58 known and 111 novel ncRNAs (including 70 snoRNAs, 27 snRNAs, 13 unclassified ncRNAs and one tRNA). Forty of the novel ncRNAs were specific to the Panicoideae, and 24% of them are located on sense-strand of the 5' or 3' terminus of protein coding genes on chromosome. Target site analysis found that 22 snoRNAs can guide to 38 2'-O-methylation and pseudouridylation modification sites of ribosomal RNAs and small nuclear RNAs. Expression analysis showed that 43 ncRNAs exhibited significantly altered expression in different tissues or developmental stages of maize seedlings, eight ncRNAs had tissue-specific expression and five ncRNAs were strictly accumulated in the early stage of leaf development. Further analysis showed that 3 of the 5 stage-specific ncRNAs (Zm-3, Zm-18, and Zm-73) can be highly induced under drought and salt stress, while one snoRNA Zm-8 can be repressed under PEG-simulated drought condition.

CONCLUSIONS

We provided a genome-wide identification and functional analysis of ncRNAs with a size range of 50-500 nt in maize. 111 novel ncRNAs were cloned and 40 ncRNAs were determined to be specific to Panicoideae. 43 ncRNAs changed significantly during maize development, three ncRNAs can be strongly induced under drought and salt stress, suggesting their roles in maize stress response. This work set a foundation for further study of intermediate-size ncRNAs in maize.

Landscape of RNAs in human lumbar disc degeneration

Accumulating evidence indicates noncoding RNAs (ncRNAs) fine-tune gene expression with mysterious machinery. We conducted a combination of mRNA, miRNA, circRNA, LncRNA microarray analyses on 10 adults' lumbar discs. Moreover, we performed additional global exploration on RNA interacting machinery in terms of in silico computational pipeline. Here we show the landscape of RNAs in human lumbar discs. In general, the RNA-abundant landscape comprises 14,635 mRNAs (37.93%), 2,059 miRNAs (5.34%), 18,995 LncRNAs (49.23%) and 2,894 (7.5%) circRNAs. Chromosome 1 contributes for RNA transcription at most (10%). Bi-directional transcription contributes evenly for RNA biogenesis, in terms of 5′ to 3′ and 3′ to 5′. Despite the majority of circRNAs are exonic, antisense (1.49%), intergenic (0.035%), intragenic (1.69%), and intronic (6.29%) circRNAs should not be ignored. A single miRNA could interact with a multitude of circRNAs. Notably, CDR1as or ciRS-7 harbors 66 consecutive binding sites for miR-7-5p (previous miR-7), evidencing our pipeline. The majority of binding sites are perfect-matched (78.95%). Collectively, global landscape of RNAs sheds novel insights on RNA interacting mechanisms in human intervertebral disc degeneration.

Identification of circular RNA in the Bombyx mori silk gland

Bombyx mori is an economically important holometabolous lepidopteran insect. In B. mori endogenous noncoding RNAs such as microRNAs (miRNAs) and Piwi-interacting RNAs play crucial biological functions in metamorphosis and sex determination. In addition, circular RNAs (circRNAs) have been recently identified as noncoding RNAs in most common model organisms and show potential as gene regulators. However, to date, there have been few studies on the circRNAs present in the B. mori genome conducted to date. Here, we identified 3916 circRNAs by deep circular transcriptome sequencing using the silk gland of B. mori. 3155 circRNAs were found to be derived from 1727 parental genes. The circRNAs displayed tissue-specific expression between the middle silk gland (MSG) and posterior silk gland (PSG), with 2532 and 880 being upregulated circRNAs in the MSG and PSG, respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that the parental genes from the MSG and PSG were generally annotated to similar categories and pathways. The interaction network of circRNAs and miRNAs showed that circRNAs might act as miRNA sponges or interact with miRNAs in some other way. Overall, the results revealed the complicated patterns of circRNAs in the B. mori silk gland providing a new angle from which to explore the mechanisms of complex gene regulation and efficient silk protein synthesis.

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.

The expression analysis of silk gland-enriched intermediate-size non-coding RNAs in silkworm Bombyx mori

Small non-protein coding RNAs (ncRNAs) play important roles in development, stress response and other cellular processes. Silkworm is an important model for studies on insect genetics and control of Lepidopterous pests. We have previously identified 189 novel intermediate-size ncRNAs in silkworm Bombyx mori, including 40 ncRNAs that showed altered expression in different developmental stages. Here we characterized the functions of these 40 ncRNAs by measuring their expressions in six tissues of the fifth instar larvae using Northern blot and real-time polymerase chain reaction assays. We identified nine ncRNAs (four small nucleolar RNAs and five unclassified ncRNAs) that were enriched in silk gland, including four ncRNAs that showed silk gland-specific expression. We further showed that three of nine silk gland-enriched ncRNAs were predominantly expressed in the anterior silk gland, whereas another three ncRNAs were highly accumulated in the posterior silk gland, suggesting that they may play different roles in fibroin synthesis. Furthermore, an unclassified ncRNA, Bm-152, exhibited converse expression pattern with its antisense host gene gartenzwerg in diverse tissues, and might regulate the expression of gartenzwerg through RNA-protein complex. In addition, two silk gland-enriched ncRNAs Bm-102 and Bm-159 can be found in histone modification complex, which indicated that they might play roles through epigenetic modifications. Taken together, we provided the first expression and preliminary functional analysis of silk gland-enriched ncRNAs, which will help understand the molecular mechanism of silk gland-development and fibroin synthesis.

Deep profiling of the novel intermediate-size noncoding RNAs in intraerythrocytic Plasmodium falciparum

Intermediate-size noncoding RNAs (is-ncRNAs) have been shown to play important regulatory roles in the development of several eukaryotic organisms. However, they have not been thoroughly explored in Plasmodium falciparum, which is the most virulent malaria parasite infecting human being. By using Illumina/Solexa paired-end sequencing of an is-ncRNA-specific library, we performed a systematic identification of novel is-ncRNAs in intraerythrocytic P. falciparum, strain 3D7. A total of 1,198 novel is-ncRNA candidates, including antisense, intergenic, and intronic is-ncRNAs, were identified. Bioinformatics analyses showed that the intergenic is-ncRNAs were the least conserved among different Plasmodium species, and antisense is-ncRNAs were more conserved than their sense counterparts. Twenty-two novel snoRNAs were identified, and eight potential novel classes of P. falciparum is-ncRNAs were revealed by clustering analysis. The expression of randomly selected novel is-ncRNAs was confirmed by RT-PCR and northern blotting assays. An obvious different expressional profile of the novel is-ncRNA between the early and late intraerythrocytic developmental stages of the parasite was observed. The expression levels of the antisense RNAs correlated with those of their cis-encoded sense RNA counterparts, suggesting that these is-ncRNAs are involved in the regulation of gene expression of the parasite. In conclusion, we accomplished a deep profiling analysis of novel is-ncRNAs in P. falciparum, analysed the conservation and structural features of these novel is-ncRNAs, and revealed their differential expression patterns during the development of the parasite. These findings provide important information for further functional characterisation of novel is-ncRNAs during the development of P. falciparum.

RIP-seq of BmAgo2-associated small RNAs reveal various types of small non-coding RNAs in the silkworm, Bombyx mori

Background

Small 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.

Results

Here, 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.

Conclusions

Using 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.

Experimental identification and analysis of macronuclear non-coding RNAs from the ciliate Tetrahymena thermophila

The ciliate Tetrahymena thermophila is an important eukaryotic model organism that has been used in pioneering studies of general phenomena, such as ribozymes, telomeres, chromatin structure and genome reorganization. Recent work has shown that Tetrahymena has many classes of small RNA molecules expressed during vegetative growth or sexual reorganization. In order to get an overview of medium-sized (40-500 nt) RNAs expressed from the Tetrahymena genome, we created a size-fractionated cDNA library from macronuclear RNA and analyzed 80 RNAs, most of which were previously unknown. The most abundant class was small nucleolar RNAs (snoRNAs), many of which are formed by an unusual maturation pathway. The modifications guided by the snoRNAs were analyzed bioinformatically and experimentally and many Tetrahymena-specific modifications were found, including several in an essential, but not conserved domain of ribosomal RNA. Of particular interest, we detected two methylations in the 5'-end of U6 small nuclear RNA (snRNA) that has an unusual structure in Tetrahymena. Further, we found a candidate for the first U8 outside metazoans, and an unusual U14 candidate. In addition, a number of candidates for new non-coding RNAs were characterized by expression analysis at different growth conditions.