scholarly journals Semiparametric methods for estimation of a non-linear exposure–outcome relationship using instrumental variables with application to Mendelian randomization

2017 ◽  
Author(s):  
James R Staley ◽  
Stephen Burgess

AbstractMendelian randomization, the use of genetic variants as instrumental variables (IV), can test for and estimate the causal effect of an exposure on an outcome. Most IV methods assume that the function relating the exposure to the expected value of the outcome (the exposure-outcome relationship) is linear. However, in practice this assumption may not hold. Indeed, often the primary question of interest is to assess the shape of this relationship. We present two novel IV methods for investigating the shape of the exposure-outcome relationship: a fractional polynomial method and a piecewise linear method. We divide the population into strata using the exposure distribution, and estimate a causal effect, referred to as a localized average causal effect (LACE), in each stratum of population. The fractional polynomial method performs meta-regression on these LACE estimates. The piecewise linear method estimates a continuous piecewise linear function, the gradient of which is the LACE estimate in each stratum. Both methods were demonstrated in a simulation study to estimate the true exposure-outcome relationship well, particularly when the relationship was a fractional polynomial (for the fractional polynomial method) or was piecewise linear (for the piecewise linear method). The methods were used to investigate the shape of relationship of body mass index with systolic blood pressure and diastolic blood pressure.Availability and implementation: https://github.com/jrs95/nlmr

2020 ◽  
Author(s):  
Carlos Cinelli ◽  
Nathan LaPierre ◽  
Brian L. Hill ◽  
Sriram Sankararaman ◽  
Eleazar Eskin

ABSTRACTMendelian Randomization (MR) exploits genetic variants as instrumental variables to estimate the causal effect of an “exposure” trait on an “outcome” trait from observational data. However, the validity of such studies is threatened by population stratification, batch effects, and horizontal pleiotropy. Although a variety of methods have been proposed to partially mitigate those problems, residual biases may still remain, leading to highly statistically significant false positives in large genetic databases. Here, we describe a suite of sensitivity analysis tools for MR that enables investigators to properly quantify the robustness of their findings against these (and other) unobserved validity threats. Specifically, we propose the routine reporting of sensitivity statistics that can be used to readily quantify the robustness of a MR result: (i) the partial R2 of the genetic instrument with the exposure and the outcome traits; and, (ii) the robustness value of both genetic associations. These statistics quantify the minimal strength of violations of the MR assumptions that would be necessary to explain away the MR causal effect estimate. We also provide intuitive displays to visualize the sensitivity of the MR estimate to any degree of violation, and formal methods to bound the worst-case bias caused by violations in terms of multiples of the observed strength of principal components, batch effects, as well as putative pleiotropic pathways. We demonstrate how these tools can aid researchers in distinguishing robust from fragile findings, by showing that the MR estimate of the causal effect of body mass index (BMI) on diastolic blood pressure is relatively robust, whereas the MR estimate of the causal effect of BMI on Townsend deprivation index is relatively fragile.


Rheumatology ◽  
2020 ◽  
Author(s):  
Yi-Lin Dan ◽  
Peng Wang ◽  
Zhongle Cheng ◽  
Qian Wu ◽  
Xue-Rong Wang ◽  
...  

Abstract Objectives Several studies have reported increased serum/plasma adiponectin levels in SLE patients. This study was performed to estimate the causal effects of circulating adiponectin levels on SLE. Methods We selected nine independent single-nucleotide polymorphisms that were associated with circulating adiponectin levels (P < 5 × 10−8) as instrumental variables from a published genome-wide association study (GWAS) meta-analysis. The corresponding effects between instrumental variables and outcome (SLE) were obtained from an SLE GWAS analysis, including 7219 cases with 15 991 controls of European ancestry. Two-sample Mendelian randomization (MR) analyses with inverse-variance weighted, MR-Egger regression, weighted median and weight mode methods were used to evaluate the causal effects. Results The results of inverse-variance weighted methods showed no significantly causal associations of genetically predicted circulating adiponectin levels and the risk for SLE, with an odds ratio (OR) of 1.38 (95% CI 0.91, 1.35; P = 0.130). MR-Egger [OR 1.62 (95% CI 0.85, 1.54), P = 0.195], weighted median [OR 1.37 (95% CI 0.82, 1.35), P = 0.235) and weighted mode methods [OR 1.39 (95% CI 0.86, 1.38), P = 0.219] also supported no significant associations of circulating adiponectin levels and the risk for SLE. Furthermore, MR analyses in using SLE-associated single-nucleotide polymorphisms as an instrumental variable showed no associations of genetically predicted risk of SLE with circulating adiponectin levels. Conclusion Our study did not find evidence for a causal relationship between circulating adiponectin levels and the risk of SLE or of a causal effect of SLE on circulating adiponectin levels.


2020 ◽  
Vol 40 (2) ◽  
pp. 156-169 ◽  
Author(s):  
Christoph F. Kurz ◽  
Michael Laxy

Causal effect estimates for the association of obesity with health care costs can be biased by reversed causation and omitted variables. In this study, we use genetic variants as instrumental variables to overcome these limitations, a method that is often called Mendelian randomization (MR). We describe the assumptions, available methods, and potential pitfalls of using genetic information and how to address them. We estimate the effect of body mass index (BMI) on total health care costs using data from a German observational study and from published large-scale data. In a meta-analysis of several MR approaches, we find that models using genetic instruments identify additional annual costs of €280 for a 1-unit increase in BMI. This is more than 3 times higher than estimates from linear regression without instrumental variables (€75). We found little evidence of a nonlinear relationship between BMI and health care costs. Our results suggest that the use of genetic instruments can be a powerful tool for estimating causal effects in health economic evaluation that might be superior to other types of instruments where there is a strong association with a modifiable risk factor.


2021 ◽  
Author(s):  
Tim T Morris ◽  
Jon Heron ◽  
Eleanor Sanderson ◽  
George Davey Smith ◽  
Kate Tilling

Background Mendelian randomization (MR) is a powerful tool through which the causal effects of modifiable exposures on outcomes can be estimated from observational data. Most exposures vary throughout the life course, but MR is commonly applied to one measurement of an exposure (e.g., weight measured once between ages 40 and 60). It has been argued that MR provides biased causal effect estimates when applied to one measure of an exposure that varies over time. Methods We propose an approach that emphasises the liability that causes the entire exposure trajectory. We demonstrate this approach using simulations and an applied example. Results We show that rather than estimating the direct or total causal effect of changing the exposure value at a given time, MR estimates the causal effect of changing the liability as induced by a specific genotype that gives rise to the exposure at that time. As such, results from MR conducted at different time points are expected to differ (unless the liability of exposure is constant over time), as we demonstrate by estimating the effect of BMI measured at different ages on systolic blood pressure. Conclusions Practitioners should not interpret MR results as timepoint-specific direct or total causal effects, but as the effect of changing the liability that causes the entire exposure trajectory. Estimates of how the effects of a genetic variant on an exposure vary over time are needed to interpret timepoint-specific causal effects.


2021 ◽  
Author(s):  
Ken Chen ◽  
Zhenhuang Zhuang ◽  
Chunli Shao ◽  
Jilin Zheng ◽  
Qing Zhou ◽  
...  

Abstract ObjectivesTo investigate the roles of cardiometabolic factors (including blood pressure, blood lipids, thyroid function, body mass, and insulin sensitivity) in mediating the causal effect of type 2 diabetes (T2DM) on cardiovascular disease (CVD) outcomes. DesignTwo-step, two-sample multivariable Mendelian randomization (MVMR) study.SettingInternational genome-wide association study (GWAS) consortia data.ExposureT2DM, blood pressure: systolic blood pressure (SBP), diastolic blood pressure (DBP); blood lipids: low-density lipoprotein (LDL), high-density lipoprotein (HDL), total cholesterol (TC), triglycerides (TG); thyroid function: hyperthyroidism, hypothyroidism; body mass index (BMI), waist-hip-ratio (WHR), and insulin sensitivity. Main outcomesCVD including coronary heart disease (CHD), myocardial infarction (MI) and stroke.MethodsSummary-level data for exposures and main outcomes were extracted from GWAS consortia. We used two-sample MR to illustrate the causal effect of T2DM on CVD subtypes and regression-based MVMR to quantify the possible mediation effects of cardiometabolic factors on CVD.ResultsEach additional unit of log odds of T2DM increased 16% risk of CHD [OR: 1.16, 95% confidence interval (CI): 1.12-1.21], 15% risk of MI (OR: 1.15, 95%CI: 1.10-1.20), and 10% risk of stroke (OR: 1.10, 95%CI: 1.06-1.13). In mediation analysis, SBP, DBP and TG were found as main mediators, while the mediation effects of other cardiometabolic factors were not significant. The proportion of total effect of T2DM on CHD mediated by SBP, DBP and TG was 16% (95%CI: 8%-24%), 7% (95%CI: 1%-13%) and 10% (95%CI: 2%-18%), respectively. Mediation effect of SBP and DBP on MI and stroke, TG on MI was also prominent, while mediation effect of TG on stroke was not significant. Combined mediation effect of all three mediators accounted for 29%, 26% and 13% of total effect of T2DM on CHD, MI and stroke, respectively.ConclusionSBP, DBP and TG mediate a substantial proportion of the causal effect of T2DM on CVD and thus interventions on these factors might reduce considerable excess risk of CVD among T2DM patients.


2020 ◽  
Vol 41 (35) ◽  
pp. 3314-3322 ◽  
Author(s):  
Marina Cecelja ◽  
Louise Keehn ◽  
Li Ye ◽  
Tim D Spector ◽  
Alun D Hughes ◽  
...  

Abstract Aims Haemodynamic determinants of blood pressure (BP) include cardiac output (CO), systemic vascular resistance (SVR), and arterial stiffness. We investigated the heritability of these phenotypes, their association with BP-related single-nucleotide polymorphisms (SNPs), and the causal association between BP and arterial stiffness. Methods and results We assessed BP, central BP components, and haemodynamic properties (during a single visit) including CO, SVR, and pulse wave velocity (PWV, measure of arterial stiffness) in 3531 (1934 monozygotic, 1586 dizygotic) female TwinsUK participants. Heritability was estimated using structural equation modelling. Association with 984 BP-associated SNP was examined using least absolute shrinkage and selection operator (LASSO) and generalized estimating equation regression. One and two-sample Mendelian randomization (MR) was used to estimate the causal direction between BP and arterial stiffness including data on 436 419 UK Biobank participants. We found high heritability for systolic and pulsatile components of BP (>50%) and PWV (65%) with overlapping genes accounting for >50% of their observed correlation. Environmental factors explained most of the variability of CO and SVR (>80%). Regression identified SNPs (n = 5) known to be associated with BP to also be associated with PWV. One-sample MR showed evidence of bi-directional causal association between BP and PWV in TwinsUK participants. Two-sample MR, confirmed a bi-directional causal effect of PWV on BP (inverse variance weighted (IVW) beta = 0.11, P < 0.02) and BP on arterial stiffness (IVW beta = 0.004, P < 0.0001). Conclusion The genetic basis of BP is mediated not only by genes regulating BP but also by genes that influence arterial stiffness. Mendelian randomization indicates a bi-directional causal association between BP and arterial stiffness.


2019 ◽  
Author(s):  
Christopher N Foley ◽  
Paul D W Kirk ◽  
Stephen Burgess

AbstractMotivationMendelian randomization is an epidemiological technique that uses genetic variants as instrumental variables to estimate the causal effect of a risk factor on an outcome. We consider a scenario in which causal estimates based on each variant in turn differ more strongly than expected by chance alone, but the variants can be divided into distinct clusters, such that all variants in the cluster have similar causal estimates. This scenario is likely to occur when there are several distinct causal mechanisms by which a risk factor influences an outcome with different magnitudes of causal effect. We have developed an algorithm MR-Clust that finds such clusters of variants, and so can identify variants that reflect distinct causal mechanisms. Two features of our clustering algorithm are that it accounts for uncertainty in the causal estimates, and it includes ‘null’ and ‘junk’ clusters, to provide protection against the detection of spurious clusters.ResultsOur algorithm correctly detected the number of clusters in a simulation analysis, outperforming the popular Mclust method. In an applied example considering the effect of blood pressure on coronary artery disease risk, the method detected four clusters of genetic variants. A hypothesis-free search suggested that variants in the cluster with a negative effect of blood pressure on coronary artery disease risk were more strongly related to trunk fat percentage and other adiposity measures than variants not in this cluster.Availability and ImplementationMR-Clust can be downloaded from https://github.com/cnfoley/[email protected] or [email protected] InformationSupplementary Material is included in the submission.


2016 ◽  
Author(s):  
Hans van Kippersluis ◽  
Cornelius A Rietveld

AbstractBackgroundThe potential of Mendelian Randomization studies is rapidly expanding due to (i) the growing power of GWAS meta-analyses to detect genetic variants associated with several exposures, and (ii) the increasing availability of these genetic variants in large-scale surveys. However, without a proper biological understanding of the pleiotropic working of genetic variants, a fundamental assumption of Mendelian Randomization (the exclusion restriction) can always be contested.MethodsWe build upon and synthesize recent advances in the econometric literature on instrumental variables (IV) estimation that test and relax the exclusion restriction. Our Pleiotropy-robust Mendelian Randomization (PRMR) method first estimates the degree of pleiotropy, and in turn corrects for it. If a sample exists for which the genetic variants do not affect the exposure, and pleiotropic effects are homogenous, PRMR obtains unbiased estimates of causal effects in case of pleiotropy.ResultsSimulations show that existing MR methods produce biased estimators for realistic forms of pleiotropy. Under the aforementioned assumptions, PRMR produces unbiased estimators. We illustrate the practical use of PRMR by estimating the causal effect of (i) cigarettes smoked per day on Body Mass Index (BMI); (ii) prostate cancer on self-reported health, and (iii) educational attainment on BMI in the UK Biobank data.ConclusionsPRMR allows for instrumental variables that violate the exclusion restriction due to pleiotropy, and corrects for pleiotropy in the estimation of the causal effect. If the degree of pleiotropy is unknown, PRMR can still be used as a sensitivity analysis.Key messagesIf genetic variants have pleiotropic effects, causal estimates of Mendelian Randomization studies will be biased.Pleiotropy-robust Mendelian Randomization (PRMR) produces unbiased causal estimates in case (i) a subsample can be identified for which the genetic variants do not affect the exposure, and (ii) pleiotropic effects are homogenous.If such a subsample does not exist, PRMR can still routinely be reported as a sensitivity analysis in any MR analysis.If pleiotropic effects are not homogenous, PRMR can be used as an informal test to gauge the exclusion restriction.


Author(s):  
Christopher N Foley ◽  
Amy M Mason ◽  
Paul D W Kirk ◽  
Stephen Burgess

Abstract Motivation Mendelian randomization is an epidemiological technique that uses genetic variants as instrumental variables to estimate the causal effect of a risk factor on an outcome. We consider a scenario in which causal estimates based on each variant in turn differ more strongly than expected by chance alone, but the variants can be divided into distinct clusters, such that all variants in the cluster have similar causal estimates. This scenario is likely to occur when there are several distinct causal mechanisms by which a risk factor influences an outcome with different magnitudes of causal effect. We have developed an algorithm MR-Clust that finds such clusters of variants, and so can identify variants that reflect distinct causal mechanisms. Two features of our clustering algorithm are that it accounts for differential uncertainty in the causal estimates, and it includes ‘null’ and ‘junk’ clusters, to provide protection against the detection of spurious clusters. Results Our algorithm correctly detected the number of clusters in a simulation analysis, outperforming methods that either do not account for uncertainty or do not include null and junk clusters. In an applied example considering the effect of blood pressure on coronary artery disease risk, the method detected four clusters of genetic variants. A post hoc hypothesis-generating search suggested that variants in the cluster with a negative effect of blood pressure on coronary artery disease risk were more strongly related to trunk fat percentage and other adiposity measures than variants not in this cluster. Availability and implementation MR-Clust can be downloaded from https://github.com/cnfoley/mrclust. Supplementary information Supplementary data are available at Bioinformatics online.


Sign in / Sign up

Export Citation Format

Share Document