Several sensitivity analyses were used to check and correct for the presence of pleiotropy in the causal estimates. Cochran’s Q was computed to quantify heterogeneity across the individual causal effects, with a P-value ≤ 0.05 indicating the presence of pleiotropy, and that consequently, a random effects IVW MR analysis should be used43 (link),48 (link). We also assessed the potential presence of horizontal pleiotropy using MR-Egger regression based on its intercept term, where deviation from zero denotes the presence of directional pleiotropy. Additionally, the slope of the MR-Egger regression provides valid MR estimates in the presence of horizontal pleiotropy when the pleiotropic effects of the genetic variants are independent from the genetic associations with the exposure49 (link),50 (link). We also computed OR estimates using the complementary weighted-median method that can give valid MR estimates under the presence of horizontal pleiotropy when up to 50% of the included instruments are invalid44 (link). The presence of pleiotropy was also assessed using the MR-PRESSO. In this, outlying SNPs are excluded from the accelerometer-measured physical activity instrument and the effect estimates are reassessed51 (link). For all of the aforementioned sensitivity analyses to identify possible pleiotropy, we considered the estimates from the extended 10 SNP instrument as the primary results due to unstable estimates from the 5 SNP instrument. A leave-one-SNP out analysis was also conducted to assess the influence of individual variants on the observed associations. We also examined the selected genetic instruments and their proxies (r2 > 0.8) and their associations with secondary phenotypes (P-value < 5 × 10−8) in Phenoscanner (http://www.phenoscanner.medschl.cam.ac.uk/) and GWAS catalog (date checked April 2019).
For the extended 10 SNP instrument, we also conducted multivariable MR analyses to adjust for potential pleiotropy due to BMI because the initial GWAS on physical activity reported several strong associations (P-value < 10−5) between the identified SNPs and BMI52 (link). The new estimates correspond to the direct causal effect of physical activity with the BMI being fixed. The genetic data on BMI were obtained from a GWAS study published by The Genetic Investigation of ANthropometric Traits (GIANT) consortium53 (link) (Supplementary Table 9). Additionally, for the extended 10 SNP instrument, we also conducted analyses with adiposity-related SNPs (i.e. those previously associated with BMI, waist circumference, weight, or body/trunk fat percentage in GWAS studies at P-value < 10−8) excluded (n = 5; rs34517439, rs6775319, rs11012732, rs1550435, rs59499656). Finally, we conducted two-sample MR analyses using BMI adjusted GWAS estimates for the 5 SNP accelerometer-measured physical activity instrument11 (link). However, the MR results using the BMI adjusted GWAS estimates should be interpreted cautiously due to the potential for collider bias11 (link).
All the analyses were conducted using the MendelianRandomisation54 (link) and TwoSampleMR55 (link) packages, and the R programming language.
Free full text: Click here