ABSTRACT
Multivariable Mendelian randomization (MVMR) is an instrumental variable technique that generalizes the MR framework for multiple exposures. Framed as a linear regression problem, it is subject to the pitfall of multi-collinearity. The bias and efficiency of MVMR estimates thus depends on the correlation of exposures. Dimensionality reduction techniques such as principal component analysis (PCA) provide transformations of all the included variables that are effectively uncorrelated. We propose the use of sparse PCA (sPCA) algorithms that create principal components of subsets of the exposures and may provide more interpretable and reliable MR estimates. The approach consists of three steps. We first apply a sparse dimension reduction method and transform the variant-exposure summary statistics to principal components. We then choose a subset of the principal components based on data-driven cutoffs, and estimate their strength as instruments with an adjusted F-statistic. Finally, we perform MR with these transformed exposures. This pipeline is demonstrated in a simulation study of highly correlated exposures and an applied example using summary data from a genome-wide association study of 118 highly correlated lipid metabolites. As a positive control, we tested the causal associations of the transformed exposures on CHD. Compared to the conventional inverse-variance weighted MVMR method and a weak-instrument robust MVMR method (MR GRAPPLE), sparse component analysis achieved a superior balance of sparsity and biologically insightful grouping of the lipid traits.
Key Messages
In multivariable MR, investigation of multiple highly correlated exposures can hinder the efficiency of the estimators and mask true associations.
Dimensionality reduction approaches such as principal component analysis (PCA) appear to be effective in summarising the variant-exposure summary statistics data in an example of correlated metabolite data.
Sparse PCA approaches have the additional benefit of providing interpretable PCs as only a few exposures contribute to each. This benefit is shown in simulation studies where there is a gain in accuracy over PCA.
In a positive control analysis, the sparse PCs were representing biologically meaningful groups of metabolites (VLDL, LDL, HDL) and in general were associated with coronary heart disease in an anticipated mannert.
Competing Interest Statement
Dr D.G. is a part-time employee of Novo Nordisk. The other authors have no relevant conflict of interest to declare.
Funding Statement
VK and Pr. JB received funding from the Expanding Excellence in England (E3) fund.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Not Applicable
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
Not applicable
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Not Applicable
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Not Applicable
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Not Applicable
Footnotes
↵* e-mail: vaskarageorg{at}hotmail.com
Data Availability
This is a computational study and no new data were generated. The code for the analyses is publicly available.
