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

Abstract

BackgroundAlthough 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).MethodsPaired 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.ResultsUsing 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.DiscussionThis 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.ConclusionsOur 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 materialThe online version of this article (doi:10.1186/s12864-015-1929-y) contains supplementary material, which is available to authori...

Genome-wide profiling of transfer RNAs and their role as novel prognostic markers for breast cancer

“…Therefore the reads captured in this study may likely be the fragments of tRNAs but it is not certain if these fragments are representative of mature full length tRNAs or if they represent actual physiological products (identified as tRFs). Nevertheless, we and others45 have confirmed that tRNA sequences can be accurately captured using small RNA sequencing and the reads captured from our study represent the known tRNAs identified and annotated to-date across all chromosomes (Supplementary Table S4). …”

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

“…Although miRNAs are only one of the many sncRNA species in smRNA-seq data sets, miRNAs remain the most popular class to study, largely because their biogenesis is relatively well understood and because the regulatory mechanism in post-transcription is known (42). However, more and more evidence shows that other non-miRNA sncRNAs also play important roles in gene regulation and certain diseases, such as cancers (5,7,43–45). For instance, there is an increasing amount of knowledge regarding the role of snoRNAs in cancer progression, and the information obtained thus far suggests that snoRNAs hold considerable potential for use as novel biomarkers and therapeutic targets in cancer treatment (4–6).…”

“…For example, a growing number of studies have reported that aberrant piRNA expression is a signature marker across distinct tumor types (3) and that snoRNAs act as oncogenes in tumorigenesis (4–6). The integration of large-scale smRNA-seq data helps researchers study the roles of these functional sncRNAs in cancer progression, but questions remain concerning the optimal methodologies for analysis, translation and utilization of such massive amounts of data (7). …”

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

“…35,36 In short, library construction from total RNA used the TruSeq Small RNA sample preparation kit (Illumina, San Diego, CA) and sequencing was performed on the Illumina HiSeq2500. Raw sequencing reads in BCL format were processed through CASAVA-1.8.2 for FASTQ conversion and demultiplexing.…”