Estimands in epigenome-wide association studies

Transcriptomics and epigenetic data integration learning module on Google Cloud

Multi-omics (genomics, transcriptomics, epigenomics, proteomics, metabolomics, etc.) research approaches are vital for understanding the hierarchical complexity of human biology and have proven to be extremely valuable in cancer research and precision medicine. Emerging scientific advances in recent years have made high-throughput genome-wide sequencing a central focus in molecular research by allowing for the collective analysis of various kinds of molecular biological data from different types of specimens in a single tissue or even at the level of a single cell. Additionally, with the help of improved computational resources and data mining, researchers are able to integrate data from different multi-omics regimes to identify new prognostic, diagnostic, or predictive biomarkers, uncover novel therapeutic targets, and develop more personalized treatment protocols for patients. For the research community to parse the scientifically and clinically meaningful information out of all the biological data being generated each day more efficiently with less wasted resources, being familiar with and comfortable using advanced analytical tools, such as Google Cloud Platform becomes imperative. This project is an interdisciplinary, cross-organizational effort to provide a guided learning module for integrating transcriptomics and epigenetics data analysis protocols into a comprehensive analysis pipeline for users to implement in their own work, utilizing the cloud computing infrastructure on Google Cloud. The learning module consists of three submodules that guide the user through tutorial examples that illustrate the analysis of RNA-sequence and Reduced-Representation Bisulfite Sequencing data. The examples are in the form of breast cancer case studies, and the data sets were procured from the public repository Gene Expression Omnibus. The first submodule is devoted to transcriptomics analysis with the RNA sequencing data, the second submodule focuses on epigenetics analysis using the DNA methylation data, and the third submodule integrates the two methods for a deeper biological understanding. The modules begin with data collection and preprocessing, with further downstream analysis performed in a Vertex AI Jupyter notebook instance with an R kernel. Analysis results are returned to Google Cloud buckets for storage and visualization, removing the computational strain from local resources. The final product is a start-to-finish tutorial for the researchers with limited experience in multi-omics to integrate transcriptomics and epigenetics data analysis into a comprehensive pipeline to perform their own biological research.This manuscript describes the development of a resource module that is part of a learning platform named ``NIGMS Sandbox for Cloud-based Learning'' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox [16] at the beginning of this Supplement. This module delivers learning materials on the analysis of bulk and single-cell ATAC-seq data in an interactive format that uses appropriate cloud resources for data access and analyses.

HIGHLIGHTS

DNA-methylation signature accurately differentiates pancreatic cancer from chronic pancreatitis in tissue and plasma

OBJECTIVE

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy. Differentiation from chronic pancreatitis (CP) is currently inaccurate in about one-third of cases. Misdiagnoses in both directions, however, have severe consequences for patients. We set out to identify molecular markers for a clear distinction between PDAC and CP.

DESIGN

Genome-wide variations of DNA-methylation, messenger RNA and microRNA level as well as combinations thereof were analysed in 345 tissue samples for marker identification. To improve diagnostic performance, we established a random-forest machine-learning approach. Results were validated on another 48 samples and further corroborated in 16 liquid biopsy samples.

RESULTS

Machine-learning succeeded in defining markers to differentiate between patients with PDAC and CP, while low-dimensional embedding and cluster analysis failed to do so. DNA-methylation yielded the best diagnostic accuracy by far, dwarfing the importance of transcript levels. Identified changes were confirmed with data taken from public repositories and validated in independent sample sets. A signature of six DNA-methylation sites in a CpG-island of the protein kinase C beta type gene achieved a validated diagnostic accuracy of 100% in tissue and in circulating free DNA isolated from patient plasma.

CONCLUSION

The success of machine-learning to identify an effective marker signature documents the power of this approach. The high diagnostic accuracy of discriminating PDAC from CP could have tremendous consequences for treatment success, once the result from still a limited number of liquid biopsy samples would be confirmed in a larger cohort of patients with suspected pancreatic cancer.

Validation of an Integrated Genetic-Epigenetic Test for the Assessment of Coronary Heart Disease

BACKGROUND

Coronary heart disease (CHD) is the leading cause of death in the world. Unfortunately, many of the key diagnostic tools for CHD are insensitive, invasive, and costly; require significant specialized infrastructure investments; and do not provide information to guide postdiagnosis therapy. In prior work using data from the Framingham Heart Study, we provided in silico evidence that integrated genetic-epigenetic tools may provide a new avenue for assessing CHD.

METHODS AND RESULTS

In this communication, we use an improved machine learning approach and data from 2 additional cohorts, totaling 449 cases and 2067 controls, to develop a better model for ascertaining symptomatic CHD. Using the DNA from the 2 new cohorts, we translate and validate the in silico findings into an artificial intelligence-guided, clinically implementable method that uses input from 6 methylation-sensitive digital polymerase chain reaction and 10 genotyping assays. Using this method, the overall average area under the curve, sensitivity, and specificity in the 3 test cohorts is 82%, 79%, and 76%, respectively. Analysis of targeted cytosine-phospho-guanine loci shows that they map to key risk pathways involved in atherosclerosis that suggest specific therapeutic approaches.

CONCLUSIONS

We conclude that this scalable integrated genetic-epigenetic approach is useful for the diagnosis of symptomatic CHD, performs favorably as compared with many existing methods, and may provide personalized insight to CHD therapy. Furthermore, given the dynamic nature of DNA methylation and the ease of methylation-sensitive digital polymerase chain reaction methodologies, these findings may pave a pathway for precision epigenetic approaches for monitoring CHD treatment response.

Epigenetic biomarkers for smoking cessation

Cigarette smoking has been associated with epigenetic alterations that may be reversible upon cessation. As the most-studied epigenetic modification, DNA methylation is strongly associated with smoking exposure, providing a potential mechanism that links smoking to adverse health outcomes. Here, we reviewed the reversibility of DNA methylation in accessible peripheral tissues, mainly blood, in relation to cigarette smoking cessation and the utility of DNA methylation as a biomarker signature to differentiate current, former, and never smokers and to quantify time since cessation. We summarized thousands of differentially methylated Cytosine-Guanine (CpG) dinucleotides and regions associated with smoking cessation from candidate gene and epigenome-wide association studies, as well as the prediction accuracy of the multi-CpG predictors for smoking status. Overall, there is robust evidence for DNA methylation signature of cigarette smoking cessation. However, there are still gaps to fill, including (1) cell-type heterogeneity in measuring blood DNA methylation; (2) underrepresentation of non-European ancestry populations; (3) limited longitudinal data to quantitatively measure DNA methylation after smoking cessation over time; and (4) limited data to study the impact of smoking cessation on other epigenetic features, noncoding RNAs, and histone modifications. Epigenetic machinery provides promising biomarkers that can improve success in smoking cessation in the clinical setting. To achieve this goal, larger and more-diverse samples with longitudinal measures of a broader spectrum of epigenetic marks will be essential to developing a robust DNA methylation biomarker assay, followed by meeting validation requirements for the assay before being implemented as a clinically useful tool.

DNA Methylation Analysis in Monozygotic Twins Discordant for ALS in Blood Cells

ALS is a fatal motor neuron disease that displays a broad variety of phenotypes ranging from early fatal courses to slowly progressing and rather benign courses. Such divergence can also be seen in genetic ALS cases with varying phenotypes bearing specific mutations, suggesting epigenetic mechanisms like DNA methylation act as disease modifiers. However, the epigenotype dictated by, in addition to other mechanisms, DNA methylation is also strongly influenced by the individual's genotype. Hence, we performed a DNA methylation study using EPIC arrays on 7 monozygotic (MZ) twin pairs discordant for ALS in whole blood, which serves as an ideal model for eliminating the effects of the genetic-epigenetic interplay to a large extent. We found one CpG site showing intra-pair hypermethylation in the affected co-twins, which maps to the Glutamate Ionotropic Receptor Kainate Type Subunit 1 gene (GRIK1). Additionally, we found 4 DMPs which were subsequently confirmed using 2 different statistical approaches. Differentially methylated regions or blocks could not be detected within the scope of this work. In conclusion, we revealed that despite a low sample size, monozygotic twin studies discordant for the disease can bring new insights into epigenetic processes in ALS, pointing to new target loci for further investigations.

Maternal blood pressure associates with placental DNA methylation both directly and through alterations in cell-type composition

BACKGROUND

Maternal blood pressure levels reflect cardiovascular adaptation to pregnancy and proper maternal-fetal exchanges through the placenta and are very sensitive to numerous environmental stressors. Maternal hypertension during pregnancy has been associated with impaired placental functions and with an increased risk for children to suffer from cardiovascular and respiratory diseases later on. Investigating changes in placental DNA methylation levels and cell-type composition in association with maternal blood pressure could help elucidate its relationships with placental and fetal development.

METHODS

Taking advantage of a large cohort of 666 participants, we investigated the association between epigenome-wide DNA methylation patterns in the placenta, measured using the Infinium HumanMethylation450 BeadChip, placental cell-type composition, estimated in silico, and repeated measurements of maternal steady and pulsatile blood pressure indicators during pregnancy.

RESULTS

At the site-specific level, no significant association was found between maternal blood pressure and DNA methylation levels after correction for multiple testing (false discovery rate < 0.05), but 5 out of 24 previously found CpG associations were replicated (p-value < 0.05). At the regional level, our analyses highlighted 64 differentially methylated regions significantly associated with at least one blood pressure component, including 35 regions associated with mean arterial pressure levels during late pregnancy. These regions were found enriched for genes implicated in lung development and diseases. Further mediation analyses show that a significant part of the association between steady blood pressure-but not pulsatile pressure-and placental methylation can be explained by alterations in placental cell-type composition. In particular, elevated blood pressure levels are associated with a decrease in the ratio between mesenchymal stromal cells and syncytiotrophoblasts, even in the absence of preeclampsia.

CONCLUSIONS

This study provides the first evidence that the association between maternal steady blood pressure during pregnancy and placental DNA methylation is both direct and partly explained by changes in cell-type composition. These results could hint at molecular mechanisms linking maternal hypertension to lung development and early origins of childhood respiratory problems and at the importance of controlling maternal blood pressure during pregnancy.

Digital methylation assessments of alcohol and cigarette consumption account for common variance in accelerated epigenetic ageing

Smoking and Heavy Alcohol Consumption (HAC) are established risk factors for myriad complex disorders of ageing. Yet many prior studies of Epigenetic Ageing (EA) have shown only modest effects of smoking and drinking on accelerated ageing. One potential reason for this conundrum might be the reliance of some prior EA studies on self-reported substance use, which may be unreliable in many samples. To test whether novel, non-self-reported indices would show a stronger association of smoking and HAC to EA, we used methylation sensitive digital PCR (MSdPCR) and data from 437 African American subjects from Wave 7 of the Family and Community Health Study Offspring Cohort to examine the effects of subjective and objective measures of smoking and HAC on 7 indices of EA. Because of limited overall correlations between the various EA indices, we examined patterns of association separately for each index. Consistent with expectations, MSdPCR assessments of smoking and HAC, but not self-reported alcohol consumption, were strongly correlated with accelerated EA. MSdPCR assessments of smoking and HAC accounted for 57% of GrimAge acceleration and the shared variance in GrimAge and DunedinPOAM accelerated EA. We conclude that MSdPCR assessments of smoking and HAC are valuable tools for understanding EA, represent directly targetable conditions for the prevention of premature ageing, and substantially improve upon self-reported assessment of smoking and HAC.

The relationship of smoking to cg05575921 methylation in blood and saliva DNA samples from several studies

Numerous studies have shown that cg05575921 methylation decreases in response to smoking. However, secondary to methodological issues, the magnitude and dose dependency of that response is as of yet unclear. This lack of certainty is a barrier to the use of DNA methylation clinically to assess and monitor smoking status. To better define this relationship, we conducted a joint analysis of methylation sensitive PCR digital (MSdPCR) assessments of cg05575921 methylation in whole blood and/or saliva DNA to smoking using samples from 421 smokers and 423 biochemically confirmed non-smokers from 4 previously published studies. We found that cg05575921 methylation manifested a curvilinear dose dependent decrease in response to increasing cigarette consumption. In whole blood DNA, the Receiver Operating Characteristic (ROC) Area Under the Curve (AUC) of cg05575921 methylation for predicting daily smoking status was 0.98. In saliva DNA, the gross AUC was 0.91 with correction for cellular heterogeneity improving the AUC to 0.94. Methylation status was significantly associated with the Fagerstrom Test for Nicotine Dependence score, but with significant sampling heterogeneity. We conclude that MSdPCR assessments of cg05575921 methylation are a potentially powerful, clinically implementable tool for the assessment and management of smoking.

Epigenetic adaptations of the masticatory mucosa to periodontal inflammation

Background

In mucosal barrier interfaces, flexible responses of gene expression to long-term environmental changes allow adaptation and fine-tuning for the balance of host defense and uncontrolled not-resolving inflammation. Epigenetic modifications of the chromatin confer plasticity to the genetic information and give insight into how tissues use the genetic information to adapt to environmental factors. The oral mucosa is particularly exposed to environmental stressors such as a variable microbiota. Likewise, persistent oral inflammation is the most important intrinsic risk factor for the oral inflammatory disease periodontitis and has strong potential to alter DNA-methylation patterns. The aim of the current study was to identify epigenetic changes of the oral masticatory mucosa in response to long-term inflammation that resulted in periodontitis.

Methods and results

Genome-wide CpG methylation of both inflamed and clinically uninflamed solid gingival tissue biopsies of 60 periodontitis cases was analyzed using the Infinium MethylationEPIC BeadChip. We validated and performed cell-type deconvolution for infiltrated immune cells using the EpiDish algorithm. Effect sizes of DMPs in gingival epithelial and fibroblast cells were estimated and adjusted for confounding factors using our recently developed “intercept-method”. In the current EWAS, we identified various genes that showed significantly different methylation between periodontitis-inflamed and uninflamed oral mucosa in periodontitis patients. The strongest differences were observed for genes with roles in wound healing (ROBO2, PTP4A3), cell adhesion (LPXN) and innate immune response (CCL26, DNAJC1, BPI). Enrichment analyses implied a role of epigenetic changes for vesicle trafficking gene sets.

Conclusions

Our results imply specific adaptations of the oral mucosa to a persistent inflammatory environment that involve wound repair, barrier integrity, and innate immune defense.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01190-7.

External validation of integrated genetic-epigenetic biomarkers for predicting incident coronary heart disease

Aim: The Framingham Risk Score (FRS) and atherosclerotic cardiovascular disease (ASCVD) Pooled Cohort Equation (PCE) for predicting risk for incident coronary heart disease (CHD) work poorly. To improve risk stratification for CHD, we developed a novel integrated genetic-epigenetic tool. Materials & methods: Using machine learning techniques and datasets from the Framingham Heart Study (FHS) and Intermountain Healthcare (IM), we developed and validated an integrated genetic-epigenetic model for predicting 3-year incident CHD. Results: Our approach was more sensitive than FRS and PCE and had high generalizability across cohorts. It performed with sensitivity/specificity of 79/75% in the FHS test set and 75/72% in the IM set. The sensitivity/specificity was 15/93% in FHS and 31/89% in IM for FRS, and sensitivity/specificity was 41/74% in FHS and 69/55% in IM for PCE. Conclusion: The use of our tool in a clinical setting could better identify patients at high risk for a heart attack.