SWAN: Subset-quantile within array normalization for illumina infinium HumanMethylation450 BeadChips

Cadmium exposure and age-associated DNA methylation changes in non-smoking women from northern Thailand

DNA methylation changes with age, and may serve as a biomarker of aging. Cadmium (Cd) modifies cellular processes that promote aging and disrupts methylation globally. Whether Cd modifies aging processes by influencing establishment of age-associated methylation marks is currently unknown. In this pilot study, we characterized methylation profiles in > 450 000 CpG sites in 40 non-smoking women (age 40-80) differentially exposed to environmental Cd from Thailand. Based on specific gravity adjusted urinary Cd, we classified them as high (HE) and low (LE) exposed and age-matched within 5 years. Urinary Cd was defined as below 2 µg/l in the LE group. We predicted epigenetic age (DNAm-age) using two published methods by Horvath and Hannum and examined the difference between epigenetic age and chronologic age (Δage). We assessed differences by Cd exposure using linear mixed models adjusted for estimated white blood cell proportions, BMI, and urinary creatinine. We identified 213 age-associated CpG sites in our population (P < 10^-4). Counterintuitively, the mean Δage was smaller in HE vs. LE (Hannum: 3.6 vs. 7.6 years, P = 0.0093; Horvath: 2.4 vs. 4.5 years, P = 0.1308). The Cd exposed group was associated with changes in methylation (P < 0.05) at 12, 8, and 20 age-associated sites identified in our population, Hannum, and Horvath. From the results of this pilot study, elevated Cd exposure is associated with methylation changes at age-associated sites and smaller differences between DNAm-age and chronologic age, in contrast to expected age-accelerating effects. Cd may modify epigenetic aging, and biomarkers of aging warrant further investigation when examining Cd and its relationship with chronic disease and mortality.

Differential DNA methylation and lymphocyte proportions in a Costa Rican high longevity region

Background

The Nicoya Peninsula in Costa Rica has one of the highest old-age life expectancies in the world, but the underlying biological mechanisms of this longevity are not well understood. As DNA methylation is hypothesized to be a component of biological aging, we focused on this malleable epigenetic mark to determine its association with current residence in Nicoya versus elsewhere in Costa Rica. Examining a population’s unique DNA methylation pattern allows us to differentiate hallmarks of longevity from individual stochastic variation. These differences may be characteristic of a combination of social, biological, and environmental contexts.

Methods

In a cross-sectional subsample of the Costa Rican Longevity and Healthy Aging Study, we compared whole blood DNA methylation profiles of residents from Nicoya (n = 48) and non-Nicoya (other Costa Rican regions, n = 47) using the Infinium HumanMethylation450 microarray.

Results

We observed a number of differences that may be markers of delayed aging, such as bioinformatically derived differential CD8+ T cell proportions. Additionally, both site- and region-specific analyses revealed DNA methylation patterns unique to Nicoyans. We also observed lower overall variability in DNA methylation in the Nicoyan population, another hallmark of younger biological age.

Conclusions

Nicoyans represent an interesting group of individuals who may possess unique immune cell proportions as well as distinct differences in their epigenome, at the level of DNA methylation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0128-2) contains supplementary material, which is available to authorized users.

A peripheral epigenetic signature of immune system genes is linked to neocortical thickness and memory

Increasing age is tightly linked to decreased thickness of the human neocortex. The biological mechanisms that mediate this effect are hitherto unknown. The DNA methylome, as part of the epigenome, contributes significantly to age-related phenotypic changes. Here, we identify an epigenetic signature that is associated with cortical thickness (P=3.86 × 10⁻⁸) and memory performance in 533 healthy young adults. The epigenetic effect on cortical thickness was replicated in a sample comprising 596 participants with major depressive disorder and healthy controls. The epigenetic signature mediates partially the effect of age on cortical thickness (P<0.001). A multilocus genetic score reflecting genetic variability of this signature is associated with memory performance (P=0.0003) in 3,346 young and elderly healthy adults. The genomic location of the contributing methylation sites points to the involvement of specific immune system genes. The decomposition of blood methylome-wide patterns bears considerable potential for the study of brain-related traits.

Cortical thickness has high heritability estimates and is known to be influenced by genetic factors. Here, Freytag and colleagues show that DNA methylation patterns of peripheral blood monocytes are also correlated with cortical thickness and memory performance in human.

Cord blood hematopoietic cells from preterm infants display altered DNA methylation patterns

Background

Premature infants are highly vulnerable to infection. This is partly attributable to the preterm immune system, which differs from that of the term neonate in cell composition and function. Multiple studies have found differential DNA methylation (DNAm) between preterm and term infants’ cord blood; however, interpretation of these studies is limited by the confounding factor of blood cell composition. This study evaluates the epigenetic impact of preterm birth in isolated hematopoietic cell populations, reducing the concern of cell composition differences.

Methods

Genome-wide DNAm was measured using the Illumina 450K array in T cells, monocytes, granulocytes, and nucleated red blood cells (nRBCs) isolated from cord blood of 5 term and 5 preterm (<31 weeks gestational age) newborns. DNAm of hematopoietic cells was compared globally across the 450K array and through site-specific linear modeling.

Results

Nucleated red blood cells (nRBCs) showed the most extensive changes in DNAm, with 9258 differentially methylated (DM) sites (FDR < 5%, |Δβ| > 0.10) discovered between preterm and term infants compared to the <1000 prematurity-DM sites identified in white blood cell populations. The direction of DNAm change with gestational age at these prematurity-DM sites followed known patterns of hematopoietic differentiation, suggesting that term hematopoietic cell populations are more epigenetically mature than their preterm counterparts. Consistent shifts in DNAm between preterm and term cells were observed at 25 CpG sites, with many of these sites located in genes involved in growth and proliferation, hematopoietic lineage commitment, and the cytoskeleton. DNAm in preterm and term hematopoietic cells conformed to previously identified DNAm signatures of fetal liver and bone marrow, respectively.

Conclusions

This study presents the first genome-wide mapping of epigenetic differences in hematopoietic cells across the late gestational period. DNAm differences in hematopoietic cells between term and <31 weeks were consistent with the hematopoietic origin of these cells during ontogeny, reflecting an important role of DNAm in their regulation. Due to the limited sample size and the high coincidence of prematurity and multiple births, the relationship between cause of preterm birth and DNAm could not be evaluated. These findings highlight gene regulatory mechanisms at both cell-specific and systemic levels that may be involved in fetal immune system maturation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0339-1) contains supplementary material, which is available to authorized users.

DNA methylation profiling in peripheral lung tissues of smokers and patients with COPD

BACKGROUND

Epigenetics changes have been shown to be affected by cigarette smoking. Cigarette smoke (CS)-mediated DNA methylation can potentially affect several cellular and pathophysiological processes, acute exacerbations, and comorbidity in the lungs of patients with chronic obstructive pulmonary disease (COPD). We sought to determine whether genome-wide lung DNA methylation profiles of smokers and patients with COPD were significantly different from non-smokers. We isolated DNA from parenchymal lung tissues of patients including eight lifelong non-smokers, eight current smokers, and eight patients with COPD and analyzed the samples using Illumina's Infinium HumanMethylation450 BeadChip.

RESULTS

Our data revealed that the differentially methylated genes were related to top canonical pathways (e.g., G beta gamma signaling, mechanisms of cancer, and nNOS signaling in neurons), disease and disorders (organismal injury and abnormalities, cancer, and respiratory disease), and molecular and cellular functions (cell death and survival, cellular assembly and organization, cellular function and maintenance) in patients with COPD. The genome-wide DNA methylation analysis identified suggestive genes, such as NOS1AP, TNFAIP2, BID, GABRB1, ATXN7, and THOC7 with DNA methylation changes in COPD lung tissues that were further validated by pyrosequencing. Pyrosequencing validation confirmed hyper-methylation in smokers and patients with COPD as compared to non-smokers. However, we did not detect significant differences in DNA methylation for TNFAIP2, ATXN7, and THOC7 genes in smokers and COPD groups despite the changes observed in the genome-wide analysis.

CONCLUSIONS

Our study suggests that DNA methylation in suggestive genes, such as NOS1AP, BID, and GABRB1 may be used as epigenetic signatures in smokers and patients with COPD if the same is validated in a larger cohort. Future studies are required to correlate DNA methylation status with transcriptomics of selective genes identified in this study and elucidate their role and involvement in the progression of COPD and its exacerbations.

Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14⁺ monocytes, CD16⁺ neutrophils, and naive CD4⁺ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

•

Genome, transcriptome, and epigenome reference panel in three human immune cell types

•

Identified 4,418 genes associated with epigenetic changes independent of genetics

•

Described genome-epigenome coordination defining cell-type-specific regulatory events

•

Functionally mapped disease mechanisms at 345 unique autoimmune disease loci

Aberrant DNA Methylation in Human iPSCs Associates with MYC-Binding Motifs in a Clone-Specific Manner Independent of Genetics

Induced pluripotent stem cells (iPSCs) show variable methylation patterns between lines, some of which reflect aberrant differences relative to embryonic stem cells (ESCs). To examine whether this aberrant methylation results from genetic variation or non-genetic mechanisms, we generated human iPSCs from monozygotic twins to investigate how genetic background, clone, and passage number contribute. We found that aberrantly methylated CpGs are enriched in regulatory regions associated with MYC protein motifs and affect gene expression. We classified differentially methylated CpGs as being associated with genetic and/or non-genetic factors (clone and passage), and we found that aberrant methylation preferentially occurs at CpGs associated with clone-specific effects. We further found that clone-specific effects play a strong role in recurrent aberrant methylation at specific CpG sites across different studies. Our results argue that a non-genetic biological mechanism underlies aberrant methylation in iPSCs and that it is likely based on a probabilistic process involving MYC that takes place during or shortly after reprogramming.

Genome-Wide Measures of Peripheral Blood Dna Methylation and Prostate Cancer Risk in a Prospective Nested Case-Control Study

BACKGROUND

Global measures of peripheral blood DNA methylation have been associated with risk of some malignancies, including breast, bladder, and gastric cancer. Here, we examined genome-wide measures of peripheral blood DNA methylation in prostate cancer and its non-aggressive and aggressive disease forms.

METHODS

We used a matched, case-control study of 687 incident prostate cancer samples, nested within a larger prospective cohort study. DNA methylation was measured in pre-diagnostic, peripheral blood samples using the Illumina Infinium HM450K BeadChip. Genome-wide measures of DNA methylation were computed as the median M-value of all CpG sites and according to CpG site location and regulatory function. We used conditional logistic regression to test for associations between genome-wide measures of DNA methylation and risk of prostate cancer and its subtypes, and by time between blood draw and diagnosis.

RESULTS

We observed no associations between the genome-wide measure of DNA methylation based on all CpG sites and risk of prostate cancer or aggressive disease. Risk of non-aggressive disease was associated with higher methylation of CpG islands (OR = 0.80; 95%CI = 0.68-0.94), promoter regions (OR = 0.79; 95%CI = 0.66-0.93), and high density CpG regions (OR = 0.80; 95%CI = 0.68-0.94). Additionally, higher methylation of all CpGs (OR = 0.66; 95%CI = 0.48-0.89), CpG shores (OR = 0.62; 95%CI = 0.45-0.84), and regulatory regions (OR = 0.68; 95% CI = 0.51-0.91) was associated with a reduced risk of overall prostate cancer within 5 years of blood draw but not thereafter.

CONCLUSIONS

A reduced risk of overall prostate cancer within 5 years of blood draw and non-aggressive prostate cancer was associated with higher genome-wide methylation of peripheral blood DNA. While these data have no immediate clinical utility, with further work they may provide insight into the early events of prostate carcinogenesis. Prostate 77:471-478, 2017. © 2017 Wiley Periodicals, Inc.

Sex differences in DNA methylation of the cord blood are related to sex-bias psychiatric diseases

Sex differences in the prevalence of psychiatric disorders are well documented, with exposure to stress during gestation differentially impacting females and males. We explored sex-specific DNA methylation in the cord blood of 39 females and 32 males born at term and with appropriate weight at birth regarding their potential connection to psychiatric outcomes. Mothers were interviewed to gather information about environmental factors (gestational exposure) that could interfere with the methylation profiles in the newborns. Bisulphite converted DNA was hybridized to Illumina HumanMethylation450 BeadChips. Excluding XYS probes, there were 2,332 differentially methylated CpG sites (DMSs) between sexes, which were enriched within brain modules of co-methylated CpGs during brain development and also differentially methylated in the brains of boys and girls. Genes associated with the DMSs were enriched for neurodevelopmental disorders, particularly for CpG sites found differentially methylated in brain tissue between patients with schizophrenia and controls. Moreover, the DMS had an overlap of 890 (38%) CpG sites with a cohort submitted to toxic exposition during gestation. This study supports the evidences that sex differences in DNA methylation of autosomes act as a primary driver of sex differences that are found in psychiatric outcomes.

Combinatorial Ranking of Gene Sets to Predict Disease Relapse: The Retinoic Acid Pathway in Early Prostate Cancer

Background

Quantitative high-throughput data deposited in consortia such as International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) present opportunities and challenges for computational analyses.

Methods

We present a computational strategy to systematically rank and investigate a large number (2¹⁰–2²⁰) of clinically testable gene sets, using combinatorial gene subset generation and disease-free survival (DFS) analyses. This approach integrates protein–protein interaction networks, gene expression, DNA methylation, and copy number data, in association with DFS profiles from patient clinical records.

Results

As a case study, we applied this pipeline to systematically analyze the role of ALDH1A2 in prostate cancer (PCa). We have previously found this gene to have multiple roles in disease and homeostasis, and here we investigate the role of the associated ALDH1A2 gene/protein networks in PCa, using our methodology in combination with PCa patient clinical profiles from ICGC and TCGA databases. Relationships between gene signatures and relapse were analyzed using Kaplan–Meier (KM) log-rank analysis and multivariable Cox regression. Relative expression versus pooled mean from diploid population was used for z-statistics calculation. Gene/protein interaction network analyses generated 11 core genes associated with ALDH1A2; combinatorial ranking of the power set of these core genes identified two gene sets (out of 2¹¹ − 1 = 2,047 combinations) with significant correlation with disease relapse (KM log rank p < 0.05). For the more significant of these two sets, referred to as the optimal gene set (OGS), patients have median survival 62.7 months with OGS alterations compared to >150 months without OGS alterations (p = 0.0248, hazard ratio = 2.213, 95% confidence interval = 1.1–4.098). Two genes comprising OGS (CYP26A1 and RDH10) are strongly associated with ALDH1A2 in the retinoic acid (RA) pathways, suggesting a major role of RA signaling in early PCa progression. Our pipeline complements human expertise in the search for prognostic biomarkers in large-scale datasets.

Epigenetic modifications and glucocorticoid sensitivity in Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS)

Background

Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) is a debilitating idiopathic disease characterized by unexplained fatigue that fails to resolve with sufficient rest. Diagnosis is based on a list of symptoms and exclusion of other fatigue-related health conditions. Despite a heterogeneous patient population, immune and hypothalamic-pituitary-adrenal (HPA) axis function differences, such as enhanced negative feedback to glucocorticoids, are recurring findings in ME/CFS studies. Epigenetic modifications, such as CpG methylation, are known to regulate long-term phenotypic differences and previous work by our group found DNA methylome differences in ME/CFS, however the relationship between DNA methylome modifications, clinical and functional characteristics associated with ME/CFS has not been examined.

Methods

We examined the DNA methylome in peripheral blood mononuclear cells (PBMCs) of a larger cohort of female ME/CFS patients using the Illumina HumanMethylation450 BeadChip Array. In parallel to the DNA methylome analysis, we investigated in vitro glucocorticoid sensitivity differences by stimulating PBMCs with phytohaemagglutinin and suppressed growth with dexamethasone. We explored DNA methylation differences using bisulfite pyrosequencing and statistical permutation. Linear regression was implemented to discover epigenomic regions associated with self-reported quality of life and network analysis of gene ontology terms to biologically contextualize results.

Results

We detected 12,608 differentially methylated sites between ME/CFS patients and healthy controls predominantly localized to cellular metabolism genes, some of which were also related to self-reported quality of life health scores. Among ME/CFS patients, glucocorticoid sensitivity was associated with differential methylation at 13 loci.

Conclusions

Our results indicate DNA methylation modifications in cellular metabolism in ME/CFS despite a heterogeneous patient population, implicating these processes in immune and HPA axis dysfunction in ME/CFS. Modifications to epigenetic loci associated with differences in glucocorticoid sensitivity may be important as biomarkers for future clinical testing. Overall, these findings align with recent ME/CFS work that point towards impairment in cellular energy production in this patient population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-017-0248-3) contains supplementary material, which is available to authorized users.

Variable DAXX gene methylation is a common feature of placental trophoblast differentiation, preeclampsia, and response to hypoxia

Placental functioning relies on the appropriate differentiation of progenitor villous cytotrophoblasts (CTBs) into extravillous cytotrophoblasts (EVCTs), including invasive EVCTs, and the multinucleated syncytiotrophoblast (ST) layer. This is accompanied by a general move away from a proliferative, immature phenotype. Genome-scale expression studies have provided valuable insight into genes that are associated with the shift to both an invasive EVCT and ST phenotype, whereas genome-scale DNA methylation analysis has shown that differentiation to ST involves widespread methylation shifts, which are counteracted by low oxygen. In the current study, we sought to identify DNA methylation variation that is associated with transition from CTB to ST in vitro and from a noninvasive to invasive EVCT phenotype after culture on Matrigel. Of the several hundred differentially methylated regions that were identified in each comparison, the majority showed a loss of methylation with differentiation. This included a large differentially methylated region (DMR) in the gene body of death domain-associated protein 6 (DAXX ), which lost methylation during both CTB syncytialization to ST and EVCT differentiation to invasive EVCT. Comparison to publicly available methylation array data identified the same DMR as among the most consistently differentially methylated genes in placental samples from preeclampsia pregnancies. Of interest, in vitro culture of CTB or ST in low oxygen increases methylation in the same region, which correlates with delayed differentiation. Analysis of combined epigenomics signatures confirmed DAXX DMR as a likely regulatory element, and direct gene expression analysis identified a positive association between methylation at this site and DAXX expression levels. The widespread dynamic nature of DAXX methylation in association with trophoblast differentiation and placenta-associated pathologies is consistent with an important role for this gene in proper placental development and function.-Novakovic, B., Evain-Brion, D., Murthi, P., Fournier, T., Saffery, R. Variable DAXX gene methylation is a common feature of placental trophoblast differentiation, preeclampsia, and response to hypoxia.

Identification of Predictive DNA Methylation Biomarkers for Chemotherapy Response in Colorectal Cancer

Resistance to 5-Fluorouracil (5-FU) is a major obstacle to the successful treatment of colorectal cancer (CRC) and posed an increased risk of recurrence. DNA methylation has been suggested as one of the underlying mechanisms for recurrent disease and its contribution to the development of drug resistance remains to be clarified. This study aimed to determine the methylation phenotype in CRC for identification of predictive markers for chemotherapy response. We performed DNA methylation profiling on 43 non-recurrent and five recurrent CRC patients using the Illumina Infinium HumanMethylation450 Beadchip assay. In addition, CRC cells with different genetic backgrounds, response to 5-FU and global methylation levels (HT29 and SW48) were treated with 5-FU and DNA methylation inhibitor 5-aza-2′-deoxycytidine (5-azadC). The singular and combined effects of these two drug classes on cell viability and global methylation profiles were investigated. Our genome-wide methylation study on the clinical specimens showed that recurrent CRCs exhibited higher methylation levels compared to non-recurrent CRCs. We identified 4787 significantly differentially methylated genes (P < 0.05); 3112 genes were hyper- while 1675 genes were hypomethylated in the recurrent group compared to the non-recurrent. Fifty eight and 47 of the significantly hypermethylated and hypomethylated genes have an absolute recurrent/non-recurrent methylation difference of ≥20%. Most of the hypermethylated genes were involved in the MAPK signaling pathway which is a key regulator for apoptosis while the hypomethylated genes were involved in the PI3K-AKT signaling pathway and proliferation process. We also demonstrate that 5-azadC treatment enhanced response to 5-FU which resulted in significant growth inhibition compared to 5-FU alone in hypermethylated cell lines SW48. In conclusion, we found the evidence of five potentially biologically important genes in recurrent CRCs that could possibly serve as a new potential therapeutic targets for patients with chemoresistance. We postulate that aberrant methylation of CCNEI, CCNDBP1, PON3, DDX43, and CHL1 in CRC might be associated with the recurrence of CRC and 5-azadC-mediated restoration of 5-FU sensitivity is mediated at least in part by MAPK signaling pathway.

Immunomethylomic approach to explore the blood neutrophil lymphocyte ratio (NLR) in glioma survival

Background

Differentially methylated regions (DMRs) within DNA isolated from whole blood can be used to estimate the proportions of circulating leukocyte subtypes. We use the term “immunomethylomics” to describe the application of these immune lineage DMRs to studying leukocyte profiles. Here, we applied this approach to peripheral blood DNA from 72 glioma patients with molecularly defined brain tumors, representing common patient groups with defined characteristic survival times and risk factors. We first estimated the proportions of leukocyte subtypes in samples using deconvolution algorithms with reference DMR libraries from isolated leukocyte populations and Illumina 450K DNA methylation data. Then, we calculated the neutrophil to lymphocyte ratio (NLR) using methylation-derived cell composition estimates (mdNLR). The NLR is considered an indicator of immunosuppressive cells in cancer patients.

Results

Elevated mdNLR scores were observed in glioma patients compared to mdNLR values of published controls. Significantly decreased survival times were associated with mdNLR ≥ 4.0 in Cox proportional hazards models adjusted for age, gender, tumor grade, and molecular subtype (HR 2.02, 95% CI, 1.11–3.69). We also identified five myeloid-related CpGs that were highly correlated with the mdNLR (adjusted R² ≥ 0.80). Each of the five myeloid CpG loci was associated with survival when adjusted for the above covariates and offer a simplified approach for utilizing fresh or archived peripheral blood samples for interrogating a very small number of methylation markers to estimate myeloid immune influences in glioma survival.

Conclusions

The mdNLR (based on DNA methylation) is a novel candidate methylation biomarker that represents immunosuppressive myeloid cells within the blood of glioma patients with potential application in clinical trials and future epidemiologic studies of glioma risk and survival.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0316-8) contains supplementary material, which is available to authorized users.

In utero exposure to maternal smoking is associated with DNA methylation alterations and reduced neuronal content in the developing fetal brain

Background

Intrauterine exposure to maternal smoking is linked to impaired executive function and behavioral problems in the offspring. Maternal smoking is associated with reduced fetal brain growth and smaller volume of cortical gray matter in childhood, indicating that prenatal exposure to tobacco may impact cortical development and manifest as behavioral problems. Cellular development is mediated by changes in epigenetic modifications such as DNA methylation, which can be affected by exposure to tobacco.

Results

In this study, we sought to ascertain how maternal smoking during pregnancy affects global DNA methylation profiles of the developing dorsolateral prefrontal cortex (DLPFC) during the second trimester of gestation. When DLPFC methylation profiles (assayed via Illumina, HM450) of smoking-exposed and unexposed fetuses were compared, no differentially methylated regions (DMRs) passed the false discovery correction (FDR ≤ 0.05). However, the most significant DMRs were hypomethylated CpG Islands within the promoter regions of GNA15 and SDHAP3 of smoking-exposed fetuses. Interestingly, the developmental up-regulation of SDHAP3 mRNA was delayed in smoking-exposed fetuses. Interaction analysis between gestational age and smoking exposure identified significant DMRs annotated to SYCE3, C21orf56/LSS, SPAG1 and RNU12/POLDIP3 that passed FDR. Furthermore, utilizing established methods to estimate cell proportions by DNA methylation, we found that exposed DLPFC samples contained a lower proportion of neurons in samples from fetuses exposed to maternal smoking. We also show through in vitro experiments that nicotine impedes the differentiation of neurons independent of cell death.

Conclusions

We found evidence that intrauterine smoking exposure alters the developmental patterning of DNA methylation and gene expression and is associated with reduced mature neuronal content, effects that are likely driven by nicotine.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0111-y) contains supplementary material, which is available to authorized users.

Genome-wide analysis of differential transcriptional and epigenetic variability across human immune cell types

Background

A healthy immune system requires immune cells that adapt rapidly to environmental challenges. This phenotypic plasticity can be mediated by transcriptional and epigenetic variability.

Results

We apply a novel analytical approach to measure and compare transcriptional and epigenetic variability genome-wide across CD14⁺CD16⁻ monocytes, CD66b⁺CD16⁺ neutrophils, and CD4⁺CD45RA⁺ naïve T cells from the same 125 healthy individuals. We discover substantially increased variability in neutrophils compared to monocytes and T cells. In neutrophils, genes with hypervariable expression are found to be implicated in key immune pathways and are associated with cellular properties and environmental exposure. We also observe increased sex-specific gene expression differences in neutrophils. Neutrophil-specific DNA methylation hypervariable sites are enriched at dynamic chromatin regions and active enhancers.

Conclusions

Our data highlight the importance of transcriptional and epigenetic variability for the key role of neutrophils as the first responders to inflammatory stimuli. We provide a resource to enable further functional studies into the plasticity of immune cells, which can be accessed from: http://blueprint-dev.bioinfo.cnio.es/WP10/hypervariability.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1156-8) contains supplementary material, which is available to authorized users.

META2: Intercellular DNA Methylation Pairwise Annotation and Integrative Analysis

Genome-wide deciphering intercellular differential DNA methylation as well as its roles in transcriptional regulation remains elusive in cancer epigenetics. Here we developed a toolkit META2 for DNA methylation annotation and analysis, which aims to perform integrative analysis on differentially methylated loci and regions through deep mining and statistical comparison methods. META2 contains multiple versatile functions for investigating and annotating DNA methylation profiles. Benchmarked with T-47D cell, we interrogated the association within differentially methylated CpG (DMC) and region (DMR) candidate count and region length and identified major transition zones as clues for inferring statistically significant DMRs; together we validated those DMRs with the functional annotation. Thus META2 can provide a comprehensive analysis approach for epigenetic research and clinical study.

Combining omics data to identify genes associated with allergic rhinitis

Allergic rhinitis is a common chronic disorder characterized by immunoglobulin E-mediated inflammation. To identify new genes associated with this trait, we performed genome- and epigenome-wide association studies and linked marginally significant CpGs located in genes or its promoter and SNPs located 1 Mb from the CpGs, by identifying cis methylation quantitative trait loci (mQTL). This approach relies on functional cellular aspects rather than stringent statistical correction. We were able to identify one gene with significant cis-mQTL for allergic rhinitis, caudal-type homeobox 1 (CDX1). We also identified 11 genes with marginally significant cis-mQTLs (p < 0.05) including one with both allergic rhinitis with or without asthma (RNF39). Moreover, most SNPs identified were not located closest to the gene they were linked to through cis-mQTLs counting the one linked to CDX1 located in a gene previously associated with asthma and atopic dermatitis. By combining omics data, we were able to identify new genes associated with allergic rhinitis and better assess the genes linked to associated SNPs.

Interleukin-2 receptor-α proximal promoter hypomethylation is associated with multiple sclerosis

Genetic studies have demonstrated association between single-nucleotide polymorphisms within the IL2RA (interleukin-2 receptor α-subunit) gene and risk of developing multiple sclerosis (MS); however, these variants do not have obvious functional consequences. DNA methylation is a source of genetic variation that could impact on autoimmune disease risk. We investigated DNA methylation of the IL2RA promoter in genomic DNA obtained from peripheral blood mononuclear cells and neural tissue using matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry. A differential methylation profile of IL2RA was identified, suggesting that IL2RA expression was regulated by DNA methylation. We extended our analysis of DNA methylation to peripheral blood mononuclear cell (PBMC) of MS cases and controls using MALDI-TOF and Illumina HumanMethylation450 arrays. Analyses of CpG sites within the proximal promoter of IL2RA in PBMC showed no differences between MS cases and controls despite an increase in IL2RA expression. In contrast, we inferred significant DNA methylation differences specific to particular leukocyte subtypes in MS cases compared with controls by deconvolution of the array data. The decrease in methylation in patients correlated with an increase in IL2RA expression in T cells from MS cases in comparison with controls. Our data suggest that differential methylation of the IL2RA promoter in T cells could be an important pathogenic mechanism in MS.

Genome-wide methylation data mirror ancestry information

Background

Genetic data are known to harbor information about human demographics, and genotyping data are commonly used for capturing ancestry information by leveraging genome-wide differences between populations. In contrast, it is not clear to what extent population structure is captured by whole-genome DNA methylation data.

Results

We demonstrate, using three large-cohort 450K methylation array data sets, that ancestry information signal is mirrored in genome-wide DNA methylation data and that it can be further isolated more effectively by leveraging the correlation structure of CpGs with cis-located SNPs. Based on these insights, we propose a method, EPISTRUCTURE, for the inference of ancestry from methylation data, without the need for genotype data.

Conclusions

EPISTRUCTURE can be used to infer ancestry information of individuals based on their methylation data in the absence of corresponding genetic data. Although genetic data are often collected in epigenetic studies of large cohorts, these are typically not made publicly available, making the application of EPISTRUCTURE especially useful for anyone working on public data. Implementation of EPISTRUCTURE is available in GLINT, our recently released toolset for DNA methylation analysis at: http://glint-epigenetics.readthedocs.io.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0108-y) contains supplementary material, which is available to authorized users.

RELIC: a novel dye-bias correction method for Illumina Methylation BeadChip

BACKGROUND

The Illumina Infinium HumanMethylation450 BeadChip and its successor, Infinium MethylationEPIC BeadChip, have been extensively utilized in epigenome-wide association studies. Both arrays use two fluorescent dyes (Cy3-green/Cy5-red) to measure methylation level at CpG sites. However, performance difference between dyes can result in biased estimates of methylation levels.

RESULTS

Here we describe a novel method, called REgression on Logarithm of Internal Control probes (RELIC) to correct for dye bias on whole array by utilizing the intensity values of paired internal control probes that monitor the two color channels. We evaluate the method in several datasets against other widely used dye-bias correction methods. Results on data quality improvement showed that RELIC correction statistically significantly outperforms alternative dye-bias correction methods. We incorporated the method into the R package ENmix, which is freely available from the Bioconductor website ( https://www.bioconductor.org/packages/release/bioc/html/ENmix.html ).

CONCLUSIONS

RELIC is an efficient and robust method to correct for dye-bias in Illumina Methylation BeadChip data. It outperforms other alternative methods and conveniently implemented in R package ENmix to facilitate DNA methylation studies.

Epigenome-Wide Tumor DNA Methylation Profiling Identifies Novel Prognostic Biomarkers of Metastatic-Lethal Progression in Men Diagnosed with Clinically Localized Prostate Cancer

PURPOSE

Aside from Gleason sum, few factors accurately identify the subset of prostate cancer patients at high risk for metastatic progression. We hypothesized that epigenetic alterations could distinguish prostate tumors with life-threatening potential.

EXPERIMENTAL DESIGN

Epigenome-wide DNA methylation profiling was performed in surgically resected primary tumor tissues from a population-based (n = 430) and a replication (n = 80) cohort of prostate cancer patients followed prospectively for at least 5 years. Metastasis was confirmed by positive bone scan, MRI, CT, or biopsy, and death certificates confirmed cause of death. AUC, partial AUC (pAUC, 95% specificity), and P value criteria were used to select differentially methylated CpG sites that robustly stratify patients with metastatic-lethal from nonrecurrent tumors, and which were complementary to Gleason sum.

RESULTS

Forty-two CpG biomarkers stratified patients with metastatic-lethal versus nonrecurrent prostate cancer in the discovery cohort, and eight of these CpGs replicated in the validation cohort based on a significant (P < 0.05) AUC (range, 0.66-0.75) or pAUC (range, 0.007-0.009). The biomarkers that improved discrimination of patients with metastatic-lethal prostate cancer include CpGs in five genes (ALKBH5, ATP11A, FHAD1, KLHL8, and PI15) and three intergenic regions. In the validation dataset, the AUC for Gleason sum alone (0.82) significantly increased with the addition of four individual CpGs (range, 0.86-0.89; all P <0.05).

CONCLUSIONS

Eight differentially methylated CpGs that distinguish patients with metastatic-lethal from nonrecurrent tumors were validated. These novel epigenetic biomarkers warrant further investigation as they may improve prognostic classification of patients with clinically localized prostate cancer and provide new insights on tumor aggressiveness. Clin Cancer Res; 23(1); 311-9. ©2016 AACR.

DNA methylation landscape of hepatoblastomas reveals arrest at early stages of liver differentiation and cancer-related alterations

Hepatoblastomas are uncommon embryonal liver tumors accounting for approximately 80% of childhood hepatic cancer. We hypothesized that epigenetic changes, including DNA methylation, could be relevant to hepatoblastoma onset. The methylomes of eight matched hepatoblastomas and non-tumoral liver tissues were characterized, and data were validated in an independent group (11 hepatoblastomas). In comparison to differentiated livers, hepatoblastomas exhibited a widespread and non-stochastic pattern of global low-level hypomethylation. The analysis revealed 1,359 differentially methylated CpG sites (DMSs) between hepatoblastomas and control livers, which are associated with 765 genes. Hypomethylation was detected in hepatoblastomas for ~58% of the DMSs with enrichment at intergenic sites, and most of the hypermethylated CpGs were located in CpG islands. Functional analyses revealed enrichment in signaling pathways involved in metabolism, negative regulation of cell differentiation, liver development, cancer, and Wnt signaling pathway. Strikingly, an important overlap was observed between the 1,359 DMSs and the CpG sites reported to exhibit methylation changes through liver development (p<0.0001), with similar patterns of methylation in both hepatoblastomas and fetal livers compared to adult livers. Overall, our results suggest an arrest at early stages of liver cell differentiation, in line with the hypothesis that hepatoblastoma ontogeny involves the disruption of liver development. This genome-wide methylation dysfunction, taken together with a relatively small number of driver genetic mutations reported for both adult and pediatric liver cancers, shed light on the relevance of epigenetic mechanisms for hepatic tumorigenesis.

DNA Methylation-Derived Neutrophil-to-Lymphocyte Ratio: An Epigenetic Tool to Explore Cancer Inflammation and Outcomes

Background: The peripheral blood neutrophil-to-lymphocyte ratio (NLR) is a cytologic marker of both inflammation and poor outcomes in patients with cancer. DNA methylation is a key element of the epigenetic program defining different leukocyte subtypes and may provide an alternative to cytology in assessing leukocyte profiles. Our aim was to create a bioinformatic tool to estimate NLR using DNA methylation, and to assess its diagnostic and prognostic performance in human populations.Methods: We developed a DNA methylation-derived NLR (mdNLR) index based on normal isolated leukocyte methylation libraries and established cell-mixture deconvolution algorithms. The method was applied to cancer case-control studies of the bladder, head and neck, ovary, and breast, as well as publicly available data on cancer-free subjects.Results: Across cancer studies, mdNLR scores were either elevated in cases relative to controls, or associated with increased hazard of death. High mdNLR values (>5) were strong indicators of poor survival. In addition, mdNLR scores were elevated in males, in nonHispanic white versus Hispanic ethnicity, and increased with age. We also observed a significant interaction between cigarette smoking history and mdNLR on cancer survival.Conclusions: These results mean that our current understanding of mature leukocyte methylomes is sufficient to allow researchers and clinicians to apply epigenetically based analyses of NLR in clinical and epidemiologic studies of cancer risk and survival.Impact: As cytologic measurements of NLR are not always possible (i.e., archival blood), mdNLR, which is computed from DNA methylation signatures alone, has the potential to expand the scope of epigenome-wide association studies. Cancer Epidemiol Biomarkers Prev; 26(3); 328-38. ©2016 AACR.

DNA methylation in lung cells is associated with asthma endotypes and genetic risk

The epigenome provides a substrate through which environmental exposures can exert their effects on gene expression and disease risk, but the relative importance of epigenetic variation on human disease onset and progression is poorly characterized. Asthma is a heterogeneous disease of the airways, for which both onset and clinical course result from interactions between host genotype and environmental exposures, yet little is known about the molecular mechanisms for these interactions. We assessed genome-wide DNA methylation using the Infinium Human Methylation 450K Bead Chip and characterized the transcriptome by RNA sequencing in primary airway epithelial cells from 74 asthmatic and 41 nonasthmatic adults. Asthma status was based on doctor's diagnosis and current medication use. Genotyping was performed using various Illumina platforms. Our study revealed a regulatory locus on chromosome 17q12-21 associated with asthma risk and epigenetic signatures of specific asthma endotypes and molecular networks. Overall, these data support a central role for DNA methylation in lung cells, which promotes distinct molecular pathways of asthma pathogenesis and modulates the effects of genetic variation on disease risk and clinical heterogeneity.

Common epigenetic variation in a European population of mentally healthy young adults

DNA methylation represents an important link between structural genetic variation and complex phenotypes. The study of genome-wide CpG methylation and its relation to traits relevant to psychiatry has become increasingly important. Here, we analyzed quality metrics of 394,043 CpG sites in two samples of 568 and 319 mentally healthy young adults. For 25% of all CpGs we observed medium to large common epigenetic variation. These CpGs were overrepresented in open sea and shore regions, as well as in intergenic regions. They also showed a strong enrichment of significant hits in association analyses. Furthermore, a significant proportion of common DNA methylation is at least partially genetically driven and thus may be observed similarly across tissues. These findings could be of particular relevance for studies of complex neuropsychiatric traits, which often rely on proxy tissues.

DNA methylation at diagnosis is associated with response to disease-modifying drugs in early rheumatoid arthritis

AIM

A proof-of-concept study to explore whether DNA methylation at first diagnosis is associated with response to disease-modifying antirheumatic drugs (DMARDs) in patients with early rheumatoid arthritis (RA).

PATIENTS & METHODS

DNA methylation was quantified in T-lymphocytes from 46 treatment-naive patients using HumanMethylation450 BeadChips. Treatment response was determined in 6 months using the European League Against Rheumatism (EULAR) response criteria.

RESULTS

Initial filtering identified 21 cytosine-phosphate-guanines (CpGs) that were differentially methylated between responders and nonresponders. After conservative adjustment for multiple testing, six sites remained statistically significant, of which four showed high sensitivity and/or specificity (≥75%) for response to treatment. Moreover, methylation at two sites in combination was the strongest factor associated with response (80.0% sensitivity, 90.9% specificity, AUC 0.85).

CONCLUSION

DNA methylation at diagnosis is associated with disease-modifying antirheumatic drug treatment response in early RA.

Differentially Methylated DNA Regions in Monozygotic Twin Pairs Discordant for Rheumatoid Arthritis: An Epigenome-Wide Study

Objectives

In an explorative epigenome-wide association study (EWAS) to search for gene independent, differentially methylated DNA positions and regions (DMRs) associated with rheumatoid arthritis (RA) by studying monozygotic (MZ) twin pairs discordant for RA.

Methods

Genomic DNA was isolated from whole blood samples from 28 MZ twin pairs discordant for RA. DNA methylation was measured using the HumanMethylation450 BeadChips. Smoking, anti-cyclic citrullinated peptide antibodies, and immunosuppressive treatment were included as covariates. Pathway analysis was performed using GREAT.

Results

Smoking was significantly associated with hypomethylation of a DMR overlapping the promoter region of the RNF5 and the AGPAT1, which are implicated in inflammation and autoimmunity, whereas DMARD treatment induced hypermethylation of the same region. Additionally, the promotor region of both S100A6 and EFCAB4B were hypomethylated, and both genes have previously been associated with RA. We replicated several candidate genes identified in a previous EWAS in treatment-naïve RA singletons. Gene-set analysis indicated the involvement of immunologic signatures and cancer-related pathways in RA.

Conclusion

We identified several differentially methylated regions associated with RA, which may represent environmental effects or consequences of the disease and plausible biological pathways pertinent to the pathogenesis of RA.

Gene-Specific DNA Methylation Changes Predict Remission in Patients with ANCA-Associated Vasculitis

ANCA-associated vasculitis is an autoimmune condition characterized by vascular inflammation and organ damage. Pharmacologically induced remission of this condition is complicated by relapses. Potential triggers of relapse are immunologic challenges and environmental insults, both of which associate with changes in epigenetic silencing modifications. Altered histone modifications implicated in gene silencing associate with aberrant autoantigen expression. To establish a link between DNA methylation, a model epigenetic gene silencing modification, and autoantigen gene expression and disease status in ANCA-associated vasculitis, we measured gene-specific DNA methylation of the autoantigen genes myeloperoxidase (MPO) and proteinase 3 (PRTN3) in leukocytes of patients with ANCA-associated vasculitis observed longitudinally (n=82) and of healthy controls (n=32). Patients with active disease demonstrated hypomethylation of MPO and PRTN3 and increased expression of the autoantigens; in remission, DNA methylation generally increased. Longitudinal analysis revealed that patients with ANCA-associated vasculitis could be divided into two groups, on the basis of whether DNA methylation increased or decreased from active disease to remission. In patients with increased DNA methylation, MPO and PRTN3 expression correlated with DNA methylation. Kaplan-Meier estimate of relapse revealed patients with increased DNA methylation at the PRTN3 promoter had a significantly greater probability of a relapse-free period (P<0.001), independent of ANCA serotype. Patients with decreased DNA methylation at the PRTN3 promoter had a greater risk of relapse (hazard ratio, 4.55; 95% confidence interval, 2.09 to 9.91). Thus, changes in the DNA methylation status of the PRTN3 promoter may predict the likelihood of stable remission and explain autoantigen gene regulation.

Racial differences in genome-wide methylation profiling and gene expression in breast tissues from healthy women

Breast cancer is more common in European Americans (EAs) than in African Americans (AAs) but mortality from breast cancer is higher among AAs. While there are racial differences in DNA methylation and gene expression in breast tumors, little is known whether such racial differences exist in breast tissues of healthy women. Genome-wide DNA methylation and gene expression profiling was performed in histologically normal breast tissues of healthy women. Linear regression models were used to identify differentially-methylated CpG sites (CpGs) between EAs (n = 61) and AAs (n = 22). Correlations for methylation and expression were assessed. Biological functions of the differentially-methylated genes were assigned using the Ingenuity Pathway Analysis. Among 485 differentially-methylated CpGs by race, 203 were hypermethylated in EAs, and 282 were hypermethylated in AAs. Promoter-related differentially-methylated CpGs were more frequently hypermethylated in EAs (52%) than AAs (27%) while gene body and intergenic CpGs were more frequently hypermethylated in AAs. The differentially-methylated CpGs were enriched for cancer-associated genes with roles in cell death and survival, cellular development, and cell-to-cell signaling. In a separate analysis for correlation in EAs and AAs, different patterns of correlation were found between EAs and AAs. The correlated genes showed different biological networks between EAs and AAs; networks were connected by Ubiquitin C. To our knowledge, this is the first comprehensive genome-wide study to identify differences in methylation and gene expression between EAs and AAs in breast tissues from healthy women. These findings may provide further insights regarding the contribution of epigenetic differences to racial disparities in breast cancer.

Differential DNA methylation in peripheral blood mononuclear cells in adolescents exposed to significant early but not later childhood adversity

Internationally adopted adolescents who are adopted as young children from conditions of poverty and deprivation have poorer physical and mental health outcomes than do adolescents conceived, born, and raised in the United States by families similar to those who adopt internationally. Using a sample of Russian and Eastern European adoptees to control for Caucasian race and US birth, and nonadopted offspring of well-educated and well-resourced parents to control for postadoption conditions, we hypothesized that the important differences in environments, conception to adoption, might be reflected in epigenetic patterns between groups, specifically in DNA methylation. Thus, we conducted an epigenome-wide association study to compare DNA methylation profiles at approximately 416,000 individual CpG loci from peripheral blood mononuclear cells of 50 adopted youth and 33 nonadopted youth. Adopted youth averaged 22 months at adoption, and both groups averaged 15 years at testing; thus, roughly 80% of their lives were lived in similar circumstances. Although concurrent physical health did not differ, cell-type composition predicted using the DNA methylation data revealed a striking difference in the white blood cell-type composition of the adopted and nonadopted youth. After correcting for cell type and removing invariant probes, 30 CpG sites in 19 genes were more methylated in the adopted group. We also used an exploratory functional analysis that revealed that 223 gene ontology terms, clustered in neural and developmental categories, were significantly enriched between groups.

The nasal methylome and childhood atopic asthma

BACKGROUND

Given the strong environmental influence on both epigenetic marks and allergic asthma in children, the epigenetic alterations in respiratory epithelia might provide insight into allergic asthma.

OBJECTIVE

We sought to identify DNA methylation and gene expression changes associated with childhood allergic persistent asthma.

METHODS

We compared genomic DNA methylation patterns and gene expression in African American children with persistent atopic asthma (n = 36) versus healthy control subjects (n = 36). Results were validated in an independent population of asthmatic children (n = 30) by using a shared healthy control population (n = 36) and in an independent population of white adult atopic asthmatic patients (n = 12) and control subjects (n = 12).

RESULTS

We identified 186 genes with significant methylation changes, differentially methylated regions or differentially methylated probes, after adjustment for age, sex, race/ethnicity, batch effects, inflation, and multiple comparisons. Genes differentially methylated included those with established roles in asthma and atopy and genes related to extracellular matrix, immunity, cell adhesion, epigenetic regulation, and airflow obstruction. The methylation changes were substantial (median, 9.5%; range, 2.6% to 29.5%). Hypomethylated and hypermethylated genes were associated with increased and decreased gene expression, respectively (P < 2.8 × 10^-6 for differentially methylated regions and P < 7.8 × 10^-10 for differentially methylated probes). Quantitative analysis in 53 differentially expressed genes demonstrated that 32 (60%) have significant methylation-expression relationships within 5 kb of the gene. Ten loci selected based on the relevance to asthma, magnitude of methylation change, and methylation-expression relationships were validated in an independent cohort of children with atopic asthma. Sixty-seven of 186 genes also have significant asthma-associated methylation changes in nasal epithelia of adult white asthmatic patients.

CONCLUSIONS

Epigenetic marks in respiratory epithelia are associated with allergic asthma and gene expression changes in inner-city children.

Statistical challenges in analyzing methylation and long-range chromosomal interaction data

With the rapid development of high throughput technologies such as array and next generation sequencing (NGS), genome-wide, nucleotide-resolution epigenomic data are increasingly available. In recent years, there has been particular interest in data on DNA methylation and 3-dimensional (3D) chromosomal organization, which are believed to hold keys to understand biological mechanisms, such as transcription regulation, that are closely linked to human health and diseases. However, small sample size, complicated correlation structure, substantial noise, biases, and uncertainties, all present difficulties for performing statistical inference. In this review, we present an overview of the new technologies that are frequently utilized in studying DNA methylation and 3D chromosomal organization. We focus on reviewing recent developments in statistical methodologies designed for better interrogating epigenomic data, pointing out statistical challenges facing the field whenever appropriate.

Characterizing the hypomethylated DNA methylation profile of nucleated red blood cells from cord blood

AIM

To provide insight into fetal nucleated red blood cell (nRBC) development using genome-wide DNA methylation (DNAm) profiling.

MATERIALS & METHODS

The DNAm profile (Illumina 450K array) of cord blood (n = 7) derived nRBCs was compared with B cells, CD4 and CD8 T cells, natural killer cells, granulocytes, monocytes and placenta (n = 5).

RESULTS

nRBCs and placenta had similarly low array-wide DNAm compared with white blood cells, but their patterns of hypomethylation differed at biologically relevant subsets of the array. High interindividual variability in nRBC DNAm was driven by a negative association between DNAm and nRBC count.

CONCLUSION

nRBC hypomethylation is likely an epigenetic signature of erythropoiesis rather than of early development. Variability in nRBC DNAm may stem from differences in the cell population's maturity or hematopoietic source.

Epigenetic signature of Gleason score and prostate cancer recurrence after radical prostatectomy

BACKGROUND

Identifying the subset of patients with clinically localized prostate cancer (PCa) at the highest risk of recurrence remains challenging, and better prognostic markers are needed. Gleason score is the best predictor of PCa aggressiveness and prognosis. In the present study, we generated an epigenetic signature based on high versus low Gleason score tumors and evaluated its ability to predict recurrence after radical prostatectomy.

METHODS

Genome-wide DNA methylation data from The Cancer Genome Atlas (TCGA; no. of patients = 333) and the elastic net method were used to generate an epigenetic signature by contrasting patients with high (8-10) versus low (≤6) Gleason score tumors. The signature was then tested in a cohort of 523 patients with clinically localized disease who had radical prostatectomy. Samples taken from the primary tumor were used for DNA methylation and mRNA expression profiling. Patients were followed for PCa recurrence on average for 8 years after diagnosis.

RESULTS

The epigenetic signature includes 52 differentially methylated CpG sites. In the testing cohort, the signature was associated with poorer recurrence-free survival (hazard ratio per 25 % increase = 1.78; 95 % confidence interval 1.48, 2.16). The signature significantly improved the area under the curve (AUC) for PCa recurrence compared to clinical-pathological parameters alone, particularly among patients diagnosed with Gleason score 7 tumors (0.64 vs. 0.76, P = 1.34E-4). Results were comparable for patients with Gleason 3 + 4 and those with 4 + 3 tumors. Gene Set Enrichment Analysis showed that higher levels of the signature were associated with increased expression of genes related to cell cycle proliferation and decreased expression of androgen-responsive genes.

CONCLUSIONS

This report shows evidence that DNA methylation patterns measured in prostate tumor cells are predictive of PCa aggressiveness. The epigenetic signature may have clinical utility to improve prognostication particularly in patients with intermediate Gleason score 7 tumors.

Comparative DNA Methylation Profiling Reveals an Immunoepigenetic Signature of HIV-related Cognitive Impairment

Monocytes/macrophages contribute to the neuropathogenesis of HIV-related cognitive impairment (CI); however, considerable gaps in our understanding of the precise mechanisms driving this relationship remain. Furthermore, whether a distinct biological profile associated with HIV-related CI resides in immune cell populations remains unknown. Here, we profiled DNA methylomes and transcriptomes of monocytes derived from HIV-infected individuals with and without CI using genome-wide DNA methylation and gene expression profiling. We identified 1,032 CI-associated differentially methylated loci in monocytes. These loci related to gene networks linked to the central nervous system (CNS) and interactions with HIV. Most (70.6%) of these loci exhibited higher DNA methylation states in the CI group and were preferentially distributed over gene bodies and intergenic regions of the genome. CI-associated DNA methylation states at 12 CpG sites associated with neuropsychological testing performance scores. CI-associated DNA methylation also associated with gene expression differences including CNS genes CSRNP1 (P = 0.017), DISC1 (P = 0.012), and NR4A2 (P = 0.005); and a gene known to relate to HIV viremia, THBS1 (P = 0.003). This discovery cohort data unveils cell type-specific DNA methylation patterns related to HIV-associated CI and provide an immunoepigenetic DNA methylation “signature” potentially useful for corroborating clinical assessments, informing pathogenic mechanisms, and revealing new therapeutic targets against CI.

Analysis of the breast cancer methylome using formalin-fixed paraffin-embedded tumour

PURPOSE

Aberrant DNA methylation occurs frequently in breast carcinogenesis. Tools for translational epigenetic studies of breast cancer involving formalin-fixed paraffin-embedded (FFPE) human tissues have now been developed. Few studies have measured genome-wide methylation in DNA derived from paraffin-embedded tumour tissues and compared the DNA methylation in corresponding adjacent non-tumour ductal epithelium (ADJ_NT). These studies are technically challenging due to the spectrum of breast cancer pathologies, the variable suitability of DNA extracted from FFPE material and the difficulties in identifying ADJ_NT. We assessed the suitability of FFPE breast cancer material for genome-wide DNA methylation assessment of tumour and ADJ_NT.

METHODS

Twenty-one archival breast tumour tissues with paired ADJ_NT obtained from separate blocks and at least 2 cm from the tumour were sourced from The Melbourne Collaborative Cohort Study (MCCS). DNA was prepared from macrodissected tissue samples and assessed for genome-wide methylation using the Infinium HumanMethylation450 Beadchip (HM450K) array.

RESULTS

The 1000 most differentially methylated probes between tumour and ADJ_NT in this FFPE-derived dataset differentiated tumour and ADJ_NT in The Cancer Genome Atlas Network data (TCGA; derived from high molecular weight DNA using the same HM450K array).

CONCLUSIONS

Large-scale studies of genome-wide DNA methylation using FFPE breast cancer specimens offer the opportunity to further refine the pathological classification of tumours, to include subtypes that are underrepresented in the TCGA data and provide the capacity to further explore intra-tumoural heterogeneity.

Genome-wide DNA methylation analysis in obsessive-compulsive disorder patients

Literatures have suggested that not only genetic but also environmental factors, interactively accounted for susceptibility of obsessive-compulsive disorder (OCD). DNA methylation may regulate expression of genes as the heritable epigenetic modification. The examination for genome-wide DNA methylation was performed on blood samples from 65 patients with OCD, as well as 96 healthy control subjects. The DNA methylation was examined at over 485,000 CpG sites using the Illumina Infinium Human Methylation450 BeadChip. As a result, 8,417 probes corresponding to 2,190 unique genes were found to be differentially methylated between OCD and healthy control subjects. Of those genes, 4,013 loci were located in CpG islands and 2,478 were in promoter regions. These included BCYRN1, BCOR, FGF13, HLA-DRB1, ARX, etc., which have previously been reported to be associated with OCD. Pathway analyses indicated that regulation of actin cytoskeleton, cell adhesion molecules (CAMs), actin binding, transcription regulator activity, and other pathways might be further associated with risk of OCD. Unsupervised clustering analysis of the top 3,000 most variable probes revealed two distinct groups with significantly more people with OCD in cluster one compared with controls (67.74% of cases v.s. 27.13% of controls, Chi-square = 26.011, df = 1, P = 3.41E-07). These results strongly suggested that differential DNA methylation might play an important role in etiology of OCD.

DMAK: A curated pan-cancer DNA methylation annotation knowledgebase

Pan-cancer analysis can identify cell- and tissue-specific genomic loci and regions with underlying biological functions. Here we present an online curated DNA Methylation Annotation Knowledgebase, DMAK, which includes the pan-cancer analysis results for differentially-methylated loci and regions by the Reduced Representation Bisulfite Sequencing profiling technology. DMAK contains 3 modules of curated information and analysis results on 688,445 CpG sites across 19 cancer and embryonic stem cell lines from ENCODE, and further analysis of survival associations with clinical sources retrieved from TCGA. The knowledgebase covers all identified differentially-methylated CpG sites and regions of interest, further annotated genomic information, together with tumor suppressor genes information and calculated methylation level. DMAK provides meaningful clues for deriving functional association network and related clinical association results based on protein-coding genes, including tumor suppressor genes, identified from differentially methylated regions of interest. Thus DMAK constitutes a comprehensive reference source for the current epigenetic research and clinical study.

Epigenetic drift in the aging genome: a ten-year follow-up in an elderly twin cohort

BACKGROUND

Current epigenetic studies on aging are dominated by the cross-sectional design that correlates subjects' ages or age groups with their measured epigenetic profiles. Such studies have been more aimed at age prediction or building up the epigenetic clock of age rather than focusing on the dynamic patterns in epigenetic changes during the aging process.

METHODS

We performed an epigenome-wide association study of intra-individual longitudinal changes in DNA methylation at CpG (cytosine-phosphate-guanine) sites measured in whole-blood samples of a cohort of 43 elderly twin pairs followed for 10 years (age at intake 73-82 years). Biological pathway analysis and survival analysis were also conducted on CpGs showing longitudinal change in their DNA-methylation levels. Classical twin models were fitted to each CpG site to estimate the genetic and environmental effects on DNA-methylation.

RESULTS

Our analysis identified 2284 CpG sites whose DNA-methylation levels changed longitudinally over the follow-up. Twin modelling revealed that the longitudinal change for 90% of these CpG sites was explained solely by individual unique environmental factors and only for 10% of these sites was it influenced by familial factors (genetic or shared environment). Over 60% of the identified CpG sites were replicated (same direction and replication P < 0.05) in an independent cross-sectional sample of 300 twins aged from 30 to 74 years. The replication rate went up to 91% for the top 53 CpGs with P < 1 × 10-07. Pathway analysis of genes linked to these CpGs identified biologically meaningful gene-sets involved in cellular-signalling events and in transmission across chemical synapses, which are important molecular underpinnings of aging-related degenerative disorders.

CONCLUSION

Our epigenome-wide association studies on a cohort of old twins followed up for 10 years identified highly replicable epigenetic biomarkers predominantly implicated in signalling pathways of degenerative disorders and survival in the elderly.

Genome-wide measures of DNA methylation in peripheral blood and the risk of urothelial cell carcinoma: a prospective nested case-control study

Background:

Global DNA methylation has been reported to be associated with urothelial cell carcinoma (UCC) by studies using blood samples collected at diagnosis. Using the Illumina HumanMethylation450 assay, we derived genome-wide measures of blood DNA methylation and assessed them for their prospective association with UCC risk.

Methods:

We used 439 case–control pairs from the Melbourne Collaborative Cohort Study matched on age, sex, country of birth, DNA sample type, and collection period. Conditional logistic regression was used to compute odds ratios (OR) of UCC risk per s.d. of each genome-wide measure of DNA methylation and 95% confidence intervals (CIs), adjusted for potential confounders. We also investigated associations by disease subtype, sex, smoking, and time since blood collection.

Results:

The risk of superficial UCC was decreased for individuals with higher levels of our genome-wide DNA methylation measure (OR=0.71, 95% CI: 0.54–0.94; P=0.02). This association was particularly strong for current smokers at sample collection (OR=0.47, 95% CI: 0.27–0.83). Intermediate levels of our genome-wide measure were associated with decreased risk of invasive UCC. Some variation was observed between UCC subtypes and the location and regulatory function of the CpGs included in the genome-wide measures of methylation.

Conclusions:

Higher levels of our genome-wide DNA methylation measure were associated with decreased risk of superficial UCC and intermediate levels were associated with reduced risk of invasive disease. These findings require replication by other prospective studies.

Genome-wide characteristics of de novo mutations in autism

De novo mutations (DNMs) are important in Autism Spectrum Disorder (ASD), but so far analyses have mainly been on the ~1.5% of the genome encoding genes. Here, we performed whole genome sequencing (WGS) of 200 ASD parent-child trios and characterized germline and somatic DNMs. We confirmed that the majority of germline DNMs (75.6%) originated from the father, and these increased significantly with paternal age only (p=4.2×10^-10). However, when clustered DNMs (those within 20kb) were found in ASD, not only did they mostly originate from the mother (p=7.7×10^-13), but they could also be found adjacent to de novo copy number variations (CNVs) where the mutation rate was significantly elevated (p=2.4×10^-24). By comparing DNMs detected in controls, we found a significant enrichment of predicted damaging DNMs in ASD cases (p=8.0×10^-9; OR=1.84), of which 15.6% (p=4.3×10^-3) and 22.5% (p=7.0×10^-5) were in the non-coding or genic non-coding, respectively. The non-coding elements most enriched for DNM were untranslated regions of genes, boundaries involved in exon-skipping and DNase I hypersensitive regions. Using microarrays and a novel outlier detection test, we also found aberrant methylation profiles in 2/185 (1.1%) of ASD cases. These same individuals carried independently identified DNMs in the ASD risk- and epigenetic- genes DNMT3A and ADNP. Our data begins to characterize different genome-wide DNMs, and highlight the contribution of non-coding variants, to the etiology of ASD.

Epigenetic dysregulation in the developing Down syndrome cortex

Using Illumina 450K arrays, 1.85% of all analyzed CpG sites were significantly hypermethylated and 0.31% hypomethylated in fetal Down syndrome (DS) cortex throughout the genome. The methylation changes on chromosome 21 appeared to be balanced between hypo- and hyper-methylation, whereas, consistent with prior reports, all other chromosomes showed 3-11 times more hyper- than hypo-methylated sites. Reduced NRSF/REST expression due to upregulation of DYRK1A (on chromosome 21q22.13) and methylation of REST binding sites during early developmental stages may contribute to this genome-wide excess of hypermethylated sites. Upregulation of DNMT3L (on chromosome 21q22.4) could lead to de novo methylation in neuroprogenitors, which then persists in the fetal DS brain where DNMT3A and DNMT3B become downregulated. The vast majority of differentially methylated promoters and genes was hypermethylated in DS and located outside chromosome 21, including the protocadherin gamma (PCDHG) cluster on chromosome 5q31, which is crucial for neural circuit formation in the developing brain. Bisulfite pyrosequencing and targeted RNA sequencing showed that several genes of PCDHG subfamilies A and B are hypermethylated and transcriptionally downregulated in fetal DS cortex. Decreased PCDHG expression is expected to reduce dendrite arborization and growth in cortical neurons. Since constitutive hypermethylation of PCDHG and other genes affects multiple tissues, including blood, it may provide useful biomarkers for DS brain development and pharmacologic targets for therapeutic interventions.

Reliability of DNA methylation measures from dried blood spots and mononuclear cells using the HumanMethylation450k BeadArray

The reliability of methylation measures from the widely used HumanMethylation450 (HM450K) microarray has not been assessed for DNA from dried blood spots (DBS) or peripheral blood mononuclear cells (PBMC), nor for combined data from different studies. Repeated HM450K methylation measures in DNA from DBS and PBMC samples were available from participants in six case-control studies nested within the Melbourne Collaborative Cohort Study. Reliability was assessed for individual CpGs by calculating the intraclass correlation coefficient (ICC) based on technical replicates (samples repeated in a single study; 126 PBMC, 136 DBS) and study duplicates (samples repeated across studies; 280 PBMC, 769 DBS) using mixed-effects models. Reliability based on technical replicates was moderate for PBMC (median ICC = 0.42), but lower for DBS (median ICC = 0.20). Study duplicates gave lower ICCs than technical replicates. CpGs that were either highly methylated or unmethylated generally had lower ICCs, which appeared to be mostly related to their lower variability. The ICCs for global methylation measures were high, typically greater than 0.70. The reliability of methylation measures determined by the HM450K microarray is wide-ranging and depends primarily on the variability in methylation at individual CpG sites. The power of association studies is low for a substantial proportion of CpGs in the HM450K assay.

CpG sites with continuously increasing or decreasing methylation from early to late human fetal brain development

Normal human brain development is dependent on highly dynamic epigenetic processes for spatial and temporal gene regulation. Recent work identified wide-spread changes in DNA methylation during fetal brain development. We profiled CpG methylation in frontal cortex of 27 fetuses from gestational weeks 12-42, using Illumina 450K methylation arrays. Sites showing genome-wide significant correlation with gestational age were compared to a publicly available data set from gestational weeks 3-26. Altogether, we identified 2016 matching developmentally regulated differentially methylated positions (m-dDMPs): 1767m-dDMPs were hypermethylated and 1149 hypomethylated during fetal development. M-dDMPs are underrepresented in CpG islands and gene promoters, and enriched in gene bodies. They appear to cluster in certain chromosome regions. M-dDMPs are significantly enriched in autism-associated genes and CpGs. Our results promote the idea that reduced methylation dynamics during fetal brain development may predispose to autism. In addition, m-dDMPs are enriched in genes with human-specific brain expression patterns and/or histone modifications. Collectively, we defined a subset of dDMPs exhibiting constant methylation changes from early to late pregnancy. The same epigenetic mechanisms involving methylation changes in cis-regulatory regions may have been adopted for human brain evolution and ontogeny.

Genome-wide pharmacologic unmasking identifies tumor suppressive microRNAs in multiple myeloma

Epigenetic alterations have emerged as an important cause of microRNA (miRNA) deregulation. In Multiple Myeloma (MM), a few tumor suppressive miRNAs silenced by DNA hypermethylation have been reported, but so far there are few systemic investigations on epigenetically silenced miRNAs. We conducted genome-wide screening for tumor suppressive miRNAs epigenetically silenced in MM. Four Human MM Cell lines were treated with demethylating agent 5'azacytidine (5'aza). Consistently upregulated miRNAs include miR-155, miR-198, miR-135a*, miR-200c, miR-125a-3p, miR-188-5p, miR-483-5p, miR-663, and miR-630. Methylation array analysis revealed increased methylation at or near miRNA-associated CpG islands in MM patients. Ectopic restoration of miR-155, miR-198, miR-135a*, miR-200c, miR-663 and miR-483-5p significantly repressed MM cell proliferation, migration and colony formation. Furthermore, we derived a 33-gene signature from predicted miRNA target genes that were also upregulated in MM patients and associated with patient survival in three independent myeloma datasets. In summary, we have revealed important, epigenetically silenced tumor suppressive miRNAs by pharmacologic reversal of epigenetic silencing.

Genome-scale methylation assessment did not identify prognostic biomarkers in oral tongue carcinomas

Background

DNA methylation profiling of heterogeneous head and neck squamous cell carcinoma (HNSCC) cohorts has been reported to predict patient outcome. We investigated if a prognostic DNA methylation profile could be found in tumour tissue from a single uniform subsite, the oral tongue. The methylation status of 109 comprehensively annotated oral tongue squamous cell carcinoma (OTSCC) formalin-fixed paraffin-embedded (FFPE) samples from a single institution were examined with the Illumina HumanMethylation450K (HM450K) array. Data pre-processing, quality control and analysis were performed using R packages. Probes mapping to SNPs, sex chromosomes and unreliable probes were accounted for prior to downstream analyses. The relationship between methylation and patient survival was examined using both agnostic approaches and feature selection. The cohort was enlarged by incorporation of 331 The Cancer Genome Atlas (TCGA) HNSCC samples, which included 91 TCGA OTSCC samples with HM450K and survival data available.

Results

Given the use of FFPE-derived DNA, we defined different cohorts for separate analyses. Overall, similar results were found between cohorts. With an unsupervised approach, no distinct hypermethylated group of samples was identified and nor was a prognostic methylation profile identified. The use of multiple downstream feature selection approaches, including a linear models for microarray data (LIMMA), centroid feature selection (CFS), and recursive feature elimination (RFE) support vector machines, similarly failed to identify a significant methylation signature informative for patient prognosis or any clinicopathological data available. Furthermore, we were unable to confirm the prognostic methylation profiles or specific prognostic loci reported within the literature for HNSCC.

Conclusions

With genome-scale assessment of DNA methylation using HM450K in one of the largest OTSCC cohorts to date, we were unable to identify a hypermethylated group of tumours or a prognostic methylation signature. This suggests that either DNA methylation in isolation is not likely to be of prognostic value or larger cohorts are required to identify such a biomarker for OTSCC.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0235-0) contains supplementary material, which is available to authorized users.

Comparison of pre-processing methodologies for Illumina 450k methylation array data in familial analyses

Background

Human methylome mapping in health and disease states has largely relied on Illumina Human Methylation 450k array (450k array) technology. Accompanying this has been the necessary evolution of analysis pipelines to facilitate data processing. The majority of these pipelines, however, cater for experimental designs where matched ‘controls’ or ‘normal’ samples are available. Experimental designs where no appropriate ‘reference’ exists remain challenging. Herein, we use data generated from our study of the inheritance of methylome profiles in families to evaluate the performance of eight normalisation pre-processing methods. Fifty individual samples representing four families were interrogated on five 450k array BeadChips. Eight normalisation methods were tested using qualitative and quantitative metrics, to assess efficacy and suitability.

Results

Stratified quantile normalisation combined with ComBat were consistently found to be the most appropriate when assessed using density, MDS and cluster plots. This was supported quantitatively by ANOVA on the first principal component where the effect of batch dropped from p < 0.01 to p = 0.97 after stratified QN and ComBat. Median absolute differences between replicated samples were the lowest after stratified QN and ComBat as were the standard error measures on known imprinted regions. Biological information was preserved after normalisation as indicated by the maintenance of a significant association between a known mQTL and methylation (p = 1.05e-05).

Conclusions

A strategy combining stratified QN with ComBat is appropriate for use in the analyses when no reference sample is available but preservation of biological variation is paramount. There is great potential for use of 450k array data to further our understanding of the methylome in a variety of similar settings. Such advances will be reliant on the determination of appropriate methodologies for processing these data such as established here.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0241-2) contains supplementary material, which is available to authorized users.

Age-associated hydroxymethylation in human bone-marrow mesenchymal stem cells

Background

Age-associated changes in genomic DNA methylation have been primarily attributed to 5-methylcytosine (5mC). However, the recent discovery of 5-hydroxymethylcytosine (5hmC) suggests that this epigenetic mark might also play a role in the process.

Methods

Here, we analyzed the genome-wide profile of 5hmc in mesenchymal stem cells (MSCs) obtained from bone-marrow donors, aged 2–89 years.

Results

We identified 10,685 frequently hydroxymethylated CpG sites in MSCs that were, as in other cell types, significantly associated with low density CpG regions, introns, the histone posttranslational modification H3k4me1 and enhancers. Study of the age-associated changes to 5hmC identified 785 hyper- and 846 hypo-hydroxymethylated CpG sites in the MSCs obtained from older individuals.

Conclusions

DNA hyper-hydroxymethylation in the advanced-age group was associated with loss of 5mC, which suggests that, at specific CpG sites, this epigenetic modification might play a role in DNA methylation changes during lifetime. Since bone-marrow MSCs have many clinical applications, and the fact that the epigenomic alterations in this cell type associated with aging identified in this study could have associated functional effects, the age of donors should be taken into account in clinical settings.

Electronic supplementary material

The online version of this article (doi:10.1186/s12967-016-0966-x) contains supplementary material, which is available to authorized users.