Graphical AbstractHighlights d exRNA sequencing complexity and reproducibility varies across isolation methods d Deconvolution shows differential access to exRNA carriers by different methods d Performance of exRNA isolation methods vary across biofluids and RNA species d miRDaR enables customized selection of optimal exRNA isolation methods SUMMARYPoor reproducibility within and across studies arising from lack of knowledge regarding the performance of extracellular RNA (exRNA) isolation methods has hindered progress in the exRNA field. A systematic comparison of 10 exRNA isolation methods across 5 biofluids revealed marked differences in the complexity and reproducibility of the resulting small RNA-seq profiles. The relative efficiency with which each method accessed different exRNA carrier subclasses was determined by estimating the proportions of extracellular vesicle (EV)-, ribonucleoprotein (RNP)-, and highdensity lipoprotein (HDL)-specific miRNA signatures in each profile. An interactive web-based application (miRDaR) was developed to help investigators select the optimal exRNA isolation method for their studies. miRDar provides comparative statistics for all expressed miRNAs or a selected subset of miRNAs in the desired biofluid for each exRNA isolation method and returns a ranked list of exRNA isolation methods prioritized by complexity, expression level, and repro-ducibility. These results will improve reproducibility and stimulate further progress in exRNA biomarker development.

Isoforms of human miRNAs (isomiRs) are constitutively expressed with tissue- and disease-subtype-dependencies. We studied 10 271 tumor datasets from The Cancer Genome Atlas (TCGA) to evaluate whether isomiRs can distinguish amongst 32 TCGA cancers. Unlike previous approaches, we built a classifier that relied solely on ‘binarized’ isomiR profiles: each isomiR is simply labeled as ‘present’ or ‘absent’. The resulting classifier successfully labeled tumor datasets with an average sensitivity of 90% and a false discovery rate (FDR) of 3%, surpassing the performance of expression-based classification. The classifier maintained its power even after a 15× reduction in the number of isomiRs that were used for training. Notably, the classifier could correctly predict the cancer type in non-TCGA datasets from diverse platforms. Our analysis revealed that the most discriminatory isomiRs happen to also be differentially expressed between normal tissue and cancer. Even so, we find that these highly discriminating isomiRs have not been attracting the most research attention in the literature. Given their ability to successfully classify datasets from 32 cancers, isomiRs and our resulting ‘Pan-cancer Atlas’ of isomiR expression could serve as a suitable framework to explore novel cancer biomarkers.

MINTbase is a repository that comprises nuclear and mitochondrial tRNA-derived fragments (‘tRFs’) found in multiple human tissues. The original version of MINTbase comprised tRFs obtained from 768 transcriptomic datasets. We used our deterministic and exhaustive tRF mining pipeline to process all of The Cancer Genome Atlas datasets (TCGA). We identified 23 413 tRFs with abundance of ≥ 1.0 reads-per-million (RPM). To facilitate further studies of tRFs by the community, we just released version 2.0 of MINTbase that contains information about 26 531 distinct human tRFs from 11 719 human datasets as of October 2017. Key new elements include: the ability to filter tRFs on-the-fly by minimum abundance thresholding; the ability to filter tRFs by tissue keywords; easy access to information about a tRF’s maximum abundance and the datasets that contain it; the ability to generate relative abundance plots for tRFs across cancer types and convert them into embeddable figures; MODOMICS information about modifications of the parental tRNA, etc. Version 2.0 of MINTbase contains 15x more datasets and nearly 4x more distinct tRFs than the original version, yet continues to offer fast, interactive access to its contents. Version 2.0 is available freely at http://cm.jefferson.edu/MINTbase/.

The fragments that derive from transfer RNAs (tRNAs) are an emerging category of regulatory RNAs. Known as tRFs, these fragments were reported for the first time only a decade ago, making them a relatively recent addition to the ever-expanding pantheon of non-coding RNAs. tRFs are short, 16–35 nucleotides (nts) in length, and produced through cleavage of mature and precursor tRNAs at various positions. Both cleavage positions and relative tRF abundance depend strongly on context, including the tissue type, tissue state, and disease, as well as the sex, population of origin, and race/ethnicity of an individual. These dependencies increase the urgency to understand the regulatory roles of tRFs. Such efforts are gaining momentum, and comprise experimental and computational approaches. System-level studies across many tissues and thousands of samples have produced strong evidence that tRFs have important and multi-faceted roles. Here, we review the relevant literature on tRF biology in higher organisms, single cell eukaryotes, and prokaryotes.

tRNA-derived fragments (tRF) are a class of potent regulatory RNAs. We mined the datasets from The Cancer Genome Atlas (TCGA) representing 32 cancer types with a deterministic and exhaustive pipeline for tRNA fragments. We found that mitochondrial tRNAs contribute disproportionally more tRFs than nuclear tRNAs. Through integrative analyses, we uncovered a multitude of statistically significant and contextdependent associations between the identified tRFs and mRNAs. In many of the 32 cancer types, these associations involve mRNAs from developmental processes, receptor tyrosine kinase signaling, the proteasome, and metabolic pathways that include glycolysis, oxidative phosphorylation, and ATP synthesis. Even though the pathways are common to multiple cancers, the association of specific mRNAs with tRFs depends on and differs from cancer to cancer. The associations between tRFs and mRNAs extend to genomic properties as well; specifically, tRFs are positively correlated with shorter genes that have a higher density in repeats, such as ALUs, MIRs, and ERVLs. Conversely, tRFs are negatively correlated with longer genes that have a lower repeat density, suggesting a possible dichotomy between cell proliferation and differentiation. Analyses of bladder, lung, and kidney cancer data indicate that the tRF-mRNA wiring can also depend on a patient's sex. Sex-dependent associations involve cyclindependent kinases in bladder cancer, the MAPK signaling pathway in lung cancer, and purine metabolism in kidney cancer. Taken together, these findings suggest diverse and wide-ranging roles for tRFs and highlight the extensive interconnections of tRFs with key cellular processes and human genomic architecture.Significance: Across 32 TCGA cancer contexts, nuclear and mitochondrial tRNA fragments exhibit associations with mRNAs that belong to concrete pathways, encode proteins with particular destinations, have a biased repeat content, and are sex dependent.