Detection of circulating extracellular mRNAs by modified small-RNA-sequencing analysis

A dataset of hidden small non-coding RNA in the testis of heat-stressed models revealed by Pandora-seq

Infertility, a worldwide reproductive health concern, impacts approximately one in five couples. Male infertility, stemming from spermatogenic dysfunction and reduced sperm quality, stands as a primary factor contributing to infertility. Given the global decrease in male fertility linked to environmental factors like the greenhouse effect, it is crucial to develop a comprehensive understanding of how increased temperatures impact both the quantity and quality of sperm. In this study, we utilized Pandora-seq technology to detect the small non-coding RNAs (sncRNAs) expression profile in the testicular tissue of heat-stressed mice. The investigation explores the dynamic shifts in sncRNAs within the mouse testis under heat stress, including miRNAs, tsRNAs, piRNAs, rsRNAs and other sncRNAs. Furthermore, we successfully identified differentially expressed sncRNAs in testicular tissues before and after heat stress. Subsequently, we conducted functional enrichment analysis on the potential predicted target genes of differentially expressed miRNAs and tsRNAs. These datasets will constitute a valuable foundational resource for further investigations into the decline in male reproductive capacity triggered by heat stress.

Depletion-assisted multiplexed cell-free RNA sequencing reveals distinct human and microbial signatures in plasma versus extracellular vesicles

BACKGROUND

Cell-free long RNAs in human plasma and extracellular vesicles (EVs) have shown promise as biomarkers in liquid biopsy, despite their fragmented nature.

METHODS

To investigate these fragmented cell-free RNAs (cfRNAs), we developed a cost-effective cfRNA sequencing method called DETECTOR-seq (depletion-assisted multiplexed cell-free total RNA sequencing). DETECTOR-seq utilised a meticulously tailored set of customised guide RNAs to remove large amounts of unwanted RNAs (i.e., fragmented ribosomal and mitochondrial RNAs) in human plasma. Early barcoding strategy was implemented to reduce costs and minimise plasma requirements.

RESULTS

Using DETECTOR-seq, we conducted a comprehensive analysis of cell-free transcriptomes in both whole human plasma and EVs. Our analysis revealed discernible distributions of RNA types in plasma and EVs. Plasma exhibited pronounced enrichment in structured circular RNAs, tRNAs, Y RNAs and viral RNAs, while EVs showed enrichment in messenger RNAs (mRNAs) and signal recognition particle RNAs (srpRNAs). Functional pathway analysis highlighted RNA splicing-related ribonucleoproteins (RNPs) and antimicrobial humoral response genes in plasma, while EVs demonstrated enrichment in transcriptional activity, cell migration and antigen receptor-mediated immune signals. Our study indicates the comparable potential of cfRNAs from whole plasma and EVs in distinguishing cancer patients (i.e., colorectal and lung cancer) from healthy donors. And microbial cfRNAs in plasma showed potential in classifying specific cancer types.

CONCLUSIONS

Our comprehensive analysis of total and EV cfRNAs in paired plasma samples provides valuable insights for determining the need for EV purification in cfRNA-based studies. We envision the cost effectiveness and efficiency of DETECTOR-seq will empower transcriptome-wide investigations in the fields of cfRNAs and liquid biopsy.

KEYPOINTS

DETECTOR-seq (depletion-assisted multiplexed cell-free total RNA sequencing) enabled efficient and specific depletion of sequences derived from fragmented ribosomal and mitochondrial RNAs in plasma. Distinct human and microbial cell-free RNA (cfRNA) signatures in whole Plasma versus extracellular vesicles (EVs) were revealed. Both Plasma and EV cfRNAs were capable of distinguishing cancer patients from normal individuals, while microbial RNAs in Plasma cfRNAs enabled better classification of cancer types than EV cfRNAs.

The impact of preanalytical variables on the analysis of cell-free DNA from blood and urine samples

Cell-free DNA (cfDNA), a burgeoning class of molecular biomarkers, has been extensively studied across a variety of biomedical fields. As a key component of liquid biopsy, cfDNA testing is gaining prominence in disease detection and management due to the convenience of sample collection and the abundant wealth of genetic information it provides. However, the broader clinical application of cfDNA is currently impeded by a lack of standardization in the preanalytical procedures for cfDNA analysis. A number of fundamental challenges, including the selection of appropriate preanalytical procedures, prevention of short cfDNA fragment loss, and the validation of various cfDNA measurement methods, remain unaddressed. These existing hurdles lead to difficulties in comparing results and ensuring repeatability, thereby undermining the reliability of cfDNA analysis in clinical settings. This review discusses the crucial preanalytical factors that influence cfDNA analysis outcomes, including sample collection, transportation, temporary storage, processing, extraction, quality control, and long-term storage. The review provides clarification on achievable consensus and offers an analysis of the current issues with the goal of standardizing preanalytical procedures for cfDNA analysis.

The tRNAVal half: A strong endogenous Toll-like receptor 7 ligand with a 5'-terminal universal sequence signature

Toll-like receptors (TLRs) are crucial components of the innate immune system. Endosomal TLR7 recognizes single-stranded RNAs, yet its endogenous ssRNA ligands are not fully understood. We previously showed that extracellular (ex-) 5'-half molecules of tRNAHisGUG (the 5'-tRNAHisGUG half) in extracellular vesicles (EVs) of human macrophages activate TLR7 when delivered into endosomes of recipient macrophages. Here, we fully explored immunostimulatory ex-5'-tRNA half molecules and identified the 5'-tRNAValCAC/AAC half, the most abundant tRNA-derived RNA in macrophage EVs, as another 5'-tRNA half molecule with strong TLR7 activation capacity. Levels of the ex-5'-tRNAValCAC/AAC half were highly up-regulated in macrophage EVs upon exposure to lipopolysaccharide and in the plasma of patients infected with Mycobacterium tuberculosis. The 5'-tRNAValCAC/AAC half-mediated activation of TLR7 effectively eradicated bacteria infected in macrophages. Mutation analyses of the 5'-tRNAValCAC/AAC half identified the terminal GUUU sequence as a determinant for TLR7 activation. We confirmed that GUUU is the optimal ratio of guanosine and uridine for TLR7 activation; microRNAs or other RNAs with the terminal GUUU motif can indeed stimulate TLR7, establishing the motif as a universal signature for TLR7 activation. These results advance our understanding of endogenous ssRNA ligands of TLR7 and offer insights into diverse TLR7-involved pathologies and their therapeutic strategies.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Immunostimulatory short non-coding RNAs in the circulation of patients with tuberculosis infection

Mycobacterium tuberculosis (Mtb) infection is among the world's deadliest infectious diseases. Developing effective treatments and biomarkers for tuberculosis requires a deeper understanding of its pathobiology and host responses. Here, we report a comprehensive characterization of circulating short non-coding RNAs (sncRNAs) in plasma samples from Mtb-infected patients. We achieved this by pre-treating plasma RNAs with T4 polynucleotide kinase to convert all RNA ends to those compatible with sncRNA sequencing. We discovered a global and drastic upregulation of plasma sncRNAs in Mtb-infected patients, with tRNA-derived sncRNAs representing the most dramatically elevated class. Most of these tRNA-derived sncRNAs originated from a limited subset of tRNAs, specifically from three tRNA isoacceptors, and exhibited skewed patterns to 5'-derived fragments, such as 5' halves, 5' tRNA fragments (tRFs), and internal tRFs (i-tRFs) from the 5' regions. Further, Mtb-infected patients displayed markedly upregulated and distinct profiles of both rRNA- and mRNA-derived sncRNAs. Some of these sncRNAs, which are abundant and specific to Mtb-infected patients, robustly activated human macrophages via Toll-like receptor 7 and induced cytokine production. This drastic accumulation of circulating, immunostimulatory sncRNAs in the plasma of Mtb-infected patients offers insights into the sncRNA-driven aspects of host immune response against infectious diseases and suggests a pool of potential therapeutic targets and biomarkers.

tsRNA: A Promising Biomarker in Breast Cancer

tRNA-derived small RNAs (tsRNAs) are a novel class of non-coding small RNAs, generated from specific cleavage sites of tRNA or pre-tRNA. tsRNAs can directly participate in RNA silencing, transcription, translation, and other processes. Their dysregulation is closely related to the occurrence and development of various cancers. Breast cancer is one of the most common and fastest-growing malignant tumors in humans. tsRNAs have been found to be dysregulated in breast cancer, serving as a new target for exploring the pathogenesis of breast cancer. They are also considered new tumor markers, providing a basis for diagnosis and treatment. This article reviews the generation, classification, mechanism of action, function of tsRNAs, and their biological effects and related mechanisms in breast cancer, in the hope of providing a new direction for the diagnosis and treatment of breast cancer.

Small RNA structural biochemistry in a post-sequencing era

High-throughput sequencing has had an enormous impact on small RNA research during the past decade. However, sequencing only offers a one-dimensional view of the transcriptome and is often highly biased. Additionally, the 'sequence, map and annotate' approach, used widely in small RNA research, can lead to flawed interpretations of the data, lacking biological plausibility, due in part to database issues. Even in the absence of technical biases, the loss of three-dimensional information is a major limitation to understanding RNA stability, turnover and function. For example, noncoding RNA-derived fragments seem to exist mainly as dimers, tetramers or as nicked forms of their parental RNAs, contrary to widespread assumptions. In this perspective, we will discuss main sources of bias during small RNA-sequencing, present several useful bias-reducing strategies and provide guidance on the interpretation of small RNA-sequencing results, with emphasis on RNA fragmentomics. As sequencing offers a one-dimensional projection of a four-dimensional reality, prior structure-level knowledge is often needed to make sense of the data. Consequently, while less-biased sequencing methods are welcomed, integration of orthologous experimental techniques is also strongly recommended.

Challenges in characterization of transcriptomes of extracellular vesicles and non-vesicular extracellular RNA carriers

Since its original discovery over a decade ago, extracellular RNA (exRNA) has been found in all biological fluids. Furthermore, extracellular microRNA has been shown to be involved in communication between various cell types. Importantly, the exRNA is protected from RNases degradation by certain carriers including membrane vesicles and non-vesicular protein nanoparticles. Each type of carrier has its unique exRNA profile, which may vary depending on cell type and physiological conditions. To clarify putative mechanisms of intercellular communication mediated by exRNA, the RNA profile of each carrier has to be characterized. While current methods of biofluids fractionation are continuously improving, they fail to completely separate exRNA carriers. Likewise, most popular library preparation approaches for RNA sequencing do not allow obtaining exhaustive and unbiased data on exRNA transcriptome. In this mini review we discuss ongoing progress in the field of exRNA, with the focus on exRNA carriers, analyze the key methodological challenges and provide recommendations on how the latter could be overcome.

Discovery of the major 15-30 nt mammalian small RNAs, their biogenesis and function

Small RNAs (sRNAs) within 15-30 nt such as miRNA, tsRNA, srRNA with 3'-OH have been identified. However, whether these sRNAs are the major 15-30 nt sRNAs is still unknown. Here we show about 90% mammalian sRNAs within 15-30 nt end with 2',3'-cyclic phosphate (3'-cP). TANT-seq was developed to simultaneously profile sRNAs with 3'-cP (sRNA-cPs) and sRNA-OHs, and huge amount of sRNA-cPs were detected. Surprisingly, sRNA-cPs and sRNA-OHs usually have distinct sequences. The data from TANT-seq were validated by a novel method termed TE-qPCR, and Northern blot. Furthermore, we found that Angiogenin and RNase 4 contribute to the biogenesis of sRNA-cPs. Moreover, much more sRNA-cPs than sRNA-OHs bind to Ago2, and can regulate gene expression. Particularly, snR-2-cP regulates Bcl2 by targeting to its 3'UTR dependent on Ago2, and subsequently regulates apoptosis. In addition, sRNA-cPs can guide the cleavage of target RNAs in Ago2 complex as miRNAs without the requirement of 3'-cP. Our discovery greatly expands the repertoire of mammalian sRNAs, and provides strategies and powerful tools towards further investigation of sRNA-cPs.

Men who inject opioids exhibit altered tRNA-Gly-GCC isoforms in semen

In addition to their role in protein translation, tRNAs can be cleaved into shorter, biologically active fragments called tRNA fragments (tRFs). Specific tRFs from spermatocytes can propagate metabolic disorders in second generations of mice. Thus, tRFs in germline cells are a mechanism of epigenetic inheritance. It has also been shown that stress and toxins can cause alterations in tRF patterns. We were therefore interested in whether injecting illicit drugs, a major stressor, impacts tRFs in germline cells. We sequenced RNA from spermatocytes and from semen-derived exosomes from people who inject illicit drugs (PWID) and from non-drug using controls, both groups of unknown fertility status. All PWID injected opioids daily, but most also used other illicit drugs. The tRF cleavage products from Gly-GCC tRNA were markedly different between spermatocytes from PWID compared to controls. Over 90% of reads in controls mapped to shorter Gly-GCC tRFs, while in PWID only 45% did. In contrast, only 4.1% of reads in controls mapped to a longer tRFs versus 45.6% in PWID. The long/short tRF ratio was significantly higher in PWID than controls (0.23 versus 0.16, P = 0.0128). We also report differential expression of a group of small nucleolar RNAs (snoRNAs) in semen-derived exosomes, including, among others, ACA14a, U19, and U3-3. Thus, PWID exhibited an altered cleavage pattern of tRNA-Gly-GCC in spermatocytes and an altered cargo of snoRNAs in semen-derived exosomes. Participants were not exclusively using opioids and were not matched with controls in terms of diet, chronic disease, or other stressors, so our finding are not conclusively linked to opioid use. However, all individuals in the PWID group did inject heroin daily. Our study indicates a potential for opioid injection and/or its associated multi-drug use habits and lifestyle changes to influence epigenetic inheritance.

Emerging roles of tRNA-derived fragments in cancer

tRNA-derived fragments (tRFs) are an emerging category of small non-coding RNAs that are generated from cleavage of mature tRNAs or tRNA precursors. The advance in high-throughput sequencing has contributed to the identification of increasing number of tRFs with critical functions in distinct physiological and pathophysiological processes. tRFs can regulate cell viability, differentiation, and homeostasis through multiple mechanisms and are thus considered as critical regulators of human diseases including cancer. In addition, increasing evidence suggest the extracellular tRFs may be utilized as promising diagnostic and prognostic biomarkers for cancer liquid biopsy. In this review, we focus on the biogenesis, classification and modification of tRFs, and summarize the multifaceted functions of tRFs with an emphasis on the current research status and perspectives of tRFs in cancer.

Integrated analysis of circulating cell free nucleic acids for cancer genotyping and immune phenotyping of tumor microenvironment

The circulating cell-free nucleic acids (ccfNAs) consist of a heterogenous cocktail of both single (ssNA) and double-stranded (dsNA) nucleic acids. These ccfNAs are secreted into the blood circulation by both healthy and malignant cells via various mechanisms including apoptosis, necrosis, and active secretion. The major source of ccfNAs are the cells of hematopoietic system under healthy conditions. These ccfNAs include fragmented circulating cell free DNA (ccfDNA), coding or messenger RNA (mRNA), long non-coding RNA (lncRNA), microRNA (miRNA), and mitochondrial DNA/RNA (mtDNA and mtRNA), that serve as prospective biomarkers in assessment of various clinical conditions. For, e.g., free fetal DNA and RNA migrate into the maternal plasma, whereas circulating tumor DNA (ctDNA) has clinical relevance in diagnostic, prognostic, therapeutic targeting, and disease progression monitoring to improve precision medicine in cancer. The epigenetic modifications of ccfDNA as well as circulating cell-free RNA (ccfRNA) such as miRNA and lncRNA show disease-related variations and hold potential as epigenetic biomarkers. The messenger RNA present in the circulation or the circulating cell free mRNA (ccf-mRNA) and long non-coding RNA (ccf-lncRNA) have gradually become substantial in liquid biopsy by acting as effective biomarkers to assess various aspects of disease diagnosis and prognosis. Conversely, the simultaneous characterization of coding and non-coding RNAs in human biofluids still poses a significant hurdle. Moreover, a comprehensive assessment of ccfRNA that may reflect the tumor microenvironment is being explored. In this review, we focus on the novel approaches for exploring ccfDNA and ccfRNAs, specifically ccf-mRNA as biomarkers in clinical diagnosis and prognosis of cancer. Integrating the detection of circulating tumor DNA (ctDNA) for cancer genotyping in conjunction with ccfRNA both quantitatively and qualitatively, may potentially hold immense promise towards precision medicine. The current challenges and future directions in deciphering the complexity of cancer networks based on the dynamic state of ccfNAs will be discussed.

Circulating messenger RNA variants as a potential biomarker for surveillance of hepatocellular carcinoma

Background and rationale

Liver derived messenger ribonucleic acid (mRNA) transcripts were reported to be elevated in the circulation of hepatocellular carcinoma (HCC) patients. We now report the detection of high-risk mRNA variants exclusively in the circulation of HCC patients. Numerous genomic alleles such as single nucleotide polymorphisms (SNPs), nucleotide insertions and deletions (called Indels), splicing variants in many genes, have been associated with elevated risk of cancer. Our findings potentially offer a novel non-invasive platform for HCC surveillance and early detection.

Approach

RNAseq analysis was carried out in the plasma of 14 individuals with a diagnosis of HCC, 8 with LC and no HCC, and 6 with no liver disease diagnosis. RNA from 6 matching tumors and 5 circulating extracellular vesicle (EV) samples from 14 of those with HCC was also analyzed. Specimens from two cholangiocarcinoma (CCA) patients were also included in our study. HCC specific SNPs and Indels referred as "variants" were identified using GATK HaplotypeCaller and annotated by SnpEff to filter out high risk variants.

Results

The variant calling on all RNA samples enabled the detection of 5.2 million SNPs, 0.91 million insertions and 0.81 million deletions. RNAseq analyses in tumors, normal liver tissue, plasma, and plasma derived EVs led to the detection of 5480 high-risk tumor specific mRNA variants in the circulation of HCC patients. These variants are concurrently detected in tumors and plasma samples or tumors and EVs from HCC patients, but none of these were detected in normal liver, plasma of LC patients or normal healthy individuals. Our results demonstrate selective detection of concordant high-risk HCC-specific mRNA variants in free plasma, plasma derived EVs and tumors of HCC patients. The variants comprise of splicing, frameshift, fusion and single nucleotide alterations and correspond to cancer and tumor metabolism pathways. Detection of these high-risk variants in matching specimens from same subjects with an enrichment in circulating EVs is remarkable. Validation of these HCC selective ctmRNA variants in larger patient cohorts is likely to identify a predictive set of ctmRNA with high diagnostic performance and thus offer a novel non-invasive serology-based biomarker for HCC.

Evaluating the Effects of Storage Conditions on Multiple Cell-Free RNAs in Plasma by High-Throughput Sequencing

Background: Plasma cell-free RNAs (cfRNAs) can serve as noninvasive biomarkers for the diagnosis and monitoring of diseases. However, the delay in blood processing may lead to unreliable results. Therefore, an unbiased evaluation based on the whole transcriptome under different storage conditions is needed. Methods: Here, blood samples were collected in ethylenediaminetetraacetic acid tubes and processed immediately (0 hour), or stored at room temperature (RT) or 4°C for different time intervals (2, 6, and 24 hours) before plasma separation. High-throughput sequencing was applied to assess the effects of storage conditions on the transcript profiles and fragment characteristics of plasma cell-free mRNA, long noncoding RNA (lncRNA), and small RNAs. Results: More genes changed their expression levels with time when blood was stored at RT compared with those at 4°C. Cell-free mRNA and lncRNA were relatively stable in blood preserved at 4°C for 6 hours, while cell-free microRNA (miRNA) and piwi-interacting RNA (piRNA) remained stable at 4°C for 24 hours. After 24 hours, more contamination of the leukocyte-derived RNAs occurred at RT, possibly due to apoptosis. Meanwhile, significant changes were also observed regarding the characteristics of the RNA fragments, including fragment size, the proportion of intron, and the pyrimidine frequency of the fragmented 3' end. Fifteen tissue-enriched genes were detected in the plasma but not expressed in leukocytes. The expression level and fragment length of these genes gradually decreased during storage, suggesting the degradation of the cfRNA and the dilution of leukocyte-derived RNA with other tissue-derived cfRNA. Conclusions: Our results suggest that the contamination of leukocyte-derived RNA and the degradation of original cfRNA contribute to the changes in the cfRNA expression profiles and the fragment characteristics during short-term storage. The storage of blood at 4°C for 6 hours allows plasma cfRNA to remain relatively stable, which will be useful for further studies or clinical applications where adequate quantification or the fragment signature of cfRNA is required.

Polyadenylation ligation-mediated sequencing (PALM-Seq) characterizes cell-free coding and non-coding RNAs in human biofluids

Background

Cell‐free messenger RNA (cf‐mRNA) and long non‐coding RNA (cf‐lncRNA) are becoming increasingly important in liquid biopsy by providing biomarkers for disease prediction, diagnosis and prognosis, but the simultaneous characterization of coding and non‐coding RNAs in human biofluids remains challenging.

Methods

Here, we developed polyadenylation ligation‐mediated sequencing (PALM‐Seq), an RNA sequencing strategy employing treatment of RNA with T4 polynucleotide kinase to generate cell‐free RNA (cfRNA) fragments with 5′ phosphate and 3′ hydroxyl and RNase H to deplete abundant RNAs, achieving simultaneous quantification and characterization of cfRNAs.

Results

Using PALM‐Seq, we successfully identified well‐known differentially abundant mRNA, lncRNA and microRNA in the blood plasma of pregnant women. We further characterized cfRNAs in blood plasma, saliva, urine, seminal plasma and amniotic fluid and found that the detected numbers of different RNA biotypes varied with body fluids. The profiles of cf‐mRNA reflected the function of originated tissues, and immune cells significantly contributed RNA to blood plasma and saliva. Short fragments (<50 nt) of mRNA and lncRNA were major in biofluids, whereas seminal plasma and amniotic fluid tended to retain long RNA. Body fluids showed distinct preferences of pyrimidine at the 3′ end and adenine at the 5′ end of cf‐mRNA and cf‐lncRNA, which were correlated with the proportions of short fragments.

Conclusion

Together, PALM‐Seq enables a simultaneous characterization of cf‐mRNA and cf‐lncRNA, contributing to elucidating the biology and promoting the application of cfRNAs.

Distinct Stress-Dependent Signatures of Cellular and Extracellular tRNA-Derived Small RNAs

The cellular response to stress is an important determinant of disease pathogenesis. Uncovering the molecular fingerprints of distinct stress responses may identify novel biomarkers and key signaling pathways for different diseases. Emerging evidence shows that transfer RNA-derived small RNAs (tDRs) play pivotal roles in stress responses. However, RNA modifications present on tDRs are barriers to accurately quantifying tDRs using traditional small RNA sequencing. Here, AlkB-facilitated methylation sequencing is used to generate a comprehensive landscape of cellular and extracellular tDR abundances in various cell types during different stress responses. Extracellular tDRs are found to have distinct fragmentation signatures from intracellular tDRs and these tDR signatures are better indicators of different stress responses than miRNAs. These distinct extracellular tDR fragmentation patterns and signatures are also observed in plasma from patients on cardiopulmonary bypass. It is additionally demonstrated that angiogenin and RNASE1 are themselves regulated by stressors and contribute to the stress-modulated abundance of sub-populations of cellular and extracellular tDRs. Finally, a sub-population of extracellular tDRs is identified for which AGO2 appears to be required for their expression. Together, these findings provide a detailed profile of stress-responsive tDRs and provide insight about tDR biogenesis and stability in response to cellular stressors.

Exploring the expanding universe of small RNAs

The world of small noncoding RNAs (sncRNAs) is ever-expanding, from small interfering RNA, microRNA and Piwi-interacting RNA to the recently emerging non-canonical sncRNAs derived from longer structured RNAs (for example, transfer, ribosomal, Y, small nucleolar, small nuclear and vault RNAs), showing distinct biogenesis and functional principles. Here we discuss recent tools for sncRNA identification, caveats in sncRNA expression analysis and emerging methods for direct sequencing of sncRNAs and systematic mapping of RNA modifications that are integral to their function.

A Novel Tissue Atlas and Online Tool for the Interrogation of Small RNA Expression in Human Tissues and Biofluids

One promising goal for utilizing the molecular information circulating in biofluids is the discovery of clinically useful biomarkers. Extracellular RNAs (exRNAs) are one of the most diverse classes of molecular cargo, easily assayed by sequencing and with expressions that rapidly change in response to subject status. Despite diverse exRNA cargo, most evaluations from biofluids have focused on small RNA sequencing and analysis, specifically on microRNAs (miRNAs). Another goal of characterizing circulating molecular information, is to correlate expression to injuries associated with specific tissues of origin. Biomarker candidates are often described as being specific, enriched in a particular tissue or associated with a disease process. Likewise, miRNA data is often reported to be specific, enriched for a tissue, without rigorous testing to support the claim. Here we provide a tissue atlas of small RNAs from 30 different tissues and three different blood cell types. We analyzed the tissues for enrichment of small RNA sequences and assessed their expression in biofluids: plasma, cerebrospinal fluid, urine, and saliva. We employed published data sets representing physiological (resting vs. acute exercise) and pathologic states (early- vs. late-stage liver fibrosis, and differential subtypes of stroke) to determine differential tissue-enriched small RNAs. We also developed an online tool that provides information about exRNA sequences found in different biofluids and tissues. The data can be used to better understand the various types of small RNA sequences in different tissues as well as their potential release into biofluids, which should help in the validation or design of biomarker studies.

exoRBase 2.0: an atlas of mRNA, lncRNA and circRNA in extracellular vesicles from human biofluids

Extracellular vesicles (EVs) are small membranous vesicles that contain an abundant cargo of different RNA species with specialized functions and clinical implications. Here, we introduce an updated online database (http://www.exoRBase.org), exoRBase 2.0, which is a repository of EV long RNAs (termed exLRs) derived from RNA-seq data analyses of diverse human body fluids. In exoRBase 2.0, the number of exLRs has increased to 19 643 messenger RNAs (mRNAs), 15 645 long non-coding RNAs (lncRNAs) and 79 084 circular RNAs (circRNAs) obtained from ∼1000 human blood, urine, cerebrospinal fluid (CSF) and bile samples. Importantly, exoRBase 2.0 not only integrates and compares exLR expression profiles but also visualizes the pathway-level functional changes and the heterogeneity of origins of circulating EVs in the context of different physiological and pathological conditions. Our database provides an attractive platform for the identification of novel exLR signatures from human biofluids that will aid in the discovery of new circulating biomarkers to improve disease diagnosis and therapy.

PANDORA-seq expands the repertoire of regulatory small RNAs by overcoming RNA modifications

Although high-throughput RNA sequencing (RNA-seq) has greatly advanced small non-coding RNA (sncRNA) discovery, the currently widely used complementary DNA library construction protocol generates biased sequencing results. This is partially due to RNA modifications that interfere with adapter ligation and reverse transcription processes, which prevent the detection of sncRNAs bearing these modifications. Here, we present PANDORA-seq (panoramic RNA display by overcoming RNA modification aborted sequencing), employing a combinatorial enzymatic treatment to remove key RNA modifications that block adapter ligation and reverse transcription. PANDORA-seq identified abundant modified sncRNAs-mostly transfer RNA-derived small RNAs (tsRNAs) and ribosomal RNA-derived small RNAs (rsRNAs)-that were previously undetected, exhibiting tissue-specific expression across mouse brain, liver, spleen and sperm, as well as cell-specific expression across embryonic stem cells (ESCs) and HeLa cells. Using PANDORA-seq, we revealed unprecedented landscapes of microRNA, tsRNA and rsRNA dynamics during the generation of induced pluripotent stem cells. Importantly, tsRNAs and rsRNAs that are downregulated during somatic cell reprogramming impact cellular translation in ESCs, suggesting a role in lineage differentiation.

Revisiting Extracellular RNA Release, Processing, and Function

It is assumed that RNAs enriched in extracellular samples were selected for release by their parental cells. However, recent descriptions of extracellular RNA (exRNA) biogenesis and their differential stabilities question this assumption, as they could produce identical outcomes. Here, we share our opinion about the importance of considering both selective and nonselective mechanisms for RNA release into the extracellular environment. In doing so, we provide new perspectives on RNA-mediated intercellular communication, including an analogy to communication through social media. We also argue that technical limitations have restricted the study of some of the most abundant exRNAs, both inside and outside extracellular vesicles (EVs). These RNAs may be better positioned to induce a response in recipient cells compared with low abundance miRNAs.

Fragmentation of extracellular ribosomes and tRNAs shapes the extracellular RNAome

A major proportion of extracellular RNAs (exRNAs) do not copurify with extracellular vesicles (EVs) and remain in ultracentrifugation supernatants of cell-conditioned medium or mammalian blood serum. However, little is known about exRNAs beyond EVs. We have previously shown that the composition of the nonvesicular exRNA fraction is highly biased toward specific tRNA-derived fragments capable of forming RNase-protecting dimers. To solve the problem of stability in exRNA analysis, we developed a method based on sequencing the size exclusion chromatography (SEC) fractions of nonvesicular extracellular samples treated with RNase inhibitors (RI). This method revealed dramatic compositional changes in exRNA population when enzymatic RNA degradation was inhibited. We demonstrated the presence of ribosomes and full-length tRNAs in cell-conditioned medium of a variety of mammalian cell lines. Their fragmentation generates some small RNAs that are highly resistant to degradation. The extracellular biogenesis of some of the most abundant exRNAs demonstrates that extracellular abundance is not a reliable input to estimate RNA secretion rates. Finally, we showed that chromatographic fractions containing extracellular ribosomes are probably not silent from an immunological perspective and could possibly be decoded as damage-associated molecular patterns.

EV-origin: Enumerating the tissue-cellular origin of circulating extracellular vesicles using exLR profile

Extracellular vesicles (EVs) are complex ecosystems that can be derived from all body cells and circulated in the body fluids. Characterizing the tissue-cellular source contributing to circulating EVs provides biological information about the cell or tissue of origin and their functional states. However, the relative proportion of tissue-cellular origin of circulating EVs in body fluid has not been thoroughly characterized. Here, we developed an approach for digital EVs quantification, called EV-origin, that enables enumerating of EVs tissue-cellular source contribution from plasma extracellular vesicles long RNA sequencing profiles. EV-origin was constructed by the input matrix of gene expression signatures and robust deconvolution algorithm, collectively used to separate the relative proportions of each tissue or cell type of interest. EV-origin respectively predicted the relative enrichment of seven types of hemopoietic cells and sixteen solid tissue subsets from exLR-seq profile. Using the EV-origin approach, we depicted an integrated landscape of the traceability system of plasma EVs for healthy individuals. We also compared the heterogenous tissue-cellular source components from plasma EVs samples with diverse disease status. Notably, the aberrant liver fraction could reflect the development and progression of hepatic disease. The liver fraction could also serve as a diagnostic indicator and effectively separate HCC patients from normal individuals. The EV-origin provides an approach to decipher the complex heterogeneity of tissue-cellular origin in circulating EVs. Our approach could inform the development of exLR-based applications for liquid biopsy.

Identification of protein-protected mRNA fragments and structured excised intron RNAs in human plasma by TGIRT-seq peak calling

Human plasma contains > 40,000 different coding and non-coding RNAs that are potential biomarkers for human diseases. Here, we used thermostable group II intron reverse transcriptase sequencing (TGIRT-seq) combined with peak calling to simultaneously profile all RNA biotypes in apheresis-prepared human plasma pooled from healthy individuals. Extending previous TGIRT-seq analysis, we found that human plasma contains largely fragmented mRNAs from > 19,000 protein-coding genes, abundant full-length, mature tRNAs and other structured small non-coding RNAs, and less abundant tRNA fragments and mature and pre-miRNAs. Many of the mRNA fragments identified by peak calling correspond to annotated protein-binding sites and/or have stable predicted secondary structures that could afford protection from plasma nucleases. Peak calling also identified novel repeat RNAs, miRNA-sized RNAs, and putatively structured intron RNAs of potential biological, evolutionary, and biomarker significance, including a family of full-length excised intron RNAs, subsets of which correspond to mirtron pre-miRNAs or agotrons.

Processing by RNase 1 forms tRNA halves and distinct Y RNA fragments in the extracellular environment

Extracellular RNAs participate in intercellular communication, and are being studied as promising minimally invasive diagnostic markers. Several studies in recent years showed that tRNA halves and distinct Y RNA fragments are abundant in the extracellular space, including in biofluids. While their regulatory and diagnostic potential has gained a substantial amount of attention, the biogenesis of these extracellular RNA fragments remains largely unexplored. Here, we demonstrate that these fragments are produced by RNase 1, a highly active secreted nuclease. We use RNA sequencing to investigate the effect of a null mutation of RNase 1 on the levels of tRNA halves and Y RNA fragments in the extracellular environment of cultured human cells. We complement and extend our RNA sequencing results with northern blots, showing that tRNAs and Y RNAs in the non-vesicular extracellular compartment are released from cells as full-length precursors and are subsequently cleaved to distinct fragments. In support of these results, formation of tRNA halves is recapitulated by recombinant human RNase 1 in our in vitro assay. These findings assign a novel function for RNase 1, and position it as a strong candidate for generation of tRNA halves and Y RNA fragments in biofluids.

Extracellular tRNAs and tRNA-derived fragments

Fragmentation of tRNAs generates a family of small RNAs collectively known as tRNA-derived fragments. These fragments vary in sequence and size but have been shown to regulate many processes involved in cell homoeostasis and adaptations to stress. Additionally, the field of extracellular RNAs (exRNAs) is rapidly growing because exRNAs are a promising source of biomarkers in liquid biopsies, and because exRNAs seem to play key roles in intercellular and interspecies communication. Herein, we review recent descriptions of tRNA-derived fragments in the extracellular space in all domains of life, both in biofluids and in cell culture. The purpose of this review is to find consensus on which tRNA-derived fragments are more prominent in each extracellular fraction (including extracellular vesicles, lipoproteins and ribonucleoprotein complexes). We highlight what is becoming clear and what is still controversial in this field, in order to stimulate future hypothesis-driven studies which could clarify the role of full-length tRNAs and tRNA-derived fragments in the extracellular space.

Profiling the circulating mRNA transcriptome in human liver disease

The human circulation contains cell-free DNA and non-coding microRNA (miRNA). Less is known about the presence of messenger RNA (mRNA). This report profiles the human circulating mRNA transcriptome in people with liver cirrhosis (LC) and hepatocellular carcinoma (HCC) to determine whether mRNA analytes can be used as biomarkers of liver disease. Using RNAseq and RT-qPCR, we investigate circulating mRNA in plasma from HCC and LC patients and demonstrate detection of transcripts representing more than 19,000 different protein coding genes. Remarkably, the circulating mRNA expression levels were similar from person to person over the 21 individuals whose samples were analyzed by RNAseq. Liver derived circulating transcripts such as albumin (ALB), apolipoprotein (APO) A1, A2 & H, serpin A1 & E1, ferritin light chain (FTL) and fibrinogen like 1 (FGL1) were significantly upregulated in HCC patient samples. Higher levels of some of these liver-specific transcripts in the plasma of HCC patients were confirmed by RT-qPCR in another cohort of 20 individuals. Several less abundant circulating transcripts associated with cancer were detected in most HCC samples, but not in healthy subjects. Liver specificity of circulating transcripts was confirmed by investigating their expression in HCC tumor and liver cancer cell lines. Liver specific mRNA sequences in the plasma were predominantly present outside circulating extracellular vesicles.

Conclusions: The circulating “mRNA” transcriptome is remarkably consistent in diversity and expression from person to person. Detection of transcripts corresponding to disease selective polypeptides suggests the possibility that circulating mRNA can work as a biomarker analyte for cancer detection.

Performance assessment of total RNA sequencing of human biofluids and extracellular vesicles

RNA profiling has emerged as a powerful tool to investigate the biomarker potential of human biofluids. However, despite enormous interest in extracellular nucleic acids, RNA sequencing methods to quantify the total RNA content outside cells are rare. Here, we evaluate the performance of the SMARTer Stranded Total RNA-Seq method in human platelet-rich plasma, platelet-free plasma, urine, conditioned medium, and extracellular vesicles (EVs) from these biofluids. We found the method to be accurate, precise, compatible with low-input volumes and able to quantify a few thousand genes. We picked up distinct classes of RNA molecules, including mRNA, lncRNA, circRNA, miscRNA and pseudogenes. Notably, the read distribution and gene content drastically differ among biofluids. In conclusion, we are the first to show that the SMARTer method can be used for unbiased unraveling of the complete transcriptome of a wide range of biofluids and their extracellular vesicles.

Comparative Analysis of Biochemical Biases by Ligation- and Template-Switch-Based Small RNA Library Preparation Protocols

BACKGROUND

Small RNAs are key players in the regulation of gene expression and differentiation. However, many different classes of small RNAs (sRNAs) have been described with distinct biogenesis pathways and, as a result, with different biochemical properties. To analyze sRNAs by deep sequencing, complementary DNA synthesis requires manipulation of the RNA molecule itself. Thus, enzymatic activities during library preparation bias the library content owing to biochemical criteria.

METHODS

We compared 4 different manipulations of RNA for library preparation: (a) a ligation-based procedure allowing only 5'-mono-phosphorylated RNA to enter the library, (b) a ligation-based procedure allowing additional 5'-triphosphates and Cap structures, (c) a ligation-independent, template-switch-based library preparation, and (d) a template-switch-based library preparation allowing 3'-phosphorylated RNAs to enter the library.

RESULTS

Our data show large differences between ligation-dependent and ligation-independent libraries in terms of their preference for individual sRNA classes such as microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs), and transfer RNA fragments. Moreover, the miRNA composition is different between both procedures, and more microRNA isoforms (isomiRs) can be identified after pyrophosphatase treatment. piRNAs are enriched in template-switch libraries, and this procedure apparently includes more different RNA species.

CONCLUSIONS

Our data indicate that miRNAomics from both methods will hardly be comparable. Ligation-based libraries enrich for canonical miRNAs, which thus may be suitable methods for miRNAomics. Template-switch libraries contain increased numbers and different compositions of fragments and long RNAs. Following different interests for other small RNA species, ligation-independent libraries appear to show a more realistic sRNA landscape with lower bias against biochemical modifications.