Processing of snoRNAs as a new source of regulatory non-coding RNAs: snoRNA fragments form a new class of functional RNAs

Small RNAs derived from tRNAs and rRNAs are highly enriched in exosomes from both old and new world Leishmania providing evidence for conserved exosomal RNA Packaging

Background

Leishmania use exosomes to communicate with their mammalian hosts and these secreted vesicles appear to contribute to pathogenesis by delivering protein virulence factors to macrophages. In other eukaryotes, exosomes were found to carry RNA cargo, such as mRNAs and small non-coding RNAs, capable of altering recipient cell phenotype. Whether leishmania exosomes also contain RNAs which they are able to deliver to bystander cells is not known. Here, we show that leishmania exosomes indeed contain RNAs and compare and contrast the RNA content of exosomes released by Leishmania donovani and Leishmania braziliensis.

Results

We purified RNA from exosomes collected from axenic amastigote culture supernatant and found that when compared with total leishmania RNA, exosomes mainly contained short RNA sequences. Exosomes with intact membranes were capable of protecting their RNA cargo from degradation by RNase. Moreover, exosome RNA cargo was delivered to host cell cytoplasm in vitro. Sequencing of exosomal RNA indicated that the majority of cargo sequences were derived from non-coding RNA species such as rRNA and tRNA. In depth analysis revealed the presence of tRNA-derived small RNAs, a novel RNA type with suspected regulatory functions. Northern blotting confirmed the specific and selective enrichment of tRNA-derived small RNAs in exosomes. We also identified a number of novel transcripts, which appeared to be specifically enriched in exosomes compared to total cell RNA. In addition, we observed the presence of sequences mapping to siRNA-coding regions in L. braziliensis , but not in L. donovani exosomes.

Conclusions

These results show that leishmania exosomes are selectively and specifically enriched in small RNAs derived almost exclusively from non-coding RNAs. These exosomes are competent to deliver their cargo of novel, potential small regulatory RNAs to macrophages where they may influence parasite-host cell interactions. The remarkably high degree of congruence in exosomal RNA content between L. donovani and L. braziliensis, argues for the presence of a conserved mechanism for exosomal RNA packaging in leishmania. These findings open up a new avenue of research on non-canonical, small RNA pathways in this trypanosomatid, which may elucidate pathogenesis and identify novel therapeutic approaches.

Electronic supplementary material

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

Analysis of 13 cell types reveals evidence for the expression of numerous novel primate- and tissue-specific microRNAs

Two decades after the discovery of the first animal microRNA (miRNA), the number of miRNAs in animal genomes remains a vexing question. Here, we report findings from analyzing 1,323 short RNA sequencing samples (RNA-seq) from 13 different human tissue types. Using stringent thresholding criteria, we identified 3,707 statistically significant novel mature miRNAs at a false discovery rate of ≤ 0.05 arising from 3,494 novel precursors; 91.5% of these novel miRNAs were identified independently in 10 or more of the processed samples. Analysis of these novel miRNAs revealed tissue-specific dependencies and a commensurate low Jaccard similarity index in intertissue comparisons. Of these novel miRNAs, 1,657 (45%) were identified in 43 datasets that were generated by cross-linking followed by Argonaute immunoprecipitation and sequencing (Ago CLIP-seq) and represented 3 of the 13 tissues, indicating that these miRNAs are active in the RNA interference pathway. Moreover, experimental investigation through stem-loop PCR of a random collection of newly discovered miRNAs in 12 cell lines representing 5 tissues confirmed their presence and tissue dependence. Among the newly identified miRNAs are many novel miRNA clusters, new members of known miRNA clusters, previously unreported products from uncharacterized arms of miRNA precursors, and previously unrecognized paralogues of functionally important miRNA families (e.g., miR-15/107). Examination of the sequence conservation across vertebrate and invertebrate organisms showed 56.7% of the newly discovered miRNAs to be human-specific whereas the majority (94.4%) are primate lineage-specific. Our findings suggest that the repertoire of human miRNAs is far more extensive than currently represented by public repositories and that there is a significant number of lineage- and/or tissue-specific miRNAs that are uncharacterized.

Mammalian introns: when the junk generates molecular diversity

Introns represent almost half of the human genome, yet their vast majority is eliminated from eukaryotic transcripts through RNA splicing. Nevertheless, they feature key elements and functions that deserve further interest. At the level of DNA, introns are genomic segments that can shelter independent transcription units for coding and non-coding RNAs which transcription may interfere with that of the host gene, and regulatory elements that can influence gene expression and splicing itself. From the RNA perspective, some introns can be subjected to alternative splicing. Intron retention appear to provide some plasticity to the nature of the protein produced, its distribution in a given cell type and timing of its translation. Intron retention may also serve as a switch to produce coding or non-coding RNAs from the same transcription unit. Conversely, splicing of introns has been directly implicated in the production of small regulatory RNAs. Hence, splicing of introns also appears to provide plasticity to the type of RNA produced from a genetic locus (coding, non-coding, short or long). We addressed these aspects to add to our understanding of mechanisms that control the fate of introns and could be instrumental in regulating genomic output and hence cell fate.

Coding and noncoding expression patterns associated with rare obesity-related disorders: Prader-Willi and Alström syndromes

Obesity is accompanied by hyperphagia in several classical genetic obesity-related syndromes that are rare, including Prader-Willi syndrome (PWS) and Alström syndrome (ALMS). We compared coding and noncoding gene expression in adult males with PWS, ALMS, and nonsyndromic obesity relative to nonobese males using readily available lymphoblastoid cells to identify disease-specific molecular patterns and disturbed mechanisms in obesity. We found 231 genes upregulated in ALMS compared with nonobese males, but no genes were found to be upregulated in obese or PWS males and 124 genes were downregulated in ALMS. The metallothionein gene (MT1X) was significantly downregulated in ALMS, in common with obese males. Only the complex SNRPN locus was disturbed (downregulated) in PWS along with several downregulated small nucleolar RNAs (snoRNAs) in the 15q11-q13 region (SNORD116, SNORD109B, SNORD109A, SNORD107). Eleven upregulated and ten downregulated snoRNAs targeting multiple genes impacting rRNA processing, developmental pathways, and associated diseases were found in ALMS. Fifty-two miRNAs associated with multiple, overlapping gene expression disturbances were upregulated in ALMS, and four were shared with obese males but not PWS males. For example, seven passenger strand microRNAs (miRNAs) (miR-93*, miR-373*, miR-29b-2*, miR-30c-1*, miR27a*, miR27b*, and miR-149*) were disturbed in association with six separate downregulated target genes (CD68, FAM102A, MXI1, MYO1D, TP53INP1, and ZRANB1). Cell cycle (eg, PPP3CA), transcription (eg, POLE2), and development may be impacted by upregulated genes in ALMS, while downregulated genes were found to be involved with metabolic processes (eg, FABP3), immune responses (eg, IL32), and cell signaling (eg, IL1B). The high number of gene and noncoding RNA disturbances in ALMS contrast with observations in PWS and males with nonsyndromic obesity and may reflect the progressing multiorgan pathology of the ALMS disease process.

Similar response patterns to topical minoxidil foam 5% in frontal and vertex scalp of men with androgenetic alopecia: a microarray analysis

BACKGROUND

There are regional variations in the scalp hair miniaturization seen in androgenetic alopecia (AGA). Use of topical minoxidil can lead to reversal of miniaturization in the vertex scalp. However, its effects on other scalp regions have been less well studied.

OBJECTIVES

To determine whether scalp biopsies from men with AGA show variable gene expression before and after 8 weeks of treatment with minoxidil topical foam 5% (MTF) vs. placebo.

METHODS

A placebo-controlled double-blinded prospective pilot study of MTF vs. placebo was conducted in 16 healthy men aged 18-49 years with Hamilton-Norwood type IV-V thinning. The subjects were asked to apply the treatment (active drug or placebo) to the scalp twice daily for 8 weeks. Stereotactic scalp photographs were taken at the baseline and final visits, to monitor global hair growth. Scalp biopsies were taken at the leading edge of hair loss from the frontal and vertex scalp before and after treatment with MTF and placebo, and microarray analysis was performed using the Affymetrix GeneChip HG U133 Plus 2.0.

RESULTS

Global stereotactic photographs showed that MTF induced hair growth in both the frontal and vertex scalp of patients with AGA. Regional differences in gene expression profiles were observed before treatment. However, MTF treatment induced the expression of hair keratin-associated genes and decreased the expression of epidermal differentiation complex and inflammatory genes in both scalp regions.

CONCLUSIONS

These data suggest that MTF is effective in the treatment of both the frontal and vertex scalp of patients with AGA.

Aging and cardiovascular diseases: the role of gene-diet interactions

In the study of longevity, increasing importance is being placed on the concept of healthy aging rather than considering the total number of years lived. Although the concept of healthy lifespan needs to be defined better, we know that cardiovascular diseases (CVDs) are the main age-related diseases. Thus, controlling risk factors will contribute to reducing their incidence, leading to healthy lifespan. CVDs are complex diseases influenced by numerous genetic and environmental factors. Numerous gene variants that are associated with a greater or lesser risk of the different types of CVD and of intermediate phenotypes (i.e., hypercholesterolemia, hypertension, diabetes) have been successfully identified. However, despite the close link between aging and CVD, studies analyzing the genes related to human longevity have not obtained consistent results and there has been little coincidence in the genes identified in both fields. The APOE gene stands out as an exception, given that it has been identified as being relevant in CVD and longevity. This review analyzes the genomic and epigenomic factors that may contribute to this, ranging from identifying longevity genes in model organisms to the importance of gene-diet interactions (outstanding among which is the case of the TCF7L2 gene).

RNA-Seq of the nucleolus reveals abundant SNORD44-derived small RNAs

Small non-coding RNAs represent RNA species that are not translated to proteins, but which have diverse and broad functional activities in physiological and pathophysiological states. The knowledge of these small RNAs is rapidly expanding in part through the use of massive parallel (deep) sequencing efforts. We present here the first deep sequencing of small RNomes in subcellular compartments with particular emphasis on small RNAs (sRNA) associated with the nucleolus. The vast majority of the cellular, cytoplasmic and nuclear sRNAs were identified as miRNAs. In contrast, the nucleolar sRNAs had a unique size distribution consisting of 19-20 and 25 nt RNAs, which were predominantly composed of small snoRNA-derived box C/D RNAs (termed as sdRNA). Sequences from 47 sdRNAs were identified, which mapped to both 5' and 3' ends of the snoRNAs, and retained conserved box C or D motifs. SdRNA reads mapping to SNORD44 comprised 74% of all nucleolar sdRNAs, and were confirmed by Northern blotting as comprising both 20 and 25 nt RNAs. A novel 120 nt SNORD44 form was also identified. The expression of the SNORD44 sdRNA and 120 nt form was independent of Dicer/Drosha-mediated processing pathways but was dependent on the box C/D snoRNP proteins/sno-ribonucleoproteins fibrillarin and NOP58. The 120 nt SNORD44-derived RNA bound to fibrillarin suggesting that C/D sno-ribonucleoproteins are involved in regulating the stability or processing of SNORD44. This study reveals sRNA cell-compartment specific expression and the distinctive unique composition of the nucleolar sRNAs.

Multiple in vivo pathways for Escherichia coli small ribosomal subunit assembly occur on one pre-rRNA

Processing of transcribed precursor ribosomal RNA (pre-rRNA) to a mature state is a conserved aspect of ribosome biogenesis in vivo. We developed an affinity-purification system to isolate and analyze in vivo-formed pre-rRNA-containing ribonucleoprotein (RNP) particles (rRNPs) from wild-type E. coli. We observed that the first processing intermediate of pre-small subunit (pre-SSU) rRNA is a platform for biogenesis. These pre-SSU-containing RNPs have differing ribosomal-protein and auxiliary factor association and rRNA folding. Each RNP lacks the proper architecture in functional regions, thus suggesting that checkpoints preclude immature subunits from entering the translational cycle. This work offers in vivo snapshots of SSU biogenesis and reveals that multiple pathways exist for the entire SSU biogenesis process in wild-type E. coli. These findings have implications for understanding SSU biogenesis in vivo and offer a general strategy for analysis of RNP biogenesis.

Timed regulation of P-element-induced wimpy testis-interacting RNA expression during rat liver regeneration

Small noncoding RNAs comprise a growing family of molecules that regulate key cellular processes, including messenger RNA (mRNA) degradation, translational repression, and transcriptional gene silencing. P-element-induced wimpy testis (PIWI)-interacting RNAs (piRNAs) represent a class of small RNAs initially identified in the germline of a variety of species, where they contribute to maintenance of genome stability, and recently found expressed also in stem and somatic cells, where their role and responsiveness to physiopathological signals remain elusive. Here, we investigated piRNA expression in rat liver and its response to the stimuli exerted by regenerative proliferation of this organ. Quantitative polymerase chain reaction analysis identify in the liver the RNAs encoding PIWIL2/HILI, PIWIL4/HIWI2, and other components of the piRNA biogenesis pathways, suggesting that this is indeed functional. RNA sequencing before, during, and after the wave of cell proliferation that follows partial hepatectomy (PH) identified ∼1,400 mammalian germline piRNAs expressed in rat liver, including 72 showing timed changes in expression 24-48 hours post-PH, a timing that corresponds to cell transition through the S phase, returning to basal levels by 168 hours, when organ regeneration is completed and hepatocytes reach quiescence.

CONCLUSION

The piRNA pathway is active in somatic cells of the liver and is subject to regulation during the pathophysiological process of organ regeneration, when these molecules are available to exert their regulatory functions on the cell genome and transcriptome, as demonstrated by the identification of several liver mRNAs representing candidate targets of these regulatory RNAs.

miRNome of Italian Large White pig subcutaneous fat tissue: new miRNAs, isomiRs and moRNAs

Small RNAs, such as micro-RNAs (miRNAs), are decisive regulators of gene expression, and they could determine adipose tissue traits. A better knowledge of porcine fat genomics is relevant given that the pig is a biomedical model for metabolic and cardiovascular human pathologies. Adipose tissue is particularly important for the meat industry. We explored the miRNome of two adult Italian Large White pig backfat samples by Illumina RNA-Seq. Using custom bioinformatic methods, the expressed miRNAs were identified and quantified and the nucleotide sequence variability of miRNA isoforms were analysed. We detected 222 known miRNAs, 68 new miRNAs and 17 miRNA-offset RNAs (moRNAs) expressed from known hairpins, and 312 new miRNAs expressed from 253 new hairpins. Porcine transcripts targeted by the most expressed miRNAs were predicted, showing that these miRNAs may have an impact on Wnt, insulin signalling and axon guidance pathways. The expression of five small RNAs, including moRNA ssc-5'-moR-21 and a miRNA from a new hairpin, was validated by a qRT-PCR assay, thus confirming the robustness of our results. The depicted miRNome complexity suggests that quantitative and qualitative features of miRNAs and non-canonical products of their precursors are worthy of further investigation to clarify their roles in the adipose tissue biology.

How does the Mediterranean diet promote cardiovascular health? Current progress toward molecular mechanisms: gene-diet interactions at the genomic, transcriptomic, and epigenomic levels provide novel insights into new mechanisms

Epidemiological evidence supports a health-promoting effect of the Mediterranean Diet (MedDiet), especially in the prevention of cardiovascular diseases. These cardiovascular benefits have been attributed to a number of components of the MedDiet such as monounsaturated fatty acids, antioxidant vitamins and phytochemicals. However, the underlying mechanisms remain unknown. Likewise, little is known about the genes that define inter-individual variation in response to the MedDiet, although the TCF7L2 gene is emerging as an illustrative candidate for determining relative risk of cardiovascular events in response to the MedDiet. Moreover, omics technologies are providing evidence supporting potential mechanisms, some of them implicating epigenetics (i.e. microRNAs, methylation), and certain data suggest that some traditional foods could contribute via microRNAs possibly acting as exogenous regulators of gene expression. Future research should aim at increasing and consolidating the nutrigenetic and nutrigenomic knowledge of the MedDiet in order to provide sound, personalized and optimized nutritional recommendations.

Identification and initial functional characterization of a human vascular cell-enriched long noncoding RNA

OBJECTIVE

Long noncoding RNAs (lncRNAs) represent a rapidly growing class of RNA genes with functions related primarily to transcriptional and post-transcriptional control of gene expression. There is a paucity of information about lncRNA expression and function in human vascular cells. Thus, we set out to identify novel lncRNA genes in human vascular smooth muscle cells and to gain insight into their role in the control of smooth muscle cell phenotypes.

APPROACH AND RESULTS

RNA sequencing (RNA-seq) of human coronary artery smooth muscle cells revealed 31 unannotated lncRNAs, including a vascular cell-enriched lncRNA (Smooth muscle and Endothelial cell-enriched migration/differentiation-associated long NonCoding RNA [SENCR]). Strand-specific reverse transcription polymerase chain reaction (PCR) and rapid amplification of cDNA ends indicate that SENCR is transcribed antisense from the 5' end of the FLI1 gene and exists as 2 splice variants. RNA fluorescence in situ hybridization and biochemical fractionation studies demonstrate SENCR is a cytoplasmic lncRNA. Consistent with this observation, knockdown studies reveal little to no cis-acting effect of SENCR on FLI1 or neighboring gene expression. RNA-seq experiments in smooth muscle cells after SENCR knockdown disclose decreased expression of Myocardin and numerous smooth muscle contractile genes, whereas several promigratory genes are increased. Reverse transcription PCR and Western blotting experiments validate several differentially expressed genes after SENCR knockdown. Loss-of-function studies in scratch wound and Boyden chamber assays support SENCR as an inhibitor of smooth muscle cell migration.

CONCLUSIONS

SENCR is a new vascular cell-enriched, cytoplasmic lncRNA that seems to stabilize the smooth muscle cell contractile phenotype.

plantDARIO: web based quantitative and qualitative analysis of small RNA-seq data in plants

High-throughput sequencing techniques have made it possible to assay an organism's entire repertoire of small non-coding RNAs (ncRNAs) in an efficient and cost-effective manner. The moderate size of small RNA-seq datasets makes it feasible to provide free web services to the research community that provide many basic features of a small RNA-seq analysis, including quality control, read normalization, ncRNA quantification, and the prediction of putative novel ncRNAs. DARIO is one such system that so far has been focussed on animals. Here we introduce an extension of this system to plant short non-coding RNAs (sncRNAs). It includes major modifications to cope with plant-specific sncRNA processing. The current version of plantDARIO covers analyses of mapping files, small RNA-seq quality control, expression analyses of annotated sncRNAs, including the prediction of novel miRNAs and snoRNAs from unknown expressed loci and expression analyses of user-defined loci. At present Arabidopsis thaliana, Beta vulgaris, and Solanum lycopersicum are covered. The web tool links to a plant specific visualization browser to display the read distribution of the analyzed sample. The easy-to-use platform of plantDARIO quantifies RNA expression of annotated sncRNAs from different sncRNA databases together with new sncRNAs, annotated by our group. The plantDARIO website can be accessed at http://plantdario.bioinf.uni-leipzig.de/.

Links between the oncoprotein YB-1 and small non-coding RNAs in breast cancer

Background

The nucleic acid-binding protein YB-1, a member of the cold-shock domain protein family, has been implicated in the progression of breast cancer and is associated with poor patient survival. YB-1 has sequence similarity to LIN28, another cold-shock protein family member, which has a role in the regulation of small noncoding RNAs (sncRNAs) including microRNAs (miRNAs). Therefore, to investigate whether there is an association between YB-1 and sncRNAs in breast cancer, we investigated whether sncRNAs were bound by YB-1 in two breast cancer cell lines (luminal A-like and basal cell-like), and whether the abundance of sncRNAs and mRNAs changed in response to experimental reduction of YB-1 expression.

Results

RNA-immunoprecipitation with an anti-YB-1 antibody showed that several sncRNAs are bound by YB-1. Some of these were bound by YB-1 in both breast cancer cell lines; others were cell-line specific. The small RNAs bound by YB-1 were derived from various sncRNA families including miRNAs such as let-7 and miR-320, transfer RNAs, ribosomal RNAs and small nucleolar RNAs (snoRNA). Reducing YB-1 expression altered the abundance of a number of transcripts encoding miRNA biogenesis and processing proteins but did not alter the abundance of mature or precursor miRNAs.

Conclusions

YB-1 binds to specific miRNAs, snoRNAs and tRNA-derived fragments and appears to regulate the expression of miRNA biogenesis and processing machinery. We propose that some of the oncogenic effects of YB-1 in breast cancer may be mediated through its interactions with sncRNAs.

Cellular functions of the microprocessor

The microprocessor is a complex comprising the RNase III enzyme Drosha and the double-stranded RNA-binding protein DGCR8 (DiGeorge syndrome critical region 8 gene) that catalyses the nuclear step of miRNA (microRNA) biogenesis. DGCR8 recognizes the RNA substrate, whereas Drosha functions as an endonuclease. Recent global analyses of microprocessor and Dicer proteins have suggested novel functions for these components independent of their role in miRNA biogenesis. A HITS-CLIP (high-throughput sequencing of RNA isolated by cross-linking immunoprecipitation) experiment designed to identify novel substrates of the microprocessor revealed that this complex binds and regulates a large variety of cellular RNAs. The microprocessor-mediated cleavage of several classes of RNAs not only regulates transcript levels, but also modulates alternative splicing events, independently of miRNA function. Importantly, DGCR8 can also associate with other nucleases, suggesting the existence of alternative DGCR8 complexes that may regulate the fate of a subset of cellular RNAs. The aim of the present review is to provide an overview of the diverse functional roles of the microprocessor.

Molecular characterization of a patient presumed to have prader-willi syndrome

Prader-Willi syndrome (PWS) is caused by the loss of RNA expression from an imprinted region on chromosome 15 that includes SNRPN, SNORD115, and SNORD116. Currently, there are no mouse models that faithfully reflect the human phenotype and investigations rely on human post-mortem material. During molecular characterization of tissue deposited in a public brain bank from a patient diagnosed with Prader-Willi syndrome, we found RNA expression from SNRPN, SNORD115, and SNORD116 which does not support a genetic diagnosis of Prader-Willi syndrome. The patient was a female, Caucasian nursing home resident with history of morbid obesity (BMI 56.3) and mental retardation. She died at age of 56 from pulmonary embolism. SNORD115 and SNORD116 are unexpectedly stable in post mortem tissue and can be used for post-mortem diagnosis. Molecular characterization of PWS tissue donors can confirm the diagnosis and identify those patients that have been misdiagnosed.

The 5' untranslated region of the serotonin receptor 2C pre-mRNA generates miRNAs and is expressed in non-neuronal cells

The serotonin receptor 2C (HTR2C) gene encodes a G protein-coupled receptor that is exclusively expressed in neurons. Here, we report that the 5′ untranslated region of the receptor pre-mRNA as well as its hosted miRNAs is widely expressed in non-neuronal cell lines. Alternative splicing of HTR2C is regulated by MBII-52. MBII-52 and the neighboring MBII-85 cluster are absent in people with Prader–Willi syndrome, which likely causes the disease. We show that MBII-52 and MBII-85 increase expression of the HTR2C 5′ UTR and influence expression of the hosted miRNAs. The data indicate that the transcriptional unit expressing HTR2C is more complex than previously recognized and likely deregulated in Prader–Willi syndrome.