- Research article
- Open Access
Serum microRNA expression as an early marker for breast cancer risk in prospectively collected samples from the Sister Study cohort
Breast Cancer Research volume 15, Article number: R42 (2013)
MicroRNAs (miRNAs) are small, non-coding, single-stranded RNAs between 18-22 nucleotides long that regulate gene expression. Expression of miRNAs is altered in tumor compared to normal tissue; there is some evidence that these changes may be reflected in the serum of cancer cases compared to healthy individuals. This has yet to be examined in a prospective study where samples are collected before diagnosis.
We used Affymetrix arrays to examine serum miRNA expression profiles in 410 participants in the Sister Study, a prospective cohort study of 50,884 women. All women in the cohort had never been diagnosed with breast cancer at the time of enrollment. We compared global miRNA expression patterns in 205 women who subsequently developed breast cancer and 205 women who remained breast cancer-free. In addition within the case group we examined the association of miRNA expression in serum with different tumor characteristics, including hormone status (ER, PR, and HER-2) and lymph node status.
Overall, 414 of 1,105 of the human miRNAs on the chip were expressed above background levels in 50 or more women. When the average expression among controls was compared to cases using conditional logistic regression, 21 miRNAs were found to be differentially expressed (P≤.05). Using qRT-PCR on a small, independent sample of 5 cases and 5 controls we verified overexpression of the 3 highest expressing miRNAs among cases, miR-18a, miR-181a, and miR-222; the differences were not statistically significant in this small set. The 21 differentially expressed miRNAs are known to target at least 82 genes; using the gene list for pathway analysis we found enrichment of genes involved in cancer-related processes. In a separate case-case analyses restricted to the 21 miRNAs, we found 7 miRNAs with differential expression for women whose breast tumors differed by HER-2 expression, and 10 miRNAs with differential expression by nodal status.
miRNA levels in serum show a number of small differences between women who later develop cancer versus those who remain cancer-free.
MicroRNAs (miRNAs) are small, non-coding, single-stranded RNAs ranging in size between 18 and 22 nucleotides; they are typically excised from longer, 60- to 110-nucleotide stem-loop precursors [1, 2]. miRNAs are involved in fundamental biological processes, including development, differentiation, apoptosis, and proliferation, and are believed to act predominately as post-transcriptional regulators that can either degrade their mRNA targets or repress their translation . A single miRNA may have multiple mRNA targets, and up to 30% of human genes may be regulated by miRNAs [4, 5].
Aberrant expression of miRNAs in cancer was initially identified in B-cell chronic lymphocytic leukemia , and miRNA dysregulation has been subsequently reported for many tumor types in which, depending on the specific target mRNA(s), they may act either as tumor suppressor genes or as oncogenes [7, 8]. In breast cancer, post-diagnosis miRNA levels have been shown to correlate with a number of tumor characteristics, including stage, vascular invasion, proliferative index, and estrogen receptor/progesterone receptor (ER/PR) status [9, 10], and may have prognostic value.
miRNAs have recently been found in human serum and plasma, where they appear to be resistant to RNAase degradation and thus relatively stable, even in stored samples . This stability has made miRNAs appealing candidates for epidemiologic studies of stored samples, particularly since miRNA profiling requires only small amounts of serum or plasma . The use of circulating miRNA profiles as potential early-detection cancer markers has generated considerable interest [13–16], although data addressing such application remain sparse. Initial studies have suggested that serum levels of miRNAs may differ between diagnosed cancer cases and controls , and several recent case control studies of breast cancer have reported evidence of differential miRNA expression levels in serum [18–21]. These studies have shown little agreement, perhaps because some have measured only a few miRNAs whereas others have used more comprehensive miRNA screens, but with a small number of subjects. None has used samples obtained prior to diagnosis. Use of such prospective samples avoids a number of important potential biases (for example, differential selection and processing of cases and controls or the possibility that the differences observed in case samples are the result of biopsy, cancer treatments, behavioral changes, stress, or other factors experienced by cases but not controls).
Here, we report on a study that prospectively collected serum samples from 205 women who subsequently developed breast cancer and 205 women who remained cancer-free and that used microarrays to comprehensively assess known miRNAs.
Materials and methods
The Sister Study  is a prospective cohort study of 50,884 women and was designed to examine the environmental and genetic determinants of breast cancer. The cohort has been previously described ; briefly, women from the US or Puerto Rico were eligible to enroll if they themselves had never had breast cancer but had a full or half-sister who had breast cancer. At baseline interview, all participants provided extensive information, including family history, reproductive history, and information about potential risk factors. Informed consent and blood samples were obtained during a home visit. For women who subsequently developed breast cancer, detailed information on diagnosis was collected from medical records and self-report. Pathology reports were abstracted for tumor grade, stage, and other information, including status for ER, PR, and HER-2 (human epidermal growth factor receptor 2) expression. The study was approved by the Institutional Review Board of the National Institute of Environmental Health Sciences, National Institutes of Health, and the Copernicus Group Institutional Review Board.
Selection of cases and controls
We designed a matched-pair nested case control study. We selected patients who had confirmed invasive breast cancer, who completed enrollment by August 2008, and whose diagnosis occurred within 18 months following blood draw (n = 242). We excluded 29 cases who lacked a serum sample or whose sample had integrity issues during collection and shipping and eight cases whose sample had limited volume, leaving 205 cases that are the focus of our study. For each case, a matched control was selected from the 50,884 participants on the basis of the following criteria: no history of cancer (other than non-melanoma skin cancer), having completed enrollment by August 2008, an available blood sample, same race (non-Hispanic white, black, Hispanic, or other), similar age at enrollment (within 5 years), and similar date of blood draw (within 2 months). Three replicate serum samples from three women (nine samples in total) who were not participants in the study but who provided blood samples that were collected and processed in the same manner as Sister Study participants were used to provide technical replicates.
Assignment to extraction batches and array chip lot
To minimize possible processing and chip lot effects, samples were assigned to processing batches of seven to nine pairs, and batches had similar distributions of age, race, and date of enrollment. For array hybridization, each batch was assigned to one of two different chip lots ('A' and 'B') in a manner designed to ensure a balance of these same characteristics. The nine replicates (described above) were assigned to the same batch and chip lot. Laboratory personnel were blind to case control status and other phenotype information.
RNA extraction, labeling, and hybridization
Total RNA was extracted in batches by using a Total RNA purification kit (cat. no. 17200; Norgen Biotek Corp., Thorold, ON, Canada). In accordance with the manufacturer's recommendation not to exceed 200 µL per column, 400 µL of total serum from each individual was split into two equal 200-µL aliquots and then processed separately following the manufacturer's recommended protocol for total RNA purification from serum. An on-column DNase digestion was added before sample elution by using an RNase-Free DNase I Kit (cat. no. 25710; Norgen Biotek Corp.), and the two aliquots were subsequently pooled. Fixed volumes rather than fixed amounts of RNA were used in accordance with other studies .
Total RNA (8 µL) was directly labeled by using Flash Tag Biotin HSR Labeling kits (cat. no. HSR30FTA; Genisphere, LLC, Hatfield, PA, USA) in accordance with the instructions of the manufacturer. RNA was heated to 80°C for 10 minutes before labeling to inactivate any residual DNase activity. RNA was hybridized for 42 hours to the GeneChip miRNA 2.0 array (cat. no. 901755; Affymetrix Inc., Santa Clara, CA, USA ). The GeneChip miRNA 2.0 arrays contain 100% miRBase version 15 coverage of 131 organisms and contain probes for 3,439 human non-coding RNAs (ncRNAs), including 1,105 miRNAs and 2,334 other ncRNAs (including scaRNAs and snoRNAs). The arrays were washed and stained by using standard Affymetrix protocols and scanned by using an Affymetrix GCS 3000 7G Scanner. Feature intensities were extracted by using miRNA 2.0 array library files. Array hybridization and scanning were completed by Precision Biomarker Resources, Inc. (Evanston, IL, USA). The average Spearman correlation coefficient values for three sets of three technical replicates were all above 0.8 (Additional file 1). Array data were deposited into the NCBI Gene Expression Omnibus (GSE44281).
Replication samples and qRT-PCR
An independent set of 10 women were used to validate selected miRNAs via quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Five women who provided consent and blood samples but who developed breast cancer prior to completing enrollment were selected as cases, along with five controls who also provided consent and blood samples and who were cancer-free but did not complete enrollment. Total RNA was extracted from serum samples of these women as described above with the addition of Synthetic C. elegans miScript miRNA Mimic (cat. no. MSY0000010; Qiagen, Valencia, CA, USA). Synthetic cel-39 was spiked-in at a final concentration of 0.25 fmol/µL prior to extraction and used as a PCR normalization control. The RNA concentration, reverse transcription, and pre-PCR steps were carried out in accordance with a previously published protocol . ExoSAP-IT (cat. no. 78250; Affymetrix Inc.) treatment followed by column purification (cat. no. 28004; Qiagen) in accordance with the protocol of the manufacturer was used to purify the pre-PCR product. Individual PCR was run in triplicate by using 1 µL of purified pre-PCR product. The reaction contained the following components: 2x Taqman universal master mix (cat. no. 4324018; ABI, Carlsbad, CA, USA), 1 µM forward primer, 1 µM universal reverse primer, and 0.2 µM probe. The reaction was run on a Bio-Rad CFX 384 Real-Time System (Bio-Rad Laboratories, Inc., Hercules, CA, USA) by using the following parameters: 55°C for 2 minutes, 95°C for 10 minutes, followed by 40 cycles of 95°C for 15 seconds and 55°C for 1 minute. PCR cycle threshold (Ct) values were recorded for each target gene and for normalization controls and were averaged across three independent runs. Primers for miR-222, miR-181a, miR-1825, and miR-18a were custom-ordered from IDT (San Diego, CA, USA) by using previously published sequences . Primers for cel-39 were designed in the same fashion as above and custom-ordered from IDT.
To determine the best candidate miRNA for PCR normalization in our data set, we ran the array expression data from the 47 miRNAs expressed in almost all individuals through the NormFinder software . NormFinder uses a model-based variance estimation approach . Using these results, we selected as a qRT-PCR normalization control miR-1825, which showed one of the highest stability values across the 410 cases and controls and had blood levels that were similar to those of the three target miRNAs. We used the average of miR-1825 and an external spike-in cel-39 control, a strategy shown to be effective for controlling both technical and biologic variability in qRT-PCR assays from serum [17, 24]. The efficiency of the four PCR assays (for miR-181a, miR-18a, miR-222, and miR-1825) was similar for all four assays (Additional file 2). Normalized relative expression was based on Ct values and calculated as 1/(Ctgene−Ctnorm).
Data processing and statistical analysis
miRNA expression intensity values were background-corrected and normalized across arrays by using the robust multichip average method . The intensity data used in all analysis were log (2)-transformed.
For each array, the miRNA probe set signals were compared with the distribution of signals for anti-genomic probes that had matching GC content (miRNA QC Tool, version 188.8.131.52), and in accordance with the recommendation of the manufacturer, Wilcoxon rank-sum test of P value of less than 0.06 was used to identify miRNAs above background. Subsequent analysis was restricted to 414 miRNAs that exceeded background levels in at least 50 women. Conditional logistic regression was used to identify differentially expressed miRNA probes between cases and controls for those 414 probes. Because analysis of circulating miRNAs in prospectively collected samples is still exploratory, we - like some other investigators of circulating miRNAs [30, 31] - regard these results as descriptive and not as tests of hypotheses and so provide P values that are unadjusted for multiple comparisons.
The association between miRNAs and the tumor characteristics of hormone receptor status (ER, PR, and HER-2) and lymph node status was tested in a case-only logistic analysis, in which race was adjusted for. Chip lot and batch were specified as random effect variables. All statistical analyses were performed by using R 2.15.
Pathway analysis with ingenuity pathway analysis
miRNAs found to be significantly associated with case control status were further analyzed with ingenuity pathway analysis (IPA) . Using IPA's microRNA target filter, we generated a list of predicted mRNA targets for each of the 21 significant miRNAs. The list was then restricted to the mRNAs listed in the IPA database as experimentally verified targets of any of the 21 miRNAs. This mRNA target list was then used to run a canonical pathway analysis.
A large number of miRNAs are detected in serum
In total, 410 serum samples from breast cancer cases (n = 205) and controls (n = 205) were analyzed in this study; baseline characteristics of the cases and controls are summarized in Table 1. Of the 1,105 human miRNAs, 414 miRNAs were detected above background threshold levels in at least 50 women. Forty-seven miRNAs were detected above background in 400 or more women (Table 2), and miR-16 showed the highest average expression. Even though expression of miRNAs showed considerable inter-individual variation, several miRNAs, including miR-1825 and miR-1228, were relatively constant among women (Figure 1).
Discovery of differentially expressed miRNAs in serum
When paired case control analysis of the 414 miRNAs expressed above background was used, 21 miRNAs showed significantly different levels in cases and controls (P ≤0.05) (Table 3). The differences were small, ranging from 4% to 19%. Higher miRNA expression in women destined to become cases was significantly more common (16 of 21 miRNAs) than would be expected by chance alone (binomial test, two-tailed P <0.05). Differential miRNA expression was not stronger in women close to their time of diagnosis, but sample size was small and all cases were diagnosed within 18 months of blood draw (data not shown). Using qRT-PCR on a small independent replication set of five cases and five controls, we further examined the three miRNAs (miR-18a, miR-181a, and miR-222) with the highest expression in cases. As predicted, all three miRNAs showed higher levels in cases, although none was statistically significant in this small set of women (Additional file 3).
The impact of miRNA alterations on regulatory pathways
To explore potential biological associations, we ran IPA on the 82 experimentally verified mRNA targets of the 21 differentially expressed miRNAs. Sixteen IPA canonical pathways, including molecular mechanisms of cancer, were enriched as were other cancer-related pathways, including p53 signaling, cyclins and cell cycle regulation, and Myc-mediated apoptosis signaling (Additional file 4).
miRNA expression association with tumor characteristics
To investigate the potential association of serum miRNA expression with tumor characteristics in the 205 women who later developed breast cancer, we subclassified them into groups based on tumor characteristics (Table 4) and performed a case-case comparison. There was no evidence of significant differences in serum miRNA levels based on tumor ER or PR staining characteristics. In comparisons of serum samples from the 25 women who developed HER-2-positive tumors with 147 samples from women who developed HER-2-negative tumors, there were seven miRNAs with significantly differential expression (P ≤0.05); one miRNA was overexpressed and six miRNAs were underexpressed in the HER-2-positive tumors (Figure 2A and Additional file 5). Case-case comparison of serum from women who subsequently developed lymph node-negative tumors (pN0, n = 153) with that of women who developed lymph node-positive tumors (pN1, pN2, or pN3, n = 52) revealed 10 differentially expressed miRNAs (P ≤0.05); five were overexpressed and five were underexpressed in node-positive tumors (Figure 2B and Additional file 6).
miRNA profiles are gaining interest as potential diagnostic or prognostic markers for breast cancer . However, existing studies have been limited by sample size or the number of miRNAs analyzed, and none has used prospectively collected samples [18, 31, 34]. Our study minimized potential biases by profiling global serum miRNA expression patterns in samples obtained from women prior to clinical diagnosis (mean time to diagnosis was 10 months). We found a set of 21 miRNAs differentially expressed in serum samples from 205 women who subsequently developed breast cancer compared with 205 women who remained cancer-free during the time of follow-up. The differences in miRNA levels were small and include both overexpression and underexpression of miRNAs in the cases, and overexpression was significantly more frequent than would be expected had the association been random. Published reports of primary breast tumors or cell lines have examined seven of the 21 differentially expressed miRNAs we found, and all seven showed agreement with the direction of change in our case serum samples (Table 3). IPA of the mRNA targets of these differentially expressed miRNAs suggested gene enrichment for cancer-related signaling pathways. Although the absolute differences in miRNA levels between serum samples of cases and controls are quite small, differences pre-date clinical diagnosis and may reflect important pathways for breast cancer development.
miR-18a, miR-181a, and miR-222 showed the highest percentage difference between cases and controls in our study; qRT-PCR of these miRNAs in a small independent replication set of cases and controls, though not statistically significant, replicated the direction of change for all three. These three miRNAs have been suggested to act as oncogenes through regulation of their potential target mRNAs. miR-18a is part of the oncogenic miR 17-92 cluster, which is often overexpressed in solid tumors, including breast . Overexpression of this cluster is believed to cooperate with c-Myc in stimulating proliferation by negatively regulating E2F1 [36, 37]. Increased expression of miR-181a in the bone marrow of patients with breast cancer has been reported to be associated with shorter disease-free survival, higher grade, and breast cancer recurrence . miR-181a is believed to target the tumor suppressor gene programmed cell death protein 4 (PDCD4) , which inhibits tumor neoplastic transformation . In breast cancer cell lines, miR-222 overexpression has been reported to be associated with tamoxifen resistance through targeting the cell cycle inhibitor p27 (Kip1) . miR-222 has also been reported to increase proliferation of ERα-negative cells while reducing the expression of various tumor suppressor proteins , and expression of miR-222 has been reported to increase cell migration in the epithelial-to-mesenchymal transition acting downstream of the RAS-RAF-MEK oncogenic pathway .
Interestingly, two recent case control studies have provided evidence that both miR-222 and miR-181a are overexpressed in the serum of patients with breast cancer. One used sequencing by oligonucleotide ligation and detection (SOLiD) of serum samples obtained prior to surgery from 13 breast cancer cases compared with samples from 10 healthy controls and found 26 miRNAs that were overexpressed in cases, including miR-222 and miR-181a; overexpression of miR-222 was validated in an independent group of 50 cases and 50 controls by using qRT-PCR . A second study used Solexa sequencing combined with Taqman low-density array chips on serum samples obtained prior to surgery from 48 breast cancer cases and 48 controls; 10 miRNAs were found to be overexpressed in the cases, and four were validated by using qRT-PCR in an independent group of 76 cases and 76 controls . That study also found overexpression of miR-222 . These studies, combined with our prospective study, provide a growing body of evidence that miR-222 measured in blood is associated with breast cancer.
Among cases, we compared the serum miRNA profiles of women with different tumor characteristics, including hormone status (ER, PR, and HER-2) and nodal status. Although there were no differences in ER or PR status, there were differences in HER-2 and lymph node status. Of the seven miRNAs differentially expressed in the serum of women who developed HER-2-overexpressing breast tumors, miR-93, miR-183, and miR-29a have been reported to be associated with breast cancer in previous studies [20, 43, 44]. In our study, miR-93 was underexpressed in the serum of women who developed HER-2-overexpressing breast tumors; interestingly, miR-93 expression was recently shown to induce a more differentiated cell phenotype in breast cancer cell lines, and expression of miR-93 in mouse mammary fat pads blocked tumor development and metastases . Of the 10 miRNAs differentially expressed in the serum of women with tumors that spread to the lymph nodes (pN1 or higher), four (miR-145, miR-124, miR-125b, and miR-320) have been reported to be associated with breast cancer in previous studies [45–48]. Of these, miR-320 is of particular interest as we found three miR-320 family members (miR-320b, miR-320d, and miR320e) to be underexpressed in the serum of women who developed lymph node-positive breast tumors. miR-320 has been reported to be decreased in breast tumor tissue and downregulation of miR-320 - through loss of phosphatase and tensin homolog (PTEN) -has been shown to promote tumor proliferation and invasiveness in mouse models; expression of miR-320 distinguished human normal breast stroma from tumor stroma and was correlated with recurrence . A study comparing miRNA expression in inflammatory breast cancer (IBC) with non-IBC also found miR-320 to be downregulated in the more aggressive IBC group of tumors . Thus, loss of miR-320 expression may be associated with a higher likelihood of lymph node involvement and a more aggressive metastatic phenotype.
Although miRNAs that are differentially expressed between tumor and normal tissue are more frequently downregulated in tumor tissue , our study (like others [20, 21]) has found that circulating miRNAs that differ in levels between breast cancer patients and controls are more frequently at higher levels in case blood samples. The mechanism underlying circulating miRNA stability is still being investigated. One model involves the active release of miRNAs from cells in membrane-bound microvesicles, including exosomes and shedding vesicles [51–53]. There is evidence that microvesicles can deliver miRNAs to recipient cells and trigger changes in target mRNA levels . A recent report has shown that vesicle-encapsulated miRNAs represent only a minor portion of circulating miRNAs but that a significant portion of circulating miRNAs are associated with Argonaute2 (Ago2) , the effector component of the miRNA-induced silencing complex . Both models support the possibility that miRNAs may be actively released into circulation and could act as signaling molecules able to regulate their target mRNA expression in recipient cells. Cancer-associated miRNAs in the circulation could also originate from immunocytes in the tumor microenvironment or from a response mediated by the body's systemic response to disease .
We find some evidence of differences in miRNA serum levels between women who subsequently developed cancer compared with women who remained clinically cancer-free. The magnitude of these differences is small, and this may limit their clinical application as circulating early-detection markers for breast cancer. This is the first study to use prospectively collected samples; limitations include a relatively short follow-up period and a sample size that is not large enough to fully explore etiologic versus diagnostic relevance. Our study was carried out within a cohort of women who each had a sister with breast cancer, putting the former at about twofold increased risk, and so the differences that we observed may not be generalizable to women without a similar family history.
human epidermal growth factor receptor 2
inflammatory breast cancer
ingenuity pathway analysis
polymerase chain reaction
quantitative reverse transcription-polymerase chain reaction.
Ambros V: The functions of animal microRNAs. Nature. 2004, 431: 350-355.
Pasquinelli AE, Hunter S, Bracht J: MicroRNAs: a developing story. Curr Opin Genet Dev. 2005, 15: 200-205.
Sayed D, Abdellatif M: MicroRNAs in development and disease. Physiol Rev. 2011, 91: 827-887.
Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005, 433: 769-773.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297.
Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, Aldler H, Rattan S, Keating M, Rai K, Rassenti L, Kipps T, Negrini M, Bullrich F, Croce CM: Frequent deletions and down-regulation of micro-RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci USA. 2002, 99: 15524-15529.
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz HR, Golub TR: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838.
Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, Prueitt RL, Yanaihara N, Lanza G, Scarpa A, Vecchione A, Negrini M, Harris CC, Croce CM: A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci USA. 2006, 103: 2257-2261.
Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M, Menard S, Palazzo JP, Rosenberg A, Musiani P, Volinia S, Nenci I, Calin GA, Querzoli P, Negrini M, Croce CM: MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005, 65: 7065-7070.
Mattie MD, Benz CC, Bowers J, Sensinger K, Wong L, Scott GK, Fedele V, Ginzinger D, Getts R, Haqq C: Optimized high-throughput microRNA expression profiling provides novel biomarker assessment of clinical prostate and breast cancer biopsies. Mol Cancer. 2006, 5: 24-
Gilad S, Meiri E, Yogev Y, Benjamin S, Lebanony D, Yerushalmi N, Benjamin H, Kushnir M, Cholakh H, Melamed N, Bentwich Z, Hod M, Goren Y, Chajut A: Serum microRNAs are promising novel biomarkers. PLoS One. 2008, 3: e3148-
Lodes MJ, Caraballo M, Suciu D, Munro S, Kumar A, Anderson B: Detection of cancer with serum miRNAs on an oligonucleotide microarray. PLoS One. 2009, 4: e6229-
Weiland M, Gao XH, Zhou L, Mi QS: Small RNAs have a large impact: circulating microRNAs as biomarkers for human diseases. RNA Biol. 2012, 9: 850-859.
Heneghan HM, Miller N, Kerin MJ: MiRNAs as biomarkers and therapeutic targets in cancer. Curr Opin Pharmacol. 2010, 10: 543-550.
Cortez MA, Calin GA: MicroRNA identification in plasma and serum: a new tool to diagnose and monitor diseases. Expert Opin Biol Ther. 2009, 9: 703-711.
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-2092.
Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, Peterson A, Noteboom J, O'Briant KC, Allen A, Lin DW, Urban N, Drescher CW, Knudsen BS, Stirewalt DL, Gentleman R, Vessella RL, Nelson PS, Martin DB, Tewari M: Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci USA. 2008, 105: 10513-10518.
Asaga S, Kuo C, Nguyen T, Terpenning M, Giuliano AE, Hoon DS: Direct serum assay for microRNA-21 concentrations in early and advanced breast cancer. Clin Chem. 2011, 57: 84-91.
van Schooneveld E, Wouters MC, Van der Auwera I, Peeters DJ, Wildiers H, Van Dam PA, Vergote I, Vermeulen PB, Dirix LY, Van Laere SJ: Expression profiling of cancerous and normal breast tissues identifies microRNAs that are differentially expressed in serum from patients with (metastatic) breast cancer and healthy volunteers. Breast Cancer Res. 2012, 14: R34-
Wu Q, Wang C, Lu Z, Guo L, Ge Q: Analysis of serum genome-wide microRNAs for breast cancer detection. Clin Chim Acta. 2012, 413: 1058-1065.
Hu Z, Dong J, Wang LE, Ma H, Liu J, Zhao Y, Tang J, Chen X, Dai J, Wei Q, Zhang C, Shen H: Serum microRNA profiling and breast cancer risk: the use of miR-484/191 as endogenous controls. Carcinogenesis. 2012, 33: 828-834.
The Sister Study homepage. [http://www.sisterstudy.org]
D'Aloisio AA, Baird DD, DeRoo LA, Sandler DP: Association of intrauterine and early-life exposures with diagnosis of uterine leiomyomata by 35 years of age in the Sister Study. Environ Health Perspect. 2010, 118: 375-381.
Kroh EM, Parkin RK, Mitchell PS, Tewari M: Analysis of circulating microRNA biomarkers in plasma and serum using quantitative reverse transcription-PCR (qRT-PCR). Methods. 2010, 50: 298-301.
Affymetrix, Inc. homepage. [http://www.affymetrix.com]
Moltzahn F, Hunkapiller N, Mir AA, Imbar T, Blelloch R: High throughput microRNA profiling: optimized multiplex qRT-PCR at nanoliter scale on the fluidigm dynamic arrayTM IFCs. J Vis Exp. 2011, 2552-54
Molecular Diagnostic Laboratory homepage. [http://www.mdl.dk/publicationsnormfinder.htm]
Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64: 5245-5250.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264.
Chen X, Hu Z, Wang W, Ba Y, Ma L, Zhang C, Wang C, Ren Z, Zhao Y, Wu S, Zhuang R, Zhang Y, Hu H, Liu C, Xu L, Wang J, Shen H, Zhang J, Zen K, Zhang CY: Identification of ten serum microRNAs from a genome-wide serum microRNA expression profile as novel noninvasive biomarkers for nonsmall cell lung cancer diagnosis. Int J Cancer. 2012, 130: 1620-1628.
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-
Ingenuity Systems, Inc. homepage. [http://www.ingenuity.com/]
Fu SW, Chen L, Man YG: miRNA biomarkers in breast cancer detection and management. J Cancer. 2011, 2: 116-122.
Heneghan HM, Miller N, Lowery AJ, Sweeney KJ, Newell J, Kerin MJ: Circulating microRNAs as novel minimally invasive biomarkers for breast cancer. Ann Surg. 2010, 251: 499-505.
Song L, Lin C, Wu Z, Gong H, Zeng Y, Wu J, Li M, Li J: miR-18a impairs DNA damage response through downregulation of ataxia telangiectasia mutated (ATM) kinase. PLoS One. 2011, 6: e25454-
He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ, Hammond SM: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435: 828-833.
O'Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005, 435: 839-843.
Ota D, Mimori K, Yokobori T, Iwatsuki M, Kataoka A, Masuda N, Ishii H, Ohno S, Mori M: Identification of recurrence-related microRNAs in the bone marrow of breast cancer patients. Int J Oncol. 2011, 38: 955-962.
Cmarik JL, Min H, Hegamyer G, Zhan S, Kulesz-Martin M, Yoshinaga H, Matsuhashi S, Colburn NH: Differentially expressed protein Pdcd4 inhibits tumor promoter-induced neoplastic transformation. Proc Natl Acad Sci USA. 1999, 96: 14037-14042.
Miller TE, Ghoshal K, Ramaswamy B, Roy S, Datta J, Shapiro CL, Jacob S, Majumder S: MicroRNA-221/222 confers tamoxifen resistance in breast cancer by targeting p27Kip1. J Biol Chem. 2008, 283: 29897-29903.
Di Leva G, Gasparini P, Piovan C, Ngankeu A, Garofalo M, Taccioli C, Iorio MV, Li M, Volinia S, Alder H, Nakamura T, Nuovo G, Liu Y, Nephew KP, Croce CM: MicroRNA cluster 221-222 and estrogen receptor alpha interactions in breast cancer. J Natl Cancer Inst. 2010, 102: 706-721.
Stinson S, Lackner MR, Adai AT, Yu N, Kim HJ, O'Brien C, Spoerke J, Jhunjhunwala S, Boyd Z, Januario T, Newman RJ, Yue P, Bourgon R, Modrusan Z, Stern HM, Warming S, de Sauvage FJ, Amler L, Yeh RF, Dornan D: TRPS1 targeting by miR-221/222 promotes the epithelial-to-mesenchymal transition in breast cancer. Sci Signal. 2011, 4: ra41-
Shimono Y, Zabala M, Cho RW, Lobo N, Dalerba P, Qian D, Diehn M, Liu H, Panula SP, Chiao E, Dirbas FM, Somlo G, Pera RA, Lao K, Clarke MF: Downregulation of miRNA-200c links breast cancer stem cells with normal stem cells. Cell. 2009, 138: 592-603.
Liu S, Patel SH, Ginestier C, Ibarra I, Martin-Trevino R, Bai S, McDermott SP, Shang L, Ke J, Ou SJ, Heath A, Zhang KJ, Korkaya H, Clouthier SG, Charafe-Jauffret E, Birnbaum D, Hannon GJ, Wicha MS: MicroRNA93 regulates proliferation and differentiation of normal and malignant breast stem cells. PLoS Genet. 2012, 8: e1002751-
Wu X, Somlo G, Yu Y, Palomares MR, Li AX, Zhou W, Chow A, Yen Y, Rossi JJ, Gao H, Wang J, Yuan YC, Frankel P, Li S, Ashing-Giwa KT, Sun G, Wang Y, Smith R, Robinson K, Ren X, Wang SE: De novo sequencing of circulating miRNAs identifies novel markers predicting clinical outcome of locally advanced breast cancer. J Transl Med. 2012, 10: 42-
Zhou M, Liu Z, Zhao Y, Ding Y, Liu H, Xi Y, Xiong W, Li G, Lu J, Fodstad O, Riker AI, Tan M: MicroRNA-125b confers the resistance of breast cancer cells to paclitaxel through suppression of pro-apoptotic Bcl-2 antagonist killer 1 (Bak1) expression. J Biol Chem. 2010, 285: 21496-21507.
Lv XB, Jiao Y, Qing Y, Hu H, Cui X, Lin T, Song E, Yu F: miR-124 suppresses multiple steps of breast cancer metastasis by targeting a cohort of pro-metastatic genes in vitro. Chin J Cancer. 2011, 30: 821-830.
Yan LX, Huang XF, Shao Q, Huang MY, Deng L, Wu QL, Zeng YX, Shao JY: MicroRNA miR-21 overexpression in human breast cancer is associated with advanced clinical stage, lymph node metastasis and patient poor prognosis. RNA. 2008, 14: 2348-2360.
Bronisz A, Godlewski J, Wallace JA, Merchant AS, Nowicki MO, Mathsyaraja H, Srinivasan R, Trimboli AJ, Martin CK, Li F, Yu L, Fernandez SA, Pecot T, Rosol TJ, Cory S, Hallett M, Park M, Piper MG, Marsh CB, Yee LD, Jimenez RE, Nuovo G, Lawler SE, Chiocca EA, Leone G, Ostrowski MC: Reprogramming of the tumour microenvironment by stromal PTEN-regulated miR-320. Nat Cell Biol. 2012, 14: 159-167.
Van der Auwera I, Limame R, van Dam P, Vermeulen PB, Dirix LY, Van Laere SJ: Integrated miRNA and mRNA expression profiling of the inflammatory breast cancer subtype. Br J Cancer. 2010, 103: 532-541.
Hunter MP, Ismail N, Zhang X, Aguda BD, Lee EJ, Yu L, Xiao T, Schafer J, Lee ML, Schmittgen TD, Nana-Sinkam SP, Jarjoura D, Marsh CB: Detection of microRNA expression in human peripheral blood microvesicles. PLoS One. 2008, 3: e3694-
Kosaka N, Iguchi H, Yoshioka Y, Takeshita F, Matsuki Y, Ochiya T: Secretory mechanisms and intercellular transfer of microRNAs in living cells. J Biol Chem. 2010, 285: 17442-17452.
Simons M, Raposo G: Exosomes--vesicular carriers for intercellular communication. Curr Opin Cell Biol. 2009, 21: 575-581.
Valadi H, Ekstrom K, Bossios A, Sjostrand M, Lee JJ, Lotvall JO: Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007, 9: 654-659.
Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, Mitchell PS, Bennett CF, Pogosova-Agadjanyan EL, Stirewalt DL, Tait JF, Tewari M: Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci USA. 2011, 108: 5003-5008.
Meister G, Landthaler M, Patkaniowska A, Dorsett Y, Teng G, Tuschl T: Human Argonaute2 mediates RNA cleavage targeted by miRNAs and siRNAs. Mol Cell. 2004, 15: 185-197.
Ma R, Jiang T, Kang X: Circulating microRNAs in cancer: origin, function and application. J Exp Clin Cancer Res. 2012, 31: 38-
Zhao JJ, Lin J, Yang H, Kong W, He L, Ma X, Coppola D, Cheng JQ: MicroRNA-221/222 negatively regulates estrogen receptor alpha and is associated with tamoxifen resistance in breast cancer. J Biol Chem. 2008, 283: 31079-31086.
Guttilla IK, Phoenix KN, Hong X, Tirnauer JS, Claffey KP, White BA: Prolonged mammosphere culture of MCF-7 cells induces an EMT and repression of the estrogen receptor by microRNAs. Breast Cancer Res Treat. 2012, 132: 75-85.
Zhang L, Huang J, Yang N, Greshock J, Megraw MS, Giannakakis A, Liang S, Naylor TL, Barchetti A, Ward MR, Yao G, Medina A, O'Brien-Jenkins A, Katsaros D, Hatzigeorgiou A, Gimotty PA, Weber BL, Coukos G: microRNAs exhibit high frequency genomic alterations in human cancer. Proc Natl Acad Sci USA. 2006, 103: 9136-9141.
Navon R, Wang H, Steinfeld I, Tsalenko A, Ben-Dor A, Yakhini Z: Novel rank-based statistical methods reveal microRNAs with differential expression in multiple cancer types. PLoS One. 2009, 4: e8003-
Leivonen SK, Makela R, Ostling P, Kohonen P, Haapa-Paananen S, Kleivi K, Enerly E, Aakula A, Hellstrom K, Sahlberg N, Kristensen VN, Borresen-Dale AL, Saviranta P, Perala M, Kallioniemi O: Protein lysate microarray analysis to identify microRNAs regulating estrogen receptor signaling in breast cancer cell lines. Oncogene. 2009, 28: 3926-3936.
Yoshimoto N, Toyama T, Takahashi S, Sugiura H, Endo Y, Iwasa M, Fujii Y, Yamashita H: Distinct expressions of microRNAs that directly target estrogen receptor alpha in human breast cancer. Breast Cancer Res Treat. 2011, 130: 331-339.
Zhang H, Su SB, Zhou QM, Lu YY: [Differential expression profiles of microRNAs between breast cancer cells and mammary epithelial cells]. Ai Zheng. 2009, 28: 493-499. Article in Chinese
Manavalan TT, Teng Y, Appana SN, Datta S, Kalbfleisch TS, Li Y, Klinge CM: Differential expression of microRNA expression in tamoxifen-sensitive MCF-7 versus tamoxifen-resistant LY2 human breast cancer cells. Cancer Lett. 2011, 313: 26-43.
Sobin LH, Gospodarowicz MK, Wittekind C: Breast tumours. TNM Online. 2003, Hoboken, NJ: John Wiley & Sons
We thank Dan Scharf at Westat (Rockville, MD, USA) for sample selection and batch and chip lot designation of the samples and Carol Messler at SSS for pulling and delivering the serum samples. We are especially grateful to Robert (Bob) Getts and Jessica Bowers at Genisphere and Keith Jones and Radharani Duttagupta at Affymetrix for their help and advice in pilot studies and troubleshooting, and we want to thank Bob Getts and Jessica Bowers for doing all of the sample labeling. All of the array hybridization and scanning was carried out by Precision Biomarkers Resources, and we would especially like to thank both David Paul and Jason Monroe at Precision Biomarkers Resources for their help with this project. This work was funded by the Intramural Research Program of the National Institutes of Health, National Institute of Environmental Health Sciences.
The authors declare that they have no competing interests.
ACG carried out sample extraction and sample preparation, did some statistical analysis, and helped to design the experiments, interpret the results, and write the manuscript. ZX performed the majority of the statistical analysis. CRW and PAW participated in study design. LDR participated in study design and provided patient data and study variables. DPS participated in study design and collected samples and data. RCG provided study reagents and processing as well as critical advice on study design, QA, and data interpretation. JAT helped to design the experiments, interpret the results, and write the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Spearman correlation coefficient values for technical replicates of arrays. Three replicate serum samples from three women (nine samples in total) were processed and hybridized to arrays as described for samples in the main study. Spearman correlation coefficients were calculated for the three pairings of replicate samples for each woman and averaged. One array from Individual 1 appeared to be an outlier but was included in the results shown above. Exclusion of this array resulted in correlation coefficients of greater than 0.97 in all three categories of probes for Individual 1. (PDF 88 KB)
Additional file 2: Efficiency of the four polymerase chain reaction (PCR) assays. The efficiency of PCR amplifications for the normalization control and three target microRNAs (miRNAs) was calculated by using DART-PCR version 1.0. The average efficiency of each of three independent PCRs was similar across all four miRNAs. (PDF 132 KB)
Additional file 3: Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) validation in 10 cases and controls. Serum levels of five cases and five controls were examined by using qRT-PCR for miR181a, miR18a, and miR-222. Box plots show normalized relative expression. Normalization was carried out by using the mean of miR-1825 and spiked in cel-39. The horizontal line represents the mean for each sample. (PDF 80 KB)
Additional file 4: . Identified ingenuity pathway analysis (IPA) canonical pathways enriched by the experimentally observed targets of the 21 miRNAs differentially expressed between the cases and non-cases. The negative log (10) false discovery rate-corrected P values are shown. Note that this test has not corrected for possible dependencies across the mRNAs considered and that statistical significance may be overstated. HER-2, human epidermal growth factor receptor 2; ILK, integrin-linked kinase; PTEN, phosphatase and tensin homolog. (PDF 87 KB)
Additional file 5: List of differentially expressed serum microRNAs (miRNAs) based on HER-2/NEU expression. (XLSX 11 KB)
About this article
Cite this article
Godfrey, A.C., Xu, Z., Weinberg, C.R. et al. Serum microRNA expression as an early marker for breast cancer risk in prospectively collected samples from the Sister Study cohort. Breast Cancer Res 15, R42 (2013). https://doi.org/10.1186/bcr3428
- Breast Cancer
- miRNA Expression
- Ingenuity Pathway Analysis
- Inflammatory Breast Cancer
- Develop Breast Cancer