Human colon cancer profiles show differential microRNA expression depending on mismatch repair status and are characteristic of undifferentiated proliferative states
- Aaron L Sarver1Email author,
- Amy J French2,
- Pedro M Borralho3,
- Venugopal Thayanithy4,
- Ann L Oberg5,
- Kevin AT Silverstein1,
- Bruce W Morlan5,
- Shaun M Riska5,
- Lisa A Boardman6,
- Julie M Cunningham2,
- Subbaya Subramanian4,
- Liang Wang2,
- Tom C Smyrk2,
- Cecilia MP Rodrigues3,
- Stephen N Thibodeau2 and
- Clifford J Steer7, 8
© Sarver et al; licensee BioMed Central Ltd. 2009
Received: 25 August 2009
Accepted: 18 November 2009
Published: 18 November 2009
Colon cancer arises from the accumulation of multiple genetic and epigenetic alterations to normal colonic tissue. microRNAs (miRNAs) are small, non-coding regulatory RNAs that post-transcriptionally regulate gene expression. Differential miRNA expression in cancer versus normal tissue is a common event and may be pivotal for tumor onset and progression.
To identify miRNAs that are differentially expressed in tumors and tumor subtypes, we carried out highly sensitive expression profiling of 735 miRNAs on samples obtained from a statistically powerful set of tumors (n = 80) and normal colon tissue (n = 28) and validated a subset of this data by qRT-PCR.
Tumor specimens showed highly significant and large fold change differential expression of the levels of 39 miRNAs including miR-135b, miR-96, miR-182, miR-183, miR-1, and miR-133a, relative to normal colon tissue. Significant differences were also seen in 6 miRNAs including miR-31 and miR-592, in the direct comparison of tumors that were deficient or proficient for mismatch repair. Examination of the genomic regions containing differentially expressed miRNAs revealed that they were also differentially methylated in colon cancer at a far greater rate than would be expected by chance. A network of interactions between these miRNAs and genes associated with colon cancer provided evidence for the role of these miRNAs as oncogenes by attenuation of tumor suppressor genes.
Colon tumors show differential expression of miRNAs depending on mismatch repair status. miRNA expression in colon tumors has an epigenetic component and altered expression that may reflect a reversion to regulatory programs characteristic of undifferentiated proliferative developmental states.
Colon cancer (CC) is the second most frequent cause of cancer-related death in the United States and in Europe [1, 2]. The development of CC is considered a stepwise process with the accumulation of different genetic and epigenetic alterations. Most tumors (~85%) are generated by chromosomal instability (CIN) and associated with high frequency aneuploidy and allelic imbalance. The remaining 15% have defective DNA mismatch repair (dMMR), which is frequently measured by either the presence of microsatellite instability (MSI) or by testing for loss of the protein products for genes involved in DNA mismatch repair, most commonly MLHI, MSH2, MSH6, and PMS2 . Sporadic CCs with dMMR have distinctive clinical and pathological features that include proximal colon predominance, poor differentiation and/or mucinous histology, intra- and peritumoral lymphocytic infiltration, diploid DNA content  and generally have a better prognosis . The presence of MSI-H tumor phenotype and loss of protein expression of MLHl and MSH2 is highly concordant .
Recent progress has been made in CC screening  and treatment protocols . Currently, the most important prognostic factor for patients is pathologic tumor staging based on the tumor-node-metastasis (TNM) system . However, several pathologic and clinical features have been associated with increased risk of tumor recurrence in resectable CC . Between 20 and 25% of stage II CC patients develop recurrence and die from the disease . Therefore, it is imperative to identify and develop accurate, reliable and sensitive biomarkers.
Genome-wide approaches have reshaped the landscape of cancer research. Meta-analysis of multi-study data has allowed the identification of overlapping sets of differentially expressed genes that may have biomarker potential . Emphasis on genome-wide gene expression analyses has been driven by the general view that alterations in protein-coding oncogenes or tumor-suppressor genes underlie tumorigenesis. However, the discovery of a growing class of small non-coding RNAs, termed miRNAs, has opened a new field of cancer research and revealed the complexity of cancer biology. miRNAs regulate gene expression post-transcriptionally by translational attenuation/repression or cleavage of target mRNAs , thereby adding a new dimension to the regulation of gene expression . Further, aberrant expression of miRNAs has been associated with a growing list of cancers .
The potential use of miRNAs in diagnosis and prognosis has been demonstrated for several forms of cancer [14, 15]. miRNA expression profiles may be better-suited targets for the discovery of novel cancer biomarkers compared to gene expression profiles. This is supported by a report demonstrating the ability of miRNA profiles to correctly classify human cancers of unknown primary origin as well as poorly differentiated tumors . A growing number of studies have addressed miRNA expression in CC [17–20]. However, comparison across studies is limited by differences in profiling platform, quantity of miRNA obtainable, methodology, in some cases sample number and a paucity of clinicopathologic data. Consequently, translation of results into clinically useful and widely applicable biomarkers is hampered. Importantly, a potentially strong contributor to the variability of data among different studies relates to the tumor resection procedures. The inadvertent collection of surrounding residual non-tumor tissue may skew experimental results, diluting quantitative estimates of particular miRNA species based on the extent and type of tissue present in the sample. Most studies do not specifically address this potential problem.
In this study, global miRNA expression was evaluated in 80 colon tumors and 28 normal mucosa samples using the BeadArray™ platform (Illumina, Inc.) to evaluate global miRNA expression of 735 miRNA targets . Our results demonstrate that, in a sufficiently statistically powered number of tumors, a larger set of miRNAs than has previously been reported is differentially expressed between normal colon and tumor tissue. Additionally, specific miRNAs are significantly differentially expressed between tumors of deficient and proficient mismatch repair status (dMMR and pMMR, respectively).
Colon Cancer Samples
CC biospecimens were obtained from Mayo Clinic Rochester patients with colorectal neoplasia beginning in 1995. Of the consented patients, US residents who had only CC were included in this study. The majority of patients in the registry are from the surrounding five-state region. Pathologic tumor staging was carried out using the tumor-node-metastasis (TNM) system. Collection of tissue material and review of patient records to obtain clinical information were performed under IRB-approved protocol.
Mismatch Repair Status
Defective DNA mismatch repair (dMMR) was defined by the presence of microsatellite instability (MSI-H) and/or the absence of protein expression for hMLH1. Proficient DNA mismatch repair (pMMR) was defined by the absence of microsatellite instability (MSS/MSI-L) and the presence of normal protein expression for hMLH1. MSI in CC cases was compared with paired normal and tumor DNA isolated (Qiagen DNA extraction kit) from formalin-fixed, paraffin-embedded (FFPE) material. Tumors were classified as MSI-H if > 30% of markers demonstrated instability and as MSS/MSI-L if < 30% demonstrated MSI [22, 23]. Immunohistochemical analysis of hMLH1 expression was performed on FFPE samples, as previously described .
Following harvest, tissue samples were immediately snap frozen and stored at -80°C. For each case, frozen tumor tissue was cut on a cryostat to generate a hematoxylin and eosin (H&E) stained slide. Tumors were then evaluated for content of tumor present, areas of tumor containing at last 70% tumor or greater were macro-dissected and tissue equivalent to 7 mm2 and 10-microns thick was then sectioned and placed in a vial containing 400 uL of RLT buffer (QIAGEN, Chatsworth, CA) including 4 μL of β-mercaptoethanol. The vial was then stored at -80°C until utilized for RNA extraction using TRIzol® LSTrizol© (Invitrogen, Corp., Carlsbad, CA), according to the manufacturer's instructions.
miRNA Expression by qRT-PCR
qRT-PCR was carried using normalization to U6 snRNA. First strand cDNA was synthesized from 1 μg of total RNA using a miScript reverse transcription kit (Qiagen, Valencia, CA). miRNA was quantified with an miRscript SYBR Green PCR kit (Catl: 218073, Qiagen, Valencia, CA) using cDNA equivalent of 50 ng total RNA per reaction. Mature miRNA-specific forward primers [see Additional file 1] were purchased from a vendor (IDT, USA) and the universal reverse primer provided by the manufacturer. Real-time PCR was performed at 55°C annealing following standard protocol of the manufacturer in 7500 Real Time PCR system (Applied Biosystem, USA) and threshold cycles (C T) were calculated using Sequence Detection Software (SDS v1.2.1, Applied Biosystem, USA). Fold expression was calculated from the triplicate of C T values following the 2-ΔΔCt method .
The chemistry involved in the miRNA BeadArray™ was similar to that used in the DASL amplification process . Total RNA [200 ng (all 108 samples) and 400 ng (38 technical replicates)] was poly-adenylated and then converted to cDNA using a biotinylated oligo-dT primer with a universal PCR sequenced at its 5'-end. This was followed by annealing of a miRNA-specific oligonucleotide pool (MSO), which consisted of three parts: a universal PCR priming site at the 5' end; an address sequence complementary to a capture sequence on the BeadArray™; and a microRNA-specific sequence at the 3'-end. Extension of MSO was facilitated by addition of a polymerase but only if their 3' bases were complementary to the cognate sequence in the cDNA template. Common primers were used to amplify the cDNA templates; the primer complimentary to the BeadArray™ was fluorescently labeled. The single-stranded PCR product was hybridized to the Sentrix Array Matrix where the labeled strand bound to the bead on the array containing the complementary address sequence. The arrays were imaged using an Illumina BeadArray™ Reader, which measures the fluorescence intensity at each addressed bead location. Intensity files were analyzed using BeadStudio version 3.1.1. Expression levels were expressed as an average Beta value.
Pre-processing and Analyses
Data were exported from Illumina's BeadStudio package, with neither background correction nor scale correction, and transformed to the log base 2 scale. Quantile normalization  was used as implemented by the R library affy [26, 27]. Samples were distributed on two Illumina SAMs; systematic effects due to SAM were evaluated and found to be minimal. A linear mixed-effects model  was used on a probe-by-probe basis to compare expression between groups accounting for the correlation between replicates or paired normal-tumor samples from a given subject. Replicate samples were averaged for further data analyses.
Support Vector Machine Classification (SVM) was carried out using the Genedata expressionist data analyses suite. Clustering was carried out on log transformed median centered data using Cluster 3.0 for Mac OS X (C Clustering Library 1.43) centroid linkage clustering and was visualized using Java TreeView Version 1.1.3. miRNA tumor expression data has been submitted to NCBI's Gene Expression Omnibus (GEO) as ID: GSE18392. Additionally, the data will be searchable on our website http://www.oncomir.umn.edu.
Genomewide CC methylation data , genes implicated in CC via forward genetic screen  and stem cell miRNA expression data  were obtained from publication. The Diana miRpath application  was used to determine pathway enrichment within the putative targets for miRNAs found to be altered in CC. miRNA::mRNA interactions were downloaded from the miRDB version 2.0 database. Custom Perl scripts were used to parse miRDB, identify predicted miRNA interactions with tumor suppressor mRNA and identify relationships between miRNA and methylated regions.
Sample counts of tissues meeting specified characteristics analyzed in this study
dMMR/MSIH tumors (sporadic)
pMMR Stage II tumor
pMMR Stage III tumor
pMMR Stage IV tumor
pMMR Stage undefined
Technical Array Replicates Show High Level of Correlation
miRNA Profiles Show Distinct Separation of Tumors from Normal Tissue
Unsupervised hierarchical clustering was carried out on the miRNA expression data (Figure 1D). Analyses of the dendrogram revealed that normal colon tissue profiles cluster together in a group separate from the colon tumor samples (Figure 1D). In addition, normal tissue gave a tighter internal correlation, signified by higher correlations, than the tumor tissue subset, which as a group showed less internal correlation.
SVM analyses were used to determine whether the miRNA profiled by the unsupervised hierarchical clustering are sufficient to classify the tumor/normal status. SVMs look for distinctive patterns in groups of data and then use them to classify group membership status for additional data. A Monte Carlo cross-validation simulation was set up to create training sets from 50% of the data and use the sets to predict tumor/normal status of the remainder of the dataset with 1000 iterations. Following combination of the results from all simulations, tumor vs normal status was correctly classified 100% of the time. The confusion matrix is included [see Additional file 4].
Using these criteria, 39 miRNAs were differentially expressed between normal colon tissue and the composite collection of colon tumor tissue (Figure 2A). Many of these miRNA show similar expression patterns and are found in close chromosomal proximity. A subset of 17 miRNAs showed increased expression in colon tumors relative to normal colon tissue. Importantly, miR-135b displayed the largest average change, a 4.55-fold increase. miR-96, miR-182, miR-182* and miR-183 were all up-regulated; expression was highly correlated and mapped to the same region of chromosome 7. Additionally, 22 miRNAs showed decreased expression in colon tumors. Among the down-regulated miRNAs, miR-133a and miR-1, as well as miR-30a-3p and miR-30a-5p, have a common chromosomal locus on chromosome 6 and chromosome 18 respectively and also display highly correlated expression levels. Further, miR-20b and miR-363 show similar levels and are expressed from a common region on chromosome X.
miRNAs Show Striking Differences with Mismatch Repair Status
Direct comparison of dMMR with pMMR tumors by linear effects mixed model using the Bonferroni correction revealed that six miRNAs were differentially expressed, with high significance, between different tumor subtypes (Figure 2B; detailed statistics provided [see Additional file 5]). miRNAs that showed decreased levels in pMMR relative to dMMR tumors included miR-552, miR-592, miR-181c and miR-196b. In contrast, miR-625 and miR-31 exhibited increased levels in dMMR relative to pMMR tumors.
As expected, unsupervised hierarchical clustering, (Figure 2C) as well as SVM analyses, were able to demonstrate separation of dMMR from pMMR tumors using statistically significant miRNAs as defined in Figure 2A-B. However, separation was not as complete as that observed between normal and tumor tissues using all probes (Figure 1D). Examination of the hierarchy generated revealed that a pMMR tumor subset gave expression patterns resembling those found in dMMR tumors, when analyzing statistically significant miRNAs. The area for which SVM was unable to correctly define the tumor MMR status was also the boundary region observed by unsupervised hierarchical clustering between dMMR and pMMR tumors. The SVM confusion matrix for this analysis is included [see Additional file 4]. Further, unsupervised clustering was not able to differentiate between pMMR stages II, III and IV using either the complete dataset or the statistically significant tumor dataset. Furthermore, no differences were observed between tumors from males and females.
dMMR vs pMMR Significant miRNAs Show Differential Levels of Response Relative to Normal Colon Tissue
Examination of miRNAs showing highly significant changes in tumor vs normal revealed that the expression levels for 6 miRNAs were different (p < 0.05 and fold change > 1.5) between dMMR and pMMR tumors (Figure 3B). miR-1, miR-133a, miR-328 and miR-9 are further decreased in dMMR tumors relative to pMMR tumors. miR-224 was induced slightly in dMMR and further increased in pMMR tumor. miR-486 was reduced only in dMMR tumors; and HS-29 was induced to a greater extent in dMMR tumors than pMMR tumors.
While the majority of miRNAs that were differentially expressed in pMMR tumors showed similar perturbations between stages, stage IV pMMR tumors displayed a trend towards increased levels of a subset of miRNAs, relative to stage II and stage III pMMR tumors. Six of the miRNAs identified as highly significant to tumors (HS-29, miR-135b, miR-32, miR-33, miR-542-5p and miR-96) displayed higher expression in pMMR stage IV as compared to stage II tumors (p < 0.05 and fold change > 1.5) (Figure 3C). The levels found in Stage III tumors trended towards those identified in Stage II tumors.
Chromosomal Regions Containing Differentially Expressed miRNAs Also Show Alterations in DNA Methylation Status
Fishers Exact test to determine probability of overlap between differentially expressed miRNAs and differentially methylated regions in colon cancer and the regions of overlap
Regions with Δ miRNA transcript
Regions with Δ methylation
Regions with no Δ methylation
Validation of miRNA Expression
Predicted Targets for miRNAs Significantly Altered in Colon Cancer are Enriched in Cancer-related Signaling Pathways
Induced miRNAs are Predicted to Target Tumor Suppressors
Recently, the results of a forward genetic screen for tumor suppressors in mice independently validated the role of a number of genes in intestinal tumorigenesis . Therefore, we sought to determine whether miRNAs up-regulated in colon tumors were predicted to interact with the 3' UTR of these tumor suppressors, potentially attenuating their transcript level/translation status. Of the 15 genes described as "drivers" of CC , interactions were found between eight genes and 8 miRNAs reported to be up-regulated in colon tumors in the current study. Furthermore, in four of these cases, five separate up-regulated miRNAs were predicted to bind to the 3' UTR of the tumor suppressor (Figure 5B).
Similarities in Colon Cancer and Stem Cell Differentiation
This study has revealed that a number of miRNAs are strongly differentially expressed during the development of CC, including miRNAs not previously reported. Using the Illumina platform, 39 miRNAs were identified that show highly statistically significant and meaningful fold change alterations in CC tumors. The core of this response showed very similar patterns in all colon tumor types studied. In addition, six miRNAs were identified that are significantly differentially expressed in dMMR tumors compared to pMMR tumors.
As previously described, high-level correlations in the miRNA data derived from both technical and biological replicate experiments were observed over a wide range of detection [21, 33]. In addition, the impact of 2-fold fluctuations in the input RNA population was less than the impact of biological variability between normal samples. As expected, therefore, distinct expression signatures were observed in the current study between normal and colon tumor tissues. Furthermore, unsupervised hierarchical clustering and SVM results showed that this was robust and highly predictive in nature.
Several groups have published miRNA profiles of colon tumors using different platforms with or without an amplification step: RT-PCR for 150 miRNAs , microRNA microarray for 389 miRNAs , and LNA-based oligonucleotide arrays for ~315 miRNAs . The present study is thus the most comprehensive, evaluating 735 miRNAs, the most statistically rigorous, using a Bonferroni correction to the multiple testing problem requiring large fold changes and used a macrodissected set of tumor tissues, therefore minimizing the effect from non-malignant cells. Most notably, 18 of the 39 miRNAs we found significantly altered in CC, have been previously reported by systematic RT-PCR .
Previous work with spotted arrays and qRT-PCR revealed the induction of miR-21 [17, 19] and the repression of miR-143  in CC compared to normal colon samples. Surprisingly, these miRNA did not show differential expression on the Illumina platform. We conducted separate qRT-PCR for these miRNAs and were able to see induction of miR-21 and repression of miR-143 in our tumor set compared to normal tissue. miR-21 showed a significant P-value (P = 2.436e-5 but had an average fold change of only 1.02 between groups. Non-responsive probes for miR-21 and miR-143 all had very high expression levels and very low variance, while responsive probes had higher variance and lower average values.
In order to understand the rationale for this discrepancy, we plotted the miRNA expression profiles relative to the average of normal tissue, alongside the absolute raw expression levels for each probe [see Additional file 6]. Our results demonstrate that miRNAs with high raw intensity values as well as a negative control, show a tight distribution of values whereas that distribution is much more variable in responsive probes. Further comparison reveals that non-responsive probes have very high expression levels and very low variance, while responsive probes have higher variance and lower average values. These effects are also evident following normalization. Filtering the expression dataset for probes with high average expression and low variance reveals 28 probes with these features [see Additional file 7].
Our observation of high expression "silent" probes is consistent with a probe cDNA hybridization model where all probe-binding sites are occupied. According to this model, further increases or decreases in miRNA levels will not be visible by array analyses due to binding saturation. This creates a ceiling above which any change, either increase or decrease in miRNA level, will remain undetectable. Additional evidence in support of this model can be found by looking at the relative concentrations of these miRNA in deep sequencing. As an example, in some cases miR-21 made up 50% of the pool in deep sequencing .
Comparison of arrays generated with both the Illumina platform and spotted arrays for 4 different commonly used cell lines showed that following removal of high average expression low variance "silent" Illumina probes, and low signal probes from the cDNA arrays, the remaining probes cluster together independently of the array platform [see Additional File 8]. This information coupled with our qRT-PCR verification of responsive miRNA, indicated that the Illumina arrays have very low false positive rates; but are potentially susceptible to false negative signal events.
dMMR and pMMR tumors showed differential expression of a small set of miRNA. The molecular etiology of those tumors involving dMMR is very heterogeneous, involving several different genes and numerous mechanisms of gene inactivation, including epigenetic, somatic and germline alterations. Among sporadic CC, the vast majority of cases with dMMR are due to inactivation of MLH1 (~95%), with MSH2 and MSH6 accounting for ~5% and < 1%, respectively . For MLH1, the most common mechanism (~90% of cases) of gene inactivation is promoter hypermethylation . In this study, analyses of dMMR cases were specifically confined to those with loss of MLH1 due to promoter hypermethylation. This was done to achieve a homogeneous group. Thus, the results derived from this study primarily reflect the biology of sporadic MLH1 CC.
Furthermore, hMLH1 methylation-associated microsatellite instability has also been strongly associated with tumors that express the CpG Island Methylator Phenotype (CIMP) . As hMLH1 methylation-associated microsatellite instability generally does not occur among sporadic cases outside the context of CIMP, it appears that the underlying basis for mismatch repair deficiency among this select group of sporadic colon cancer is a broader epigenetic control defect that affects hMLH1 in some, but not all CIMP tumors. CIMP tumors represent another subset of all CC. Thus, the few pMMR samples that closely resemble the dMMR subset as a group, with respect to the miRNA profile, may have an underlying CIMP phenotype, which would be common to both groups. In fact, all but one of the cases, for which the CIMP phenotype was available, cluster within a single group containing both dMMR and pMMR.
Collectively, our data supports a model for colon tumorigenesis encompassing miRNA::mRNA interactions. In support of this model, the predicted mRNA target lists compiled from the 39 altered miRNAs in CC are enriched in tumorigenesis and cancer-related signaling pathways. The simplest explanation of this phenomenon is that up-regulated miRNAs directly or indirectly decrease expression of tumor suppressor proteins in contrast to down-regulated miRNAs that may lead to increased oncogene expression.
The predicted interactions were further explored between up-regulated miRNAs and a set of "drivers" of CC as defined by a forward genetic screen in mice . Many of the induced miRNA found in this study were predicted to interact with the tumor suppressors. In addition, PTEN , SMAD4 , and NOTCH1 , are known tumor suppressors whose transcript or protein levels are decreased in tumors. miRNA-mediated decreases in tumor suppressors provide an elegant explanation for the observed tumor expression patterns.
In this analysis, highly significant increases were observed in miR-135b in CC and the interaction between miR-135b and APC was identified as relevant. It has recently been shown that increased levels of miR-135a/b lead directly to decreased protein expression of the CC tumor suppressor APC via a direct binding interaction between miR-135b and APC mRNA 3' UTR . Although APC mutations are found in a majority of CC, deregulation of miR-135b may have an adverse effect on APC in the remainder of cases.
Inadvertent expression and activation of tumor suppressors could easily disrupt normal growth. High levels of intestinal cell proliferation are required to offset the very high turnover rate of intestinal tissue . We propose that miRNA mediated attenuation of transcription/translation of tumor suppressors is a necessary step in normal intestinal development allowing for cell proliferation. With differentiation, the intestinal epithelial cells no longer replicate, which coincides with a change in miRNA expression and increased levels of tumor suppressors. miRNA expression profiles found in CC show striking similarities with those miRNA profiles found in undifferentiated embryonic stem cells relative to differentiated stem cells (Figure 6A) and cancer cell lines . Embryonic stem cells and colon tumor cells are both capable of almost unlimited mitotic division. A potential interpretation of this observation is that the intestinal crypt cells fail to properly differentiate and may instead continue to actively divide leading to tumor formation (Figure 6B). This is further supported by the observed decreases in miR-1 and miR-133a, which are involved in maintaining the differentiation status of muscle cells . In addition, miR-1 over-expression has been shown to cause expression of differentiated muscle cell mRNA . The decreased levels of miR-1 and miR-133a in dMMR tumors relative to pMMR tumors may also explain why sporadic CC with dMMR show poor differentiation.
Temporal patterns of gene expression during cell-cell adhesion-initiated polarization of cultured human Caco-2 cells found similar transcript changes to those observed during migration and differentiation of intestinal epithelial cells in vivo, despite the absence of morphogen gradients and interactions with stromal cells characteristic of enterocyte differentiation in situ . Here, we propose that miRNAs may be involved in these changes and that the miRNA profile within each cell is modified as cells stop dividing, undergo chromatin remodeling and become further differentiated with limited ability to replicate.
Taken together, these findings suggest that the inability to properly differentiate due to loss of epigenetic control may be an important factor in the development of colon cancer. The role of induced miRNA found in tumors may be to (i) drive tumorigenesis by attenuating the translation of tumor suppressors thereby maintaining a state capable of cell division; or (ii) act as passengers whose levels reflect that they are locked into a state of uncontrolled cell replication.
Several basic questions arise from these findings. How and why are miRNA expression levels altered and how can they be modified? Our analyses show that miR-135b was highly up-regulated in colon tumors. PDRM5 is a tumor suppressor and a target of epigenetic silencing in CC . Of interest, PDRM5 has been shown to bind to the promoter of miR-135b and PDRM5 silencing leads to increased levels of miR-135b. Mechanistically, PRDM5 recruits HDAC1 and G9a to the miR-135b promoter, and this chromatin remodeling results in decreased miR-135b expression. Decreases in PDRM5 expression levels lead to increased expression of miR-135b presumably due to loss of recruitment of this chromatin-remodeling complex. Intriguingly, PDRM5 was also found by CHIP-CHIP to be present at the promoter regions of miR-9, miR-378, miR-196b, miR-96-182-183, as well as additional colon cancer relevant miRNA miR-21 and miR-143-145 . Decreases in PDRM5 were associated with both increases and decreases in miRNA expression. miR-9 was also down-regulated in invasive breast cancers via promoter hypermethylation . Therefore, PDRM5-mediated chromatin remodeling may be a general epigenetic mechanism related to a wide range of tumors. Drugs that modify chromatin structure and methylation status in a sequence specific fashion may be useful in colon cancer treatment due to the reversible nature of histone modification and DNA methylation. Those that modify miRNA transcript level or mimic/inhibit miRNA function may also be useful.
The data presented here unifies miRNA expression platforms, broadens the miRNAome associated with CC and defines highly significant differences in miRNA expression between dMMR and pMMR tumors. We hypothesize that the differential miRNA expression patterns observed may modify oncogene and tumor suppressor protein expression and thereby may promote tumorigenesis by blocking the proper differentiation of intestinal progenitor cells, consequently leading to improper cell division and tumor formation. By identifying the miRNA alterations associated with CC at a high level of definition in a statistically powerful dataset, we have provided an important template on which to discover critical biomarkers and develop more effective treatment protocols for colon cancer at all stages.
We thank the Minnesota Supercomputing Institute for providing access to computational resources.
Funding Supported in part by grant PTDC/SAU-GMG/099161/2008 from Fundação para a Ciência e a Tecnologia (FCT), Lisbon, Portugal (to C.M.P.R.), and by Ph.D. fellowship SFRH/BD/24165/2005 (to P.M.B.).
This work was supported by a grant from the Minnesota Partnership for Biotechnology and Medical Genomics [Medica Foundation 2006-100412].
- Ferlay J, Autier P, Boniol M, Heanue M, Colombet M, Boyle P: Estimates of the cancer incidence and mortality in Europe in 2006. Ann Oncol. 2007, 18: 581-592. 10.1093/annonc/mdl498.View ArticlePubMedGoogle Scholar
- Jemal A, Siegel R, Ward E, Hao Y, Xu J, Murray T, Thun MJ: Cancer statistics. CA: Cancer J Clin. 2008, 58: 71-96. 10.3322/CA.2007.0010.Google Scholar
- Baudhuin LM, Burgart LJ, Leontovich O, Thibodeau SN: Use of microsatellite instability and immunohistochemistry testing for the identification of individuals at risk for Lynch Syndrome. Fam Cancer. 2005, 4 (3): 255-265. 10.1007/s10689-004-1447-6.View ArticlePubMedGoogle Scholar
- Grady WM, Carethers JM: Genomic and epigenetic instability in colorectal cancer pathogenesis. Gastroenterology. 2008, 135: 1079-1099. 10.1053/j.gastro.2008.07.076.View ArticlePubMedPubMed CentralGoogle Scholar
- Lindor NM, Burgart LJ, Leontovich O, Goldberg RM, Cunningham JM, Sargent DJ, Walsh-Vockley C, Petersen GM, Walsh MD, Leggett BA, et al: Immunohistochemistry versus microsatellite instability testing in phenotyping colorectal tumors. J Clin Onc. 2002, 20: 1043-1048. 10.1200/JCO.20.4.1043.View ArticleGoogle Scholar
- Zauber AG, Levin TR, Jaffe CC, Galen BA, Ransohoff DF, Brown ML: Implications of new colorectal cancer screening technologies for primary care practice. Med Care. 2008, 46: S138-S146. 10.1097/MLR.0b013e31818192ef.View ArticlePubMedPubMed CentralGoogle Scholar
- Wolpin BM, Mayer RJ: Systemic treatment of colorectal cancer. Gastroenterology. 2008, 134: 1296-1310. 10.1053/j.gastro.2008.02.098.View ArticlePubMedPubMed CentralGoogle Scholar
- Kehoe J, Khatri VP: Staging and prognosis of colon cancer. Surg Oncol Clin N Am. 2006, 15: 129-146. 10.1016/j.soc.2005.08.006.View ArticlePubMedGoogle Scholar
- Benson AB: New approaches to the adjuvant therapy of colon cancer. Oncologist. 2006, 11: 973-980. 10.1634/theoncologist.11-9-973.View ArticlePubMedGoogle Scholar
- Chan SK, Griffith OL, Tai IT, Jones SJ: Meta-analysis of colorectal cancer gene expression profiling studies identifies consistently reported candidate biomarkers. Cancer Epidemiol Biomarkers Prev. 2008, 17: 543-552. 10.1158/1055-9965.EPI-07-2615.View ArticlePubMedGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- Carthew RW: Gene regulation by microRNAs. Curr Opin Genet Dev. 2006, 16: 203-208. 10.1016/j.gde.2006.02.012.View ArticlePubMedGoogle Scholar
- Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006, 6: 259-269. 10.1038/nrc1840.View ArticlePubMedGoogle Scholar
- Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, et al: A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci USA. 2006, 103: 2257-2261. 10.1073/pnas.0510565103.View ArticlePubMedPubMed CentralGoogle Scholar
- Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838. 10.1038/nature03702.View ArticlePubMedGoogle Scholar
- Rosenfeld N, Aharonov R, Meiri E, Rosenwald S, Spector Y, Zepeniuk M, Benjamin H, Shabes N, Tabak S, Levy A, et al: MicroRNAs accurately identify cancer tissue origin. Nat Biotechnol. 2008, 26: 462-469. 10.1038/nbt1392.View ArticlePubMedGoogle Scholar
- Bandres E, Cubedo E, Agirre X, Malumbres R, Zarate R, Ramirez N, Abajo A, Navarro A, Moreno I, Monzo M, et al: Identification by real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues. Mol Cancer. 2006, 5: 29-10.1186/1476-4598-5-29.View ArticlePubMedPubMed CentralGoogle Scholar
- Schepeler T, Reinert JT, Ostenfeld MS, Christensen LL, Silahtaroglu AN, Dyrskjot L, Wiuf C, Sorensen FJ, Kruhoffer M, Laurberg S, et al: Diagnostic and prognostic microRNAs in stage II colon cancer. Cancer Res. 2008, 68: 6416-6424. 10.1158/0008-5472.CAN-07-6110.View ArticlePubMedGoogle Scholar
- Schetter AJ, Leung SY, Sohn JJ, Zanetti KA, Bowman ED, Yanaihara N, Yuen ST, Chan TL, Kwong DL, et al: MicroRNA expression profiles associated with prognosis and therapeutic outcome in colon adenocarcinoma. JAMA. 2008, 299: 425-436. 10.1001/jama.299.4.425.PubMedPubMed CentralGoogle Scholar
- Slaby O, Svoboda M, Fabian P, Serdova T, Knoflickova D, Bednarikova M, Nenutil R, Vyzula R: Altered expression of miR-21, miR-31, miR-143 and miR-145 is related to clinicopathologic features of colorectal cancer. Oncology. 2007, 72: 397-402. 10.1159/000113489.View ArticlePubMedGoogle Scholar
- Chen J, Lozach J, Garcia EW, Barnes B, Luo S, Mikoulitch I, Zhou L, Schroth G, Fan JB: Highly sensitive and specific microRNA expression profiling using BeadArray technology. Nucleic Acids Res. 2008, 36: e87-10.1093/nar/gkn387.View ArticlePubMedPubMed CentralGoogle Scholar
- Boland CR, Thibodeau SN, Hamilton SR, Sidransky D, Eshleman JR, Burt RW, Meltzer SJ, Rodriguez-Bigas MA, Fodde R, Ranzani GN, et al: A National Cancer Institute Workshop on Microsatellite Instability for cancer detection and familial predisposition: development of international criteria for the determination of microsatellite instability in colorectal cancer. Cancer Res. 1998, 58: 5248-5257.PubMedGoogle Scholar
- Thibodeau SN, French AJ, Cunningham JM, Tester D, Burgart LJ, Roche PC, McDonnell SK, Schaid DJ, Vockley CW, Michels VV, et al: Microsatellite instability in colorectal cancer: different mutator phenotypes and the principal involvement of hMLH1. Cancer Res. 1998, 58: 1713-1718.PubMedGoogle Scholar
- Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008, 3: 1101-1108. 10.1038/nprot.2008.73.View ArticlePubMedGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.View ArticlePubMedGoogle Scholar
- Gautier L, Cope L, Bolstad BM, Irizarry RA: affy--analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 307-315. 10.1093/bioinformatics/btg405.View ArticlePubMedGoogle Scholar
- R Development Core Team R: A language and environment for statistical computing. 2008, R Foundation for Statistical Computing, Vienna, AustriaGoogle Scholar
- Pinheiro JC, Douglas MB: Mixed-effects models in S and S-PLUS. 2002, Springer-Verlag New York, LLCGoogle Scholar
- Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, et al: The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nature Genetics. 2009, 41: 178-186. 10.1038/ng.298.View ArticlePubMedPubMed CentralGoogle Scholar
- Starr TK, Allaei R, Silverstein KA, Staggs RA, Sarver AL, Bergemann TL, Gupta M, O'Sullivan MG, Matise I, Dupuy AJ, et al: A transposon-based genetic screen in mice identifies genes altered in colorectal cancer. Science. 2009, 323: 1747-1750. 10.1126/science.1163040.View ArticlePubMedPubMed CentralGoogle Scholar
- Laurent LC, Chen J, Ulitsky I, Mueller FJ, Lu C, Shamir R, Fan JB, Loring JF: Comprehensive microRNA profiling reveals a unique human embryonic stem cell signature dominated by a single seed sequence. Stem Cells. 2008, 26: 1506-1516. 10.1634/stemcells.2007-1081.View ArticlePubMedGoogle Scholar
- Papadopoulos GL, Alexiou P, Maragkakis M, Reczko M, Hatzigeorgiou AG: DIANA-mirPath: Integrating human and mouse microRNAs in pathways. Bioinformatics. 2009, 25: 1991-1993. 10.1093/bioinformatics/btp299.View ArticlePubMedGoogle Scholar
- Cunningham JM, Oberg AL, Borralho PM, Kren BT, French AJ, Wang L, Bots BM, Morlan BW, Silverstein KAT, Staggs R, et al: Evaluation of a new high-dimensional miRNA profiling platform. BMC Med Genomics. 2009, 2: 57-10.1186/1755-8794-2-57.View ArticlePubMedPubMed CentralGoogle Scholar
- Cunningham JM, Kim CY, Christensen ER, Tester DJ, Parc Y, Burgart LJ, Halling KC, McDonnell SK, Schaid DJ, Vockley CW, et al: The frequency of hereditary defective mismatch repair in a prospective series of unselected colorectal carcinomas. Am J Human Genet. 2001, 69: 780-790. 10.1086/323658.View ArticleGoogle Scholar
- Cunningham JM, Kim CY, Christensen ER, Tester DJ, Parc Y, Burgart LJ, Halling KC, McDonnell SK, Schaid DJ, Walsh Vockley C, et al: Hypermethylation of the hMLH1 promoter in colon cancer with microsatellite instability. Cancer Res. 1998, 58: 3455-3460.PubMedGoogle Scholar
- Weisenberger DJ, Siegmund KD, Campan M, Young J, Long TI, Faasse MA, Kang GH, Widschwendter M, Weener D, Buchanan D, et al: CpG island methylator phenotype underlies sporadic microsatellite instability and is tightly associated with BRAF mutation in colorectal cancer. Nat Genet. 2006, 38: 787-793. 10.1038/ng1834.View ArticlePubMedGoogle Scholar
- Zhou X, Loukola A, Salovaara R, Nystrom-Lahti M, Peltomaki P, de la Chapelle A, Aaltonen LA, Eng C: PTEN mutational spectra, expression levels, and subcellular localization in microsatellite stable and unstable colorectal cancers. Am J Pathol. 2002, 161: 439-447.View ArticlePubMedPubMed CentralGoogle Scholar
- Salovaara R, Roth S, Loukola A, Launonen V, Sistonen P, Avizienyte E, Kristo P, Jarvinen H, Souchelnytskyi S, Sarlomo-Rikala M, et al: Frequent loss of SMAD4/DPC4 protein in colorectal cancers. Gut. 2002, 51: 56-59. 10.1136/gut.51.1.56.View ArticlePubMedPubMed CentralGoogle Scholar
- Nicolas M, Wolfer A, Raj K, Kummer JA, Mill P, van Noort M, Hui CC, Clevers H, Dotto GP, Radtke F: Notch1 functions as a tumor suppressor in mouse skin. Nat Genet. 2003, 33: 416-421. 10.1038/ng1099.View ArticlePubMedGoogle Scholar
- Nagel R, le Sage C, Diosdado B, Waal van der M, Oude Vrielink JA, Bolijn A, Meijer GA, Agami R: Regulation of the adenomatous polyposis coli gene by the miR-135 family in colorectal cancer. Cancer Res. 2008, 68: 5795-5802. 10.1158/0008-5472.CAN-08-0951.View ArticlePubMedGoogle Scholar
- Radtke F, Clevers H: Self-renewal and cancer of the gut: Two sides of a coin. Science. 2005, 307: 1904-1909. 10.1126/science.1104815.View ArticlePubMedGoogle Scholar
- Sweetman D, Goljanek K, Rathjen T, Oustanina S, Braun T, Dalmay T, Munsterberg A: Specific requirements of MRFs for the expression of muscle specific microRNAs, miR-1, miR-206 and miR-133. Dev Biol. 2008, 321: 491-499. 10.1016/j.ydbio.2008.06.019.View ArticlePubMedGoogle Scholar
- Lim LP, Lau NC, Garret-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs down-regulate large numbers of target mRNAs. Nature. 2005, 433: 769-773. 10.1038/nature03315.View ArticlePubMedGoogle Scholar
- Sääf AM, Halbleib JM, Chen X, Yuen ST, Leung SY, Nelson WJ, Brown PO: Parallels between global transcriptional programs of polarizing Caco-2 intestinal epithelial cells in vitro and gene expression programs in normal colon and colon cancer. Mol Biol Cell. 2007, 18: 4245-4260. 10.1091/mbc.E07-04-0309.View ArticlePubMedPubMed CentralGoogle Scholar
- Watanabe Y, Toyota M, Kondo Y, Suzuki H, Imai T, Ohe-Toyota M, Maruyaa R, Nohima M, Sasaki Y, Sekido Y, et al: PRDM5 identified as a target of epigenetic silencing in colorectal and gastric cancer. Clin Cancer Res. 2007, 13: 4786-4794. 10.1158/1078-0432.CCR-07-0305.View ArticlePubMedGoogle Scholar
- Duan Z, Person RE, Lee HH, Huang S, Donadieu J, Badolato R, Grimes HL, Papayannopoullou T, Horwitz MS: Epigenetic regulation of protein-coding and microRNA genes by the Gfi1-interacting tumor suppressor PRDM5. Mol Cell Biol. 2007, 27: 6889-6902. 10.1128/MCB.00762-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Lehmann U, Hasemeier B, Christgen M, Muller M, Romermann D, Langer F, Kreipe H: Epigenetic inactivation of microRNA gene hsa-mir-9-1 in human breast cancer. J Pathol. 2008, 214: 17-24. 10.1002/path.2251.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/9/401/prepub
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.