Methods
We used published large-scale genome-wide association data and applied univariate and multivariate Mendelian randomization methods. The main analysis model was inverse variance weighting, and the reliability of the results was tested using MR-Egger intercept analysis, Cochran's 