Multi-Platform Whole-Genome Microarray Analyses Refine the Epigenetic Signature of Breast Cancer Metastasis with Gene Expression and Copy Number

Background

We have previously identified genome-wide DNA methylation changes in a cell line model of breast cancer metastasis. These complex epigenetic changes that we observed, along with concurrent karyotype analyses, have led us to hypothesize that complex genomic alterations in cancer cells (deletions, translocations and ploidy) are superimposed over promoter-specific methylation events that are responsible for gene-specific expression changes observed in breast cancer metastasis.

Methodology/Principal Findings

We undertook simultaneous high-resolution, whole-genome analyses of MDA-MB-468GFP and MDA-MB-468GFP-LN human breast cancer cell lines (an isogenic, paired lymphatic metastasis cell line model) using Affymetrix gene expression (U133), promoter (1.0R), and SNP/CNV (SNP 6.0) microarray platforms to correlate data from gene expression, epigenetic (DNA methylation), and combination copy number variant/single nucleotide polymorphism microarrays. Using Partek Software and Ingenuity Pathway Analysis we integrated datasets from these three platforms and detected multiple hypomethylation and hypermethylation events. Many of these epigenetic alterations correlated with gene expression changes. In addition, gene dosage events correlated with the karyotypic differences observed between the cell lines and were reflected in specific promoter methylation patterns. Gene subsets were identified that correlated hyper (and hypo) methylation with the loss (or gain) of gene expression and in parallel, with gene dosage losses and gains, respectively. Individual gene targets from these subsets were also validated for their methylation, expression and copy number status, and susceptible gene pathways were identified that may indicate how selective advantage drives the processes of tumourigenesis and metastasis.

Conclusions/Significance

Our approach allows more precisely profiling of functionally relevant epigenetic signatures that are associated with cancer progression and metastasis.     

5-Aza-2′-deoxycytidine Leads to Reduced Embryo Implantation and Reduced Expression of DNA Methyltransferases and Essential Endometrial Genes

Background

The DNA demethylating agent 5-aza-2′-deoxycytidine (5-aza-CdR) incorporates into DNA and decreases DNA methylation, sparking interest in its use as a potential therapeutic agent. We aimed to determine the effects of maternal 5-aza-CdR treatment on embryo implantation in the mouse and to evaluate whether these effects are associated with decreased levels of DNA methyltransferases (Dnmts) and three genes (estrogen receptor α [Esr1], progesterone receptor [Pgr], and homeobox A10 [Hoxa10]) that are vital for control of endometrial changes during implantation.

Methods and Principal Findings

Mice treated with 5-aza-CdR had a dose-dependent decrease in number of implantation sites, with defected endometrial decidualization and stromal cell proliferation. Western blot analysis on pseudo-pregnant day 3 (PD3) showed that 0.1 mg/kg 5-aza-CdR significantly repressed Dnmt3a protein level, and 0.5 mg/kg 5-aza-CdR significantly repressed Dnmt1, Dnmt3a, and Dnmt3b protein levels in the endometrium. On PD5, mice showed significantly decreased Dnmt3a protein level with 0.1 mg/kg 5-aza-CdR, and significantly decreased Dnmt1 and Dnmt3a with 0.5 mg/kg 5-aza-CdR. Immunohistochemical staining showed that 5-aza-CdR repressed DNMT expression in a cell type–specific fashion within the uterus, including decreased expression of Dnmt1 in luminal and/or glandular epithelium and of Dnmt3a and Dnmt3b in stroma. Furthermore, the 5′ flanking regions of the Esr1, Pgr, and Hoxa10 were hypomethylated on PD5. Interestingly, the higher (0.5 mg/kg) dose of 5-aza-CdR decreased protein expression of Esr1, Pgr, and Hoxa10 in the endometrium on PD5 in both methylation-dependent and methylation-independent manners.

Conclusions

The effects of 5-aza-CdR on embryo implantation in mice were associated with altered expression of endometrial Dnmts and genes controlling endometrial changes, suggesting that altered gene methylation, and not cytotoxicity alone, contributes to implantation defects induced by 5-aza-CdR.     

Reproducibility and concordance of differential DNA methylation and gene expression in cancer

Background

Hundreds of genes with differential DNA methylation of promoters have been identified for various cancers. However, the reproducibility of differential DNA methylation discoveries for cancer and the relationship between DNA methylation and aberrant gene expression have not been systematically analysed.

Methodology/Principal Findings

Using array data for seven types of cancers, we first evaluated the effects of experimental batches on differential DNA methylation detection. Second, we compared the directions of DNA methylation changes detected from different datasets for the same cancer. Third, we evaluated the concordance between methylation and gene expression changes. Finally, we compared DNA methylation changes in different cancers. For a given cancer, the directions of methylation and expression changes detected from different datasets, excluding potential batch effects, were highly consistent. In different cancers, DNA hypermethylation was highly inversely correlated with the down-regulation of gene expression, whereas hypomethylation was only weakly correlated with the up-regulation of genes. Finally, we found that genes commonly hypomethylated in different cancers primarily performed functions associated with chronic inflammation, such as ‘keratinization’, ‘chemotaxis’ and ‘immune response’.

Conclusions

Batch effects could greatly affect the discovery of DNA methylation biomarkers. For a particular cancer, both differential DNA methylation and gene expression can be reproducibly detected from different studies with no batch effects. While DNA hypermethylation is significantly linked to gene down-regulation, hypomethylation is only weakly correlated with gene up-regulation and is likely to be linked to chronic inflammation.     

Dietary Flavanols Modulate the Transcription of Genes Associated with Cardiovascular Pathology without Changes in Their DNA Methylation State

Background

In a recent intervention study, the daily supplementation with 200 mg monomeric and oligomeric flavanols (MOF) from grape seeds for 8 weeks revealed a vascular health benefit in male smokers. The objective of the present study was to determine the impact of MOF consumption on the gene expression profile of leukocytes and to assess changes in DNA methylation.

Methodology/Principal Findings

Gene expression profiles were determined using whole genome microarrays (Agilent) and DNA methylation was assessed using HumanMethylation450 BeadChips (Illumina). MOF significantly modulated the expression of 864 genes. The majority of the affected genes are involved in chemotaxis, cell adhesion, cell infiltration or cytoskeleton organisation, suggesting lower immune cell adhesion to endothelial cells. This was corroborated by in vitro experiments showing that MOF exposure of monocytes attenuates their adhesion to TNF-α-stimulated endothelial cells. Nuclear factor kappa B (NF-κB) reporter gene assays confirmed that MOF decrease the activity of NF-κB. Strong inter-individual variability in the leukocytes'' DNA methylation was observed. As a consequence, on group level, changes due to MOF supplementation could not be found.

Conclusion

Our study revealed that an 8 week daily supplementation with 200 mg MOF modulates the expression of genes associated with cardiovascular disease pathways without major changes of their DNA methylation state. However, strong inter-individual variation in leukocyte DNA methylation may obscure the subtle epigenetic response to dietary flavanols. Despite the lack of significant changes in DNA methylation, the modulation of gene expression appears to contribute to the observed vascular health effect of MOF in humans.     

Association of Childhood Chronic Physical Aggression with a DNA Methylation Signature in Adult Human T Cells

Background

Chronic physical aggression (CPA) is characterized by frequent use of physical aggression from early childhood to adolescence. Observed in approximately 5% of males, CPA is associated with early childhood adverse environments and long-term negative consequences. Alterations in DNA methylation, a covalent modification of DNA that regulates genome function, have been associated with early childhood adversity.

Aims

To test the hypothesis that a trajectory of chronic physical aggression during childhood is associated with a distinct DNA methylation profile during adulthood.

Methods

We analyzed genome-wide promoter DNA methylation profiles of T cells from two groups of adult males assessed annually for frequency of physical aggression between 6 and 15 years of age: a group with CPA and a control group. Methylation profiles covering the promoter regions of 20 000 genes and 400 microRNAs were generated using MeDIP followed by hybridization to microarrays.

Results

In total, 448 distinct gene promoters were differentially methylated in CPA. Functionally, many of these genes have previously been shown to play a role in aggression and were enriched in biological pathways affected by behavior. Their locations in the genome tended to form clusters spanning millions of bases in the genome.

Conclusions

This study provides evidence of clustered and genome-wide variation in promoter DNA methylation in young adults that associates with a history of chronic physical aggression from 6 to 15 years of age. However, longitudinal studies of methylation during early childhood will be necessary to determine if and how this methylation variation in T cells DNA plays a role in early development of chronic physical aggression.     

Global DNA hypomethylation in peripheral blood leukocytes as a biomarker for cancer risk: a meta-analysis

Background

Good biomarkers for early detection of cancer lead to better prognosis. However, harvesting tumor tissue is invasive and cannot be routinely performed. Global DNA methylation of peripheral blood leukocyte DNA was evaluated as a biomarker for cancer risk.

Methods

We performed a meta-analysis to estimate overall cancer risk according to global DNA hypomethylation levels among studies with various cancer types and analytical methods used to measure DNA methylation. Studies were systemically searched via PubMed with no language limitation up to July 2011. Summary estimates were calculated using a fixed effects model.

Results

The subgroup analyses by experimental methods to determine DNA methylation level were performed due to heterogeneity within the selected studies (p<0.001, I²: 80%). Heterogeneity was not found in the subgroup of %5-mC (p = 0.393, I²: 0%) and LINE-1 used same target sequence (p = 0.097, I²: 49%), whereas considerable variance remained in LINE-1 (p<0.001, I²: 80%) and bladder cancer studies (p = 0.016, I²: 76%). These results suggest that experimental methods used to quantify global DNA methylation levels are important factors in the association study between hypomethylation levels and cancer risk. Overall, cancer risks of the group with the lowest DNA methylation levels were significantly higher compared to the group with the highest methylation levels [OR (95% CI): 1.48 (1.28–1.70)].

Conclusions

Global DNA hypomethylation in peripheral blood leukocytes may be a suitable biomarker for cancer risk. However, the association between global DNA methylation and cancer risk may be different based on experimental methods, and region of DNA targeted for measuring global hypomethylation levels as well as the cancer type. Therefore, it is important to select a precise and accurate surrogate marker for global DNA methylation levels in the association studies between global DNA methylation levels in peripheral leukocyte and cancer risk.     

Pre-Diagnostic Leukocyte Genomic DNA Methylation and the Risk of Colorectal Cancer in Women

Background

Abnormal one-carbon metabolism may lead to general genomic (global) hypomethylation, which may predispose an individual to the development of colorectal neoplasia.

Methods

We evaluated the association between pre-diagnostic leukocyte genomic DNA methylation level and the risk of colorectal cancer in a nested case-control study of 358 colorectal cancer cases and 661 matched controls within the all-female cohort of the Nurses’ Health Study (NHS). Among control subjects, we further examined major plasma components in the one-carbon metabolism pathway in relation to genomic DNA methylation level. Liquid chromatography/tandem mass spectrometry was used to examine leukocyte genomic DNA methylation level. We calculated odds ratios (ORs) and 95% confidence intervals (95% CIs) using logistic regression.

Results

Overall genomic DNA methylation level was not associated with the risk of colorectal cancer (p for trend, 0.45). Compared with women in the lowest quintile of methylation, the multivariate OR of colorectal cancer risk was 1.32 (95% CI, 0.82–2.13) for those in the highest quintile. We did not find significant associations between major plasma components of one-carbon metabolism or risk factors for colorectal cancer and genomic DNA methylation level (all p for trend >0.05). Also, neither one-carbon metabolism-related plasma components nor well-known risk factors for colorectal cancer modified the association between genomic DNA methylation level and the risk of colorectal cancer (all p for interaction >0.05).

Conclusions

We found no evidence that hypomethylation of leukocyte genomic DNA increases risk of colorectal cancer among women. Additional studies are needed to investigate the association between pre-diagnostic genomic DNA methylation level and colorectal cancer risk among diverse populations.     

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.     

Distinct functional patterns of gene promoter hypomethylation and hypermethylation in cancer genomes

Background

Aberrant DNA methylation plays important roles in carcinogenesis. However, the functional significance of genome-wide hypermethylation and hypomethylation of gene promoters in carcinogenesis currently remain unclear.

Principal Findings

Based on genome-wide methylation data for five cancer types, we showed that genes with promoter hypermethylation were highly consistent in function across different cancer types, and so were genes with promoter hypomethylation. Functions related to “developmental processes” and “regulation of biology processes” were significantly enriched with hypermethylated genes but were depleted of hypomethylated genes. In contrast, functions related to “cell killing” and “response to stimulus”, including immune and inflammatory response, were associated with an enrichment of hypomethylated genes and depletion of hypermethylated genes. We also observed that some families of cytokines secreted by immune cells, such as IL10 family cytokines and chemokines, tended to be hypomethylated in various cancer types. These results provide new hints for understanding the distinct functional roles of genome-wide hypermethylation and hypomethylation of gene promoters in carcinogenesis.

Conclusions

Genes with promoter hypermethylation and hypomethylation are highly consistent in function across different cancer types, respectively, but these two groups of genes tend to be enriched in different functions associated with cancer. Especially, we speculate that hypomethylation of gene promoters may play roles in inducing immunity and inflammation disorders in precancerous conditions, which may provide hints for improving epigenetic therapy and immunotherapy of cancer.     

Quantitative reconstruction of leukocyte subsets using DNA methylation

Background

Cell lineage-specific DNA methylation patterns distinguish normal human leukocyte subsets and can be used to detect and quantify these subsets in peripheral blood. We have developed an approach that uses DNA methylation to simultaneously quantify multiple leukocyte subsets, enabling investigation of immune modulations in virtually any blood sample including archived samples previously precluded from such analysis. Here we assess the performance characteristics and validity of this approach.

Results

Using Illumina Infinium HumanMethylation27 and VeraCode GoldenGate Methylation Assay microarrays, we measure DNA methylation in leukocyte subsets purified from human whole blood and identify cell lineage-specific DNA methylation signatures that distinguish human T cells, B cells, NK cells, monocytes, eosinophils, basophils and neutrophils. We employ a bioinformatics-based approach to quantify these cell types in complex mixtures, including whole blood, using DNA methylation at as few as 20 CpG loci. A reconstruction experiment confirms that the approach could accurately measure the composition of mixtures of human blood leukocyte subsets. Applying the DNA methylation-based approach to quantify the cellular components of human whole blood, we verify its accuracy by direct comparison to gold standard immune quantification methods that utilize physical, optical and proteomic characteristics of the cells. We also demonstrate that the approach is not affected by storage of blood samples, even under conditions prohibiting the use of gold standard methods.

Conclusions

Cell mixture distributions within peripheral blood can be assessed accurately and reliably using DNA methylation. Thus, precise immune cell differential estimates can be reconstructed using only DNA rather than whole cells.