Exploring genome-wide DNA methylation profiles altered in hepatocellular carcinoma using Infinium HumanMethylation 450 BeadChips

Hepatocellular carcinoma (HCC) incidence has increased in the US and also has one of the fastest growing death rates of any cancer. The purpose of the current study was to discover novel genome-wide aberrant DNA methylation patterns in HCC tumors that are predominantly HCV-related. Infinium HumanMethylation 450K BeadChip arrays were used to examine genome-wide DNA methylation profiles in 66 pairs of HCC tumor and adjacent non-tumor tissues. After Bonferroni adjustment, a total of 130,512 CpG sites significantly differed in methylation level in tumor compared with non-tumor tissues, with 28,017 CpG sites hypermethylated and 102,495 hypomethylated in tumor tissues. Absolute tumor/non-tumor methylation differences ≥ 20% were found in 24.9% of the hypermethylated and 43.1% of the hypomethylated CpG sites; almost 10,000 CpG sites have ≥ 30% DNA methylation differences. Most (60.1%) significantly hypermethylated CpG sites are located in CpG islands, with 21.6% in CpG shores and 3.6% in shelves. In contrast, only a small proportion (8.2%) of significantly hypomethylated CpG sites are situated in islands, while most are found in open sea (60.2%), shore (17.3%) or shelf (14.3%) regions. A total of 2,568 significant CpG sites (2,441 hypermethylated and 127 hypomethylated) covering 589 genes are located within 684 differentially methylated regions defined as regions with at least two significant CpG sites displaying > 20% methylation differences in the same direction within 250-bp. The top 500 significant CpG sites can significantly distinguish HCC tumor from adjacent tissues with one misclassification. Within adjacent non-tumor tissues, we also identified 75 CpG sites significantly associated with gender, 228 with HCV infection, 17,207 with cirrhosis, and 56 with both HCV infection and cirrhosis after multiple comparisons adjustment. Aberrant DNA methylation profiles across the genome were identified in tumor tissues from US HCC cases that are predominantly related to HCV infection. These results demonstrate the significance of aberrant DNA methylation in HCC tumorigenesis.     

DNA Methylation Profile Distinguishes Clear Cell Sarcoma of the Kidney from Other Pediatric Renal Tumors

A number of specific, distinct neoplastic entities occur in the pediatric kidney, including Wilms’ tumor, clear cell sarcoma of the kidney (CCSK), congenital mesoblastic nephroma (CMN), rhabdoid tumor of the kidney (RTK), and the Ewing’s sarcoma family of tumors (ESFT). By employing DNA methylation profiling using Illumina Infinium HumanMethylation27, we analyzed the epigenetic characteristics of the sarcomas including CCSK, RTK, and ESFT in comparison with those of the non-neoplastic kidney (NK), and these tumors exhibited distinct DNA methylation profiles in a tumor-type-specific manner. CCSK is the most frequently hypermethylated, but least frequently hypomethylated, at CpG sites among these sarcomas, and exhibited 490 hypermethylated and 46 hypomethylated CpG sites in compared with NK. We further validated the results by MassARRAY, and revealed that a combination of four genes was sufficient for the DNA methylation profile-based differentiation of these tumors by clustering analysis. Furthermore, THBS1 CpG sites were found to be specifically hypermethylated in CCSK and, thus, the DNA methylation status of these THBS1 sites alone was sufficient for the distinction of CCSK from other pediatric renal tumors, including Wilms’ tumor and CMN. Moreover, combined bisulfite restriction analysis could be applied for the detection of hypermethylation of a THBS1 CpG site. Besides the biological significance in the pathogenesis, the DNA methylation profile should be useful for the differential diagnosis of pediatric renal tumors.     

A Comprehensive Microarray-Based DNA Methylation Study of 367 Hematological Neoplasms

Background

Alterations in the DNA methylation pattern are a hallmark of leukemias and lymphomas. However, most epigenetic studies in hematologic neoplasms (HNs) have focused either on the analysis of few candidate genes or many genes and few HN entities, and comprehensive studies are required.

Methodology/Principal Findings

Here, we report for the first time a microarray-based DNA methylation study of 767 genes in 367 HNs diagnosed with 16 of the most representative B-cell (n = 203), T-cell (n = 30), and myeloid (n = 134) neoplasias, as well as 37 samples from different cell types of the hematopoietic system. Using appropriate controls of B-, T-, or myeloid cellular origin, we identified a total of 220 genes hypermethylated in at least one HN entity. In general, promoter hypermethylation was more frequent in lymphoid malignancies than in myeloid malignancies, being germinal center mature B-cell lymphomas as well as B and T precursor lymphoid neoplasias those entities with highest frequency of gene-associated DNA hypermethylation. We also observed a significant correlation between the number of hypermethylated and hypomethylated genes in several mature B-cell neoplasias, but not in precursor B- and T-cell leukemias. Most of the genes becoming hypermethylated contained promoters with high CpG content, and a significant fraction of them are targets of the polycomb repressor complex. Interestingly, T-cell prolymphocytic leukemias show low levels of DNA hypermethylation and a comparatively large number of hypomethylated genes, many of them showing an increased gene expression.

Conclusions/Significance

We have characterized the DNA methylation profile of a wide range of different HNs entities. As well as identifying genes showing aberrant DNA methylation in certain HN subtypes, we also detected six genes—DBC1, DIO3, FZD9, HS3ST2, MOS, and MYOD1—that were significantly hypermethylated in B-cell, T-cell, and myeloid malignancies. These might therefore play an important role in the development of different HNs.     

Characterization of the IGF2 Imprinted Gene Methylation Status in Bovine Oocytes during Folliculogenesis

DNA methylation reprogramming occurs during mammalian gametogenesis and embryogenesis. Sex-specific DNA methylation patterns at specific CpG islands controlling imprinted genes are acquired during this window of development. Characterization of the DNA methylation dynamics of imprinted genes acquired by oocytes during folliculogenesis is essential for understanding the physiological and genetic aspects of female gametogenesis and to determine the parameters for oocyte competence. This knowledge can be used to improve in vitro embryo production (IVP), specifically because oocyte competence is one of the most important aspects determining the success of IVP. Imprinted genes, such as IGF2, play important roles in embryo development, placentation and fetal growth. The aim of this study was to characterize the DNA methylation profile of the CpG island located in IGF2 exon 10 in oocytes during bovine folliculogenesis. The methylation percentages in oocytes from primordial follicles, final secondary follicles, small antral follicles, large antral follicles, MII oocytes and spermatozoa were 73.74 ± 2.88%, 58.70 ± 7.46%, 56.00 ± 5.58%, 65.77 ± 5.10%, 56.35 ± 7.45% and 96.04 ± 0.78%, respectively. Oocytes from primordial follicles showed fewer hypomethylated alleles (15.5%) than MII oocytes (34.6%) (p = 0.039); spermatozoa showed only hypermethylated alleles. Moreover, MII oocytes were less methylated than spermatozoa (p<0.001). Our results showed that the methylation pattern of this region behaves differently between mature oocytes and spermatozoa. However, while this region has a classical imprinted pattern in spermatozoa that is fully methylated, it was variable in mature oocytes, showing hypermethylated and hypomethylated alleles. Furthermore, our results suggest that this CpG island may have received precocious reprogramming, considering that the hypermethylated pattern was already found in growing oocytes from primordial follicles. These results may contribute to our understanding of the reprogramming of imprinted genes during bovine oogenesis.     

Epigenetic Basis of Regeneration: Analysis of Genomic DNA Methylation Profiles in the MRL/MpJ Mouse

Epigenetic regulation plays essential role in cell differentiation and dedifferentiation, which are the intrinsic processes involved in regeneration. To investigate the epigenetic basis of regeneration capacity, we choose DNA methylation as one of the most important epigenetic mechanisms and the MRL/MpJ mouse as a model of mammalian regeneration known to exhibit enhanced regeneration response in different organs. We report the comparative analysis of genomic DNA methylation profiles of the MRL/MpJ and the control C57BL/6J mouse. Methylated DNA immunoprecipitation followed by microarray analysis using the Nimblegen ‘3 × 720 K CpG Island Plus RefSeq Promoter’ platform was applied in order to carry out genome-wide DNA methylation profiling covering 20 404 promoter regions. We identified hundreds of hypo- and hypermethylated genes and CpG islands in the heart, liver, and spleen, and 37 of them in the three tissues. Decreased inter-tissue diversification and the shift of DNA methylation balance upstream the genes distinguish the genomic methylation patterns of the MRL/MpJ mouse from the C57BL/6J. Homeobox genes and a number of other genes involved in embryonic morphogenesis are significantly overrepresented among the genes hypomethylated in the MRL/MpJ mouse. These findings indicate that epigenetic patterning might be a likely molecular basis of regeneration capability in the MRL/MpJ mouse.     

High-density array analysis of DNA methylation in Tamoxifen-resistant breast cancer cell lines

Roughly two-thirds of all breast cancers are ERα-positive and can be treated with the antiestrogen, Tamoxifen, however resistance occurs in 33% of women who take the drug for more than 5 y. Aberrant DNA methylation, an epigenetic mechanism that alters gene expression in cancer, is thought to play a role in this resistance. To develop an understanding of Tamoxifen-resistance and identify novel pathways and targets of aberrant methylation, DNA from MCF-7 breast cancer cells and Tamoxifen-resistant derivatives, TMX2–11 and TMX2–28, were analyzed using the Illumina HumanMethylation450 BeadChip. Normalizing against MCF-7 values, ERα-positive TMX2–11 had 4000 hypermethylated sites and ERα-negative TMX2–28 had over 33 000. Analysis of CpG sites altered in both TMX2–11 and TMX2–28 revealed that the Tamoxifen-resistant cell lines share 3000 hypermethylated and 200 hypomethylated CpGs. ZNF350 and MAGED1, two genes hypermethylated in both cell lines, were examined in greater detail. Treatment with 5-aza-2′deoxycitidine caused a significant reduction in promoter methylation of both ZNF350 and MAGED1 and a corresponding increase in expression in TMX2–28. A similar relationship between methylation and expression was not detected in TMX2–11. Our findings are indicative of the variable responses to methylation-targeted breast cancer therapy and highlight the need for biomarkers that accurately predict treatment outcome.     

Short‐term calorie restriction ameliorates genomewide,age‐related alterations in DNA methylation

DNA methylation plays major roles in many biological processes, including aging, carcinogenesis, and development. Analyses of DNA methylation using next‐generation sequencing offer a new way to profile and compare methylomes across the genome in the context of aging. We explored genomewide DNA methylation and the effects of short‐term calorie restriction (CR) on the methylome of aged rat kidney. Whole‐genome methylation of kidney in young (6 months old), old (25 months old), and OCR (old with 4‐week, short‐term CR) rats was analyzed by methylated DNA immunoprecipitation and next‐generation sequencing (MeDIP‐Seq). CpG islands and repetitive regions were hypomethylated, but 5′‐UTR, exon, and 3′‐UTR hypermethylated in old and OCR rats. The methylation in the promoter and intron regions was decreased in old rats, but increased in OCR rats. Pathway enrichment analysis showed that the hypermethylated promoters in old rats were associated with degenerative phenotypes such as cancer and diabetes. The hypomethylated promoters in old rats related significantly to the chemokine signaling pathway. However, the pathways significantly enriched in old rats were not observed from the differentially methylated promoters in OCR rats. Thus, these findings suggest that short‐term CR could partially ameliorate age‐related methylation changes in promoters in old rats. From the epigenomic data, we propose that the hypermethylation found in the promoter regions of disease‐related genes during aging may indicate increases in susceptibility to age‐related diseases. Therefore, the CR‐induced epigenetic changes that ameliorate age‐dependent aberrant methylation may be important to CR's health‐ and life‐prolonging effects.     

DNA methylome profiling identifies novel methylated genes in African American patients with colorectal neoplasia

The identification of genes that are differentially methylated in colorectal cancer (CRC) has potential value for both diagnostic and therapeutic interventions specifically in high-risk populations such as African Americans (AAs). However, DNA methylation patterns in CRC, especially in AAs, have not been systematically explored and remain poorly understood. Here, we performed DNA methylome profiling to identify the methylation status of CpG islands within candidate genes involved in critical pathways important in the initiation and development of CRC. We used reduced representation bisulfite sequencing (RRBS) in colorectal cancer and adenoma tissues that were compared with DNA methylome from a healthy AA subject’s colon tissue and peripheral blood DNA. The identified methylation markers were validated in fresh frozen CRC tissues and corresponding normal tissues from AA patients diagnosed with CRC at Howard University Hospital. We identified and validated the methylation status of 355 CpG sites located within 16 gene promoter regions associated with CpG islands. Fifty CpG sites located within CpG islands—in genes ATXN7L1 (2), BMP3 (7), EID3 (15), GAS7 (1), GPR75 (24), and TNFAIP2 (1)—were significantly hypermethylated in tumor vs. normal tissues (P < 0.05). The methylation status of BMP3, EID3, GAS7, and GPR75 was confirmed in an independent, validation cohort. Ingenuity pathway analysis mapped three of these markers (GAS7, BMP3 and GPR) in the insulin and TGF-β1 network—the two key pathways in CRC. In addition to hypermethylated genes, our analysis also revealed that LINE-1 repeat elements were progressively hypomethylated in the normal-adenoma-cancer sequence. We conclude that DNA methylome profiling based on RRBS is an effective method for screening aberrantly methylated genes in CRC. While previous studies focused on the limited identification of hypermethylated genes, ours is the first study to systematically and comprehensively identify novel hypermethylated genes, as well as hypomethylated LINE-1 sequences, which may serve as potential biomarkers for CRC in African Americans. Our discovered biomarkers were intimately linked to the insulin/TGF-B1 pathway, further strengthening the association of diabetic disorders with colon oncogenic transformation.     

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 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.     

Fetal DNA Methylation Associates with Early Spontaneous Preterm Birth and Gestational Age

Spontaneous preterm birth (PTB, <37 weeks gestation) is a major public health concern, and children born preterm have a higher risk of morbidity and mortality throughout their lives. Recent studies suggest that fetal DNA methylation of several genes varies across a range of gestational ages (GA), but it is not yet clear if fetal epigenetic changes associate with PTB. The objective of this study is to interrogate methylation patterns across the genome in fetal leukocyte DNA from African Americans with early PTB (24^1/7–34^0/7 weeks; N = 22) or term births (39^0/7–40^6/7weeks; N = 28) and to evaluate the association of each CpG site with PTB and GA. DNA methylation was assessed across the genome with the HumanMethylation450 BeadChip. For each individual sample and CpG site, the proportion of DNA methylation was estimated. The associations between methylation and PTB or GA were evaluated by fitting a separate linear model for each CpG site, adjusting for relevant covariates. Overall, 29 CpG sites associated with PTB (FDR<.05; 5.7×10⁻¹⁰<p<2.9×10⁻⁶) independent of GA. Also, 9637 sites associated with GA (FDR<.05; 9.5×10⁻¹⁶<p<1.0×10⁻³), with 61.8% decreasing in methylation with shorter GA. GA-associated CpG sites were depleted in the CpG islands of their respective genes (p<2.2×10⁻¹⁶). Gene set enrichment analysis (GSEA) supported enrichment of GA-associated CpG sites in genes that play a role in embryonic development as well as the extracellular matrix. Additionally, this study replicated the association of several CpG sites associated with gestational age in other studies (CRHBP, PIK3CD and AVP). Dramatic differences in fetal DNA methylation are evident in fetuses born preterm versus at term, and the patterns established at birth may provide insight into the long-term consequences associated with PTB.     

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.