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.     

Sperm methylation profiles reveal features of epigenetic inheritance and evolution in primates

During germ cell and preimplantation development, mammalian cells undergo nearly complete reprogramming of DNA methylation patterns. We profiled the methylomes of human and chimp sperm as a basis for comparison to methylation patterns of ESCs. Although the majority of promoters escape methylation in both ESCs and sperm, the corresponding hypomethylated regions show substantial structural differences. Repeat elements are heavily methylated in both germ and somatic cells; however, retrotransposons from several subfamilies evade methylation more effectively during male germ cell development, whereas other subfamilies show the opposite trend. Comparing methylomes of human and chimp sperm revealed a subset of differentially methylated promoters and strikingly divergent methylation in retrotransposon subfamilies, with an evolutionary impact that is apparent in the underlying genomic sequence. Thus, the features that determine DNA methylation patterns differ between male germ cells and somatic cells, and elements of these features have diverged between humans and chimpanzees.     

Normal early pregnancy

The objective of this study was to analyze genome-wide differential methylation patterns in maternal leukocyte DNA in early pregnant and non-pregnant states. This is an age and body mass index matched case-control study comparing the methylation patterns of 27,578 cytosine-guanine (CpG) sites in 14,495 genes in maternal leukocyte DNA in early pregnancy (n = 14), in the same women postpartum (n = 14), and in nulligravid women (n = 14) on a BeadChip platform. Transient widespread hypomethylation was found in early pregnancy as compared with the non-pregnant states. Methylation of nine genes was significantly different in early pregnancy compared with both postpartum and nulligravid states (< 10% False Discovery Rate). Early pregnancy may be characterized by widespread hypomethylation compared with non-pregnant states; there is no apparent permanent methylation imprint after a normal term gestation. Nine potential candidate genes were identified as differentially methylated in early pregnancy and may play a role in the maternal adaptation to pregnancy.     

Normal early pregnancy: A transient state of epigenetic change favoring hypomethylation

The objective of this study was to analyze genome-wide differential methylation patterns in maternal leukocyte DNA in early pregnant and non-pregnant states. This is an age and body mass index matched case-control study comparing the methylation patterns of 27,578 cytosine-guanine (CpG) sites in 14,495 genes in maternal leukocyte DNA in early pregnancy (n = 14), in the same women postpartum (n = 14), and in nulligravid women (n = 14) on a BeadChip platform. Transient widespread hypomethylation was found in early pregnancy as compared with the non-pregnant states. Methylation of nine genes was significantly different in early pregnancy compared with both postpartum and nulligravid states (< 10% False Discovery Rate). Early pregnancy may be characterized by widespread hypomethylation compared with non-pregnant states; there is no apparent permanent methylation imprint after a normal term gestation. Nine potential candidate genes were identified as differentially methylated in early pregnancy and may play a role in the maternal adaptation to pregnancy.     

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.     

In utero arsenic exposure and epigenome-wide associations in placenta,umbilical artery,and human umbilical vein endothelial cells

Exposure to arsenic early in life has been associated with increased risk of several chronic diseases and is believed to alter epigenetic programming in utero. In the present study, we evaluate the epigenome-wide association of arsenic exposure in utero and DNA methylation in placenta (n = 37), umbilical artery (n = 45) and human umbilical vein endothelial cells (HUVEC) (n = 52) in a birth cohort using the Infinium HumanMethylation450 BeadChip array. Unadjusted and cell mixture adjusted associations for each tissue were examined along with enrichment analyses relative to CpG island location and omnibus permutation tests of association among biological pathways. One CpG in artery (cg26587014) and 4 CpGs in placenta (cg12825509; cg20554753; cg23439277; cg21055948) reached a Bonferroni adjusted level of significance. Several CpGs were differentially methylated in artery and placenta when controlling the false discovery rate (q-value<0.05), but none in HUVEC. Enrichment of hypomethylated CpG islands was observed for artery while hypermethylation of open sea regions were present in placenta relative to prenatal arsenic exposure. The melanogenesis pathway was differentially methylated in artery (Max F P < 0.001), placenta (Max F P < 0.001), and HUVEC (Max F P = 0.02). Similarly, the insulin-signaling pathway was differentially methylated in artery (Max F P = 0.02), placenta (Max F P = 0.02), and HUVEC (Max F P = 0.02). Our results show that prenatal arsenic exposure can alter DNA methylation in artery and placenta but not in HUVEC. Further studies are needed to determine if these alterations in DNA methylation mediate the effect of prenatal arsenic exposure and health outcomes later in life.     

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.     

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.     

Maternal genome-wide DNA methylation patterns and congenital heart defects

The majority of congenital heart defects (CHDs) are thought to result from the interaction between multiple genetic, epigenetic, environmental, and lifestyle factors. Epigenetic mechanisms are attractive targets in the study of complex diseases because they may be altered by environmental factors and dietary interventions. We conducted a population based, case-control study of genome-wide maternal DNA methylation to determine if alterations in gene-specific methylation were associated with CHDs. Using the Illumina Infinium Human Methylation27 BeadChip, we assessed maternal gene-specific methylation in over 27,000 CpG sites from DNA isolated from peripheral blood lymphocytes. Our study sample included 180 mothers with non-syndromic CHD-affected pregnancies (cases) and 187 mothers with unaffected pregnancies (controls). Using a multi-factorial statistical model, we observed differential methylation between cases and controls at multiple CpG sites, although no CpG site reached the most stringent level of genome-wide statistical significance. The majority of differentially methylated CpG sites were hypermethylated in cases and located within CpG islands. Gene Set Enrichment Analysis (GSEA) revealed that the genes of interest were enriched in multiple biological processes involved in fetal development. Associations with canonical pathways previously shown to be involved in fetal organogenesis were also observed. We present preliminary evidence that alterations in maternal DNA methylation may be associated with CHDs. Our results suggest that further studies involving maternal epigenetic patterns and CHDs are warranted. Multiple candidate processes and pathways for future study have been identified.     

Genome wide analysis of acute myeloid leukemia reveal leukemia specific methylome and subtype specific hypomethylation of repeats

Methylated DNA immunoprecipitation followed by high-throughput sequencing (MeDIP-seq) has the potential to identify changes in DNA methylation important in cancer development. In order to understand the role of epigenetic modulation in the development of acute myeloid leukemia (AML) we have applied MeDIP-seq to the DNA of 12 AML patients and 4 normal bone marrows. This analysis revealed leukemia-associated differentially methylated regions that included gene promoters, gene bodies, CpG islands and CpG island shores. Two genes (SPHKAP and DPP6) with significantly methylated promoters were of interest and further analysis of their expression showed them to be repressed in AML. We also demonstrated considerable cytogenetic subtype specificity in the methylomes affecting different genomic features. Significantly distinct patterns of hypomethylation of certain interspersed repeat elements were associated with cytogenetic subtypes. The methylation patterns of members of the SINE family tightly clustered all leukemic patients with an enrichment of Alu repeats with a high CpG density (P<0.0001). We were able to demonstrate significant inverse correlation between intragenic interspersed repeat sequence methylation and gene expression with SINEs showing the strongest inverse correlation (R(2) = 0.7). We conclude that the alterations in DNA methylation that accompany the development of AML affect not only the promoters, but also the non-promoter genomic features, with significant demethylation of certain interspersed repeat DNA elements being associated with AML cytogenetic subtypes. MeDIP-seq data were validated using bisulfite pyrosequencing and the Infinium array.     

Genome-wide DNA methylation profiles of maternal peripheral blood and placentas: potential risk factors for preeclampsia and validation of <Emphasis Type="Italic">GRK5</Emphasis>

Hypoxic placentation has been considered as a key step for the development of preeclampsia (PE); however, the underlying epigenetic mechanisms are still not fully understood. The purpose of this study is to investigate the whole genome DNA methylation status of PE. A microarray analysis using the Infinium HumanMethylation450 BeadChip assay in the placentas and maternal peripheral blood (PB) from PE patients and normal controls was performed. For validation, a quantitative RT-PCR analysis was used. Maternal PB showed 71 differentially methylated CpG loci (44 hypermethylated and 27 hypomethylated), while placenta revealed 365 loci (37 hypermethylated and 328 hypomethylated) at the statistical significance level of |Δβ| ≥ 0.17 and P ≤ 0.01. Notably, among the candidates showing significant signals, GRK5 (a member of G protein-coupled receptor kinase family that has previously been known to be associated with PE) showed a significantly hypomethylated level in the placentas of PE patients (Δβ = ?0.176, P = 2.8 × 10^?5). In the validation for the potential effect of GRK5 methylation on the gene regulation, GRK5 expression was significantly increased in the placentas from PE patients compared to those from controls (P = 0.027). In further GO analysis, genes of MHC class II protein complex showed the most significant differential methylation in the maternal PB of PE patients, while genes of palate development were differentially methylated in the placenta. Although further replication and functional studies are required, our preliminary results suggest that PE has distinct DNA methylation profiles in the maternal PB and placentas, which may provide insight into future research.     

MIRA-SNuPE,a quantitative,multiplex method for measuring allele-specific DNA methylation

5-methyl-C (5mC) and 5-hydroxymethyl-C (5hmC) are epigenetic marks with well-known and putative roles in gene regulation, respectively. These two DNA covalent modifications cannot be distinguished by bisulfite sequencing or restriction digestion, the standard methods of 5mC detection. The methylated CpG island recovery assay (MIRA), however, specifically detects 5mC but not 5hmC. We further developed MIRA for the analysis of allele-specific CpG methylation at differentially methylated regions (DMRs) of imprinted genes. MIRA specifically distinguished between the parental alleles by capturing the paternally methylated H19/Igf2 DMR and maternally methylated KvDMR1 in mouse embryo fibroblasts (MEFs) carrying paternal and maternal duplication of mouse distal Chr7, respectively. MIRA in combination with multiplex single nucleotide primer extension (SNuPE) assays specifically captured the methylated parental allele from normal cells at a set of maternally and paternally methylated DMRs. The assay correctly recognized aberrant biallelic methylation in a case of loss of imprinting. The MIRA-SNuPE assays revealed that placenta exhibited less DNA methylation bias at DMRs compared to yolk sac, amnion, brain, heart, kidney, liver and muscle. This method should be useful for the analysis of allele-specific methylation events related to genomic imprinting, X chromosome inactivation and for verifying and screening haplotype-associated methylation differences in the human population.Key words: epigenetics, imprinting, DMR, MIRA, MBD, DNA methylation, SNuPE     

A novel imprinted transgene located near a repetitive element that exhibits allelic imbalance in DNA methylation during early development

A mouse line carrying a lacZ transgene driven by the human EEF1A1/EF1alpha promoter was established. Although the promoter is known to show ubiquitous activity, only paternal transgene alleles were expressed, resulting in a transgene imprinting. At mid‐gestation, the promoter sequence was differentially methylated, hypomethylated for paternal and hypermethylated for maternal alleles. In germline, the promoter was a typical differentially methylated region. After fertilization, however, both alleles were hypermethylated. Thus, the differential methylation of the promoter required for transgene imprinting was re‐established during later embryonic development independently of the germline differential methylation. Furthermore, also a retroelement promoter closely‐flanking imprinted transgene and its wild type counterpart displayed similar differential methylation during early development. The retroelement promoter was methylated differentially also in germline, but in an opposite pattern to the embryonic differential methylation. These results suggest that there might be an unknown epigenetic regulation inducing transgene imprinting independently of DNA methylation in the transgene insertion site. Then, besides CpG dinucleotides, non‐CpG cytosines of the retroelement promoter were highly methylated especially in the transgene‐active mid‐gestational embryos, suggesting that an unusual epigenetic regulation might protect the active transgene against de novo methylation occurring generally in mid‐gestational embryo.