Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2015 Sep 24:16:727.
doi: 10.1186/s12864-015-1929-y.

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

Affiliations

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

Yan Guo et al. BMC Genomics. .

Abstract

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.

PubMed Disclaimer

Figures

Fig. 1
Fig. 1
Read count and alignment distribution example taken from one sample. The figures were produced using all read counts per category, not just unique reads per category. The other samples in this study follow a similar pattern. a. Read count distribution after trimming adaptors. The smaller peak at 22 base pairs indicates the abundance of miRNA and the larger peak at 33 base pairs indicates the abundance primarily of tRNA. b. The reads alignment distribution by RNA type. The majority of the reads aligned to tRNA instead of miRNA
Fig. 2
Fig. 2
a. Cluster analysis and heatmap using tRNA expression of all samples. Three phenotype bars are drawn below the dendrogram: pre-treatment, post-treatment and normal controls. Two clusters are visible (light green and light red). These two clusters do not separate pre- and post-treatment, but distinguish MDS and normal samples reasonably well. b. The six differentially expressed tRNA between disease and normal

Similar articles

Cited by

References

    1. Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63. doi: 10.1038/nrg2484. - DOI - PMC - PubMed
    1. Shendure J. The beginning of the end for microarrays? Nat Methods. 2008;5(7):585–7. doi: 10.1038/nmeth0708-585. - DOI - PubMed
    1. Guo Y, Sheng Q, Li J, Ye F, Samuels DC, Shyr Y. Large Scale Comparison of Gene Expression Levels by Microarrays and RNAseq Using TCGA Data. PLoS One. 2013;8(8):e71462. doi: 10.1371/journal.pone.0071462. - DOI - PMC - PubMed
    1. Guo Y, Li CI, Ye F, Shyr Y. Evaluation of read count based RNAseq analysis methods. BMC Genomics. 2013;14(Suppl 8):S2. doi: 10.1186/1471-2164-14-S8-S2. - DOI - PMC - PubMed
    1. Asmann YW, Klee EW, Thompson EA, Perez EA, Middha S, Oberg AL, et al. 3′ tag digital gene expression profiling of human brain and universal reference RNA using Illumina Genome Analyzer. BMC Genomics. 2009;10:531. doi: 10.1186/1471-2164-10-531. - DOI - PMC - PubMed

Publication types