Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

The human miRNA repertoire of different blood compounds

Petra Leidinger1, Christina Backes2, Benjamin Meder3, Eckart Meese1 and Andreas Keller2*

Author Affiliations

1 Institute of Human Genetics, Saarland University, Homburg, Germany

2 Chair for Clinical Bioinformatics, Saarland University, Saarbrücken, Germany

3 Internal Medicine III, University Hospital Heidelberg, Heidelberg, Germany

For all author emails, please log on.

BMC Genomics 2014, 15:474  doi:10.1186/1471-2164-15-474


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/15/474


Received:8 April 2014
Accepted:6 June 2014
Published:14 June 2014

© 2014 Leidinger et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.

Abstract

Background

MiRNAs from body fluids gain more and more attraction as biomarker candidates. Besides serum, patterns from whole blood are increasingly considered as markers for human pathologies. Usually, the contribution of different cell types to the respective signature remains however unknown. In this study we provide insights into the human miRNome of different compounds of the blood including CD3, CD14, CD15, CD19, CD56 positive cells as well as exosomes.

Methods

We measured the miRNA repertoire for each cell type and whole blood for two individuals at three time points over the course of one year in order to provide evidence that the cell type miRNomes can be reproducibly detected.

Results

For measurements repeated after 24 hours we found on average correlation of 0.97, even after one year profiles still correlated with 0.96, demonstrating the enormous stability of the cell type specific miRNomes. Highest correlation was found for CD15 positive cells, exceeding Pearson correlation of 0.99. For exosomes a significantly higher variability of miRNA expression was detected. In order to estimate the complexity and variability of the cell type specific miRNomes, we generated profiles for all considered cell types in a total of seven unaffected individuals. While CD15 positive cells showed the most complex miRNome consisting of 328 miRNAs, we detected significantly less miRNAs (186, p = 1.5*10-5) in CD19 positive cells. Moreover, our analysis showed functional enrichment in many relevant categories such as onco-miRNAs and tumor miRNA suppressors. Interestingly, exosomes were enriched just for onco-miRNAs but not for miRNA tumor suppressors.

Conclusion

In sum, our results provide evidence that blood cell type specific miRNomes are very consistent between individuals and over time.

Keywords:
miRNA; Exosome; Cell separation

Background

Small non-coding RNAs (miRNAs) are involved in the majority of biological processes, including proliferation, apoptosis, differentiation, and development [1-3]. Moreover, the majority of human genes are regulated by miRNAs. Not surprisingly, miRNAs also play a crucial role in pathogenic processes.

In the beginning of miRNA research mostly tissue has been in the focus of research, thus for all human tissue types specific profiles have been discovered, e.g. from patients with lung cancer [4], breast cancer [5] or glioblastoma [6]. More recently body fluids as source for non-invasive or minimal-invasive markers have become more important. Besides serum (biomarker discovery for non-ischaemic systolic heart failure [7], pulmonary tuberculosis [8], non small-cell lung cancer [9,10], breast cancer [11], prostate cancer [12], or ovarian cancer [13]), one of the most frequently applied approaches remains to measure miRNAs from whole blood (biomarkers for myocardial infarction [14], lung cancer [15], multiple sclerosis [16,17], melanoma [18], ovarian cancer [19], COPD [20], glioblastoma [21] and Alzheimer’s Disease [22]). As miRNAs are known to be very specific for different tissues, they are likewise known to be specific for different compounds of the blood. The respective knowledge on the origin of blood-borne miRNA patterns has however been hardly considered.

In our study we aimed to provide miRNomes for common compounds of the human blood, including CD3, CD14, CD15, CD19, CD56 positive cells and exosomes. For all samples we also provided the measurement of whole blood using PAXgene tubes as comparison. In a first stage, we measured the miRNA repertoire for each cell type and whole blood for two individuals at three time points (t0; t1 = t0 + 1 year; t2 = t1 + 1 day) in order to get an estimation how stable cell type specific miRNomes are over time. Here, we detected a very high reproducibility and stability for all tested cell types.

In a second stage we extended the measurement to 7 individuals in order to understand the complexity and variability of cell type specific miRNomes. Our results showed that the cell type miRNomes are highly specific but comparably constant between different individuals. Thus, our study represents a standard repertoire of miRNAs in the considered cell types.

Results and discussion

Stability of the human blood compound miRNomes

A central question of our study was how reproducible and reliable can miRNomes of individuals be measured. To this end we picked one male and one female individual and collected blood at three time points: t0, t1, which was one year after t0 and t2, which was one day after t1. Blood cell type extraction and detection of the respective miRNomes was done immediately after blood collection. This set-up allows us to test for the reproducibility of cell type specific miRNomes over a very short period of time but at the same time offers to report the stability of the miRNA repertoire in individuals over a longer period.

For measurements repeated within 24 hours we found on average a high correlation of 0.97, indicating a reasonable reproducibility of the blood collection, cell type separation and miRNA profiling approach. Astonishingly, profiles correlated even after one year still very well. On average, Pearson correlation was as high as 0.96, demonstrating the enormous stability of the cell type specific miRNomes over a longer period of time.

Considering the two consecutive time points, the highest correlation for a specific cell type was reached by CD15 positive cells, exceeding 0.99. Just CD3 positive cells reached almost the same level of reproducibility. A similar behavior was also observed for the correlation within cell types in general. Here, CD3, CD15 and CD56 positive cells reached correlation of 0.98, CD14 positive cells and PAXgene whole blood reached average correlation of 0.97, CD19 positive cells showed lowest but still substantial correlation of 0.96.To understand how the different cell populations are related to each other we applied hierarchical clustering on the 50 most variable miRNAs. As the dendrogram in Figure 1 presents, different cell types cluster perfectly to each other. PAXgene blood shows the most different patterns, CD14 and CD15 positive cells clustered well together as well as CD3, CD19 and CD56 positive cells. Notably, for PAXgene blood and CD3 positive cells, individuals likewise cluster together. Regarding the three time points we just observed a tendency of the latter two time points to cluster well together, again demonstrating the high stability of blood cell type specific miRNomes over a longer period of time.

thumbnailFigure 1. Cluster Heatmap. The figure presents clustering in different blood cell types (represented by different colors on top of the heatmap) for two individuals (represented by two different shapes on top of heatmap) at three different timepoints (represented by numbers on top of heatmap).

Complexity and specificity of cell type specific miRNomes

In order to improve our understanding of the complexity and specificity of cell type miRNomes we extended our cohort to 7 unaffected individuals and measured the same cell fractions. The bar diagram in Figure 2 presents the average number of detected miRNAs per cell fraction. While CD15 positive cells showed the most complex miRNome consisting of as much as 328 of the measured 1,205 markers, we detected significantly less miRNAs (186, p = 1.5*10-5) in CD19 positive cells. CD56, CD14 and CD3 positive cells showed comparable number of miRNAs with 262, 247 and 223 detected miRNAs on average. The respective miRNomes per cell fraction are available in Supplemental Table 1.Based on the different complexity in miRNome we concluded that the respective profiles as such may also show a high degree of variability. To validate this hypothesis we focused on the miRNAs with highest variance among all cell types and calculated the expression as percentage of average expression over all cell types for each miRNA. As the spider diagram and heatmap in Figure 3 detail, we discovered 6 miRNAs being significantly higher expressed in CD19, CD3 and CD56 positive cells as compared to CD14 and CD15 positive cells: miR-150, miR-342-5p, miR-146a, miR-342-3p, miR-155 and miR-151-5p. 7 miRNAs showed the opposite behavior: miR-223*, miR-199b-5p, miR-582-5p, miR-424, miR-223, miR-338-3p and miR-340. These results are in line with our initial findings on the two individuals measured over 1 year, overall CD14 and CD15 positive cells have similar patterns as well as CD19, CD3 and CD56 positive cells. However, we also detect miRNAs that are highly expressed just in single cell types, such as miR-145, miR-143 and miR-3663-3p for CD15 positive cells. Vice versa, miR-363 and miR-181b are over expressed in CD56 positive cells and show significantly decreased expression in CD15 positive cells.

thumbnailFigure 2. Bar diagram detaining the number of different miRNAs in all tested cell compounds.

Table 1. Enriched clusters, families and functions

thumbnailFigure 3. Spider diagram showing the abundance of selected miRNAs in different cell types. Total blood cells have been used to normalize the data (100%).

Notably, we discovered for each cell fraction at least a single specific miRNA, i.e. a miRNA that was just expressed in this cell type (in two experiments) while not in any other cell type. For CD56 positive cells highest number of specific miRNAs was detected (25), in contrast, the other fractions showed substantially less specific markers: CD14 (3), CD56 (5), CD3 (4) CD19 (1). The respective specific markers are summarized in Figure 4. The full miRNome profile of the respective individuals for all cell types is available in Supplemental Table 2.

thumbnailFigure 4. miRNAs that are specific for the different cell compounds. The miRNAs below the respective compounds have been detected just in the respective compounds while not being available in any other type.

Table 2. Information for all blood donors

The miRNome profile of exosomes

Besides measuring whole miRNome profiles from different tissues or body fluids like whole blood, serum or urine, miRNA profiles of exosomes that are released by cells and circulate in body fluids become likewise more central for current research. Thus, we also extracted exosomes from serum of the same two individuals at t1 and t2 and measured the respective exosome miRNomes. Repeated extraction of exosomes of the same individuals correlated well to each other, on average, pearson correlation was above 0.95. While this still represents a reasonable stability of exosome profiles, the values were clearly below measurements for the different cell types, corresponding to a higher variability of miRNAs in exosomes.Altogether, we found just a small number of miRNAs being expressed in these microvesicles. On average, 128 miRNAs were detected, a number well below the average number of detected miRNAs in the different analyzed blood cell types from the same blood donors. Despite the lower complexity of exosome miRNA profiles we searched for miRNAs being present just in exosomes while not in any other cell fraction. Indeed, we detected markers present at least in two different exosome measurements while not in any blood cell type or whole blood of the same individuals. The respective 16 markers include hsa-miR-1299, hsa-miR-3124, hsa-miR-4290, hsa-miR-2278, hsa-miR-32*, hsa-miR-3149, hsa-miR-877*, hsa-miR-584, hsa-miR-3148, hsa-miR-122, hsa-miR-718, hsa-miR-670, hsa-miR-3680*, hsa-miR-193b*, hsa-miR-3911 and hsa-miR-1228* (see Figure 4).

Family and function of blood cell compound miRNAs

In order to find potentially enriched miRNA families or clusters as well as biological functions we used the freely available TAM tool [23]. Specifically, we analyzed all miRNAs for each blood compound that has been detected at least in 2 measurements. Our analysis revealed that in many cases specific miRNA families are enriched in selected cell types, for CD14 positive cells, miR-106a cluster was significant while for CD19 positive cells miR-106b cluster was found to be enriched. Interestingly, miR-17 cluster was significant for all cell types but not for exosomes. The same holds for the Akt pathway, Cell proliferation, Apoptosis, Granulopoiesis and others. Most interestingly, miRNA tumor suppressors were enriched in all categories besides exomes, however, in the case of onco-miRNAs, exosomes became likewise significant. All results are summarized in Table 1.

Conclusion

While for gene expression various studies have been carried out in order to systematically explore patterns in different cell types, for miRNAs few studies have been published [24,25]. As an example, Yu and co-workers demonstrated that positive and negative selection do not have a significant influence of the miRNA pattern [26]. In order to reduce the risk of contamination with other cells and to reach high purity of the isolated cell fractions we decided to apply positive cell selection. In two stages we then systematically explored the specific miRNomes of 5 cell fractions including CD3, CD14, CD15, CD19 and CD56. We also performed profiling of the same individuals using PAXgene blood tubes, representing one of the gold standard approaches for detecting circulating disease markers for various human pathologies. Even more importantly, we likewise generated profiles from exosomes of the respective individuals.

In summary our study reveals that 1) miRNA patterns of specific compounds of human blood can be reliably measured and are stable even over the course of one year; 2) the different cell types show substantial differences with respect to the complexity of their respective miRNomes; 3) CD14 and CD15 positive cells have a common myeloid progenitor and show a similar expression pattern, likewise CD3, CD19 and CD56 positive cells (common lymphoid progenitor) are generally related to each other; 4) Exosomes show a less complex and more variable miRNA profile as compared to other cell types.

Methods

Study set-up and miRNA profiling

The study set-up is presented in Figure 5. All participants in the study gave written informed consent. The local ethics committee (Ärztekammer des Saarlandes) approved the study (reference: 01/08). From all blood donors (details see Table 2), blood was collected in PAXgene RNA blood tubes (Becton Dickinson) and in S-Monovette® with EDTA K2 (Sarstedt). RNA from PAXgene whole blood was extracted according to manufacturers instruction using the PAXgene Blood miRNA Kit (Qiagen). For different cell types, positive selection was applied. While for CD3, CD19, CD15, and CD14 positive cells non-magnetic beads along with a sieve (pluriSelect) were used, for CD56 positive cells magnetic beads (Miltenyi Biotech) were used. The selection has been carried out using the leukocyte fraction. Exosome extraction has been carried out using ExoQuick Exosome Precipitation Solution (System Biosciences). Total RNA including small RNAs was isolated from different cell types and from exosomes using miRNeasy Micro Kit (Qiagen) according to manufacturer’s instructions.

thumbnailFigure 5. Study set-up. The diagram shows the different time points and cohort sizes for which the blood cell compounds were calculated.

Expression analysis was performed as previously described using Agilent’s Sure Print G3 Human v16 8x60K miRNA microarrays [27]. These arrays contain 40 replicates of each of the 1,205 miRNAs as annotated in miRBase v16 [28].

Bioinformatics analysis

The initial data analysis was performed using Agilent’s Feature Extraction software as described in the manufacturers instructions. For further down-stream analysis, the 40 replicates of background corrected miRNA measurements were summarized using Agilent’s Feature Extraction software. These measurements are available as raw data at GEO with the reference GSE56590. The further down-stream analysis has been carried out using the freely available statistical programming language R [29] in version 3.0.2.

In order to carry out hierarchical clustering and calculate heatmaps source code from the heatmap.2 function, provided as part of the “gplots” CRAN package (version 2.12.1) has been used. In more detail, hierarchical clustering relying on the Euclidian distance has been carried out on quantile normalized data (normalization has been done by the “preprocessCore” package using standard parameters). In order to estimate the linear dependencies between the different replicates the Pearson correlation coefficient has been calculated (standard “stats” package). If not mentioned explicitly, significance vales have been calculated using t-test (unpaired, two-tailed, homoscedastic variance). As alpha level, 0.05 has been used through the manuscript. In order to evaluate whether homoscedastic t-test is applicable to the data, normal distribution was verified using the Shapiro Wilk test and the variance has been tested using Bartletts test (both provided by the R “stats” package.

Finally, functional enrichment analysis has been carried out using TAM, the tool for annotations of microRNAs, which is freely accessible online (http://202.38.126.151/hmdd/tools/tam.html/ webcite[23]). TAM relies on the hypergeometric distribution to calculate p-values.

Availability of supporting data

The array data have been deposited in the Gene expression Omnibus GEO and are freely accessible under reference GSE56590.

Competing interest

The authors declare that they have no competing interests.

Authors’ contributions

PL carried out the experiments, CB contributed in data analysis, BM contributed in experimental design and study set up, EM wrote the manuscript, AK wrote the manuscript and contributed in data analysis. All authors read and approved the final manuscript.

Acknowledgement

This work has been funded in parts by Deutsche Forschungsgemeinschaft (DFG, LE2783/1-1) and the EU as part of the BestAgeing project.

References

  1. Cheng AM, Byrom MW, Shelton J, Ford LP: Antisense inhibition of human miRNAs and indications for an involvement of miRNA in cell growth and apoptosis.

    Nucleic Acids Res 2005, 33(4):1290-1297. OpenURL

  2. Pogue AI, Cui JG, Li YY, Zhao Y, Culicchia F, Lukiw WJ: Micro RNA-125b (miRNA-125b) function in astrogliosis and glial cell proliferation.

    Neurosci Lett 2010, 476(1):18-22. OpenURL

  3. Yoshizawa M, Taguchi YH, Yasuda J: Inference of Gene Regulation via miRNAs During ES Cell Differentiation Using MiRaGE Method.

    Int J Mol Sci 2012, 12(12):9265-9276. OpenURL

  4. Yu SL, Chen HY, Chang GC, Chen CY, Chen HW, Singh S, Cheng CL, Yu CJ, Lee YC, Chen HS, Su TJ, Chiang CC, Li HN, Hong QS, Su HY, Chen CC, Chen WJ, Liu CC, Chan WK, Li KC, Chen JJ, Yang PC: MicroRNA signature predicts survival and relapse in lung cancer.

    Cancer Cell 2008, 13(1):48-57. OpenURL

  5. Ma L, Teruya-Feldstein J, Weinberg RA: Tumour invasion and metastasis initiated by microRNA-10b in breast cancer.

    Nature 2007, 449(7163):682-688. OpenURL

  6. Chan JA, Krichevsky AM, Kosik KS: MicroRNA-21 is an antiapoptotic factor in human glioblastoma cells.

    Cancer Res 2005, 65(14):6029-6033. OpenURL

  7. Vogel B, Keller A, Frese KS, Leidinger P, Sedaghat-Hamedani F, Kayvanpour E, Kloos W, Backe C, Thanaraj A, Brefort T, Beier M, Hardt S, Meese E, Katus HA, Meder B: Multivariate miRNA signatures as biomarkers for non-ischaemic systolic heart failure.

    Eur Heart J 2013, 34(36):2812-2823. OpenURL

  8. Fu Y, Yi Z, Wu X, Li J, Xu F: Circulating microRNAs in patients with active pulmonary tuberculosis.

    Journal of clinical microbiology 2011, 49(12):4246-4251. OpenURL

  9. Keller A, Leidinger P, Gislefoss R, Haugen A, Langseth H, Staehler P, Lenhof HP, Meese E: Stable serum miRNA profiles as potential tool for non-invasive lung cancer diagnosis.

    RNA Biol 2011, 8(3):506-516. OpenURL

  10. Wang ZX, Bian HB, Wang JR, Cheng ZX, Wang KM, De W: Prognostic significance of serum miRNA-21 expression in human non-small cell lung cancer.

    J Surg Oncol 2011, 104(7):847-851. OpenURL

  11. Gotte M: MicroRNAs in breast cancer pathogenesis.

    Minerva ginecologica 2010, 62(6):559-571. OpenURL

  12. Yin Y, Li M, Li H, Jiang Y, Cao LY, Zhang HF, Xu XC: [Expressions of 6 microRNAs in prostate cancer].

    Zhonghua nan ke xue = National journal of andrology 2010, 16(7):599-605. OpenURL

  13. Resnick KE, Alder H, Hagan JP, Richardson DL, Croce CM, Cohn DE: The detection of differentially expressed microRNAs from the serum of ovarian cancer patients using a novel real-time PCR platform.

    Gynecol Oncol 2009, 112(1):55-59. OpenURL

  14. Meder B, Keller A, Vogel B, Haas J, Sedaghat-Hamedani F, Kayvanpour E, Just S, Borries A, Rudloff J, Leidinger P, Meese E, Katus HA, Rottbauer W: MicroRNA signatures in total peripheral blood as novel biomarkers for acute myocardial infarction.

    Basic Res Cardiol 2011, 106(1):13-23. OpenURL

  15. Keller A, Leidinger P, Borries A, Wendschlag A, Wucherpfennig F, Scheffler M, Huwer H, Lenhof HP, Meese E: miRNAs in lung cancer - studying complex fingerprints in patient's blood cells by microarray experiments.

    BMC Cancer 2009, 9:353. OpenURL

  16. Keller A, Leidinger P, Lange J, Borries A, Schroers H, Scheffler M, Lenhof HP, Ruprecht K, Meese E: Multiple sclerosis: microRNA expression profiles accurately differentiate patients with relapsing-remitting disease from healthy controls.

    PLoS One 2009, 4(10):e7440. OpenURL

  17. Keller A, Leidinger P, Steinmeyer F, Stahler C, Franke A, Hemmrich-Stanisak G, Kappel A, Wright I, Dorr J, Paul F, Diem R, Tocariu-Krick B, Meder B, Backes C, Meese E, Ruprecht K: Comprehensive analysis of microRNA profiles in multiple sclerosis including next-generation sequencing.

    Mult Scler 2014, 20(3):295-303. OpenURL

  18. Leidinger P, Keller A, Borries A, Reichrath J, Rass K, Jager SU, Lenhof HP, Meese E: High-throughput miRNA profiling of human melanoma blood samples.

    BMC Cancer 2010, 10:262. OpenURL

  19. Hausler SF, Keller A, Chandran PA, Ziegler K, Zipp K, Heuer S, Krockenberger M, Engel JB, Honig A, Scheffler M, Dietl J, Wischhusen J: Whole blood-derived miRNA profiles as potential new tools for ovarian cancer screening.

    Br J Cancer 2010, 103(5):693-700. OpenURL

  20. Leidinger P, Keller A, Borries A, Huwer H, Rohling M, Huebers J, Lenhof HP, Meese E: Specific peripheral miRNA profiles for distinguishing lung cancer from COPD.

    Lung Cancer 2011, 74(1):41-47. OpenURL

  21. Roth P, Wischhusen J, Happold C, Chandran PA, Hofer S, Eisele G, Weller M, Keller A: A specific miRNA signature in the peripheral blood of glioblastoma patients.

    J Neurochem 2011, 118(3):449-457. OpenURL

  22. Leidinger P, Backes C, Deutscher S, Schmitt K, Mueller SC, Frese K, Haas J, Ruprecht K, Paul F, Stahler C, Lang CJ, Meder B, Bartfai T, Meese E, Keller A: A blood based 12-miRNA signature of Alzheimer disease patients.

    Genome Biol 2013, 14(7):R78. OpenURL

  23. Lu M, Shi B, Wang J, Cao Q, Cui Q: TAM: a method for enrichment and depletion analysis of a microRNA category in a list of microRNAs.

    BMC Bioinformatics 2010, 11:419. OpenURL

  24. Beliakova-Bethell N, Massanella M, White C, Lada SM, Du P, Vaida F, Blanco J, Spina CA, Woelk CH: The effect of cell subset isolation method on gene expression in leukocytes.

    Cytometry A 2014, 85(1):94-104. OpenURL

  25. Debey S, Schoenbeck U, Hellmich M, Gathof BS, Pillai R, Zander T, Schultze JL: Comparison of different isolation techniques prior gene expression profiling of blood derived cells: impact on physiological responses, on overall expression and the role of different cell types.

    Pharmacogenomics J 2004, 4(3):193-207. OpenURL

  26. Yu IW, Espinoza DA, McAlexander MA, Witwer KW: OpenArray profiling reveals no differential modulation of miRNA by positive and negative CD4+ T cell immunoselection.

    Exp Hematol 2014, 42(1):11-13. OpenURL

  27. Leidinger P, Keller A, Backes C, Huwer H, Meese E: MicroRNA expression changes after lung cancer resection: A follow-up study.

    RNA Biol 2012, 9(6):900-910. OpenURL

  28. Griffiths-Jones S: miRBase: the microRNA sequence database.

    Methods Mol Biol 2006, 342:129-138. OpenURL

  29. Team R: R: A Language and Environment for Statistical Computing. In. Vienna: R Foundation for Statistical Computing; 2008. OpenURL