Abstract
Colorectal cancer (CRC) is the third most common malignancy and the second most common cause of cancer death worldwide. Early detection of CRC can improve patient survival rates; thus, the identification of noninvasive diagnostic markers is urgently needed. MicroRNAs (miRNAs) have extensive potential to diagnose several diseases, including cancer. In this study, we compared the expression pattern of miRNAs from plasma and stool samples of patients with early stages of CRC (I, II) with that of healthy subjects. We performed miRNA profiling using microarrays on plasma and stool samples of eight patients with CRC and four healthy subjects. Seven miRNAs were found to be underexpressed in both plasma and stool samples of patients with CRC versus healthy subjects. Then, we aimed to verify two out of these seven differentially expressed miRNAs (let-7a-5p and let-7f-5p) by quantitative reverse transcriptase polymerase chain reaction on a larger set of plasma and stool samples of 51 patients with CRC and 26 healthy subjects. We confirmed the results of microarray analysis since their expression was significantly lower in stool and plasma samples of patients with CRC. Moreover, receiver operating characteristic curve analysis demonstrated that fecal let-7f expression levels have significant sensitivity and specificity to distinguish between patients with CRC and healthy subjects. In conclusion, if the results are confirmed in larger series of patients, underexpressed let-7a-5p and let-7f-5p miRNAs in both plasma and stool samples of patients with CRC may serve potentially as noninvasive molecular biomarkers for the early detection of CRC.
Introduction
Despite the existence of screening, preventive, and therapeutic strategies, colorectal cancer (CRC) is still the third most common malignancy and the second most common cause of cancer death, with an estimated 1.2 million new cases and more than 600,000 deaths annually, worldwide.1–3 Initiation and progression of CRC is a slow process lasting several years, and it usually begins with cryptic lesions and gradually develops to adenoma and finally to malignant carcinoma, with the potential for invasion and metastasis. 4 Early detection of this neoplasm is critical because patient survival rate and response to treatment are strongly associated with tumor stage and disease progression at the time of diagnosis. It has been estimated that almost 95% of patients would be cured if diagnosis was made at an early or precursor lesion stage. 5
Therefore, new, accurate, and reliable diagnostic methods that improve early detection and monitoring of CRC are necessary.
MicroRNAs (miRNAs) are short (18–24 nucleotides), highly conserved noncoding endogenous RNAs, which regulate gene expression posttranscriptionally.6,7 Numerous studies have shown that the alteration of miRNAs expression is involved in the pathogenesis of several types of cancers, such as CRC, by regulating the expression of tumor suppressors and oncogenes.8,9 Differential expression of miRNAs has also been detected in serum, plasma,10–15 and stool of patients with CRC.16–22
In this study, we compared the expression pattern of miRNAs in paired plasma and stool samples of patients with early stages (I, II) of CRC and of neoplasm-free individuals as controls.
Materials and Methods
Patients
In this study, a total of 51 patients with CRC and 26 healthy controls were enrolled for the evaluation of miRNAs from plasma and stool samples at the Digestive Disease Research Institute, Shariati Hospital, Tehran, Iran, from August 2011 to October 2012. Patients with CRC and healthy controls had similar age (P-value = 0.776, independent samples t-test) and gender (P-value = 0.567, chi-square test) distribution. Clinical and demographical data of participants are presented in Table 1.
Characteristics of patients with CRC and healthy subjects.
TNM: tumor-node-metastasis staging system according to the Union for international Cancer Control (UICC) guidelines.
P-value for age was calculated by independent sample t-test; P-value for gender was calculated by chi-square test.
Plasma and stool samples were collected from 51 newly diagnosed, pathologically confirmed patients with early stages (I, II) of CRC before receiving any treatment and from 26 healthy subjects who were free of CRC and colorectal polyps as examined by colonoscopy. Tumor was staged according to the tumor node metastasis staging system of the Union for International Cancer Control guidelines. The patients with a history of concurrent malignant diseases were excluded from this study. Plasma and stool samples were stored at –-80°C until use. Informed consent was obtained from all participants for the use of their blood and stool samples in this study. This research complied with the principles of the Declaration of Helsinki and was approved by the Institutional Review Board of Shariati Hospital.
RNA isolation
Extraction of total RNA, including miRNAs, from plasma and stool samples was carried out using miRNAeasy Serum/Plasma Kit (Qiagen GmbH) and miRNAeasy Mini Kit (Qiagen GmbH), respectively. For stool samples, the protocol was slightly modified to increase the RNA yield. Briefly, ~100 mg of stool was homogenized with 1 mL of RNase-free water and 200 μL of this was lysed with QIAzol lysis reagent in a proportion of 1:6. RNA was precipitated with chloroform, and the aqueous phase was mixed with 1.5 volumes of 100% ethanol. 17 The concentration of extracted RNA was measured using the Nano-Drop-1000 Spectrophotometer (Thermo Fisher Scientific, Inc.). Agilent's Bioanalyzer (Agilent Technologies) was used for checking the quality of total RNA with the RNA 6000 chip and miRNA with the small RNA chip (Agilent Technologies).
MicroRNA microarray expression profiling and data analysis
Eight paired stool and plasma samples of patients with CRC and four paired stool and plasma samples of healthy individuals were randomly selected to identify the miRNA expression signature. MiRNA profiling was done with Agilent's miRNA microarray system (V19), containing 2,006 human miRNAs (Agilent Technologies). RNA labeling and hybridization were performed according to Agilent's protocol version 2.4. Briefly, 100 ng of the total RNA was dephosphory-lated with calf intestinal alkaline phosphatase by incubation at 37°C for 30 minutes. In the denaturation step, dimethyl sulfoxide (100%) was added to each sample, and after incubation at 100°C for eight minutes, the samples were immediately transferred to an ice–water bath. Then, samples were labeled with Cyanine3-pCp by incubation with T4 RNA ligase for two hours at 16°C. After labeling, samples were completely dried in a vacuum concentrator at 45°C-55°C and resuspended in nuclease-free water. The samples were hybridized to array active side in Agilent SureHyb Chambers at 55°C for 20 hours. The microarrays were then washed with washing buffers 1 and 2 and scanned on an Agilent's Microarray Scanner. The information of probe features was extracted from the microarray scan data using Agilent Feature Extraction Software. Finally, the miRNA expression analysis was performed using GeneSpring GX Analysis software (Agilent Technologies).
Statistical significance of differential expression between the CRC and control groups was determined by t-test. The false discovery rate method of Benjamini-Hochberg procedure was used to correct for multiple comparisons. miRNAs with a fold-change of >2, and P-values <0.05 were considered differentially expressed.
Validation of microarray results with quantitative reverse transcriptase polymerase chain reaction (qRT-PCR)
A set of 51 paired plasma and stool samples of CRC patients (stage I, n = 12; stage II, n = 39) and 26 paired plasma and stool samples of healthy controls (including the patients and healthy subjects used for miRNA microarray analysis) were used for validation of the miRNA microarray findings by SYBR-Green-based quantitative reverse transcriptase polymerase chain reaction (qRT-PCR). Two miRNAs, let-7a-5p and let-7f-5p, were selected for validation based on miRNA array results and also because of the known role of let-7 family of miRNAs in human cancer. Reverse transcription of the extracted RNA was performed by the miScript Reverse Transcription Kit (Qiagen GmbH), according to the manufacturer's recommendations. Quantification of miRNA was performed using SYBR Green miScript PCR system (Qiagen GmbH) on Light Cycler v.3.5 software (Roche Applied Science). Each reaction was performed in a 20 μL volume with 10 ng of template cDNA, and all reactions were run in duplicates. RNU6B miScript Primer Assay (Qiagen GmbH) was used as an endogenous control for normalization of the data. Relative Expression Software Tool (REST 2009) was used to evaluate statistically significant differences in miRNA expression between the two groups. 23 P-values <0.05 were considered statistically significant.
Results
Plasma and stool miRNA microarray profiling
RNA concentrations ranged from 20 to 51 ng/μL and from 62 to 1,430 ng/μL in plasma and stool samples, respectively. Quality of total RNA and miRNA was acceptable to perform subsequent miRNAs microarrays and the validation study. Analysis of our array results from comparisons of plasma and stool samples of eight patients with CRC versus four healthy subjects using GeneSpring GX software showed several significantly differentially expressed miRNAs. At the fold change cutoff of 2.0 and P < 0.05, q < 0.3 results of CRC plasma miRNA arrays in comparison with normal samples revealed 78 significantly downregulated and two significantly upreg-ulated miRNAs, and for stool samples, our results revealed 148 significantly downregulated and two significantly upreg-ulated miRNAs. All data of microarray results are listed in Supplementary Table 1. Interestingly, seven miRNAs were downregulated in both plasma and stool specimens of patients with CRC compared with those of healthy controls (P < 0.05, q < 0.3; Fig. 1).

Diagramatic representation of the miRNAs differentially expressed in plasma and stool of patients with CRC compared to healthy subjects. The venn diagram shows 150 differentially expressed miRNAs in stool and 80 in plasma compared to healthy individuals. The intersection depicts seven miRNAs underexpressed in both plasma and stool of patients with CRC, which are shown in detail in a table.
Validation of miRNA expression differences by qRT-PCR
Two miRNAs (let-7a-5p and let-7f-5p) out of the seven differentially expressed miRNAs in both plasma and stool of patients with CRC by microarray analysis were validated by qRT-PCR on plasma and stool samples of 51 patients with CRC and 26 healthy subjects. We observed significantly lower expression by Mann-Whitney test of let-7a-5p (P = 0.02) and let-7f-5p (P = 0.02) in plasma samples of patients with CRC and of let-7a-5p (P = 0.03) and let-7f-5p (P = 0.003) in stool samples of patients with CRC compared with that of the normal group samples, confirming the results obtained with microarray analysis (Figs. 2 and 3). There was no correlation between let-7a or let-7f expression and cancer stage or tumor site in patients with CRC.

Plasma let-7a-5p and let-7f-5p expression levels in 51 patients with CRC versus controls. Expression levels of the miRNAs were normalized to RNU6B. (A) Let-7a is significantly downregulated in patient samples compared to normal ones (P-value = 0.02); (B) Let-7f is significantly downregulated in plasma samples of patients compared to normal group (P-value = 0.02). Data were analyzed using Mann-Whitney test. Graphs represent the mean value and standard error of mean of each group.

Stool let-7a-5p and let-7f-5p expression level in 51 patients with CRC versus controls. Expression levels of the miRNAs were normalized to RNU6B. (A) Let-7a is significantly downregulated in stool samples of patients compared with the normal group (P-value = 0.03); (B) Let-7f is significantly downregulated in samples of patients in comparison with the normal group (P-value = 0.003). Data were analyzed using Mann-Whitney test. Graphs represent the mean value and standard error of mean of each group.
Evaluation of specificity and sensitivity of the let-7 miRNAs
Receiver operating characteristic (ROC) curve analyses were performed on plasma and fecal let-7a-5p and let-7f-5p separately. The ROC curve analysis demonstrated that let-7f-5p expression level in the stool samples has a significant specificity and sensitivity to discriminate patients with CRC from normal individuals. The area under curve was 70.9% with a standard deviation error of 0.06 (P-value = 0.003; Fig. 4).

Area under curve (AUC) of receiver operating characteristic (ROC) analysis for let-7f relative expression in fecal samples. ROC curve analysis revealed that fecal let-7f relative expression has significant specificity and sensitivity to distinguish between patients with CRC and normal individuals. The AUC value is 70.9% with the standard error of 0.06 (P-value = 0.003).
Discussion
To our knowledge, no study has been done on simultaneous evaluation of miRNA expression in both plasma and stool samples of patients with CRC to date. In this study, for the first time, we investigated simultaneously the plasma and stool miRNA expression in patients with early stages (I, II) of CRC and in neoplasm-free individuals as controls. We believed that aberrant expression of specific miRNAs on these types of samples might help to identify potential CRC biomarkers.
When expression profiles of miRNA from plasma and stool samples of the patients with CRC were compared with those of the healthy subjects, 80 and 150 differentially expressed miRNAs, respectively, were detected. From the differentially expressed miRNAs, only seven were seen in both plasma and stool samples. With regard to cancer, let-7a-5p and let-7f-5p were the most intriguing from those seven miRNAs, and we wanted to confirm their downregulation in a bigger cohort of patients by qRT-PCR. Our validation study agreed with the results from microarrays; expression levels of the let-7a-5p and let-7f-5p miRNAs were significantly decreased in both plasma and stool samples of patients with early stages of CRC compared with those of healthy subjects. Moreover, ROC curve analysis demonstrated that fecal let-7f expression levels have significant sensitivity and specificity to distinguish between patients with CRC and healthy subjects.
Reduced expression of let-7 family members has been reported in many tumors.24–28 For example, lower expression of let-7f has been reported as a marker for early stage ovarian cancer detection. 29 Moreover, there are numerous reports that indicate that this miRNA family is differentially expressed in CRC.22,30–33 However, there is very little information of downregulation of the let-7 miRNA family members in CRC plasma or stool samples. Wang et al performed miRNA profiling on two pooled plasma samples of 10 patients with CRC and 10 normal controls targeting 742 miRNAs by qRT-PCR. 13 Interstingly, let-7a miRNA was downregulated, among others, with a fold change >5.
It is well known that let-7f miRNAs are located at chromosome 9q22.3, forming a cluster along with let-7d miRNA. However, the mechanisms of action of these let-7 family members on mRNA are not fully understood, even though some of their targets are known. Experimentally verified let-7a target genes with strong evidence include some known cancer genes, such as BCL2, CASP8, CCND2, DICER1, EWSR1, FOXA1, HMGA1, HMGA2, HRAS, KRAS, MYC, MPL, NF2, NRAS, and PRDM1 (according to the databases miR-TarBase 34 and COSMIC, 35 accessed on September 10, 2015).
Experimentally verified let-7f-5p target cancer genes include CCND1, MYC, MPL, NKX2–1, and PRDM1 (according to the same databases). However, it is important to note that all the let-7 variants share an identical seed sequence that is critical for target recognition. 36 These target genes of let-7 miRNAs are involved in carcinogenesis regulating directly or indirectly cell cycle, apoptosis, and cell adhesion. For example, reduced expression of let-7 is related to cell proliferation and oncogenic transformation by increasing expression of HMGA2.37,38 Moreover, restoring let-7 levels in pancreatic cancer-derived cell lines strongly inhibits cell proliferation. 24 More specifically with regard to let-7a and let-7f, several studies have connected these miRNAs with oncogenesis and/or cell growth. Let-7f has been reported to modulate β-catenin/Wnt signaling that is related to growth and development of stem cells via its binding to Axin 2 and subsequent degradation of β-catenin in human mesenchymal stem cells. 39 It has been also reported that the miRNA let-7a can revert cell growth in lymphoma and in lung cancer.40,41 In addition, let-7a plays a role in STAT3 signaling at least in cervical carcinogenesis. 42 Interestingly, regarding CRC, previous results have shown that transfection of let-7a-1 precursor miRNA to human colon cancer cells resulted in lower levels of RAS and c-myc proteins and subsequent significant growth suppression. 33
It is important to note that hsa-miR-21–5p was downregu-lated in both plasma and stool of patients with CRC compared with that of healthy subjects after the microarray analysis. Surprisingly, this miRNA is upregulated in patients with CRC, even in their serum, according to previous studies. 43 However, it has been associated with metastasis and poor survival, and our patients with CRC were all in early and curable stage. It seems that miR-21 has significant potential as a prognostic marker, rather than as an early detection marker. 43
Our validation study results confirmed that the expression level of the let-7a-5p and let-7f-5p was significantly decreased in both plasma and stool samples of patients with stages I and II CRC compared with those of healthy subjects. Moreover, ROC curve analysis demonstrated that fecal let-7f-5p expression levels have significant sensitivity and specificity to distinguish between patients with CRC and healthy subjects. We suggest that the downregulation of fecal let-7f-5p in patients with CRC may be useful in early diagnosis of this cancer. Nevertheless, further studies on larger series of patients, especially with different stages of CRC, adenoma, and other precancer-ous lesions, are needed to clarify whether these miRNAs have the potential to be used for early detection of CRC.
Author Contributions
Conceived and designed the experiments: RG and RM. Performed the experiments: RG and NM. Analyzed the data: NN, AM, S. Khorrami, EA, and MP. Contributed reagents/materials/analysis tools: NM, VKS, and S. Knuutila. Prepared and wrote the paper: RG, NM, GA, and S. Knuutila. All the authors reviewed and approved the final manuscript.
Supplementary Material
Differentially expressed miRNAs in patients with CRC versus normal subjects.
