High-dimensional multivariate mediation with application to neuroimaging data

Optogenetic fMRI reveals therapeutic circuits of subthalamic nucleus deep brain stimulation

While deep brain stimulation (DBS) is widely employed for managing motor symptoms in Parkinson's disease (PD), its exact circuit mechanisms remain controversial. To identify the neural targets affected by therapeutic DBS in PD, we analyzed DBS-evoked whole brain activity in female hemi-parkinsonian rats using function magnetic resonance imaging (fMRI). We delivered subthalamic nucleus (STN) DBS at various stimulation pulse repetition rates using optogenetics, allowing unbiased examinations of cell-type specific STN feed-forward neural activity. Unilateral STN optogenetic stimulation elicited pulse repetition rate-dependent alterations of blood-oxygenation-level-dependent (BOLD) signals in SNr (substantia nigra pars reticulata), GP (globus pallidus), and CPu (caudate putamen). Notably, these manipulations effectively ameliorated pathological circling behavior in animals expressing the kinetically faster Chronos opsin, but not in animals expressing ChR2. Furthermore, mediation analysis revealed that the pulse repetition rate-dependent behavioral rescue was significantly mediated by optogenetically induced activity changes in GP and CPu, but not in SNr. This suggests that the activation of GP and CPu are critically involved in the therapeutic mechanisms of STN DBS.

A Bayesian nonparametric approach for multiple mediators with applications in mental health studies

Mediation analysis with contemporaneously observed multiple mediators is a significant area of causal inference. Recent approaches for multiple mediators are often based on parametric models and thus may suffer from model misspecification. Also, much of the existing literature either only allow estimation of the joint mediation effect or estimate the joint mediation effect just as the sum of individual mediator effects, ignoring the interaction among the mediators. In this article, we propose a novel Bayesian nonparametric method that overcomes the two aforementioned drawbacks. We model the joint distribution of the observed data (outcome, mediators, treatment, and confounders) flexibly using an enriched Dirichlet process mixture with three levels. We use standardization (g-computation) to compute all possible mediation effects, including pairwise and all other possible interaction among the mediators. We thoroughly explore our method via simulations and apply our method to a mental health data from Wisconsin Longitudinal Study, where we estimate how the effect of births from unintended pregnancies on later life mental depression (CES-D) among the mothers is mediated through lack of self-acceptance and autonomy, employment instability, lack of social participation, and increased family stress. Our method identified significant individual mediators, along with some significant pairwise effects.

Causal mediation analysis using high-dimensional image mediator bounded in irregular domain with an application to breast cancer

Mammography is the primary breast cancer screening strategy. Recent methods have been developed using the mammogram image to improve breast cancer risk prediction. However, it is unclear on the extent to which the effect of risk factors on breast cancer risk is mediated through tissue features summarized in mammogram images and the extent to which it is through other pathways. While mediation analysis has been conducted using mammographic density (a summary measure within the image), the mammogram image is not necessarily well described by a single summary measure and, in addition, such a measure provides no spatial information about the relationship between the exposure risk factor and the risk of breast cancer. Thus, to better understand the role of the mammogram images that provide spatial information about the state of the breast tissue that is causally predictive of the future occurrence of breast cancer, we propose a novel method of causal mediation analysis using mammogram image mediator while accommodating the irregular shape of the breast. We apply the proposed method to data from the Joanne Knight Breast Health Cohort and leverage new insights on the decomposition of the total association between risk factor and breast cancer risk that was mediated by the texture of the underlying breast tissue summarized in the mammogram image.

Methods for mediation analysis with high-dimensional DNA methylation data: Possible choices and comparisons

Epigenetic researchers often evaluate DNA methylation as a potential mediator of the effect of social/environmental exposures on a health outcome. Modern statistical methods for jointly evaluating many mediators have not been widely adopted. We compare seven methods for high-dimensional mediation analysis with continuous outcomes through both diverse simulations and analysis of DNAm data from a large multi-ethnic cohort in the United States, while providing an R package for their seamless implementation and adoption. Among the considered choices, the best-performing methods for detecting active mediators in simulations are the Bayesian sparse linear mixed model (BSLMM) and high-dimensional mediation analysis (HDMA); while the preferred methods for estimating the global mediation effect are high-dimensional linear mediation analysis (HILMA) and principal component mediation analysis (PCMA). We provide guidelines for epigenetic researchers on choosing the best method in practice and offer suggestions for future methodological development.

A high-dimensional mediation model for a neuroimaging mediator: Integrating clinical, neuroimaging, and neurocognitive data to mitigate late effects in pediatric cancer

Pediatric cancer treatment, especially for brain tumors, can have profound and complicated late effects. With the survival rates increasing because of improved detection and treatment, a more comprehensive understanding of the impact of current treatments on neurocognitive function and brain structure is critically needed. A frontline medulloblastoma clinical trial (SJMB03) has collected data, including treatment, clinical, neuroimaging, and cognitive variables. Advanced methods for modeling and integrating these data are critically needed to understand the mediation pathway from the treatment through brain structure to neurocognitive outcomes. We propose an integrative Bayesian mediation analysis approach to model jointly a treatment exposure, a high-dimensional structural neuroimaging mediator, and a neurocognitive outcome and to uncover the mediation pathway. The high-dimensional imaging-related coefficients are modeled via a binary Ising-Gaussian Markov random field prior (BI-GMRF), addressing the sparsity, spatial dependency, and smoothness and increasing the power to detect brain regions with mediation effects. Numerical simulations demonstrate the estimation accuracy, power, and robustness. For the SJMB03 study, the BI-GMRF method has identified white matter microstructure that is damaged by cancer-directed treatment and impacts late neurocognitive outcomes. The results provide guidance on improving treatment planning to minimize long-term cognitive sequela for pediatric brain tumor patients.

Mediation Analysis for High-Dimensional Mediators and Outcomes with an Application to Multimodal Imaging Data

Multimodal neuroimaging data have attracted increasing attention for brain research. An integrated analysis of multimodal neuroimaging data and behavioral or clinical measurements provides a promising approach for comprehensively and systematically investigating the underlying neural mechanisms of different phenotypes. However, such an integrated data analysis is intrinsically challenging due to the complex interactive relationships between the multimodal multivariate imaging variables. To address this challenge, a novel multivariate-mediator and multivariate-outcome mediation model (MMO) is proposed to simultaneously extract the latent systematic mediation patterns and estimate the mediation effects based on a dense bi-cluster graph approach. A computationally efficient algorithm is developed for dense bicluster structure estimation and inference to identify the mediation patterns with multiple testing correction. The performance of the proposed method is evaluated by an extensive simulation analysis with comparison to the existing methods. The results show that MMO performs better in terms of both the false discovery rate and sensitivity compared to existing models. The MMO is applied to a multimodal imaging dataset from the Human Connectome Project to investigate the effect of systolic blood pressure on whole-brain imaging measures for the regional homogeneity of the blood oxygenation level-dependent signal through the cerebral blood flow.

Brain mediators of negative affect-induced physical symptom reporting in patients with functional somatic syndromes

Functional somatic syndromes (FSS) include fibromyalgia, irritable bowel syndrome (IBS), and others. In FSS patients, merely viewing negative affective pictures can elicit increased physical symptoms. Our aim was to investigate the neural mechanisms underlying such negative affect-induced physical symptoms in FSS patients. Thirty patients with fibromyalgia and/or IBS and 30 healthy controls (all women) watched neutral, positive and negative affective picture blocks during functional MRI scanning and rated negative affect and physical symptoms after every block. We compared brain-wide activation during negative versus neutral picture viewing in FSS patients versus controls using robust general linear model analysis. Further, we compared neurologic pain signature (NPS), stimulus intensity-independent pain signature (SIIPS) and picture-induced negative emotion signature (PINES) responses to the negative versus neutral affect contrast and investigated whether they mediated between-group differences in affective picture-induced physical symptom reporting. More physical symptoms were reported after viewing negative compared to neutral pictures, and this effect was larger in patients than controls (p = 0.025). Accordingly, patients showed stronger activation in somatosensory regions during negative versus neutral picture viewing. NPS, but not SIIPS nor PINES, responses were higher in patients than controls during negative versus neutral pictures (p = 0.026). These differential NPS responses partially mediated between-group differences in physical symptoms. In conclusion, picture-induced negative affect elicits physical symptoms in FSS patients as a result of activation of somatosensory and nociceptive brain patterns, supporting the idea that affect-driven alterations in processing of somatic signals is a critical mechanism underlying FSS.

The DNA methylation landscape of the human oxytocin receptor gene (OXTR): data-driven clusters and their relation to gene expression and childhood adversity

The oxytocin receptor gene (OXTR) is of interest when investigating the effects of early adversity on DNA methylation. However, there is heterogeneity regarding the selection of the most promising CpG sites to target for analyses. The goal of this study was to determine functionally relevant clusters of CpG sites within the OXTR CpG island in 113 mother-infant dyads, with 58 of the mothers reporting childhood maltreatment (CM). OXTR DNA methylation was analyzed in peripheral/umbilical blood mononuclear cells. Different complexity reduction approaches were used to reduce the 188 CpG sites into clusters of co-methylated sites. Furthermore, associations between OXTR DNA methylation (cluster- and site-specific level) and OXTR gene expression and CM were investigated in mothers. Results showed that, first, CpG sections differed strongly regarding their statistical utility for research of individual differences in DNA methylation. Second, cluster analyses and Partial Least Squares (PLS) suggested two clusters consisting of intron1/exon2 and the protein-coding region of exon3, respectively, as most strongly associated with outcome measures. Third, cross-validated PLS regression explained 7% of variance in CM, with low cross-validated variance explained for the prediction of gene expression. Fourth, substantial mother-child correspondence was observed in correlation patterns within the identified clusters, but only modest correspondence outside these clusters. This study makes an important contribution to the mapping of the DNA methylation landscape of the OXTR CpG island by highlighting clusters of CpG sites that show desirable statistical properties and predictive value. We provide a Companion Web Application to facilitate the choice of CpG sites.

Relation between white matter integrity, perfusion, and processing speed in early-stage schizophrenia

OBJECTIVE

Cerebral blood flow (CBF) plays a critical role in the maintenance of neuronal integrity, and CBF alterations have been linked to deleterious white matter changes. Several studies report CBF and white matter structural alterations individually. However, whether and how these pathological changes relate to each other remains elusive. By using our cohort of individuals with early-stage schizophrenia, we investigated the relationship between CBF and white matter structure.

METHOD

We studied 51 early-stage schizophrenia patients and age- and sex-matched healthy controls. We investigated the relationship among tissue structure (assessed with diffusion weighted imaging), perfusion (accessed by pseudo-continuous arterial labeling imaging), and neuropsychological indices (focusing on processing speed). We focused on the corpus callosum, due to its major role in associative functions and directness on revealing the architecture of a major white matter bundle. We performed mediation analysis to identify the possible mechanism underlay the relationship among cognition and white matter integrity and perfusion.

RESULTS

The CBF and the fractional anisotropy (FA) were inversely correlated in the corpus callosum of early-stage schizophrenia patients. While CBF negatively correlated with processing speed, FA correlated positively with this cognitive measure. These results were not observed in controls. Mediation analysis revealed that the effect of FA on processing speed was mediated via the CBF.

CONCLUSIONS

We provide evidence of a relationship between brain perfusion and white matter integrity in the corpus callosum in early-stage schizophrenia. These findings may shed the light on underlying metabolic support for structural changes with cognitive impact in schizophrenia.

A modality-independent proto-organization of human multisensory areas

The processing of multisensory information is based upon the capacity of brain regions, such as the superior temporal cortex, to combine information across modalities. However, it is still unclear whether the representation of coherent auditory and visual events requires any prior audiovisual experience to develop and function. Here we measured brain synchronization during the presentation of an audiovisual, audio-only or video-only version of the same narrative in distinct groups of sensory-deprived (congenitally blind and deaf) and typically developed individuals. Intersubject correlation analysis revealed that the superior temporal cortex was synchronized across auditory and visual conditions, even in sensory-deprived individuals who lack any audiovisual experience. This synchronization was primarily mediated by low-level perceptual features, and relied on a similar modality-independent topographical organization of slow temporal dynamics. The human superior temporal cortex is naturally endowed with a functional scaffolding to yield a common representation across multisensory events.

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.

A machine learning based approach towards high-dimensional mediation analysis

Mediation analysis is used to investigate the role of intermediate variables (mediators) that lie in the path between an exposure and an outcome variable. While significant research has focused on developing methods for assessing the influence of mediators on the exposure-outcome relationship, current approaches do not easily extend to settings where the mediator is high-dimensional. These situations are becoming increasingly common with the rapid increase of new applications measuring massive numbers of variables, including brain imaging, genomics, and metabolomics. In this work, we introduce a novel machine learning based method for identifying high dimensional mediators. The proposed algorithm iterates between using a machine learning model to map the high-dimensional mediators onto a lower-dimensional space, and using the predicted values as input in a standard three-variable mediation model. Hence, the machine learning model is trained to maximize the likelihood of the mediation model. Importantly, the proposed algorithm is agnostic to the machine learning model that is used, providing significant flexibility in the types of situations where it can be used. We illustrate the proposed methodology using data from two functional Magnetic Resonance Imaging (fMRI) studies. First, using data from a task-based fMRI study of thermal pain, we combine the proposed algorithm with a deep learning model to detect distributed, network-level brain patterns mediating the relationship between stimulus intensity (temperature) and reported pain at the single trial level. Second, using resting-state fMRI data from the Human Connectome Project, we combine the proposed algorithm with a connectome-based predictive modeling approach to determine brain functional connectivity measures that mediate the relationship between fluid intelligence and working memory accuracy. In both cases, our multivariate mediation model links exposure variables (thermal pain or fluid intelligence), high dimensional brain measures (single-trial brain activation maps or resting-state brain connectivity) and behavioral outcomes (pain report or working memory accuracy) into a single unified model. Using the proposed approach, we are able to identify brain-based measures that simultaneously encode the exposure variable and correlate with the behavioral outcome.

Methods for Mediation Analysis with High-Dimensional DNA Methylation Data: Possible Choices and Comparison

Epigenetic researchers often evaluate DNA methylation as a mediator between social/environmental exposures and disease, but modern statistical methods for jointly evaluating many mediators have not been widely adopted. We compare seven methods for high-dimensional mediation analysis with continuous outcomes through both diverse simulations and analysis of DNAm data from a large national cohort in the United States, while providing an R package for their implementation. Among the considered choices, the best-performing methods for detecting active mediators in simulations are the Bayesian sparse linear mixed model by Song et al. (2020) and high-dimensional mediation analysis by Gao et al. (2019); while the superior methods for estimating the global mediation effect are high-dimensional linear mediation analysis by Zhou et al. (2021) and principal component mediation analysis by Huang and Pan (2016). We provide guidelines for epigenetic researchers on choosing the best method in practice and offer suggestions for future methodological development.

Associations between neighborhood socioeconomic status, parental education, and executive system activation in youth

Socioeconomic status (SES) can impact cognitive performance, including working memory (WM). As executive systems that support WM undergo functional neurodevelopment during adolescence, environmental stressors at both individual and community levels may influence cognitive outcomes. Here, we sought to examine how SES at the neighborhood and family level impacts task-related activation of the executive system during adolescence and determine whether this effect mediates the relationship between SES and WM performance. To address these questions, we studied 1,150 youths (age 8-23) that completed a fractal n-back WM task during functional magnetic resonance imaging at 3T as part of the Philadelphia Neurodevelopmental Cohort. We found that both higher neighborhood SES and parental education were associated with greater activation of the executive system to WM load, including the bilateral dorsolateral prefrontal cortex, posterior parietal cortex, and precuneus. The association of neighborhood SES remained significant when controlling for task performance, or related factors like exposure to traumatic events. Furthermore, high-dimensional multivariate mediation analysis identified distinct patterns of brain activity within the executive system that significantly mediated the relationship between measures of SES and task performance. These findings underscore the importance of multilevel environmental factors in shaping executive system function and WM in youth.

Persistent organic pollutant exposure contributes to Black/White differences in leukocyte telomere length in the National Health and Nutrition Examination Survey

Despite racial disparities in diseases of aging and premature mortality, non-Hispanic Black Americans tend to have longer leukocyte telomere length (LTL), a biomarker of cellular aging, than non-Hispanic White Americans. Previous findings suggest that exposure to certain persistent organic pollutants (POPs) is both racially-patterned and associated with longer LTL. We examine whether Black/White differences in LTL are explained by differences in exposure to 15 POPs by estimating the indirect effect (IE) of self-reported race on LTL that is mediated through nine polychlorinated biphenyls (PCBs), three furans, and three dioxins, as well as their mixtures. Our study population includes 1,251 adults from the 1999-2000 and 2001-2002 cycles of the cross-sectional National Health and Nutrition Examination Survey. We characterized single-pollutant mediation effects by constructing survey-weighted linear regression models. We also implemented various approaches to quantify a global mediation effect of all POPs, including unpenalized linear regression, ridge regression, and examination of three summary exposure scores. We found support for the hypothesis that exposure to PCBs partially mediates Black/White differences in LTL. In single-pollutant models, there were significant IEs of race on LTL through six individual PCBs (118, 138, 153, 170, 180, and 187). Ridge regression (0.013, CI 0.001, 0.023; 26.0% mediated) and models examining summative exposure scores with linear combinations derived from principal components analysis (0.019, CI 0.009, 0.029; 34.8% mediated) and Toxic Equivalency Quotient (TEQ) scores (0.016, CI 0.005, 0.026; 28.8% mediated) showed significant IEs when incorporating survey weights. Exposures to individual POPs and their mixtures, which may arise from residential and occupational segregation, may help explain why Black Americans have longer LTL than their White counterparts, providing an environmental explanation for counterintuitive race differences in cellular aging.

Cerebral blood flow and cardiovascular risk effects on resting brain regional homogeneity

Regional homogeneity (ReHo) is a measure of local functional brain connectivity that has been reported to be altered in a wide range of neuropsychiatric disorders. Computed from brain resting-state functional MRI time series, ReHo is also sensitive to fluctuations in cerebral blood flow (CBF) that in turn may be influenced by cerebrovascular health. We accessed cerebrovascular health with Framingham cardiovascular risk score (FCVRS). We hypothesize that ReHo signal may be influenced by regional CBF; and that these associations can be summarized as FCVRS→CBF→ReHo. We used three independent samples to test this hypothesis. A test-retest sample of N = 30 healthy volunteers was used for test-retest evaluation of CBF effects on ReHo. Amish Connectome Project (ACP) sample (N = 204, healthy individuals) was used to evaluate association between FCVRS and ReHo and testing if the association diminishes given CBF. The UKBB sample (N = 6,285, healthy participants) was used to replicate the effects of FCVRS on ReHo. We observed strong CBF→ReHo links (p<2.5 × 10^-3) using a three-point longitudinal sample. In ACP sample, marginal and partial correlations analyses demonstrated that both CBF and FCVRS were significantly correlated with the whole-brain average (p<10^-6) and regional ReHo values, with the strongest correlations observed in frontal, parietal, and temporal areas. Yet, the association between ReHo and FCVRS became insignificant once the effect of CBF was accounted for. In contrast, CBF→ReHo remained significantly linked after adjusting for FCVRS and demographic covariates (p<10^-6). Analysis in N = 6,285 replicated the FCVRS→ReHo effect (p = 2.7 × 10^-27). In summary, ReHo alterations in health and neuropsychiatric illnesses may be partially driven by region-specific variability in CBF, which is, in turn, influenced by cardiovascular factors.

Bayesian network mediation analysis with application to the brain functional connectome

The brain functional connectome, the collection of interconnected neural circuits along functional networks, facilitates a cutting-edge understanding of brain functioning, and has a potential to play a mediating role within the effect pathway between an exposure and an outcome. While existing mediation analytic approaches are capable of providing insight into complex processes, they mainly focus on a univariate mediator or mediator vector, without considering network-variate mediators. To fill the methodological gap and accomplish this exciting and urgent application, in the article, we propose an integrative mediation analysis under a Bayesian paradigm with networks entailing the mediation effect. To parameterize the network measurements, we introduce individually specified stochastic block models with unknown block allocation, and naturally bridge effect elements through the latent network mediators induced by the connectivity weights across network modules. To enable the identification of truly active mediating components, we simultaneously impose a feature selection across network mediators. We show the superiority of our model in estimating different effect components and selecting active mediating network structures. As a practical illustration of this approach's application to network neuroscience, we characterize the relationship between a therapeutic intervention and opioid abstinence as mediated by brain functional sub-networks.

Migraine disability, pain catastrophizing, and headache severity are associated with evoked pain and targeted by mind-body therapy

ABSTRACT

Meta-analysis suggests that migraine patients are no more sensitive to experimentally evoked pain than healthy control subjects. At the same time, studies have linked some migraine symptoms to quantitative sensory testing (QST) profiles. Unfortunately, previous studies associating migraine symptoms and QST have important methodological shortcomings, stemming from small sample sizes, and frequent use of univariate statistics for multivariate research questions. In the current study, we seek to address these limitations by using a large sample of episodic migraine patients (n = 103) and a multivariate analysis that associates pain ratings from many thermal intensities simultaneously with 12 clinical measures ranging from headache frequency to sleep abnormalities. We identified a single dimension of association between thermal QST and migraine symptoms that relates to pain ratings for all stimulus intensities and a subset of migraine symptoms relating to disability (Headache Impact Test 6 and Brief Pain Inventory interference), catastrophizing (Pain Catastrophizing Scale), and pain severity (average headache pain, Brief Pain Inventory severity, and Short-Form McGill Pain Questionnaire 2). Headache frequency, allodynia, affect, and sleep disturbances were unrelated to this dimension. Consistent with previous research, we did not observe any difference in QST ratings between migraine patients and healthy control subjects. Additionally, we found that the linear combination of symptoms related to QST was modified by the mind-body therapy enhanced mindfulness-based stress reduction (MBSR+). These results suggest that QST has a selective relationship with pain symptoms even in the absence of between-subjects differences between chronic pain patients and healthy control subjects.

HIMA2: high-dimensional mediation analysis and its application in epigenome-wide DNA methylation data

Mediation analysis plays a major role in identifying significant mediators in the pathway between environmental exposures and health outcomes. With advanced data collection technology for large-scale studies, there has been growing research interest in developing methodology for high-dimensional mediation analysis. In this paper we present HIMA2, an extension of the HIMA method (Zhang in Bioinformatics 32:3150-3154, 2016). First, the proposed HIMA2 reduces the dimension of mediators to a manageable level based on the sure independence screening (SIS) method (Fan in J R Stat Soc Ser B 70:849-911, 2008). Second, a de-biased Lasso procedure is implemented for estimating regression parameters. Third, we use a multiple-testing procedure to accurately control the false discovery rate (FDR) when testing high-dimensional mediation hypotheses. We demonstrate its practical performance using Monte Carlo simulation studies and apply our method to identify DNA methylation markers which mediate the pathway from smoking to reduced lung function in the Coronary Artery Risk Development in Young Adults (CARDIA) Study.

Causal Inference in Radiomics: Framework, Mechanisms, and Algorithms

The widespread use of machine learning algorithms in radiomics has led to a proliferation of flexible prognostic models for clinical outcomes. However, a limitation of these techniques is their black-box nature, which prevents the ability for increased mechanistic phenomenological understanding. In this article, we develop an inferential framework for estimating causal effects with radiomics data. A new challenge is that the exposure of interest is latent so that new estimation procedures are needed. We leverage a multivariate version of partial least squares for causal effect estimation. The methodology is illustrated with applications to two radiomics datasets, one in osteosarcoma and one in glioblastoma.

On optimal two-stage testing of multiple mediators

Mediation analysis in high‐dimensional settings often involves identifying potential mediators among a large number of measured variables. For this purpose, a two‐step familywise error rate procedure called ScreenMin has been recently proposed. In ScreenMin, variables are first screened and only those that pass the screening are tested. The proposed data‐independent threshold for selection has been shown to guarantee asymptotic familywise error rate. In this work, we investigate the impact of the threshold on the finite‐sample familywise error rate. We derive a power maximizing threshold and show that it is well approximated by an adaptive threshold of Wang et al. (2016, arXiv preprint arXiv:1610.03330). We illustrate the investigated procedures on a case‐control study examining the effect of fish intake on the risk of colorectal adenoma. We also apply our procedure in the context of replicability analysis to identify single nucleotide polymorphisms (SNP) associated with crop yield in two distinct environments.

Multimodal data integration via mediation analysis with high-dimensional exposures and mediators

Motivated by an imaging proteomics study for Alzheimer's disease (AD), in this article, we propose a mediation analysis approach with high‐dimensional exposures and high‐dimensional mediators to integrate data collected from multiple platforms. The proposed method combines principal component analysis with penalized least squares estimation for a set of linear structural equation models. The former reduces the dimensionality and produces uncorrelated linear combinations of the exposure variables, whereas the latter achieves simultaneous path selection and effect estimation while allowing the mediators to be correlated. Applying the method to the AD data identifies numerous interesting protein peptides, brain regions, and protein–structure–memory paths, which are in accordance with and also supplement existing findings of AD research. Additional simulations further demonstrate the effective empirical performance of the method.

Co-sparse Non-negative Matrix Factorization

Non-negative matrix factorization, which decomposes the input non-negative matrix into product of two non-negative matrices, has been widely used in the neuroimaging field due to its flexible interpretability with non-negativity property. Nowadays, especially in the neuroimaging field, it is common to have at least thousands of voxels while the sample size is only hundreds. The non-negative matrix factorization encounters both computational and theoretical challenge with such high-dimensional data, i.e., there is no guarantee for a sparse and part-based representation of data. To this end, we introduce a co-sparse non-negative matrix factorization method to high-dimensional data by simultaneously imposing sparsity in both two decomposed matrices. Instead of adding some sparsity induced penalty such as l ₁ norm, the proposed method directly controls the number of non-zero elements, which can avoid the bias issues and thus yield more accurate results. We developed an alternative primal-dual active set algorithm to derive the co-sparse estimator in a computationally efficient way. The simulation studies showed that our method achieved better performance than the state-of-art methods in detecting the basis matrix and recovering signals, especially under the high-dimensional scenario. In empirical experiments with two neuroimaging data, the proposed method successfully detected difference between Alzheimer's patients and normal person in several brain regions, which suggests that our method may be a valuable toolbox for neuroimaging studies.

A new Mendelian Randomization method to estimate causal effects of multivariable brain imaging exposures

The advent of simultaneously collected imaging-genetics data in large study cohorts provides an unprecedented opportunity to assess the causal effect of brain imaging traits on externally measured experimental results (e.g., cognitive tests) by treating genetic variants as instrumental variables. However, classic Mendelian Randomization methods are limited when handling high-throughput imaging traits as exposures to identify causal effects. We propose a new Mendelian Randomization framework to jointly select instrumental variables and imaging exposures, and then estimate the causal effect of multivariable imaging data on the outcome. We validate the proposed method with extensive data analyses and compare it with existing methods. We further apply our method to evaluate the causal effect of white matter microstructure integrity (WM) on cognitive function. The findings suggest that our method achieved better performance regarding sensitivity, bias, and false discovery rate compared to individually assessing the causal effect of a single exposure and jointly assessing the causal effect of multiple exposures without dimension reduction. Our application results indicated that WM measures across different tracts have a joint causal effect that significantly impacts the cognitive function among the participants from the UK Biobank.

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.

High-throughput mediation analysis of human proteome and metabolome identifies mediators of post-bariatric surgical diabetes control

To improve the power of mediation in high-throughput studies, here we introduce High-throughput mediation analysis (Hitman), which accounts for direction of mediation and applies empirical Bayesian linear modeling. We apply Hitman in a retrospective, exploratory analysis of the SLIMM-T2D clinical trial in which participants with type 2 diabetes were randomized to Roux-en-Y gastric bypass (RYGB) or nonsurgical diabetes/weight management, and fasting plasma proteome and metabolome were assayed up to 3 years. RYGB caused greater improvement in HbA1c, which was mediated by growth hormone receptor (GHR). GHR’s mediation is more significant than clinical mediators, including BMI. GHR decreases at 3 months postoperatively alongside increased insulin-like growth factor binding proteins IGFBP1/BP2; plasma GH increased at 1 year. Experimental validation indicates (1) hepatic GHR expression decreases in post-bariatric rats; (2) GHR knockdown in primary hepatocytes decreases gluconeogenic gene expression and glucose production. Thus, RYGB may induce resistance to diabetogenic effects of GH signaling.

Trial Registration: Clinicaltrials.gov NCT01073020.

Factors underlying the effects of gastric bypass surgery on glucose homeostasis are incompletely understood. Here the authors developed and applied high-throughput mediation analysis to identify proteome/metabolome mediators of improved glucose homeostasis after to gastric bypass surgery, and report that improved glycemia was mediated by the growth hormone receptor.

Bayesian Sparse Mediation Analysis with Targeted Penalization of Natural Indirect Effects

Causal mediation analysis aims to characterize an exposure's effect on an outcome and quantify the indirect effect that acts through a given mediator or a group of mediators of interest. With the increasing availability of measurements on a large number of potential mediators, like the epigenome or the microbiome, new statistical methods are needed to simultaneously accommodate high-dimensional mediators while directly target penalization of the natural indirect effect (NIE) for active mediator identification. Here, we develop two novel prior models for identification of active mediators in high-dimensional mediation analysis through penalizing NIEs in a Bayesian paradigm. Both methods specify a joint prior distribution on the exposure-mediator effect and mediator-outcome effect with either (a) a four-component Gaussian mixture prior or (b) a product threshold Gaussian prior. By jointly modeling the two parameters that contribute to the NIE, the proposed methods enable penalization on their product in a targeted way. Resultant inference can take into account the four-component composite structure underlying the NIE. We show through simulations that the proposed methods improve both selection and estimation accuracy compared to other competing methods. We applied our methods for an in-depth analysis of two ongoing epidemiologic studies: the Multi-Ethnic Study of Atherosclerosis (MESA) and the LIFECODES birth cohort. The identified active mediators in both studies reveal important biological pathways for understanding disease mechanisms.

Penalized mediation models for multivariate data

Statistical methods to integrate multiple layers of data, from exposures to intermediate traits to outcome variables, are needed to guide interpretation of complex data sets for which variables are likely contributing in a causal pathway from exposure to outcome. Statistical mediation analysis based on structural equation models provide a general modeling framework, yet they can be difficult to apply to high-dimensional data and they are not automated to select the best fitting model. To overcome these limitations, we developed novel algorithms and software to simultaneously evaluate multiple exposure variables, multiple intermediate traits, and multiple outcome variables. Our penalized mediation models are computationally efficient and simulations demonstrate that they produce reliable results for large data sets. Application of our methods to a study of vascular disease demonstrates their utility to identify novel direct effects of single-nucleotide polymorphisms (SNPs) on coronary heart disease and peripheral artery disease, while disentangling the effects of SNPs on the intermediate risk factors including lipids, cigarette smoking, systolic blood pressure, and type 2 diabetes.

Varying Coefficient Mediation Model and Application to Analysis of Behavioral Economics Data

This article is concerned with causal mediation analysis with varying indirect and direct effects. We propose a varying coefficient mediation model, which can also be viewed as an extension of moderation analysis on a causal diagram. We develop a new estimation procedure for the direct and indirect effects based on B-splines. Under mild conditions, rates of convergence and asymptotic distributions of the resulting estimates are established. We further propose a F-type test for the direct effect. We conduct simulation study to examine the finite sample performance of the proposed methodology, and apply the new procedures for empirical analysis of behavioral economics data.

Pathway Lasso: Pathway Estimation and Selection with High-Dimensional Mediators

In many scientific studies, it becomes increasingly important to delineate the pathways through a large number of mediators, such as genetic and brain mediators. Structural equation modeling (SEM) is a popular technique to estimate the pathway effects, commonly expressed as the product of coefficients. However, it becomes unstable and computationally challenging to fit such models with high-dimensional mediators. This paper proposes a sparse mediation model using a regularized SEM approach, where sparsity means that a small number of mediators have a nonzero mediation effect between a treatment and an outcome. To address the model selection challenge, we innovate by introducing a new penalty called Pathway Lasso. This penalty function is a convex relaxation of the non-convex product function for the mediation effects, and it enables a computationally tractable optimization criterion to estimate and select pathway effects simultaneously. We develop a fast ADMM-type algorithm to compute the model parameters, and we show that the iterative updates can be expressed in closed form. We also prove the asymptotic consistency of our Pathway Lasso estimator for the mediation effect. On both simulated data and an fMRI data set, the proposed approach yields higher pathway selection accuracy and lower estimation bias than competing methods.

Mediation analysis for survival data with High-Dimensional mediators

MOTIVATION

Mediation analysis has become a prevalent method to identify causal pathway(s) between an independent variable and a dependent variable through intermediate variable(s). However, little work has been done when the intermediate variables (mediators) are high-dimensional and the outcome is a survival endpoint. In this paper, we introduce a novel method to identify potential mediators in a causal framework of high-dimensional Cox regression.

RESULTS

We first reduce the data dimension through a mediation-based sure independence screening (SIS) method. A de-biased Lasso inference procedure is used for Cox's regression parameters. We adopt a multiple-testing procedure to accurately control the false discovery rate (FDR) when testing high-dimensional mediation hypotheses. Simulation studies are conducted to demonstrate the performance of our method. We apply this approach to explore the mediation mechanisms of 379,330 DNA methylation markers between smoking and overall survival among lung cancer patients in the TCGA lung cancer cohort. Two methylation sites (cg08108679 and cg26478297) are identified as potential mediating epigenetic markers.

AVAILABILITY

Our proposed method is available with the R package HIMA at https://cran.r-project.org/web/packages/HIMA/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Performance of approaches relying on multidimensional intermediary data to decipher causal relationships between the exposome and health: A simulation study under various causal structures

Challenges in the assessment of the health effects of the exposome, defined as encompassing all environmental exposures from the prenatal period onwards, include a possibly high rate of false positive signals. It might be overcome using data dimension reduction techniques. Data from the biological layers lying between the exposome and its possible health consequences, such as the methylome, may help reducing exposome dimension. We aimed to quantify the performances of approaches relying on the incorporation of an intermediary biological layer to relate the exposome and health, and compare them with agnostic approaches ignoring the intermediary layer. We performed a Monte-Carlo simulation, in which we generated realistic exposome and intermediary layer data by sampling with replacement real data from the Helix exposome project. We generated a Gaussian outcome assuming linear relationships between the three data layers, in 2381 scenarios under five different causal structures, including mediation and reverse causality. We tested 3 agnostic methods considering only the exposome and the health outcome: ExWAS (for Exposome-Wide Association study), DSA, LASSO; and 3 methods relying on an intermediary layer: two implementations of our new oriented Meet-in-the-Middle (oMITM) design, using ExWAS and DSA, and a mediation analysis using ExWAS. Methods' performances were assessed through their sensitivity and FDP (False-Discovery Proportion). The oMITM-based methods generally had lower FDP than the other approaches, possibly at a cost in terms of sensitivity; FDP was in particular lower under a structure of reverse causality and in some mediation scenarios. The oMITM-DSA implementation showed better performances than oMITM-ExWAS, especially in terms of FDP. Among the agnostic approaches, DSA showed the highest performance. Integrating information from intermediary biological layers can help lowering FDP in studies of the exposome health effects; in particular, oMITM seems less sensitive to reverse causality than agnostic exposome-health association studies.

Testing for mediation effect with application to human microbiome data

Mediation analysis has been commonly used to study the effect of an exposure on an outcome through a mediator. In this paper, we are interested in exploring the mediation mechanism of microbiome, whose special features make the analysis challenging. We consider the isometric logratio transformation of the relative abundance as the mediator variable. Then, we present a de-biased Lasso estimate for the mediator of interest and derive its standard error estimator, which can be used to develop a test procedure for the interested mediation effect. Extensive simulation studies are conducted to assess the performance of our method. We apply the proposed approach to test the mediation effect of human gut microbiome between the dietary fiber intake and body mass index.

High-Dimensional Mediation Analysis With Confounders in Survival Models

Mediation analysis is a common statistical method for investigating the mechanism of environmental exposures on health outcomes. Previous studies have extended mediation models with a single mediator to high-dimensional mediators selection. It is often assumed that there are no confounders that influence the relations among the exposure, mediator, and outcome. This is not realistic for the observational studies. To accommodate the potential confounders, we propose a concise and efficient high-dimensional mediation analysis procedure using the propensity score for adjustment. Results from simulation studies demonstrate the proposed procedure has good performance in mediator selection and effect estimation compared with methods that ignore all confounders. Of note, as the sample size increases, the performance of variable selection and mediation effect estimation is as well as the results shown in the method which include all confounders as covariates in the mediation model. By applying this procedure to a TCGA lung cancer data set, we find that lung cancer patients who had serious smoking history have increased the risk of death via the methylation markers cg21926276 and cg20707991 with significant hazard ratios of 1.2093 (95% CI: 1.2019-1.2167) and 1.1388 (95% CI: 1.1339-1.1438), respectively.

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.

Application of an analytical framework for multivariate mediation analysis of environmental data

Diverse toxicological mechanisms may mediate the impact of environmental toxicants (phthalates, phenols, polycyclic aromatic hydrocarbons, and metals) on pregnancy outcomes. In this study, we introduce an analytical framework for multivariate mediation analysis to identify mediation pathways (q = 61 mediators) in the relationship between environmental toxicants (p = 38 analytes) and gestational age at delivery. Our analytical framework includes: (1) conducting pairwise mediation for unique exposure-mediator combinations, (2) exposure dimension reduction by estimating environmental risk scores, and (3) multivariate mediator analysis using either Bayesian shrinkage mediation analysis, population value decomposition, or mediation pathway penalization. Dimension reduction demonstrates that a one-unit increase in phthalate risk score is associated with a total effect of 1.07 lower gestational age (in weeks) at delivery (95% confidence interval: 0.48-1.67) and eicosanoids from the cytochrome p450 pathway mediated 26% of this effect (95% confidence interval: 4-63%). Eicosanoid products derived from the cytochrome p450 pathway may be important mediators of phthalate toxicity.

Identifying neural signatures mediating behavioral symptoms and psychosis onset: High-dimensional whole brain functional mediation analysis

Along the pathway from behavioral symptoms to the development of psychotic disorders sits the multivariate mediating brain. The functional organization and structural topography of large-scale multivariate neural mediators among patients with brain disorders, however, are not well understood. Here, we design a high-dimensional brain-wide functional mediation framework to investigate brain regions that intermediate between baseline behavioral symptoms and future conversion to full psychosis among individuals at clinical high risk (CHR). Using resting-state functional magnetic resonance imaging (fMRI) data from 263 CHR subjects, we extract an α brain atlas and a β brain atlas: the former underlines brain areas associated with prodromal symptoms and the latter highlights brain areas associated with disease onset. In parallel, we identify and separate mediators that potentially positively and negatively mediate symptoms and psychosis, respectively, and quantify the effect of each neural mediator on disease development. Taken together, these results paint a brain-wide picture of neural markers that are potentially mediating behavioral symptoms and the development of psychotic disorders; additionally, they underscore a statistical framework that is useful to uncover large-scale intermediating variables in a regulatory biological system.

Bayesian shrinkage estimation of high dimensional causal mediation effects in omics studies

Causal mediation analysis aims to examine the role of a mediator or a group of mediators that lie in the pathway between an exposure and an outcome. Recent biomedical studies often involve a large number of potential mediators based on high-throughput technologies. Most of the current analytic methods focus on settings with one or a moderate number of potential mediators. With the expanding growth of -omics data, joint analysis of molecular-level genomics data with epidemiological data through mediation analysis is becoming more common. However, such joint analysis requires methods that can simultaneously accommodate high-dimensional mediators and that are currently lacking. To address this problem, we develop a Bayesian inference method using continuous shrinkage priors to extend previous causal mediation analysis techniques to a high-dimensional setting. Simulations demonstrate that our method improves the power of global mediation analysis compared to simpler alternatives and has decent performance to identify true nonnull contributions to the mediation effects of the pathway. The Bayesian method also helps us to understand the structure of the composite null cases for inactive mediators in the pathway. We applied our method to Multi-Ethnic Study of Atherosclerosis and identified DNA methylation regions that may actively mediate the effect of socioeconomic status on cardiometabolic outcomes.

Multimodal neuroimaging data integration and pathway analysis

With advancements in technology, the collection of multiple types of measurements on a common set of subjects is becoming routine in science. Some notable examples include multimodal neuroimaging studies for the simultaneous investigation of brain structure and function and multi-omics studies for combining genetic and genomic information. Integrative analysis of multimodal data allows scientists to interrogate new mechanistic questions. However, the data collection and generation of integrative hypotheses is outpacing available methodology for joint analysis of multimodal measurements. In this article, we study high-dimensional multimodal data integration in the context of mediation analysis. We aim to understand the roles that different data modalities play as possible mediators in the pathway between an exposure variable and an outcome. We propose a mediation model framework with two data types serving as separate sets of mediators and develop a penalized optimization approach for parameter estimation. We study both the theoretical properties of the estimator through an asymptotic analysis and its finite-sample performance through simulations. We illustrate our method with a multimodal brain pathway analysis having both structural and functional connectivity as mediators in the association between sex and language processing.

Estimation and inference for the indirect effect in high-dimensional linear mediation models

Mediation analysis is difficult when the number of potential mediators is larger than the sample size. In this paper we propose new inference procedures for the indirect effect in the presence of high-dimensional mediators for linear mediation models. We develop methods for both incomplete mediation, where a direct effect may exist, and complete mediation, where the direct effect is known to be absent. We prove consistency and asymptotic normality of our indirect effect estimators. Under complete mediation, where the indirect effect is equivalent to the total effect, we further prove that our approach gives a more powerful test compared to directly testing for the total effect. We confirm our theoretical results in simulations, as well as in an integrative analysis of gene expression and genotype data from a pharmacogenomic study of drug response. We present a novel analysis of gene sets to understand the molecular mechanisms of drug response, and also identify a genome-wide significant noncoding genetic variant that cannot be detected using standard analysis methods.

Challenges Raised by Mediation Analysis in a High-Dimension Setting

BACKGROUND

Mediation analysis is used in epidemiology to identify pathways through which exposures influence health. The advent of high-throughput (omics) technologies gives opportunities to perform mediation analysis with a high-dimension pool of covariates.

OBJECTIVE

We aimed to highlight some biostatistical issues of this expanding field of high-dimension mediation.

DISCUSSION

The mediation techniques used for a single mediator cannot be generalized in a straightforward manner to high-dimension mediation. Causal knowledge on the relation between covariates is required for mediation analysis, and it is expected to be more limited as dimension and system complexity increase. The methods developed in high dimension can be distinguished according to whether mediators are considered separately or as a whole. Methods considering each potential mediator separately do not allow efficient identification of the indirect effects when mutual influences exist among the mediators, which is expected for many biological (e.g., epigenetic) parameters. In this context, methods considering all potential mediators simultaneously, based, for example, on data reduction techniques, are more adapted to the causal inference framework. Their cost is a possible lack of ability to single out the causal mediators. Moreover, the ability of the mediators to predict the outcome can be overestimated, in particular because many machine-learning algorithms are optimized to increase predictive ability rather than their aptitude to make causal inference. Given the lack of overarching validated framework and the generally complex causal structure of high-dimension data, analysis of high-dimension mediation currently requires great caution and effort to incorporate a priori biological knowledge. https://doi.org/10.1289/EHP6240.

High-dimensional mediation analysis in survival models

Mediation analysis with high-dimensional DNA methylation markers is important in identifying epigenetic pathways between environmental exposures and health outcomes. There have been some methodology developments of mediation analysis with high-dimensional mediators. However, high-dimensional mediation analysis methods for time-to-event outcome data are still yet to be developed. To address these challenges, we propose a new high-dimensional mediation analysis procedure for survival models by incorporating sure independent screening and minimax concave penalty techniques for variable selection, with the Sobel and the joint method for significance test of indirect effect. The simulation studies show good performance in identifying correct biomarkers, false discovery rate control, and minimum estimation bias of the proposed procedure. We also apply this approach to study the causal pathway from smoking to overall survival among lung cancer patients potentially mediated by 365,307 DNA methylations in the TCGA lung cancer cohort. Mediation analysis using a Cox proportional hazards model estimates that patients who have serious smoking history increase the risk of lung cancer through methylation markers including cg21926276, cg27042065, and cg26387355 with significant hazard ratios of 1.2497(95%CI: 1.1121, 1.4045), 1.0920(95%CI: 1.0170, 1.1726), and 1.1489(95%CI: 1.0518, 1.2550), respectively. The three methylation sites locate in the three genes which have been showed to be associated with lung cancer event or overall survival. However, the three CpG sites (cg21926276, cg27042065 and cg26387355) have not been reported, which are newly identified as the potential novel epigenetic markers linking smoking and survival of lung cancer patients. Collectively, the proposed high-dimensional mediation analysis procedure has good performance in mediator selection and indirect effect estimation.

Multiple Brain Networks Mediating Stimulus-Pain Relationships in Humans

The brain transforms nociceptive input into a complex pain experience comprised of sensory, affective, motivational, and cognitive components. However, it is still unclear how pain arises from nociceptive input and which brain networks coordinate to generate pain experiences. We introduce a new high-dimensional mediation analysis technique to estimate distributed, network-level patterns that formally mediate the relationship between stimulus intensity and pain. We applied the model to a large-scale analysis of functional magnetic resonance imaging data (N = 284), focusing on brain mediators of the relationship between noxious stimulus intensity and trial-to-trial variation in pain reports. We identify mediators in both traditional nociceptive pathways and in prefrontal, midbrain, striatal, and default-mode regions unrelated to nociception in standard analyses. The whole-brain mediators are specific for pain versus aversive sounds and are organized into five functional networks. Brain mediators predicted pain ratings better than previous brain measures, including the neurologic pain signature (Wager et al. 2013). Our results provide a broader view of the networks underlying pain experience, as well as novel brain targets for interventions.

Complier stochastic direct effects: identification and robust estimation

Mediation analysis is critical to understanding the mechanisms underlying exposure-outcome relationships. In this paper, we identify the instrumental variable-direct effect of the exposure on the outcome not through the mediator, using randomization of the instrument. We call this estimand the complier stochastic direct effect (CSDE). To our knowledge, such an estimand has not previously been considered or estimated. We propose and evaluate several estimators for the CSDE: a ratio of inverse-probability of treatment-weighted estimators (IPTW), a ratio of estimating equation estimators (EE), a ratio of targeted minimum loss-based estimators (TMLE), and a TMLE that targets the CSDE directly. These estimators are applicable for a variety of study designs, including randomized encouragement trials, like the Moving to Opportunity housing voucher experiment we consider as an illustrative example, treatment discontinuities, and Mendelian randomization. We found the IPTW estimator to be the most sensitive to finite sample bias, resulting in bias of over 40% even when all models were correctly specified in a sample size of N=100. In contrast, the EE estimator and TMLE that targets the CSDE directly were far less sensitive. The EE and TML estimators also have advantages in terms of efficiency and reduced reliance on correct parametric model specification.

FWER and FDR control when testing multiple mediators

Motivation

The biological pathways linking exposures and disease risk are often poorly understood. To gain insight into these pathways, studies may try to identify biomarkers that mediate the exposure/disease relationship. Such studies often simultaneously test hundreds or thousands of biomarkers.

Results

We consider a set of m biomarkers and a corresponding set of null hypotheses, where the jth null hypothesis states that biomarker j does not mediate the exposure/disease relationship. We propose a Multiple Comparison Procedure (MCP) that rejects a set of null hypotheses or, equivalently, identifies a set of mediators, while asymptotically controlling the Family-Wise Error Rate (FWER) or False Discovery Rate (FDR). We use simulations to show that, compared to currently available methods, our proposed method has higher statistical power to detect true mediators. We then apply our method to a breast cancer study and identify nine metabolites that may mediate the known relationship between an increased BMI and an increased risk of breast cancer.

Availability and implementation

R package MultiMed on https://github.com/SiminaB/MultiMed.

Supplementary information

Supplementary data are available at Bioinformatics online.

Sparse Principal Component based High-Dimensional Mediation Analysis

Causal mediation analysis aims to quantify the intermediate effect of a mediator on the causal pathway from treatment to outcome. When dealing with multiple mediators, which are potentially causally dependent, the possible decomposition of pathway effects grows exponentially with the number of mediators. An existing approach incorporated the principal component analysis (PCA) to address this challenge based on the fact that the transformed mediators are conditionally independent given the orthogonality of the principal components (PCs). However, the transformed mediator PCs, which are linear combinations of original mediators, can be difficult to interpret. A sparse high-dimensional mediation analysis approach is proposed which adopts the sparse PCA method to the mediation setting. The proposed approach is applied to a task-based functional magnetic resonance imaging study, illustrating its ability to detect biologically meaningful results related to an identified mediator.