3.1 Instrumental variables
For lipoprotein(a) as the exposure, IVs were chosen as SNPs linked to lipoprotein(a) (4 SNPs for CeD, 8 SNPs for CD, 7 SNPs for IBD, 52 SNPs for MS, 65 SNPs for Pso, 46 SNPs for RA, 57 SNPs for SLE, 77 SNPs for T1D, 8 SNPs for UC). For lipoprotein(a) as the outcome, IVs were chosen as SNPs linked to IMIDs (10 SNPs for CeD, 44 SNPs for CD, 40 SNPs for IBD,12 SNPs for MS, 6 SNPs for Pso,15 SNPs for RA, 9 SNPs for SLE, 21 SNPs for T1D, 29 SNPs for UC). The absence of weak instrument bias was indicated by all F-statistics > 10. The comprehensive details regarding the instrumental variables can be found in Table S1-S3.
3.2 Causal estimates of genetic susceptibility to lipoprotein(a) and IMIDs risk
The findings from the MR analysis exploring the causal association between lipoprotein(a) and nine IMIDs traits are depicted in Fig. 2. Lipoprotein(a) exhibited no causal association with CeD (OR = 0.797, 95% CI: 0.498 to 1.274), CD (OR = 1.231, 95% CI: 0.519 to 2.920), IBD (OR = 1.026, 95% CI: 0.804 to 1.309), MS (OR = 1.007, 95% CI: 0.822 to 1.232), Pso (OR = 1.059, 95% CI: 0.910 to 1.231), RA (OR = 0.936, 95% CI: 0.782 to 1.120), SLE (OR = 1.039, 95% CI: 0.801 to 1.349), T1D (OR = 1.065, 95% CI: 0.931 to 1.218), and UC (OR = 1.016, 95% CI: 0.707 to 1.461). This discovery aligns with the outcomes derived from alternative MR techniques, including MR Egger and weighted median.
There was noticeable heterogeneity observed in our instrumental variables for Lp(a) in relation to CD (Q P.val = 7.65E-11), SLE (Q P.val = 3.66E-5), MS (Q P.val = 0.0001), Pso (Q P.val = 0.0043) and T1D (Q P.val = 0.0051) as outcome (Table S4). We observed that all MR-Egger regression intercepts were not significantly different from zero, indicating no indication of horizontal pleiotropy between the Lp(a) instrumental variables and IMIDs (intercept p > 0.05), except for T1D where a marginal deviation was found (intercept p = 0.021) (Table S4). However, the MR-PRESSO analysis revealed significant horizontal pleiotropy in certain analysis. Nevertheless, the causal estimates of Lp(a) with MS, Pso, SLE, and T1D remained consistent even after conducting outlier-corrected analyses (Table S5).
In addition, the sensitivity analysis plots indicated that no individual SNP was expected to have a substantial impact on the causal relationship, thus affirming the reliability of our findings (Figure S1-4). Taken collectively, these findings provide compelling evidence supporting the absence of a causal association between Lp(a) and IMIDs.
3.3 Causal estimates of genetic susceptibility to IMIDs and lipoprotein(a) levels
Furthermore, conducting reverse studies investigating the association between exposure to the risk of 9 IBIDs and the outcome of Lp(a) levels, we found no significant association between CeD (OR = 0.998, 95% CI: 0.995 to 1.002), CD (OR = 1.000, 95% CI: 0.991 to 1.008), IBD (OR = 0.995, 95% CI: 0.987 to 1.003), MS (OR = 0.992, 95% CI: 0.983 to 1.001), Pso (OR = 1.003, 95% CI: 0.989 to 1.018), RA (OR = 1.003, 95% CI: 0.998 to 1.009), SLE (OR = 0.998, 95% CI: 0.993 to 1.004), T1D (OR = 0.998, 95% CI: 0.994 to 1.002), UC (OR = 1.004, 95% CI: 0.994 to 1.013) and Lp(a) in the IVW analysis results (Fig. 3). The outcomes obtained from each of the three MR techniques exhibited concurrence. There was noticeable diversity observed in our instrumental variables for CD (Q P.val = 1.60E-6), and UC (Q P.val = 0.0024 (Table S6). The presence of imbalanced horizontal pleiotropy was not indicated by the MR-Egger intercept, as it exhibited a central tendency around zero in all MR analyses (Table S6). Although in certain analyses, MR-PRESSO revealed the presence of substantial horizontal pleiotropy, the causal estimates of Lp(a) with UC and CD remained robust after outlier-corrected analyses (Table S7). Additionally, the sensitivity analysis plots, which employed a leave-one-out approach, indicated that the individual impact of each SNP on the causal association was not significant. This finding further strengthens our conclusions (Figure S5-8).
3.4 Multivariable MR
To account for potential pleiotropic pathways arising from the relationship between different lipid traits, we employed a multivariable Mendelian randomization model incorporating Lp(a), HDL-C, LDL-C, and TG as joint exposures for each IMIDs outcome. Following adjustment for HDL-C, LDL-C, and TG, genetically elevated Lp(a) showed no causal association with the onset of IMIDs, consistent with the findings of univariable MR analysis. However after adjusting for other lipid traits, genetically predicted HDL-C (ORMVMR = 0.80, 95% CI: 0.68–0.95; P = 0.011) and TG (ORMVMR = 0.80, 95% CI: 0.66–0.98; P = 0.033) were negatively associated with type 1 diabetes. The association between genetically predicted LDL and Pso remained marginally significant even after adjusting for multiple lipid traits. (ORMVMR = 0.80, 95% CI: 0.64–0.99; P = 0.045) (Table S8).
Incorporating a distinct panel of 75 and 42 SNPs that have demonstrated strong and separate associations with CD and UC, correspondingly, multivariable Mendelian randomization analysis provided compelling evidence suggesting no causal relationship between the genetic liability to these autoimmune traits and lipoprotein(a) levels (CD (OR = 1.00, 95% CI: 0.99–1.01, P = 0.997), UC (OR = 1.00, 95% CI: 0.99–1.01, P = 0.679)).