We first conducted univariable MR analyses for each lipid-related trait. For this, we harmonised SNPs identified from our GWASs of lipoprotein lipid traits in the UKBB to those SNPs available in CARDIoGRAMplusC4D by either matching the SNP directly or by selecting proxy SNPs in high LD (r2 > 0.8). This led to a small drop in the number of SNPs being available for MR, with a median of 93% SNPs identified in GWASs available for MR (the numbers available for each trait are provided in Table 1 ). We used the inverse variance weighted approach, which, in brief, takes the form of a linear regression of the SNP–outcome association regressed on the SNP–exposure association weighted by the inverse of the square of the standard error of the SNP–outcome association, with the intercept constrained at the origin.
We next conducted multivariable MR, which is a statistical approach that allows for the association of SNPs with multiple phenotypes to be incorporated into the analysis, permitting an estimation of the direct effect of each phenotype on the outcome (i.e., an effect that is not mediated by any other factor in the model [28 (link)]); seeS1 Fig for further details. In this manuscript, we use the term ‘adjusted’ in the context of multivariable MR to mean ‘direct’ effects, i.e., the effect of a lipid trait on CHD whilst accounting for either mediation or confounding by another trait included in the model. For the multivariable MR analyses, we fitted a model with apolipoprotein B, LDL cholesterol, and triglycerides to identify which one or more traits appeared to be responsible for the effect of ‘atherogenic’ lipid-related traits on risk of CHD. We then took the atherogenic trait(s) that retained an effect on CHD in the multivariable MR model forward and further adjusted for apolipoprotein A-I and HDL cholesterol to assess the potential causal roles of HDL-related phenotypes in the development of CHD. In the setting of multivariable MR, we included all GWAS-associated SNPs for all traits in the model. This meant that there were differing numbers of SNPs in the 2 multivariable models tested.
We characterised instrument strengths in both the univariable and multivariable MR settings as follows: for the univariable estimates, we generated the mean F-statistic, using the approximation described by Bowden and colleagues [44 (link)]. For the multivariable estimate, we generated the conditional F-statistic [28 (link),45 ]. Further details are provided inS1 Text .
We next conducted multivariable MR, which is a statistical approach that allows for the association of SNPs with multiple phenotypes to be incorporated into the analysis, permitting an estimation of the direct effect of each phenotype on the outcome (i.e., an effect that is not mediated by any other factor in the model [28 (link)]); see
We characterised instrument strengths in both the univariable and multivariable MR settings as follows: for the univariable estimates, we generated the mean F-statistic, using the approximation described by Bowden and colleagues [44 (link)]. For the multivariable estimate, we generated the conditional F-statistic [28 (link),45 ]. Further details are provided in
Full text: Click here