- Methodology article
- Open Access
Discovery pipeline for epigenetically deregulated miRNAs in cancer: integration of primary miRNA transcription
BMC Genomics volume 12, Article number: 54 (2011)
Cancer is commonly associated with widespread disruption of DNA methylation, chromatin modification and miRNA expression. In this study, we established a robust discovery pipeline to identify epigenetically deregulated miRNAs in cancer.
Using an integrative approach that combines primary transcription, genome-wide DNA methylation and H3K9Ac marks with microRNA (miRNA) expression, we identified miRNA genes that were epigenetically modified in cancer. We find miR-205, miR-21, and miR-196b to be epigenetically repressed, and miR-615 epigenetically activated in prostate cancer cells.
We show that detecting changes in primary miRNA transcription levels is a valuable method for detection of local epigenetic modifications that are associated with changes in mature miRNA expression.
MiRNA genes are typically transcribed by RNA polymerase II into primary miRNA (pri-miRNA) transcripts, and transcription appears to be regulated in a similar process as traditional coding genes . Pri-miRNA transcripts are long non-coding RNAs (ncRNA) with stem-loop secondary structures that contain precursor miRNAs (pre-miRNA). Pri-miRNA transcripts are cleaved co-transcriptionally by the enzyme Drosha into stem loop pre-miRNAs, which are exported to the cytoplasm by Exportin 5, and further processed by Dicer into mature miRNA that can be loaded in the RNA induced silencing complex (RISC). Pre-miRNA levels are low relative to pri- and mature miRNAs, suggesting that dicing is an efficient mechanism with little regulation . During development, many pri-miRNAs are expressed but not efficiently processed to mature miRNA . In healthy tissues, the ratio of pri- to mature miRNAs has been shown to be close to one, while in cancer cells, a large number of miRNA genes are transcribed but not processed to mature miRNA . An expanding body of evidence supports a role for miRNAs in disease progression and the potential for epigenetic mechanisms to regulate miRNA expression [5, 6]. Epigenetics is the heritable modification of gene expression without changes in the DNA sequence. The importance of regulated epigenetic information is highlighted by the disruption of multiple epigenetic marks in various disease states, including cancer, which is commonly associated with deregulation of DNA methylation, histone modifications and miRNA expression . Cytosine methylation of CpG islands associated with gene promoters is a well-studied epigenetic mark, and aberrant DNA hypermethylation in combination with altered histone architecture is a common hallmark of neoplastic cells . Histone acetylation at lysine 9 (H3K9Ac) is associated with an active chromatin state and changes in its global patterning are linked to poor prognosis in multiple cancers . MiRNAs have been shown to function as tumor suppressors and oncogenes, through the targeting of essential elements of cellular growth, proliferation and apoptotic pathways . Genomic loci encoding miRNAs can become epigenetically remodelled in cancer, however the extent and effect of these changes in cancer remains to be elucidated .
We argue that a minority of miRNA genes, critical to regulation of cellular growth, proliferation, or survival, gain epigenetic modifications to enable persistent enhancement or inhibition of their activity in order to maintain a neoplastic phenotype. Using an integrative approach combining primary transcription, genome-wide DNA methylation and H3K9Ac patterns with mature miRNA levels, we identified a number of miRNA genes that showed epigenetic remodeling in cancer in combination with de-regulation of both the primary and mature miRNA. This demonstrates that integration of multiple methods of analysis may provide a reliable method for the identification of epigenetically regulated miRNAs.
Results and Discussion
We proposed that epigenetically de-regulated miRNAs in cancer would show three characteristics: first, differential expression at the primary transcript level, second, differential mature miRNA expression, and finally an associated change in one or more epigenetic marks. To eliminate post-transcriptionally regulated miRNAs that showed accumulation of epigenetic marks independent of their transcription, an integrative pipeline that assessed all three criteria was developed (Figure 1).
Normal prostate epithelial (PrEC) and prostate cancer (LNCaP) cells were used as a model system to test this integrative pipeline in the discovery of epigenetically de-regulated miRNAs in cancer. Custom high-density NimbleGen tiling arrays covering miRNA loci from miRBase v12.0 were designed and used to probe total RNA to determine transcription levels, and DNA from methylated-DNA immunoprecipitation (MeDIP) and H3K9Ac immunoprecipitation to investigate epigenetic changes. The oligonucleotide probes of these arrays have an average length of 60nt, precluding interference of the mature 20-23nt miRNAs, thus assaying the full-length primary transcripts. Plotting the array probe intensities against chromosomal coordinates generated high-resolution maps of expression and epigenetic mark enrichment across each locus (plots for miR-205 are shown in Figure 1A & 1B as examples). To generate lists of loci significantly changed in LNCaP cells for each assay, we selected those that were consistently changed in the cancer cell line within 500 bp of the center of the miRNA gene. For our comparisons, we limited the number of miRNAs to those that were assayed by both TaqMan® Low Density Arrays (TLDA) (Figure 1C) and tiling array: 341 miRNAs in total. All changes are articulated as LNCaP with respect to PrEC.
Looking at those loci whose primary transcript repression or activation in LNCaP was associated with a significant change in the mature miRNA, we identified 29 upregulated and 33 downregulated miRNAs (Figure 2A,B). We then looked for concomitant epigenetic changes at these loci in the tiling array data. A high proportion of the upregulated (31%; 9/29) and downregulated loci (33%; 11/33) showed an overlap with at least one of these epigenetic marks (Figure 2C,D). The majority (18/20) of these miRNAs have previously been identified as deregulated in various cancers, and five with prior evidence for epigenetic regulation (Table 1). A concordant change in all assays; expression, DNA methylation, and H3K9Ac restricted the total number of candidate miRNAs to 2/29 up- and 4/33 downregulated (Figure 2C,D). Given the large number of possible epigenetic modifications (for a recent review see ), it is possible that many of the differentially expressed loci that show no change in DNA methylation or H3K9Ac could undergo changes in other epigenetic marks. Furthermore, miRNAs may be regulated by DNA methylation at CpG islands outside the 2000 bp regions interrogated in this analysis.
In order to perform an independent test for miRNAs regulated by DNA methylation, we treated LNCaP and PrEC cells with the drug 5-Aza-CdR, an inhibitor of DNA methylation, and analysed their expression by TLDA. Restricting our list to miRNAs that were significantly deregulated between PrEC and LNCaP we found 30 upregulated by 5-Aza-CdR in PrEC and 16 in LNCaP. Of the miRNA loci showing epigenetic activation by the three criteria of primary transcript expression, DNA methylation loss and H3K9Ac gain, only miR-615 showed 5-Aza-CdR induced reactivation in the normal prostate PrEC cells (Figure 2E). MiR-205, miR-196b and miR-21 satisfied all three criteria for epigenetic repression, and LNCaP 5-Aza-CdR reactivation (Figure 2F).
MiRNA genes showing non-canonical epigenetic regulation were also observed. That is, miRNAs showing up regulation of primary and mature miRNA expression, with a concomitant gain of DNA methylation or loss of H3K9Ac, or vice versa. Twelve miRNAs showed activation and gain of repressive epigenetic marks and six (including the chromosome 14 miRNA cluster) showed repression in combination with active chromatin (Additional file 1, Table S1). None of these miRNAs showed significant up regulation upon treatment with 5-Aza-CdR. Several miRNAs were repressed in the presence of the demethylating drug, suggesting the repression of these loci may be dependent on genes that are themselves subject to regulation by DNA methylation.
To validate the changes in primary miRNA expression, we performed pri-miRNA TaqMan® qPCR on the four miRNAs that met all our criteria. These assays confirmed the repression of pri-miR-205 and pri-miR-21, and over expression of pri-miR-615 in LNCaP cells (Figure 3A,B,D). There was no significant change in pri-miR-196b between LNCaP and PrEC, however treatment with 5-Aza-CdR led to overexpression in PrEC cells (Figure 3C). 5-Aza-CdR led to overexpression of pri-miR-21 and pri-miR-615 in LNCaP, as would be predicted by their DNA methylation status (Figure 3B,D). However, reactivation of pri-miR-205 by 5-Aza-CdR was not significant (Figure 3A).
To verify changes in DNA methylation we analysed recently published methyl-CpG binding domain-based capture (MBDCap) data from LNCaP and PrEC cells . MBDCap preferentially captures DNA methylated double-stranded DNA which is then eluted for high-throughput sequencing . We examined the number of sequencing reads that mapped to within 300 bp of the selected miRNA loci. In agreement with the tiling array data, there was enrichment of sequencing reads at the MIR196b locus (Additional file 1, Figure S1). MIR205 and MIR21 had increased reads in LNCaP over PrEC, however there was also an increase in the input DNA controls. MIR615 showed a low number of sequencing events in both LNCaP and PrEC and input controls (Additional file 1, Figure S1).
The identification of miR-205 as an epigenetically repressed miRNA is consistent with a previous study that examined chromatin modifications in prostate cells using genome-wide promoter arrays, which revealed a gain of the repressive H3K27-trimethylation mark and loss of the active H3K4-trimethylation mark at the MIR205 locus in prostate cancer cells . Furthermore, we recently reported that MIR205 locus is silenced and gains promoter hypermethylation and repressive chromatin marks in invasive bladder tumors and undifferentiated bladder cell lines . Two expression array studies have also identified miR-205 as repressed in prostate cancer [15, 16]. While the repression of pri-miR-205 was validated by qPCR, the lack of significant change in DNA methylation observed in the MBDCap high-throughput sequencing data could be due to differences between the techniques. Both MIR205 and MIR21 reside in regions of low CpG density, whereas MIR196b and MIR615 are located in regions of high CpG density, which may bias the purification efficiencies. Variations between pri-miRNA and array results, such as with pri-miR196b, could be accounted for by differences in assay design. An advantage of the tiling array approach is the high-density coverage of all loci. PCR typically provides improved dynamic range over array analysis, however, the efficiency of a reaction can be limited by location of the primers - as may be the case with pri-miRNA transcripts that have yet to be exhaustively mapped. MiR-21 shows significant epigenetic repression in our array and qPCR analyses, contradictory to its well-established role as an oncogenic miRNA that is overexpressed in cancer , however recent evidence has also shown that miR-21 knockdown elicited no change in proliferative or invasive properties in prostate cancer cells .
The potential for epigenetic regulation of MIR21 has previously been shown in ovarian cancer, where treatment with 5-Aza-CdR leads to miR-21 overexpression , consistent with hypermethylation of the locus. It is possible that epigenetic mechanisms may play a role in context dependent regulation of miR-21 expression, and it would be interesting to examine the prevalence of epigenetic marks at in a cohort of clinical prostate cancer samples. The epigenetic regulation of both MIR615 and MIR196b is particularly intriguing as both genes are encoded from CpG islands located within HOX gene clusters. MIR196b has recently been shown to be hypomethylated in acute lymphoblastic leukemia  and gastric cancer  but to the best of our knowledge, the epigenetic regulation of MIR615 is a novel discovery. MIR615 is located within the HOXC cluster on chromosome 13, which we previously have shown to be epigenetically deregulated in LNCaP cells in long range epigenetically silenced domains . There is evidence that HOX associated miRNAs are essential for development and can function by targeting HOX messenger RNAs . Furthermore, HOX genes have been reported to be involved in abnormal development and malignancy in numerous cancers via both genetic and epigenetic mechanisms .
There are inherent limitations to the various available global epigenetic discovery techniques . Immunoprecipitation based techniques can be affected by antibody sensitivity and MeDIP can be biased towards regions of high CpG content, leading to increased false-negatives at GC-poor sites . In addition, treatment with demethylating drugs such as 5-Aza-CdR can uncover false positives through activation or repression of downstream targets . A multi-platform approach can help rule out many of these false positive results, but the problem of false negatives remains. Here we show that using miRNA primary transcript data enables leverage of the known mechanism of post-transcriptional regulation of miRNAs as an additional means of epigenetic marker discovery. Epigenetic marks such as DNA methylation possess several advantages as biomarkers; DNA is more robust than RNA, can be isolated from formalin fixed paraffin-embedded archival samples and DNA methylation provides a positive readout for genes or miRNAs that are downregulated .
We found that up to 33% of transcriptionally deregulated miRNA loci displayed concordant DNA methylation and H3K9 acetylation patterns. Extending the analysis to additional epigenetic marks, such as the methylation of H3K4, H3K9 and H3K27 and acetylation at H3K14 will increase the number of identified loci. The use of tiling arrays in this study enabled analysis within 2000 bp of a miRNA gene and the increasing availability of high-throughput sequencing technology will allow this approach to be extended genome wide. The prevalence of miRNAs that show contradicting expression and epigenetic marks suggests that there is a subset of genes that are subject to a complex pattern of regulation, where specific marks may recruit activators or repressors to a promoter. We propose that in cancer, a change in primary transcript expression at miRNA loci provides an important first step in the identification of those miRNAs that display variant epigenetic marks. Overlapping these loci with change in mature miRNA restricts analysis to miRNAs that are functionally relevant and that provide possible selective advantage to the cancer cell. Further discovery of DNA hyper- and hypo-methylated miRNAs in cancer using this type of integrative discovery pipeline should assist in identifying loci that may provide reliable novel diagnostic tools.
Cell Lines and Culture
LNCaP cells were propagated in T-media as described previously . PrEC cells were propagated in PrEGM (Lonza, Switzerland). 5-Aza-2'-deoxycytidine (5-Aza-CdR) (Fluka, MO, USA) treatment was performed for 72 hours with fresh media and drug added at 24 and 48 hours. We optimized 5-Aza-CdR treatment of LNCaP and PrEC cells at 3 μM and 1 μM, respectively, based on reactivation of the imprinted gene IGF2 and the prostate cancer hypermethylated gene GSTP1 in response to treatment, as described previously . Total RNA was isolated using Trizol reagent according to supplier's protocol (Invitrogen, CA, USA). Genomic DNA was extracted from untreated LNCaP and PrEC using the Mammalian Blood and Tissue Kit (Qiagen, Germany).
TaqMan® Low-Density Arrays (TLDA)
TLDA v1.0 Early Access (Applied Biosystems, CA, USA) assays for 368 miRNAs were performed as per manufacturer's protocols. In brief, 1 μg total RNA from 5-Aza-CdR treated and untreated LNCaP and PrEC biological duplicates was reverse transcribed with pools of miRNA specific RT-primers, and the resultant cDNA was subjected to real-time reverse transcription quantitative (qRT)-PCR using specific forward primers and TaqMan® fluorescent probes in an array format. RNU48 was selected as the most reliable control for expression normalization. Initial analyses were performed using SDS software (Applied Biosystems). Statistical tests for changes in expression (via Cts) were performed using moderated t-tests without multiple testing corrections using the R limma package . Mature miRNAs were deemed differentially expressed in LNCaP compared to PrEC at p <0.05 and fold change of ≥1.5.
Expression Tiling Arrays
Total RNA (6 μg) from untreated LNCaP and PrEC duplicates was reverse transcribed with random primers and aminoallyl- conjugated dUTP. cDNA products were precipitated and resuspended with mono-reactive Cy3/5 dye (GE Healthcare, UK), and the products were purified with Fairplay columns (Stratagene, TX, USA). Labeling efficiency was determined spectrophotometrically. Hybridizations for untreated LNCaP and PrEC were repeated in two formats to test a method to minimize probe GC bias: competitive hybridization with fluorescently labeled gDNA in the Cy5 channel, and in a dye swap arrangement. The biological repeats were highly reproducible (Pearson correlations: LNCaP 0.94 and PrEC 0.96) and co-hybridization showed no correction to probe binding efficiency (data not shown).
Tiling Array MeDIP-chip/H3K9Ac-ChIP-chip
4 μg of genomic DNA from untreated LNCaP and PrEC biological duplicates was sonicated, then denatured followed by incubation with Anti-5-Methylcytosine antibody (Calibiochem) overnight at 4°C. The antibody/DNA complexes were affinity purified using protein A/G agarose beads (Santa Cruz Biotechnology). Products were washed, eluted and ethanol precipitated. Efficiency of immunoprecipitation was confirmed by PCR at three loci; centromeric chromosome 16, GSTP1 and GAPDH, and compared with no-antibody and input controls, as described previously . Biological duplicates were highly reproducible (Pearson correlation: LNCaP and PrEC 0.96).
LNCaP and PrEC cells were formaldehyde fixed in duplicate prior to SDS lysis and sonication. Immunoprecipitation (IP) was performed overnight at 4°C with an anti-acetyl H3K9 antibody (Millipore, CA, USA). DNA was eluted with 0.1M NaHCO3, 1% SDS elution buffer and reverse cross-linked. IPs and input controls were purified by phenol chloroform extraction and ethanol precipitation, as described previously , duplicate IP reactions were then pooled to ensure sufficient material for array analysis. IP DNA from MeDIP and H3K9Ac were labeled as per the manufacturers protocol (NimbleGen, WI, USA). Briefly, 1 μg IP or input DNA was labeled with Cy-dye random nonamers (Trilink Biotechnologies, CA, USA) using Klenow fragment (New England Biolabs, MA, USA), precipitated, and resuspended in hybridization buffer.
We designed probes for 384 K feature tiling arrays (NimbleGen) to cover the positive strand for the following genomic regions: all human miRNAs listed in mirBASE v12.0 +/-2000 bp; HOX A,B,C and D gene clusters; the chromosome 14 Dlk1 - Dio3 - miRNA cluster region; and hybridization control regions. Array hybridization was performed for 16 hours with 6 μg labeled cDNA with hybribization buffer and alignment oligonucleotides (NimbleGen). Arrays were washed as per the manufacturers protocol and scanned using a GenePix 4000B scanner (Molecular Devices, CA, USA). To screen for miRNA loci that were significantly changed between LNCaP and PrEC cells, we fitted a linear model to the quantile-normalized probe-level data for all probes within 500 bp of the center of a miRNA gene. From the linear model fit, the effect of interest is the average change between PrEC and LNCaP across all probes, after adjusting for probe-specific effects. The t-statistics shown in Table 1 represent the significance of this average difference for each miRNA. MiRNAs were deemed statistically significant using a t-statistic cutoff of ≤-2 or ≥2. Venn diagrams were drawn using Venny , and restricted to overlaps for those miRNAs analysed by TLDA. Results of all array experiments are listed in Additional file 2.
Pri-miRNA Real-Time PCR
We used Applied Biosystems pri-miRNA stem-loop specific TaqMan assays to quantify pri-miRNAs, as per the manufacturers protocol. Briefly, first-strand cDNA was synthesized using 10 ng DNase (RNase free DNase, Qiagen) treated total RNA, or no-template controls, and Superscript reverse transcriptase (Invitrogen). Real-time PCR reactions for each miRNA (10 μl volume) included 4 μl of RT product, were performed in triplicate and included no-template and PCR controls using an ABI7900 PCR system (Applied Biosystems). Reactions performed under the following conditions: 95°C for 10 minutes, followed by 40 cycles at 95°C for 15 seconds and 60°C for 1 minute. Relative quantitative levels of individual miRNAs were determined using 18S for normalization, and statistically significant differences were analysed using the students t-test.
TaqMan Low Density Array
Wiklund ED, Kjems J, Clark SJ: Epigenetic architecture and microRNA: reciprocal regulators. Epigenomics. 2010
Siomi H, Siomi MC: Posttranscriptional regulation of microRNA biogenesis in animals. Molecular cell. 2010, 38: 323-332. 10.1016/j.molcel.2010.03.013.
Thomson JM, Newman M, Parker JS, Morin-Kensicki EM, Wright T, Hammond SM: Extensive post-transcriptional regulation of microRNAs and its implications for cancer. Genes & development. 2006, 20: 2202-2207.
Lee EJ, Baek M, Gusev Y, Brackett DJ, Nuovo GJ, Schmittgen TD: Systematic evaluation of microRNA processing patterns in tissues, cell lines, and tumors. RNA (New York, NY. 2008, 14: 35-42.
Negrini M, Nicoloso MS, Calin GA: MicroRNAs and cancer--new paradigms in molecular oncology. Current opinion in cell biology. 2009, 21: 470-479. 10.1016/j.ceb.2009.03.002.
Guil S, Esteller M: DNA methylomes, histone codes and miRNAs: tying it all together. The international journal of biochemistry & cell biology. 2009, 41: 87-95.
Egger G, Liang G, Aparicio A, Jones PA: Epigenetics in human disease and prospects for epigenetic therapy. Nature. 2004, 429: 457-463. 10.1038/nature02625.
Wilson IM, Davies JJ, Weber M, Brown CJ, Alvarez CE, MacAulay C, Schubeler D, Lam WL: Epigenomics: mapping the methylome. Cell cycle (Georgetown, Tex. 2006, 5: 155-158. 10.4161/cc.5.2.2367.
Cooper CS, Foster CS: Concepts of epigenetics in prostate cancer development. British journal of cancer. 2009, 100: 240-245. 10.1038/sj.bjc.6604771.
Bartels CL, Tsongalis GJ: MicroRNAs: novel biomarkers for human cancer. Clinical chemistry. 2009, 55: 623-631. 10.1373/clinchem.2008.112805.
Robinson MD, Stirzaker C, Statham AL, Coolen MW, Song JZ, Nair SS, Strbenac D, Speed TP, Clark SJ: Evaluation of affinity-based genome-wide DNA methylation data: Effects of CpG density, amplification bias, and copy number variation. Genome research. 2010
Serre D, Lee BH, Ting AH: MBD-isolated Genome Sequencing provides a high-throughput and comprehensive survey of DNA methylation in the human genome. Nucleic acids research. 2010, 38: 391-399. 10.1093/nar/gkp992.
Ke XS, Qu Y, Rostad K, Li WC, Lin B, Halvorsen OJ, Haukaas SA, Jonassen I, Petersen K, Goldfinger N, et al: Genome-wide profiling of histone h3 lysine 4 and lysine 27 trimethylation reveals an epigenetic signature in prostate carcinogenesis. PLoS ONE. 2009, 4: e4687-10.1371/journal.pone.0004687.
Wiklund ED, Bramsen JB, Hulf T, Dyrskjot L, Ramanathan R, Hansen TB, Villadsen SB, Gao S, Ostenfeld MS, Borre M, et al: Coordinated epigenetic repression of the miR-200 family and miR-205 in invasive bladder cancer. International journal of cancer. 2010
Schaefer A, Jung M, Mollenkopf HJ, Wagner I, Stephan C, Jentzmik F, Miller K, Lein M, Kristiansen G, Jung K: Diagnostic and prognostic implications of microRNA profiling in prostate carcinoma. International journal of cancer. 2010, 126: 1166-1176.
Gandellini P, Folini M, Longoni N, Pennati M, Binda M, Colecchia M, Salvioni R, Supino R, Moretti R, Limonta P, et al: miR-205 Exerts tumor-suppressive functions in human prostate through down-regulation of protein kinase Cepsilon. Cancer Res. 2009, 69: 2287-2295. 10.1158/0008-5472.CAN-08-2894.
Ribas J, Lupold SE: The transcriptional regulation of miR-21, its multiple transcripts, and their implication in prostate cancer. Cell Cycle. 2010, 9: 923-929. 10.4161/cc.9.5.10930.
Folini M, Gandellini P, Longoni N, Profumo V, Callari M, Pennati M, Colecchia M, Supino R, Veneroni S, Salvioni R, et al: miR-21: an oncomir on strike in prostate cancer. Molecular cancer. 2010, 9: 12-10.1186/1476-4598-9-12.
Iorio MV, Visone R, Di Leva G, Donati V, Petrocca F, Casalini P, Taccioli C, Volinia S, Liu CG, Alder H, et al: MicroRNA Signatures in Human Ovarian Cancer. Cancer Res. 2007, 67: 8699-8707. 10.1158/0008-5472.CAN-07-1936.
Schotte D, Lange-Turenhout EA, Stumpel DJ, Stam RW, Buijs-Gladdines JG, Meijerink JP, Pieters R, Den Boer ML: Expression of miR-196b is not exclusively MLL-driven but is especially linked to activation of HOXA genes in pediatric acute lymphoblastic leukemia. Haematologica. 2010, 95: 1675-1682. 10.3324/haematol.2010.023481.
Tsai KW, Hu LY, Wu CW, Li SC, Lai CH, Kao HW, Fang WL, Lin Wc: Epigenetic regulation of miR-196b expression in gastric cancer. Genes, Chromosomes and Cancer. 2010, 49: 969-980. 10.1002/gcc.20804.
Coolen MW, Stirzaker C, Song JZ, Statham AL, Kassir Z, Moreno CS, Young AN, Varma V, Speed TP, Cowley M, et al: Consolidation of the cancer genome into domains of repressive chromatin by long-range epigenetic silencing (LRES) reduces transcriptional plasticity. Nature cell biology. 2010, 12: 235-246.
Yekta S, Tabin CJ, Bartel DP: MicroRNAs in the Hox network: an apparent link to posterior prevalence. Nature reviews. 2008, 9: 789-796. 10.1038/nrg2400.
Shah N, Sukumar S: The Hox genes and their roles in oncogenesis. Nat Rev Cancer. 2010
Hayashi H, Nagae G, Tsutsumi S, Kaneshiro K, Kozaki T, Kaneda A, Sugisaki H, Aburatani H: High-resolution mapping of DNA methylation in human genome using oligonucleotide tiling array. Hum Genet. 2007, 120: 701-711. 10.1007/s00439-006-0254-6.
Shi H, Wei SH, Leu YW, Rahmatpanah F, Liu JC, Yan PS, Nephew KP, Huang TH: Triple analysis of the cancer epigenome: an integrated microarray system for assessing gene expression, DNA methylation, and histone acetylation. Cancer Res. 2003, 63: 2164-2171.
Meiers I, Shanks JH, Bostwick DG: Glutathione S-transferase pi (GSTP1) hypermethylation in prostate cancer: review 2007. Pathology. 2007, 39: 299-304. 10.1080/00313020701329906.
Stirzaker C, Song JZ, Davidson B, Clark SJ: Transcriptional gene silencing promotes DNA hypermethylation through a sequential change in chromatin modifications in cancer cells. Cancer Res. 2004, 64: 3871-3877. 10.1158/0008-5472.CAN-03-3690.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-
Song JZ, Stirzaker C, Harrison J, Melki JR, Clark SJ: Hypermethylation trigger of the glutathione-S-transferase gene (GSTP1) in prostate cancer cells. Oncogene. 2002, 21: 1048-1061. 10.1038/sj.onc.1205153.
Coolen MW, Stirzaker C, Song JZ, Statham AL, Kassir Z, Moreno CS, Young AN, Varma V, Speed TP, Cowley M: Consolidation of the cancer genome into domains of repressive chromatin by long-range epigenetic silencing (LRES) reduces transcriptional plasticity. Nature cell biology. 12: 235-246.
Oliveros JC: VENNY. An interactive tool for comparing lists with Venn Diagrams. 2007, [http://bioinfogp.cnb.csic.es/tools/venny/index.html]
Ozen M, Creighton CJ, Ozdemir M, Ittmann M: Widespread deregulation of microRNA expression in human prostate cancer. Oncogene. 2008, 27: 1788-1793. 10.1038/sj.onc.1210809.
Vrba L, Jensen TJ, Garbe JC, Heimark RL, Cress AE, Dickinson S, Stampfer MR, Futscher BW: Role for DNA methylation in the regulation of miR-200c and miR-141 expression in normal and cancer cells. PLoS ONE. 2010, 5: e8697-10.1371/journal.pone.0008697.
Satzger I, Mattern A, Kuettler U, Weinspach D, Voelker B, Kapp A, Gutzmer R: MicroRNA-15b represents an independent prognostic parameter and is correlated with tumor cell proliferation and apoptosis in malignant melanoma. International journal of cancer. 2010, 126: 2553-2562.
Hu Y, Correa AM, Hoque A, Guan B, Ye F, Huang J, Swisher SG, Wu TT, Ajani JA, Xu XC: Prognostic significance of differentially expressed miRNAs in esophageal cancer. International journal of cancer. 2010
Hui AB, Lenarduzzi M, Krushel T, Waldron L, Pintilie M, Shi W, Perez-Ordonez B, Jurisica I, O'Sullivan B, Waldron J, et al: Comprehensive MicroRNA profiling for head and neck squamous cell carcinomas. Clin Cancer Res. 2010, 16: 1129-1139. 10.1158/1078-0432.CCR-09-2166.
Nakano H, Miyazawa T, Kinoshita K, Yamada Y, Yoshida T: Functional screening identifies a microRNA, miR-491 that induces apoptosis by targeting Bcl-X(L) in colorectal cancer cells. International journal of cancer. 2009
Cicatiello L, Mutarelli M, Grober OM, Paris O, Ferraro L, Ravo M, Tarallo R, Luo S, Schroth GP, Seifert M, et al: Estrogen receptor alpha controls a gene network in luminal-like breast cancer cells comprising multiple transcription factors and microRNAs. The American journal of pathology. 2010, 176: 2113-2130. 10.2353/ajpath.2010.090837.
Zhang H, Luo XQ, Zhang P, Huang LB, Zheng YS, Wu J, Zhou H, Qu LH, Xu L, Chen YQ: MicroRNA patterns associated with clinical prognostic parameters and CNS relapse prediction in pediatric acute leukemia. PLoS ONE. 2009, 4: e7826-10.1371/journal.pone.0007826.
Chen Y, Gorski DH: Regulation of angiogenesis through a microRNA (miR-130a) that down-regulates antiangiogenic homeobox genes GAX and HOXA5. Blood. 2008, 111: 1217-1226. 10.1182/blood-2007-07-104133.
Sarver AL, French AJ, Borralho PM, Thayanithy V, Oberg AL, Silverstein KA, Morlan BW, Riska SM, Boardman LA, Cunningham JM, et al: Human colon cancer profiles show differential microRNA expression depending on mismatch repair status and are characteristic of undifferentiated proliferative states. BMC Cancer. 2009, 9: 401-10.1186/1471-2407-9-401.
Bandres E, Agirre X, Bitarte N, Ramirez N, Zarate R, Roman-Gomez J, Prosper F, Garcia-Foncillas J: Epigenetic regulation of microRNA expression in colorectal cancer. International journal of cancer. 2009, 125: 2737-2743. 10.1002/ijc.24638.
Popovic R, Riesbeck LE, Velu CS, Chaubey A, Zhang J, Achille NJ, Erfurth FE, Eaton K, Lu J, Grimes HL, et al: Regulation of mir-196b by MLL and its overexpression by MLL fusions contributes to immortalization. Blood. 2009, 113: 3314-3322. 10.1182/blood-2008-04-154310.
Qin W, Shi Y, Zhao B, Yao C, Jin L, Ma J, Jin Y: miR-24 regulates apoptosis by targeting the open reading frame (ORF) region of FAF1 in cancer cells. PLoS ONE. 2010, 5: e9429-10.1371/journal.pone.0009429.
Valastyan S, Reinhardt F, Benaich N, Calogrias D, Szasz AM, Wang ZC, Brock JE, Richardson AL, Weinberg RA: A pleiotropically acting microRNA, miR-31, inhibits breast cancer metastasis. Cell. 2009, 137: 1032-1046. 10.1016/j.cell.2009.03.047.
Lee EJ, Gusev Y, Jiang J, Nuovo GJ, Lerner MR, Frankel WL, Morgan DL, Postier RG, Brackett DJ, Schmittgen TD: Expression profiling identifies microRNA signature in pancreatic cancer. International journal of cancer. 2007, 120: 1046-1054. 10.1002/ijc.22394.
We thank Warren Kaplan, Mark Cowley and Bruce Tabor for bioinformatic support; Jenny Song and Clare Stirzaker for technical assistance; and Kate Patterson for critical reading of the manuscript. This work was supported by Cure Cancer Australia and New South Wales Cancer Institute Fellowships (to TH), NH&MRC project grants (SJC) and the SIROCCO small RNA EU consortium (EDW).
TH and SJC conceived the project. TH wrote the manuscript, cultured cells, performed array experiments and analysis. TS cultured cells and performed ChIP, MeDIP and TLDA experiments. EDW performed tiling array experiments. SB performed MeDIP experiments. ALS, DS and MDR performed data analyses. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1:Table S1: MiRNAs showing non-canonical epigenetic regulation. Table listing miRNA loci and genes showing significant non-canonical epigenetic regulation. Values are expressed as changes in LNCaP relative to PrEC. Tiling array t-statistics represent the significance of the average difference for each miRNA, with increasing value indicating more significant change. Mature miRNA levels are given as fold change with the associated p-value. Chr: Chromosome, Str: Strand, meth: DNA-methylation, L: LNCaP, P: PrEC, 5Aza: treatment with 5-Aza-CdR, p *: p-value <0.05, Epi./Cancer: publication indicating epigenetic/cancer regulation. Chr14 cluster data includes: MIR127, MIR323, MIR376a-2, MIR432, MIR376a-1, MIR654, MIR382, MIR487b, MIR134, MIR487b, and MIR411. (DOC 122 KB)
Additional file 2:Figure S1: Validation of DNA methylation in selected miRNA loci. Analysis of published DNA methylation enrichment data at selected miRNA loci from Robinson et al.. Graph displays absolute number of MBDCap high throughput sequencing reads within 300 bp of the miRNA locus. (XLS 238 KB)
About this article
Cite this article
Hulf, T., Sibbritt, T., Wiklund, E.D. et al. Discovery pipeline for epigenetically deregulated miRNAs in cancer: integration of primary miRNA transcription. BMC Genomics 12, 54 (2011). https://doi.org/10.1186/1471-2164-12-54
- LNCaP Cell
- miRNA Gene
- Epigenetic Mark
- Tiling Array
- PrEC Cell