Mining diverse small RNA species in the deep transcriptome

MEG3 Expression Indicates Lymph Node Metastasis and Presence of Cancer-Associated Fibroblasts in Papillary Thyroid Cancer

Papillary thyroid cancer is the most common endocrine malignancy, occurring at an incidence rate of 12.9 per 100,000 in the US adult population. While the overall 10-year survival of PTC nears 95%, the presence of lymph node metastasis (LNM) or capsular invasion indicates the need for extensive neck dissection with possible adjuvant radioactive iodine therapy. While imaging modalities such as ultrasound and CT are currently in use for the detection of suspicious cervical lymph nodes, their sensitivities for tumor-positive nodes are low. Therefore, advancements in preoperative detection of LNM may optimize the surgical and medical management of patients with thyroid cancer. To this end, we analyzed bulk RNA-sequencing datasets to identify candidate markers highly predictive of LNM. We identified MEG3, a long-noncoding RNA previously described as a tumor suppressor when expressed in malignant cells, as highly associated with LNM tissue. Furthermore, the expression of MEG3 was highly predictive of tumor infiltration with cancer-associated fibroblasts, and single-cell RNA-sequencing data revealed the expression of MEG3 was isolated to cancer-associated fibroblasts (CAFs) in the most aggressive form of thyroid cancers. Our findings suggest that MEG3 expression, specifically in CAFs, is highly associated with LNM and may be a driver of aggressive disease.

Current Status of Regulatory Non-Coding RNAs Research in the Tritryp

Trypanosomatids are protozoan parasites that cause devastating vector-borne human diseases. Gene expression regulation of these organisms depends on post-transcriptional control in responding to diverse environments while going through multiple developmental stages of their complex life cycles. In this scenario, non-coding RNAs (ncRNAs) are excellent candidates for a very efficient, quick, and economic strategy to regulate gene expression. The advent of high throughput RNA sequencing technologies show the presence and deregulation of small RNA fragments derived from canonical ncRNAs. This review seeks to depict the ncRNA landscape in trypanosomatids, focusing on the small RNA fragments derived from functional RNA molecules observed in RNA sequencing studies. Small RNA fragments derived from canonical ncRNAs (tsRNAs, snsRNAs, sdRNAs, and sdrRNAs) were identified in trypanosomatids. Some of these RNAs display changes in their levels associated with different environments and developmental stages, demanding further studies to determine their functional characterization and potential roles. Nevertheless, a comprehensive and detailed ncRNA annotation for most trypanosomatid genomes is still needed, allowing better and more extensive comparative and functional studies.

Noncoding RNAs in Papillary Thyroid Cancer: Interaction with Cancer-Associated Fibroblasts (CAFs) in the Tumor Microenvironment (TME) and Regulators of Differentiation and Lymph Node Metastasis

A large majority of all thyroid cancers are papillary thyroid carcinomas (PTC), named for the specific papillary architecture observed histologically. Despite the high rate of success with modern diagnostic and therapeutic algorithms, there are significant areas where the management of PTC can be improved. Aggressive PTC subtypes that are refractory to radioactive iodine (RAI) therapy carry a more severe prognosis and account for most of PTC-related deaths. As lymph node metastasis is present in roughly 40% of all adult PTC cases, higher specificity in these tests is a clinical need, especially since lymph node metastases are associated with reduced survival and higher recurrence rates. Additionally, this cancer can progress to more dedifferentiated and aggressive variants, such as poorly differentiated papillary thyroid cancer (PDPTC) and anaplastic thyroid cancer (ATC). Therefore, development of more sensitive and specific detection methods that allow unnecessary surgeries to be avoided is of the utmost importance. The body of large-scale, unbiased gene expression analysis in PTC has focused on the coding transcriptome, specifically mRNAs and microRNAs. However, there have been implications for the potential use of long noncoding RNAs (lncRNAs) in PTC diagnosis, prognosis, and treatment via the utilization of genome-wide studies of patient samples. lncRNAs have diverse regulatory potential in gene expression, alternative splicing, posttranscriptional mRNA modification, and epigenomic alterations. Many lncRNAs have tissue-specific expression and are demonstrated to play key roles in cancer progression and prognosis. However, lncRNAs are not being exploited as biomarkers or therapeutic targets currently, despite their elucidated effects on oncogenesis. These potent biomarkers would be revolutionary in detection at early stages, as this significantly increases the chances of survival. Their aberrant expression in cancer and correlation with steps in tumorigenesis as well as their role in differentiation would allow for a promising role as a prognostic and diagnostic biomarker in thyroid cancer. This would help prevent the more aggressive ATC that derives from dedifferentiation of the less aggressive PTC and FTC. The targeting of the specific lncRNAs could also pose a valuable treatment option via preventing or reversing this dedifferentiation process and making this usually refractory form of thyroid cancer more responsive to standard treatment options.

Transfer RNA-derived non-coding RNAs (tncRNAs): Hidden regulation of plants' transcriptional regulatory circuits

The emergence of distinct classes of non-coding RNAs has led to better insights into the eukaryotic gene regulatory networks. Amongst them, the existence of transfer RNA (tRNA)-derived non-coding RNAs (tncRNAs) demands exploration in the plant kingdom. We have designed a methodology to uncover the entire perspective of tncRNAome in plants. Using this pipeline, we have identified diverse tncRNAs with a size ranging from 14 to 50 nucleotides (nt) by utilizing 2448 small RNA-seq samples from six angiosperms, and studied their various features, including length, codon-usage, cleavage pattern, and modified tRNA nucleosides. Codon-dependent generation of tncRNAs suggests that the tRNA cleavage is highly specific rather than random tRNA degradation. The nucleotide composition analysis of tncRNA cleavage positions indicates that they are generated through precise endoribonucleolytic cleavage machinery. Certain nucleoside modifications detected on tncRNAs were found to be conserved across the plants, and hence may influence tRNA cleavage, as well as tncRNA functions. Pathway enrichment analysis revealed that common tncRNA targets are majorly enriched during metabolic and developmental processes. Further distinct tissue-specific tncRNA clusters highlight their role in plant development. Significant number of tncRNAs differentially expressed under abiotic and biotic stresses highlights their potential role in stress resistance. In summary, this study has developed a platform that will help in the understanding of tncRNAs and their involvement in growth, development, and response to various stresses. The workflow, software package, and results are freely available at http://nipgr.ac.in/tncRNA.

Identification of Abundant and Functional dodecaRNAs (doRNAs) Derived from Ribosomal RNA

Using a modified RNA-sequencing (RNA-seq) approach, we discovered a new family of unusually short RNAs mapping to ribosomal RNA 5.8S, which we named dodecaRNAs (doRNAs), according to the number of core nucleotides (12 nt) their members contain. Using a new quantitative detection method that we developed, we confirmed our RNA-seq data and determined that the minimal core doRNA sequence and its 13-nt variant C-doRNA (doRNA with a 5' Cytosine) are the two most abundant doRNAs, which, together, may outnumber microRNAs. The C-doRNA/doRNA ratio is stable within species but differed between species. doRNA and C-doRNA are mainly cytoplasmic and interact with heterogeneous nuclear ribonucleoproteins (hnRNP) A0, A1 and A2B1, but not Argonaute 2. Reporter gene activity assays suggest that C-doRNA may function as a regulator of Annexin II receptor (AXIIR) expression. doRNAs are differentially expressed in prostate cancer cells/tissues and may control cell migration. These findings suggest that unusually short RNAs may be more abundant and important than previously thought.

Small non-coding RNAs in human cancer: function, clinical utility, and characterization

Small non-coding RNAs (sncRNAs) play critical roles in multiple regulatory processes, including transcription, post-transcription, and translation. Emerging evidence reveals the critical roles of sncRNAs in cancer development and their potential role as biomarkers and/or therapeutic targets. In this paper, we review recent research on four sncRNA species with functional significance in cancer: small nucleolar RNAs, transfer RNA, small nuclear RNAs, and piwi-interacting RNAs. We introduce their functional roles in tumorigenesis and discuss the potential utility of sncRNAs as prognostic and diagnostic biomarkers and therapeutic targets. We further summarize approaches to characterize sncRNAs in a high-throughput manner, including the specific library construction and computational framework. Our review provides a perspective of the functions, clinical utility, and characterization of sncRNAs in cancer.

Alternative Applications of Genotyping Array Data Using Multivariant Methods

One of the forerunners that pioneered the revolution of high-throughput genomic technologies is the genotyping microarray technology, which can genotype millions of single-nucleotide variants simultaneously. Owing to apparent benefits, such as high speed, low cost, and high throughput, the genotyping array has gained lasting applications in genome-wide association studies (GWAS) and thus accumulated an enormous amount of data. Empowered by continuous manufactural upgrades and analytical innovation, unconventional applications of genotyping array data have emerged to address more diverse genetic problems, holding promise of boosting genetic research into human diseases through the re-mining of the rich accumulated data. Here, we review several unconventional genotyping array analysis techniques that have been built on the idea of large-scale multivariant analysis and provide empirical application examples. These unconventional outcomes of genotyping arrays include polygenic score, runs of homozygosity (ROH)/heterozygosity ratio, distant pedigree computation, and mitochondrial DNA (mtDNA) copy number inference.

Physiological RNA dynamics in RNA-Seq analysis

Physiological RNA dynamics cause problems in transcriptome analysis. Physiological RNA accumulation affects the analysis of RNA quantification, and physiological RNA degradation affects the analysis of the RNA sequence length, feature site and quantification. In the present article, we review the effects of physiological degradation and accumulation of RNA on analysing RNA sequencing data. Physiological RNA accumulation and degradation probably led to such phenomena as incorrect estimations of transcription quantification, differential expressions, co-expressions, RNA decay rates, alternative splicing, boundaries of transcription, novel genes, new single-nucleotide polymorphisms, small RNAs and gene fusion. Thus, the transcriptomic data obtained up to date warrant further scrutiny. New and improved techniques and bioinformatics software are needed to produce accurate data in transcriptome research.

Pervasive Small RNAs in Cardiometabolic Research: Great Potential Accompanied by Biological and Technical Barriers

Advances in small RNA sequencing have revealed the enormous diversity of small noncoding RNA (sRNA) classes in mammalian cells. At this point, most investigators in diabetes are aware of the success of microRNA (miRNA) research and appreciate the importance of posttranscriptional gene regulation in glycemic control. Nevertheless, miRNAs are just one of multiple classes of sRNAs and likely represent only a minor fraction of sRNA sequences in a given cell. Despite the widespread appreciation of sRNAs, very little research into non-miRNA sRNA function has been completed, likely due to some major barriers that present unique challenges for study. To emphasize the importance of sRNA research in cardiometabolic diseases, we highlight the success of miRNAs and competitive endogenous RNAs in cholesterol and glucose metabolism. Moreover, we argue that sequencing studies have demonstrated that miRNAs are just the tip of the iceberg for sRNAs. We are likely standing at the precipice of immense discovery for novel sRNA-mediated gene regulation in cardiometabolic diseases. To realize this potential, we must first address critical barriers with an open mind and refrain from viewing non-miRNA sRNA function through the lens of miRNAs, as they likely have their own set of distinct regulatory factors and functional mechanisms.

Insect genomes: progress and challenges

In the wake of constant improvements in sequencing technologies, numerous insect genomes have been sequenced. Currently, 1219 insect genome-sequencing projects have been registered with the National Center for Biotechnology Information, including 401 that have genome assemblies and 155 with an official gene set of annotated protein-coding genes. Comparative genomics analysis showed that the expansion or contraction of gene families was associated with well-studied physiological traits such as immune system, metabolic detoxification, parasitism and polyphagy in insects. Here, we summarize the progress of insect genome sequencing, with an emphasis on how this impacts research on pest control. We begin with a brief introduction to the basic concepts of genome assembly, annotation and metrics for evaluating the quality of draft assemblies. We then provide an overview of genome information for numerous insect species, highlighting examples from prominent model organisms, agricultural pests and disease vectors. We also introduce the major insect genome databases. The increasing availability of insect genomic resources is beneficial for developing alternative pest control methods. However, many opportunities remain for developing data-mining tools that make maximal use of the available insect genome resources. Although rapid progress has been achieved, many challenges remain in the field of insect genomics.

sRNA Profiling Combined With Gene Function Analysis Reveals a Lack of Evidence for Cross-Kingdom RNAi in the Wheat - Zymoseptoria tritici Pathosystem

Cross-kingdom small RNA (sRNA) silencing has recently emerged as a mechanism facilitating fungal colonization and disease development. Here we characterized RNAi pathways in Zymoseptoria tritici, a major fungal pathogen of wheat, and assessed their contribution to pathogenesis. Computational analysis of fungal sRNA and host mRNA sequencing datasets was used to define the global sRNA populations in Z. tritici and predict their mRNA targets in wheat. 389 in planta-induced sRNA loci were identified. sRNAs generated from some of these loci were predicted to target wheat mRNAs including those potentially involved in pathogen defense. However, molecular approaches failed to validate targeting of selected wheat mRNAs by fungal sRNAs. Mutant strains of Z. tritici carrying deletions of genes encoding key components of RNAi such as Dicer-like (DCL) and Argonaute (AGO) proteins were generated, and virulence bioassays suggested that these are dispensable for full infection of wheat. Nonetheless, our results did suggest the existence of non-canonical DCL-independent pathway(s) for sRNA biogenesis in Z. tritici. dsRNA targeting essential fungal genes applied in vitro or generated from an RNA virus vector in planta in a procedure known as HIGS (Host-Induced Gene Silencing) was ineffective in preventing Z. tritici growth or disease. We also demonstrated that Z. tritici is incapable of dsRNA uptake. Collectively, our data suggest that RNAi approaches for gene function analyses in this fungal species and potentially also as a control measure may not be as effective as has been demonstrated for some other plant pathogenic fungi.

Bioinformatic analysis of endogenous and exogenous small RNAs on lipoproteins

To comprehensively study extracellular small RNAs (sRNA) by sequencing (sRNA-seq), we developed a novel pipeline to overcome current limitations in analysis entitled, "Tools for Integrative Genome analysis of Extracellular sRNAs (TIGER)". To demonstrate the power of this tool, sRNA-seq was performed on mouse lipoproteins, bile, urine and livers. A key advance for the TIGER pipeline is the ability to analyse both host and non-host sRNAs at genomic, parent RNA and individual fragment levels. TIGER was able to identify approximately 60% of sRNAs on lipoproteins and >85% of sRNAs in liver, bile and urine, a significant advance compared to existing software. Moreover, TIGER facilitated the comparison of lipoprotein sRNA signatures to disparate sample types at each level using hierarchical clustering, correlations, beta-dispersions, principal coordinate analysis and permutational multivariate analysis of variance. TIGER analysis was also used to quantify distinct features of exRNAs, including 5' miRNA variants, 3' miRNA non-templated additions and parent RNA positional coverage. Results suggest that the majority of sRNAs on lipoproteins are non-host sRNAs derived from bacterial sources in the microbiome and environment, specifically rRNA-derived sRNAs from Proteobacteria. Collectively, TIGER facilitated novel discoveries of lipoprotein and biofluid sRNAs and has tremendous applicability for the field of extracellular RNA.

CleanTag Adapters Improve Small RNA Next-Generation Sequencing Library Preparation by Reducing Adapter Dimers

Next-generation small RNA sequencing is a valuable tool which is increasing our knowledge regarding small noncoding RNAs and their function in regulating genetic information. Library preparation protocols for small RNA have thus far been restricted due to higher RNA input requirements (>10 ng), long workflows, and tedious manual gel purifications. Small RNA library preparation methods focus largely on the prevention or depletion of a side product known as adapter dimer that tends to dominate the reaction. Adapter dimer is the ligation of two adapters to one another without an intervening library RNA insert or any useful sequencing information. The amplification of this side reaction is favored over the amplification of tagged library since it is shorter. The small size discrepancy between these two species makes separation and purification of the tagged library very difficult. Adapter dimer hinders the use of low input samples and the ability to automate the workflow so we introduce an improved library preparation protocol which uses chemically modified adapters (CleanTag) to significantly reduce the adapter dimer. CleanTag small RNA library preparation workflow decreases adapter dimer to allow for ultra-low input samples (down to approx. 10 pg total RNA), elimination of the gel purification step, and automation. We demonstrate how to carry out this streamlined protocol to improve NGS data quality and allow for the use of sample types with limited RNA material.

Multi-perspective quality control of Illumina RNA sequencing data analysis

Quality control (QC) is a critical step in RNA sequencing (RNA-seq). Yet, it is often ignored or conducted on a limited basis. Here, we present a multi-perspective strategy for QC of RNA-seq experiments. The QC of RNA-seq can be divided into four related stages: (1) RNA quality, (2) raw read data (FASTQ), (3) alignment and (4) gene expression. We illustrate the importance of conducting QC at each stage of an RNA-seq experiment and demonstrate our recommended RNA-seq QC strategy. Furthermore, we discuss the major and often neglected quality issues associated with the three major types of RNA-seq: mRNA, total RNA and small RNA. This RNA-seq QC overview provides comprehensive guidance for researchers who conduct RNA-seq experiments.

Target-enrichment sequencing for detailed characterization of small RNAs

Identification of important, functional small RNA (sRNA) species is currently hampered by the lack of reliable and sensitive methods to isolate and characterize them. We have developed a method, termed target-enrichment of sRNAs (TEsR), that enables targeted sequencing of rare sRNAs and diverse precursor and mature forms of sRNAs not detectable by current standard sRNA sequencing methods. It is based on the amplification of full-length sRNA molecules, production of biotinylated RNA probes, hybridization to one or multiple targeted RNAs, removal of nontargeted sRNAs and sequencing. By this approach, target sRNAs can be enriched by a factor of 500-30,000 while maintaining strand specificity. TEsR enriches for sRNAs irrespective of length or different molecular features, such as the presence or absence of a 5' cap or of secondary structures or abundance levels. Moreover, TEsR allows the detection of the complete sequence (including sequence variants, and 5' and 3' ends) of precursors, as well as intermediate and mature forms, in a quantitative manner. A well-trained molecular biologist can complete the TEsR procedure, from RNA extraction to sequencing library preparation, within 4-6 d.

The miRNAome dynamics during developmental and metabolic reprogramming of tomato root infected with potato cyst nematode

Cyst-forming plant-parasitic nematodes are pests threatening many crops. By means of their secretions cyst nematodes induce the developmental and metabolic reprogramming of host cells that lead to the formation of a syncytium, which is the sole food source for growing nematodes. The in depth micro RNA (miRNA) dynamics in the syncytia induced by Globodera rostochiensis in tomato roots was studied. The miRNAomes were obtained from syncytia covering the early and intermediate developmental stages, and were the subject of differential expression analysis. The expression of 1235 miRNAs was monitored. The fold change (log₂FC) ranged from -7.36 to 8.38, indicating that this transcriptome fraction was very variable. Moreover, we showed that the DE (differentially expressed) miRNAs do not fully overlap between the selected time points, suggesting infection stage specific regulation by miRNA. The correctness of RNA-seq expression profiling was confirmed by qRT-PCR (quantitative Real Time Polymerase Chain Reaction) for seven miRNA species. Down- and up-regulated miRNA species, including their isomiRs, were further used to identify their potential targets. Among them there are a large number of transcription factors linked to different aspects of plant development belonging to gene families, such as APETALA2 (AP2), SQUAMOSA (MADS-box), MYB, GRAS, and AUXIN RESPONSE FACTOR (ARF). The substantial portion of potential target genes belong to the NB-LRR and RLK (RECEPTOR-LIKE KINASE) families, indicating the involvement of miRNA mediated regulation in defense responses. We also collected the evidence for target cleavage in the case of 29 miRNAs using one of three alternative methods: 5' RACE (5' Rapid Amplification of cDNA Ends), a search of tasiRNA within our datasets, and the meta-analysis of tomato degradomes in the GEO (Gene Expression Omnibus) database. Eight target transcripts showed a negative correlation with their respective miRNAs at two or three time points. These results indicate a large regulatory potential for miRNAs in tuning the development and defense responses.

Data Mining of Small RNA-Seq Suggests an Association Between Prostate Cancer and Altered Abundance of 5' Transfer RNA Halves in Seminal Fluid and Prostatic Tissues

Extracellular RNAs are gaining clinical interest as biofluid-based noninvasive markers for diseases, especially cancer. In particular, derivatives of transfer RNA (tRNA) are emerging as a new class of small-noncoding RNAs with high biomarker potential. We and others previously reported alterations in serum levels of specific tRNA halves in disease states including cancer. Here, we explored seminal fluid for tRNA halves as potential markers of prostate cancer. We found that 5′ tRNA halves are abundant in seminal fluid and are elevated in prostate cancer relative to noncancer patients. Importantly, most of these tRNA halves are also detectable in prostatic tissues, and a subset were increased in malignant relative to adjacent normal tissue. These findings emphasize the potential of 5′ tRNA halves as noninvasive markers for prostate cancer screening and diagnosis and provide leads for future work to elucidate a putative role of the 5′ tRNA halves in carcinogenesis.

The discovery potential of RNA processing profiles

Small non-coding RNAs (sncRNAs) are highly abundant molecules that regulate essential cellular processes and are classified according to sequence and structure. Here we argue that read profiles from size-selected RNA sequencing capture the post-transcriptional processing specific to each RNA family, thereby providing functional information independently of sequence and structure. We developed SeRPeNT, a new computational method that exploits reproducibility across replicates and uses dynamic time-warping and density-based clustering algorithms to identify, characterize and compare sncRNAs by harnessing the power of read profiles. We applied SeRPeNT to: (i) generate an extended human annotation with 671 new sncRNAs from known classes and 131 from new potential classes, (ii) show pervasive differential processing of sncRNAs between cell compartments and (iii) predict new molecules with miRNA-like behaviour from snoRNA, tRNA and long non-coding RNA precursors, potentially dependent on the miRNA biogenesis pathway. Furthermore, we validated experimentally four predicted novel non-coding RNAs: a miRNA, a snoRNA-derived miRNA, a processed tRNA and a new uncharacterized sncRNA. SeRPeNT facilitates fast and accurate discovery and characterization of sncRNAs at an unprecedented scale. SeRPeNT code is available under the MIT license at https://github.com/comprna/SeRPeNT.

The Landscape of Small Non-Coding RNAs in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is an operational term for breast cancers lacking targetable estrogen receptor expression and HER2 amplifications. TNBC is, therefore, inherently heterogeneous, and is associated with worse prognosis, greater rates of metastasis, and earlier onset. TNBC displays mutational and transcriptional diversity, and distinct mRNA transcriptional subtypes exhibiting unique biology. High-throughput sequencing has extended cancer research far beyond protein coding regions that include non-coding small RNAs, such as miRNA, isomiR, tRNA, snoRNAs, snRNA, yRNA, 7SL, and 7SK. In this study, we performed small RNA profiling of 26 TNBC cell lines, and compared the abundance of non-coding RNAs among the transcriptional subtypes of triple negative breast cancer. We also examined their co-expression pattern with corresponding mRNAs. This study provides a detailed description of small RNA expression in triple-negative breast cancer cell lines that can aid in the development of future biomarker and novel targeted therapies.

Shortcuts to a functional adipose tissue: The role of small non-coding RNAs

Metabolic diseases such as type 2 diabetes are a major public health issue worldwide. These diseases are often linked to a dysfunctional adipose tissue. Fat is a large, heterogenic, pleiotropic and rather complex tissue. It is found in virtually all cavities of the human body, shows unique plasticity among tissues, and harbors many cell types in addition to its main functional unit - the adipocyte. Adipose tissue function varies depending on the localization of the fat depot, the cell composition of the tissue and the energy status of the organism. While the white adipose tissue (WAT) serves as the main site for triglyceride storage and acts as an important endocrine organ, the brown adipose tissue (BAT) is responsible for thermogenesis. Beige adipocytes can also appear in WAT depots to sustain heat production upon certain conditions, and it is becoming clear that adipose tissue depots can switch phenotypes depending on cell autonomous and non-autonomous stimuli. To maintain such degree of plasticity and respond adequately to changes in the energy balance, three basic processes need to be properly functioning in the adipose tissue: i) adipogenesis and adipocyte turnover, ii) metabolism, and iii) signaling. Here we review the fundamental role of small non-coding RNAs (sncRNAs) in these processes, with focus on microRNAs, and demonstrate their importance in adipose tissue function and whole body metabolic control in mammals.

The Utilization of Formalin Fixed-Paraffin-Embedded Specimens in High Throughput Genomic Studies

High throughput genomic assays empower us to study the entire human genome in short time with reasonable cost. Formalin fixed-paraffin-embedded (FFPE) tissue processing remains the most economical approach for longitudinal tissue specimen storage. Therefore, the ability to apply high throughput genomic applications to FFPE specimens can expand clinical assays and discovery. Many studies have measured the accuracy and repeatability of data generated from FFPE specimens using high throughput genomic assays. Together, these studies demonstrate feasibility and provide crucial guidance for future studies using FFPE specimens. Here, we summarize the findings of these studies and discuss the limitations of high throughput data generated from FFPE specimens across several platforms that include microarray, high throughput sequencing, and NanoString.

MicroRNAs and tRNA-derived fragments predict the transformation of myelodysplastic syndromes to acute myeloid leukemia

Myelodysplastic syndromes (MDS) are clonal hematopoietic disorders of the elderly that carry an increased risk of progression to acute myeloid leukemia (AML). Since small non-coding RNAs (sRNAs), including microRNA (miRNAs), act as regulators of cellular differentiation, we hypothesized that changes to sRNAs might be implicated in the progression of MDS to AML. We conducted sRNA sequencing on three sets of patients: Group A (MDS patients who never progressed to AML); Group B (MDS patients who later progressed to an AML); and Group C (AML patients with myelodysplasia-related changes, including patients with a known preceding diagnosis of MDS). We identified five miRNAs that differentiated Groups A and B, independent of bone marrow blast percentage, including three members of the miR-181 family, as well as differential patterns of miRNA isoforms (isomiRs) and tDRs. Thus, we have identified sRNA biomarkers that predict MDS cases that are likely to progress to AML.

Comprehensive evaluation of extracellular small RNA isolation methods from serum in high throughput sequencing

Background

DNA and RNA fractions from whole blood, serum and plasma are increasingly popular analytes that are currently under investigation for their utility in the diagnosis and staging of disease. Small non-coding ribonucleic acids (sRNAs), specifically microRNAs (miRNAs) and their variant isoforms (isomiRs), and transfer RNA (tRNA)-derived small RNAs (tDRs) comprise a repertoire of molecules particularly promising in this regard.

Results

In this designed study, we compared the performance of various methods and kits for isolating circulating extracellular sRNAs (ex-sRNAs). ex-sRNAs from one healthy individual were isolated using five different isolation kits: Qiagen Circulating Nucleic Acid Kit, ThermoFisher Scientific Ambion TRIzol LS Reagent, Qiagen miRNEasy, QiaSymphony RNA extraction kit and the Exiqon MiRCURY RNA Isolation Kit. Each isolation method was repeated four times. A total of 20 small RNA sequencing (sRNAseq) libraries were constructed, sequenced and compared using a rigorous bioinformatics approach. The Circulating Nucleic Acid Kit had the greatest miRNA isolation variability, but had the lowest isolation variability for other RNA classes (isomiRs, tDRs, and other miscellaneous sRNAs (osRNA). However, the Circulating Nucleic Acid Kit consistently generated the fewest number of reads mapped to the genome, as compared to the best-performing method, Ambion TRIzol, which mapped 10% of the miRNAs, 7.2% of the tDRs and 23.1% of the osRNAs. The other methods performed intermediary, with QiaSymphony mapping 14% of the osRNAs, and miRNEasy mapping 4.6% of the tDRs and 2.9% of the miRNAs, achieving the second best kit performance rating overall.

Conclusions

In summary, each isolation kit displayed different performance characteristics that could be construed as biased or advantageous, depending upon the downstream application and number of samples that require processing.

Electronic supplementary material

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

YM500v3: a database for small RNA sequencing in human cancer research

We previously presented the YM500 database, which contains >8000 small RNA sequencing (smRNA-seq) data sets and integrated analysis results for various cancer miRNome studies. In the updated YM500v3 database (http://ngs.ym.edu.tw/ym500/) presented herein, we not only focus on miRNAs but also on other functional small non-coding RNAs (sncRNAs), such as PIWI-interacting RNAs (piRNAs), tRNA-derived fragments (tRFs), small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs). There is growing knowledge of the role of sncRNAs in gene regulation and tumorigenesis. We have also incorporated >10 000 cancer-related RNA-seq and >3000 more smRNA-seq data sets into the YM500v3 database. Furthermore, there are two main new sections, ‘Survival' and ‘Cancer', in this updated version. The ‘Survival’ section provides the survival analysis results in all cancer types or in a user-defined group of samples for a specific sncRNA. The ‘Cancer’ section provides the results of differential expression analyses, miRNA–gene interactions and cancer miRNA-related pathways. In the ‘Expression’ section, sncRNA expression profiles across cancer and sample types are newly provided. Cancer-related sncRNAs hold potential for both biotech applications and basic research.

Isolation of High-density Lipoproteins for Non-coding Small RNA Quantification

The diversity of small non-coding RNAs (sRNA) is rapidly expanding and their roles in biological processes, including gene regulation, are emerging. Most interestingly, sRNAs are also found outside of cells and are stably present in all biological fluids. As such, extracellular sRNAs represent a novel class of disease biomarkers and are likely involved in cell signaling and intercellular communication networks. To assess their potential as biomarkers, sRNAs can be quantified in plasma, urine, and other fluids. Nevertheless, to fully understand the impact of extracellular sRNAs as endocrine signals, it is important to determine which carriers are transporting and protecting them in biological fluids (e.g., plasma), which cells and tissues contribute to extracellular sRNA pools, and cells and tissues capable of accepting and utilizing extracellular sRNA. To accomplish these goals, it is critical to isolate highly pure populations of extracellular carriers for sRNA profiling and quantification. We have previously demonstrated that lipoproteins, particularly high-density lipoproteins (HDL), transport functional microRNAs (miRNA) between cells and HDL-miRNAs are significantly altered in disease. Here, we detail a new protocol that utilizes tandem HDL isolation with density-gradient ultracentrifugation (DGUC) and fast-protein-liquid chromatography (FPLC) to obtain highly pure HDL for downstream profiling and quantification of all sRNAs, including miRNAs, using both high-throughput sequencing and real-time PCR approaches. This protocol will be a valuable resource for the investigation of sRNAs on HDL.

Small RNA Library Preparation Method for Next-Generation Sequencing Using Chemical Modifications to Prevent Adapter Dimer Formation

For most sample types, the automation of RNA and DNA sample preparation workflows enables high throughput next-generation sequencing (NGS) library preparation. Greater adoption of small RNA (sRNA) sequencing has been hindered by high sample input requirements and inherent ligation side products formed during library preparation. These side products, known as adapter dimer, are very similar in size to the tagged library. Most sRNA library preparation strategies thus employ a gel purification step to isolate tagged library from adapter dimer contaminants. At very low sample inputs, adapter dimer side products dominate the reaction and limit the sensitivity of this technique. Here we address the need for improved specificity of sRNA library preparation workflows with a novel library preparation approach that uses modified adapters to suppress adapter dimer formation. This workflow allows for lower sample inputs and elimination of the gel purification step, which in turn allows for an automatable sRNA library preparation protocol.

lncRNA expression in the auditory forebrain during postnatal development

The biological processes governing brain development and maturation depend on complex patterns of gene and protein expression, which can be influenced by many factors. One of the most overlooked is the long noncoding class of RNAs (lncRNAs), which are known to play important regulatory roles in an array of biological processes. Little is known about the distribution of lncRNAs in the sensory systems of the brain, and how lncRNAs interact with other mechanisms to guide the development of these systems. In this study, we profiled lncRNA expression in the mouse auditory forebrain during postnatal development at time points before and after the onset of hearing (P7, P14, P21, adult). First, we generated lncRNA profiles of the primary auditory cortex (A1) and medial geniculate body (MG) at each age. Then, we determined the differential patterns of expression by brain region and age. These analyses revealed that the lncRNA expression profile was distinct between both brain regions and between each postnatal age, indicating spatial and temporal specificity during maturation of the auditory forebrain. Next, we explored potential interactions between functionally-related lncRNAs, protein coding RNAs (pcRNAs), and associated proteins. The maturational trajectories (P7 to adult) of many lncRNA - pcRNA pairs were highly correlated, and predictive analyses revealed that lncRNA-protein interactions tended to be strong. A user-friendly database was constructed to facilitate inspection of the expression levels and maturational trajectories for any lncRNA or pcRNA in the database. Overall, this study provides an in-depth summary of lncRNA expression in the developing auditory forebrain and a broad-based foundation for future exploration of lncRNA function during brain development.

What is an RNA? A top layer for RNA classification

Every ribonucleic acid begins its cellular life as a transcript. If the transcript or its processing product has a function it should be regarded an RNA. Nonfunctional transcripts, by-products from processing, degradation intermediates, even those originating from (functional) RNAs, and non-functional products of transcriptional gene regulation accomplished via the act of transcription, as well as stochastic (co)transcripts could simply be addressed as transcripts (class 0). The copious functional RNAs (class I), often maturing after one or more processing steps, already are systematized into ever expanding sub-classifications ranging from micro RNAs to rRNAs. Established sub-classifications addressing a wide functional diversity remain unaffected. mRNAs (class II) are distinct from any other RNA by virtue of their potential to be translated into (poly)peptide(s) on ribosomes. We are not proposing a novel RNA classification, but wish to add a basic concept with existing terminology (transcript, RNA, and mRNA) that should serve as an additional framework for carefully delineating RNA function from an avalanche of RNA sequencing data. At the same time, this top level hierarchical model should illuminate important principles of RNA evolution and biology thus heightening our awareness that in biology boundaries and categorizations are typically fuzzy.

Small RNA Modifications: Integral to Function and Disease

Small RNAs have the potential to store a secondary layer of labile biological information in the form of modified nucleotides. Emerging evidence has shown that small RNAs including microRNAs (miRNAs), PIWI-interacting RNAs (piRNAs) and tRNA-derived small RNAs (tsRNAs) harbor a diversity of RNA modifications. These findings highlight the importance of RNA modifications in the modulation of basic properties such as RNA stability and other complex physiological processes involved in stress responses, metabolism, immunity, and epigenetic inheritance of environmentally acquired traits, among others. High-resolution, high-throughput methods for detecting, mapping and screening these small RNA modifications now provide opportunities to uncover their diagnostic potential as sensitive disease markers.

A micro-RNA expression signature for human NAFLD progression

BACKGROUND

The spectrum of nonalcoholic fatty liver disease (NAFLD) describes disease conditions deteriorating from nonalcoholic fatty liver (NAFL) to nonalcoholic steatohepatitis (NASH) to cirrhosis (CIR) to hepatocellular carcinoma (HCC). From a molecular and biochemical perspective, our understanding of the etiology of this disease is limited by the broad spectrum of disease presentations, the lack of a thorough understanding of the factors contributing to disease susceptibility, and ethical concerns related to repeat sampling of the liver. To better understand the factors associated with disease progression, we investigated by next-generation RNA sequencing the altered expression of microRNAs (miRNAs) in liver biopsies of class III obese subjects (body mass index ≥40 kg/m(2)) biopsied at the time of elective bariatric surgery.

METHODS

Clinical characteristics and unbiased RNA expression profiles for 233 miRs, 313 transfer RNAs (tRNAs), and 392 miscellaneous small RNAs (snoRNAs, snRNAs, rRNAs) were compared among 36 liver biopsy specimens stratified by disease severity.

RESULTS

The abundances of 3 miRNAs that were found to be differentially regulated (miR-301a-3p and miR-34a-5p increased and miR-375 decreased) with disease progression were validated by RT-PCR. No tRNAs or miscellaneous RNAs were found to be associated with disease severity. Similar patterns of increased miR-301a and decreased miR-375 expression were observed in 134 hepatocellular carcinoma (HCC) samples deposited in The Cancer Genome Atlas (TCGA).

CONCLUSIONS

Our analytical results suggest that NAFLD severity is associated with a specific pattern of altered hepatic microRNA expression that may drive the hallmark of this disorder: altered lipid and carbohydrate metabolism. The three identified miRNAs can potentially be used as biomarkers to access the severity of NAFLD. The persistence of this miRNA expression pattern in an external validation cohort of HCC samples suggests that specific microRNA expression patterns may permit and/or sustain NAFLD development to HCC.

RNA Sequencing of Formalin-Fixed, Paraffin-Embedded Specimens for Gene Expression Quantification and Data Mining

Background. Proper rRNA depletion is crucial for the successful utilization of FFPE specimens when studying gene expression. We performed a study to evaluate two major rRNA depletion methods: Ribo-Zero and RNase H. RNAs extracted from 4 samples were treated with the two rRNA depletion methods in duplicate and sequenced (N = 16). We evaluated their reducibility, ability to detect RNA, and ability to molecularly subtype these triple negative breast cancer specimens. Results. Both rRNA depletion methods produced consistent data between the technical replicates. We found that the RNase H method produced higher quality RNAseq data as compared to the Ribo-Zero method. In addition, we evaluated the RNAseq data generated from the FFPE tissue samples for noncoding RNA, including lncRNA, enhancer/super enhancer RNA, and single nucleotide variation (SNV). We found that the RNase H is more suitable for detecting high-quality, noncoding RNAs as compared to the Ribo-Zero and provided more consistent molecular subtype identification between replicates. Unfortunately, neither method produced reliable SNV data. Conclusions. In conclusion, for FFPE specimens, the RNase H rRNA depletion method performed better than the Ribo-Zero. Neither method generates data sufficient for SNV detection.

Practicability of mitochondrial heteroplasmy detection through an Illumina genotyping array

MOTIVATION

High-throughput genomic data often contain unexpected information that can be mined for alternative applications. Despite the rise of high-throughput sequencing, Illumina genotyping arrays remain a driving force in large scale genetic and epidemiology studies. By processing and analyzing genotyping data of over 100,000 samples genotyped on Illumina genotyping arrays, we discovered evidence that indicates that mitochondrial heteroplasmy can be estimated from the fluorescence intensity data of the array. To verify our hypothesis, we conducted a sequencing validation study.

RESULT

Mitochondrial DNA targeted sequencing was performed on three samples that had been genotyped using the Illumina exome genotyping array. In each sample chosen, one heteroplasmy target was identified from the genotyping array, and sequencing data verified all three putative heteroplasmic sites. The estimated heteroplasmy level difference between that estimated from the genotyping fluorescence intensity and that directly measured from sequencing was 3.2% on average. Our analysis showed that an Illumina genotyping array can accurately and reliably estimate high-level heteroplasmy (>40%); however, intensity data from a genotyping array is not suitable for estimating low level heteroplasmy (<25%).

HDL and microRNA therapeutics in cardiovascular disease

microRNAs (miRNA) are small non-coding RNAs (sRNA) that post-transcriptionally regulate gene (mRNA) expression and are implicated in many biological processes and diseases. Many miRNAs have been reported to be altered in cardiovascular disease (CVD); both cellular and extracellular miRNA levels are affected by hypercholesterolemia and atherosclerosis. We and other groups have reported that lipoproteins transport miRNAs in circulation and these lipoprotein signatures are significantly altered in hypercholesterolemia and coronary artery disease (CAD). Extracellular miRNAs are a new class of potential biomarkers for CVD; however, they may also be new drug targets as high-density lipoproteins (HDL) transfer functional miRNAs to recipient cells in an endocrine-like form of intercellular communication that likely suppresses vascular inflammation. Recently, RNA-based drugs have emerged as the next frontier in drug therapy, and there are many miRNA inhibitors and mimics in clinical development. Here, we discuss specific miRNA drug targets and how their manipulation may impact CVD. We also address the potential for manipulating HDL-miRNA levels to treat CVD and the use of HDL as a delivery vehicle for RNA and chemical drugs. Finally, we outline the current and future challenges for HDL and miRNA-based therapeutics for the prevention and treatment of CVD.

Sources and Functions of Extracellular Small RNAs in Human Circulation

Various biotypes of endogenous small RNAs (sRNAs) have been detected in human circulation, including microRNAs, transfer RNAs, ribosomal RNA, and yRNA fragments. These extracellular sRNAs (ex-sRNAs) are packaged and secreted by many different cell types. Ex-sRNAs exhibit differences in abundance in several disease states and have, therefore, been proposed for use as effective biomarkers. Furthermore, exosome-borne ex-sRNAs have been reported to elicit physiological responses in acceptor cells. Exogenous ex-sRNAs derived from diet (most prominently from plants) and microorganisms have also been reported in human blood. Essential issues that remain to be conclusively addressed concern the (a) presence and sources of exogenous ex-sRNAs in human bodily fluids, (b) detection and measurement of ex-sRNAs in human circulation, (c) selectivity of ex-sRNA export and import, (d) sensitivity and specificity of ex-sRNA delivery to cellular targets, and (e) cell-, tissue-, organ-, and organism-wide impacts of ex-sRNA-mediated cell-to-cell communication. We survey the present state of knowledge of most of these issues in this review.

Comparative transcriptomic analysis of human and Drosophila extracellular vesicles

Extracellular vesicles (EVs) are membrane-enclosed nanoparticles containing specific repertoires of genetic material. In mammals, EVs can mediate the horizontal transfer of various cargos and signaling molecules, notably miRNA and mRNA species. Whether this form of intercellular communication prevails in other metazoans remains unclear. Here, we report the first parallel comparative morphologic and transcriptomic characterization of EVs from Drosophila and human cellular models. Electronic microscopy revealed that human and Drosophila cells release similar EVs with diameters ranging from 30 to 200 nm, which contain complex populations of transcripts. RNA-seq identified abundant ribosomal RNAs, related pseudogenes and retrotransposons in human and Drosophila EVs. Vault RNAs and Y RNAs abounded in human samples, whereas small nucleolar RNAs involved in pseudouridylation were most prevalent in Drosophila EVs. Numerous mRNAs were identified, largely consisting of exonic sequences displaying full-length read coverage and enriched for translation and electronic transport chain functions. By analogy with human systems, these sizeable similarities suggest that EVs could potentially enable RNA-mediated intercellular communication in Drosophila.

Nutrigenomics, the Microbiome, and Gene-Environment Interactions: New Directions in Cardiovascular Disease Research, Prevention, and Treatment: A Scientific Statement From the American Heart Association

Cardiometabolic diseases are the leading cause of death worldwide and are strongly linked to both genetic and nutritional factors. The field of nutrigenomics encompasses multiple approaches aimed at understanding the effects of diet on health or disease development, including nutrigenetic studies investigating the relationship between genetic variants and diet in modulating cardiometabolic risk, as well as the effects of dietary components on multiple "omic" measures, including transcriptomics, metabolomics, proteomics, lipidomics, epigenetic modifications, and the microbiome. Here, we describe the current state of the field of nutrigenomics with respect to cardiometabolic disease research and outline a direction for the integration of multiple omics techniques in future nutrigenomic studies aimed at understanding mechanisms and developing new therapeutic options for cardiometabolic disease treatment and prevention.

Identification of Differentially Expressed microRNAs between the Fenpropathrin Resistant and Susceptible Strains in Tetranychus cinnabarinus

The carmine spider mite (Tetranychus cinnabarinus) is one of the most serious pests on crops and its control mainly depends on chemical acaricides. The excessive and improper acaricides use has resulted in mite resistance to many acaricides, including fenpropathrin. Previous studies have indicated fenpropathrin resistance is a complex development process involving many genes, but information on resistance mechanism of post-transcription regulation is rare. Using Illumina sequencing, several categories of sRNAs were identified from susceptible (TS) and fenpropathrin-resistant strains (TR) of T. cinnabarinus, including 75 known microRNAs (miRNAs) and 64 novel miRNAs, whose target genes containing 78592 miRNA-target pairs were predicted by 6 algorithms. Also, 12 significantly differently expressed miRNAs were identified between the TS and TR libraries and RT-qPCR validation also performed a well consistency with sequencing. The targets of significantly differentially expressed miRNAs included 7 glutathione S-transferase, 7 cytochrome P450 and 16 carboxyl/choline esterase genes, their function in fenpropathrin resistance were further analyzed. The present study provides the firstly large-scale characterization of miRNAs in T. cinnabarinus and the comparison between TS and TR strains gives a clue on how miRNA involves in fenpropathrin resistance.

Mitochondria sequence mapping strategies and practicability of mitochondria variant detection from exome and RNA sequencing data

The rapid progress in high-throughput sequencing has significantly enriched our capacity for studying the mitochondrial DNA (mtDNA). In addition to performing specific mitochondrial targeted sequencing, an increasingly popular alternative approach is using the off-target reads from exome sequencing to infer mtDNA variants, including single nucleotide polymorphisms (SNPs) and heteroplasmy. However, the effectiveness and practicality of this approach have not been tested. Recently, RNAseq data have also been suggested as a good source for alternative data mining, but whether mitochondrial variants can be detected from RNAseq data has not been validated. We designed a study to evaluate the practicability of mtDNA variant detection using exome and RNA sequencing data. Five breast cancer cell lines were sequenced through mitochondrial targeted, exome, and RNA sequencing. Mitochondrial targeted sequencing was used as the gold standard to compute the validation and false discovery rates of SNP and heteroplasmy detection in exome and RNAseq data. We found that exome and RNA sequencing can accurately detect mitochondrial SNPs. However, the lower false discovery rate makes exome sequencing a better choice for heteroplasmy detection than RNAseq. Furthermore, we examined three alignment strategies and found that aligning reads directly to the mitochondrial reference genome or aligning reads to the nuclear and mitochondrial references genomes simultaneously produced the best results, and that aligning to the nuclear genome first and afterwards to the mitochondrial genome performed poorly. In conclusion, our study provides important guidelines for future studies that intend to use either exome sequencing or RNAseq data to infer mitochondrial SNPs and heteroplasmy.

Piwi-interacting RNAs in cancer: emerging functions and clinical utility

PIWI-interacting RNAs (piRNAs) are emerging players in cancer genomics. Originally described in the germline, there are over 20,000 piRNA genes in the human genome. In contrast to microRNAs, piRNAs interact with PIWI proteins, another member of the Argonaute family, and function primarily in the nucleus. There, they are involved in the epigenetic silencing of transposable elements in addition to the transcriptional regulation of genes. It has recently been demonstrated that piRNAs are also expressed across a variety of human somatic tissue types in a tissue-specific manner. An increasing number of studies have shown that aberrant piRNA expression is a signature feature across multiple tumour types; however, their specific tumorigenic functions remain unclear. In this article, we discuss the emerging functional roles of piRNAs in a variety of cancers, and highlight their potential clinical utilities.

tRNA-Derived Fragments (tRFs): Emerging New Roles for an Ancient RNA in the Regulation of Gene Expression

This review will summarise the recent discoveries and current state of research on short noncoding RNAs derived from tRNAs—known as tRNA-derived fragments (tRFs). It will describe the features of the known subtypes of these RNAs; including sequence characteristics, protein interactors, expression characteristics, biogenesis, and similarity to canonical miRNA pathways. Also their role in regulating gene expression; including mediating translational suppression, will be discussed. We also highlight their potential use as biomarkers, functions in gene regulation and links to disease. Finally, this review will speculate as to the origin and rationale for the conservation of this novel class of noncoding RNAs amongst both prokaryotes and eukaryotes.

Circulating Extracellular microRNA in Systemic Autoimmunity

MicroRNAs (miRNAs) are differentially regulated in healthy, activated, inflamed, neoplastic, or otherwise pathological cells and tissues. While their main functions are executed intracellularly, many miRNAs can reproducibly be detected extracellularly in plasma and serum. This circulating, extracellular miRNA is protected against degradation by complexation with carrier proteins and/or by being enclosed in subcellular membrane vesicles. This, together with their tissue- and disease-specific expression, has fuelled the interest in using circulating microRNA profiles as harbingers of disease, i.e., as diagnostic analytes and as clues to dysregulated pathways in disease. Many studies show that inflammation and immune dysregulation, e.g., in autoimmune diseases, are associated with distinct miRNA expression changes in targeted tissues and in innate and adaptive immunity cells such as lymphocytes, natural killer cells, neutrophil granulocytes, and monocyte-macrophages. Exploratory studies (only validated in a few cases) also show that specific profiles of circulating miRNAs are associated with different systemic autoimmune diseases including systemic lupus erythematosus (SLE), systemic sclerosis, and rheumatoid arthritis. Even though the link between cellular alterations and extracellular profiles is still unpredictable, the data suggest that circulating miRNAs in autoimmunity may become diagnostically useful. Here, we review important circulating miRNAs in animal models of inflammation and in systemic autoimmunity and summarize some proposed functions of miRNAs in immune regulation and dysregulation. We conclude that the studies suggest new hypotheses and additional experiments, and that further diagnostic development is highly dependent on analytical method development and on obtaining sufficient numbers of uniformly processed samples from clinically well-characterized patients and controls.

RNAseq by Total RNA Library Identifies Additional RNAs Compared to Poly(A) RNA Library

The most popular RNA library used for RNA sequencing is the poly(A) captured RNA library. This library captures RNA based on the presence of poly(A) tails at the 3' end. Another type of RNA library for RNA sequencing is the total RNA library which differs from the poly(A) library by capture method and price. The total RNA library costs more and its capture of RNA is not dependent on the presence of poly(A) tails. In practice, only ribosomal RNAs and small RNAs are washed out in the total RNA library preparation. To evaluate the ability of detecting RNA for both RNA libraries we designed a study using RNA sequencing data of the same two breast cancer cell lines from both RNA libraries. We found that the RNA expression values captured by both RNA libraries were highly correlated. However, the number of RNAs captured was significantly higher for the total RNA library. Furthermore, we identify several subsets of protein coding RNAs that were not captured efficiently by the poly(A) library. One of the most noticeable is the histone-encode genes, which lack the poly(A) tail.

Transfer RNA detection by small RNA deep sequencing and disease association with myelodysplastic syndromes

Background

Although advances in sequencing technologies have popularized the use of microRNA (miRNA) sequencing (miRNA-seq) for the quantification of miRNA expression, questions remain concerning the optimal methodologies for analysis and utilization of the data. The construction of a miRNA sequencing library selects RNA by length rather than type. However, as we have previously described, miRNAs represent only a subset of the species obtained by size selection. Consequently, the libraries obtained for miRNA sequencing also contain a variety of additional species of small RNAs. This study looks at the prevalence of these other species obtained from bone marrow aspirate specimens and explores the predictive value of these small RNAs in the determination of response to therapy in myelodysplastic syndromes (MDS).

Methods

Paired pre and post treatment bone marrow aspirate specimens were obtained from patients with MDS who were treated with either azacytidine or decitabine (24 pre-treatment specimens, 23 post-treatment specimens) with 22 additional non-MDS control specimens. Total RNA was extracted from these specimens and submitted for next generation sequencing after an additional size exclusion step to enrich for small RNAs. The species of small RNAs were enumerated, single nucleotide variants (SNVs) identified, and finally the differential expression of tRNA-derived species (tDRs) in the specimens correlated with diseasestatus and response to therapy.

Results

Using miRNA sequencing data generated from bone marrow aspirate samples of patients with known MDS (N = 47) and controls (N = 23), we demonstrated that transfer RNA (tRNA) fragments (specifically tRNA halves, tRHs) are one of the most common species of small RNA isolated from size selection. Using tRNA expression values extracted from miRNA sequencing data, we identified six tRNA fragments that are differentially expressed between MDS and normal samples. Using the elastic net method, we identified four tRNAs-derived small RNAs (tDRs) that together can explain 67 % of the variation in treatment response for MDS patients. Similar analysis of specifically mitochondrial tDRs (mt-tDRs) identified 13 mt-tDRs which distinguished disease status in the samples and a single mt-tDR which predited response. Finally, 14 SNVs within the tDRs were found in at least 20 % of the MDS samples and were not observed in any of the control specimens.

Discussion

This study highlights the prevalence of tDRs in RNA-seq studies focused on small RNAs. The potential etiologies of these species, both technical and biologic, are discussed as well as important challenges in the interpretation of tDR data.

Conclusions

Our analysis results suggest that tRNA fragments can be accurately detected through miRNA sequencing data and that the expression of these species may be useful in the diagnosis of MDS and the prediction of response to therapy.

Electronic supplementary material

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

Does every transcript originate from a gene?

Outdated gene definitions favored regions corresponding to mature messenger RNAs, in particular, the open reading frame. In eukaryotes, the intergenic space was widely regarded nonfunctional and devoid of RNA transcription. Original concepts were based on the assumption that RNA expression was restricted to known protein-coding genes and a few so-called structural RNA genes, such as ribosomal RNAs or transfer RNAs. With the discovery of introns and, more recently, sensitive techniques for monitoring genome-wide transcription, this view had to be substantially modified. Tiling microarrays and RNA deep sequencing revealed myriads of transcripts, which cover almost entire genomes. The tremendous complexity of non-protein-coding RNA transcription has to be integrated into novel gene definitions. Despite an ever-growing list of functional RNAs, questions concerning the mass of identified transcripts are under dispute. Here, we examined genome-wide transcription from various angles, including evolutionary considerations, and suggest, in analogy to novel alternative splice variants that do not persist, that the vast majority of transcripts represent raw material for potential, albeit rare, exaptation events.