Skip to main content
Download PDF

Diversity and signature of small RNA in different bodily fluids using next generation sequencing

BMC Genomics volume 19, Article number: 408 (2018) Cite this article

Abstract

Background

Small RNAs are critical components in regulating various cellular pathways. These molecules may be tissue-associated or circulating in bodily fluids and have been shown to associate with different tumors. Next generation sequencing (NGS) on small RNAs is a powerful tool for profiling and discovery of microRNAs (miRNAs).

Results

In this study, we isolated total RNA from various bodily fluids: blood, leukocytes, serum, plasma, saliva, cell-free saliva, urine and cell-free urine. Next, we used Illumina’s NGS platform and intensive bioinformatics analysis to investigate the distribution and signature of small RNAs in the various fluids. Successful NGS was accomplished despite the variations in RNA concentrations among the different fluids. Among the fluids studied, blood and plasma were found to be the most promising fluids for small RNA profiling as well as novel miRNA prediction. Saliva and urine yielded lower numbers of identifiable molecules and therefore were less reliable in small RNA profiling and less useful in predicting novel molecules. In addition, all fluids shared many molecules, including 139 miRNAs, the most abundant tRNAs, and the most abundant piwi-interacting RNAs (piRNAs). Fluids of similar origin (blood, urine or saliva) displayed closer clustering, while each fluid still retains its own characteristic signature based on its unique molecules and its levels of the common molecules. Donor urine samples showed sex-dependent differential clustering, which may prove useful for future studies.

Conclusions

This study shows the successful clustering and unique signatures of bodily fluids based on their miRNA, tRNA and piRNA content. With this information, cohorts may be differentiated based on multiple molecules from each small RNA class by a multidimensional assessment of the overall molecular signature.

Background

Small RNAs are a class of mainly non-coding RNAs (ncRNAs) characterized by their small nucleotide length of less than 200 nt [1]. Within this class there are key RNA types with a size range of 14–35 nt that are highly important for diagnostic biomarker discovery and the development of therapeutic agents [2,3,4]. These include microRNAs (miRNAs), transfer RNA-derived RNAs (tDRs) and Piwi-interacting RNAs (piRNAs). MicroRNAs are non-coding molecules of about 19–23 nt that bind to and downregulate messenger RNAs (mRNAs) [5]. They down regulate gene expression, playing a major role in essential biological pathways, such as differentiation, proliferation, metastasis and apoptosis. [6,7,8,9,10,11]. MicroRNAs represent an entire layer of gene expression regulation, regulating more than 50% of protein coding mRNAs in mammalian cells [12]. To date, 2588 human mature miRNAs have been identified and are currently included in miRBase 21 [13]. Aside from being found in tissues and cells, miRNAs are found in bodily fluids in extracellular vesicles or in complexes with argonaute or lipoproteins [14,15,16,17]. They have been reported in bodily fluids such as blood, plasma, serum, urine, tears, saliva, breast milk, amniotic fluid, seminal fluid and colostrum [18,19,20,21,22,23]. MicroRNAs have been linked to many diseases and are highly promising molecular biomarkers [24,25,26,27].

Mature tRNAs and nascent pre-tRNA transcripts are processed enzymatically to produce well defined tDRs in a regulated process, suggesting that they are not random degradation products [28, 29]. Sizes of tDRs range from 30 to 35 nt for tRNA halves and 14 to 26 nt for the shorter fragments [4]. Various studies suggest that tDRs are involved in different functions including stress responses in human diseases, where they act as inhibitors of global translation and transcription regulation [30,31,32,33,34,35,36]. Like miRNAs, they can conduct specific gene silencing and have potential as cancer biomarkers [7, 37,38,39,40,41,42,43,44,45,46,47]. Finally, piRNAs have a nucleotide size range between 26 to 31. They modulate different gene expression pathways by interacting with Piwi proteins [48, 49]. Currently there are 23,439 piRNA molecules in the piRNABank [50]. They are abundant in gonads and mediate transposon repression to conserve genome integrity [51,52,53].

The use of next generation sequencing (NGS) technology in small RNA detection has advanced research in the field at unprecedented speed. NGS shines light on the key role of ncRNAs in transcriptome regulation in healthy and disease conditions and accelerates the profiling and discovery of molecules [54]. The technology enables the analysis of multiple samples in parallel and provides precise quantification of each molecule, making it superior to previous genomic technologies. The demonstrated capacities of NGS have led to advances in biological and medical genomics and transcriptomics [55,56,57].

Solid tumor miRNAs are well represented in bodily fluids, indicating their importance as cancer biomarkers [58,59,60]. Almost all bodily fluids from healthy individuals contain miRNAs. Therefore, bodily fluids represent an excellent candidate for non-invasive detection of miRNAs and have been used in applications such as cancer biomarker discovery [8]. Transfer RNA-derived small RNAs are thought to have a dual function, as they act as suppressors and oncogenes [42]. In addition, altered piRNAs levels were found to be associated with lung, breast, gastric and colon cancers [61,62,63,64]. However, no comprehensive study has been reported on tDRs or piRNAs in different bodily fluids.

The field of small RNA is expanding with the profiling and discovery of molecules in various disease conditions and treatments. Therefore, it is important to explore the small RNA content in normal individuals to better understand the small RNA profile in each fluid as well as their relative distribution among the different fluids. To gain insights into the distribution and signature of small RNA in bodily fluids, we carried out a comparative study on RNA from different fluids collected from the same donors and used NGS to explore and describe their small RNA content.

Methods

Sample collection, preparation and RNA purification

Blood, saliva and urine were collected from 4 healthy donors, 2 females and 2 males between the ages of 20 to 30. The various bodily fluids were collected from each individual within a 2-h period. Collection and sample de-identification was performed under an IRB approved protocol (16198–16:02:416–11-2017). Three 10 mL blood samples were collected from each donor. Two of these samples were collected in Vacutainer® plastic EDTA tubes (BD, USA) and used for RNA isolation from whole blood, leukocytes and plasma. RNA was isolated directly from 0.2 mL of blood using the Total RNA Purification Plus Kit (Norgen Biotek Corp., Canada). Leukocytes were prepared from 0.5 mL of blood by using the Leukocyte RNA Purification Kit (Norgen Biotek Corp., Canada). One entire tube was used to prepare plasma and was centrifuged at 200 RCF for 10 min at room temperature. Plasma was collected and stored at − 70 °C until isolation. The last blood sample was collected in a Vacutainer® glass serum tube with silicon coated interior (BD, USA) and used for RNA isolation from serum. The tube was left to stand at room temperature for 45 min, and then it was centrifuged at 1300 RCF for 15 min. Serum was collected and stored at − 70 °C until isolation. Both plasma and serum RNA were isolated from 0.2 mL using the Plasma/Serum RNA Isolation Mini Kit (Norgen Biotek Corp., Canada). All kits were used according to the manufacturer’s instructions.

Two milliliters of saliva was collected from each donor into Falcon 50 mL centrifuge tubes (BD, USA) and 0.3 mL was used directly for saliva RNA isolation using the Total RNA Purification Kit (Norgen Biotek Corp., Canada), according to the manufacturer’s instructions. Another 0.3 mL of saliva was transferred into a 1.5 mL tube (Eppendorf, Germany) and spun down at 200 RCF for 10 min to remove cells, and then the supernatant was used for cell-free saliva RNA purification by the same kit. A similar approach was used for urine sample preparation, where 100 mL of urine was collected from each donor into disposable cups (Sarstedt, Germany). RNA was isolated from 30 mL of urine using the Urine Cell-Free Circulating RNA Purification Maxi Kit (Norgen Biotek Corp., Canada). The kit’s procedure was modified by skipping the initial centrifugation steps to purify RNA from total urine. Another urine sample was processed by the same kit without any modifications to the manufacturer’s procedures to isolate RNA from 30 mL of cell-free urine. Purified RNA from all samples were tested for positive amplification by miR-21 stem-loop RT-PCR [65]. RNA concentration was then estimated by the Agilent 2100 Bioanalyzer System (Agilent Technologies, USA) using the RNA 6000 Nano Total RNA chip.

Small RNA library construction and high-throughput sequencing

The small RNA libraries were prepared from the RNA isolated from each sample using the Small RNA Library Prep Kit for Illumina (Norgen Biotek Corp., Canada) according to the manufacturer’s instructions. Briefly, 6 μL of purified RNA was mixed with the 3′ adapter and incubated at 70 °C for 2 min before being used in a ligation step by adding T4 RNA ligase 2 (truncated), buffer and RNase inhibitor. The reaction was incubated at 28 °C for 1 h then heat inactivated at 70 °C for 10 min. The excess 3′ adapters were removed by the addition of the reverse primer and incubating the reaction at 75 °C for 5 min, 37 °C for 15 min and 25 °C for 5 min. The 5′ adapter was denatured at 70 °C for 2 min and then added together with 10 mM ATP and T4 RNA ligase 1 to the reaction and incubated at 28 °C for 1 h followed by heat inactivation at 70 °C for 10 min. The two adapters were diluted 1:1 before being added to the reactions and all the incubation steps were performed in a thermocycler with cooling on ice between the different steps. Reverse transcription was performed on the ligation reaction product by adding a mixture containing 10 mM dNTPs, first strand buffer and TruScript reverse transcriptase, and incubating the reaction at 50 °C for 1 h before heat inactivation at 70 °C for 15 min. Finally, the reverse transcription reaction product was amplified and indexed in a 15 cycle PCR reaction by adding the NGS PCR master mix, PCR reverse primer and the unique index primer for each sample.

The PCR reaction product was cleaned and separated on a 6% Novex® TBE PAGE gel (Life Technologies, USA). The gel was stained with SYBR® Gold Nucleic Acid Gel Stain (Life Technologies, USA) and a library size range from 125 bp to 170 bp was excised from the gel and placed in a Gel Breaker Tube (IST Engineering, USA), then centrifuged at 14000 RCF for 2 min. The prepared libraries were then eluted overnight in nuclease-free water (Ambion, USA) and cleaned. The library was quantified by the High Sensitivity DNA Analysis Kit on the Agilent 2100 Bioanalyzer System (Agilent Technologies, USA). Libraries were diluted to 4 nM, pooled, and sequenced on the Illumina HiSeq 4000 at The McGill University and Génome Québec Innovation Centre (Montreal, Canada), using the HiSeq 3000/4000 SBS Kit (50 cycles).

Read mapping and small RNA annotation

The sequence raw data from the Illumina HiSeq 4000 were converted to fastq format. Files were then used in the Genboree Workbench’s exceRpt small RNA-seq pipeline (version 4.6.2) for read mapping to the hg38 human genome version [66]. This allowed for a single mismatched base down to 18 nucleotides. After adapter trimming, read quality was assessed by FASTQC to filter out reads with a quality score lower than 30 on the PHRED scale. Reads were first mapped to the UniVec and human ribosomal RNA (rRNA) sequences to exclude them before mapping to databases of miRBase version 21, gtRNAdb and piRNABank to assign reads to miRNAs, tRNAs and piRNAs, respectively. Identified tRNAs are tRNA-derived RNA fragments due to the fact that the library insert size is below 50 nt. Remaining sequences were then annotated to gencode version 24 (hg38) which includes protein coding transcripts (protein_coding), mitochondrial rRNA (Mt_rRNA), mitochondrial tRNA (Mt_tRNA), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), long intergenic noncoding RNA (lincRNA) and miscellaneous RNA (misc_RNA).

Data analysis

Raw read counts obtained from the Genboree Workbench’s exceRpt small RNA-seq pipeline were further analyzed using R (version 3.4.0) and R studio (version 1.0.143). The following R packages were used in the analysis: RnaSeqGeneEdgeRQL (version 1.0.0) for counts per million (CPM) filtration and normalization by using trimmed mean of M-values (TMM) [67], ggfortify (version 0.4.1) and ComplexHeatmap (version 1.14.0) for principal component analysis (PCA) plot and heatmaps based on the filtered and normalized data, respectively. VennDiagram (version 1.6.17) was used to illustrate Venn diagrams. miRDeep2 (version 2.0.0.8) was used to predict novel miRNA candidates and tDRmapper was used to identify tDRs.

Results

Small RNA profiles in the various bodily fluids used in this study provide an atlas of miRNAs, tRNAs and piRNAs relative distribution. They also provide in depth molecular analysis and a guide for NGS-based small RNA expression studies that employ one or more of these bodily fluids as a source of biological data. It is important to look at the normal characteristics of small RNA molecules in each fluid in terms of abundance and representation. The origin and nature of these fluids can pose a significant effect on their use in certain studies that might require specific handling during preparation and sequencing to ensure the validity of results.

RNA concentration variations in the different bodily fluids

Concentration of RNA from each bodily fluid tested was measured using an Agilent 2100 Bioanalyzer. The average range of RNA content in 1 L of bodily fluids was as low as 0.01 mg in urine to as high as 11.2 mg in saliva. Bodily fluids can be categorized based on their RNA content; significantly higher amounts of RNA can be recovered from saliva, blood and cell-free saliva (4.2 to 11.2 mg/L). Leukocytes, serum and plasma had moderate yields of 0.8 to 1.8 mg/L, while urine and cell-free urine had significantly lower RNA content of 0.01 mg/L. Blood, leukocytes, saliva and cell-free saliva had lower standard deviations in their RNA content (< 50% of average), whereas plasma, serum, urine and cell-free urine had higher deviations between their samples (70–85% of average). The average concentration of the isolated RNA from all bodily fluids ranged from 67.2 ng/μL to 3.4 ng/μL. They can be classified into high (> 20 ng/μL) from blood, leukocytes, saliva and cell-free saliva and low (< 10 ng/μL) from plasma, serum, urine and cell-free urine. The RNA integrity number (RIN) was more than 7 for RNA from leukocytes and lower for RNA from blood, saliva and cell-free saliva (about 2–3). RNA from serum, urine and cell-free urine had a low RIN of 1 or less, while RNA from plasma had no measurable RIN from any sample (Table 1).

Table 1 Variations in RNA concentration, RIN value, and yield among the different bodily fluids
Full size table

Input read alignment

Reads obtained from sequencing were used for alignment and mapping to the human genome after adapter clipping and quality filtering. The range of average input reads from the various bodily fluids was between 9.5 million reads (serum) to 15.7 million reads (blood). The average input reads from all the bodily fluid samples tested was 12.57 ± 3.54 million reads. The descending order of samples based on their number of input reads was: blood, cell-free urine, urine, leukocyte, plasma, cell-free saliva, saliva and serum. The percentage of successfully clipped reads was more than 60% from all sample types, with a minimum percentage of reads failing quality filters. Reads were mapped to human rRNA to exclude rRNA sequences before mapping to human genome. The average percentage of reads aligned to human rRNA was less than 12% in all bodily fluids except saliva, cell-free saliva and leukocytes, which had average percentages of 16.6 ± 7.5, 13.0 ± 6.2 and 36.0 ± 1.5, respectively. This is indirectly proportional to reads used for alignment to the human genome. More than 50% of reads used for alignment were mapped to the human genome in blood, plasma, serum, urine and cell-free urine. The percentage was lower in the leukocytes as well as total and cell-free samples of saliva.

In saliva and cell-free saliva, the average percentage of unmapped reads was about 50% of reads used for alignment (48.9 ± 19.7% and 50.3 ± 10.1% of input reads, respectively). Conversely, urine and cell-free urine had an average percentage of unmapped reads relative to reads used in alignment of 21.2 ± 17.8% and 25.5 ± 23.9%, respectively. The percentage of input reads alignment from each bodily fluid can be found in Table 2.

Table 2 Percentage of input reads alignment from each bodily fluid
Full size table

Small RNA biotype mapping

Reads that were mapped to the human genome were then mapped and classified to the various small RNA biotypes. The average total reads mapped to small RNA biotypes within each bodily fluid ranged from 1.3 to 12.8 million reads. Blood, plasma, cell-free urine and urine had more than 8 million reads mapped to biotypes (12.8, 9.7, 8.7 and 8.0 million reads, respectively). Serum, leukocytes, cell-free saliva and saliva had 5.4 million reads or less (5.4, 2.6, 1.3 and 1.3 million reads, respectively). The distribution of biotypes within each bodily fluid showed distinct patterns. Plasma had a high percentage of miscellaneous RNA (misc_RNA; 58.0 ± 39.4), while urine and cell-free urine had high amounts of tRNAs (91.3 ± 77.5% and 91.3 ± 90.3%, respectively). The other bodily fluids had a more diverse pattern with no single biotype exceeding 50% of the content. MicroRNAs represented more than 85% of blood biotypes, 25% of leukocytes, and 15–25% of plasma, serum and cell-free saliva. Saliva, cell-free urine and urine contained the lowest miRNA content (5.3–12.0%). Transfer RNA was the predominant biotype in urine and cell-free urine (> 90%), while serum, saliva and cell-free saliva contained moderate tRNA content (20–50%). Leukocytes had 18.4 ± 4.0 tRNA while plasma and blood contained the lowest tRNA fractions of 5.8 ± 2.5 and 2.1 ± 0.7, respectively. Finally, piRNAs represented less than 2% of the reads in blood, serum, saliva, cell-free saliva, urine and cell-free urine, while more than 5% of the reads were piRNAs in leukocytes and plasma. The percentages of the various biotypes in each bodily fluid are listed in Table 3 and illustrated in Fig. 1. The biotype distribution in each donor is illustrated in an additional file (Additional file 1: Figure S2) and shows a relatively similar pattern for each bodily fluid between the donors. However, in urine samples there was a difference in patterns between male and female donors.

Table 3 Percentage of biotype counts in the various bodily fluids
Full size table
Fig. 1
figure1

Relative biotype distribution among the various bodily fluids. The graph was generated using average percentage of biotype counts of each fluid. MicroRNAs are the largest biotype of blood, while tRNAs are the major biotype of urine. Saliva has the highest biotype diversity among all fluids

Full size image

miRNA

The most variable 50 miRNAs were calculated based on TMM-normalized miRNA counts at a CPM corresponding to a minimum of 5 counts in a library to achieve a high confidence level. These miRNAs were then used for PCA which showed the separation of various bodily fluids based on their miRNA expression (Fig. 2). The analysis revealed closeness between saliva and cell-free saliva samples and between urine and cell-free urine samples. Close clustering was seen between blood and leukocytes, and between plasma and serum. The data disclosed consistent biological origin and miRNA expression-based separation of bodily fluid profiles. The Z scores of the most variable 50 miRNAs were used to generate a heatmap that illustrated the pattern of expression as well as the relationship between samples (Fig. 3). It clustered bodily fluids based on their biology. The dendrogram showed that invasive bodily fluids (blood, leukocytes, plasma and serum) branched apart from non-invasive fluids (saliva, cell-free saliva, urine and cell-free urine). Furthermore, it showed that cell-depletion from saliva or urine did not have a major effect on their clustering. A large set of miRNAs appeared to be highly abundant or severely depleted in the various bodily fluids. Particularly, urine and cell-free urine as well as plasma and serum had different sets of upregulated miRNAs.

Fig. 2
figure2

Principal component analysis of the most variable miRNAs in each bodily fluid. Analysis of the most variable 50 miRNAs was calculated based on TMM-normalized miRNA counts. Four pairs of fluids show close clustering: blood/leukocyte, plasma/serum, saliva/cell-free saliva and urine/cell-free urine

Full size image
Fig. 3
figure3

Heatmap clustering of the most variable miRNAs in each of the bodily fluids. The sex of the sample donor is indicated as (F) for female donors or (M) for male donors. The analysis was generated using Z scores of the most variable 50 miRNAs. The dendrogram shows distinct clustering of invasive fluids (blood, leukocyte, plasma and serum) and non-invasive fluids (saliva, cell-free saliva, urine and cell-free urine)

Full size image

Bodily fluids can be classified into two groups based on their collection procedures: invasive (blood, leukocytes, plasma and serum) and non-invasive (saliva and urine). Variations can be seen in the detectable number of miRNAs, at a minimum of 5 counts in 3 or more individuals, between the two groups. The range of detected miRNAs from invasive fluids was 307–440 while the range from the non-invasive fluids was 178–233. Blood had the largest number of detected miRNAs and saliva had the lowest. Plasma contained more miRNAs compared to serum. Similarly, more miRNAs were detected from the cell-free preparation of saliva than saliva. Almost no differences could be seen between urine and cell-free urine (Table 4). The overlap between detected miRNAs in the different bodily fluids was illustrated in Venn diagrams (Fig. 4). About 97% of serum miRNAs were shared with plasma. More than 90% of miRNAs in saliva and cell-free saliva were shared with blood, leukocytes and plasma. Saliva had 98.9% of its miRNA identical to cell-free saliva, while the latter had only 75.5% of its miRNAs overlapped with saliva. About 85–91% of urine and cell-free urine miRNAs were overlapped with blood, leukocytes, plasma and serum. Cell-free urine had 91.3% similarity with urine, while urine had 92.7% similarity to cell-free urine. Saliva and urine shared more than 77% of their miRNAs. In addition, Venn diagrams were used to demonstrate the overlap between invasive and non-invasive bodily fluids (Fig. 5). The invasive fluids had 230 common miRNAs (Additional file 2: Table S1). Blood had 98 unique miRNAs, which was 2 to 3-fold higher than plasma and leukocytes. In contrast, non-invasive fluids had 148 common miRNAs (Additional file 3: Table S2), and had lower numbers of unique miRNAs. The non-invasive fluids shared 144 common miRNAs with blood (Additional file 4: Table S3), while the latter had 209 unique miRNAs that were absent from all non-invasive fluids (Additional file 5: Figure S1). Common miRNAs between all fluids were 139 (Additional file 6: Table S4).

Table 4 Detected number of miRNAs in the various bodily fluids
Full size table
Fig. 4
figure4

Overlap of miRNA content between various bodily fluids. MicroRNAs of each fluid were filtered to keep molecules that have a minimum of 5 counts in 3 or more individuals. The highest overlap is seen between fluids within the same category: invasive fluids (blood, leukocyte, plasma and serum) or non-invasive fluids (saliva, cell-free saliva, urine and cell-free urine)

Full size image
Fig. 5
figure5

Venn diagram showing the overlap between (a) invasive and (b) non-invasive bodily fluids. Only miRNAs that are present at a minimum of 5 counts in 3 or more individuals were included in the comparison. Invasive fluids have higher number of shared and unique miRNAs compared to non-invasive fluids

Full size image

We ran a detailed analysis on the miRNA composition of each bodily fluid. We looked at the 20 most abundant miRNAs and calculated their fractions of the total miRNA content of each sample (Table 5). These 20 most abundant miRNAs covered about 74 to 94% of all miRNA counts. They represented 94 and 91% of blood and cell-free urine, respectively. They were only 75% of all miRNA counts in saliva and cell-free saliva, while representing 80–88% of the rest of the other bodily fluids. The deep analysis of miRNAs revealed that 2 specific miRNAs were dominating the counts of 4 bodily fluids. Hsa-miR-486-5p made up 53.6 ± 1.9% and 43.2 ± 13.2% of miRNA counts of blood and serum, respectively, while hsa-miR-10b-5p represented 38.4 ± 8.3% and 45.6 ± 5.6% of miRNA counts of urine and cell-free urine, respectively. Other miRNAs that represented 10–20% of total miRNA counts included let-7f-5p (11.11 ± 1.12%) in blood, miR-146b-5p (11.23 ± 1.85) and let-7f-5p (10.58 ± 1.33) in leukocytes, miR-486-5p (15.16 ± 1.63) and miR-191-5p (10.84 ± 0.25) in plasma, miR-143-3p (10.65 ± 4.16) in saliva, miR-143-3p (14.92 ± 4.47) and miR-191-5p (11.59 ± 1.64) in cell-free saliva, and miR-10a-5p in urine (11.56 ± 2.67) and cell-free urine (13.6 ± 0.97). The proportions of the top 10 most abundant miRNAs in each fluid are illustrated in Fig. 6. Among the top 20 most abundant miRNAs of each fluid, 5 miRNAs were found common to all fluids. These include: hsa-let-7a-5p, hsa-let-7f-5p, hsa-miR-191-5p, hsa-miR-26a-5p and hsa-miR-486-5p. These five miRNAs represent more than 50% of blood and serum miRNA counts, 25 to 45% of plasma and saliva (total and cell-free) and less than 11% of urine and cell-free urine miRNA counts (Fig. 7).

Table 5 Twenty most abundant miRNAs detected in each bodily fluid
Full size table
Fig. 6
figure6

Top 10 most abundant miRNAs relative to all miRNA counts in each fluid. The 10 miRNAs that have the highest read counts in each fluid were illustrated relative to the total miRNA read counts of the fluid. Counts of the remaining miRNAs were summed up and illustrated as “other”

Full size image
Fig. 7
figure7

Relative abundance of the top 5 common miRNAs between the different fluids. Counts of the 5 common miRNAs are presented relative to the total miRNA counts of each fluid. The five common miRNAs represent a large fraction of the invasive fluids with the highest percentage in blood. They represent lower fractions in the non-invasive fluids

Full size image

Analysis of unique and novel miRNAs are valuable in evaluating the usefulness of using a specific specimen as a source of information. We ran the analysis of unique miRNAs in our study using 3 comparison groups: all fluids, invasive fluids and non-invasive fluids. The number of unique miRNAs detected in each comparison group are listed in Table 6. Blood, leukocytes and plasma had significantly higher numbers of unique miRNAs compared to the rest of fluids. Blood had the highest number of unique molecules (94 miRNAs), while plasma and leukocytes had 42 and 30 unique miRNAs, respectively. The comparison within the invasive group showed a similar trend. The comparison within the non-invasive group showed that cell-free saliva had the most unique miRNAs (32 miRNAs) followed by cell-free urine (11 miRNAs). Saliva and urine had minimal numbers of unique miRNAs compared to their cell-free preparations. The list of unique miRNAs from the three comparison groups are listed in three additional files (Additional file 7: Table S5, Additional file 8: Table S6, Additional file 9: Table S7).

Table 6 Number of unique miRNAs in each bodily fluid at a minimum of 5 counts
Full size table

All bodily fluids were analyzed for their novel miRNA candidates using miRDeep2. Signal-to-noise ratio of more than 10 was used to select for miRDeep2 score cutoff [68]. In bodily fluids where signal-to-noise ratio was less than 10 (saliva, urine and cell-free urine), we selected the score cutoff that corresponds to the highest signal-to-noise ratio (Table 7). Invasive bodily fluids had higher numbers of novel miRNA candidates than non-invasive fluids. The highest number of novel miRNA candidates was observed in plasma and blood, with 50 and 48 candidates, respectively. Serum had 20 novel candidates while leukocytes had 16 candidates. All the non-invasive fluids had 7 or less novel candidates. Sequences of novel miRNA candidates were matched to the miRCarta database (v1.0) of newly predicted human miRNAs [69]. More than 50% of blood and leukocytes novel miRNA candidates were present in miRCarta (66 and 56%, respectively). Plasma and serum had less miRNA candidates matching miRCarta database (30 and 45%, respectively). No miRCarta matches were found for the novel miRNA candidates of the non-invasive fluids. An additional file containing the list of novel miRNA candidates in each bodily fluid and their sequences, as well as the matching results to the miRCarta database is provided (Additional file 10: Table S8).

Table 7 Number of novel miRNA candidates in the different bodily fluids
Full size table

tRNA

Mapped tRNAs represented tDRs down to 18 nucleotides. The predominant tRNA fragments in all the bodily fluids was tRNAGly. This tRNA composed 86.5 and 87.6% of the total tRNA content in urine and cell-free urine, respectively. For the remaining bodily fluids, it made up 72.0 to 84.1% of the total tRNA content. The second most abundant tRNA was tRNAGlu, with a range of 6.7 to 21.4% of the tRNA content. Further analysis of tRNAs by tDRmapper to look at the exact tDRs composition showed that tRNAGly-GCC and tRNAGlu-CTC were the predominant fragments in all fluids. All samples, regardless of the fluid origin, shared similar tDRs composition. The quantification and coverage of the top 50 tDRs in each fluid are presented in Additional files 11 and 12. None of the remaining tRNAs in any of the fluids exceeded 3.1% of the tRNA content (Table 8). For blood, plasma, saliva and cell-free saliva, there was a higher diversity of tRNAs that represented 1% or more of the total tRNA content of the sample (5 tRNAs or more). However, for leukocytes, serum, urine and cell-free urine, there was a lower diversity (3 to 4 tRNAs).

Table 8 List of tRNAs that represents ≥1% of all tRNA counts in each bodily fluid
Full size table

Read counts of tRNAs were normalized using TMM at a CPM corresponding to a minimum of 5 counts in a library. The normalized reads were then used to generate principal component analysis which showed the separation of various bodily fluids based on their tRNA levels (Fig. 8). The analysis revealed closeness between saliva and cell-free saliva as well as closeness between the invasive fluids. However, urine and cell-free urine were dispersed between both saliva and serum. The Z scores of these tRNAs were used to generate a heatmap that indicates the levels of various tRNAs in the different samples (Fig. 9). Serum and the male urine/cell-free urine samples showed distant clustering from the rest of samples. The female urine/cell-free urine clustered with saliva/cell-free saliva. Blood, leukocytes and plasma showed similar clustering. The data shows clustering patterns based on sample biology and no difference between cell-depleted and non-depleted conditions.

Fig. 8
figure8

Principal component analysis of tRNAs in each bodily fluid. The plot was generated based on TMM-normalized tRNA counts. Samples of the same origin clustered closer to each other. However, urine samples are more dispersed from each other

Full size image
Fig. 9
figure9

Heatmap clustering of tRNAs in the various bodily fluids. The sex of sample donor is indicated as (F) for female donors or (M) for male donors. The analysis was generated using Z scores of TMM-normalized tRNA counts. The dendrogram shows distant clustering of urine samples (total and cell-free) based on sex. Male urine samples (total and cell-free) clustered close to serum while female urine samples (total and cell-free) are clustered with saliva

Full size image

piRNA

All bodily fluids had piR-016658 at different levels. The highest levels were seen in blood and serum (92.3 ± 1.8% and 94.0 ± 2.7%, respectively), followed by plasma and leukocytes (81.8 ± 33.1% and 73.5 ± 3.8%, respectively). It was the highest piRNA in cell-free saliva (40.3 ± 17.0%). It had lower concentrations, yet more than 10%, in saliva, urine and urine-cell free (14.9 ± 1.7%, 15.6 ± 12.9% and 14.5 ± 12.5, respectively). Urine and urine-cell free had piR-019825 as the highest piRNA (46.0 ± 40.4% and 58.7 ± 32.1%, respectively). Interestingly, piR-019825 was the second highest piRNA in plasma where it represented 15.2 ± 30.4% of the piRNA content (Table 9). An additional file contains the list of piRNAs at an average of 1% or more of the entire piRNAs counts of each bodily fluid (Additional file 13: Table S9).

Table 9 Predominant piRNAs in each bodily fluid
Full size table

Read counts of piRNAs were TMM-normalized at a CPM corresponding to a minimum of 5 counts in a library. The normalized counts were used to generate a PCA plot (Fig. 10). Close clustering was obtained between blood and leukocytes, saliva and cell-free saliva, and urine and cell-free urine. Plasma and serum showed relatively distant clustering from each other and from the other fluids. A similar clustering pattern was obtained from the heatmap. However, urine samples and most of the cell-free urine samples showed sex-dependent clustering (Fig. 11). Female urine samples clustered close to saliva samples, while male urine samples clustered close to plasma and serum samples. While clustering reflects the close biology of the samples, it has a distinct trend when compared to that of miRNA as it showed more overlap between invasive and non-invasive fluids.

Fig. 10
figure10

Principal component analysis of piRNAs in each bodily fluid. Analysis was generated based on TMM-normalized piRNA counts. Samples are clustered based on biology and fluids that share similar origin have close clustering. Close clustering is seen between the following fluid pairs: blood/leukocytes, saliva/cell-free saliva and urine/cell-free urine. Serum and plasma show distant clustering from the other fluids

Full size image
Fig. 11
figure11

Heatmap clustering of piRNAs in each bodily fluid. The sex of sample donor is indicated as (F) for female donors or (M) for male donors. Heatmap was generated using Z scores of TMM-normalized tRNA counts. The dendrogram shows distant clustering of urine samples (total and cell-free) based on sex. Male urine and cell-free urine samples clustered close to serum and plasma while female urine and most of cell-free urine samples are clustered with saliva

Full size image

Discussion

In this study, we investigated small RNA profiles in various bodily fluids using NGS in order to understand the distribution of the various biotypes between fluids as well as the molecular signature of each fluid. Purified total RNA from each fluid showed large variations in the RNA content and integrity. Saliva, cell-free saliva and blood have the highest RNA content. These elevated levels in both blood and saliva are due to their high number of cells. However, the high RNA levels in the cell-free saliva preparation were most likely due to the high bacterial content, as the cell-free preparation steps utilized in this study were aimed at removing mammalian but not bacterial cells. The lowest RNA yields were found in both total and cell-free urine preparations. This indicated very low cell content as well as minimal cell-free RNA content of urine samples collected from healthy individuals. Intra-fluid RNA yields are more consistent from fluids that have high cell content (blood, saliva and leukocytes). RNA integrity as measured by RIN value was also dependent on the cellular content of the fluid. This severely affected the RNA integrity of urine, cell-free urine, plasma and serum where they almost have no measurable integrity. In fact, while many small RNA-Seq library preparation pipelines require RNA with a minimum RIN value, this study showed that many bodily fluids of low cell content may not meet such requirements. However, in this study libraries were successfully constructed from RNA samples with low or no RIN value, suggesting that using RIN value as a sole determination of RNA quality may not be universally applicable.

For better representation of the actual library preparation workflow, we used a standard NGS preparation protocol based on equal input volumes from each bodily fluid preparation. An average of 12.57 ± 3.54 million reads were obtained from all samples, despite the large inter-fluid and intra-fluid variations in RNA yield and integrity. This indicates that regardless of sample type, a clean purification and robust library preparation can yield similar sequencing read outcomes. The critical parameter that would then define suitability of a sample to be used in small RNA profiling and discovery would be its actual biotype content. Read alignment to human genome varied between the different fluids based on their molecular composition. The lower percentage of leukocytes reads successfully aligned to the human genome is a result of their higher rRNA content. The lower percentage of saliva and cell-free saliva reads successfully aligned to the human genome is due to the high percentage of unmapped reads. We conducted exogenous mapping analysis on the saliva and cell-free saliva samples and 85–90% of the unmapped reads were mapped to bacteria (Additional file 14: Figure S3 and Additional file 15: Figure S4). This agrees with results of a recent study by Yeri et al. [70]. This in turn reduces the amount of valuable human RNA molecules that can be used in profiling or discovery. Efficient removal of salivary bacteria can be achieved by centrifugation [71]. However, removal of bacteria in this manner would reduce the number of reads directed towards bacterial sequences, thereby hindering the study of bacterial communities and/or pathogens that might be contained within these fluids.

While the profile of leukocytes showed a fairly even distribution of biotypes, all the other bodily fluids showed predominant reads from one or more biotypes. This can be of a significant value if these predominant biotypes are of known biological importance such as miRNAs, tRNAs or piRNAs. Blood has the highest levels of miRNAs (86.6 ± 12.3), which were about 3-fold the levels of leukocytes miRNAs (29.9 ± 3.3). A large portion of these blood miRNAs are lost from plasma and serum, as the separation and coagulation processes might be of the factors that affect miRNA distribution and recovery. Variations between plasma and serum miRNA content results from the stress during coagulation [18]. Non-invasive fluids had lower miRNA fractions, keeping in mind that they had a large percentage of unmapped reads (about 50% of saliva and 20% of urine reads). The relatively high bacterial content of saliva as well as the filtered and diluted nature of urine were key factors in this result. Recent analysis of urine, saliva and plasma miRNAs from NGS data showed lower miRNA counts from urine and saliva [70]. Profiling of miRNA in bodily fluids by RT-PCR in an earlier study showed similar low urine RNA concentrations and low numbers of detected miRNAs, while saliva had the highest number of miRNAs among the studied fluids [72]. The information obtained in our results could be used to guide methods for targeting specific biotypes in bodily fluids (via enrichment, separation, or depletion, for example).

The various bodily fluids have unique miRNA, tRNA and piRNA profiles that characterize the type and origin of the fluid as seen from the PCA plots and heatmaps. The different samples are well clustered based on miRNAs according to sample type and their biology. In addition, invasive and non-invasive fluids have distinct profiles and less variations between the fluids within the same group. PCA plots and heatmaps generated for tRNAs and piRNAs show a biology-related clustering, with overlap between invasive and non-invasive fluids. An interesting observation was the differential clustering of tRNAs and piRNAs from urine and cell-free urine samples based on the sex of the donor, indicating sex-related expression of these molecules. In addition, urine showed close clustering to serum, indicating that the latter might be the true liquid part of blood. The differentiation of the three biotypes (miRNA, tRNA and piRNA) in the different bodily fluids might be a result of fluid origin and biological functions. The higher impact of origin of fluid on miRNA distribution may refer to more specialized functions of miRNA in comparison to piRNAs and tRNAs, which might be involved in more general biological roles.

Non-invasive fluids had almost half the number of identified miRNAs. For urine, this may result from the filtering process by the kidneys. For saliva, the lower circulating nucleic acid content, relative to blood-related fluids may be the cause [73,74,75,76]. The larger miRNA fraction and lower numbers of unmapped reads for the invasive fluids explained their higher number of identified miRNAs. Almost every fluid had unique miRNAs that are specific only to that fluid, providing a specific signature for each fluid. The number of common and unique miRNAs between two fluids varied depending on biological relatedness. The invasive fluids collected in this study were more similar to the other invasive fluids, and the non-invasive fluids were more similar to the other non-invasive fluids collected. This may be due to the fact that the invasive fluids collected were all blood derived. The fluid-specific unique miRNAs can result from different cells secreting the different fluids. They may also result from the natural filtration process of some fluids, where some molecules are enriched, while others are depleted. None of the unique miRNAs were found among the top 20 most abundant molecules of each fluid. This indicated the need for higher read depth to detect miRNAs that might have specific functions. Five miRNAs: hsa-let-7a-5p, hsa-let-7f-5p, hsa-miR-191-5p, hsa-miR-26a-5p and hsa-miR-486-5p were common among the top 20 most abundant molecules in all fluids, indicating shared origin or function. These five abundant common miRNAs represented a large portion of the miRNA counts of invasive fluids (more than 40%), while they were relatively lower in non-invasive fluids (less than 30%). There was a set of 139 core miRNAs that are common among the different fluids and a set of 144 miRNAs that were shared between non-invasive fluids and blood. While the levels of these molecules may vary between fluid types, they might be promising biomarker candidates that can be detected from multiple sources, including non-invasive fluids.

An interesting observation was the variation between the most expressed miRNAs in the different fluids. In blood, plasma and serum, hsa-miR-486-5p was the most expressed, while hsa-miR-143-3p was the most expressed in saliva and cell-free saliva and hsa-miR-10b-5p was the predominant miRNA in urine and cell-free urine. We searched these 3 miRNAs on the human miRNA tissue atlas [77] to identify the tissue of origin. High quantile normalized expression levels of hsa-miR-486-5p were found in vein and muscle specimens, while hsa-miR-143-3p was highly elevated in esophagus and relatively high in colon, bladder and prostate specimens. The expression levels of hsa-miR-10b-5p were very high in the epididymis and elevated in kidney, colon and muscle specimens which explains the relatively higher expression levels of this miRNA in male urine samples. Recent studies indicate the importance of hsa-miR-486-5P as a cancer biomarker in non-small cell lung cancer [78], gastric cancer [79] and oral tongue squamous cell carcinoma [80]. It may act as a tumor suppressor miRNA [81] and may also be used to predict the efficacy of cancer vaccine treatment for colorectal cancer [82]. However, in many other cancer studies, this miRNA was not deregulated.

The large breadth of unique miRNAs found in blood, combined with an abundance of predicted novel miRNA candidates demonstrates the superiority of blood for miRNA profiling and discovery. However, blood has high levels of hsa-miR-486-5p, representing over 50% of its miRNA content. Other bodily fluids that had a relatively high miRNA content are plasma and, to some extent, serum. Plasma did not suffer from the presence of a predominant molecule as did blood and serum. This resulted in a high number of predicted novel miRNA candidates compared to all the other bodily fluids, making plasma a good alternative to blood. However, depletion of hsa-miR-486-5p from blood and serum could be a useful tool to direct a greater proportion of reads to other miRNA sequences.

Both saliva and urine did not offer the same advantage as the invasive fluids. They had lower miRNA content and this affected their molecular diversity. The most expressed miRNA in saliva and cell-free saliva was hsa-miR-143-3p (10–15%). It is also the second most expressed miRNA in saliva exosomes [83]. It is differentially expressed in senescence [84] and as a tumor suppressor in gliomas [85]. MicroRNA hsa-miR-10b-5p represented about 38–45% of urine and cell-free urine miRNAs. It has been recently reported to be the most expressed miRNA in urine samples [86]. hsa-miR-10b-5p plays a role in carcinoma metastasis and is overexpressed in colorectal cancer [87,88,89,90]. Due to its high expression, a lower proportion of reads will map to other miRNA sequences and its depletion should be considered as a priority for improving the diversity of miRNAs within urine specimens.

Novel miRNA prediction was not as efficient when dealing with non-invasive fluids. Their low signal-to-noise ratio made it hard to obtain accurate prediction. The only exception was the cell-free saliva, where a fair signal-to-noise ratio was achieved, and 7 novel miRNA candidates had been identified. It also had double the number of unique miRNAs compared to the other saliva and urine samples. Due to the removal of mammalian cells by centrifugation, cell-free saliva usually captures more circulating miRNAs than total saliva or the cellular fraction of saliva [20, 91]. This made the cell free saliva sample superior in terms of discovering unique and novel miRNAs candidates compared to the other non-invasive fluids. Matches were found between the miRCarta database of newly predicted human miRNAs and the predicted novel miRNA candidates of the invasive fluids (30 to 66%), while no matches were found between miRCarta and the non-invasive fluids. This might be due to the large number of miRNA studies from the invasive fluids as well as the higher counts and diversity of miRNAs from these fluids compared to the non-invasive fluids. This further indicates the higher potential of the invasive fluids in novel miRNA prediction. Although our findings are based on prediction of candidate miRNA and have not been validated by another technique, they showed the potential of these various fluids in novel miRNA discovery.

While tRNA fragments are a minor portion of blood and plasma small RNAs, they were well represented in serum and saliva preparations (39–46%) and were the major small RNA species of urine and cell-free urine (> 90%). The main component of these tRNAs in all the fluids was tRNAGly (72.0 to 87.6%), followed by tRNAGlu (6.7 to 21.4%). Urine samples, unlike the other fluids, had high sample-to-sample variations, with tRNAs ranging from 47 to 98% of small RNAs. Similar variations have been reported in a recent study on urine from ovarian cancer patients [86]. However, these variations may be correlated to the sex of the individual, where male urine has over 90% and female urine has about 70% or less. A larger study is needed to validate these findings. It is also interesting that the specific tRNA molecular composition of these tRNA fractions is consistent. Despite the overall fluctuations in urine tRNA fractions, changes in the molecular signature of tRNA molecules might still be valid for potential biomarker discovery. However, it may be limited by the lower urine tRNA molecular diversity, compared to blood, plasma and saliva. These variations in diversity were also observed between plasma, saliva and urine in a recent study [70]. However, the percentage abundance of molecules was different.

Plasma and leukocytes contain relatively high amounts of piRNA (8 and 5.8%, respectively). All the other bodily fluids contain less than 2%, which is related to the small fraction of piRNAs that are consistently being expressed in normal and cancer cells [92]. It is interesting that a single piRNA molecule, hsa-piR-016658, was the most expressed in all bodily fluids except in saliva, urine and cell-free urine, where it was the second most abundant. This molecule is associated with patients with prostate cancer [93]. The most abundant piRNA molecule in urine and cell-free urine was hsa-piR-019825, which is deregulated in colorectal cancer patients [93]. Given the high number of human piRNAs, they may play a role as an important small RNA species with functional targets that are yet to be elucidated and correlated with various disease conditions. Recent studies have identified differentially expressed piRNA molecules as potential biomarkers of various cancers [94,95,96,97]. The relatively high levels of these molecules in plasma might prove important as potential biomarkers. The low piRNAs levels in the other fluids can be overcome by size selection methods, albeit not easily as they overlap with other small RNA species.

Plasma and serum had a large fraction of reads mapping to miscellaneous RNA (misc_RNA) (58 and 35%, respectively). Only 4 YRNA-derived small RNAs (s-RNYs) sequences were elevated within the misc_RNA fractions of these two fluids: RNY4, RNY4P10, RNY4P7 and YRNA.295. It has been previously reported that s-RNYs are abundant in human serum and plasma [98]. They are potential cancer biomarkers and regulators of inflammation and cell death [99, 100].

Conclusions

Our study showed that it is possible to successfully accomplish NGS of the different bodily fluids (blood, plasma, serum, saliva, cell-free saliva, urine and cell-free urine), even with the high variations in the volumes used for RNA purification as well as the high variations in concentrations of the isolated RNA. Despite the ease of collection and handling of non-invasive fluids, they did not provide the same small RNA diversity and sample consistency as invasive fluids. However, this study showed that these samples can still be routinely profiled. Furthermore, the signatures of these non-invasive fluids are very likely linked to their origin. For example, urine may be a good candidate for studying diseases related to organs such as kidney and bladder, although careful result interpretation should be considered when investigating male and female urine, as their biotypes may be sex-dependent. This observation is limited by the sample size of our study and is yet to be investigated on a large sample size study. An organ and fluid small RNA index might be needed to track and correlate origins and functions of the various molecules. Processing of larger volumes of urine, and bacterial removal from saliva preparations might improve their NGS mapping to human targets. In addition, depletion of specific molecules or selection/enrichment of target molecules from almost every bodily fluid may significantly increase flow cell capacity for target molecules and in turn provide a meaningful read depth. Successful clustering of bodily fluids based on their miRNA distribution can be expanded to cohorts that can be differentiated according to their miRNA, and possibly in combination with other small RNAs. Therefore, a biomarker within these fluids would be the overall biotype distribution and the molecular signature within these biotypes, rather than a single molecule.

Abbreviations

CPM:

Counts per million

lincRNA:

Long intergenic noncoding RNA

miRNA:

MicroRNA

misc_RNA:

Miscellaneous RNA

mRNA:

Messenger RNA

Mt_rRNA:

Mitochondrial rRNA

Mt_tRNA:

Mitochondrial tRNA

ncRNA:

Non-coding RNA

NGS:

Next generation sequencing

PCA:

Principal component analysis

piRNA:

Piwi-interacting RNA

RIN:

RNA integrity number

rRNA:

Ribosomal RNA

snoRNA:

Small nucleolar RNA

snRNA:

Small nuclear RNA

s-RNYs:

Y RNA-derived small RNAs

TMM:

Trimmed mean of M-values

References

  1. 1.

    Storz G. An expanding universe of noncoding RNAs. Science. 2002;296:1260–1263. Available from: http://www.sciencemag.org/cgi/doi/10.1126/science.1072249

  2. 2.

    Gambari R, Brognara E, Spandidos DA, Fabbri E. Targeting oncomiRNAs and mimicking tumor suppressor miRNAs: Ew trends in the development of miRNA therapeutic strategies in oncology (review). Int. J. Oncol. [internet]. 2016;49:5–32. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27175518

  3. 3.

    Keller A, Meese E. Can circulating miRNAs live up to the promise of being minimal invasive biomarkers in clinical settings? Wiley Interdiscip. Rev. RNA [internet]. 2016;7:148–156. Available from: http://www.ncbi.nlm.nih.gov/pubmed/26670867

  4. 4.

    Cristodero M, Polacek N. The multifaceted regulatory potential of tRNA-derived fragments. Non-coding RNA Investig. [Internet]. 2017;7–7. Available from: http://ncri.amegroups.com/article/view/3820/4459

  5. 5.

    Guay C, Regazzi R. Circulating microRNAs as novel biomarkers for diabetes mellitus. Nat Rev Endocrinol [Internet] 2013;9:513–521. Available from: http://www.nature.com/doifinder/10.1038/nrendo.2013.86

  6. 6.

    Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell [internet]. 2004;116:281–297. Available from: http://www.ncbi.nlm.nih.gov/pubmed/14744438

  7. 7.

    Kim VN, Han J, Siomi MC. Biogenesis of small RNAs in animals. Nat Rev Mol Cell Biol [Internet] 2009;10:126–139. Available from: http://www.nature.com/doifinder/10.1038/nrm2632

  8. 8.

    Ortiz-Quintero B. Cell-free microRNAs in blood and other body fluids, as cancer biomarkers. Cell Prolif. [internet]. 2016;49:281–303. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27218664

  9. 9.

    Brennecke J, Hipfner DR, Stark A, Russell RB, Cohen SM. Bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell [internet]. 2003;113:25–36. Available from: http://www.ncbi.nlm.nih.gov/pubmed/12679032

  10. 10.

    Scapoli L, Palmieri A, Lo Muzio L, Pezzetti F, Rubini C, Girardi A, et al. MicroRNA expression profiling of oral carcinoma identifies new markers of tumor progression. Int. J. Immunopathol. Pharmacol. 2010;23:1229–34. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21244772

    Article  PubMed  CAS  Google Scholar 

  11. 11.

    Liang Z, Bian X, Shim H. Downregulation of microRNA-206 promotes invasion and angiogenesis of triple negative breast cancer. Biochem. Biophys. Res. Commun. 2016;477:461–6. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27318091

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. 12.

    Friedman RC, Farh KKH, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009;19:92–105. Available from: http://www.ncbi.nlm.nih.gov/pubmed/18955434

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. 13.

    miRBase: the microRNA database [Internet]. [cited 2018 Mar 7]. Available from: http://www.mirbase.org/cgi-bin/browse.pl?org=hsa

  14. 14.

    Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. [cited 2018 Mar 7]; Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3064324/pdf/pnas.201019055.pdf

  15. 15.

    Hunter MP, Ismail N, Zhang X, Aguda BD, Lee EJ, Yu L, et al. Detection of microRNA expression in human peripheral blood microvesicles. PLoS One. 2008 [cited 2018 Mar 7];3. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2577891/pdf/pone.0003694.pdf

  16. 16.

    Valadi H, Ekström K, Bossios A, Sjöstrand M, Lee JJ, Lötvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat. Cell Biol; 2007 [cited 2018 Mar 7];9:654–659. Available from: http://www.nature.com/articles/ncb1596. Nature Publishing Group

  17. 17.

    Vickers KC, Remaley AT. Lipid-based carriers of microRNAs and intercellular communication. [cited 2018 Mar 7]; Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5570485/pdf/nihms892052.pdf

  18. 18.

    Wang K, Zhang S, Weber J, Baxter D, Galas DJ. Export of microRNAs and microRNA-protective protein by mammalian cells. Nucleic Acids Res. 2010;38:7248–59. Available from: http://www.ncbi.nlm.nih.gov/pubmed/20615901

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. 19.

    Dumache R, Ciocan V, Muresan C, Rogobete AF, Enache A. Circulating microRNAs as promising biomarkers in forensic body fluids identification. Clin. Lab. 2015;61:1129–35. Available from: http://www.ncbi.nlm.nih.gov/pubmed/26554231

    PubMed  CAS  Google Scholar 

  20. 20.

    Park NJ, Zhou H, Elashoff D, Henson BS, Kastratovic DA, Abemayor E, et al. Salivary microRNA: Discovery, characterization, and clinical utility for oral cancer detection. Clin. Cancer Res. 2009;15:5473–7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19706812

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. 21.

    Tölle A, Jung M, Rabenhorst S, Kilic E, Jung K, Weikert S. Identification of microRNAs in blood and urine as tumour markers for the detection of urinary bladder cancer. Oncol. Rep. 2013;30:1949–56. Available from: http://www.ncbi.nlm.nih.gov/pubmed/23877086

    Article  PubMed  CAS  Google Scholar 

  22. 22.

    Suryawanshi S, Vlad AM, Lin HM, Mantia-Smaldone G, Laskey R, Lee M, et al. Plasma MicroRNAs as novel biomarkers for endometriosis and endometriosis-associated ovarian cancer. Clin. Cancer Res. 2013;19:1213–24. Available from: http://www.ncbi.nlm.nih.gov/pubmed/23362326

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. 23.

    Hu Z, Chen X, Zhao Y, Tian T, Jin G, Shu Y, et al. Serum microRNA signatures identified in a genome-wide serum microRNA expression profiling predict survival of non-small-cell lung cancer. J. Clin. Oncol. 2010;28:1721–6. Available from: http://www.ncbi.nlm.nih.gov/pubmed/20194856

    Article  PubMed  Google Scholar 

  24. 24.

    Gilad S, Meiri E, Yogev Y, Benjamin S, Lebanony D, Yerushalmi N, et al. Serum microRNAs are promising novel biomarkers. PLoS One. 2008;3:e3148. Available from: http://www.ncbi.nlm.nih.gov/pubmed/18773077

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. 25.

    Xie Z, Chen G, Zhang X, Li D, Huang J, Yang C, et al. Salivary MicroRNAs as Promising Biomarkers for Detection of Esophageal Cancer. Lo AWI, editor. PLoS One. 2013 [cited 2016 Aug 10];8:e57502. Available from: http://dx.plos.org/10.1371/journal.pone.0057502. Public Library of Science

  26. 26.

    Kosaka N, Iguchi H, Ochiya T. Circulating microRNA in body fluid: a new potential biomarker for cancer diagnosis and prognosis. Cancer Sci. 2010;101:2087–92. Available from: http://www.ncbi.nlm.nih.gov/pubmed/20624164

    Article  PubMed  CAS  Google Scholar 

  27. 27.

    Zhao N, Jin L, Fei G, Zheng Z, Zhong C. Serum microRNA-133b is associated with low ceruloplasmin levels in Parkinson’s disease. Park. Relat. Disord, Available from. 2014;20:1177–80. http://www.ncbi.nlm.nih.gov/pubmed/25218846

  28. 28.

    Kirchner S, Ignatova Z. Emerging roles of tRNA in adaptive translation, signalling dynamics and disease. Nat. Rev. Genet. 2015;16:98–112. Available from: http://www.nature.com/doifinder/10.1038/nrg3861

    Article  PubMed  CAS  Google Scholar 

  29. 29.

    Selitsky SR, Baran-Gale J, Honda M, Yamane D, Masaki T, Fannin EE, et al. Small tRNA-derived RNAs are increased and more abundant than microRNAs in chronic hepatitis B and C. Sci. Rep. 2015;5:7675. Available from: http://www.nature.com/articles/srep07675

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. 30.

    Ivanov P, Emara MM, Villen J, Gygi SP, Anderson P. Angiogenin-induced tRNA fragments inhibit translation initiation. Mol. Cell. 2011;43:613–23. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21855800

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. 31.

    Yamasaki S, Ivanov P, Hu GF, Anderson P. Angiogenin cleaves tRNA and promotes stress-induced translational repression. J. Cell Biol. 2009;185:35–42. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19332886

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. 32.

    Saikia M, Jobava R, Parisien M, Putnam A, Krokowski D, Gao X-H, et al. Angiogenin-Cleaved tRNA Halves Interact with Cytochrome c, Protecting Cells from Apoptosis during Osmotic Stress. Mol. Cell. Biol. 2014;34:2450–2463. Available from: http://mcb.asm.org/cgi/doi/10.1128/MCB.00136-14

  33. 33.

    Schaffer AE, Eggens VRC, Caglayan AO, Reuter MS, Scott E, Coufal NG, et al. CLP1 founder mutation links tRNA splicing and maturation to cerebellar development and neurodegeneration. Cell. 2014;157:651–63. Available from: http://www.ncbi.nlm.nih.gov/pubmed/24766810

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. 34.

    Gebetsberger J, Zywicki M, Künzi A, Polacek N. TRNA-derived fragments target the ribosome and function as regulatory non-coding RNA in Haloferax volcanii. Archaea. 2012;2012:260909. Available from: http://www.ncbi.nlm.nih.gov/pubmed/23326205

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. 35.

    Gebetsberger J, Wyss L, Mleczko AM, Reuther J, Polacek N. A tRNA-derived fragment competes with mRNA for ribosome binding and regulates translation during stress. RNA Biol. 2017;14:1364–73. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27892771

    Article  PubMed  Google Scholar 

  36. 36.

    Sharma U, Conine CC, Shea JM, Boskovic A, Derr AG, Bing XY, et al. Biogenesis and function of tRNA fragments during sperm maturation and fertilization in mammals. Science. 2016;351:391–396. Available from: http://www.sciencemag.org/cgi/doi/10.1126/science.aad6780

  37. 37.

    Venkatesh T, Suresh PS, Tsutsumi R. TRFs: miRNAs in disguise. Gene. 2016;579:133–8. Available from: http://www.ncbi.nlm.nih.gov/pubmed/26743126

    Article  PubMed  CAS  Google Scholar 

  38. 38.

    Haussecker D, Huang Y, Lau A, Parameswaran P, Fire AZ, Kay MA. Human tRNA-derived small RNAs in the global regulation of RNA silencing. Rna. 2010;16:673–695. Available from: http://rnajournal.cshlp.org/cgi/doi/10.1261/rna.2000810

  39. 39.

    Elbarbary RA, Takaku H, Uchiumi N, Tamiya H, Abe M, Takahashi M, et al. Modulation of gene expression by human cytosolic tRNase Z(L) through 5′-half-tRNA. PLoS One. 2009;4:e5908. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19526060

  40. 40.

    Ghildiyal M, Zamore PD. Small silencing RNAs: an expanding universe. Nat Rev Genet [Internet] 2009;10:94–108. Available from: http://www.nature.com/doifinder/10.1038/nrg2504

  41. 41.

    Kumar A, Karmarkar AM, Tan A, Graham JE, Arcari CM, Ottenbacher KJ, et al. The effect of obesity on incidence of disability and mortality in Mexicans aged 50 years and older. Salud Publica Mex. 2015 [cited 2017 Jul 12];57:S31–S38. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4572697/pdf/nihms673180.pdf

  42. 42.

    Green D, Fraser WD, Dalmay T. Transfer RNA-derived small RNAs in the cancer transcriptome. Pflugers Arch. Eur. J. Physiol. 2016;468:1041–7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27095039

    Article  CAS  Google Scholar 

  43. 43.

    Garcia-Silva MR, Cabrera-Cabrera F, Güida MC, Cayota A. Hints of tRNA-derived small RNAs role in RNA silencing mechanisms. Genes (Basel), Available from. 2012;3:603–14. http://www.ncbi.nlm.nih.gov/pubmed/24705078

  44. 44.

    Goodarzi H, Liu X, Nguyen HCB, Zhang S, Fish L, Tavazoie SF. Endogenous tRNA-derived fragments suppress breast cancer progression via YBX1 displacement. Cell. 2015;161:790–802. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25957686

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. 45.

    Maute RL, Schneider C, Sumazin P, Holmes A, Califano A, Basso K, et al. tRNA-derived microRNA modulates proliferation and the DNA damage response and is down-regulated in B cell lymphoma. Proc. Natl. Acad. Sci. 2013;110:1404–1409. Available from: http://www.pnas.org/lookup/doi/10.1073/pnas.1206761110

  46. 46.

    Atala A. Re: sex hormone-dependent tRNA halves enhance cell proliferation in breast and prostate cancers. J. Urol. 2016;195:1168–9. Available from: http://www.ncbi.nlm.nih.gov/pubmed/26124144

    CAS  PubMed  Google Scholar 

  47. 47.

    Balatti V, Pekarsky Y, Croce CM. Role of the tRNA-derived small RNAs in Cancer: new potential biomarkers and target for therapy [internet]. 1st ed. Elsevier Inc.; 2017. Available from: https://doi.org/10.1016/bs.acr.2017.06.007. Adv. Cancer Res

  48. 48.

    Hatfield SD, Shcherbata HR, Fischer KA, Nakahara K, Carthew RW, Ruohola-Baker H. Stem cell division is regulated by the microRNA pathway. Nature. 2005;435:974–978. Available from: http://www.nature.com/doifinder/10.1038/nature03816

  49. 49.

    Rouget C, Papin C, Boureux A, Meunier AC, Franco B, Robine N, et al. Maternal mRNA deadenylation and decay by the piRNA pathway in the early Drosophila embryo. Nature. 2010;467:1128–1132. Available from: http://www.nature.com/doifinder/10.1038/nature09465

  50. 50.

    piRNABank: : a web resource on classified and clustered Piwi-interacting RNAs [Internet]. [cited 2018 Mar 7]. Available from: http://pirnabank.ibab.ac.in/stats.html

  51. 51.

    Aravin AA, Sachidanandam R, Bourc’his D, Schaefer C, Pezic D, Toth KF, et al. A piRNA Pathway Primed by Individual Transposons Is Linked to De Novo DNA Methylation in Mice. Mol. Cell. 2008;31:785–99. Available from: http://www.ncbi.nlm.nih.gov/pubmed/18922463

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. 52.

    Aravin AA, Bourc’his D. Small RNA guides for de novo DNA methylation in mammalian germ cells. Genes Dev. 2008;22:970–5. Available from: http://www.ncbi.nlm.nih.gov/pubmed/18413711

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. 53.

    Hirakata S, Siomi MC. piRNA biogenesis in the germline: From transcription of piRNA genomic sources to piRNA maturation. Biochim. Biophys. Acta - Gene Regul. Mech. 2016;1859:82–92. Available from: https://doi.org/10.1016/j.bbagrm.2015.09.002. Elsevier B.V

  54. 54.

    Esteller M. Non-coding RNAs in human disease. Nat. Rev. Genet. [Internet]. 2011;12:861–874. Available from: http://www.nature.com/doifinder/10.1038/nrg3074

  55. 55.

    Zhang J, Chiodini R, Badr A, Zhang G. The impact of next-generation sequencing on genomics [Internet]. J. Genet. Genomics. 2011 [cited 2017 Jul 16]. p. 95–109. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3076108/pdf/nihms-282401.pdf

  56. 56.

    Byron SA, Van Keuren-Jensen KR, Engelthaler DM, Carpten JD, Craig DW. Translating RNA sequencing into clinical diagnostics: Opportunities and challenges. Nat. Rev. Genet. 2016 [cited 2017 Jul 16];17:257–271. Available from: http://www.nature.com/doifinder/10.1038/nrg.2016.10

  57. 57.

    Shore S, Henderson JM, Lebedev A, Salcedo MP, Zon G, McCaffrey AP, et al. Small RNA library preparation method for next-generation sequencing using chemical modifications to prevent adapter dimer formation. PLoS One. 2016 [cited 2017 Jul 16];11. Available from: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0167009&type=printable

  58. 58.

    Chen X, Ba Y, Ma L, Cai X, Yin Y, Wang K, et al. Characterization of microRNAs in serum: A novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008;18:997–1006. Available from: http://www.nature.com/doifinder/10.1038/cr.2008.282

  59. 59.

    Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc. Natl. Acad. Sci. 2008;105:10513–10518. Available from: http://www.pnas.org/cgi/doi/10.1073/pnas.0804549105

  60. 60.

    Wang J, Zhang KY, Liu SM, Sen S. Tumor-associated circulating micrornas as biomarkers of cancer. Molecules. 2014;19:1912–38. Available from: http://www.ncbi.nlm.nih.gov/pubmed/24518808

    Article  PubMed  CAS  Google Scholar 

  61. 61.

    Cheng J, Guo JM, Xiao BX, Miao Y, Jiang Z, Zhou H, et al. PiRNA, the new non-coding RNA, is aberrantly expressed in human cancer cells. Clin. Chim. Acta. 2011;412:1621–5. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21616063

    Article  PubMed  CAS  Google Scholar 

  62. 62.

    Hashim A, Rizzo F, Marchese G, Ravo M, Tarallo R, Nassa G, et al. RNA sequencing identifies specific PIWI-interacting small non-coding RNA expression patterns in breast cancer. Oncotarget. 2014;5:9901–10. Available from: http://www.oncotarget.com/fulltext/2476

    Article  PubMed  PubMed Central  Google Scholar 

  63. 63.

    Li Y, Wu X. Piwi-Interacting RNAs (piRNAs) Are Dysregulated in Renal Cell Carcinoma and Associated with Tumor Metastasis and Cancer-Specific Survival. Mol. Med. 2015;21:1. Available from: http://www.molmed.org/content/pdfstore/14_203_Li.pdf

    Google Scholar 

  64. 64.

    Reeves ME, Firek M, Jliedi A, Amaar YG. Identification and characterization of RASSF1C piRNA target genes in lung cancer cells. Oncotarget. 2017; Available from: http://www.oncotarget.com/abstract/15965

  65. 65.

    Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005 [cited 2017 Jan 1];33:e179. Available from: https://academic.oup.com/nar/article-lookup/doi/10.1093/nar/gni178

  66. 66.

    Subramanian SL, Kitchen RR, Alexander R, Carter BS, Cheung KH, Laurent LC, et al. Integration of extracellular RNA profiling data using metadata, biomedical ontologies and Linked Data technologies. J. Extracell. Vesicles. 2015;4:27497. Available from: http://www.ncbi.nlm.nih.gov/pubmed/26320941

    Article  PubMed  CAS  Google Scholar 

  67. 67.

    Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11:R25. Available from: http://genomebiology.biomedcentral.com/articles/10.1186/gb-2010-11-3-r25

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  68. 68.

    Dhahbi JM, Atamna H, Boffelli D, Magis W, Spindler SR, Martin DIK. Deep sequencing reveals novel micrornas and regulation of microRNA expression during cell senescence. PLoS One. 2011;6:e20509. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21637828

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  69. 69.

    Backes C, Fehlmann T, Kern F, Kehl T, Lenhof HP, Meese E, et al. MiRCarta: A central repository for collecting miRNA candidates. Nucleic Acids Res. 2018 [cited 2018 Apr 12];46:D160–D167. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5753177/pdf/gkx851.pdf

  70. 70.

    Yeri A, Courtright A, Reiman R, Carlson E, Beecroft T, Janss A, et al. Total extracellular small RNA profiles from plasma, saliva, and urine of healthy subjects. Sci. Rep. 2017 [cited 2017 Jul 11];7:44061. Available from: http://www.nature.com/articles/srep44061

  71. 71.

    Park NJ, Zhou X, Yu T, Brinkman BMN, Zimmermann BG, Palanisamy V, et al. Characterization of salivary RNA by cDNA library analysis. Arch. Oral Biol. 2007 [cited 2017 Oct 11];52:30–35. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2743855/pdf/nihms15843.pdf

  72. 72.

    Weber JA, Baxter DH, Zhang S, Huang DY, Huang KH, Lee MJ, et al. The microRNA spectrum in 12 body fluids. Clin Chem. 2010;56:1733–41.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  73. 73.

    Gallo A, Tandon M, Alevizos I, Illei GG. The majority of microRNAs detectable in serum and saliva is concentrated in exosomes. Afarinkia K, editor. PLoS One. 2012 [cited 2016 Aug 10];7:e30679. Available from: http://dx.plos.org/10.1371/journal.pone.0030679. Public Library of Science

  74. 74.

    Majem B, Rigau M, Reventós J, Wong DT. Non-coding RNAs in saliva: emerging biomarkers for molecular diagnostics. Int. J. Mol. Sci. 2015;16:8676–98. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25898412

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  75. 75.

    Spielmann N, Ilsley D, Gu J, Lea K, Brockman J, Heater S, et al. The human salivary RNA transcriptome revealed by massively parallel sequencing. Clin. Chem. 2012;58:1314–21. Available from: http://www.ncbi.nlm.nih.gov/pubmed/22773539

    Article  PubMed  CAS  Google Scholar 

  76. 76.

    Li M, Zeringer E, Barta T, Schageman J, Cheng A, Vlassov A V. Analysis of the RNA content of the exosomes derived from blood serum and urine and its potential as biomarkers. Philos. Trans. R. Soc. B Biol. Sci. 2014;369:20130502. Available from: http://rstb.royalsocietypublishing.org/cgi/doi/10.1098/rstb.2013.0502

  77. 77.

    Ludwig N, Leidinger P, Becker K, Backes C, Fehlmann T, Pallasch C, et al. Distribution of miRNA expression across human tissues. Nucleic Acids Res. 2016 [cited 2018 Apr 13];44:3865–3877. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4856985/pdf/gkw116.pdf

  78. 78.

    Sromek M, Glogowski M, Chechlinska M, Kulinczak M, Szafron L, Zakrzewska K, et al. Changes in plasma miR-9, miR-16, miR-205 and miR-486 levels after non-small cell lung cancer resection. Cell. Oncol. 2017;40:529–36. Available from: http://www.ncbi.nlm.nih.gov/pubmed/28634901

    Article  CAS  Google Scholar 

  79. 79.

    Sierzega M, Kaczor M, Kolodziejczyk P, Kulig J, Sanak M, Richter P. Evaluation of serum microRNA biomarkers for gastric cancer based on blood and tissue pools profiling: The importance of MIR-21 and MIR-331. Br. J. Cancer. 2017;117:266–73. Available from: http://www.nature.com/doifinder/10.1038/bjc.2017.190

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  80. 80.

    Chen Z, Yu T, Cabay RJ, Jin Y, Mahjabeen I, Luan X, et al. miR-486-3p, miR-139-5p, and miR-21 as Biomarkers for the Detection of Oral Tongue Squamous Cell Carcinoma. Biomark. Cancer [Internet]. 2017;9:1–8. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5224348/

  81. 81.

    Ye H, Yu X, Xia J, Tang X, Tang L, Chen F. MiR-486-3p targeting ECM1 represses cell proliferation and metastasis in cervical cancer. Biomed. Pharmacother. [Internet], Available from. 2016;80:109–14. http://www.ncbi.nlm.nih.gov/pubmed/27133046

  82. 82.

    Shindo Y, Hazama S, Nakamura Y, Inoue Y, Kanekiyo S, Suzuki N, et al. miR-196b, miR-378a and miR-486 are predictive biomarkers for the efficacy of vaccine treatment in colorectal cancer. Oncol. Lett. [Internet]. 2017;14:1355–62. Available from: http://www.ncbi.nlm.nih.gov/pubmed/28789351

    Article  Google Scholar 

  83. 83.

    Ogawa Y, Taketomi Y, Murakami M, Tsujimoto M, Yanoshita R. Small RNA transcriptomes of two types of exosomes in human whole saliva determined by next generation sequencing. Biol Pharm Bull [Internet]. 2013;36:66–75. Available from: http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=23302638

    Article  CAS  Google Scholar 

  84. 84.

    Bonifacio LN, Jarstfer MB. MiRNA profile associated with replicative senescence, extended cell culture, and ectopic telomerase expression in human foreskin fibroblasts. PLoS One [Internet]. 2010;5:1–8. Available from: http://www.ncbi.nlm.nih.gov/pubmed/20824140

    Google Scholar 

  85. 85.

    Wang L, Shi Z, Jiang C, Liu X, Chen Q, Qian X, et al. MiR-143 acts as a tumor suppressor by targeting N-RAS and enhances temozolomide-induced apoptosis in glioma. Oncotarget [Internet]. 2014;5:5416–27. Available from: http://www.oncotarget.com/fulltext/2116

    Google Scholar 

  86. 86.

    Zhou K, Spillman MA, Behbakht K, Komatsu JM, Abrahante JE, Hicks D, et al. A method for extracting and characterizing RNA from urine: For downstream PCR and RNAseq analysis. Anal. Biochem. [Internet]. 2017;536:8–15. Available from: https://doi.org/10.1016/j.ab.2017.08.003. Elsevier Inc

  87. 87.

    Zhang L, Sun J, Wang B, Ren JC, Su W, Zhang T. MicroRNA-10b triggers the epithelial–mesenchymal transition (EMT) of laryngeal carcinoma Hep-2 cells by directly targeting the E-cadherin. Appl. Biochem. Biotechnol. [internet]. 2015;176:33–44. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25875782

    Article  CAS  Google Scholar 

  88. 88.

    Xiao H, Li H, Yu G, Xiao W, Hu J, Tang K, et al. MicroRNA-10b promotes migration and invasion through KLF4 and HOXD10 in human bladder cancer. Oncol. Rep. [Internet]. 2014;31:1832–8. Available from: http://www.ncbi.nlm.nih.gov/pubmed/26311318

    Article  CAS  Google Scholar 

  89. 89.

    Ma Z, Chen Y, Min L, Li L, Huang H, Li J, et al. Augmented miR-10b expression associated with depressed expression of its target gene KLF4 involved in gastric carcinoma. Int. J. Clin. Exp. Pathol. [Internet]. 2015;8:5071–9. Available from: http://www.ncbi.nlm.nih.gov/pubmed/26191201

    CAS  Google Scholar 

  90. 90.

    Abdelmaksoud-Dammak R, Chamtouri N, Triki M, Saadallah-Kallel A, Ayadi W, Charfi S, et al. Overexpression of miR-10b in colorectal cancer patients: Correlation with TWIST-1 and E-cadherin expression. Tumor Biol. [Internet]. 2017 [cited 2017 Oct 20];39:101042831769591. Available from: http://journals.sagepub.com/doi/10.1177/1010428317695916. SAGE PublicationsSage UK: London, England

  91. 91.

    Lin X, Lo H-C, Wong DTW, Xiao X. Noncoding RNAs in human saliva as potential disease biomarkers. Front Genet [Internet] 2015;6:1–6. Available from: http://www.frontiersin.org/RNA/10.3389/fgene.2015.00175/full

  92. 92.

    Martinez VD, Vucic EA, Thu KL, Hubaux R, Enfield KSS, Pikor LA, et al. Unique somatic and malignant expression patterns implicate PIWI-interacting RNAs in cancer-type specific biology. Sci. Rep. [Internet]. 2015 [cited 2017 Oct 10];5. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4444957/pdf/srep10423.pdf

  93. 93.

    Yuan T, Huang X, Woodcock M, Du M, Dittmar R, Wang Y, et al. Plasma extracellular RNA profiles in healthy and cancer patients. Sci. Rep. [Internet]. 2016;6:19413. Available from: http://www.nature.com/articles/srep19413

    Article  CAS  Google Scholar 

  94. 94.

    Lim SL, Ricciardelli C, Oehler MK, De Arao Tan IMD, Russell D, Grützner F. Overexpression of piRNA pathway genes in epithelial ovarian cancer. PLoS One [Internet]. 2014 [cited 2017 Oct 10];9. Available from: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0099687&type=printable

  95. 95.

    Müller S, Raulefs S, Bruns P, Afonso-Grunz F, Plötner A, Thermann R, et al. Next-generation sequencing reveals novel differentially regulated mRNAs, lncRNAs, miRNAs, sdRNAs and a piRNA in pancreatic cancer. Mol. Cancer [Internet]. 2015;14:94. Available from: http://molecular-cancer.biomedcentral.com/articles/10.1186/s12943-015-0358-5

  96. 96.

    Martinez VD, Enfield KSS, Rowbotham DA, Lam WL. An atlas of gastric PIWI-interacting RNA transcriptomes and their utility for identifying signatures of gastric cancer recurrence. Gastric Cancer [internet]. 2016;19:660–5. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25779424

    Article  CAS  Google Scholar 

  97. 97.

    Koduru S V, Tiwari AK, Hazard SW, Mahajan M, Ravnic DJ. Exploration of small RNA-seq data for small non-coding RNAs in Human Colorectal Cancer. J. Genomics [Internet]. 2017 [cited 2017 Oct 10];5:16–31. Available from: http://www.jgenomics.com/v05p0016.htm

  98. 98.

    Dhahbi JM, Spindler SR, Atamna H, Boffelli D, Mote P, Martin DIK. 5’-YRNA fragments derived by processing of transcripts from specific YRNA genes and pseudogenes are abundant in human serum and plasma. Physiol. Genomics [Internet]. 2013;45:990–8. Available from: http://physiolgenomics.physiology.org/cgi/doi/10.1152/physiolgenomics.00129.2013

    Article  CAS  Google Scholar 

  99. 99.

    Dhahbi, Spinder S, Atamna H, Boffelli D, Martin D. Deep Sequencing of Serum Small RNAs Identifies Patterns of 5&amp;#39; tRNA Half and YRNA Fragment Expression Associated with Breast Cancer. Biomark. Cancer [Internet]. 2014 [cited 2016 Apr 11];6:37. Available from: http://www.la-press.com/deep-sequencing-of-serum-small-rnas-identifies-patterns-of-5-trna-half-article-a4553

  100. 100.

    Hizir Z, Bottini S, Grandjean V, Trabucchi M, Repetto E. RNY (YRNA)-derived small RNAs regulate cell death and inflammation in monocytes/macrophages. Cell Death Dis. [Internet]. 2017 [cited 2017 Oct 11];8:e2530. Available from: http://www.nature.com/doifinder/10.1038/cddis.2016.429

Download references

Availability of data and materials

The raw sequencing files of the study are available on the NCBI Sequence Read Archive (SRA) through the following link: https://www.ncbi.nlm.nih.gov/sra/SRP136264, SRA accession: SRP136264.

Author information

Affiliations

  1. Norgen Biotek Corp, Thorold, ON, L2V 4Y6, Canada

    Mohamed El-Mogy, Bernard Lam, Taha A. Haj-Ahmad, Darrick Yu, Lucas Nosal, Nezar Rghei, Pam Roberts & Yousef Haj-Ahmad

  2. Molecular Biology Department, National Research Centre, Dokki, Giza, Egypt

    Mohamed El-Mogy

  3. Department of Biological Sciences, Brock University, St. Catharines, ON, L2S 3A1, Canada

    Shannon McGowan & Yousef Haj-Ahmad

Authors
  1. Mohamed El-Mogy
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Bernard Lam
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Taha A. Haj-Ahmad
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Shannon McGowan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Darrick Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Lucas Nosal
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Nezar Rghei
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Pam Roberts
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Yousef Haj-Ahmad
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

ME, BL, NR and YHA conceived, designed and initiated the project. ME, BL, SM performed experiments. THA and DY contributed to the analysis tools. THA, LN and ME carried out data analysis. ME and PR drafted the manuscript. BL, THA, NR and YHA helped to revise the manuscript. LN re-edited the language of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mohamed El-Mogy.

Ethics declarations

Ethics approval and consent to participate

The study protocol and consent were reviewed and approved by Veritas IRB Ethics Review Board (Veritas IRB, Montreal, Canada. IRB tracking number: 16198–16:02:416–11-2017). Healthy volunteer donors were recruited by advertising in local communities and all participants gave their written informed consent.

Competing interests

The following authors are employees at Norgen Biotek Corp.: ME, BL, THA, DY, LN, PR and NR. YHA is the President and CEO of Norgen Biotek Corp. SM received an Industrial Undergraduate Student Research Award from the Natural Sciences and Engineering Research Council of Canada (NSERC) at Norgen Biotek Corp. Some of Norgen Biotek’s products have been used in the study, however the study is for basic scientific exploratory purposes and is not intended to promote or test any of Norgen Biotek’s products.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S2. Relative biotype distribution among the various bodily fluids of each donor. (TIF 566 kb)

Additional file 2:

Table S1. Common miRNAs between invasive fluids. (DOCX 14 kb)

Additional file 3:

Table S2. Common miRNAs between non-invasive fluids. (DOCX 13 kb)

Additional file 4:

Table S3. Common miRNAs between non-invasive fluids and blood. (DOCX 13 kb)

Additional file 5:

Figure S1. Venn diagram showing the overlap between blood and the non-invasive bodily fluids. (TIF 536 kb)

Additional file 6:

Table S4. Common miRNAs between all fluids. (DOCX 13 kb)

Additional file 7:

Table S5. Unique miRNAs detected in each bodily fluid. (DOCX 14 kb)

Additional file 8:

Table S6. Unique miRNAs detected in the invasive bodily fluids. (DOCX 14 kb)

Additional file 9:

Table S7. Unique miRNAs detected in the non-invasive bodily fluids. (DOCX 13 kb)

Additional file 10:

Table S8. Candidate novel miRNAs detected by miRDeep2 in the different bodily fluids and their matching result to the miRCarta database. (XLSX 43 kb)

Additional file 11:

Top 50 tDRs in each fluid with their quantification and coverage percentage. (XLSX 4259 kb)

Additional file 12:

Profiles of the top 50 mature tDR in each fluid. (PDF 4737 kb)

Additional file 13:

Table S9. piRNAs that represent an average of 1% or more of the entire piRNA counts of each bodily fluid. (DOCX 18 kb)

Additional file 14:

Figure S3. Exogenous mapping of unmapped saliva reads. (PDF 24 kb)

Additional file 15:

Figure S4. Exogenous mapping of unmapped cell-free saliva reads. (PDF 24 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

El-Mogy, M., Lam, B., Haj-Ahmad, T.A. et al. Diversity and signature of small RNA in different bodily fluids using next generation sequencing. BMC Genomics 19, 408 (2018). https://doi.org/10.1186/s12864-018-4785-8

Download citation

Keywords

  • miRNA
  • tRNA
  • piRNA
  • Next generation sequencing
  • Blood
  • Plasma
  • Serum
  • Saliva
  • Urine

BMC Genomics

ISSN: 1471-2164

Contact us