- Open Access
Diagnostic and prognostic microRNAs in the serum of breast cancer patients measured by droplet digital PCR
Biomarker Research volume 3, Article number: 12 (2015)
Breast cancer circulating biomarkers include carcinoembryonic antigen and carbohydrate antigen 15–3, which are used for patient follow-up. Since sensitivity and specificity are low, novel and more useful biomarkers are needed. The presence of stable circulating microRNAs (miRNAs) in serum or plasma suggested a promising role for these tiny RNAs as cancer biomarkers. To acquire an absolute concentration of circulating miRNAs and reduce the impact of preanalytical and analytical variables, we used the droplet digital PCR (ddPCR) technique.
We investigated a panel of five miRNAs in the sera of two independent cohorts of breast cancer patients and disease-free controls. The study showed that miR-148b-3p and miR-652-3p levels were significantly lower in the serum of breast cancer patients than that in controls in both cohorts. For these two miRNAs, the stratification of breast cancer patients versus controls was confirmed by receiver operating characteristic curve analyses. In addition, we showed that higher levels of serum miR-10b-5p were associated with clinicobiological markers of poor prognosis.
The study revealed the usefulness of the ddPCR approach for the quantification of circulating miRNAs. The use of the ddPCR quantitative approach revealed very good agreement between two independent cohorts in terms of comparable absolute miRNA concentrations and consistent trends of dysregulation in breast cancer patients versus controls. Overall, this study supports the use of the quantitative ddPCR approach for monitoring the absolute levels of diagnostic and prognostic tumor-specific circulating miRNAs.
Breast cancer is the most frequent cancer and the second leading cause of cancer death among women in industrialized countries. Approximately 1.3 million women develop breast cancer every year . Advances in early diagnosis and treatments have contributed to the decrease of mortality rates over the years. The overall 5-year survival is 90 % when breast cancer is diagnosed at an early stage as opposed to 20 % if disease has spread to distant organs . Physical examination, mammography, and biopsy are the current approaches to breast cancer diagnosis . Carcinoembryonic antigen and carbohydrate antigen 15–3 are circulating tumor markers that are mainly used for patient follow-up . However, the sensitivity of these markers is low, thus calling attention to the need for novel and more accurate noninvasive diagnostic biomarkers.
Studies on circulating microRNAs (miRNAs) opened potential opportunities for the discovery of new tumor biomarkers. miRNAs are a class of small noncoding RNAs that regulate gene expression at the post-transcriptional level . They play a crucial role in the regulation of most, if not all, human genes and their involvement in the deregulation of pathological states such as cancer has been well established . Moreover, miRNAs can be detected in serum or plasma, and their levels may be specifically altered in pathological conditions. Because of their remarkable stability in plasma and serum and the possibility of measuring their levels using noninvasive methods, various studies have suggested a role for circulating miRNAs as novel cancer biomarkers [7–11].
Despite promising results, however, it became evident that several variables (sample collection and storage, RNA purification methods, quantification and normalization methods) could affect final results . In this study, we took advantage of the droplet digital PCR (ddPCR) technique for assessing circulating miRNA levels, an approach that allows absolute quantification without the need for internal/external normalization. Using ddPCR, we investigated five miRNAs in the sera of two independent cohorts of breast cancer patients and disease-free controls to verify whether miRNAs could represent useful diagnostic biomarkers of breast cancer.
Circulating miRNAs in sera of breast cancer patients versus healthy controls
We selected five miRNAs (miR-10b-5p, miR-145-5p, miR-148b-3p, miR-425-5p, miR-652-3p) derived from microarray experiments  or described in recently published scientific literature as being potential circulating biomarkers (Additional file 1: Table S1). Two independent sets of serum samples from breast cancer patients and disease-free controls were analyzed (Table 1). One group of samples was collected at the University Hospital of Ferrara, Italy, from 2012 to 2014 (cohort A), while the second group of samples was collected at the Mercy’s Woman Center in Oklahoma City, OK, USA, from 2005 to 2013 (cohort B). Serum samples from both cohorts were collected and processed according to the same protocol, and the levels of circulating miRNAs were assessed by ddPCR. This technique allows the measurement of the absolute concentration of circulating miRNAs with no need for a reference gene, a condition particularly important for this type of sample. Hence, miRNA levels were expressed as copies per microliter of serum.
miR-148b-3p and miR-652-3p levels were significantly lower in breast cancer patients than in controls in both cohorts (p = 0.0042 and p < 0.0001, respectively, in cohort A; p = 0.0115 and p = 0.0043, respectively, in cohort B) (Fig. 1). miR-145-5p and miR-425-5p were also down-regulated in breast cancer patients compared with controls in both cohorts, but the differences were statistically significant in cohort B only (miR-145-5p: p = 0.0257, miR-425-5p: p = 0.0226) (Fig. 1). Conversely, miR-10b-5p exhibited a weak increase in cancer patients compared with controls. This trend was statistically significant in cohort B only (p = 0.016) (Fig. 1). The reduced representation of high tumor stages was likely responsible for the lack of statistical significance for the results of cohort A (see next section). When combined, the two cohorts produced a highly significant discrimination between cancer patients and controls for all investigated miRNAs, with the strongest discrimination achieved by miR-652 and miR-148b (Additional file 1: Figure S1). To validate these results in comparison with the most commonly used method based on Real Time PCR (RT-PCR) approach, we investigated miR-652 and miR-10b. As RT-PCR normalizer, we employed the non-human Cel-miR-39, which we routinely add at a defined concentration to any serum sample. Results of these analyses confirmed the significant discrimination between samples from breast cancer patients versus controls in the same direction shown by ddPCR-based analyses, thus showing that ddPCR did not introduce any experimental bias if compared to quantitative RT-PCR (Additional file 1: Figure S2).
Receiver operating characteristic (ROC) curve analysis was performed to evaluate the diagnostic value of the five miRNAs (Fig. 2). miR-652 and miR-148b appeared to represent valuable diagnostic biomarkers. miR-652 was of particular interest because of the highly significant ROC curves in both cohorts. Here, we confirmed a significant lower level of miR-652-3p not only in Luminal A cancer patients (estrogen receptor [ER]/progesterone receptor [PR] positive and human epidermal growth factor receptor 2 [HER2] negative), as previously reported ], but also in non-Luminal A cancer patients, versus controls (p = 0.020 for Luminal A and p = 0.004 for non-Luminal A) (Additional file 1: Figure S3).
Association of miR-10b with prognostic parameters
Associations between each miRNA and clinicopathological features were investigated in patients of cohort B (Table 2). Since cohort A included only patients with stage I or II tumors, it was excluded from these analyses. In cohort B, the level of serum miR-10b-5p revealed a concordant increase with tumor stage (Fig. 3). Patients with stage II to IV cancers exhibited significantly higher levels of miR-10b in comparison with patients with a stage I tumor (p = 0.0047) or with controls (p = 0.0028). Conversely, no significant difference was found between stage I patients and controls.
Notably, miR-10b-5p was also significantly up-regulated in association with other clinicopathological features of prognostic significance, including higher tumor grading and lymph node metastases (Fig. 3). Albeit not statistically significant, the average level of miR-10b-5p was also higher in patients carrying HER2-positive or ER/PR-negative cancers. The number of triple negative breast cancer were too few to lead to any meaningful result.
These data indicate that miR-10b-5p is significantly associated with parameters associated with a poor prognosis. No significant association was found between miR-145-5p, miR-148b-3p, miR-425-5p, or miR-652-3p and any clinicopathological feature.
The presence of stable miRNAs circulating in plasma or serum suggested their potential use as noninvasive biomarkers in cancer patients. In the last few years, several authors have demonstrated that a number of circulating miRNAs could discriminate breast cancer patients from healthy individuals [15, 16] or could be linked to breast cancer subtypes [14, 17]. Unfortunately, the combined effects of several variables made results poorly reproducible and their translation into clinically useful applications not feasible [12, 18–20]. For example, the heterogeneity of investigated populations (age, tumor features), or the differences in sample type (plasma or serum) or sample processing protocols could have been responsible for the apparent discrepancies among the various studies. Moreover, a variety of methods were used to normalize data, thereby producing non-comparable or difficult-to-compare results .
Here, we used the ddPCR technique to measure circulating miRNAs. ddPCR is a technique that can achieve absolute quantification of nucleic acids by combining limiting dilutions, end-point PCR, and Poisson statistics. In fact, the partitioning of the PCR reaction into up to 20,000 separate droplets mimics a binary distribution of the target. More important, being an end-point PCR, ddPCR can tolerate wide variations in amplification efficiencies without affecting copy number estimation of the target [21–24].
Using ddPCR, we analyzed the levels of five miRNAs in the serum of two independent cohorts of breast cancer patients and disease-free controls. Blood samples were collected at two independent institutions and processed separately, with no differences in procedures for obtaining serum samples. We analyzed serum, instead of plasma, as it is the most commonly available patient material, and the procedure used to collect serum is homogeneous at different institutions, thus helping to reduce uncertainties in preanalytical procedures. Notably, all of the analyzed miRNAs showed comparable absolute levels in the sera of the two cohorts (see Fig. 1). Most important, both cohorts exhibited consistent trends of dysregulation in breast cancer patients versus controls. The differences between breast cancer cases and controls in cohort B were statistically significant for all five miRNAs, whereas only miR-148b-3p and miR-652-3p reached statistical significance in cohort A, possibly because of differences in clinicopathological characteristics. We performed ROC curve analyses in which we evaluated the possible diagnostic potential of the five circulating miRNAs. Areas under the curves and P-values were significant for miR-148b-3p and miR-652-3p in both cohorts (see Fig. 2), suggesting the potential value of these two miRNAs as breast cancer biomarkers.
Our data on circulating miRNAs show both similarities to and discrepancies from those of previous reports (listed in Additional file 1: Table S1). However, differences in experimental settings and technical approaches make it difficult to compare the findings from various reports. Similar to our findings, three previous reports found a decreased level of miR-145-5p in breast cancer patients [17, 25, 26]. The up-regulation of miR-10b-5p in breast cancer patients in our study is also in agreement with several published reports [16, 25, 27]. We did not find significant differences in cohort A. However, since this group consisted of patients carrying stage I or II cancers, this finding does not contradict a positive correlation between miR-10b and more advanced disease; it is also consistent with the results described in a report by Roth et al., who found a higher level of miR-10b in patients with metastatic disease . For mir-652-3p, our results are in agreement with one report that indicated a decreased level in patients with Luminal A-like breast cancer in comparison with controls . Conversely, Cuk et al. showed increased levels of miR-652-3p, likely because the study was performed on plasma instead of serum samples . In support of this suggestion, analysis of a small number (n = 20) of plasma samples from our breast cancer cohort revealed that miR-652-3p was indeed increased in breast cancer patients compared with that in controls (data not shown). Finding an opposite trend in serum or plasma samples is not new [12, 13] and raises the question about the different genesis of circulating miRNAs . Concerning miR-148b-3p and miR-425-5p, the published literature reports opposite results to ours, as both miRNAs were found to be higher in breast cancer patients than in healthy controls. However, studies on miR-148b were all performed by using plasma [15, 29, 30], whereas the miR-425 study was performed with serum, but normalization was based on the mean of assays performed on all samples ; thus, these results are not directly comparable with ours.
In the search for possible correlations with clinicopathological features, we found that miR-10b-5p levels were increased in the serum of patients with a high cancer stage or grade, or with the presence of lymph node metastases. Albeit not statistically significant, the average expression of miR-10b-5p was also higher in cases with ER/PR-negative tumors and in cases with HER2-positive tumors. Together, these findings indicate that miR-10b-5p represents a biomarker of tumor aggressiveness. This suggestion is also supported by other studies that indicated higher levels of circulating miR-10b in patients with metastatic breast cancer and worsening clinical stage [10, 27, 31]. It is notable that the role of miR-10b in invasion and metastasis has been thoroughly investigated and its importance proven [32–35]. miR-10b has also been found to be highly expressed in the vasculature of high-grade breast cancer . Taken together, these findings suggest that tumor and microenvironment features may be directly responsible for the increased levels of circulating miR-10b in the bloodstream of patients with advanced breast cancer. The possibility of assessing the absolute levels of miR-10b in the serum of patients by using a robust technique such as ddPCR could represent a potential new approach for monitoring disease behavior in breast cancer patients.
Overall, this study supports the use of the quantitative ddPCR approach for monitoring the absolute levels of specific miRNAs as diagnostic and prognostic serum biomarkers in breast cancer patients.
Two cohorts of patients were investigated. Serum samples from cohort A (n = 55) were collected at the General Surgery Unit of the University Hospital of Ferrara, Italy, from 28 breast cancer patients and 27 age-matched disease-free controls. Serum samples from cohort B (n = 94) were collected at the Mercy Women’s Center in Oklahoma City, OK, USA, from 59 breast cancer patients and 35 age-matched controls (Additional file 1: Table S2). Ethical approval was granted by the ethics committees of the respective institutions. Written informed consent was obtained from all individuals enrolled in the study. Table 1 summarizes the clinicopathological features of the patients.
Sample preparation and RNA purification from serum
Blood samples from cohort A and B were collected in red stopper clot tubes (Greiner Bio-One VACUETTE in cohort A, BD Vacutainer in cohort B) and processed within 1 h; they were centrifuged at 1000 g for 10 min at room temperature, and serum was stored at −80 °C in 200 μL aliquots until use. 3 μL of 4.16 nM solution of synthetic miRNA cel-miR-39-3p (ucaccggguguaaaucagcuug) from C.Elegans (synthesized by IDT) was added to each aliquot. Total RNA was isolated from 200 μL of serum by using the MiRNeasy kit (Qiagen). RNA was eluted from spin columns in 35 μL of nuclease-free water.
Reverse transcription, ddPCR and Real-Time PCR
cDNA was synthesized in a 20 μL reaction by using the Universal cDNA synthesis kit II (Exiqon), starting from 3 μL of RNA according to the manufacturer’s guidelines for serum and plasma samples. Synthesized cDNA was diluted 50-fold, and 8 μL was assayed in a 20 μL PCR reaction volume according to the manufacturer’s protocol for miRCURY LNA Universal RT microRNA PCR (Exiqon) with EvaGreen (Bio-Rad). Each PCR reaction was mixed with 70 μL of droplet generator oil for EvaGreen in a disposable cartridge and applied to the QX200 droplet generator device (Bio-Rad) that portioned each sample into 20,000 nanoliter-sized droplets. Each sample was then transferred into a 96-well PCR plate and PCR was performed to the end point according to the manufacturer’s protocol (Bio-Rad). The same procedure was applied to all test samples and negative controls. At the end of the PCR reaction, the QX200 droplet reader (Bio-Rad) was used to count PCR-positive and PCR-negative droplets: the “singulator” unpacks the emulsified droplets and streams them in a single line past a two-color optical detection system. Positive droplets, which contain at least one copy of the target miRNA, exhibit increased fluorescence compared with negative droplets. The fraction of PCR-positive droplets enables the target to be quantified according to Poisson distribution. Quantitative Real-Time PCR was performed on miR-10b-5p and miR-652-3p using the ExiLent CyberGreen mastermix (Exiqon). Cel-miR-39 was used to normalize miRNAs levels with the 2^-ΔΔCt method.
Statistical analysis was performed by using Prism software version 5.0 (GraphPad, La Jolla, CA). An unpaired t-test with Welch’s correction was performed to assess the significance of differences between data distribution. A p-value of less than 0.05 was deemed to be significant. Binary logistic regression analysis was performed and ROC curves were generated to evaluate the ability of the chosen miRNAs to discriminate cancer cases versus controls.
Droplet digital PCR
Human epidermal growth factor receptor 2
Receiver operating characteristic
Area under the curve
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61:69–90.
Howlader N, Noone A, Krapcho M, Garshell J, Miller D, Altekruse S, et al. SEER Cancer Statistics Review, 1975–2011. http://seer.cancer.gov/csr/1975_2011/. National Cancer Institute: Bethesda, MD; 2014.
Checka CM, Chun JE, Schnabel FR, Lee J, Toth H. The relationship of mammographic density and age: implications for breast cancer screening. AJR Am J Roentgenol. 2012;198:W292–5.
Harris L, Fritsche H, Mennel R, Norton L, Ravdin P, Taube S, et al. American Society of Clinical Oncology 2007 update of recommendations for the use of tumor markers in breast cancer. J Clin Oncol. 2007;25:5287–312.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Negrini M, Ferracin M, Sabbioni S, Croce CM. MicroRNAs in human cancer: from research to therapy. J Cell Sci. 2007;120:1833–40.
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 U S A. 2008;105:10513–8.
Turchinovich A, Weiz L, Langheinz A, Burwinkel B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011;39:7223–33.
Lawrie CH, Gal S, Dunlop HM, Pushkaran B, Liggins AP, Pulford K, et al. Detection of elevated levels of tumour-associated microRNAs in serum of patients with diffuse large B-cell lymphoma. Br J Haematol. 2008;141:672–5.
Roth C, Rack B, Muller V, Janni W, Pantel K, Schwarzenbach H. Circulating microRNAs as blood-based markers for patients with primary and metastatic breast cancer. Breast Cancer Res. 2010;12:R90.
Zhao H, Shen J, Medico L, Wang D, Ambrosone CB, Liu S. A pilot study of circulating miRNAs as potential biomarkers of early stage breast cancer. PLoS One. 2010;5, e13735.
Jarry J, Schadendorf D, Greenwood C, Spatz A, van Kempen LC. The validity of circulating microRNAs in oncology: five years of challenges and contradictions. Mol Oncol. 2014;8:819–29.
Ferracin M, Lupini L, Salamon I, Saccenti E, Musa G, Zagatti B, et al. Absolute quantification of cell-free microRNAs in cancer patients. Oncotarget 2015, in press.
McDermott AM, Miller N, Wall D, Martyn LM, Ball G, Sweeney KJ, et al. Identification and validation of oncologic miRNA biomarkers for luminal A-like breast cancer. PLoS One. 2014;9:e87032.
Cuk K, Zucknick M, Heil J, Madhavan D, Schott S, Turchinovich A, et al. Circulating microRNAs in plasma as early detection markers for breast cancer. Int J Cancer. 2013;132:1602–12.
Mar-Aguilar F, Mendoza-Ramirez JA, Malagon-Santiago I, Espino-Silva PK, Santuario-Facio SK, Ruiz-Flores P, et al. Serum circulating microRNA profiling for identification of potential breast cancer biomarkers. Dis Markers. 2013;34:163–9.
Kodahl AR, Lyng MB, Binder H, Cold S, Gravgaard K, Knoop AS, et al. Novel circulating microRNA signature as a potential non-invasive multi-marker test in ER-positive early-stage breast cancer: a case control study. Mol Oncol. 2014;8:874–83.
Cheng HH, Yi HS, Kim Y, Kroh EM, Chien JW, Eaton KD, et al. Plasma processing conditions substantially influence circulating microRNA biomarker levels. PLoS One. 2013;8:e64795.
Mestdagh P, Hartmann N, Baeriswyl L, Andreasen D, Bernard N, Chen C, et al. Evaluation of quantitative miRNA expression platforms in the microRNA quality control (miRQC) study. Nat Methods. 2014;11:809–15.
Sourvinou IS, Markou A, Lianidou ES. Quantification of circulating miRNAs in plasma: effect of preanalytical and analytical parameters on their isolation and stability. J Mol Diagn. 2013;15:827–34.
Dingle TC, Sedlak RH, Cook L, Jerome KR. Tolerance of droplet-digital PCR vs real-time quantitative PCR to inhibitory substances. Clin Chem. 2013;59:1670–2.
Hindson CM, Chevillet JR, Briggs HA, Gallichotte EN, Ruf IK, Hindson BJ, et al. Absolute quantification by droplet digital PCR versus analog real-time PCR. Nat Methods. 2013;10:1003–5.
Miotto E, Saccenti E, Lupini L, Callegari E, Negrini M, Ferracin M. Quantification of circulating mirnas by droplet digital PCR: comparison of evagreen- and taqman-based chemistries. Cancer Epidemiol Biomarkers Prev. 2014;23:2638–42.
Pinheiro LB, Coleman VA, Hindson CM, Herrmann J, Hindson BJ, Bhat S, et al. Evaluation of a droplet digital polymerase chain reaction format for DNA copy number quantification. Anal Chem. 2012;84:1003–11.
Chan M, Liaw CS, Ji SM, Tan HH, Wong CY, Thike AA, et al. Identification of circulating microRNA signatures for breast cancer detection. Clin Cancer Res. 2013;19:4477–87.
Ng EK, Li R, Shin VY, Jin HC, Leung CP, Ma ES, et al. Circulating microRNAs as specific biomarkers for breast cancer detection. PLoS One. 2013;8:e53141.
Zhao FL, Hu GD, Wang XF, Zhang XH, Zhang YK, Yu ZS. Serum overexpression of microRNA-10b in patients with bone metastatic primary breast cancer. J Int Med Res. 2012;40:859–66.
Etheridge A, Gomes CP, Pereira RW, Galas D, Wang K. The complexity, function and applications of RNA in circulation. Front Genet. 2013;4:115.
Cuk K, Zucknick M, Madhavan D, Schott S, Golatta M, Heil J, et al. Plasma microRNA panel for minimally invasive detection of breast cancer. PLoS One. 2013;8:e76729.
Shen J, Hu Q, Schrauder M, Yan L, Wang D, Medico L, et al. Circulating miR-148b and miR-133a as biomarkers for breast cancer detection. Oncotarget. 2014;5:5284–94.
Chen W, Cai F, Zhang B, Barekati Z, Zhong XY. The level of circulating miRNA-10b and miRNA-373 in detecting lymph node metastasis of breast cancer: potential biomarkers. Tumour Biol. 2013;34:455–62.
Ma L, Reinhardt F, Pan E, Soutschek J, Bhat B, Marcusson EG, et al. Therapeutic silencing of miR-10b inhibits metastasis in a mouse mammary tumor model. Nat Biotechnol. 2010;28:341–7.
Ma L, Teruya-Feldstein J, Weinberg RA. Tumour invasion and metastasis initiated by microRNA-10b in breast cancer. Nature. 2007;449:682–8.
Biagioni F, Bossel Ben-Moshe N, Fontemaggi G, Yarden Y, Domany E, Blandino G. The locus of microRNA-10b: a critical target for breast cancer insurgence and dissemination. Cell Cycle. 2013;12:2371–5.
Liu Y, Zhao J, Zhang PY, Zhang Y, Sun SY, Yu SY, et al. MicroRNA-10b targets E-cadherin and modulates breast cancer metastasis. Med Sci Monit. 2012;18:BR299–308.
Plummer PN, Freeman R, Taft RJ, Vider J, Sax M, Umer BA, et al. MicroRNAs regulate tumor angiogenesis modulated by endothelial progenitor cells. Cancer Res. 2013;73:341–52.
This work was supported by a grant from the Italian Ministry of Research and the University, FIRB 2011 (Project RBAPIIBYNP), and from the University of Ferrara (FAR programs 2012–14) to MN, and by a grant from the Italian Association for Cancer Research MFAG (Project 11676) to MF.
The authors have no interest, sponsorship, or affiliation that could be perceived as a conflict of interest in the context of this research.
AM, ES, SOE, and VVP performed laboratory experiments; MP and PQ produced and collected histopathological information at the University of Ferrara Hospital; AF clinically followed up patients at the University of Ferrara Hospital; MVZ and PC collected serum samples of breast cancer patients at the University of Ferrara Hospital; JMS and AP collected serum samples from age-matched healthy controls; AH collected serum samples and the associated clinicopathological information at the Mercy Hospital in Oklahoma City; AM, SS, MN, and MF designed the study and critically evaluated the results; and AM and MN wrote the manuscript.
Distribution of miRNA levels in sera of the two combined cohorts of breast cancer and disease-free patients. Figure S2. Validation of ddPCR results using Real-Time PCR. Figure S3. Serum miR-652-3p is significantly reduced in patients affected by either Luminal A or non-Luminal A breast cancers. Table S1. Published data on circulating microRNAs in human breast cancer patients. Table S2. Demographic characteristics of study populations.