Skip to main content
  • Research article
  • Open access
  • Published:

Circulating miRNAs miR-34a and miR-150 associated with colorectal cancer progression

Abstract

Background

Screening for the early detection of colorectal cancer is important to improve patient survival. The aim of this study was to investigate the potential of circulating cell-free miRNAs as biomarkers of CRC, and their efficiency at delineating patients with polyps and benign adenomas from normal and cancer patient groups.

Methods

The expression of 667 miRNAs was assessed in a discovery set of 48 plasma samples comprising normal, polyp, adenoma, and early and advanced cancer samples. Three miRNAs (miR-34a, miR-150, and miR-923) were further examined in a validation cohort of 97 subjects divided into the same five groups, and in an independent public dataset of 40 CRC samples and paired normal tissues.

Results

High levels of circulating miR-34a and low miR-150 levels distinguished groups of patients with polyps from those with advanced cancer (AUC = 0.904), and low circulating miR-150 levels separated patients with adenomas from those with advanced cancer (AUC = 0.875). In addition, the altered expression of miR-34a and miR-150 in an independent public dataset of forty CRC samples and paired normal tissues was confirmed.

Conclusion

We identified two circulating miRNAs capable of distinguishing patient groups with different diseases of the colon from each other, and patients with advanced cancer from benign disease groups.

Peer Review reports

Background

Colorectal cancer (CRC) poses a significant threat to the health of global populations; it is the second most commonly diagnosed cancer in females and the third in males [1]. CRC develops in a progressive fashion during which normal colon epithelial cells transform to form benign growths such as polyps. These polyps may then progress to benign adenomas, and ultimately to invasive cancer lesions. The progression of the cancer has also been associated with sequential genetic changes in genes such as K-RAS, APC, DCC, and P53 [2]. However CRC is a heterogeneous disease with various patient-related confounding factors such as the anatomic location of the tumour, race/ethnicity of the patient, and genetic and dietary interactions influencing the development of the disease [3].

Screening at risk populations for CRC has significantly improved the outcome for patients, for instance diagnosis while the disease remains localised to the colon dramatically improves patient survival, and removal of early lesions such as adenomatous polyps may prevent disease formation [4]. There are currently several potential screening tests available to detect CRC including the faecal occult blood test (FOBT), flexible sigmoidoscopy (FS), optical colonoscopy (OC) and computed tomography colonography (CTC). FOBT is a simple, cheap and safe test that relies on the assumption that large adenomas and cancerous lesions may bleed, and that these blood products are detectable in the faecal matter of patients. Although cheap and non-invasive, this test is vulnerable to false positive and negative results due to incorrect sample storage, or confounding medical complaints such as haemorrhoids. The other examinations involve more costly and invasive procedures which although allow direct access to colorectal lesions also suffer from low patient acceptance and procedural risks such as perforation of the colon [4].

The focus of the scientific community has thus shifted to exploring the identification of non-invasive biomarkers of disease from bio-fluids such as saliva, urine, and blood. MicroRNAs (miRNAs) are nucleic acid markers that have been recently investigated in this context. MiRNAs are short (20-22nt) non-coding RNAs that negatively regulate gene expression through either mRNA degradation or translational repression [5]. MiRNA expression has been shown to be altered in cancerous tissue compared to normal tissue and different miRNAs have been attributed oncogenic and tumour suppressor qualities [6]. In 2008, Chen et al. detected miRNAs in the serum and plasma blood components of humans and other animals. This primary study illustrated that miRNAs remain stable in serum after being subject to severe conditions such as extremely low or high pH, 10 freeze-thaw cycles, extended storage, boiling, and RNase digestion [7]. In addition to their presence in serum and plasma, miRNAs have also been detected in other body fluids such as urine, saliva, and amniotic fluid making them ideal potential candidates as non-invasive biomarkers of disease [8].

Expression levels of circulating miRNAs have shown some potential at distinguishing cancer patients and healthy controls for prostate [9], ovarian [10], lung [11,12], and breast cancers [13]. Several studies have also investigated circulating miRNA levels for the detection of CRC. Initial approaches analysed small numbers of circulating miRNAs in CRC patient samples compared to normal controls [14]. Other groups performed miRNA profiling on pooled plasma samples and validated candidate biomarkers on additional individual samples [15], and others performed profiling on a small number of CRC tissue/serum/plasma samples before validation in a larger sample set [16]. These studies have produced conflicting results [17] and so recently, groups have begun to perform profiling on larger sample sets and included plasma from patients with adenomas in addition to CRC to improve the specificity of disease detection [18].

In 2008, a guideline was released from the American Cancer Society which highlighted the importance for patients to have access to screening tests that will facilitate cancer prevention through the early detection of cancer, and the detection and removal of polyps [4]. A clear deficit in the search for circulating biomarkers for the early detection of CRC to date is the lack of adenomatous polyp samples and the lack of separation of advanced and early stage cancers represented in studies [14-16,18-20]. The aim of this study was therefore to investigate the potential of circulating cell-free miRNAs not only as biomarkers of CRC, but also their efficiency at delineating patients presenting with polyps and benign adenomas from normal and cancer groups. To facilitate this we performed miRNA profiling for 667 miRNAs on a discovery set of 48 plasma samples comprising 8 normal, 8 polyp, 16 adenoma samples, 8 early stage cancer samples (stage I/II), and 8 advanced cancer samples (stage III/IV). Three candidate miRNAs; miR-34a, miR-150, and miR-923 were then further examined in a validation cohort of 97 independent plasma samples comprising 20 normal, 20 polyp, 20 adenoma samples, 23 early stage cancer samples, and 14 advanced cancer samples. In addition, we confirmed the altered expression of two of the miRNAs in an independent dataset of 40 CRC samples and their paired normal tissues. We found circulating levels of miR-34a and miR-150 to be capable of distinguishing patients groups with benign and malignant diseases of the colon from each other, and sets of miRNAs that distinguish patients with advanced cancer from benign disease groups. Specifically, we found high levels of circulating miR-34a and low miR-150 levels to distinguish patients with polyps from those with advanced cancer, and low circulating miR-150 levels to separate patients with adenomas from those with advanced cancer.

Methods

Patients selection and sample collection

Cases with positive colonoscopy results for malignancy, confirmed by histology as colon or rectal carcinomas, were recruited between December 2007 and December 2010 at the Department of Surgery, Adelaide and Meath Hospital and at the Thomayer Hospital in Prague, Czech Republic. Control subjects or subjects diagnosed with polyps or adenomatous polyps were selected during the same period from individuals undergoing colonoscopy for various gastrointestinal complaints (macroscopic bleeding, positive faecal occult blood test or abdominal pain of unknown origin). The participating subjects gave written informed consent in accordance with the Declaration of Helsinki at the precipitating site that was approved by Tallaght Hospital/St. James’s Hospital Joint Research Ethics Committee, The Adelaide and Meath Hospital, Dublin, Incorporating The National Children’s Hospital, Tallaght, Dublin 24, Ireland and the Ethical Committee of the Institute of Experimental Medicine, Prague, Czech Republic. See Table 1 for clinical information on samples used.

Table 1 Clinical information on the discovery and validation plasma sample cohorts

Two separate patient cohorts were identified, a discovery set (n = 48) comprising 8 normal, 8 polyp, 16 adenoma samples, 8 early stage cancer samples (stage I/II), and 8 advanced cancer samples (stage III/IV), and a validation set (n = 97) comprising 20 normal, 20 polyp, 20 adenoma samples, 23 early stage cancer samples, and 14 advanced cancer samples. In addition, an independent public dataset [21] of quantitative real-time PCR (qRT-PCR) raw data was downloaded from the NCBI GEO archive (accession no: GSE28364) which contains information on 40 CRC samples and their paired normal tissues.

Plasma samples were collected according to standard phlebotomy procedures. 10 ml of blood sample was collected into EDTA plasma tubes and immediately placed in ice. The tubes were centrifuged at 1000 x g for 10 minutes at 4°C. Plasma was denuded by pipette from the cellular material, aliquoted into cryovial tubes, labelled and stored at -80°C until the time of analysis. The time from sample procurement to storage at -80°C was less than 3 hours. Each plasma sample underwent no more than 3 freeze/thaw cycles prior to analysis.

RNA extraction

Total RNA was isolated from 60 μl of each plasma sample using the miRNeasy mini kit (Cat no 217004, Qiagen). The Qiagen supplementary protocol (Purification of total RNA, including small RNAs, from serum or plasma) was utilised with the following modifications: thawed plasma samples were centrifuged at 1000 x g for 5 minutes at 4°C to remove excess debris from samples, RNA was extracted from the upper 50 μl of each sample. To elute the RNA, 50 μl of nuclease-free water was added to each spin column and incubated for 1 minute at room temperature before centrifuging into non-stick RNase-free microfuge tubes (Cat no AM12350, Ambion) to elute the RNA.

MiRNA profiling of plasma with TaqMan® low-density arrays

TaqMan® Array Human MicroRNA A and B Cards v2.0 (Cat no 4400238, Applied Biosystems) were employed to examine the expression of 667 miRNAs in 48 plasma samples in the discovery cohort. Reverse transcription and quantitative PCR (qPCR) were performed on equal volumes of RNA from each sample according to the manufacturer’s instructions using TaqMan® MicroRNA Reverse Transcription Kit (Cat no 4366596, Applied Biosystems) and Megaplex RT Primers to convert the miRNAs to cDNA, TaqMan® PreAmp Master Mix (Cat no 4391128, Applied Biosystems) and Megaplex PreAmp Primers for a preamplification step before real-time analysis. qPCR was performed using TaqMan® Universal Master Mix II, no UNG (Cat no 4440048, Applied Biosystems) on the 7900HT Fast Real-Time PCR system (Applied Biosystems). The Sequence Detector System software version 2.2.2 was utilised to generate study files using a fixed threshold value of 0.1 for statistical analysis (accession no: GSE67075).

Validation of miRNA expression using qRT-PCR

Individual TaqMan® miRNA assays were used for miRNA quantification in the 97 plasma samples in the validation cohort. To improve reverse transcription efficiency a miRNA multiplex RT primer pool was made from the singleplex RT primers of the four miRNAs to be analysed; miR-34a, miR-150, miR-923, and miR-let7e (this miRNA was used as the endogenous control as it showed very little variation in the discovery cohort, ΔCt SD = 0.865). 100 μl of each 20X RT primer were added to an RNase-free microfuge tube. The tube was dried in a speed vacuum (MAXI dry plus, Medical Supply Company, Ireland) at 50°C for 1 hour. The primers were re-suspended in 100 μl of nuclease-free water and 300 μl of 0.1X TE buffer was added to yield a 5X multiplex RT primer pool. The TaqMan® MicroRNA Reverse Transcription Kit (Cat no 4366596, Applied Biosystems) was used to perform reverse transcription reactions. Each reaction contained 1.8 μl of RT buffer (10X), 0.18 μl of dNTPs (25 mM), 3.6 μl of miRNA multiplex RT primer pool (5X), 1.2 μl of Multiscribe RT enzyme (50 U/μl), 5.22 μl of nuclease-free water and 6 μl of extracted total RNA. The reactions were incubated at 16°C for 30 minutes, 42°C for 30 minutes and 85°C for 5 minutes (G-STORM, GS1, Somerton Biotechnology Centre, UK).

Real-time PCR analysis was performed on 96 well plates (Cat no 4346906, Applied Biosystems). Technical triplicate PCRs were performed for each sample, and no template controls and a pooled sample containing cDNA from all 97 samples were included on each plate to ensure inter-plate reproducibility. Each reaction contained 1 μl of TaqMan miRNA assay (20X), 10 μl of TaqMan® Universal Master Mix II, no UNG (Cat no 4440048, Applied Biosystems), 7.67 μl of nuclease-free water, and 1.33 μl of cDNA. The reactions were incubated at 95°C for 10 minutes, and 40 cycles of 95°C for 15 seconds and 60°C for 15 seconds on the 7900HT Fast Real-Time PCR system (Applied Biosystems). The Sequence Detector System software version 2.2.2 was utilised to generate study files using a fixed threshold value of 0.1 for statistical analysis.

Statistical analysis

In the discovery cohort (n = 48), each miRNA was normalised by the ΔΔCt method using the average within sample Ct value [22]. This technique involves the use of the mean expression value of all expressed microRNAs in a given sample as a normalisation factor for microRNA real-time quantitative PCR data. Thus the average within sample Ct value for each card is calculated by averaging all miRNA Ct values for each individual sample. This was performed using the Bioconductor package HTqPCR (www.bioconductor.org). The non-parametric Kruskal-Wallis test was used to determine between group variations by rank as the data was not normally distributed. A Wilcoxon rank sum test was subsequently used to perform pair-wise comparisons between the 5 groups for the significant miRNAs identified by the Kruskal-Wallis test.

As an alternative to spiking un-related miRNA constructs into our samples we utilised the miRNA profiling data of the discovery cohort of samples to choose an appropriate endogenous control for use in the validation cohort. This involved analysing the expression of all 667 miRNAs across all 48 samples in the discovery cohort allowing us to choose one of the least variant miRNAs. As MammU6 showed highly variant expression in the discovery cohort, miR-let7e was chosen for use as an endogenous control for the validation set as it was one of the least variant miRNAs in the discovery phase experiment (ΔCt standard deviation of 0.86). When the let7e Cts were examined across all samples in the validation cohort this miRNA proved an appropriate endogenous control with a Ct standard deviation of 1.64. Statistically significant differences were determined using the non-parametric Wilcoxon rank sum test. The p-values for the validation set were adjusted using the Benjamini and Hochberg method [23] to account for multiple testing.

For consistency, the independent public dataset from Reid et al. [21] (accession no: GSE28364) was normalised using the same approach used to analyse the discovery cohort qRT-PCR data. This independent study used TaqMan® Array Human MicroRNA Cards v2.0 to analyse miRNA expression in 40 CRC tumour samples and their paired normal tissues. In order to mimic this structure in our validation plasma sample cohort, we grouped samples into ‘non-malignant’ and ‘malignant’ groups. As there were only two groups (normal versus cancer) in this analysis, the Wilcoxon rank sum test was used to determine significantly differentially regulated miRNAs. For this analysis of the validation cohort, miR-34a, miR-150 and miR-923 were first normalised against the endogenous control (miR-let7e) and the Wilcoxon rank sum test was used to determine significance between the groups.

Logistic regression (LR) and receiver operator characteristic (ROC) curve analysis were performed on miR-34a, miR-150 and miR-923 in the validation cohort. The markers were combined using LR and the ROC curves were used for interpretation of the models generated. The area under the curve (AUC) from the ROC curve for a given model was used to determine the probability of a correct prediction. The LR model for single miRNAs or combinations of miRNAs which gave the highest AUC was considered the most discriminating model and therefore the best marker at distinguishing between the groups of interest. All calculations were carried out in the R statistical environment (http://cran.r-project.org/) using the HTqPCR and stats packages.

Results

Differential expression of miRNAs in the discovery cohort

This study examined the expression of 667 miRNAs in plasma samples of a discovery cohort of 48 patients with benign and malignant disease of the colon compared to age and sex matched disease-free controls (Table 1). Statistical analysis revealed 73 miRNAs that have significantly different levels (p-value <0.05) in at least one of the disease groups (polyp n = 8, adenoma n = 16, early stage cancer (stage I/II) n = 8, and advanced cancer (stage III/IV) n = 8) compared to the healthy controls n = 8 (Table 2). Forty miRNAs were significantly different in advanced cancer, 22 in early stage cancer, 7 in adenoma, and 22 in the polyp group compared to normal controls. Ten miRNAs were significantly altered in both the early stage and advanced cancer groups compared to the normal controls; miR-923, miR-801, miR-144*, miR-135a*, miR-500, miR-497, miR-150, miR-30c and RNU48 showed lower levels while miR-532-3p was more abundant in the cancer groups compared to the controls. MiR-34a was significantly higher in early stage cancer compared to healthy controls.

Table 2 miRNAs with significantly different levels in disease groups compared to normal in the discovery cohort

miRNAs were prioritised for subsequent confirmation in the validation sample set if they showed consistently altered levels between the control group and each of the disease groups, and if all or most of these changes were deemed to be statistically significant. MiR-34a was chosen for validation as it was increased in the plasma of diseased patients compared to controls. MiR-150 and miR-923 were chosen for validation as their plasma levels progressively decreased as the sample/disease groups progress toward malignancy (Table 2).

Altered levels of miR-34a, miR-150, and miR-923 in the validation cohort

The three candidate miRNAs; miR-34a, miR-150, and miR-923 were analysed in a validation cohort of 97 independent plasma samples comprising 20 normal, 20 polyp, 20 adenoma samples, 23 early stage cancer samples, and 14 advanced cancer samples (see Table 3).

Table 3 Fold changes and associated p-values for miR-34a, miR-150, and miR-923 in the validation cohort

Among the three miRNAs analysed in the validation cohort, only miR-34a distinguished the normal and precancerous lesion groups from the disease samples (Figure 1A). MiR-34a expression was significantly increased in the adenoma (FC 2.09, p-value = 0.028) and early stage cancer (FC 2.84, p-value = 0.002) groups compared to healthy controls, and moderately in the advanced cancer group (FC 1.80, p-value = 0.081). The levels of this miRNA were also significantly higher in the adenoma (FC 2.71, p-value = 0.002), early stage cancer (FC 3.69, p-value = 3.90e-5) and advanced cancer (FC 2.34, p-value = 0.006) groups compared to the polyp samples. Alternatively, miR-150 distinguished the polyp and adenoma groups from the advanced cancer ((FC -2.31, p-value = 0.007, FC -2.55, p-value = 0.001 respectively) (Figure 1B). Following adjustment for multiple testing, miR-923 was not found to have significantly altered levels in the validation cohort (Figure 1C).

Figure 1
figure 1

MiR-34a, miR-150, and miR-923 in the validation cohort. Beeswarm plots of plasma levels of (A) miR-34, (B) miR-150, and (C) miR-923 in the five sample groups; normal (n = 20), polyp (n = 20), adenoma (n = 20), early stage cancer (stage I/II) (n = 23), and advanced cancer (stage III/IV) (n = 14). Expression levels of the miRNAs (log10 scale at y-axis) are normalised to let-7e. Statistically significant differences between groups were determined using Wilcoxon rank sum tests and are represented as follows; p-value <0.05 = *, p-value <0.01 = **, p-value < 0.001 = ***.

Validation of altered miR-34a, miR-150, and miR-923 expression in an independent dataset of matched colon tumour and normal tissues

In an effort to determine whether the levels of circulating cell free miRNAs that were observed in our study reflected the biology of the tumour, we investigated their expression in a publically available qPCR dataset of CRC tumours and matched normal tissues [21]. All miRNAs found to be significantly differentially expressed in the independent dataset are listed in Additional file 1. As this independent dataset only contains information on cancer and adjacent normal samples, the normal, polyp and adenoma plasma samples in our validation cohort were combined into a ‘non-malignant’ group and the early stage and advanced cancer plasma samples were combined into a ‘cancer’ group to facilitate comparison of the miRNA expression changes.

Figure 2 illustrates that the altered levels of all three circulating miRNAs observed in the plasma samples of the validation cohort mirrors the expression changes in the CRC tissues in the independent tissue sample set. Mir-34a was found to be significantly up-regulated in both the validation and the independent sample sets, 2.62 fold up-regulated (p-value = <0.001) and 1.71 fold up-regulated (p-value = 8.50e-6) respectively (Figure 2A). In further examining the expression of miR-34a in the combined plasma sample groups in the validation cohort, it was observed that it was also able to distinguish the non-malignant group (normal, polyp & adenoma groups) from the early stage cancer group (FC = 3.07, p-value = <0.001), and the normal and polyp groups from the benign adenoma group (FC = 1.90, p-value = 0.001). Additionally, miR-150 was significantly down-regulated in both the cancer plasma (FC = -1.29, p-value = 0.013) and cancer tissue (FC = -2.55, p-value = 4.104e-08) samples. Although not significantly, miR-923 was also down-regulated in both the cancer plasma (FC = -1.11, p-value = 0.779) and cancer tissue (FC = -2.26, p-value = 0.353) samples. See Table 3 for all FCs and associated p-values for miR-34a, miR-150, and miR-923 in the validation cohort.

Figure 2
figure 2

Expression levels of miR-34a, miR-150, and miR-923 in cancerous and non-malignant sample groups in the validation cohort and independent data set. (A) Cancer samples show significantly higher miR-34a expression (p-value < 0.001) compared to non-malignant samples in both the validation plasma samples and independent tumour tissue samples. (B) Cancer samples show significantly lower miR-150 expression compared to non-malignant samples in both the validation plasma samples (p-value =0.003) and independent tumour tissue samples (p-value < 0.001). (C) No significant change in miR-923 expression is observed in either the validation plasma samples (p-value = 0.56) or the independent tumour tissue samples (p-value = 0.35). The non-malignant group in the validation cohort represents the normal, polyp, and adenoma samples. Cancer groups in both studies represent cancers of all stages I-IV. Expression levels of the miRNAs (log10 scale at y-axis) are normalised to let-7e in the validation cohort and mean normalised in the independent dataset. Statistically significant differences between groups were determined using Wilcoxon rank sum tests in both the validation and independent cohorts.

Diagnostic potential of circulating miR-34a, miR-150, and miR-923 for the detection of disease of the colon

Once the altered expression of miR-34a, miR-150, and miR-923 was confirmed in our 97 validation plasma samples, and in the independent tissue samples, LR and ROC analyses were used to evaluate the potential of these miRNAs to distinguish between the disease and control blood plasma samples. First LR was used to identify the linear model with the best discriminatory power between sample groups, and the quality of this model was depicted by the area under the curve (AUC) of the ROC curve. For individual miRNAs, AUCs ranged from 0.488 to 0.875 (see Additional file 2). In an effort to identify the most powerful candidates as diagnostic markers of disease, we focused on comparing sample groups that showed significant differential expression of more than one of the three miRNAs. Thus we have identified sets of miRNAs that could distinguish advanced cancer from benign disease groups (Figure 3). Individually, plasma levels of miR-34a, miR-150, and miR-923 discriminate polyp samples from advanced cancer samples with AUC = 0.796 (CI:0.646-0.947), 0.825 (CI:0.681-0.969), and 0.746 (CI:0.412-0.817) respectively. The discriminatory power of the analysis was improved by combining markers miR-34a (p-value = 0.016, CI:-2.911- -0.314) and miR-150 (p-value = 0.031, CI:0.392-2.640) which increased the AUC to 0.904. The Akaike’s Information Criterion (AIC) was used as a measure of the quality of the model. The model with the lowest AIC was considered the best fit i.e. the combination of miRNAs that produced the lowest AIC. Where the addition of a miRNA did not lower the AIC it was excluded as the simplest model that best describes the data is preferable. MiR-923 expression was found not to improve the fit of the model and thus was not included in the combination (Figure 3A). In searching for an optimal model to distinguish adenoma samples from advanced cancer samples, it was found that miR-150 alone (rather than in combination with miR-923) was sufficient to with an AUC of 0.875 (CI:0.754-0.996) (Figure 3B).

Figure 3
figure 3

ROC curve analysis combining plasma levels of the three miRNAs to distinguish advanced cancer from benign disease groups. (A) Plasma levels of miR-34a, miR-150, and miR-923 discriminate polyp samples (n = 20) from advanced cancer (stage III/IV) samples (n = 14) individually; AUC = 0.796, 0.825, and 0.746 respectively. Combining markers miR-34a and miR-150 increases the AUC to 0.904. miR-923 expression does not improve the fit of the model and thus is not included in the combination. (B) Plasma levels of miR-150 and miR-923 discriminate adenoma samples (n = 20) from advanced cancer (stage III/IV) samples (n = 14) individually; AUC = 0.875 and 0.732 respectively. Combining these markers does not improve the fit of the model, thus miR-150 alone distinguishes adenoma samples from advanced cancer samples with an AUC of 0.875.

Discussion

Screening for the early detection of CRC is important to improve patient survival and facilitate cancer prevention through the detection and removal of polyps. The aim of this study therefore was to investigate the potential of circulating cell-free miRNAs not only as biomarkers of CRC, but also their efficiency at delineating patients presenting with precancerous lesions, i.e. polyps and benign adenomas from normal and cancer patient groups. MiRNA profiling was performed in the discovery sample cohort consisting of five groups; normal, polyps, benign adenomas, early stage cancer, and advanced cancer (Table 2) and identified three candidate miRNAs (miR-34a, miR-150, and miR-923) which were then further examined in a validation cohort of 97 samples divided into the same five groups as before (Figure 1 and Table 3). In addition, we confirmed that the altered circulating levels of miR-34a and miR-150 mirror the expression changes evaluated in the tumours of an independent dataset of 40 CRC samples and their paired normal tissues (Figure 2).

miR-34a is a p53-regulated miRNA that has been shown to influence both cellular senescence and apoptosis [24,25]. Different studies have demonstrated its up or down regulation in CRC compared to normal tissue (as reviewed in [26]). Wu et al. [27] demonstrated the involvement of this miRNA in CRC invasion and metastasis through targeting FRA1, a FOS transcription factor that is capable of forming activator protein-1 (AP-1) heterodimers. Increased levels of this circulating miRNA have been detected in patients with chronic hepatitis C and non-alcoholic fatty liver disease [28] and levels have been found to be decreased in whole blood samples of patients with CRC compared to healthy controls [29]. Brunet and colleagues studied miRNA expression in stage III CRC tissue samples compared to normal controls and found miR-34a to be significantly up-regulated [30].

In one of the first studies to investigate the altered expression of miRNAs in cancer [31] miR-150 was shown to be up-regulated in colorectal tissue compared to normal tissue. However several subsequent studies have shown this miRNA to be down-regulated in CRC tissue compared to normal tissue [21,32]. Indeed a recent study on 239 samples from Ma et al. found that miR-150 was down-regulated in adenoma and CRC tissues compared to normal tissue, and this down-regulation was associated with decreased overall survival and a worse response to adjuvant chemotherapy [33]. Decreased circulating levels of miR-150 have been identified in patients with acute myeloid leukaemia [34], and have been associated with poor prognosis for critically ill patients [35]. Furthermore, Wang and colleagues found miR-150 expression to be down-regulated in their 10 pooled CRC plasma samples compared to 10 pooled control samples, although its altered expression was not validated in their additional individual samples [15]. To our knowledge there are no studies outlining miR-923 expression in CRC or detecting circulating levels of this miRNA, however it has been shown to be down-regulated in chronic lymphocytic leukaemia patients [36] and up-regulated in taxol resistant breast cancer cells [37].

Of the three miRNAs analysed in our validation cohort, only miR-34a distinguished the normal group from the disease groups (Figure 1A). A large amount of inter-individual variability was noted in the normal samples assessed for miR-150 and miR-923 expression which may account for the fact that they do not significantly separate the normal group from the disease groups (Figure 1B&C). The reason for this variability may lie in the fact that these subjects had sufficient medical complaints to present themselves for colonoscopy, but although they do not present polyps, adenomas or cancer of the colon they may have had other conditions such as irritable bowel disease which may influence the results. In addition, the high number of adenoma samples (n = 16) compared to the other sample group numbers (n = 8) in the discovery cohort may explain why we observed statistically significant alterations in miR-150 and miR-923 in the initial analysis (Table 2) but not in the validation cohort (Table 3). Despite the variability within sample groups, we found circulating levels of miR-34a and miR-150 to be capable of distinguishing cancer patients from the non-malignant group of patients (Figure 2), in addition they were also capable of delineating patient groups with different diseases of the colon from each other (Table 3). Moreover, the discovery miRNA profiling results (Table 2) provide additional miRNA candidates, for instance miR-144-5p that may have potential as circulating miRNA biomarkers of CRC which can be exploited and independently validated by other research groups. In our opinion this is an important step towards the identification of specific biomarkers for early stages of disease.

There have been several publications examining the potential for miRNAs to act as circulating biomarkers for the detection of CRC. Recently, Faltejskova and colleagues attempted to validate the serum levels of four miRNAs (miR-17-3p, miR-29a, miR-92a and miR-135b) that had been proposed by other groups as potential circulating biomarkers of CRC. They used qPCR to assess the miRNA expression levels in 100 CRC patients and 30 healthy controls, and did not detect any significant changes in the expression of any of the miRNAs evaluated [17]. We examined the lists of significantly differentially expressed miRNAs in our discovery cohort to determine whether we also identified biomarkers found by other groups. We did not find miR-21 [19], miR-141 [38], miR-29a, miR-17-3p or miR-92 [14,16], miR-601 or miR-760 [15] to be significantly differentially expressed in any of the disease groups compared to our healthy controls. We did, however find miR-19a, miR-19b, and miR-15b significantly altered in some of our comparisons (see Table 2). These miRNAs were among those found by GirĂ¡ldez et al. to be significantly up-regulated in plasma samples of CRC patients in their study in 2012 [18].

Although there is some concordance among the results of different groups in the search for biomarker miRNAs, uncertainty remains as to which miRNAs are the most appropriate markers of disease. Disparities in patient age and time of sample collection (i.e. before or after surgery/treatment) in different studies may impact on the reproducibility of results. In addition to these variables, it has been noted previously that miRNA profiles vary between different ethnic groups [39], male and female patients [40], and that blood cell contaminants can contribute to circulating miRNA profiles [40,41]. Blood cell contaminants of plasma and serum samples may be of particular importance in evaluating the potential of circulating miRNAs as biomarkers of disease. In fact, Pritchard and colleagues suggest that the elevated miR-92a levels detected in the plasma of patients with colon cancer are due the higher levels of red blood cell haemolysis in patients with this disease [41]. This poses the question as to whether all of these putative biomarkers should be discarded due to their expression in blood cells, or whether extensive validation and perhaps additional profiling on larger more diverse patient cohorts will confirm the most reliable biomarker candidates. We would argue against discarding biomarkers because of their detection in hematopoietic cells, particularly if, as we have shown, their expression reflects that found within the tumour (Figure 2). If we were to discard markers for their presence in blood cells, we would also have to discount miR-21 as a valid marker of disease as it was also detected by Duttagupta and colleagues [40]. This miRNA is commonly up-regulated in cancer, has been identified in the serum and stool samples of cancer patients, and multiple studies have linked its expression to advanced disease and worse outcome for patients [42]. In an effort to control confounding factors in this study, all samples were age and sex matched, blood was taken at the time of colonoscopy before treatment commenced, and an additional centrifugation step to remove cellular debris prior to RNA extraction recommended by Duttagupta and colleagues [40] was included in the sample processing.

In order to examine the diagnostic potential of our three candidate miRNAs in detecting different diseases of the colon we employed ROC curve analysis. To identify the most powerful candidate combinations we focused our analysis on comparing sample groups that showed significant differential expression of more than one of the three miRNAs. This approach allowed us to identify marker combinations which distinguish patient groups with benign disease of the colon from those with advanced stage cancer (Figure 3). Specifically, miR-34a and miR-150 abundance were capable of differentiating patients with polyps from those with advanced cancer, AUC = 0.904, and miR-150 abundance separates patients with adenomas from those with advanced cancer, AUC = 0.875. To further confirm the true diagnostic potential of these circulating cell-free miRNAs it is now important for these results to be independently replicated in additional samples by another group. If this independent validation were successful, a prospective validation of the miRNA candidate biomarkers would be warranted.

Conclusions

To our knowledge this is the first study to examine circulating miRNA levels in samples from patients with polyps in addition to normal, adenoma, and early stage and advanced cancer samples. We identified two circulating miRNAs capable of distinguishing patient groups with different diseases of the colon from each other, and sets of miRNAs that distinguish patients with advanced cancer from benign disease groups. We also found miR-34a expression to be significantly increased in early stage cancer samples compared to the non-malignant patient samples, and in adenoma samples compared to normal and polyp samples.

Abbreviations

CRC:

Colorectal cancer

LR:

Logistic regression

ROC:

Receiver operator curve

AUC:

Area under the curve

FC:

Fold change

References

  1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA: Cancer J Clin. 2011;61(2):69–90.

    Google Scholar 

  2. Arends JW. Molecular interactions in the Vogelstein model of colorectal carcinoma. J Pathol. 2000;190(4):412–6.

    Article  CAS  PubMed  Google Scholar 

  3. Manne U, Shanmugam C, Katkoori VR, Bumpers HL, Grizzle WE. Development and progression of colorectal neoplasia. Cancer Biomark. 2010;9(1-6):235–65.

    PubMed  PubMed Central  Google Scholar 

  4. Levin B, Lieberman DA, McFarland B, Smith RA, Brooks D, Andrews KS, et al. Screening and surveillance for the early detection of colorectal cancer and adenomatous polyps, 2008: a joint guideline from the American Cancer Society, the US Multi-Society Task Force on Colorectal Cancer, and the American College of Radiology. CA: Cancer J Clin. 2008;58(3):130–60.

    Google Scholar 

  5. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Esquela-Kerscher A, Slack FJ. Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006;6(4):259–69.

    Article  CAS  PubMed  Google Scholar 

  7. Chen X, Ba Y, Ma L, Cai X, Yin Y, Wang K, et al. Characterization of microRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008;18(10):997–1006.

    Article  CAS  PubMed  Google Scholar 

  8. Gilad S, Meiri E, Yogev Y, Benjamin S, Lebanony D, Yerushalmi N, et al. Serum microRNAs are promising novel biomarkers. PLoS One. 2008;3(9), e3148.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Brase JC, Johannes M, Schlomm T, Falth M, Haese A, Steuber T, et al. Circulating miRNAs are correlated with tumor progression in prostate cancer. Int J Cancer. 2011;128(3):608–16.

    Article  CAS  PubMed  Google Scholar 

  10. Kan CW, Hahn MA, Gard GB, Maidens J, Huh JY, Marsh DJ, et al. Elevated levels of circulating microRNA-200 family members correlate with serous epithelial ovarian cancer. BMC Cancer. 2012;12:627.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Heegaard NH, Schetter AJ, Welsh JA, Yoneda M, Bowman ED, Harris CC. Circulating micro-RNA expression profiles in early stage nonsmall cell lung cancer. Int J Cancer. 2012;130(6):1378–86.

    Article  CAS  PubMed  Google Scholar 

  12. Bianchi F, Nicassio F, Marzi M, Belloni E, Dall’olio V, Bernard L, et al. A serum circulating miRNA diagnostic test to identify asymptomatic high-risk individuals with early stage lung cancer. EMBO Mol Med. 2011;3(8):495–503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. 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(1), e53141.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Huang Z, Huang D, Ni S, Peng Z, Sheng W, Du X. Plasma microRNAs are promising novel biomarkers for early detection of colorectal cancer. Int J Cancer. 2010;127(1):118–26.

    Article  CAS  PubMed  Google Scholar 

  15. Wang Q, Huang Z, Ni S, Xiao X, Xu Q, Wang L, et al. Plasma miR-601 and miR-760 are novel biomarkers for the early detection of colorectal cancer. PLoS One. 2012;7(9), e44398.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Ng EK, Chong WW, Jin H, Lam EK, Shin VY, Yu J, et al. Differential expression of microRNAs in plasma of patients with colorectal cancer: a potential marker for colorectal cancer screening. Gut. 2009;58(10):1375–81.

    Article  CAS  PubMed  Google Scholar 

  17. Faltejskova P, Bocanek O, Sachlova M, Svoboda M, Kiss I, Vyzula R, et al. Circulating miR-17-3p, miR-29a, miR-92a and miR-135b in serum: Evidence against their usage as biomarkers in colorectal cancer. Cancer Biomark. 2012;12(4):199–204.

    PubMed  Google Scholar 

  18. Giraldez MD, Lozano JJ, Ramirez G, Hijona E, Bujanda L, Castells A, et al. Circulating microRNAs as biomarkers of colorectal cancer: results from a genome-wide profiling and validation study. Clin Gastroenterol Hepatol. 2013;11(6):681–8. e683.

    Article  CAS  PubMed  Google Scholar 

  19. Kanaan Z, Rai SN, Eichenberger MR, Roberts H, Keskey B, Pan J, et al. Plasma miR-21: a potential diagnostic marker of colorectal cancer. Ann Surg. 2012;256(3):544–51.

    Article  PubMed  Google Scholar 

  20. Pu XX, Huang GL, Guo HQ, Guo CC, Li H, Ye S, et al. Circulating miR-221 directly amplified from plasma is a potential diagnostic and prognostic marker of colorectal cancer and is correlated with p53 expression. J Gastroenterol Hepatol. 2010;25(10):1674–80.

    Article  CAS  PubMed  Google Scholar 

  21. Reid JF, Sokolova V, Zoni E, Lampis A, Pizzamiglio S, Bertan C, et al. miRNA profiling in colorectal cancer highlights miR-1 involvement in MET-dependent proliferation. Mol Cancer Res. 2012;10(4):504–15.

    Article  CAS  PubMed  Google Scholar 

  22. Mestdagh P, Van Vlierberghe P, De Weer A, Muth D, Westermann F, Speleman F, et al. A novel and universal method for microRNA RT-qPCR data normalization. Genome Biol. 2009;10(6):R64.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.

    Google Scholar 

  24. Chang TC, Wentzel EA, Kent OA, Ramachandran K, Mullendore M, Lee KH, et al. Transactivation of miR-34a by p53 broadly influences gene expression and promotes apoptosis. Mol Cell. 2007;26(5):745–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Fujita K, Mondal AM, Horikawa I, Nguyen GH, Kumamoto K, Sohn JJ, et al. p53 isoforms Delta133p53 and p53beta are endogenous regulators of replicative cellular senescence. Nat Cell Biol. 2009;11(9):1135–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Schetter AJ, Okayama H, Harris CC. The role of microRNAs in colorectal cancer. Cancer J. 2012;18(3):244–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Wu J, Wu G, Lv L, Ren YF, Zhang XJ, Xue YF, et al. MicroRNA-34a inhibits migration and invasion of colon cancer cells via targeting to Fra-1. Carcinogenesis. 2012;33(3):519–28.

    Article  CAS  PubMed  Google Scholar 

  28. Cermelli S, Ruggieri A, Marrero JA, Ioannou GN, Beretta L. Circulating microRNAs in patients with chronic hepatitis C and non-alcoholic fatty liver disease. PLoS One. 2011;6(8), e23937.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Nugent M, Miller N, Kerin MJ. Circulating miR-34a levels are reduced in colorectal cancer. J Surg Oncol. 2012;106(8):947–52.

    Article  CAS  PubMed  Google Scholar 

  30. Brunet Vega A, Pericay C, Moya I, Ferrer A, Dotor E, Pisa A, et al. microRNA expression profile in stage III colorectal cancer: circulating miR-18a and miR-29a as promising biomarkers. Oncol Rep. 2013;30(1):320–6.

    PubMed  Google Scholar 

  31. Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, et al. A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A. 2006;103(7):2257–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Slattery ML, Wolff E, Hoffman MD, Pellatt DF, Milash B, Wolff RK. MicroRNAs and colon and rectal cancer: differential expression by tumor location and subtype. Genes Chromosomes Cancer. 2011;50(3):196–206.

    Article  CAS  PubMed  Google Scholar 

  33. Ma Y, Zhang P, Wang F, Zhang H, Yang J, Peng J, et al. miR-150 as a potential biomarker associated with prognosis and therapeutic outcome in colorectal cancer. Gut. 2012;61(10):1447–53.

    Article  PubMed  Google Scholar 

  34. Fayyad-Kazan H, Bitar N, Najar M, Lewalle P, Fayyad-Kazan M, Badran R, et al. Circulating miR-150 and miR-342 in plasma are novel potential biomarkers for acute myeloid leukemia. J Transl Med. 2013;11:31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Roderburg C, Luedde M, Vargas Cardenas D, Vucur M, Scholten D, Frey N, et al. Circulating microRNA-150 serum levels predict survival in patients with critical illness and sepsis. PLoS One. 2013;8(1), e54612.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Zhu DX, Zhu W, Fang C, Fan L, Zou ZJ, Wang YH, et al. miR-181a/b significantly enhances drug sensitivity in chronic lymphocytic leukemia cells via targeting multiple anti-apoptosis genes. Carcinogenesis. 2012;33(7):1294–301.

    Article  CAS  PubMed  Google Scholar 

  37. Zhou M, Liu Z, Zhao Y, Ding Y, Liu H, Xi Y, et al. 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(28):21496–507.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Cheng H, Zhang L, Cogdell DE, Zheng H, Schetter AJ, Nykter M, et al. Circulating plasma MiR-141 is a novel biomarker for metastatic colon cancer and predicts poor prognosis. PLoS One. 2011;6(3), e17745.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. 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(10), e13735.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Duttagupta R, Jiang R, Gollub J, Getts RC, Jones KW. Impact of cellular miRNAs on circulating miRNA biomarker signatures. PLoS One. 2011;6(6), e20769.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Pritchard CC, Kroh E, Wood B, Arroyo JD, Dougherty KJ, Miyaji MM, et al. Blood cell origin of circulating microRNAs: a cautionary note for cancer biomarker studies. Cancer Prev Res. 2012;5(3):492–7.

    Article  CAS  Google Scholar 

  42. Mazeh H, Mizrahi I, Ilyayev N, Halle D, Brucher B, Bilchik A, et al. The diagnostic and prognostic role of microRNA in colorectal cancer - a comprehensive review. J Cancer. 2013;4(3):281–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Funding is acknowledged from the Science Foundation Ireland Strategic Research Cluster, Molecular Therapeutics for Cancer Ireland (www.mtci.ie) and from the Grant Agency of the Czech Republic (grants CZ GA CR: GAP304/10/1286 and 1585).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sinéad T Aherne.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

STA performed the RNA extraction, miRNA profiling of plasma with TaqMan Low-Density Arrays, validation of miRNA expression using qRT-PCR, participated in analysis and interpretation of the results and prepared the manuscript. SFM performed all of the bioinformatic/statistical analysis. STA, SFM, DJH, BP, AC, ML, PD, PN, PD and MC contributed to the result interpretation and manuscript preparation. DJH, BP, AC, ML, PD, PN assisted in the collection and provision of clinical samples. PD and MC conceived the study, participated in its design, coordination and interpretation of the results and finalized the manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1:

Significantly differentially expressed miRNAs in the independent dataset.

Additional file 2:

AUC values for miR-34a, miR-150, and miR-923 in the validation cohort.

Rights and permissions

This article is published under license to 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Aherne, S.T., Madden, S.F., Hughes, D.J. et al. Circulating miRNAs miR-34a and miR-150 associated with colorectal cancer progression. BMC Cancer 15, 329 (2015). https://doi.org/10.1186/s12885-015-1327-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-015-1327-5

Keywords