Correcting for cell-type heterogeneity in epigenome-wide association studies: revisiting previous analyses

An improved epigenetic counter to track mitotic age in normal and precancerous tissues

The cumulative number of stem cell divisions in a tissue, known as mitotic age, is thought to be a major determinant of cancer-risk. Somatic mutational and DNA methylation (DNAm) clocks are promising tools to molecularly track mitotic age, yet their relationship is underexplored and their potential for cancer risk prediction in normal tissues remains to be demonstrated. Here we build and validate an improved pan-tissue DNAm counter of total mitotic age called stemTOC. We demonstrate that stemTOC's mitotic age proxy increases with the tumor cell-of-origin fraction in each of 15 cancer-types, in precancerous lesions, and in normal tissues exposed to major cancer risk factors. Extensive benchmarking against 6 other mitotic counters shows that stemTOC compares favorably, specially in the preinvasive and normal-tissue contexts. By cross-correlating stemTOC to two clock-like somatic mutational signatures, we confirm the mitotic-like nature of only one of these. Our data points towards DNAm as a promising molecular substrate for detecting mitotic-age increases in normal tissues and precancerous lesions, and hence for developing cancer-risk prediction strategies.

A meta-analysis of immune-cell fractions at high resolution reveals novel associations with common phenotypes and health outcomes

BACKGROUND

Changes in cell-type composition of tissues are associated with a wide range of diseases and environmental risk factors and may be causally implicated in disease development and progression. However, these shifts in cell-type fractions are often of a low magnitude, or involve similar cell subtypes, making their reliable identification challenging. DNA methylation profiling in a tissue like blood is a promising approach to discover shifts in cell-type abundance, yet studies have only been performed at a relatively low cellular resolution and in isolation, limiting their power to detect shifts in tissue composition.

METHODS

Here we derive a DNA methylation reference matrix for 12 immune-cell types in human blood and extensively validate it with flow-cytometric count data and in whole-genome bisulfite sequencing data of sorted cells. Using this reference matrix, we perform a directional Stouffer and fixed effects meta-analysis comprising 23,053 blood samples from 22 different cohorts, to comprehensively map associations between the 12 immune-cell fractions and common phenotypes. In a separate cohort of 4386 blood samples, we assess associations between immune-cell fractions and health outcomes.

RESULTS

Our meta-analysis reveals many associations of cell-type fractions with age, sex, smoking and obesity, many of which we validate with single-cell RNA sequencing. We discover that naïve and regulatory T-cell subsets are higher in women compared to men, while the reverse is true for monocyte, natural killer, basophil, and eosinophil fractions. Decreased natural killer counts associated with smoking, obesity, and stress levels, while an increased count correlates with exercise and sleep. Analysis of health outcomes revealed that increased naïve CD4 + T-cell and N-cell fractions associated with a reduced risk of all-cause mortality independently of all major epidemiological risk factors and baseline co-morbidity. A machine learning predictor built only with immune-cell fractions achieved a C-index value for all-cause mortality of 0.69 (95%CI 0.67-0.72), which increased to 0.83 (0.80-0.86) upon inclusion of epidemiological risk factors and baseline co-morbidity.

CONCLUSIONS

This work contributes an extensively validated high-resolution DNAm reference matrix for blood, which is made freely available, and uses it to generate a comprehensive map of associations between immune-cell fractions and common phenotypes, including health outcomes.

Cell-Type Deconvolution of Bulk DNA Methylation Data with EpiSCORE

DNA methylation data generated from bulk tissue represents a mixture of many different cell types. Variation in the cell-type composition of tissues is thus a major confounder when inferring differential DNA methylation. Due to the high cost of single-cell methylome sequencing, computational methods that can dissect the cell-type heterogeneity of bulk DNA methylomes offer an efficient and cost-effective solution, especially in the context of large-scale EWAS. In this chapter, we present a step-by-step tutorial of Epigenetic cell-type deconvolution using Single-Cell Omic References (EpiSCORE), a reference-based method that leverages the high-resolution nature of single-cell RNA-Seq datasets to facilitate microdissection of bulk-tissue DNA methylomes.

Targeted Methylation Profiling of Single Laser-Capture Microdissected Post-Mortem Brain Cells by Adapted Limiting Dilution Bisulfite Pyrosequencing (LDBSP)

A reoccurring issue in neuroepigenomic studies, especially in the context of neurodegenerative disease, is the use of (heterogeneous) bulk tissue, which generates noise during epigenetic profiling. A workable solution to this issue is to quantify epigenetic patterns in individually isolated neuronal cells using laser capture microdissection (LCM). For this purpose, we established a novel approach for targeted DNA methylation profiling of individual genes that relies on a combination of LCM and limiting dilution bisulfite pyrosequencing (LDBSP). Using this approach, we determined cytosine-phosphate-guanine (CpG) methylation rates of single alleles derived from 50 neurons that were isolated from unfixed post-mortem brain tissue. In the present manuscript, we describe the general workflow and, as a showcase, demonstrate how targeted methylation analysis of various genes, in this case, RHBDF2, OXT, TNXB, DNAJB13, PGLYRP1, C3, and LMX1B, can be performed simultaneously. By doing so, we describe an adapted data analysis pipeline for LDBSP, allowing one to include and correct CpG methylation rates derived from multi-allele reactions. In addition, we show that the efficiency of LDBSP on DNA derived from LCM neurons is similar to the efficiency obtained in previously published studies using this technique on other cell types. Overall, the method described here provides the user with a more accurate estimation of the DNA methylation status of each target gene in the analyzed cell pools, thereby adding further validity to this approach.

Challenges in Analyzing Functional Epigenetic Data in Perspective of Adolescent Psychiatric Health

The formative period of adolescence plays a crucial role in the development of skills and abilities for adulthood. Adolescents who are affected by mental health conditions are at risk of suicide and social and academic impairments. Gene–environment complementary contributions to the molecular mechanisms involved in psychiatric disorders have emphasized the need to analyze epigenetic marks such as DNA methylation (DNAm) and non-coding RNAs. However, the large and diverse bioinformatic and statistical methods, referring to the confounders of the statistical models, application of multiple-testing adjustment methods, questions regarding the correlation of DNAm across tissues, and sex-dependent differences in results, have raised challenges regarding the interpretation of the results. Based on the example of generalized anxiety disorder (GAD) and depressive disorder (MDD), we shed light on the current knowledge and usage of methodological tools in analyzing epigenetics. Statistical robustness is an essential prerequisite for a better understanding and interpretation of epigenetic modifications and helps to find novel targets for personalized therapeutics in psychiatric diseases.

A pan-tissue DNA methylation atlas enables in silico decomposition of human tissue methylomes at cell-type resolution

Bulk-tissue DNA methylomes represent an average over many different cell types, hampering our understanding of cell-type-specific contributions to disease development. As single-cell methylomics is not scalable to large cohorts of individuals, cost-effective computational solutions are needed, yet current methods are limited to tissues such as blood. Here we leverage the high-resolution nature of tissue-specific single-cell RNA-sequencing datasets to construct a DNA methylation atlas defined for 13 solid tissue types and 40 cell types. We comprehensively validate this atlas in independent bulk and single-nucleus DNA methylation datasets. We demonstrate that it correctly predicts the cell of origin of diverse cancer types and discovers new prognostic associations in olfactory neuroblastoma and stage 2 melanoma. In brain, the atlas predicts a neuronal origin for schizophrenia, with neuron-specific differential DNA methylation enriched for corresponding genome-wide association study risk loci. In summary, the DNA methylation atlas enables the decomposition of 13 different human tissue types at a high cellular resolution, paving the way for an improved interpretation of epigenetic data.

This resource presents an in silico generated DNA methylation atlas that can be used for cell-type deconvolution of human tissues.

Cell-type heterogeneity: Why we should adjust for it in epigenome and biomarker studies

Most studies aiming to identify epigenetic biomarkers do so from complex tissues that are composed of many different cell-types. By definition, these cell-types vary substantially in terms of their epigenetic profiles. This cell-type specific variation among healthy cells is completely independent of the variation associated with disease, yet it dominates the epigenetic variability landscape. While cell-type composition of tissues can change in disease and this may provide accurate and reproducible biomarkers, not adjusting for the underlying cell-type heterogeneity may seriously limit the sensitivity and precision to detect disease-relevant biomarkers or hamper our understanding of such biomarkers. Given that computational and experimental tools for tackling cell-type heterogeneity are available, we here stress that future epigenetic biomarker studies should aim to provide estimates of underlying cell-type fractions for all samples in the study, and to identify biomarkers before and after adjustment for cell-type heterogeneity, in order to obtain a more complete and unbiased picture of the biomarker-landscape. This is critical, not only to improve reproducibility and for the eventual clinical application of such biomarkers, but importantly, to also improve our molecular understanding of disease itself.

Novel epigenetic network biomarkers for early detection of esophageal cancer

BACKGROUND

Early detection of esophageal cancer is critical to improve survival. Whilst studies have identified biomarkers, their interpretation and validity is often confounded by cell-type heterogeneity.

RESULTS

Here we applied systems-epigenomic and cell-type deconvolution algorithms to a discovery set encompassing RNA-Seq and DNA methylation data from esophageal adenocarcinoma (EAC) patients and matched normal-adjacent tissue, in order to identify robust biomarkers, free from the confounding effect posed by cell-type heterogeneity. We identify 12 gene-modules that are epigenetically deregulated in EAC, and are able to validate all 12 modules in 4 independent EAC cohorts. We demonstrate that the epigenetic deregulation is present in the epithelial compartment of EAC-tissue. Using single-cell RNA-Seq data we show that one of these modules, a proto-cadherin module centered around CTNND2, is inactivated in Barrett's Esophagus, a precursor lesion to EAC. By measuring DNA methylation in saliva from EAC cases and controls, we identify a chemokine module centered around CCL20, whose methylation patterns in saliva correlate with EAC status.

CONCLUSIONS

Given our observations that a CCL20 chemokine network is overactivated in EAC tissue and saliva from EAC patients, and that in independent studies CCL20 has been found to be overactivated in EAC tissue infected with the bacterium F. nucleatum, a bacterium that normally inhabits the oral cavity, our results highlight the possibility of using DNAm measurements in saliva as a proxy for changes occurring in the esophageal epithelium. Both the CTNND2/CCL20 modules represent novel promising network biomarkers for EAC that merit further investigation.

DNA Methylation Profiles of Purified Cell Types in Bronchoalveolar Lavage: Applications for Mixed Cell Paediatric Pulmonary Studies

In epigenome-wide association studies analysing DNA methylation from samples containing multiple cell types, it is essential to adjust the analysis for cell type composition. One well established strategy for achieving this is reference-based cell type deconvolution, which relies on knowledge of the DNA methylation profiles of purified constituent cell types. These are then used to estimate the cell type proportions of each sample, which can then be incorporated to adjust the association analysis. Bronchoalveolar lavage is commonly used to sample the lung in clinical practice and contains a mixture of different cell types that can vary in proportion across samples, affecting the overall methylation profile. A current barrier to the use of bronchoalveolar lavage in DNA methylation-based research is the lack of reference DNA methylation profiles for each of the constituent cell types, thus making reference-based cell composition estimation difficult. Herein, we use bronchoalveolar lavage samples collected from children with cystic fibrosis to define DNA methylation profiles for the four most common and clinically relevant cell types: alveolar macrophages, granulocytes, lymphocytes and alveolar epithelial cells. We then demonstrate the use of these methylation profiles in conjunction with an established reference-based methylation deconvolution method to estimate the cell type composition of two different tissue types; a publicly available dataset derived from artificial blood-based cell mixtures and further bronchoalveolar lavage samples. The reference DNA methylation profiles developed in this work can be used for future reference-based cell type composition estimation of bronchoalveolar lavage. This will facilitate the use of this tissue in studies examining the role of DNA methylation in lung health and disease.

Detecting cord blood cell type-specific epigenetic associations with gestational diabetes mellitus and early childhood growth

BACKGROUND

Epigenome-wide association studies (EWAS) have provided opportunities to understand the role of epigenetic mechanisms in development and pathophysiology of many chronic diseases. However, an important limitation of conventional EWAS is that profiles of epigenetic variability are often obtained in samples of mixed cell types. Here, we aim to assess whether changes in cord blood DNA methylation (DNAm) associated with gestational diabetes mellitus (GDM) exposure and early childhood growth markers occur in a cell type-specific manner.

RESULTS

We analyzed 275 cord blood samples collected at delivery from a prospective pre-birth cohort with genome-wide DNAm profiled by the Illumina MethylationEPIC array. We estimated proportions of seven common cell types in each sample using a cord blood-specific DNAm reference panel. Leveraging a recently developed approach named CellDMC, we performed cell type-specific EWAS to identify CpG loci significantly associated with GDM, or 3-year-old body mass index (BMI) z-score. A total of 1410 CpG loci displayed significant cell type-specific differences in methylation level between 23 GDM cases and 252 controls with a false discovery rate < 0.05. Gene Ontology enrichment analysis indicated that LDL transportation emerged from CpG specifically identified from B-cells DNAm analyses and the mitogen-activated protein kinase pathway emerged from CpG specifically identified from natural killer cells DNAm analyses. In addition, we identified four and six loci associated with 3-year-old BMI z-score that were specific to CD8+ T-cells and monocytes, respectively. By performing genome-wide permutation tests, we validated that most of our detected signals had low false positive rates.

CONCLUSION

Compared to conventional EWAS adjusting for the effects of cell type heterogeneity, the proposed approach based on cell type-specific EWAS could provide additional biologically meaningful associations between CpG methylation, prenatal maternal GDM or 3-year-old BMI. With careful validation, these findings may provide new insights into the pathogenesis, programming, and consequences of related childhood metabolic dysregulation. Therefore, we propose that cell type-specific analyses are worth cautious explorations.

Determining subpopulation methylation profiles from bisulfite sequencing data of heterogeneous samples using DXM

Epigenetic changes, such as aberrant DNA methylation, contribute to cancer clonal expansion and disease progression. However, identifying subpopulation-level changes in a heterogeneous sample remains challenging. Thus, we have developed a computational approach, DXM, to deconvolve the methylation profiles of major allelic subpopulations from the bisulfite sequencing data of a heterogeneous sample. DXM does not require prior knowledge of the number of subpopulations or types of cells to expect. We benchmark DXM's performance and demonstrate improvement over existing methods. We further experimentally validate DXM predicted allelic subpopulation-methylation profiles in four Diffuse Large B-Cell Lymphomas (DLBCLs). Lastly, as proof-of-concept, we apply DXM to a cohort of 31 DLBCLs and relate allelic subpopulation methylation profiles to relapse. We thus demonstrate that DXM can robustly find allelic subpopulation methylation profiles that may contribute to disease progression using bisulfite sequencing data of any heterogeneous sample.

Cell Types in Environmental Epigenetic Studies: Biological and Epidemiological Frameworks

PURPOSE OF REVIEW

This article introduces the roles of perinatal DNA methylation in human health and disease, highlights the challenges of tissue and cellular heterogeneity to studying DNA methylation, summarizes approaches to overcome these challenges, and offers recommendations in conducting research in environmental epigenetics.

RECENT FINDINGS

Epigenetic modifications are essential for human development and are labile to environmental influences, especially during gestation. Epigenetic dysregulation is also a hallmark of multiple diseases. Environmental epigenetic studies routinely measure DNA methylation in readily available tissues. However, tissues and cell types exhibit specific epigenetic patterning and heterogeneity between samples complicates epigenetic studies. Failure to account for cell-type heterogeneity limits identification of biological mechanisms and biases study results. Tissue-level epigenetic measures represent a convolution of epigenetic signals from individual cell types. Tissue-specific epigenetics is an evolving field and the use of disease-affected target, surrogate, or multiple tissues has inherent trade-offs and affects inference. Likewise, experimental and bioinformatic approaches to accommodate cell-type heterogeneity have varying assumptions and inherent trade-offs that affect inference. The relationships between exposure, disease, tissue-level DNA methylation, cell type-specific DNA methylation, and cell-type heterogeneity must be carefully considered in study design and analysis. Causal diagrams can inform study design and analytic strategies. Properly addressing cell-type heterogeneity limits sources of potential bias, avoids misinterpretation of study results, and allows investigators to distinguish shifts in cell-type proportions from direct changes to cellular epigenetic programming, both of which provide insights into environmental disease etiology and aid development of novel methods for prevention and treatment.

A cell-type deconvolution meta-analysis of whole blood EWAS reveals lineage-specific smoking-associated DNA methylation changes

Highly reproducible smoking-associated DNA methylation changes in whole blood have been reported by many Epigenome-Wide-Association Studies (EWAS). These epigenetic alterations could have important implications for understanding and predicting the risk of smoking-related diseases. To this end, it is important to establish if these DNA methylation changes happen in all blood cell subtypes or if they are cell-type specific. Here, we apply a cell-type deconvolution algorithm to identify cell-type specific DNA methylation signals in seven large EWAS. We find that most of the highly reproducible smoking-associated hypomethylation signatures are more prominent in the myeloid lineage. A meta-analysis further identifies a myeloid-specific smoking-associated hypermethylation signature enriched for DNase Hypersensitive Sites in acute myeloid leukemia. These results may guide the design of future smoking EWAS and have important implications for our understanding of how smoking affects immune-cell subtypes and how this may influence the risk of smoking related diseases.

Smoking-associated DNA methylation changes in whole blood have been reported by many EWAS. Here, the authors use a cell-type deconvolution algorithm to identify cell-type specific DNA methylation signals in seven EWAS, identifying lineage-specific smoking-associated DNA methylation changes.

A comparison of epigenetic mitotic-like clocks for cancer risk prediction

Background

DNA methylation changes that accrue in the stem cell pool of an adult tissue in line with the cumulative number of cell divisions may contribute to the observed variation in cancer risk among tissues and individuals. Thus, the construction of epigenetic “mitotic” clocks that can measure the lifetime number of stem cell divisions is of paramount interest.

Methods

Building upon a dynamic model of DNA methylation gain in unmethylated CpG-rich regions, we here derive a novel mitotic clock (“epiTOC2”) that can directly estimate the cumulative number of stem cell divisions in a tissue. We compare epiTOC2 to a different mitotic model, based on hypomethylation at solo-WCGW sites (“HypoClock”), in terms of their ability to measure mitotic age of normal adult tissues and predict cancer risk.

Results

Using epiTOC2, we estimate the intrinsic stem cell division rate for different normal tissue types, demonstrating excellent agreement (Pearson correlation = 0.92, R² = 0.85, P = 3e−6) with those derived from experiment. In contrast, HypoClock’s estimates do not (Pearson correlation = 0.30, R² = 0.09, P = 0.29). We validate these results in independent datasets profiling normal adult tissue types. While both epiTOC2 and HypoClock correctly predict an increased mitotic rate in cancer, epiTOC2 is more robust and significantly better at discriminating preneoplastic lesions characterized by chronic inflammation, a major driver of tissue turnover and cancer risk. Our data suggest that DNA methylation loss at solo-WCGWs is significant only when cells are under high replicative stress and that epiTOC2 is a better mitotic age and cancer risk prediction model for normal adult tissues.

Conclusions

These results have profound implications for our understanding of epigenetic clocks and for developing cancer risk prediction or early detection assays. We propose that measurement of DNAm at the 163 epiTOC2 CpGs in adult pre-neoplastic lesions, and potentially in serum cell-free DNA, could provide the basis for building feasible pre-diagnostic or cancer risk assays. epiTOC2 is freely available from 10.5281/zenodo.2632938

Capturing functional epigenomes for insight into metabolic diseases

Background

Metabolic diseases such as obesity are known to be driven by both environmental and genetic factors. Although genome-wide association studies of common variants and their impact on complex traits have provided some biological insight into disease etiology, identified genetic variants have been found to contribute only a small proportion to disease heritability, and to map mainly to non-coding regions of the genome. To link variants to function, association studies of cellular traits, such as epigenetic marks, in disease-relevant tissues are commonly applied.

Scope of the review

We review large-scale efforts to generate genome-wide maps of coordinated epigenetic marks and their utility in complex disease dissection with a focus on DNA methylation. We contrast DNA methylation profiling methods and discuss the advantages of using targeted methods for single-base resolution assessments of methylation levels across tissue-specific regulatory regions to deepen our understanding of contributing factors leading to complex diseases.

Major conclusions

Large-scale assessments of DNA methylation patterns in metabolic disease-linked study cohorts have provided insight into the impact of variable epigenetic variants in disease etiology. In-depth profiling of epigenetic marks at regulatory regions, particularly at tissue-specific elements, will be key to dissect the genetic and environmental components contributing to metabolic disease onset and progression.

•

Changes in epigenetic marks have been linked to metabolic disease phenotypes.

•

Disease-linked sites of variable DNA methylation status are enriched in distal regulatory regions of disease-linked tissues.

•

Distal regulatory elements remain underrepresented in popular array-based methylation profiling technologies.

•

Novel next-generation capture methods provide cost-effective solutions to assess the impact of DNA methylation in metabolic diseases specifically at regulatory elements.

•

Improvements in methodologies to account for tissue heterogeneity and causality will be crucial in future epigenome-wide association studies.

methylCC: technology-independent estimation of cell type composition using differentially methylated regions

A major challenge in the analysis of DNA methylation (DNAm) data is variability introduced from intra-sample cellular heterogeneity, such as whole blood which is a convolution of DNAm profiles across a unique cell type. When this source of variability is confounded with an outcome of interest, if unaccounted for, false positives ensue. Current methods to estimate the cell type proportions in whole blood DNAm samples are only appropriate for one technology and lead to technology-specific biases if applied to data generated from other technologies. Here, we propose the technology-independent alternative: methylCC, which is available at https://github.com/stephaniehicks/methylCC.

TOAST: improving reference-free cell composition estimation by cross-cell type differential analysis

In the analysis of high-throughput data from complex samples, cell composition is an important factor that needs to be accounted for. Except for a limited number of tissues with known pure cell type profiles, a majority of genomics and epigenetics data relies on the "reference-free deconvolution" methods to estimate cell composition. We develop a novel computational method to improve reference-free deconvolution, which iteratively searches for cell type-specific features and performs composition estimation. Simulation studies and applications to six real datasets including both DNA methylation and gene expression data demonstrate favorable performance of the proposed method. TOAST is available at https://bioconductor.org/packages/TOAST .

Detection of cell-type-specific risk-CpG sites in epigenome-wide association studies

In epigenome-wide association studies, the measured signals for each sample are a mixture of methylation profiles from different cell types. Current approaches to the association detection claim whether a cytosine-phosphate-guanine (CpG) site is associated with the phenotype or not at aggregate level and can suffer from low statistical power. Here, we propose a statistical method, HIgh REsolution (HIRE), which not only improves the power of association detection at aggregate level as compared to the existing methods but also enables the detection of risk-CpG sites for individual cell types.

Cellular heterogeneity is one of the major confounding factors in EWAS studies. Here the authors present a statistical method, HIgh REsolution (HIRE), which enables the detection of risk-CpG sites for individual cell types.

CONFINED: distinguishing biological from technical sources of variation by leveraging multiple methylation datasets

Methylation datasets are affected by innumerable sources of variability, both biological (cell-type composition, genetics) and technical (batch effects). Here, we propose a reference-free method based on sparse canonical correlation analysis to separate the biological from technical sources of variability. We show through simulations and real data that our method, CONFINED, is not only more accurate than the state-of-the-art reference-free methods for capturing known, replicable biological variability, but it is also considerably more robust to dataset-specific technical variability than previous approaches. CONFINED is available as an R package as detailed at https://github.com/cozygene/CONFINED.

Epigenetic marks of prenatal air pollution exposure found in multiple tissues relevant for child health

BACKGROUND

Prenatal air pollution exposure has been linked to many adverse health conditions in the offspring. However, little is known about the mechanisms underlying these associations. Epigenetics may be one plausible biologic link. Here, we sought to identify site-specific and global DNA methylation (DNAm) changes, in developmentally relevant tissues, associated with prenatal exposure to nitrogen dioxide (NO₂) and ozone (O₃). Additionally, we assessed whether sex-specific changes in methylation exist and whether DNAm changes are consistently observed across tissues.

METHODS

Genome-scale DNAm measurements were obtained using the Infinium HumanMethylation450k platform for 133 placenta and 175 cord blood specimens from Early Autism Risk Longitudinal Investigation (EARLI) neonates. Ambient NO₂ and O₃ exposure levels were based on prenatal address locations of EARLI mothers and the Environmental Protection Agency's AirNOW monitoring network using inverse distance weighting. We computed sample-level aggregate methylation measures for each of 5 types of genomic regions including genome-wide, open sea, shelf, shore, and island regions. Linear regression was performed for each genomic region; per-sample aggregate methylation measures were modeled as a function of quantitative exposure level with covariate adjustment. In addition, bumphunting was performed to identify differentially methylated regions (DMRs) associated with prenatal O₃ and NO₂ exposures in each tissue and by sex, with adjustment for technical and biological sources of variation.

RESULTS

We identified global and locus-specific changes in DNA methylation related to prenatal exposure to NO₂ and O₃ in 2 developmentally relevant tissues. Neonates with increased prenatal O₃ exposure had lower aggregate levels of DNAm at CpGs located in open sea and shelf regions of the genome. We identified 6 DMRs associated with prenatal NO₂ exposure, including 3 sex-specific. An additional 3 sex-specific DMRs were associated with prenatal O₃ exposure levels. DMRs initially detected in cord blood samples (n = 4) showed consistent exposure-related changes in DNAm in placenta. However, the DMRs initially detected in placenta (n = 5) did not show DNAm differences in cord blood and, thus, they appear to be tissue-specific.

CONCLUSIONS

We observed global, locus, and sex-specific methylation changes associated with prenatal NO₂ and O₃ exposures. Our findings support DNAm is a biologic target of prenatal air pollutant exposures and highlight epigenetic involvement in sex-specific differential susceptibility to environmental exposure effects in 2 developmentally relevant tissues.

Genome-wide DNA methylation profiling shows a distinct epigenetic signature associated with lung macrophages in cystic fibrosis

Background

Lung macrophages are major participants in the pulmonary innate immune response. In the cystic fibrosis (CF) lung, the inability of lung macrophages to successfully regulate the exaggerated inflammatory response suggests dysfunctional innate immune cell function. In this study, we aim to gain insight into innate immune cell dysfunction in CF by investigating alterations in DNA methylation in bronchoalveolar lavage (BAL) cells, composed primarily of lung macrophages of CF subjects compared with healthy controls. All analyses were performed using primary alveolar macrophages from human subjects collected via bronchoalveolar lavage. Epigenome-wide DNA methylation was examined via Illumina MethylationEPIC (850 K) array. Targeted next-generation bisulfite sequencing was used to validate selected differentially methylated CpGs. Methylation-based sample classification was performed using the recursively partitioned mixture model (RPMM) and was tested against sample case-control status. Differentially methylated loci were identified by fitting linear models with adjustment of age, sex, estimated cell type proportions, and repeat measurement.

Results

RPMM class membership was significantly associated with the CF disease status (P = 0.026). One hundred nine CpG loci were differentially methylated in CF BAL cells (all FDR ≤ 0.1). The majority of differentially methylated loci in CF were hypo-methylated and found within non-promoter CpG islands as well as in putative enhancer regions and DNase hyper-sensitive regions.

Conclusions

These results support a hypothesis that epigenetic changes, specifically DNA methylation at a multitude of gene loci in lung macrophages, may participate, at least in part, in driving dysfunctional innate immune cells in the CF lung.

Electronic supplementary material

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

Identification of differentially methylated cell types in epigenome-wide association studies

An outstanding challenge of epigenome-wide association studies (EWASs) performed in complex tissues is the identification of the specific cell type(s) responsible for the observed differential DNA methylation. Here we present a statistical algorithm called CellDMC ( https://github.com/sjczheng/EpiDISH ), which can identify differentially methylated positions and the specific cell type(s) driving the differential methylation. We validated CellDMC on in silico mixtures of DNA methylation data generated with different technologies, as well as on real mixtures from epigenome-wide association and cancer epigenome studies. CellDMC achieved over 90% sensitivity and specificity in scenarios where current state-of-the-art methods did not identify differential methylation. By applying CellDMC to an EWAS performed in buccal swabs, we identified smoking-associated differentially methylated positions occurring in the epithelial compartment, which we validated in smoking-related lung cancer. CellDMC may be useful in the identification of causal DNA-methylation alterations in disease.

Cumulative lifetime maternal stress and epigenome-wide placental DNA methylation in the PRISM cohort

Evolving evidence links maternal stress exposure to changes in placental DNA methylation of specific genes regulating placental function that may have implications for the programming of a host of chronic disorders. Few studies have implemented an epigenome-wide approach. Using the Infinium HumanMethylation450 BeadChip (450K), we investigated epigenome-wide placental DNA methylation in relation to maternal experiences of traumatic and non-traumatic stressors over her lifetime assessed using the Life Stressor Checklist-Revised (LSC-R) survey (n = 207). We found differential DNA methylation at epigenome-wide statistical significance (FDR = 0.05) for 112 CpGs. Additionally, we observed three clusters that exhibited differential methylation in response to high maternal lifetime stress. Enrichment analyses, conducted at an FDR = 0.20, revealed lysine degradation to be the most significant pathway associated with maternal lifetimes stress exposure. Targeted enrichment analyses of the three largest clusters of probes, identified using the gap statistic, were enriched for genes associated with endocytosis (i.e., SMAP1, ANKFY1), tight junctions (i.e., EPB41L4B), and metabolic pathways (i.e., INPP5E, EEF1B2). These pathways, also identified in the top 10 KEGG pathways associated with maternal lifetime stress exposure, play important roles in multiple physiological functions necessary for proper fetal development. Further, two genes were identified to exhibit multiple probes associated with maternal lifetime stress (i.e., ANKFY1, TM6SF1). The methylation status of the probes belonging to each cluster and/or genes exhibiting multiple hits, may play a role in the pathogenesis of adverse health outcomes in children born to mothers with increased lifetime stress exposure.

DNA methylation analysis on purified neurons and glia dissects age and Alzheimer's disease-specific changes in the human cortex

Background

Epigenome-wide association studies (EWAS) based on human brain samples allow a deep and direct understanding of epigenetic dysregulation in Alzheimer’s disease (AD). However, strong variation of cell-type proportions across brain tissue samples represents a significant source of data noise. Here, we report the first EWAS based on sorted neuronal and non-neuronal (mostly glia) nuclei from postmortem human brain tissues.

Results

We show that cell sorting strongly enhances the robust detection of disease-related DNA methylation changes even in a relatively small cohort. We identify numerous genes with cell-type-specific methylation signatures and document differential methylation dynamics associated with aging specifically in neurons such as CLU, SYNJ2 and NCOR2 or in glia RAI1,CXXC5 and INPP5A. Further, we found neuron or glia-specific associations with AD Braak stage progression at genes such as MCF2L, ANK1, MAP2, LRRC8B, STK32C and S100B. A comparison of our study with previous tissue-based EWAS validates multiple AD-associated DNA methylation signals and additionally specifies their origin to neuron, e.g., HOXA3 or glia (ANK1). In a meta-analysis, we reveal two novel previously unrecognized methylation changes at the key AD risk genes APP and ADAM17.

Conclusions

Our data highlight the complex interplay between disease, age and cell-type-specific methylation changes in AD risk genes thus offering new perspectives for the validation and interpretation of large EWAS results.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0211-3) contains supplementary material, which is available to authorized users.

A novel cell-type deconvolution algorithm reveals substantial contamination by immune cells in saliva, buccal and cervix

Aim: An outstanding challenge in epigenome studies is the estimation of cell-type proportions in complex epithelial tissues. Materials & methods: Here, we construct and validate a DNA methylation reference and algorithm for complex tissues that contain epithelial, immune and nonimmune stromal cells. Results: Using this reference, we show that easily accessible tissues such as saliva, buccal and cervix exhibit substantial variation in immune cell (IC) contamination. We further validate our reference in the context of oral cancer, where it correctly predicts an increased IC infiltration in cancer but suppressed in patients with highest smoking exposure. Finally, our method can improve the specificity of differentially methylated CpG calls in epithelial cancer. Conclusion: The degree and variation of IC contamination in complex epithelial tissues is substantial. We provide a valuable resource and tool for assessing the epithelial purity and IC contamination of samples and for identifying differential methylation in such complex tissues.

Distinct chromatin signatures of DNA hypomethylation in aging and cancer

Cancer is an aging‐associated disease, but the underlying molecular links between these processes are still largely unknown. Gene promoters that become hypermethylated in aging and cancer share a common chromatin signature in ES cells. In addition, there is also global DNA hypomethylation in both processes. However, the similarity of the regions where this loss of DNA methylation occurs is currently not well characterized, and it is unknown if such regions also share a common chromatin signature in aging and cancer. To address this issue, we analyzed TCGA DNA methylation data from a total of 2,311 samples, including control and cancer cases from patients with breast, kidney, thyroid, skin, brain, and lung tumors and healthy blood, and integrated the results with histone, chromatin state, and transcription factor binding site data from the NIH Roadmap Epigenomics and ENCODE projects. We identified 98,857 CpG sites differentially methylated in aging and 286,746 in cancer. Hyper‐ and hypomethylated changes in both processes each had a similar genomic distribution across tissues and displayed tissue‐independent alterations. The identified hypermethylated regions in aging and cancer shared a similar bivalent chromatin signature. In contrast, hypomethylated DNA sequences occurred in very different chromatin contexts. DNA hypomethylated sequences were enriched at genomic regions marked with the activating histone posttranslational modification H3K4me1 in aging, while in cancer, loss of DNA methylation was primarily associated with the repressive H3K9me3 mark. Our results suggest that the role of DNA methylation as a molecular link between aging and cancer is more complex than previously thought.

An optimized library for reference-based deconvolution of whole-blood biospecimens assayed using the Illumina HumanMethylationEPIC BeadArray

Genome-wide methylation arrays are powerful tools for assessing cell composition of complex mixtures. We compare three approaches to select reference libraries for deconvoluting neutrophil, monocyte, B-lymphocyte, natural killer, and CD4+ and CD8+ T-cell fractions based on blood-derived DNA methylation signatures assayed using the Illumina HumanMethylationEPIC array. The IDOL algorithm identifies a library of 450 CpGs, resulting in an average R2 = 99.2 across cell types when applied to EPIC methylation data collected on artificial mixtures constructed from the above cell types. Of the 450 CpGs, 69% are unique to EPIC. This library has the potential to reduce unintended technical differences across array platforms.

Epigenome-wide association studies for systemic autoimmune diseases: The road behind and the road ahead

Epigenetics is known to be an important mechanism in the pathogenesis of autoimmune diseases. Epigenetic variations can act as integrators of environmental and genetic exposures and propagate activated states in immune cells. Studying epigenetic alterations by means of genome-wide approaches promises to unravel novel molecular mechanisms related to disease etiology, disease progression, clinical manifestations and treatment responses. This paper reviews what we have learned in the last five years from epigenome-wide studies for three systemic autoimmune diseases, namely systemic lupus erythematosus, primary Sjögren's syndrome, and rheumatoid arthritis. We examine the degree of epigenetic sharing between different diseases and the possible mediating role of epigenetic associations in genetic and environmental risks. Finally, we also shed light into the use of epigenetic markers towards a better precision medicine regarding disease prediction, prevention and personalized treatment in systemic autoimmunity.

Analysis of DNA modifications in aging research

As geroscience research extends into the role of epigenetics in aging and age-related disease, researchers are being confronted with unfamiliar molecular techniques and data analysis methods that can be difficult to integrate into their work. In this review, we focus on the analysis of DNA modifications, namely cytosine methylation and hydroxymethylation, through next-generation sequencing methods. While older techniques for modification analysis performed relative quantitation across regions of the genome or examined average genome levels, these analyses lack the desired specificity, rigor, and genomic coverage to firmly establish the nature of genomic methylation patterns and their response to aging. With recent methodological advances, such as whole genome bisulfite sequencing (WGBS), bisulfite oligonucleotide capture sequencing (BOCS), and bisulfite amplicon sequencing (BSAS), cytosine modifications can now be readily analyzed with base-specific, absolute quantitation at both cytosine-guanine dinucleotide (CG) and non-CG sites throughout the genome or within specific regions of interest by next-generation sequencing. Additional advances, such as oxidative bisulfite conversion to differentiate methylation from hydroxymethylation and analysis of limited input/single-cells, have great promise for continuing to expand epigenomic capabilities. This review provides a background on DNA modifications, the current state-of-the-art for sequencing methods, bioinformatics tools for converting these large data sets into biological insights, and perspectives on future directions for the field.

Systems-epigenomics inference of transcription factor activity implicates aryl-hydrocarbon-receptor inactivation as a key event in lung cancer development

Background

Diverse molecular alterations associated with smoking in normal and precursor lung cancer cells have been reported, yet their role in lung cancer etiology remains unclear. A prominent example is hypomethylation of the aryl hydrocarbon-receptor repressor (AHRR) locus, which is observed in blood and squamous epithelial cells of smokers, but not in lung cancer.

Results

Using a novel systems-epigenomics algorithm, called SEPIRA, which leverages the power of a large RNA-sequencing expression compendium to infer regulatory activity from messenger RNA expression or DNA methylation (DNAm) profiles, we infer the landscape of binding activity of lung-specific transcription factors (TFs) in lung carcinogenesis. We show that lung-specific TFs become preferentially inactivated in lung cancer and precursor lung cancer lesions and further demonstrate that these results can be derived using only DNAm data. We identify subsets of TFs which become inactivated in precursor cells. Among these regulatory factors, we identify AHR, the aryl hydrocarbon-receptor which controls a healthy immune response in the lung epithelium and whose repressor, AHRR, has recently been implicated in smoking-mediated lung cancer. In addition, we identify FOXJ1, a TF which promotes growth of airway cilia and effective clearance of the lung airway epithelium from carcinogens.

Conclusions

We identify TFs, such as AHR, which become inactivated in the earliest stages of lung cancer and which, unlike AHRR hypomethylation, are also inactivated in lung cancer itself. The novel systems-epigenomics algorithm SEPIRA will be useful to the wider epigenome-wide association study community as a means of inferring regulatory activity.

Electronic supplementary material

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

Statistical and integrative system-level analysis of DNA methylation data

Epigenetics plays a key role in cellular development and function. Alterations to the epigenome are thought to capture and mediate the effects of genetic and environmental risk factors on complex disease. Currently, DNA methylation is the only epigenetic mark that can be measured reliably and genome-wide in large numbers of samples. This Review discusses some of the key statistical challenges and algorithms associated with drawing inferences from DNA methylation data, including cell-type heterogeneity, feature selection, reverse causation and system-level analyses that require integration with other data types such as gene expression, genotype, transcription factor binding and other epigenetic information.

Epigenetics in pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children. ALL arises from the malignant transformation of progenitor B- and T-cells in the bone marrow into leukemic cells, but the mechanisms underlying this transformation are not well understood. Recent technical advances and decreasing costs of methods for high-throughput DNA sequencing and SNP genotyping have stimulated systematic studies of the epigenetic changes in leukemic cells from pediatric ALL patients. The results emerging from these studies are increasing our understanding of the epigenetic component of leukemogenesis and have demonstrated the potential of DNA methylation as a biomarker for lineage and subtype classification, prognostication, and disease progression in ALL. In this review, we provide a concise examination of the epigenetic studies in ALL, with a focus on DNA methylation and mutations perturbing genes involved in chromatin modification, and discuss the future role of epigenetic analyses in research and clinical management of ALL.

Maximizing ecological and evolutionary insight in bisulfite sequencing data sets

Genome-scale bisulfite sequencing approaches have opened the door to ecological and evolutionary studies of DNA methylation in many organisms. These approaches can be powerful. However, they introduce new methodological and statistical considerations, some of which are particularly relevant to non-model systems. Here, we highlight how these considerations influence a study's power to link methylation variation with a predictor variable of interest. Relative to current practice, we argue that sample sizes will need to increase to provide robust insights. We also provide recommendations for overcoming common challenges and an R Shiny app to aid in study design.

Is cellular heterogeneity merely a confounder to be removed from epigenome-wide association studies?

Excitement about DNA methylation biomarkers has been tempered by a growing appreciation of the complex causal relations with cell fate. Intersample differences in DNA methylation can be partitioned into those that are independent of cellular heterogeneity and those that are caused by differential mixtures of cell types. Generally, the field has assumed that the former are more likely to be causative of disease. The latter has been considered a likely consequence of disease and a confounder to be removed. We argue that the conceptual separation of these signals is artificial and not necessarily informative about causation. DNA methylation is a very sensitive measure of cell fate mix and therefore reveals much about underlying disease etiology including aspects of causation.