- Research article
- Open Access
Diversity and signature of small RNA in different bodily fluids using next generation sequencing
BMC Genomics volume 19, Article number: 408 (2018)
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).
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.
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.
Small RNAs are a class of mainly non-coding RNAs (ncRNAs) characterized by their small nucleotide length of less than 200 nt . 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) . 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 . To date, 2588 human mature miRNAs have been identified and are currently included in miRBase 21 . 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 . 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 . 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 . 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 . Transfer RNA-derived small RNAs are thought to have a dual function, as they act as suppressors and oncogenes . 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.
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 . 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 . 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).
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) , 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 220.127.116.11) was used to predict novel miRNA candidates and tDRmapper was used to identify tDRs.
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).
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.
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.
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.
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).
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).
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).
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 . 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 . 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).
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).
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.
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).
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.
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. . 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 . 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 . 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 . 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 . 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  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 , gastric cancer  and oral tongue squamous cell carcinoma . It may act as a tumor suppressor miRNA  and may also be used to predict the efficacy of cancer vaccine treatment for colorectal cancer . 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 . It is differentially expressed in senescence  and as a tumor suppressor in gliomas . 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 . 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 . 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 . 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 . 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 . The most abundant piRNA molecule in urine and cell-free urine was hsa-piR-019825, which is deregulated in colorectal cancer patients . 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 . They are potential cancer biomarkers and regulators of inflammation and cell death [99, 100].
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.
Counts per million
Long intergenic noncoding RNA
Next generation sequencing
Principal component analysis
RNA integrity number
Small nucleolar RNA
Small nuclear RNA
Y RNA-derived small RNAs
Trimmed mean of M-values
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
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
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
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
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
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell [internet]. 2004;116:281–297. Available from: http://www.ncbi.nlm.nih.gov/pubmed/14744438
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
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
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
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
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
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
miRBase: the microRNA database [Internet]. [cited 2018 Mar 7]. Available from: http://www.mirbase.org/cgi-bin/browse.pl?org=hsa
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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.
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
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
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
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
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
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
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
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/
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Dhahbi, Spinder S, Atamna H, Boffelli D, Martin D. Deep Sequencing of Serum Small RNAs Identifies Patterns of 5&#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
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
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.
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.
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.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S2. Relative biotype distribution among the various bodily fluids of each donor. (TIF 566 kb)
Table S1. Common miRNAs between invasive fluids. (DOCX 14 kb)
Table S2. Common miRNAs between non-invasive fluids. (DOCX 13 kb)
Table S3. Common miRNAs between non-invasive fluids and blood. (DOCX 13 kb)
Figure S1. Venn diagram showing the overlap between blood and the non-invasive bodily fluids. (TIF 536 kb)
Table S4. Common miRNAs between all fluids. (DOCX 13 kb)
Table S5. Unique miRNAs detected in each bodily fluid. (DOCX 14 kb)
Table S6. Unique miRNAs detected in the invasive bodily fluids. (DOCX 14 kb)
Table S7. Unique miRNAs detected in the non-invasive bodily fluids. (DOCX 13 kb)
Table S8. Candidate novel miRNAs detected by miRDeep2 in the different bodily fluids and their matching result to the miRCarta database. (XLSX 43 kb)
Top 50 tDRs in each fluid with their quantification and coverage percentage. (XLSX 4259 kb)
Profiles of the top 50 mature tDR in each fluid. (PDF 4737 kb)
Table S9. piRNAs that represent an average of 1% or more of the entire piRNA counts of each bodily fluid. (DOCX 18 kb)
Figure S3. Exogenous mapping of unmapped saliva reads. (PDF 24 kb)
Figure S4. Exogenous mapping of unmapped cell-free saliva reads. (PDF 24 kb)
About this article
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
- Next generation sequencing