Abstract
BACKGROUND: Numerous observational studies have investigated on the correlation of whole, semi-skimmed, and skimmed milk with coronary artery disease (CAD) and myocardial infarction (MI) risk; However, no consensus has been reached and evidence on any causal links between these exposures and outcomes remains unclear. This study aimed to conduct univariate and multivariate Mendelian randomization (MR) analyses, using publicly released genome-wide association study summary statistics (GWAS) from the IEU GWAS database, to ascertain the causal association of milk with various fat content with CAD and MI risk. METHODS: For the exposure data, 29, 15, and 30 single-nucleotide polymorphisms for whole milk, semi-skimmed milk, and skimmed milk, respectively, obtained from 360,806 Europeans, were used as instrumental variables. CAD and MI comprised 141,217 and 395,795 samples, respectively. We used inverse variance weighted (IVW), weighted median, MR-Egger regression, and MR Pleiotropy Residual Sum and Outlier analyses to determine whether pleiotropy and heterogeneity could skew the MR results. Sensitivity tests were conducted to verify the robustness of the results. RESULTS: After adjusting for false discovery rates (FDR), we discovered proof that skimmed milk intake is a genetically predicted risk factor for CAD (odds ratio [OR] = 5.302; 95% confidence interval [CI] 2.261-12.432; P < 0.001; FDR-corrected P < 0.001) and MI (OR = 2.287; 95% CI 1.218-4.300; P = 0.010; FDR-corrected P = 0.009). Most sensitivity assessments yielded valid results. Multivariable MR for CAD and MI produced results consistent with those obtained using the IVW method. There was no causal relationship between whole or semi-skimmed milk, and CAD or MI. CONCLUSION: Our findings indicate that the consumption of skimmed milk may increase the risk of CAD and MI. This evidence may help inform dietary recommendations for preventing cardiovascular disease. Further studies are required to elucidate the underlying mechanisms.