Estrogen imprinting is used to describe a phenomenon in which early developmental exposure to endocrine disruptors increases breast cancer risk later in adult life. We propose that long-lived, self-regenerating stem and progenitor cells are more susceptible to the exposure injury than terminally differentiated epithelial cells in the breast duct. Mammospheres, containing enriched breast progenitors, were used as an exposure system to simulate this imprinting phenomenon in vitro. Using MeDIP-chip, a methylation microarray screening method, we found that 0.5% (120 loci) of human CpG islands were hypermethylated in epithelial cells derived from estrogen-exposed progenitors compared with the non–estrogen-exposed control cells. This epigenetic event may lead to progressive silencing of tumor suppressor genes, including RUNX3, in these epithelial cells, which also occurred in primary breast tumors. Furthermore, normal tissue in close proximity to the tumor site also displayed RUNX3 hypermethylation, suggesting that this aberrant event occurs in early breast carcinogenesis. The high prevalence of estrogen-induced epigenetic changes in primary tumors and the surrounding histologically normal tissues provides the first empirical link between estrogen injury of breast stem/progenitor cells and carcinogenesis. This finding also offers a mechanistic explanation as to why a tumor suppressor gene, such as RUNX3, can be heritably silenced by epigenetic mechanisms in breast cancer.
The dysregulation of transforming growth factor-β (TGF-β) signaling plays a crucial role in ovarian carcinogenesis and in maintaining cancer stem cell properties. Classified as a member of the ATP-binding cassette (ABC) family, ABCA1 was previously identified by methylated DNA immunoprecipitation microarray (mDIP-Chip) to be methylated in ovarian cancer cell lines, A2780 and CP70. By microarray, it was also found to be upregulated in immortalized ovarian surface epithelial (IOSE) cells following TGF-β treatment. Thus, we hypothesized that ABCA1 may be involved in ovarian cancer and its initiation.
We first compared the expression level of ABCA1 in IOSE cells and a panel of ovarian cancer cell lines and found that ABCA1 was expressed in HeyC2, SKOV3, MCP3, and MCP2 ovarian cancer cell lines but downregulated in A2780 and CP70 ovarian cancer cell lines. The reduced expression of ABCA1 in A2780 and CP70 cells was associated with promoter hypermethylation, as demonstrated by bisulfite pyro-sequencing. We also found that knockdown of ABCA1 increased the cholesterol level and promoted cell growth in vitro and in vivo. Further analysis of ABCA1 methylation in 76 ovarian cancer patient samples demonstrated that patients with higher ABCA1 methylation are associated with high stage (P = 0.0131) and grade (P = 0.0137). Kaplan-Meier analysis also found that patients with higher levels of methylation of ABCA1 have shorter overall survival (P = 0.019). Furthermore, tissue microarray using 55 ovarian cancer patient samples revealed that patients with a lower level of ABCA1 expression are associated with shorter progress-free survival (P = 0.038).
ABCA1 may be a tumor suppressor and is hypermethylated in a subset of ovarian cancer patients. Hypermethylation of ABCA1 is associated with poor prognosis in these patients.
Electronic supplementary material
The online version of this article (doi:10.1186/s13148-014-0036-2) contains supplementary material, which is available to authorized users.
Ovarian cancer; Epigenetics; ABCA1
Ulcerative colitis (UC) is characterized by chronic inflammation of the colon. During inflammation, NF-κB is increased in colonic epithelial cells and in immune cells, leading to increases in proinflammatory cytokines. These events then increase DNA methyltransferases (DNMTs), which silence a subset of tumor suppressor genes by promoter methylation. Negative regulators of the Wnt pathway are frequently methylated in UC, leading to dysregulation of the pathway and, potentially, to colorectal cancer. We determined if black raspberries (BRBs) influence promoter methylation of suppressors in the Wnt pathway in dextran sodium sulfate (DSS)-induced UC. C57BL/6J mice received 1% DSS and were fed either control or 5% BRB diets. Mice were euthanized on days 7, 14 and 28, and their colons, spleen and bone marrow were collected. Berries reduced ulceration at day 28. This was accompanied by decreased staining of macrophages and neutrophils and decreased NF-κB p65 nuclear localization in the colon at all time points. At day 7, BRBs demethylated the promoter of dkk3, leading to its increased messenger RNA (mRNA) expression in colon, spleen and bone marrow. β-Catenin nuclear localization, c-Myc staining as well as protein expression of DNMT3B, histone deacetylases 1 and 2 (HDAC1 and HDAC2) and methyl-binding domain 2 (MBD2) were all decreased in colon; mRNA expression of these four proteins was decreased in bone marrow cells by BRBs. These results suggest that BRBs suppress colonic ulceration by correcting promoter hypermethylation of suppressor genes in the colon, as well as in the spleen and bone marrow that systematically regulate inflammation.
Summary: Our results suggest that dietary BRBs suppress colonic ulceration by correcting promoter hypermethylation of suppressor genes in the colon, as well as in the spleen and bone marrow that systematically regulate inflammation in DSS-induced UC.
Ulcerative colitis (UC) is frequently an intermediate step to colon cancer. The interleukin-10 knock-out (KO) mouse is a genetic model of this progression. We report that KO mice fed 5% black raspberries (BRBs) had significantly less colonic ulceration as compared to KO mice that consumed the control diet. Dysfunction of the Wnt signaling pathway is a key event in UC-associated colon carcinogenesis. Therefore, we investigated the effects of BRBs on the Wnt pathway and found that the BRB-fed KO mice exhibited a significantly lower level of β-catenin nuclear translocation. We followed-up this observation by evaluating the effect of BRBs on selected Wnt pathway antagonists. The mRNA expression levels of wif1, sox17, and qki were diminished in the KO mice, while they were expressed at normal levels in KO mice fed BRBs. The lower mRNA expression of these genes in colon from the KO mice correlated with hypermethylation of their promoter regions; BRBs decreased their promoter methylation and increased mRNA expression of these genes. This hypomethylation was associated with elevated protein expression of key proteins/enzymes that augment methylation, e.g., dnmt3b, hdac1, hdac2, and mbd2 in the KO mice; and BRBs decreased protein expression of these proteins/enzymes. The KO mouse model recapitulates what occurs in human UC. Promoter methylation of CDH1 and SFRP1 was significantly higher in human UC tissues compared to their adjacent normal tissues. In conclusion, our results suggest that BRBs inhibit colonic ulceration and ultimately colon cancer partly through inhibiting aberrant epigenetic events that dysregulate Wnt signaling.
black raspberries; Wnt pathway; methylation; colon carcinogenesis; IL-10 KO mice
There is emerging evidence that Wnt pathway activity may increase during the progression from colorectal adenoma to carcinoma and that this increase is potentially an important step towards the invasive stage. Here, we investigated whether epigenetic silencing of Wnt antagonists is the biological driver for this increased Wnt activity in human tissues and how these methylation changes correlate with MSI (Microsatelite Instability) and CIMP (CpG Island Methylator Phenotype) statuses as well as known mutations in genes driving colorectal neoplasia.
We conducted a systematic analysis by pyrosequencing, to determine the promoter methylation of CpG islands associated with 17 Wnt signaling component genes. Methylation levels were correlated with MSI and CIMP statuses and known mutations within the APC, BRAF and KRAS genes in 264 matched samples representing the progression from normal to pre-invasive adenoma to colorectal carcinoma.
We discovered widespread hypermethylation of the Wnt antagonists SFRP1, SFRP2, SFRP5, DKK2, WIF1 and SOX17 in the transition from normal to adenoma with only the Wnt antagonists SFRP1, SFRP2, DKK2 and WIF1 showing further significant increase in methylation from adenoma to carcinoma. We show this to be accompanied by loss of expression of these Wnt antagonists, and by an increase in nuclear Wnt pathway activity. Mixed effects models revealed that mutations in APC, BRAF and KRAS occur at the transition from normal to adenoma stages whilst the hypermethylation of the Wnt antagonists continued to accumulate during the transitions from adenoma to carcinoma stages.
Our study provides strong evidence for a correlation between progressive hypermethylation and silencing of several Wnt antagonists with stepping-up in Wnt pathway activity beyond the APC loss associated tumour-initiating Wnt signalling levels.
Electronic supplementary material
The online version of this article (doi:10.1186/1471-2407-14-891) contains supplementary material, which is available to authorized users.
DNA methylation is a common epigenetic marker that regulates gene expression. A robust and cost-effective way for measuring whole genome methylation is Methyl-CpG binding domain-based capture followed by sequencing (MBDCap-seq). In this study, we proposed BIMMER, a Hidden Markov Model (HMM) for differential Methylation Regions (DMRs) identification, where HMMs were proposed to model the methylation status in normal and cancer samples in the first layer and another HMM was introduced to model the relationship between differential methylation and methylation statuses in normal and cancer samples. To carry out the prediction for BIMMER, an Expectation-Maximization algorithm was derived. BIMMER was validated on the simulated data and applied to real MBDCap-seq data of normal and cancer samples. BIMMER revealed that 8.83% of the breast cancer genome are differentially methylated and the majority are hypo-methylated in breast cancer.
DNA methylation; differential methylation; MBDCap-seq; Hidden Markov Model (HMM)
Genetic amplification, mutation, and translocation are known to play a causal role in the up-regulation of an oncogene in cancer cells. Here we report an emerging role of microRNA, the epigenetic deregulation of which may also lead to this oncogenic activation. SOX4, an oncogene belonging to the SRY-related high-mobility-group-box family, was found to be over-expressed (P<0.005) in endometrial tumors (n=74) compared with uninvolved controls (n=20). This gene is computationally predicted to be the target of a microRNA, miR-129-2. When compared to the matched endometria, the expression of miR-129-2 was lost in 27 of 31 primary endometrial tumors that also showed a concomitant gain of SOX4 expression (P<0.001). This inverse relationship is associated with hypermethylation of the miR-129-2 CpG island, which was observed in endometrial cancer cell lines (n=6) and 68% of 117 endometrioid endometrial tumors analyzed. Reactivation of miR-129-2 in cancer cells by pharmacological induction of histone acetylation and DNA demethylation resulted in decreased SOX4 expression. In addition, restoration of miR-129-2 by cell transfection led to decreased SOX4 expression and reduced proliferation of cancer cells. Further analysis found a significant correlation of hypermethylated miR-129-2 with microsatellite instability and MLH1 methylation status (P<0.001), and poor overall survival (P<0.039) in patients. Therefore, these results imply that the aberrant expression of SOX4 is, in part, caused by epigenetic repression of miR-129-2 in endometrial cancer. Unlike the notion that promoter hypomethylation may up-regulate an oncogene, we present a new paradigm in which hypermethylation-mediated silencing of a microRNA de-represses its oncogenic target in cancer cells.
DNA methylation; epigenetics; miR-129-2; SOX4
Early exposure to xenoestrogens may predispose to breast cancer risk later in adult life. It is likely that long-lived, self-regenerating epithelial progenitor cells are more susceptible to these exposure injuries over time and transmit the injured memory through epigenetic mechanisms to their differentiated progeny. Here, we used progenitor-containing mammospheres as an in vitro exposure model to study this epigenetic effect. Expression profiling identified that, relative to control cells, 9.1% of microRNAs (82 of 898 loci) were altered in epithelial progeny derived from mammospheres exposed to a synthetic estrogen, diethylstilbestrol. Repressive chromatin marks, trimethyl Lys27 of histone H3 (H3K27me3) and dimethyl Lys9 of histone H3 (H3K9me2), were found at a down-regulated locus, miR-9-3, in epithelial cells preexposed to diethylstilbestrol. This was accompanied by recruitment of DNA methyltransferase 1 that caused an aberrant increase in DNA methylation of its promoter CpG island in mammosphere-derived epithelial cells on diethylstilbestrol preexposure. Functional analyses suggest that miR-9-3 plays a role in the p53-related apoptotic pathway. Epigenetic silencing of this gene, therefore, reduces this cellular function and promotes the proliferation of breast cancer cells. Promoter hypermethylation of this microRNA may be a hallmark for early breast cancer development, and restoration of its expression by epigenetic and microRNA-based therapies is another viable option for future treatment of this disease.
Ovarian cancer ranks the most lethal among gynecologic neoplasms in women. To develop potential bio-markers for diagnosis, we have identified five novel genes (CYP39A1, GTF2A1, FOXD4L4, EBP, and HAAO) that are hypermethylated in ovarian tumors, compared with the non-malignant normal ovarian surface epithelia, using the quantitative methylation-specific polymerase chain reactions. Interestingly enough, multivariate Cox regression analysis has identified hypermethylation of CYP39A1 correlated with an increase rate of relapsing (P=0.032, hazard ratio >1). Concordant hypermethylation in at least three loci was observed in 50 out of 55 (91%) of ovarian tumors examined. The test sensitivity and specificity were assessed to be 96 and 67% for CYP39A1; 95 and 88% for GTF2A1; 93 and 67% for FOXD4L4; 81 and 67% for EBP; 89 and 82% for HAAO, respectively. Our data have identified, for the first time, GTF2A1 alone, or GTF2A1 plus HAAO are excellent candidate biomarkers for detecting this disease. Moreover, the known functions of these gene products further implicate dys-regulated transcriptional control, cholesterol metabolism, or synthesis of quinolinic acids, may play important roles in attributing to ovarian neoplasm. Molecular therapies, by reversing the aberrant epigenomes using inhibitory agents or by abrogating the upstream signaling pathways that convey the epigenomic perturbations, may be developed into promising treatment regimens.
ovarian cancer; epigenetics; DNA methylation; biomarkers; quantitative methylation-specific polymerase chain reaction
The interplay between histone modifications and promoter hypermethylation provides a causative explanation for epigenetic gene silencing in cancer. Less is known about the upstream initiators that direct this process. Here, we report that the Cystatin M (CST6) tumor suppressor gene is concurrently down-regulated with other loci in breast epithelial cells co-cultured with cancer-associated fibroblasts (CAFs). Promoter hypermethylation of CST6 is associated with aberrant AKT1 activation in epithelial cells, as well as the disabled INNP4B regulator resulted from the suppression by CAFs. Repressive chromatin, marked by trimethyl-H3K27 and dimethyl-H3K9, and de novo DNA methylation is established at the promoter. The findings suggest that microenvironmental stimuli are triggers in this epigenetic cascade, leading to the long-term silencing of CST6 in breast tumors. Our present findings implicate a causal mechanism defining how tumor stromal fibroblasts support neoplastic progression by manipulating the epigenome of mammary epithelial cells. The result also highlights the importance of direct cell-cell contract between epithelial cells and the surrounding fibroblasts that confer this epigenetic perturbation. Since this two-way interaction is anticipated, the described co-culture system can be used to determine the effect of epithelial factors on fibroblasts in future studies.
p57Kip2, a cyclin-dependent kinase inhibitor, is considered to be a candidate tumor suppressor gene that has been implicated in Beckwith-Wiedemann syndrome and sporadic cancers. In addition, decreased expression of p57Kip2 protein has been frequently observed in pancreatic, lung, breast, bladder, gastrointestinal tract and prostate cancers. However, p57Kip2 gene mutations are rare in these cancers suggesting that other unknown mechanisms might be at play in reducing its expression. The aim of this study was to investigate the molecular mechanism of down-regulation of p57Kip2 in prostate cancer.
We observed a significant negative correlation between the expression of p57Kip2 and microRNA-21 (miR-21) in prostate cancer samples and after androgen deprivation with castration in the CWR22 human prostate cancer xenograft model. We report that miR-21 targeted the coding region and decreased p57Kip2 mRNA and protein levels in prostate cancer cells. Conversely, inhibition of endogenous miR-21 by an anti-miR-21 inhibitor strongly induced p57Kip2 expression. Furthermore, we found that knockdown of p57Kip2 reversed the effects of the anti-miR-21 inhibitor on cell migration and anchorage-independent cell growth.
Our results indicate that miR-21 is able to downregulate p57Kip2 expression by targeting the coding region of the gene and is also able to attenuate p57Kip2 mediated functional responses. This is the first report demonstrating that p57Kip2 is a novel target of miR-21 in prostate cancer and revealing a novel oncogenic function of this microRNA.
Electronic supplementary material
The online version of this article (doi:10.1186/1476-4598-13-212) contains supplementary material, which is available to authorized users.
p57Kip2; microRNA-21 and prostate cancer
DNA motifs are short sequences varying from 6 to 25 bp and can be highly variable and degenerated. One major approach for predicting transcription factor (TF) binding is using position weight matrix (PWM) to represent information content of regulatory sites; however, when used as the sole means of identifying binding sites suffers from the limited amount of training data available and a high rate of false-positive predictions. ChIPMotifs program is a de novo motif finding tool developed for ChIP-based high-throughput data, and W-ChIPMotifs is a Web application tool for ChIPMotifs. It composes various ab initio motif discovery tools such as MEME, MaMF, Weeder and optimizes the significance of the detected motifs by using bootstrap re-sampling error estimation and a Fisher test. Using these techniques, we determined a PWM for OCT4 which is similar to canonical OCT4 consensus sequence. In a separate study, we also use de novo motif discovery to suggest that ZNF263 binds to a 24-nt site that differs from the motif predicted by the zinc finger code in several positions.
Motif; ChIP; Position weight matrix; OCT4; ZNF263
A causal role of gene amplification in tumorigenesis is well-known, while amplification of DNA regulatory elements as an oncogenic driver remains unclear. In this study, we integrated next-generation sequencing approaches to map distant estrogen response elements (DEREs) that remotely control transcription of target genes through chromatin proximity. Two densely mapped DERE regions located on chromosomes 17q23 and 20q13 were frequently amplified in ERα-positive luminal breast cancer. These aberrantly amplified DEREs deregulated target gene expression potentially linked to cancer development and tamoxifen resistance. Progressive accumulation of DERE copies was observed in normal breast progenitor cells chronically exposed to estrogenic chemicals. These findings may extend to other DNA regulatory elements, the amplification of which can profoundly alter target transcriptome during tumorigenesis.
The enhancer pioneer transcription factor FoxA1 is a global mediator of steroid receptor (SR) action in hormone-dependent cancers. In castration-resistant prostate cancer (CRPC), FoxA1 acts as an androgen receptor co-factor to drive G2-M phase cell-cycle transit. Here we describe a mechanistically distinct SR-independent role for FoxA1 in driving G1-S phase cell-cycle transit in CRPC. By comparing FoxA1 binding sites in prostate cancer cell genomes, we defined a co-dependent set of FoxA1-MYBL2 and FoxA1-CREB1 binding sites within the regulatory regions of the Cyclin E2 and E2F1 genes that are critical for CRPC growth. Binding at these sites upregulate the Cyclin E2 and Cyclin A2 genes in CRPC but not in earlier stage androgen-dependent prostate cancer (ADPC), establishing a stage-specific role for this pathway in CRPC growth. Mechanistic investigations indicated that FoxA1, MYBL2 or CREB1 induction of histone H3 acetylation facilitated nucleosome disruption as the basis for co-dependent transcriptional activation and G1-S phase cell-cycle transit. Our findings establish FoxA1 as a pivotal driver of the cell-cycle in CRPC which promotes G1-S phase transit as well as G2-M phase transit through two distinct mechanisms.
FoxA1; cistrome; cell-cycle G1-S progression; castration-resistant prostate cancer
Aberrant expression of SOX4 in endometrial cancer has been identified and partially was contributed to hypermethylation of miR-129-2. Other miRNAs are suspected to influence SOX 4 as well. The current study seeks to identify other hypermethylated miRNAs that regulate SOX4 in endometrial carcinomas.
Methylation levels of miRNA promoter regions were measured by combined bisulfite restriction analysis (COBRA) and pyrosequencing assays. Gene expression was determined by RT-qPCR. Methylation level of a miRNA locus was corrected with clinicopathologic factors for 252 gynecological specimens.
In silico analysis identified 13 miRNA loci bound on the 3′-UTR of SOX4. Using COBRA assays, increased methylation of miR-203, miR-219-2, miR-596, and miR-618 was detected in endometrial cancer cells relative to those seen in a normal cell line and in normal endometrium. Transfection of a miR-203 mimic decreased SOX4 gene expression. Hypermethylation of miR-203 was detected in 52% of type I endometrioid endometrial carcinomas (n=131) but was not seen in any of 10 uninvolved normal endometria (P<0.001). Methylation status of miR-203 was significantly associated with microsatellite instability and MLH1 methylation in endometrial tumors (P<0.001). Furthermore, hypermethylation of miR-203 was found in endometrioid and clear endometrial subtype tumors, but not in cervical squamous cell and ovarian carcinomas.
Hypermethylation of miR-203 is a frequent event in endometrial carcinomas and is strongly associated with microsatellite instability and MLH1 methylation status. Thus, miR-203 methylation level might represent a marker for patients with endometrioid and clear endometrial sub-cancers.
Endometrial carcinoma; DNA methylation; SOX4; miR-203
Promoter hypermethylation associated tumor suppressor genes (TSGs) silencing has been explored as a therapeutic target for hypomethylating agents. Promoter methylation change may serve as a pharmacodynamic endpoint for evaluation of the efficacy of these agents and predict the patient’s clinical response. Herein, a LC-MS/MS assay has been developed for quantitative regional DNA methylation analysis using the molar ratio of 5-methyl-2′-deoxycytidine (5mdC) to 2′-deoxycytidine (2dC) in the enzymatic hydrolysate of fully methylated bisulfite-converted PCR amplicons as the methylation indicator. The assay can differentiate 5% of promoter methylation level with an intra-day precision ranging from 3.00 to 16.0% using two TSGs: HIN-1 and RASSF1A. This method was applied to characterize decitabine-induced promoter DNA methylation changes of these two TSGs in a breast cancer MCF-7 cell line. Promoter methylation of these TSGs was found to decrease in a dose-dependent manner. Correspondingly, the expression of these TSGs was enhanced. The sensitivity and reproducibility of the method make it a valuable tool for specific gene methylation analysis, which could aid characterization of hypomethylating activity on specific genes by hypomethylating agents in a clinical setting.
Regional DNA Methylation; LC-MS/MS; Quantification
In prostate cancer, androgen receptor (AR) binding and androgen-responsive gene expression are defined by hormone-independent binding patterns of the pioneer factors FoxA1 and GATA2. Insufficient evidence of the mechanisms by which GATA2 contributes to this process precludes complete understanding of a key determinant of tissue-specific AR activity. Our observations suggest that GATA2 facilitates androgen-responsive gene expression by three distinct modes of action. By occupying novel binding sites within the AR gene locus, GATA2 positively regulates AR expression before and after androgen stimulation. Additionally, GATA2 engages AR target gene enhancers prior to hormone stimulation, producing an active and accessible chromatin environment via recruitment of the histone acetyltransferase p300. Finally, GATA2 functions in establishing and/or sustaining basal locus looping by recruiting the Mediator subunit MED1 in the absence of androgen. These mechanisms may contribute to the generally positive role of GATA2 in defining AR genome-wide binding patterns that determine androgen-responsive gene expression profiles. We also find that GATA2 and FoxA1 exhibit both independent and codependent co-occupancy of AR target gene enhancers. Identifying these determinants of AR transcriptional activity may provide a foundation for the development of future prostate cancer therapeutics that target pioneer factor function.
Alternative splicing (AS), in higher eukaryotes, is one of the mechanisms of post-transcriptional regulation that generate multiple transcripts from the same gene. One particular mode of AS is the skipping event where an exon may be alternatively excluded or constitutively included in the resulting mature mRNA. Both transcript isoforms from this skipping event site, i.e. in which the exon is either included (inclusion isoform) or excluded (skipping isoform), are typically present in one cell, and maintain a subtle balance that is vital to cellular function and dynamics. However, how the prevailing conditions dictate which isoform is expressed and what biological factors might influence the regulation of this process remain areas requiring further exploration. In this study, we have developed a novel computational method, graph-based exon-skipping scanner (GESS), for de novo detection of skipping event sites from raw RNA-seq reads without prior knowledge of gene annotations, as well as for determining the dominant isoform generated from such sites. We have applied our method to publicly available RNA-seq data in GM12878 and K562 cells from the ENCODE consortium and experimentally validated several skipping site predictions by RT-PCR. Furthermore, we integrated other sequencing-based genomic data to investigate the impact of splicing activities, transcription factors (TFs) and epigenetic histone modifications on splicing outcomes. Our computational analysis found that splice sites within the skipping-isoform-dominated group (SIDG) tended to exhibit weaker MaxEntScan-calculated splice site strength around middle, ‘skipping’, exons compared to those in the inclusion-isoform-dominated group (IIDG). We further showed the positional preference pattern of splicing factors, characterized by enrichment in the intronic splice sites immediately bordering middle exons. Finally, our analysis suggested that different epigenetic factors may introduce a variable obstacle in the process of exon–intron boundary establishment leading to skipping events.
Targeting tumor metabolism by energy restriction-mimetic agents (ERMAs) has emerged as a strategy for cancer therapy/prevention. Evidence suggests a mechanistic link between ERMA-mediated antitumor effects and epigenetic gene regulation.
Microarray analysis showed that a novel thiazolidinedione-derived ERMA, CG-12, and glucose deprivation could suppress DNA methyltransferase (DNMT)1 expression and reactivate DNA methylation-silenced tumor suppressor genes in LNCaP prostate cancer cells. Thus, we investigated the effects of a potent CG-12 derivative, CG-5, vis-à-vis 2-deoxyglucose, glucose deprivation and/or 5-aza-deoxycytidine, on DNMT isoform expression (Western blotting, RT-PCR), DNMT1 transcriptional activation (luciferase reporter assay), and expression of genes frequently hypermethylated in prostate cancer (quantitative real-time PCR). Promoter methylation was assessed by pyrosequencing analysis. SiRNA-mediated knockdown and ectopic expression of DNMT1 were used to validate DNMT1 as a target of CG-5.
CG-5 and glucose deprivation upregulated the expression of DNA methylation-silenced tumor suppressor genes, including GADD45a, GADD45b, IGFBP3, LAMB3, BASP1, GPX3, and GSTP1, but also downregulated methylated tumor/invasion-promoting genes, including CD44, S100A4, and TACSTD2. In contrast, 5-aza-deoxycytidine induced global reactivation of these genes. CG-5 mediated these epigenetic effects by transcriptional repression of DNMT1, which was associated with reduced expression of Sp1 and E2F1. SiRNA-mediated knockdown and ectopic expression of DNMT1 corroborated DNMT1's role in the modulation of gene expression by CG-5. Pyrosequencing revealed differential effects of CG-5 versus 5-aza-deoxycytidine on promoter methylation in these genes.
These findings reveal a previously uncharacterized epigenetic effect of ERMAs on DNA methylation-silenced tumor suppressor genes, which may foster novel strategies for prostate cancer therapy.
energy restriction-mimetic agent; prostate cancer; energy restriction; DNA methyltransferases; epigenetics
Aberrant DNA methylation of CpG islands, CpG island shores and first exons is known to play a key role in the altered gene expression patterns in all human cancers. To date, a systematic study on the effect of DNA methylation on gene expression using high resolution data has not been reported. In this study, we conducted an integrated analysis of MethylCap-sequencing data and Affymetrix gene expression microarray data for 30 breast cancer cell lines representing different breast tumor phenotypes. As well-developed methods for the integrated analysis do not currently exist, we created a series of four different analysis methods. On the computational side, our goal is to develop methylome data analysis protocols for the integrated analysis of DNA methylation and gene expression data on the genome scale. On the cancer biology side, we present comprehensive genome-wide methylome analysis results for differentially methylated regions and their potential effect on gene expression in 30 breast cancer cell lines representing three molecular phenotypes, luminal, basal A and basal B. Our integrated analysis demonstrates that methylation status of different genomic regions may play a key role in establishing transcriptional patterns in molecular subtypes of human breast cancer.
One big limitation of computational tools for analyzing ChIP-seq data is that most of them ignore non-unique tags (NUTs) that match the human genome even though NUTs comprise up to 60% of all raw tags in ChIP-seq data. Effectively utilizing these NUTs would increase the sequencing depth and allow a more accurate detection of enriched binding sites, which in turn could lead to more precise and significant biological interpretations. In this study, we have developed a computational tool, LOcating Non-Unique matched Tags (LONUT), to improve the detection of enriched regions from ChIP-seq data. Our LONUT algorithm applies a linear and polynomial regression model to establish an empirical score (ES) formula by considering two influential factors, the distance of NUTs to peaks identified using uniquely matched tags (UMTs) and the enrichment score for those peaks resulting in each NUT being assigned to a unique location on the reference genome. The newly located tags from the set of NUTs are combined with the original UMTs to produce a final set of combined matched tags (CMTs). LONUT was tested on many different datasets representing three different characteristics of biological data types. The detected sites were validated using de novo motif discovery and ChIP-PCR. We demonstrate the specificity and accuracy of LONUT and show that our program not only improves the detection of binding sites for ChIP-seq, but also identifies additional binding sites.
DNA methylation of promoter CpG islands is associated with gene suppression, and its unique genome-wide profiles have been linked to tumor progression. Coupled with high-throughput sequencing technologies, it can now efficiently determine genome-wide methylation profiles in cancer cells. Also, experimental and computational technologies make it possible to find the functional relationship between cancer-specific methylation patterns and their clinicopathological parameters.
Cancer methylome system (CMS) is a web-based database application designed for the visualization, comparison and statistical analysis of human cancer-specific DNA methylation. Methylation intensities were obtained from MBDCap-sequencing, pre-processed and stored in the database. 191 patient samples (169 tumor and 22 normal specimen) and 41 breast cancer cell-lines are deposited in the database, comprising about 6.6 billion uniquely mapped sequence reads. This provides comprehensive and genome-wide epigenetic portraits of human breast cancer and endometrial cancer to date. Two views are proposed for users to better understand methylation structure at the genomic level or systemic methylation alteration at the gene level. In addition, a variety of annotation tracks are provided to cover genomic information. CMS includes important analytic functions for interpretation of methylation data, such as the detection of differentially methylated regions, statistical calculation of global methylation intensities, multiple gene sets of biologically significant categories, interactivity with UCSC via custom-track data. We also present examples of discoveries utilizing the framework.
CMS provides visualization and analytic functions for cancer methylome datasets. A comprehensive collection of datasets, a variety of embedded analytic functions and extensive applications with biological and translational significance make this system powerful and unique in cancer methylation research. CMS is freely accessible at: http://cbbiweb.uthscsa.edu/KMethylomes/.
An emerging Hi-C protocol has the ability to probe three-dimensional (3D) architecture and capture chromatin interactions in a genome-wide scale. It provides informative results to address how chromatin organization changes contribute to disease/tumor occurrence and progression in response to stimulation of environmental chemicals or hormones.
In this study, using MCF7 cells as a model system, we found estrogen stimulation significantly impact chromatin interactions, leading to alteration of gene regulation and the associated histone modification states. Many chromosomal interaction regions at different levels of interaction frequency were identified. In particular, the top 10 hot regions with the highest interaction frequency are enriched with breast cancer specific genes. Furthermore, four types of E2-mediated strong differential (gain- or loss-) chromosomal (intra- or inter-) interactions were classified, in which the number of gain-chromosomal interactions is less than the number of loss-chromosomal interactions upon E2 stimulation. Finally, by integrating with eight histone modification marks, DNA methylation, regulatory elements regions, ERα and Pol-II binding activities, associations between epigenetic patterns and high chromosomal interaction frequency were revealed in E2-mediated gene regulation.
The work provides insight into the effect of chromatin interaction on E2/ERα regulated downstream genes in breast cancer cells.
Recent genome-wide profiling reveals highly complex regulation networks among ERα and its targets. We integrated estrogen (E2)-stimulated time-series ERα ChIP-seq and gene expression data to identify the ERα-centered transcription factor (TF) hubs and their target genes, and inferred the time-variant hierarchical network structures using a Bayesian multivariate modeling approach. With its recurrent motif patterns, we determined three embedded regulatory modules from the ERα core transcriptional network. The GO analyses revealed the distinct biological function associated with each of three embedded modules. The survival analysis showed the genes in each module were able to render a significant survival correlation in breast cancer patient cohorts. In summary, our Bayesian statistical modeling and modularity analysis not only reveals the dynamic properties of the ERα-centered regulatory network and associated distinct biological functions, but also provides a reliable and effective genomic analytical approach for the analysis of dynamic regulatory network for any given TF.
The transcriptional response driven by Hypoxia-inducible factor (HIF) is central to the adaptation to oxygen restriction. Despite recent characterization of genome-wide HIF DNA binding locations and hypoxia-regulated transcripts in different cell types, the molecular bases of HIF target selection remain unresolved. Herein, we combined multi-level experimental data and computational predictions to identify sequence motifs that may contribute to HIF target selectivity. We obtained a core set of bona fide HIF binding regions by integrating multiple HIF1 DNA binding and hypoxia expression profiling datasets. This core set exhibits evolutionarily conserved binding regions and is enriched in functional responses to hypoxia. Computational prediction of enriched transcription factor binding sites identified sequence motifs corresponding to several stress-responsive transcription factors, such as activator protein 1 (AP1), cAMP response element-binding (CREB), or CCAAT-enhancer binding protein (CEBP). Experimental validations on HIF-regulated promoters suggest a functional role of the identified motifs in modulating HIF-mediated transcription. Accordingly, transcriptional targets of these factors are over-represented in a sorted list of hypoxia-regulated genes. Altogether, our results implicate cooperativity among stress-responsive transcription factors in fine-tuning the HIF transcriptional response.