S2 File -

Background
Multimorbidity has become an important health challenge in the aging population. Accumulated evidence has shown that multimorbidity has complex association patterns, but the further mechanisms underlying the association patterns are largely unknown.
Methods
Summary statistics of 14 conditions/diseases were available from the genome-wide association study (GWAS). Linkage disequilibrium score regression analysis (LDSC) was applied to estimate the genetic correlations. Pleiotropic SNPs between two genetically correlated traits were detected using pleiotropic analysis under the composite null hypothesis (PLACO). PLACO-identified SNPs were mapped to genes by Functional Mapping and Annotation of Genome-Wide Association Studies (FUMA), and gene set enrichment analysis and tissue differential expression were performed for the pleiotropic genes. Two-sample Mendelian randomization analyses assessed the bidirectional causality between conditions/diseases.
Results
LDSC analyses revealed the genetic correlations for 20 pairs based on different two-disease combinations of 14 conditions/diseases, and genetic correlations for 10 pairs were significant after Bonferroni adjustment (P<0.05/91 = 5.49E-04). Significant pleiotropic SNPs were detected for 11 pairs of correlated conditions/diseases. The corresponding pleiotropic genes were differentially expressed in the brain, nerves, heart, and blood vessels and enriched in gluconeogenesis and drug metabolism, biotransformation, and neurons. Comprehensive causal analyses showed strong causality between hypertension, stroke, and high cholesterol, which drive the development of multiple diseases.
Conclusions
This study highlighted the complex mechanisms underlying the association patterns that include the shared genetic components and causal effects among the 14 conditions/diseases. These findings have important implications for guiding the early diagnosis, management, and treatment of comorbidities.