Genes suppressed by DNA methylation in non-small cell lung cancer reveal the epigenetics of epithelial–mesenchymal transition

Background

DNA methylation is associated with aberrant gene expression in cancer, and has been shown to correlate with therapeutic response and disease prognosis in some types of cancer. We sought to investigate the biological significance of DNA methylation in lung cancer.

Results

We integrated the gene expression profiles and data of gene promoter methylation for a large panel of non-small cell lung cancer cell lines, and identified 578 candidate genes with expression levels that were inversely correlated to the degree of DNA methylation. We found these candidate genes to be differentially methylated in normal lung tissue versus non-small cell lung cancer tumors, and segregated by histologic and tumor subtypes. We used gene set enrichment analysis of the genes ranked by the degree of correlation between gene expression and DNA methylation to identify gene sets involved in cellular migration and metastasis. Our unsupervised hierarchical clustering of the candidate genes segregated cell lines according to the epithelial-to-mesenchymal transition phenotype. Genes related to the epithelial-to-mesenchymal transition, such as AXL, ESRP1, HoxB4, and SPINT1/2, were among the nearly 20% of the candidate genes that were differentially methylated between epithelial and mesenchymal cells. Greater numbers of genes were methylated in the mesenchymal cells and their expressions were upregulated by 5-azacytidine treatment. Methylation of the candidate genes was associated with erlotinib resistance in wild-type EGFR cell lines. The expression profiles of the candidate genes were associated with 8-week disease control in patients with wild-type EGFR who had unresectable non-small cell lung cancer treated with erlotinib, but not in patients treated with sorafenib.

Conclusions

Our results demonstrate that the underlying biology of genes regulated by DNA methylation may have predictive value in lung cancer that can be exploited therapeutically.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1079) contains supplementary material, which is available to authorized users.     

Epigenetic natural variation in Arabidopsis thaliana

Cytosine methylation of repetitive sequences is widespread in plant genomes, occurring in both symmetric (CpG and CpNpG) as well as asymmetric sequence contexts. We used the methylation-dependent restriction enzyme McrBC to profile methylated DNA using tiling microarrays of Arabidopsis Chromosome 4 in two distinct ecotypes, Columbia and Landsberg erecta. We also used comparative genome hybridization to profile copy number polymorphisms. Repeated sequences and transposable elements (TEs), especially long terminal repeat retrotransposons, are densely methylated, but one third of genes also have low but detectable methylation in their transcribed regions. While TEs are almost always methylated, genic methylation is highly polymorphic, with half of all methylated genes being methylated in only one of the two ecotypes. A survey of loci in 96 Arabidopsis accessions revealed a similar degree of methylation polymorphism. Within-gene methylation is heritable, but is lost at a high frequency in segregating F₂ families. Promoter methylation is rare, and gene expression is not generally affected by differences in DNA methylation. Small interfering RNA are preferentially associated with methylated TEs, but not with methylated genes, indicating that most genic methylation is not guided by small interfering RNA. This may account for the instability of gene methylation, if occasional failure of maintenance methylation cannot be restored by other means.     

Folic acid supplementation alters the DNA methylation profile and improves insulin resistance in high-fat-diet-fed mice

Folic acid (FA) supplementation may protect from obesity and insulin resistance, the effects and mechanism of FA on chronic high-fat-diet-induced obesity-related metabolic disorders are not well elucidated. We adopted a genome-wide approach to directly examine whether FA supplementation affects the DNA methylation profile of mouse adipose tissue and identify the functional consequences of these changes. Mice were fed a high-fat diet (HFD), normal diet (ND) or an HFD supplemented with folic acid (20 μg/ml in drinking water) for 10 weeks, epididymal fat was harvested, and genome-wide DNA methylation analyses were performed using methylated DNA immunoprecipitation sequencing (MeDIP-seq). Mice exposed to the HFD expanded their adipose mass, which was accompanied by a significant increase in circulating glucose and insulin levels. FA supplementation reduced the fat mass and serum glucose levels and improved insulin resistance in HFD-fed mice. MeDIP-seq revealed distribution of differentially methylated regions (DMRs) throughout the adipocyte genome, with more hypermethylated regions in HFD mice. Methylome profiling identified DMRs associated with 3787 annotated genes from HFD mice in response to FA supplementation. Pathway analyses showed novel DNA methylation changes in adipose genes associated with insulin secretion, pancreatic secretion and type 2 diabetes. The differential DNA methylation corresponded to changes in the adipose tissue gene expression of Adcy3 and Rapgef4 in mice exposed to a diet containing FA. FA supplementation improved insulin resistance, decreased the fat mass, and induced DNA methylation and gene expression changes in genes associated with obesity and insulin secretion in obese mice fed a HFD.     

Hierarchical Clustering of Breast Cancer Methylomes Revealed Differentially Methylated and Expressed Breast Cancer Genes

Oncogenic transformation of normal cells often involves epigenetic alterations, including histone modification and DNA methylation. We conducted whole-genome bisulfite sequencing to determine the DNA methylomes of normal breast, fibroadenoma, invasive ductal carcinomas and MCF7. The emergence, disappearance, expansion and contraction of kilobase-sized hypomethylated regions (HMRs) and the hypomethylation of the megabase-sized partially methylated domains (PMDs) are the major forms of methylation changes observed in breast tumor samples. Hierarchical clustering of HMR revealed tumor-specific hypermethylated clusters and differential methylated enhancers specific to normal or breast cancer cell lines. Joint analysis of gene expression and DNA methylation data of normal breast and breast cancer cells identified differentially methylated and expressed genes associated with breast and/or ovarian cancers in cancer-specific HMR clusters. Furthermore, aberrant patterns of X-chromosome inactivation (XCI) was found in breast cancer cell lines as well as breast tumor samples in the TCGA BRCA (breast invasive carcinoma) dataset. They were characterized with differentially hypermethylated XIST promoter, reduced expression of XIST, and over-expression of hypomethylated X-linked genes. High expressions of these genes were significantly associated with lower survival rates in breast cancer patients. Comprehensive analysis of the normal and breast tumor methylomes suggests selective targeting of DNA methylation changes during breast cancer progression. The weak causal relationship between DNA methylation and gene expression observed in this study is evident of more complex role of DNA methylation in the regulation of gene expression in human epigenetics that deserves further investigation.     

Methylome sequencing for fibrolamellar hepatocellular carcinoma depicts distinctive features

With the goal of studying epigenetic alterations in fibrolamellar hepatocellular carcinoma (FLC) and establish an associated DNA methylation signature, we analyzed LINE-1 methylation in a cohort of FLC and performed next-generation sequencing of DNA methylation in a training set of pure-FLCs and non-cirrhotic hepatocellular carcinomas (nc-HCC). DNA methylation was correlated with gene expression. Furthermore, we established and validated an epigenetic signature differentiating pure-FLC from other HCCs. LINE-1 methylation correlated with shorter recurrence-free survival and overall survival in resected pure-FLC patients. Unsupervised clustering using CG sites located in islands distinguished pure-FLC from nc-HCC. Major DNA methylation changes occurred outside promoters, mainly in gene bodies and intergenic regions located in the vicinity of liver developmental genes (i.e., SMARCA4 and RXRA). Partially methylated domains were more prone to DNA methylation changes. Furthermore, we identified several putative tumor suppressor genes (e.g., DLEU7) and oncogenes (e.g., DUSP4). While ∼70% of identified gene promoters gaining methylation were marked by bivalent histone marks (H3K4me3/H3K27me3) in embryonic stem cells, ∼70% of those losing methylation were marked by H3K4me3. Finally, we established a pure FLC DNA methylation signature and validated it in an independent dataset. Our analysis reveals a distinct epigenetic signature of pure FLC as compared to nc-HCC, with DNA methylation changes occurring in the vicinity of liver developmental genes. These data suggest new options for targeting FLC based on cancer epigenome aberrations.     

Identification of New Differentially Methylated Genes That Have Potential Functional Consequences in Prostate Cancer

Many differentially methylated genes have been identified in prostate cancer (PCa), primarily using candidate gene-based assays. Recently, several global DNA methylation profiles have been reported in PCa, however, each of these has weaknesses in terms of ability to observe global DNA methylation alterations in PCa. We hypothesize that there remains unidentified aberrant DNA methylation in PCa, which may be identified using higher resolution assay methods. We used the newly developed Illumina HumanMethylation450 BeadChip in PCa (n = 19) and adjacent normal tissues (n = 4) and combined these with gene expression data for identifying new DNA methylation that may have functional consequences in PCa development and progression. We also confirmed our methylation results in an independent data set. Two aberrant DNA methylation genes were validated among an additional 56 PCa samples and 55 adjacent normal tissues. A total 28,735 CpG sites showed significant differences in DNA methylation (FDR adjusted P<0.05), defined as a mean methylation difference of at least 20% between PCa and normal samples. Furthermore, a total of 122 genes had more than one differentially methylated CpG site in their promoter region and a gene expression pattern that was inverse to the direction of change in DNA methylation (e.g. decreased expression with increased methylation, and vice-versa). Aberrant DNA methylation of two genes, AOX1 and SPON2, were confirmed via bisulfate sequencing, with most of the respective CpG sites showing significant differences between tumor samples and normal tissues. The AOX1 promoter region showed hypermethylation in 92.6% of 54 tested PCa samples in contrast to only three out of 53 tested normal tissues. This study used a new BeadChip combined with gene expression data in PCa to identify novel differentially methylated CpG sites located within genes. The newly identified differentially methylated genes may be used as biomarkers for PCa diagnosis.     

Genome-wide DNA methylation profiling of chronic lymphocytic leukemia allows identification of epigenetically repressed molecular pathways with clinical impact

We performed a genome-wide analysis of aberrant DNA methylation in chronic lymphocytic leukemia (CLL) using methylated CpG island amplification (MCA) coupled with a promoter microarray. We identified 280 potential targets of aberrant DNA methylation in CLL. These genes were located more frequently in chromosomes 19 (16%, p = 0.001), 16 (11%, p = 0.001), 17 (10%, p = 0.02) and 11 (9%, p = 0.02) and could be grouped in several functional networks. Methylation status was confirmed for 22 of these genes (SOX11, DLX1, FAM62C, SOX14, RSPO1, ADCY5, HAND2, SPOCK, MLL, ING1, PRIMA1, BCL11B, LTBP2, BNC1, NR2F2, SALL1, GALGT2, LHX1, DLX4, KLK10, TFAP2 and APP) in 78 CLL patients by pyrosequencing. As a proof of principle, we analyzed the expression of 2 genes, PRIMA1 and APP, in primary cells and of GALGT2, TFAP2C and PRIMA1 in leukemia cells. There was an inverse association between methylation and gene expression. This could be reversed by treatment with 5-aza-2′-deoxycytidine in cell lines. Treatment in a clinical trial with 5-azacitidine resulted in decreased methylation of LINE, DLX4 and SALL1 in the peripheral blood B-cells of patients with CLL. IgV_H mutational status or ZAP-70 expression were not associated with specific methylation profiles. By multivariate analysis, methylation of LINE and APP was associated with shorter overall survival (p = 0.045 and 0.0035, respectively). This study demonstrates that aberrant DNA methylation is common and has potential prognostic and therapeutic value in CLL.Key words: chronic lymphocytic leukemia, DNA methylation, MCA/promoter microarray, epigenetics     

Intronic Parent-of-Origin Dependent Differential Methylation at the Actn1 Gene Is Conserved in Rodents but Is Not Associated with Imprinted Expression

Parent-of-origin differential DNA methylation has been associated with regulation of the preferential expression of paternal or maternal alleles of imprinted genes. Based on this association, recent studies have searched for parent-of-origin dependent differentially methylated regions in order to identify new imprinted genes in their vicinity. In a previous genome-wide analysis of mouse brain DNA methylation, we found a novel differentially methylated region in a CpG island located in the last intron of the alpha 1 Actinin (Actn1) gene. In this region, preferential methylation of the maternal allele was observed; however, there were no reports of imprinted expression of Actn1. Therefore, we have tested if differential methylation of this region is common to other tissues and species and affects the expression of Actn1. We have found that Actn1 differential methylation occurs in diverse mouse tissues. Moreover, it is also present in other murine rodents (rat), but not in the orthologous human region. In contrast, we have found no indication of an imprinted effect on gene expression of Actn1 in mice: expression is always biallelic regardless of sex, tissue type, developmental stage or isoform. Therefore, we have identified a novel parent-of-origin dependent differentially methylated region that has no apparent association with imprinted expression of the closest genes. Our findings sound a cautionary note to genome-wide searches on the use of differentially methylated regions for the identification of imprinted genes and suggest that parent-of-origin dependent differential methylation might be conserved for functions other that the control of imprinted expression.     

Detection of MGMT, RASSF1A, p15INK4B, and p14ARF promoter methylation in circulating tumor-derived DNA of central nervous system cancer patients

Despite the growing understanding of the mechanisms of carcinogenesis, cancers of the central nervous system are usually associated with unfavorable prognosis. The use of an appropriate molecular marker may improve the treatment outcome by allowing early diagnosis and treatment susceptibility monitoring. Since methylation of tumor-derived DNA can be detected in the serum of cancer patients, this makes DNA methylation-based biomarkers one of the most promising diagnostic strategies. In this study, the methylation profiles of MGMT, RASSF1A, p15INK4B, and p14ARF genes were evaluated in serum free-circulating DNA and the corresponding tumor tissue in a group of 33 primary or metastatic central nervous system cancer patients. Gene promoter methylation was assessed using methylation-specific polymerase chain reaction (PCR). All the tested genes were found to be methylated to a different extent in both serum and tumor samples. In comparison to metastatic brain tumor patients, the patients with glial tumors were characterized by a higher frequency of gene hypermethylation. The hypermethylation of RASSF1A differentiated primary from metastatic brain cancers. Moreover, the gene methylation profiles observed in serum, in most cases, matched the methylation profiles detected in paired tumor samples.     

Lung Fibroblasts from Patients with Idiopathic Pulmonary Fibrosis Exhibit Genome-Wide Differences in DNA Methylation Compared to Fibroblasts from Nonfibrotic Lung

Excessive fibroproliferation is a central hallmark of idiopathic pulmonary fibrosis (IPF), a chronic, progressive disorder that results in impaired gas exchange and respiratory failure. Fibroblasts are the key effector cells in IPF, and aberrant expression of multiple genes contributes to their excessive fibroproliferative phenotype. DNA methylation changes are critical to the development of many diseases, but the DNA methylome of IPF fibroblasts has never been characterized. Here, we utilized the HumanMethylation 27 array, which assays the DNA methylation level of 27,568 CpG sites across the genome, to compare the DNA methylation patterns of IPF fibroblasts (n = 6) with those of nonfibrotic patient controls (n = 3) and commercially available normal lung fibroblast cell lines (n = 3). We found that multiple CpG sites across the genome are differentially methylated (as defined by P value less than 0.05 and fold change greater than 2) in IPF fibroblasts compared to fibroblasts from nonfibrotic controls. These methylation differences occurred both in genes recognized to be important in fibroproliferation and extracellular matrix generation, as well as in genes not previously recognized to participate in those processes (including organ morphogenesis and potassium ion channels). We used bisulfite sequencing to independently verify DNA methylation differences in 3 genes (CDKN2B, CARD10, and MGMT); these methylation changes corresponded with differences in gene expression at the mRNA and protein level. These differences in DNA methylation were stable throughout multiple cell passages. DNA methylation differences may thus help to explain a proportion of the differences in gene expression previously observed in studies of IPF fibroblasts. Moreover, significant variability in DNA methylation was observed among individual IPF cell lines, suggesting that differences in DNA methylation may contribute to fibroblast heterogeneity among patients with IPF. These results demonstrate that IPF fibroblasts exhibit global differences in DNA methylation that may contribute to the excessive fibroproliferation associated with this disease.     

Abnormal Hypermethylation at Imprinting Control Regions in Patients with S-Adenosylhomocysteine Hydrolase (AHCY) Deficiency

S-adenosylhomocysteine hydrolase (AHCY) deficiency is a rare autosomal recessive disorder in methionine metabolism caused by mutations in the AHCY gene. Main characteristics are psychomotor delay including delayed myelination and myopathy (hypotonia, absent tendon reflexes etc.) from birth, mostly associated with hypermethioninaemia, elevated serum creatine kinase levels and increased genome wide DNA methylation. The prime function of AHCY is to hydrolyse and efficiently remove S-adenosylhomocysteine, the by-product of transmethylation reactions and one of the most potent methyltransferase inhibitors. In this study, we set out to more specifically characterize DNA methylation changes in blood samples from patients with AHCY deficiency. Global DNA methylation was increased in two of three analysed patients. In addition, we analysed the DNA methylation levels at differentially methylated regions (DMRs) of six imprinted genes (MEST, SNRPN, LIT1, H19, GTL2 and PEG3) as well as Alu and LINE1 repetitive elements in seven patients. Three patients showed a hypermethylation in up to five imprinted gene DMRs. Abnormal methylation in Alu and LINE1 repetitive elements was not observed. We conclude that DNA hypermethylation seems to be a frequent but not a constant feature associated with AHCY deficiency that affects different genomic regions to different degrees. Thus AHCY deficiency may represent an ideal model disease for studying the molecular origins and biological consequences of DNA hypermethylation due to impaired cellular methylation status.     

Identification of differentially methylated markers among cytogenetic risk groups of acute myeloid leukemia

Aberrant DNA methylation is known to occur in cancer, including hematological malignancies such as acute myeloid leukemia (AML). However, less is known about whether specific methylation profiles characterize specific subcategories of AML. We examined this issue by using comprehensive high-throughput array-based relative methylation analysis (CHARM) to compare methylation profiles among patients in different AML cytogenetic risk groups. We found distinct profiles in each group, with the high-risk group showing overall increased methylation compared with low- and mid-risk groups. The differentially methylated regions (DMRs) distinguishing cytogenetic risk groups of AML were enriched in the CpG island shores. Specific risk-group associated DMRs were located near genes previously known to play a role in AML or other malignancies, such as MN1, UHRF1, HOXB3, and HOXB4, as well as TRIM71, the function of which in cancer is not well characterized. These findings were verified by quantitative bisulfite pyrosequencing and by comparison with results available at the TCGA cancer genome browser. To explore the potential biological significance of the observed methylation changes, we correlated our findings with gene expression data available through the TCGA database. The results showed that decreased methylation at HOXB3 and HOXB4 was associated with increased gene expression of both HOXB genes specific to the mid-risk AML, while increased DNA methylation at DCC distinctive to the high-risk AML was associated with increased gene expression. Our results suggest that the differential impact of cytogenetic changes on AML prognosis may, in part, be mediated by changes in methylation.     

Systematic profiling identifies PDLIM2 as a novel prognostic predictor for oesophageal squamous cell carcinoma (ESCC)

Till now, no appropriate biomarkers for high‐risk population screening and prognosis prediction have been identified for patients with oesophageal squamous cell carcinoma (ESCC). In this study, by the combined use of data from the Gene Expression Omnibus (GEO) datasets and The Cancer Genome Atlas (TCGA)‐oesophageal carcinoma (ESCA), we aimed to screen dysregulated genes with prognostic value in ESCC and the genetic and epigenetic alterations underlying the dysregulation. About 222 genes that had at least fourfold change in ESCC compared with adjacent normal tissues were identified using the microarray data in GDS3838. Among these genes, only PDLIM2 was associated with nodal invasion and overall survival (OS) at the same time. The high PDLIM2 expression group had significantly longer OS and its expression was independently associated with better OS (HR: 0.64, 95% CI: 0.43‐0.95, P = 0.03), after adjustment for gender and pathologic stages. The expression of its exon 7/8/9/10 had the highest AUC value (0.724) and better prognostic value (HR: 0.43, 95% CI: 0.22‐0.83, P = 0.01) than total PDLIM2 expression. PDLIM2 DNA copy deletion was common in ESCC and was associated with decreased gene expression. The methylation status of two CpG sites (cg23696886 and cg20449614) in the proximal promoter region of PDLIM2 showed a moderate negative correlation with the gene expression in PDLIM2 copy neutral/amplification group. In conclusion, we infer that PDLIM2 expression might be a novel prognostic indicator for ESCC patients. Its exon 7/8/9/10 expression had the best prognostic value. Its down‐regulation might be associated with gene‐level copy deletion and promoter hypermethylation.