ALL-tRNAseq enables robust tRNA profiling in tissue samples.

Transfer RNAs (tRNAs) are small adaptor RNAs essential for mRNA translation. Alterations in the cellular tRNA population can directly affect mRNA decoding rates and translational efficiency during cancer development and progression. To evaluate changes in the composition of the tRNA pool, multiple sequencing approaches have been developed to overcome reverse transcription blocks caused by the stable structures of these molecules and their numerous base modifications. However, it remains unclear whether current sequencing protocols faithfully capture tRNAs existing in cells or tissues. This is specifically challenging for clinical tissue samples that often present variable RNA qualities. For this reason, we developed ALL-tRNAseq, which combines the highly processive MarathonRT and RNA demethylation for the robust assessment of tRNA expression, together with a randomized adapter ligation strategy prior to reverse transcription to assess tRNA fragmentation levels in both cell lines and tissues. Incorporation of tRNA fragments not only informed on sample integrity but also significantly improved tRNA profiling of tissue samples. Our data showed that our profiling strategy effectively improves classification of oncogenic signatures in glioblastoma and diffuse large B-cell lymphoma tissues, particularly for samples presenting higher levels of RNA fragmentation, further highlighting the utility of ALL-tRNAseq for translational research.

miEAA 2023: updates, new functional microRNA sets and improved enrichment visualizations.

MicroRNAs (miRNAs) are small non-coding RNAs that play a critical role in regulating diverse biological processes. Extracting functional insights from a list of miRNAs is challenging, as each miRNA can potentially interact with hundreds of genes. To address this challenge, we developed miEAA, a flexible and comprehensive miRNA enrichment analysis tool based on direct and indirect miRNA annotation. The latest release of miEAA includes a data warehouse of 19 miRNA repositories, covering 10 different organisms and 139 399 functional categories. We have added information on the cellular context of miRNAs, isomiRs, and high-confidence miRNAs to improve the accuracy of the results. We have also improved the representation of aggregated results, including interactive Upset plots to aid users in understanding the interaction among enriched terms or categories. Finally, we demonstrate the functionality of miEAA in the context of ageing and highlight the importance of carefully considering the miRNA input list. MiEAA is free to use and publicly available at https://www.ccb.uni-saarland.de/mieaa/.

isomiRdb: microRNA expression at isoform resolution.

A significant fraction of mature miRNA transcripts carries sequence and/or length variations, termed isomiRs. IsomiRs are differentially abundant in cell types, tissues, body fluids or patients' samples. Not surprisingly, multiple studies describe a physiological and pathophysiological role. Despite their importance, systematically collected and annotated isomiR information available in databases remains limited. We thus developed isomiRdb, a comprehensive resource that compiles miRNA expression data at isomiR resolution from various sources. We processed 42 499 human miRNA-seq datasets (5.9 × 1011 sequencing reads) and consistently analyzed them using miRMaster and sRNAbench. Our database provides online access to the 90 483 most abundant isomiRs (>1 RPM in at least 1% of the samples) from 52 tissues and 188 cell types. Additionally, the full set of over 3 million detected isomiRs is available for download. Our resource can be queried at the sample, miRNA or isomiR level so users can quickly answer common questions about the presence/absence of a particular miRNA/isomiR in tissues of interest. Further, the database facilitates to identify whether a potentially interesting new isoform has been detected before and its frequency. In addition to expression tables, isomiRdb can generate multiple interactive visualisations including violin plots and heatmaps. isomiRdb is free to use and publicly available at: https://www.ccb.uni-saarland.de/isomirdb.

The RNA binding protein IGF2BP2/IMP2 alters the cargo of cancer cell-derived extracellular vesicles supporting tumor-associated macrophages.

BACKGROUND: Tumor cells release extracellular vesicles (EVs) that contribute to the polarization of macrophages towards tumor-associated macrophages (TAMs). High expression levels of the RNA binding protein IGF2BP2/IMP2 are correlated with increased tumor cell proliferation, invasion, and poor prognosis in the clinic. However, there is a lack of understanding of whether IMP2 affects the cargo of cancer cell-derived EVs, thereby modulating macrophage polarization. METHODS: EVs were isolated from IMP2-expressing HCT116 parental cells (WT) and CRISPR/Cas9 IMP2 knockout (KO) cells. EVs were characterized according to MISEV guidelines, microRNA cargo was assessed by microRNA-Seq, and the protein cargo was analyzed by proteomics. Primary human monocyte-derived macrophages (HMDMs) were polarized by EVs, and the expression of genes and surface markers was assessed using qPCR and flow cytometry, respectively. Morphological changes of macrophages, as well as the migratory potential of cancer cells, were assessed by the Incucyte® system and macrophage matrix degradation potential by zymography. Changes in the metabolic activity of macrophages were quantified using a Seahorse® analyzer. For in vivo studies, EVs were injected into the yolk sac of zebrafish larvae, and macrophages were isolated by fluorescence-activated cell sorting. RESULTS: EVs from WT and KO cells had a similar size and concentration and were positive for 25 vesicle markers. The expression of tumor-promoting genes was higher in macrophages polarized with WT EVs than KO EVs, while the expression of TNF and IL6 was reduced. A similar pattern was observed in macrophages from zebrafish larvae treated in vivo. WT EV-polarized macrophages showed a higher abundance of TAM-like surface markers, higher matrix degrading activity, as well as a higher promotion of cancer cell migration. MicroRNA-Seq revealed a significant difference in the microRNA composition of WT and KO EVs, particularly a high abundance of miR-181a-5p in WT EVs, which was absent in KO EVs. Inhibitors of macropinocytosis and phagocytosis antagonized the delivery of miR-181a-5p into macrophages and the downregulation of the miR-181a-5p target DUSP6. Proteomics data showed differences in protein cargo in KO vs. WT EVs, with the differentially abundant proteins mainly involved in metabolic pathways. WT EV-treated macrophages exhibited a higher basal oxygen consumption rate and a lower extracellular acidification rate than KO EV-treated cells. CONCLUSION: Our results show that IMP2 determines the cargo of EVs released by cancer cells, thereby modulating the EVs' actions on macrophages. Expression of IMP2 is linked to the secretion of EVs that polarize macrophages towards a tumor-promoting phenotype.

Skin treatment with non-thermal plasma modulates the immune system through miR-223-3p and its target genes.

Non-thermal plasma, a partially ionized gas, holds significant potential for clinical applications, including wound-healing support, oral therapies, and anti-tumour treatments. While its applications showed promising outcomes, the underlying molecular mechanisms remain incompletely understood. We thus apply non-thermal plasma to mouse auricular skin and conducted non-coding RNA sequencing, as well as single-cell blood sequencing. In a time-series analysis (five timepoints spanning 2 hours), we compare the expression of microRNAs in the plasma-treated left ears to the unexposed right ears of the same mice as well as to the ears of unexposed control mice. Our findings indicate specific effects in the treated ears for a set of five miRNAs: mmu-miR-144-5p, mmu-miR-144-3p, mmu-miR-142a-5p, mmu-miR-223-3p, and mmu-miR-451a. Interestingly, mmu-miR-223-3p also exhibits an increase over time in the right non-treated ear of the exposed mice, suggesting systemic effects. Notably, this miRNA, along with mmu-miR-142a-5p and mmu-miR-144-3p, regulates genes and pathways associated with wound healing and tissue regeneration (namely ErbB, FoxO, Hippo, and PI3K-Akt signalling). This co-regulation is particularly remarkable considering the significant seed dissimilarities among the miRNAs. Finally, single-cell sequencing of PBMCs reveals the downregulation of 12 from 15 target genes in B-cells, Cd4+ and Cd8+ T-cells. Collectively, our data provide evidence for a systemic effect of non-thermal plasma.

sRNAbench and sRNAtoolbox 2022 update: accurate miRNA and sncRNA profiling for model and non-model organisms.

The NCBI Sequence Read Archive currently hosts microRNA sequencing data for over 800 different species, evidencing the existence of a broad taxonomic distribution in the field of small RNA research. Simultaneously, the number of samples per miRNA-seq study continues to increase resulting in a vast amount of data that requires accurate, fast and user-friendly analysis methods. Since the previous release of sRNAtoolbox in 2019, 55 000 sRNAbench jobs have been submitted which has motivated many improvements in its usability and the scope of the underlying annotation database. With this update, users can upload an unlimited number of samples or import them from Google Drive, Dropbox or URLs. Micro- and small RNA profiling can now be carried out using high-confidence Metazoan and plant specific databases, MirGeneDB and PmiREN respectively, together with genome assemblies and libraries from 441 Ensembl species. The new results page includes straightforward sample annotation to allow downstream differential expression analysis with sRNAde. Unassigned reads can also be explored by means of a new tool that performs mapping to microbial references, which can reveal contamination events or biologically meaningful findings as we describe in the example. sRNAtoolbox is available at: https://arn.ugr.es/srnatoolbox/.

MirGeneDB 2.1: toward a complete sampling of all major animal phyla.

We describe an update of MirGeneDB, the manually curated microRNA gene database. Adhering to uniform and consistent criteria for microRNA annotation and nomenclature, we substantially expanded MirGeneDB with 30 additional species representing previously missing metazoan phyla such as sponges, jellyfish, rotifers and flatworms. MirGeneDB 2.1 now consists of 75 species spanning over ∼800 million years of animal evolution, and contains a total number of 16 670 microRNAs from 1549 families. Over 6000 microRNAs were added in this update using ∼550 datasets with ∼7.5 billion sequencing reads. By adding new phylogenetically important species, especially those relevant for the study of whole genome duplication events, and through updating evolutionary nodes of origin for many families and genes, we were able to substantially refine our nomenclature system. All changes are traceable in the specifically developed MirGeneDB version tracker. The performance of read-pages is improved and microRNA expression matrices for all tissues and species are now also downloadable. Altogether, this update represents a significant step toward a complete sampling of all major metazoan phyla, and a widely needed foundation for comparative microRNA genomics and transcriptomics studies. MirGeneDB 2.1 is part of RNAcentral and Elixir Norway, publicly and freely available at http://www.mirgenedb.org/.

Ageing-associated small RNA cargo of extracellular vesicles.

Previous work on murine models and humans demonstrated global as well as tissue-specific molecular ageing trajectories of RNAs. Extracellular vesicles (EVs) are membrane vesicles mediating the horizontal transfer of genetic information between different tissues. We sequenced small regulatory RNAs (sncRNAs) in two mouse plasma fractions at five time points across the lifespan from 2-18 months: (1) sncRNAs that are free-circulating (fc-RNA) and (2) sncRNAs bound outside or inside EVs (EV-RNA). Different sncRNA classes exhibit unique ageing patterns that vary between the fcRNA and EV-RNA fractions. While tRNAs showed the highest correlation with ageing in both fractions, rRNAs exhibited inverse correlation trajectories between the EV- and fc-fractions. For miRNAs, the EV-RNA fraction was exceptionally strongly associated with ageing, especially the miR-29 family in adipose tissues. Sequencing of sncRNAs and coding genes in fat tissue of an independent cohort of aged mice up to 27 months highlighted the pivotal role of miR-29a-3p and miR-29b-3p in ageing-related gene regulation that we validated in a third cohort by RT-qPCR.

miRTargetLink 2.0-interactive miRNA target gene and target pathway networks.

Which genes, gene sets or pathways are regulated by certain miRNAs? Which miRNAs regulate a particular target gene or target pathway in a certain physiological context? Answering such common research questions can be time consuming and labor intensive. Especially for researchers without computational experience, the integration of different data sources, selection of the right parameters and concise visualization can be demanding. A comprehensive analysis should be central to present adequate answers to complex biological questions. With miRTargetLink 2.0, we develop an all-in-one solution for human, mouse and rat miRNA networks. Users input in the unidirectional search mode either a single gene, gene set or gene pathway, alternatively a single miRNA, a set of miRNAs or an miRNA pathway. Moreover, genes and miRNAs can jointly be provided to the tool in the bidirectional search mode. For the selected entities, interaction graphs are generated from different data sources and dynamically presented. Connected application programming interfaces (APIs) to the tailored enrichment tools miEAA and GeneTrail facilitate downstream analysis of pathways and context-annotated categories of network nodes. MiRTargetLink 2.0 is freely accessible at https://www.ccb.uni-saarland.de/mirtargetlink2.

Validation of human microRNA target pathways enables evaluation of target prediction tools.

MicroRNAs are regulators of gene expression. A wide-spread, yet not validated, assumption is that the targetome of miRNAs is non-randomly distributed across the transcriptome and that targets share functional pathways. We developed a computational and experimental strategy termed high-throughput miRNA interaction reporter assay (HiTmIR) to facilitate the validation of target pathways. First, targets and target pathways are predicted and prioritized by computational means to increase the specificity and positive predictive value. Second, the novel webtool miRTaH facilitates guided designs of reporter assay constructs at scale. Third, automated and standardized reporter assays are performed. We evaluated HiTmIR using miR-34a-5p, for which TNF- and TGFB-signaling, and Parkinson's Disease (PD)-related categories were identified and repeated the pipeline for miR-7-5p. HiTmIR validated 58.9% of the target genes for miR-34a-5p and 46.7% for miR-7-5p. We confirmed the targeting by measuring the endogenous protein levels of targets in a neuronal cell model. The standardized positive and negative targets are collected in the new miRATBase database, representing a resource for training, or benchmarking new target predictors. Applied to 88 target predictors with different confidence scores, TargetScan 7.2 and miRanda outperformed other tools. Our experiments demonstrate the efficiency of HiTmIR and provide evidence for an orchestrated miRNA-gene targeting.

Comprehensive, integrative genomic analysis of microRNA expression profiles in different tissues of two wheat cultivars with different traits.

Wheat is one of the most important food sources on Earth. MicroRNAs (miRNA) play important roles in wheat productivity. To identify wheat miRNAs, we constructed and sequenced sRNA libraries from leaves and roots of two wheat cultivars (RAC875 and Kukri) with many different traits. Given that available miRNA wheat complement in the plant-specific database PmiREN ( https://pmiren.com ) does not include root tissues and root-associated miRNAs might thus be missing, we performed first the prediction of novel miRNAs using the sRNAbench tool. We found a total of 150 putatively novel miRNA genes with expression of both arms from 289 unique mature sequences and nearly 30% of all miRNA reads in roots corresponded to novel miRNAs. In contrast, this figure in leaves dropped to under 3%, confirming the undersampling of roots in the complement of known miRNAs. By using 120 publicly available wheat datasets, 598 Zea mays small RNA libraries, 64 plant species genomes, wheat degradome library, and functional enrichment analysis, a subset of novel miRNAs were confirmed as bona-fide miRNAs. Of the total 605 miRNAs identified in this study inclusive of 316 known miRNAs, 528 miRNAs were shared by both cultivars, 429 miRNAs were shared by both root tissues and 329 miRNAs were shared by both leaf tissues. In addition, 32 miRNAs were specific to Kukri while 45 miRNAs were specific to RAC875. These miRNAs had diverse functions, such as regulation of gene transcription, protein translation, energy metabolism, and cell cycle progression. Our data provide a genome-wide miRNA expression profile in these two wheat cultivars and help functional studies of wheat genomics.

Quantitative proteomics analysis reveals core and variable tick salivary proteins at the tick-vertebrate host interface.

Few studies have examined tick proteomes, how they adapt to their environment, and their roles in the parasite-host interactions that drive tick infestation and pathogen transmission. Here we used a proteomics approach to screen for biologically and immunologically relevant proteins acting at the tick-host interface during tick feeding and, as proof of principle, measured host antibody responses to some of the discovered candidates. We used a label-free quantitative proteomic workflow to study salivary proteomes of (i) wild Ixodes ricinus ticks fed on different hosts, (ii) wild or laboratory ticks fed on the same host, and (iii) adult ticks cofed with nymphs. Our results reveal high and stable expression of several protease inhibitors and other tick-specific proteins under different feeding conditions. Most pathways functionally enriched in sialoproteomes were related to proteolysis, endopeptidase, and amine-binding activities. The generated catalogue of tick salivary proteins enabled the selection of six candidate secreted immunogenic peptides for rabbit immunizations, three of which induced strong and durable antigen-specific antibody responses in rabbits. Furthermore, rabbits exposed to ticks mounted immune responses against the candidate peptides/proteins, confirming their expression at the tick-vertebrate interface. Our approach provides insights into tick adaptation strategies to different feeding conditions and promising candidates for developing antitick vaccines or markers of exposure of vertebrate hosts to tick bites.

mirnaQC: a webserver for comparative quality control of miRNA-seq data.

Although miRNA-seq is extensively used in many different fields, its quality control is frequently restricted to a PhredScore-based filter. Other important quality related aspects like microRNA yield, the fraction of putative degradation products (such as rRNA fragments) or the percentage of adapter-dimers are hard to assess using absolute thresholds. Here we present mirnaQC, a webserver that relies on 34 quality parameters to assist in miRNA-seq quality control. To improve their interpretability, quality attributes are ranked using a reference distribution obtained from over 36 000 publicly available miRNA-seq datasets. Accepted input formats include FASTQ and SRA accessions. The results page contains several sections that deal with putative technical artefacts related to library preparation, sequencing, contamination or yield. Different visualisations, including PCA and heatmaps, are available to help users identify underlying issues. Finally, we show the usefulness of this approach by analysing two publicly available datasets and discussing the different quality issues that can be detected using mirnaQC.

MirGeneDB 2.0: the metazoan microRNA complement.

Small non-coding RNAs have gained substantial attention due to their roles in animal development and human disorders. Among them, microRNAs are special because individual gene sequences are conserved across the animal kingdom. In addition, unique and mechanistically well understood features can clearly distinguish bona fide miRNAs from the myriad other small RNAs generated by cells. However, making this distinction is not a common practice and, thus, not surprisingly, the heterogeneous quality of available miRNA complements has become a major concern in microRNA research. We addressed this by extensively expanding our curated microRNA gene database - MirGeneDB - to 45 organisms, encompassing a wide phylogenetic swath of animal evolution. By consistently annotating and naming 10,899 microRNA genes in these organisms, we show that previous microRNA annotations contained not only many false positives, but surprisingly lacked >2000 bona fide microRNAs. Indeed, curated microRNA complements of closely related organisms are very similar and can be used to reconstruct ancestral miRNA repertoires. MirGeneDB represents a robust platform for microRNA-based research, providing deeper and more significant insights into the biology and evolution of miRNAs as well as biomedical and biomarker research. MirGeneDB is publicly and freely available at http://mirgenedb.org/.

Unification of miRNA and isomiR research: the mirGFF3 format and the mirtop API.

MOTIVATION: MicroRNAs (miRNAs) are small RNA molecules (∼22 nucleotide long) involved in post-transcriptional gene regulation. Advances in high-throughput sequencing technologies led to the discovery of isomiRs, which are miRNA sequence variants. While many miRNA-seq analysis tools exist, the diversity of output formats hinders accurate comparisons between tools and precludes data sharing and the development of common downstream analysis methods. RESULTS: To overcome this situation, we present here a community-based project, miRNA Transcriptomic Open Project (miRTOP) working towards the optimization of miRNA analyses. The aim of miRTOP is to promote the development of downstream isomiR analysis tools that are compatible with existing detection and quantification tools. Based on the existing GFF3 format, we first created a new standard format, mirGFF3, for the output of miRNA/isomiR detection and quantification results from small RNA-seq data. Additionally, we developed a command line Python tool, mirtop, to create and manage the mirGFF3 format. Currently, mirtop can convert into mirGFF3 the outputs of commonly used pipelines, such as seqbuster, isomiR-SEA, sRNAbench, Prost! as well as BAM files. Some tools have also incorporated the mirGFF3 format directly into their code, such as, miRge2.0, IsoMIRmap and OptimiR. Its open architecture enables any tool or pipeline to output or convert results into mirGFF3. Collectively, this isomiR categorization system, along with the accompanying mirGFF3 and mirtop API, provide a comprehensive solution for the standardization of miRNA and isomiR annotation, enabling data sharing, reporting, comparative analyses and benchmarking, while promoting the development of common miRNA methods focusing on downstream steps of miRNA detection, annotation and quantification. AVAILABILITY AND IMPLEMENTATION: https://github.com/miRTop/mirGFF3/ and https://github.com/miRTop/mirtop. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

liqDB: a small-RNAseq knowledge discovery database for liquid biopsy studies.

MiRNAs are important regulators of gene expression and are frequently deregulated under pathologic conditions. They are highly stable in bodily fluids which makes them feasible candidates to become minimally invasive biomarkers. In fact, several studies already proposed circulating miRNA-based biomarkers for different types of neoplastic, cardiovascular and degenerative diseases. However, many of these studies rely on small RNA sequencing experiments that are based on different RNA extraction and processing protocols, rendering results incomparable. We generated liqDB, a database for liquid biopsy small RNA sequencing profiles that provides users with meaningful information to guide their small RNA liquid biopsy research and to overcome technical and conceptual problems. By means of a user-friendly web interface, miRNA expression profiles from 1607 manually annotated samples can be queried and explored at different levels. Result pages include downloadable expression matrices, differential expression analysis, most stably expressed miRNAs, cluster analysis and relevant visualizations by means of boxplots and heatmaps. We anticipate that liqDB will be a useful tool in liquid biopsy research as it provides a consistently annotated large compilation of experiments together with tools for reproducible analysis, comparison and hypothesis generation. LiqDB is available at http://bioinfo5.ugr.es/liqdb.

sRNAbench and sRNAtoolbox 2019: intuitive fast small RNA profiling and differential expression.

Since the original publication of sRNAtoolbox in 2015, small RNA research experienced notable advances in different directions. New protocols for small RNA sequencing have become available to address important issues such as adapter ligation bias, PCR amplification artefacts or to include internal controls such as spike-in sequences. New microRNA reference databases were developed with different foci, either prioritizing accuracy (low number of false positives) or completeness (low number of false negatives). Additionally, other small RNA molecules as well as microRNA sequence and length variants (isomiRs) have continued to gain importance. Finally, the number of microRNA sequencing studies deposited in GEO nearly triplicated from 2014 (280) to 2018 (764). These developments imply that fast and easy-to-use tools for expression profiling and subsequent downstream analysis of miRNA-seq data are essential to many researchers. Key features in this sRNAtoolbox release include addition of all major RNA library preparation protocols to sRNAbench and improvements in sRNAde, a tool that summarizes several aspects of small RNA sequencing studies including the detection of consensus differential expression. A special emphasis was put on the user-friendliness of the tools, for instance sRNAbench now supports parallel launching of several jobs to improve reproducibility and user time efficiency.

miRglmm: a generalized linear mixed model of isomiR-level counts improves estimation of miRNA-level differential expression and uncovers variable differential expression between isomiRs.

MicroRNA-seq data is produced by aligning small RNA sequencing reads of different miRNA transcript isoforms, called isomiRs, to known microRNAs. Aggregation to microRNA-level counts discards information and violates core assumptions of differential expression (DE) methods developed for mRNA-seq data. We establish miRglmm, a DE method for microRNA-seq data, that uses a generalized linear mixed model of isomiR-level counts, facilitating detection of miRNA with differential expression or differential isomiR usage. We demonstrate that miRglmm outperforms current DE methods in estimating DE for miRNA, whether or not there is significant isomiR variability, and simultaneously provides estimates of isomiR-level DE.

Full-length tRNAs lacking a functional CCA tail are selectively sorted into the lumen of extracellular vesicles.

Small extracellular vesicles (sEVs) are heterogenous lipid membrane particles typically less than 200 nm in size and secreted by most cell types either constitutively or upon activation signals. sEVs isolated from biofluids contain RNAs, including small non-coding RNAs (ncRNAs), that can be either encapsulated within the EV lumen or bound to the EV surface. EV-associated microRNAs (miRNAs) are, despite a relatively low abundance, extensively investigated for their selective incorporation and their role in cell-cell communication. In contrast, the sorting of highly-structured ncRNA species is understudied, mainly due to technical limitations of traditional small RNA sequencing protocols. Here, we adapted ALL-tRNAseq to profile the relative abundance of highly structured and potentially methylated small ncRNA species, including transfer RNAs (tRNAs), small nucleolar RNAs (snoRNAs), and Y RNAs in bulk EV preparations. We determined that full-length tRNAs, typically 75 to 90 nucleotides in length, were the dominant small ncRNA species (>60% of all reads in the 18-120 nucleotides size-range) in all cell culture-derived EVs, as well as in human plasma-derived EV samples, vastly outnumbering 21 nucleotides-long miRNAs. Nearly all EV-associated tRNAs were protected from external RNAse treatment, indicating a location within the EV lumen. Strikingly, the vast majority of luminal-sorted, full-length, nucleobase modification-containing EV-tRNA sequences, harbored a dysfunctional 3' CCA tail, 1 to 3 nucleotides truncated, rendering them incompetent for amino acid loading. In contrast, in non-EV associated extracellular particle fractions (NVEPs), tRNAs appeared almost exclusively fragmented or 'nicked' into tRNA-derived small RNAs (tsRNAs) with lengths between 18 to 35 nucleotides. We propose that in mammalian cells, tRNAs that lack a functional 3' CCA tail are selectively sorted into EVs and shuttled out of the producing cell, offering a new perspective into the physiological role of secreted EVs and luminal cargo-selection.