Abstract:Causal discovery is a powerful tool to disclose underlying structures by analyzing purely observational data. Genetic variants can provide useful complementary information for structure learning. Here, we propose a novel algorithm MRSL (Mendelian Randomization (MR)-based Structure Learning algorithm), which combines the graph theory with univariable and multivariable MR to learn the true structure using only GWAS summary statistics. Specifically, MRSL also utilizes topological sorting to improve the precision of structure learning and provides three adjusting categories for multivariable MR. Results of simulation reveal that MRSL has up to two-fold higher F1 score than other eight competitive methods. Additionally, the computing time of MRSL is 100 times faster than other methods. Furthermore, we apply MRSL to 26 biomarkers and 44 ICD10-defined diseases from UK Biobank. The results cover most of expected causal links which have biological interpretations and several new links supported by clinical case reports or previous observational literatures. ### Competing Interest Statement The authors have declared no competing interest. ### Funding Statement This study was funded by the National Natural Science Foundation of China (Grant 82173625,82003557) and the Shandong Provincial Key Research and Development project (2018CXGC1210). ### Author Declarations I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained. Yes The details of the IRB/oversight body that provided approval or exemption for the research described are given below: The study used (or will use) ONLY openly available human data that were originally located at: http://www.nealelab.is/uk-biobank. 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. Yes 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). Yes 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. Yes All data produced are available online at http://www.nealelab.is/uk-biobank. <http://www.nealelab.is/uk-biobank>

Mrmlm V4.0: an R Platform for Multi-locus Genome-wide Association Studies

Mrmlm V4.0.2: an R Platform for Multi-locus Genome-wide Association Studies.

BLUPmrMLM: A Fast mrMLM Algorithm in Genome-wide Association Studies

A fast mrMLM algorithm for multi-locus genome-wide association studies

Improving Power and Accuracy of Genome-Wide Association Studies Via a Multi-Locus Mixed Linear Model Methodology

A SUPER powerful method for genome wide association study.

Methodological Implementation of Mixed Linear Models in Multi-Locus Genome-Wide Association Studies.

The Improved FASTmrEMMA and GCIM Algorithms for Genome-Wide Association and Linkage Studies in Large Mapping Populations

Editorial: the Applications of New Multi-Locus GWAS Methodologies in the Genetic Dissection of Complex Traits.

Iterative Usage of Fixed and Random Effect Models for Powerful and Efficient Genome-Wide Association Studies

The MR-Base platform supports systematic causal inference across the human phenome

A Fast Multi-Locus Ridge Regression Algorithm for High-Dimensional Genome-Wide Association Studies

A fast multi-locus random-SNP-effect EMMA for genome-wide association studies

Principle,Optimization and Application of Mixed Models in Genome-Wide Association Study

Iterative Sure Independence Screening Em-Bayesian Lasso Algorithm For Multi-Locus Genome-Wide Association Studies

Powerful Identification of Large Quantitative Trait Loci Using Genome-Wide R/glmnet-Based Regression

Efficient molecular mendelian randomization screens with LaScaMolMR.jl

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

FastBiCmrMLM: a fast and powerful compressed variance component mixed logistic model for big genomic case-control genome-wide association study

A Multi-Locus Linear Mixed Model Methodology for Detecting Small-Effect QTLs for Quantitative Traits in MAGIC, NAM, and ROAM Populations.

MRSL: A phenome-wide causal discovery algorithm based on GWAS summary data