Mediation analysis of alcohol consumption, DNA methylation, and epithelial ovarian cancer

Methods for large-scale single mediator hypothesis testing: Possible choices and comparisons

Mediation hypothesis testing for a large number of mediators is challenging due to the composite structure of the null hypothesis,H 0 : α β = 0 ${H}_{0}:\alpha \beta =0$ (α $\alpha $ : effect of the exposure on the mediator after adjusting for confounders;β $\beta $ : effect of the mediator on the outcome after adjusting for exposure and confounders). In this paper, we reviewed three classes of methods for large-scale one at a time mediation hypothesis testing. These methods are commonly used for continuous outcomes and continuous mediators assuming there is no exposure-mediator interaction so that the productα β $\alpha \beta $ has a causal interpretation as the indirect effect. The first class of methods ignores the impact of different structures under the composite null hypothesis, namely, (1)α = 0 , β ≠ 0 $\alpha =0,\beta \ne 0$ ; (2)α ≠ 0 , β = 0 $\alpha \ne 0,\beta =0$ ; and (3)α = β = 0 $\alpha =\beta =0$ . The second class of methods weights the reference distribution under each case of the null to form a mixture reference distribution. The third class constructs a composite test statistic using the three p values obtained under each case of the null so that the reference distribution of the composite statistic is approximatelyU ( 0 , 1 ) $U(0,1)$ . In addition to these existing methods, we developed the Sobel-comp method belonging to the second class, which uses a corrected mixture reference distribution for Sobel's test statistic. We performed extensive simulation studies to compare all six methods belonging to these three classes in terms of the false positive rates (FPRs) under the null hypothesis and the true positive rates under the alternative hypothesis. We found that the second class of methods which uses a mixture reference distribution could best maintain the FPRs at the nominal level under the null hypothesis and had the greatest true positive rates under the alternative hypothesis. We applied all methods to study the mediation mechanism of DNA methylation sites in the pathway from adult socioeconomic status to glycated hemoglobin level using data from the Multi-Ethnic Study of Atherosclerosis (MESA). We provide guidelines for choosing the optimal mediation hypothesis testing method in practice and develop an R package medScan available on the CRAN for implementing all the six methods.

Mediation analysis of multiple mediators with incomplete omics data

There is an increasing interest in using multiple types of omics features (e.g., DNA sequences, RNA expressions, methylation, protein expressions, and metabolic profiles) to study how the relationships between phenotypes and genotypes may be mediated by other omics markers. Genotypes and phenotypes are typically available for all subjects in genetic studies, but typically, some omics data will be missing for some subjects, due to limitations such as cost and sample quality. In this article, we propose a powerful approach for mediation analysis that accommodates missing data among multiple mediators and allows for various interaction effects. We formulate the relationships among genetic variants, other omics measurements, and phenotypes through linear regression models. We derive the joint likelihood for models with two mediators, accounting for arbitrary patterns of missing values. Utilizing computationally efficient and stable algorithms, we conduct maximum likelihood estimation. Our methods produce unbiased and statistically efficient estimators. We demonstrate the usefulness of our methods through simulation studies and an application to the Metabolic Syndrome in Men study.

High-Dimensional DNA Methylation Mediates the Effect of Smoking on Crohn's Disease

Epigenome-wide mediation analysis aims to identify high-dimensional DNA methylation at cytosine–phosphate–guanine (CpG) sites that mediate the causal effect of linking smoking with Crohn’s disease (CD) outcome. Studies have shown that smoking has significant detrimental effects on the course of CD. So we assessed whether DNA methylation mediates the association between smoking and CD. Among 103 CD cases and 174 controls, we estimated whether the effects of smoking on CD are mediated through DNA methylation CpG sites, which we referred to as causal mediation effect. Based on the causal diagram, we first implemented sure independence screening (SIS) to reduce the pool of potential mediator CpGs from a very large to a moderate number; then, we implemented variable selection with de-sparsifying the LASSO regression. Finally, we carried out a comprehensive mediation analysis and conducted sensitivity analysis, which was adjusted for potential confounders of age, sex, and blood cell type proportions to estimate the mediation effects. Smoking was significantly associated with CD under odds ratio (OR) of 2.319 (95% CI: 1.603, 3.485, p < 0.001) after adjustment for confounders. Ninety-nine mediator CpGs were selected from SIS, and then, seven candidate CpGs were obtained by de-sparsifying the LASSO regression. Four of these CpGs showed statistical significance, and the average causal mediation effects (ACME) were attenuated from 0.066 to 0.126. Notably, three significant mediator CpGs had absolute sensitivity parameters of 0.40, indicating that these mediation effects were robust even when the assumptions were slightly violated. Genes (BCL3 and FKBP5) harboring these four CpGs were related to CD. These findings suggest that changes in methylation are involved in the mechanism by which smoking increases risk of CD.

A Review of High-Dimensional Mediation Analyses in DNA Methylation Studies

DNA methylation alterations have been widely studied as mediators of environmentally induced disease risks. With new advances in technique, epigenome-wide DNA methylation data (EWAS) have become the new standard for epigenetic studies in human populations. However, to date most epigenetic studies of mediation effects only involve selected (gene-specific) candidate methylation markers. There is an urgent need for appropriate analytical methods for EWAS mediation analysis. In this chapter, we provide an overview of recent advances on high-dimensional mediation analysis, with application to two DNA methylation data.

DNA methylation as a mediator of associations between the environment and chronic diseases: A scoping review on application of mediation analysis

DNA methylation (DNAm) is one of the most studied epigenetic modifications. DNAm has emerged as a key biological mechanism and biomarkers to test associations between environmental exposure and outcomes in epidemiological studies. Although previous studies have focused on associations between DNAm and either exposure/outcomes, it is useful to test for mediation of the association between exposure and outcome by DNAm. The purpose of this scoping review is to introduce the methodological essence of statistical mediation analysis and to examine emerging epidemiological research applying mediation analyses. We conducted this scoping review for published peer-reviewed journals on this topic using online databases (PubMed, Scopus, Cochrane, and CINAHL) ending in December 2020. We extracted a total of 219 articles by initial screening. After reviewing titles, abstracts, and full texts, a total of 69 articles were eligible for this review. The breakdown of studies assigned to each category was 13 for smoking (18.8%), 8 for dietary intake and famine (11.6%), 6 for other lifestyle factors (8.7%), 8 for clinical endpoints (11.6%), 22 for environmental chemical exposures (31.9%), 2 for socioeconomic status (SES) (2.9%), and 10 for genetic factors and race (14.5%). In this review, we provide an exposure-wide summary for the mediation analysis using DNAm levels. However, we found heterogenous methods and interpretations in mediation analysis with typical issues such as different cell compositions and tissue-specificity. Further accumulation of evidence with diverse exposures, populations and with rigorous methodology will be expected to provide further insight in the role of DNAm in disease susceptibility.

Review of Mendelian Randomization Studies on Ovarian Cancer

Ovarian cancer (OC) is one of the deadliest gynecological cancers worldwide. Previous observational epidemiological studies have revealed associations between modifiable environmental risk factors and OC risk. However, these studies are prone to confounding, measurement error, and reverse causation, undermining robust causal inference. Mendelian randomization (MR) analysis has been established as a reliable method to investigate the causal relationship between risk factors and diseases using genetic variants to proxy modifiable exposures. Over recent years, MR analysis in OC research has received extensive attention, providing valuable insights into the etiology of OC as well as holding promise for identifying potential therapeutic interventions. This review provides a comprehensive overview of the key principles and assumptions of MR analysis. Published MR studies focusing on the causality between different risk factors and OC risk are summarized, along with comprehensive analysis of the method and its future applications. The results of MR studies on OC showed that higher BMI and height, earlier age at menarche, endometriosis, schizophrenia, and higher circulating β-carotene and circulating zinc levels are associated with an increased risk of OC. In contrast, polycystic ovary syndrome; vitiligo; higher circulating vitamin D, magnesium, and testosterone levels; and HMG-CoA reductase inhibition are associated with a reduced risk of OC. MR analysis presents a2 valuable approach to understanding the causality between different risk factors and OC after full consideration of its inherent assumptions and limitations.

Biochemical Mechanisms Associating Alcohol Use Disorders with Cancers

Simple Summary

Of all yearly deaths attributable to alcohol consumption globally, approximately 12% are due to cancers, representing approximately 0.4 million deceased individuals. Ethanol metabolism disturbs cell biochemistry by targeting the structure and function of essential biomolecules (proteins, nucleic acids, and lipids) and by provoking alterations in cell programming that lead to cancer development and cancer malignancy. A better understanding of the metabolic and cell signaling realm affected by ethanol is paramount to designing effective treatments and preventive actions tailored to specific neoplasias.

Abstract

The World Health Organization identifies alcohol as a cause of several neoplasias of the oropharynx cavity, esophagus, gastrointestinal tract, larynx, liver, or female breast. We review ethanol’s nonoxidative and oxidative metabolism and one-carbon metabolism that encompasses both redox and transfer reactions that influence crucial cell proliferation machinery. Ethanol favors the uncontrolled production and action of free radicals, which interfere with the maintenance of essential cellular functions. We focus on the generation of protein, DNA, and lipid adducts that interfere with the cellular processes related to growth and differentiation. Ethanol’s effects on stem cells, which are responsible for building and repairing tissues, are reviewed. Cancer stem cells (CSCs) of different origins suffer disturbances related to the expression of cell surface markers, enzymes, and transcription factors after ethanol exposure with the consequent dysregulation of mechanisms related to cancer metastasis or resistance to treatments. Our analysis aims to underline and discuss potential targets that show more sensitivity to ethanol’s action and identify specific metabolic routes and metabolic realms that may be corrected to recover metabolic homeostasis after pharmacological intervention. Specifically, research should pay attention to re-establishing metabolic fluxes by fine-tuning the functioning of specific pathways related to one-carbon metabolism and antioxidant processes.

Longitudinal data reveal strong genetic and weak non-genetic components of ethnicity-dependent blood DNA methylation levels

Epigenetic architecture is influenced by genetic and environmental factors, but little is known about their relative contributions or longitudinal dynamics. Here, we studied DNA methylation (DNAm) at over 750,000 CpG sites in mononuclear blood cells collected at birth and age 7 from 196 children of primarily self-reported Black and Hispanic ethnicities to study race-associated DNAm patterns. We developed a novel Bayesian method for high-dimensional longitudinal data and showed that race-associated DNAm patterns at birth and age 7 are nearly identical. Additionally, we estimated that up to 51% of all self-reported race-associated CpGs had race-dependent DNAm levels that were mediated through local genotype and, quite surprisingly, found that genetic factors explained an overwhelming majority of the variation in DNAm levels at other, previously identified, environmentally-associated CpGs. These results indicate that race-associated blood DNAm patterns in particular, and blood DNAm levels in general, are primarily driven by genetic factors, and are not as sensitive to environmental exposures as previously suggested, at least during the first 7 years of life.

Statistical methods for mediation analysis in the era of high-throughput genomics: Current successes and future challenges

Mediation analysis investigates the intermediate mechanism through which an exposure exerts its influence on the outcome of interest. Mediation analysis is becoming increasingly popular in high-throughput genomics studies where a common goal is to identify molecular-level traits, such as gene expression or methylation, which actively mediate the genetic or environmental effects on the outcome. Mediation analysis in genomics studies is particularly challenging, however, thanks to the large number of potential mediators measured in these studies as well as the composite null nature of the mediation effect hypothesis. Indeed, while the standard univariate and multivariate mediation methods have been well-established for analyzing one or multiple mediators, they are not well-suited for genomics studies with a large number of mediators and often yield conservative p-values and limited power. Consequently, over the past few years many new high-dimensional mediation methods have been developed for analyzing the large number of potential mediators collected in high-throughput genomics studies. In this work, we present a thorough review of these important recent methodological advances in high-dimensional mediation analysis. Specifically, we describe in detail more than ten high-dimensional mediation methods, focusing on their motivations, basic modeling ideas, specific modeling assumptions, practical successes, methodological limitations, as well as future directions. We hope our review will serve as a useful guidance for statisticians and computational biologists who develop methods of high-dimensional mediation analysis as well as for analysts who apply mediation methods to high-throughput genomics studies.

Gene-based mediation analysis in epigenetic studies

Mediation analysis has been a useful tool for investigating the effect of mediators that lie in the path from the independent variable to the outcome. With the increasing dimensionality of mediators such as in (epi)genomics studies, high-dimensional mediation model is needed. In this work, we focus on epigenetic studies with the goal to identify important DNA methylations that act as mediators between an exposure disease outcome. Specifically, we focus on gene-based high-dimensional mediation analysis implemented with kernel principal component analysis to capture potential nonlinear mediation effect. We first review the current high-dimensional mediation models and then propose two gene-based analytical approaches: gene-based high-dimensional mediation analysis based on linearity assumption between mediators and outcome (gHMA-L) and gene-based high-dimensional mediation analysis based on nonlinearity assumption (gHMA-NL). Since the underlying true mediation relationship is unknown in practice, we further propose an omnibus test of gene-based high-dimensional mediation analysis (gHMA-O) by combing gHMA-L and gHMA-NL. Extensive simulation studies show that gHMA-L performs better under the model linear assumption and gHMA-NL does better under the model nonlinear assumption, while gHMA-O is a more powerful and robust method by combining the two. We apply the proposed methods to two datasets to investigate genes whose methylation levels act as important mediators in the relationship: (1) between alcohol consumption and epithelial ovarian cancer risk using data from the Mayo Clinic Ovarian Cancer Case-Control Study and (2) between childhood maltreatment and comorbid post-traumatic stress disorder and depression in adulthood using data from the Gray Trauma Project.

A region-based method for causal mediation analysis of DNA methylation data

Exposure to environmental factors can affect DNA methylation at a 5'-cytosine-phosphate-guanine-3' (CpG) site or a genomic region, which can then affect an outcome. In other words, environmental effects on an outcome could be mediated by DNA methylation. To date, single CpG-site-based mediation analysis has been employed extensively. More recently, however, there has been considerable interest in studying differentially methylated regions (DMRs), both because DMRs are more likely to have functional effects than single CpG sites and because testing DMRs reduces multiple testing. In this report, we propose a novel causal mediation approach under the counterfactual framework to test the significance of total (TE), direct (DE), and indirect effects (IE) of predictors on response variable with a methylated region (MR) as the mediator (denoted as MR-Mediation). Functional linear transformation is used to reduce the possible high dimension of the CpG sites in a predefined MR and to account for their location information. In our simulation studies, MR-Mediation retained the desired Type I error rates for TE, DE, and IE tests. Furthermore, MR-Mediation had better power performance than testing mean methylation level as the mediator in most considered scenarios, especially for IE (i.e., mediated effect) test, which could be more interesting than the other two effect tests. We further illustrate our proposed method by analysing the methylation mediated effect of exposure to gun violence on total immunoglobulin E or atopic asthma among participants in the Epigenetic Variation and Childhood Asthma in Puerto Ricans study.

DNA Methylation in Ovarian Cancer Susceptibility

Simple Summary

It is well established that ovarian cancer “runs in families”, where ovarian and other cancers (commonly breast cancer) occur at early ages at onset and in multiple generations. After decades of genetic studies, rare high-risk genetic mutations in cancer susceptibility genes and over 40 common genetic variants with much smaller risks have been identified. However, based on familial studies, we know that additional heritable genetic risk factors exist. It is possible that epigenetic variation—differences in how DNA is read, and which genes are actively expressed (or not) —also contributes to ovarian cancer susceptibility. This review summarizes the current collection of epidemiological studies that have investigated the role of DNA methylation—one type of epigenetic mechanism—in the risk of ovarian cancer.

Abstract

Epigenetic alterations are somatically acquired over the lifetime and during neoplastic transformation but may also be inherited as widespread ‘constitutional’ alterations in normal tissues that can cause cancer predisposition. Epithelial ovarian cancer (EOC) has an established genetic susceptibility and mounting epidemiological evidence demonstrates that DNA methylation (DNAm) intermediates as well as independently contributes to risk. Targeted studies of known EOC susceptibility genes (CSGs) indicate rare, constitutional BRCA1 promoter methylation increases familial and sporadic EOC risk. Blood-based epigenome-wide association studies (EWAS) for EOC have detected a total of 2846 differentially methylated probes (DMPs) with 71 genes replicated across studies despite significant heterogeneity. While EWAS detect both symptomatic and etiologic DMPs, adjustments and analytic techniques may enrich risk associations, as evidenced by the detection of dysregulated methylation of BNC2—a known CSG identified by genome-wide associations studies (GWAS). Integrative genetic–epigenetic approaches have mapped methylation quantitative trait loci (meQTL) to EOC risk, revealing DNAm variations that are associated with nine GWAS loci and, further, one novel risk locus. Increasing efforts to mapping epigenome variation across populations and cell types will be key to decoding both the genomic and epigenomic causal pathways to EOC.

Testing Mediation Effects in High-Dimensional Epigenetic Studies

Mediation analysis has been a powerful tool to identify factors mediating the association between exposure variables and outcomes. It has been applied to various genomic applications with the hope to gain novel insights into the underlying mechanism of various diseases. Given the high-dimensional nature of epigenetic data, recent effort on epigenetic mediation analysis is to first reduce the data dimension by applying high-dimensional variable selection techniques, then conducting testing in a low dimensional setup. In this paper, we propose to assess the mediation effect by adopting a high-dimensional testing procedure which can produce unbiased estimates of the regression coefficients and can properly handle correlations between variables. When the data dimension is ultra-high, we first reduce the data dimension from ultra-high to high by adopting a sure independence screening (SIS) method. We apply the method to two high-dimensional epigenetic studies: one is to assess how DNA methylations mediate the association between alcohol consumption and epithelial ovarian cancer (EOC) status; the other one is to assess how methylation signatures mediate the association between childhood maltreatment and post-traumatic stress disorder (PTSD) in adulthood. We compare the performance of the method with its counterpart via simulation studies. Our method can be applied to other high-dimensional mediation studies where high-dimensional mediation variables are collected.

Genetic Data from Nearly 63,000 Women of European Descent Predicts DNA Methylation Biomarkers and Epithelial Ovarian Cancer Risk

DNA methylation is instrumental for gene regulation. Global changes in the epigenetic landscape have been recognized as a hallmark of cancer. However, the role of DNA methylation in epithelial ovarian cancer (EOC) remains unclear. In this study, high-density genetic and DNA methylation data in white blood cells from the Framingham Heart Study (N = 1,595) were used to build genetic models to predict DNA methylation levels. These prediction models were then applied to the summary statistics of a genome-wide association study (GWAS) of ovarian cancer including 22,406 EOC cases and 40,941 controls to investigate genetically predicted DNA methylation levels in association with EOC risk. Among 62,938 CpG sites investigated, genetically predicted methylation levels at 89 CpG were significantly associated with EOC risk at a Bonferroni-corrected threshold of P < 7.94 × 10-7. Of them, 87 were located at GWAS-identified EOC susceptibility regions and two resided in a genomic region not previously reported to be associated with EOC risk. Integrative analyses of genetic, methylation, and gene expression data identified consistent directions of associations across 12 CpG, five genes, and EOC risk, suggesting that methylation at these 12 CpG may influence EOC risk by regulating expression of these five genes, namely MAPT, HOXB3, ABHD8, ARHGAP27, and SKAP1. We identified novel DNA methylation markers associated with EOC risk and propose that methylation at multiple CpG may affect EOC risk via regulation of gene expression. SIGNIFICANCE: Identification of novel DNA methylation markers associated with EOC risk suggests that methylation at multiple CpG may affect EOC risk through regulation of gene expression.