Abstract:Background The success of a Mendelian randomization (MR) study critically depends on the validity of the assumptions underlying MR. We focus on detecting heterogeneity (also known as horizontal pleiotropy) in two‐sample summary‐data MR. A popular approach is to apply Cochran's Q statistic method, developed for meta‐analysis. However, Cochran's Q statistic, including its modifications, is known to lack power when its degrees of freedom are large. Furthermore, there is no theoretical justification for the claimed null distribution of the minimum of the modified Cochran's Q statistic with exact weighting (Qmin ), although it seems to perform well in simulation studies. Method The principle of our proposed method is straightforward: if a set of variables are valid instruments, then any linear combination of these variables is still a valid instrument. Specifically, this principle holds when these linear combinations are formed using eigenvectors derived from a variance matrix. Each linear combination follows a known normal distribution from which a p value can be calculated. We use the minimum p value for these eigenvector‐based linear combinations as the test statistic. Additionally, we explore a modification of the modified Cochran's Q statistic by replacing the weighting matrix with a truncated singular value decomposition. Results Extensive simulation studies reveal that the proposed methods outperform Cochran's Q statistic, including those with modified weights, and MR‐PRESSO, another popular method for detecting heterogeneity, in cases where the number of instruments is not large or the Wald ratios take two values. We also demonstrate these methods using empirical examples. Furthermore, we show that Qmin does not follow, but is dominated by, the claimed null chi‐square distribution. The proposed methods are implemented in an R package iGasso. Conclusions Dimension reduction techniques are useful for generating powerful tests of heterogeneity in MR.

Random effect based tests for multinomial logistic regression in genetic association studies

Random-effect Based Test for Multinomial Logistic Regression: Choice of the Reference Level and Its Impact on the Testing

Genetic association tests in family samples for multi-category phenotypes

Regression-based Approach for Testing the Association Between Multi-Region Haplotype Configuration and Complex Trait

Association Test Between Haplotypes and Longitudinal Traits in Complex Pedigrees.

Efficient penalized generalized linear mixed models for variable selection and genetic risk prediction in high-dimensional data

Single-Locus Genetic Association Analysis By Ordinal Tests

The Testing of Random Effect in Mixed Generalized Single Index Models

Maximal Conditional Chi-Square Importance in Random Forests

Testing Gene-Gene Interactions Based on a Neighborhood Perspective in Genome-wide Association Studies

Analysis of Case-Control Association Studies: SNPs, Imputation and Haplotypes

An iterative approach to detect pleiotropy and perform Mendelian Randomization analysis using GWAS summary statistics

Summary statistics-based association test for identifying the pleiotropic effects with set of genetic variants

The eigen higher criticism and eigen Berk-Jones tests for multiple trait association studies based on GWAS summary statistics.

Testing and correcting for weak and pleiotropic instruments in two‐sample multivariable Mendelian randomization

Powerful Test of Heterogeneity in Two‐Sample Summary‐Data Mendelian Randomization

Statistical Inference for Genetic Relatedness Based on High-Dimensional Logistic Regression

Mendelian randomization with fine-mapped genetic data: Choosing from large numbers of correlated instrumental variables

Linear Models for Analysis of Multiple Single Nucleotide Polymorphisms with Quantitative Traits in Unrelated Individuals

Integrating summarized data from multiple genetic variants in Mendelian randomization: bias and coverage properties of inverse-variance weighted methods

Almost exact Mendelian randomization