Genome-wide DNA methylation studies suggest distinct DNA methylation patterns in pediatric embryonal and alveolar rhabdomyosarcomas

Rhabdomyosarcoma is the most common soft-tissue sarcoma in children. While cytogenetic abnormalities have been well characterized in this disease, aberrant epigenetic events such as DNA hypermethylation have not been described in genome-wide studies. We have analyzed the methylation status of 25,500 promoters in normal skeletal muscle, and in cell lines and tumor samples of embryonal and alveolar rhabdomyosarcoma from pediatric patients. We identified over 1,900 CpG islands that are hypermethylated in rhabdomyosarcomas relative to skeletal muscle. Genes involved in tissue development, differentiation, and oncogenesis such as DNAJA4, HES5, IRX1, BMP8A, GATA4, GATA6, ALX3, and P4HTM were hypermethylated in both RMS cell lines and primary samples, implicating aberrant DNA methylation in the pathogenesis of rhabdomyosarcoma. Furthermore, cluster analysis revealed embryonal and alveolar subtypes had distinct DNA methylation patterns, with the alveolar subtype being enriched in DNA hypermethylation of polycomb target genes. These results suggest that DNA methylation signatures may aid in the diagnosis and risk stratification of pediatric rhabdomyosarcoma and help identify new targets for therapy.     

DNA methylation changes associated with risk factors in tumors of the upper aerodigestive tract

Cancers of the upper aerodigestive tract (UADT) are common forms of malignancy associated with tobacco and alcohol exposures, although human papillomavirus and nutritional deficiency are also important risk factors. While somatically acquired DNA methylation changes have been associated with UADT cancers, what triggers these events and precise epigenetic targets are poorly understood. In this study, we applied quantitative profiling of DNA methylation states in a panel of cancer-associated genes to a case-control study of UADT cancers. Our analyses revealed a high frequency of aberrant hypermethylation of several genes, including MYOD1, CHRNA3 and MTHFR in UADT tumors, whereas CDKN2A was moderately hypermethylated. Among differentially methylated genes, we identified a new gene (the nicotinic acetycholine receptor gene) as target of aberrant hypermethylation in UADT cancers, suggesting that epigenetic deregulation of nicotinic acetycholine receptors in non-neuronal tissues may promote the development of UADT cancers. Importantly, we found that sex and age is strongly associated with the methylation states, whereas tobacco smoking and alcohol intake may also influence the methylation levels in specific genes. This study identifies aberrant DNA methylation patterns in UADT cancers and suggests a potential mechanism by which environmental factors may deregulate key cellular genes involved in tumor suppression and contribute to UADT cancers.     

DNA methylation abnormalities in congenital heart disease

Congenital heart defects represent the most common malformation at birth, occurring also in ∼50% of individuals with Down syndrome. Congenital heart defects are thought to have multifactorial etiology, but the main causes are largely unknown. We have explored the global methylation profile of fetal heart DNA in comparison to blood DNA from control subjects: an absolute correlation with the type of tissue was detected. Pathway analysis revealed a significant enrichment of differential methylation at genes related to muscle contraction and cardiomyopathies in the developing heart DNA. We have also searched for abnormal methylation profiles on developing heart-tissue DNA of syndromic and non-syndromic congenital heart defects. On average, 3 regions with aberrant methylation were detected per sample and 18 regions were found differentially methylated between groups. Several epimutations were detected in candidate genes involved in growth regulation, apoptosis and folate pathway. A likely pathogenic hypermethylation of several intragenic sites at the MSX1 gene, involved in outflow tract morphogenesis, was found in a fetus with isolated heart malformation. In addition, hypermethylation of the GATA4 gene was present in fetuses with Down syndrome with or without congenital heart defects, as well as in fetuses with isolated heart malformations. Expression deregulation of the abnormally methylated genes was detected. Our data indicate that epigenetic alterations of relevant genes are present in developing heart DNA in fetuses with both isolated and syndromic heart malformations. These epimutations likely contribute to the pathogenesis of the malformation by cis-acting effects on gene expression.     

DNA methylation analyses of the connexin gene family reveal silencing of GJC1 (Connexin45) by promoter hypermethylation in colorectal cancer

Gap junctions are specialized plasma membrane domains consisting of channels formed by members of the connexin protein family. Gap junctional intercellular communication is often lost in cancers due to aberrant localization or downregulation of connexins, and connexins are therefore suggested to act as tumor suppressor genes in various tissues. The aim of this study was to investigate the expression pattern and DNA promoter methylation status of connexins in colorectal cancer. Expression of six (GJA1, GJA9, GJB1, GJB2, GJC1 and GJD3) connexin genes was detected in normal colonic tissue samples. GJC1 expression was reduced in colorectal carcinomas compared to normal tissue samples. All analyzed connexins were hypermethylated in colon cancer cell lines, although at various frequencies. GJA4, GJB6 and GJD2 were hypermethylated in 60% (29/48), 25% (12/48) and 96% (46/48) of primary colorectal carcinomas, respectively. However, the methylation status was not associated with gene expression. GJC1 has two alternative promoters, which were methylated in 42% (32/76) and 38% (25/65) of colorectal tumors, and in none of the normal mucosa samples. Expression of GJC1 was significantly lower in methylated compared with unmethylated samples (p < 0.01) and was restored in cell lines treated with the demethylating drug 5-aza-2'deoxycytidine. Taken together, DNA hypermethylation of the promoter region of GJC1, encoding connexin45, is an important mechanism in silencing gene expression in colorectal cancer.     

Genome‐wide mapping of cytosine methylation revealed dynamic DNA methylation patterns associated with genes and centromeres in rice

We conducted genome‐wide mapping of cytosine methylation using methylcytosine immunoprecipitation combined with Illumina sequencing. The chromosomal distribution pattern of methylated DNA is similar to the heterochromatin distribution pattern on rice chromosomes. The DNA methylation patterns of rice genes are similar to those in Arabidopsis thaliana, including distinct methylation patterns asssociated with gene bodies and promoters. The DNA sequences in the core domains of rice Cen4, Cen5 and Cen8 showed elevated methylation levels compared with sequences in the pericentromeric regions. In addition, elevated methylation levels were associated with the DNA sequences in the CENH3‐binding subdomains, compared with the sequences in the flanking H3 subdomains. In contrast, the centromeric domain of Cen11, which is composed exclusively of centromeric satellite DNA, is hypomethylated compared with the pericentromeric domains. Thus, the DNA sequences associated with functional centromeres can be either hypomethylated or hypermethylated. The methylation patterns of centromeric DNA appear to be correlated with the composition of the associated DNA sequences. We propose that both hypomethylation and hypermethylation of CENH3‐associated DNA sequences can serve as epigenetic marks to distinguish where CENH3 deposition will occur within the surrounding H3 chromatin.     

Global and gene-specific promoter methylation analysis in primary hyperparathyroidism

Epigenetic mechanisms involved in primary hyperparathyroidism are poorly understood as studies are limited. In order to understand the role of aberrant DNA promoter methylation in the pathogenesis of parathyroid tumors, we have quantified the CpG island promoter methylation density of several candidate genes including APC (promoter 1A and 1B), β-catenin (CTNNB1), CASR, CDC73/HRPT2, MEN1, P16 (CDKN2A), PAX1, RASSF1A, SFRP1 and VDR in 72 parathyroid tumors and 3 normal parathyroid references using bisulfite pyrosequencing. Global methylation levels were assessed for LINE-1. We also compared methylation levels with gene expression levels measured by qRT-PCR for genes showing frequent hypermethylation. The adenomas displayed frequent hypermethylation of APC 1A (37/66; 56%), RASSF1A (34/66; 52%) and β-catenin (19/66; 29%). One of the three atypical adenomas was hypermethylated for APC 1A. The three carcinomas were hypermethylated for RASSF1A and SFRP1, and the latter was only observed in this subtype. The global methylation density was similar in tumors (mean 70%) and parathyroid reference samples (mean 70%). In general, hypermethylated genes had reduced expression in the parathyroid adenomas using qRT-PCR. Among the adenomas, methylation of APC 1A correlated with adenoma weight (r = 0.306, p < 0.05). Furthermore, the methylation status of RASSF1A correlated with each of APC 1A (r = 0.289, p < 0.05) and β-catenin (r = 0.315, p < 0.01). Our findings suggest a role for aberrant DNA promoter methylation of APC 1A, β-catenin and RASSF1A in a subset of parathyroid tumors.     

DNA methylation changes associated with risk factors in tumors of the upper aerodigestive tract

Cancers of the upper aerodigestive tract (UADT) are common forms of malignancy associated with tobacco and alcohol exposures, although human papillomavirus and nutritional deficiency are also important risk factors. While somatically acquired DNA methylation changes have been associated with UADT cancers, what triggers these events and precise epigenetic targets are poorly understood. In this study, we applied quantitative profiling of DNA methylation states in a panel of cancer-associated genes to a case-control study of UADT cancers. Our analyses revealed a high frequency of aberrant hypermethylation of several genes, including MYOD1, CHRNA3 and MTHFR in UADT tumors, whereas CDKN2A was moderately hypermethylated. Among differentially methylated genes, we identified a new gene (the nicotinic acetycholine receptor gene) as target of aberrant hypermethylation in UADT cancers, suggesting that epigenetic deregulation of nicotinic acetycholine receptors in non-neuronal tissues may promote the development of UADT cancers. Importantly, we found that sex and age is strongly associated with the methylation states, whereas tobacco smoking and alcohol intake may also influence the methylation levels in specific genes. This study identifies aberrant DNA methylation patterns in UADT cancers and suggests a potential mechanism by which environmental factors may deregulate key cellular genes involved in tumor suppression and contribute to UADT cancers.Key words: DNA methylation, upper aerodigestive tract, cancer, risk factors, biomarkers     

Fetal DNA hypermethylation in tight junction pathway is associated with neural tube defects: A genome-wide DNA methylation analysis

Neural tube defects (NTDs) are a spectrum of severe congenital malformations of fusion failure of the neural tube during early embryogenesis. Evidence on aberrant DNA methylation in NTD development remains scarce, especially when exposure to environmental pollutant is taken into consideration. DNA methylation profiling was quantified using the Infinium HumanMethylation450 array in neural tissues from 10 NTD cases and 8 non-malformed controls (stage 1). Subsequent validation was performed using a Sequenom MassARRAY system in neural tissues from 20 NTD cases and 20 non-malformed controls (stage 2). Correlation analysis of differentially methylated CpG sites in fetal neural tissues and polycyclic aromatic hydrocarbons concentrations in fetal neural tissues and maternal serum was conducted. Differentially methylated CpG sites of neural tissues were further validated in fetal mice with NTDs induced by benzo(a)pyrene given to pregnant mice. Differentially hypermethylated CpG sites in neural tissues from 17 genes and 6 pathways were identified in stage 1. Subsequently, differentially hypermethylated CpG sites in neural tissues from 6 genes (BDKRB2, CTNNA1, CYFIP2, MMP7, MYH2, and TIAM2) were confirmed in stage 2. Correlation analysis showed that methylated CpG sites in CTNNA1 and MYH2 from NTD cases were positively correlated to polycyclic aromatic hydrocarbon level in fetal neural tissues and maternal serum. The correlation was confirmed in NTD-affected fetal mice that were exposed to benzo(a)pyrene in utero. In conclusion, hypermethylation of the CTNNA1 and MYH2 genes in tight junction pathway is associated with the risk for NTDs, and the DNA methylation aberration may be caused by exposure to benzo(a)pyrene.     

DNA methylation changes in atypical adenomatous hyperplasia, adenocarcinoma in situ, and lung adenocarcinoma

Background

Aberrant DNA methylation is common in lung adenocarcinoma, but its timing in the phases of tumor development is largely unknown. Delineating when abnormal DNA methylation arises may provide insight into the natural history of lung adenocarcinoma and the role that DNA methylation alterations play in tumor formation.

Methodology/Principal Findings

We used MethyLight, a sensitive real-time PCR-based quantitative method, to analyze DNA methylation levels at 15 CpG islands that are frequently methylated in lung adenocarcinoma and that we had flagged as potential markers for non-invasive detection. We also used two repeat probes as indicators of global DNA hypomethylation. We examined DNA methylation in 249 tissue samples from 93 subjects, spanning the putative spectrum of peripheral lung adenocarcinoma development: histologically normal adjacent non-tumor lung, atypical adenomatous hyperplasia (AAH), adenocarcinoma in situ (AIS, formerly known as bronchioloalveolar carcinoma), and invasive lung adenocarcinoma. Comparison of DNA methylation levels between the lesion types suggests that DNA hypermethylation of distinct loci occurs at different time points during the development of lung adenocarcinoma. DNA methylation at CDKN2A ex2 and PTPRN2 is already significantly elevated in AAH, while CpG islands at 2C35, EYA4, HOXA1, HOXA11, NEUROD1, NEUROD2 and TMEFF2 are significantly hypermethylated in AIS. In contrast, hypermethylation at CDH13, CDX2, OPCML, RASSF1, SFRP1 and TWIST1 and global DNA hypomethylation appear to be present predominantly in invasive cancer.

Conclusions/Significance

The gradual increase in DNA methylation seen for numerous loci in progressively more transformed lesions supports the model in which AAH and AIS are sequential stages in the development of lung adenocarcinoma. The demarcation of DNA methylation changes characteristic for AAH, AIS and adenocarcinoma begins to lay out a possible roadmap for aberrant DNA methylation events in tumor development. In addition, it identifies which DNA methylation changes might be used as molecular markers for the detection of preinvasive lesions.     

Genome-wide aberrant DNA methylation of microRNA host genes in hepatocellular carcinoma

Mature microRNAs (miRNAs) are a class of small non-coding RNAs involved in posttranslational gene silencing. Previous studies found that downregulation of miRNAs is a common feature observed in solid tumors, including hepatocellular carcinoma (HCC). We employed a genome-wide approach to test the hypothesis that DNA methylation alterations in miRNA host genes may cause deregulated miRNA expression in HCC. We analyzed tumor and adjacent non-tumor tissues from 62 Taiwanese HCC cases using Infinium HumanMethylation27 DNA Analysis BeadChips that include 254 CpG sites covering 110 miRNAs from 64 host genes. Expression levels of three identified miRNAs (miR-10a, miR-10b and miR-196b) were measured in a subset of 37 HCC tumor and non-tumor tissues. After Bonferroni adjustment, a total of 54 CpG sites from 27 host genes significantly differed in DNA methylation levels between tumor and adjacent non-tumor tissues with 53 sites significantly hypermethylated in tumor tissues. Among the 54 significant CpG sites, 15 sites had more than 2-fold tumor/non-tumor changes, 17 sites had differences > 10%, and 10 sites had both features [including 8 significantly hypermethylated CpG sites in the host genes of miR-10a, miR-10b and miR-196b (HOXB4, HOXD4 and HOXA9, respectively)]. Significant downregulation of miR-10a was observed in tumor compared with non-tumor tissues (0.50 vs. 1.73, p = 0.031). The concordance for HOXB4 methylation alteration and dysregulation of miR-10a was 73.5%. No significant change was observed for miR-10b expression. Unexpectedly, miR-196b was significantly upregulated in tumor compared with non-tumor tissues (p = 0.0001). These data suggest that aberrant DNA methylation may lead to dysregulation of miR-10a in HCC tumor tissues.     

Comprehensive methylome analysis of ovarian tumors reveals hedgehog signaling pathway regulators as prognostic DNA methylation biomarkers

Women with advanced stage ovarian cancer (OC) have a five-year survival rate of less than 25%. OC progression is associated with accumulation of epigenetic alterations and aberrant DNA methylation in gene promoters acts as an inactivating “hit” during OC initiation and progression. Abnormal DNA methylation in OC has been used to predict disease outcome and therapy response. To globally examine DNA methylation in OC, we used next-generation sequencing technology, MethylCap-sequencing, to screen 75 malignant and 26 normal or benign ovarian tissues. Differential DNA methylation regions (DMRs) were identified, and the Kaplan–Meier method and Cox proportional hazard model were used to correlate methylation with clinical endpoints. Functional role of specific genes identified by MethylCap-sequencing was examined in in vitro assays. We identified 577 DMRs that distinguished (p < 0.001) malignant from non-malignant ovarian tissues; of these, 63 DMRs correlated (p < 0.001) with poor progression free survival (PFS). Concordant hypermethylation and corresponding gene silencing of sonic hedgehog pathway members ZIC1 and ZIC4 in OC tumors was confirmed in a panel of OC cell lines, and ZIC1 and ZIC4 repression correlated with increased proliferation, migration and invasion. ZIC1 promoter hypermethylation correlated (p < 0.01) with poor PFS. In summary, we identified functional DNA methylation biomarkers significantly associated with clinical outcome in OC and suggest our comprehensive methylome analysis has significant translational potential for guiding the design of future clinical investigations targeting the OC epigenome. Methylation of ZIC1, a putative tumor suppressor, may be a novel determinant of OC outcome.     

Genome‐wide DNA methylation changes with age in disease‐free human skeletal muscle

A decline in skeletal muscle mass and function with aging is well recognized, but remains poorly characterized at the molecular level. Here, we report for the first time a genome‐wide study of DNA methylation dynamics in skeletal muscle of healthy male individuals during normal human aging. We predominantly observed hypermethylation throughout the genome within the aged group as compared to the young subjects. Differentially methylated CpG (dmCpG) nucleotides tend to arise intragenically and are underrepresented in promoters and are overrepresented in the middle and 3′ end of genes. The intragenic methylation changes are overrepresented in genes that guide the formation of the junction of the motor neuron and myofibers. We report a low level of correlation of gene expression from previous studies of aged muscle with our current analysis of DNA methylation status. For those genes that had both changes in methylation and gene expression with age, we observed a reverse correlation, with the exception of intragenic hypermethylated genes that were correlated with an increased gene expression. We suggest that a minimal number of dmCpG sites or select sites are required to be altered in order to correlate with gene expression changes. Finally, we identified 500 dmCpG sites that perform well in discriminating young from old samples. Our findings highlight epigenetic links between aging postmitotic skeletal muscle and DNA methylation.     

Genome-wide DNA methylation profiling of stomach cancer in the ethnic population of Mizoram,North East India

Stomach cancer is the fifth most common cancer in terms of prevalence and incidence and the fourth leading cause of mortality in men and women worldwide. It is well-established that aberrant DNA methylation in cells can lead to carcinogenesis. The primary objective of our study was to investigate the aberrant DNA methylation status of genes associated with stomach cancer with a particular reference to the ethnic population of Mizoram, North East India. The site-level analysis identified 2883 CpG sites differentially methylated, representing ～922 genes. Out of which 476 Differentially Methylated Positions (DMPs) were promoter-associated, 452 DMPs were hypermethylated, and 24 were hypomethylated. The region-level analysis identified 462 Differentially Methylated Regions (DMRs) corresponding to ～320 genes, of which ～281 genes were hypermethylated and ～40 genes were hypomethylated. TCGA analysis showed that some of the genes had been previously implicated in other cancers including stomach cancer. Five hypermethylated genes were selected as candidate genes for further investigations and they have shown to be novel and could serve as candidate hypermethylation biomarkers for stomach cancer in this particular ethnic group.     

Genome-wide signatures of differential DNA methylation in pediatric acute lymphoblastic leukemia

Background

Although aberrant DNA methylation has been observed previously in acute lymphoblastic leukemia (ALL), the patterns of differential methylation have not been comprehensively determined in all subtypes of ALL on a genome-wide scale. The relationship between DNA methylation, cytogenetic background, drug resistance and relapse in ALL is poorly understood.

Results

We surveyed the DNA methylation levels of 435,941 CpG sites in samples from 764 children at diagnosis of ALL and from 27 children at relapse. This survey uncovered four characteristic methylation signatures. First, compared with control blood cells, the methylomes of ALL cells shared 9,406 predominantly hypermethylated CpG sites, independent of cytogenetic background. Second, each cytogenetic subtype of ALL displayed a unique set of hyper- and hypomethylated CpG sites. The CpG sites that constituted these two signatures differed in their functional genomic enrichment to regions with marks of active or repressed chromatin. Third, we identified subtype-specific differential methylation in promoter and enhancer regions that were strongly correlated with gene expression. Fourth, a set of 6,612 CpG sites was predominantly hypermethylated in ALL cells at relapse, compared with matched samples at diagnosis. Analysis of relapse-free survival identified CpG sites with subtype-specific differential methylation that divided the patients into different risk groups, depending on their methylation status.

Conclusions

Our results suggest an important biological role for DNA methylation in the differences between ALL subtypes and in their clinical outcome after treatment.