A Two-sample Mendelian Randomization Study of Oily Nut Intake and Primary Lung Cancer
-
摘要:
目的 以两样本孟德尔随机化方法探讨油性坚果摄取量与原发性肺癌的因果关系。 方法 从全基因组关联研究openGWAS公开数据库中筛选出与油性坚果摄入量相关的(P < 5×10-7)单核苷酸多态性位点作为工具变量,以逆方差加权法为该研究评估因果效应的主要分析方法,并采用MR-Egger回归和加权中位数法用于补充说明,以OR值及95%CI评价油性坚果摄入量与原发性肺癌之间是否存在因果关系。采用mr_pleiotropy_test法及MR-PRESSO法检验水平多效性,leave-one-out留一法进行敏感性检验。 结果 从ukb-b-16576数据集共提取100个SNP作为工具变量,MR-Egger回归、加权中位数和逆方差加权法计算得OR值和95%可信区间分别[0.63,95%CI(0.11~3.65),P = 0.611],[0.34,95%CI(0.19~0.61),P = 0.001],[0.32,95%CI(0.21~0.48),P = 0.001)],所有b值结果方向均小于0。mr_pleiotropy_test检验显示MR-Egger截距为-7.66E-03,接近于0,P = 0.441,同样MR-PRESSO包GlobalTest检验P = 0.162,所选的工具变量不存在水平多效性,经过leave-one-out敏感性分析,未发现任何SNP位点对结果产生显著影响,表明结果具有较高的稳定性和可靠性。 结论 两样本孟德尔随机化分析表明油性坚果摄入量与原发性肺癌的发病风险呈负向因果关系,但需要更大样本量和不同种族背景的人的数据来进一步验证。 Abstract:Objective To investigate the potential causal relationship between oily nut consumption and primary lung cancer by two-sample Mendelian randomization method. Methods Single nucleotide polymorp hisms (P < 5×10-7) associated with oily nut intake were screened from the openGWAS public database of genome - wide association study as instrumental variables. The inverse variance weighting method was the main analysis method to evaluate the causal effect of this study, and MR-Egger regression and weighted median method were used to supplement the explanation. OR value and 95%CI were used to evaluate whether there was a causal relationship between oily nut intake and primary lung cancer. The mr_pleiotropy_test and MR-PRESSO were used to test the level pleiotropy, and the leave-one-out method was used to test the sensitivity. Results A total of 100 SNPs were extracted from the ukb-b-16576 dataset as instrumental variables. The MR-Egger regression, weighted median, and inverse variance weighting methods yielded OR values and 95% confidence intervals [0.63, 95%CI (0.11-3.65), P = 0.611], [0.34, 95%CI (0.19-0.61), P = 0.001], [0.32, 95%CI (0.21-0.48), P = 0.001], rerespectively. The mr_pleiotropy_test showed that the MR-Egger intercept was close to 0 (-7.66E-03) with a p-value of 0.441, indicating no significant horizontal pleiotropy. Similarly, the MR-PRESSO package GlobalTest yielded a p-value of 0.162, further confirming the absence of horizontal pleiotropy among the selected instrumental variables.After conducting the leave-one-out sensitivity analysis, no significant influence on the results was observed from any specific SNP locus, indicating that the results are highly stable and reliable. Conclusions The two-sample Mendelian randomization analysis revealed a negative causal relationship between oily nut intake and the risk of primary lung cancer. However, further verification is needed with larger sample sizes and data from individuals of different ethnic backgrounds to strengthen the findings. -
Key words:
- Oily nuts /
- Primary lung cancer /
- Mendel randomization /
- Causal inference
-
表 1 工具变量SNP信息表
Table 1. Instrumental variable SNP information table
SNP chr EA OA EAF B S R2 F P rs3101339 1 C A 0.6 1.43E−02 1.71E−03 3.36E−03 69.91 6.20E−17 rs10740991 10 C G 0.72 1.67E−02 1.86E−03 3.35E−03 81.23 2.00E−19 rs746868 6 G C 0.61 −1.29E−02 1.71E−03 3.03E−03 56.65 5.20E−14 rs4269101 3 G T 0.72 −1.38E−02 1.86E−03 2.75E−03 55.17 1.10E−13 rs9385269 6 T C 0.52 1.21E−02 1.68E−03 2.65E−03 51.48 7.20E−13 rs4800488 18 A C 0.49 1.20E−02 1.67E−03 2.49E−03 51.36 7.70E−13 rs4149513 2 A G 0.49 1.17E−02 1.67E−03 2.46E−03 49.25 2.20E−12 rs1582322 16 G A 0.6 9.94E−03 1.72E−03 2.32E−03 33.58 6.80E−09 。。。。。。。。。。。。。。。。。。。。。。。。。 。。。。。。。。。。。。。。。。。。。。。。。。。 rs7829800 8 G A 0.67 −1.04E−02 1.79E−03 2.30E−03 34.16 5.10E−09 rs261809 1 G A 0.54 −9.63E−03 1.68E−03 2.17E−03 32.88 9.80E−09 rs6533183 4 T C 0.66 −9.54E−03 1.76E−03 2.16E−03 29.42 5.80E−08 rs4140799 14 A G 0.53 9.46E−03 1.68E−03 2.11E−03 31.74 1.80E−08 rs10087278 8 A G 0.59 8.99E−03 1.70E−03 2.11E−03 27.92 1.30E−07 rs10129747 14 G A 0.53 9.36E−03 1.68E−03 2.07E−03 30.98 2.60E−08 rs3757323 6 C T 0.63 8.90E−03 1.74E−03 2.06E−03 26.28 2.90E−07 rs11073575 15 C T 0.58 −8.90E−03 1.70E−03 2.05E−03 27.28 1.80E−07 表 2 5种MR分析方法OR估计值及95%CI
Table 2. OR estimates and 95% CI of five MR analysis methods
method b se P OR 95%CI Q值 Q_pval MR Egger −0 .47 0.89 0.611 0.63(0.11,3.65) 99.32 0.15 Weighted median −1.07 0.29 0.001 0.34(0.19,0.61) 100.01 0.16 Inverse variance weighted −1.14 0.21 0.001 0.32(0.21,0.48) Simple mode −0.99 0.76 0.192 0.37(0.08,1.64) Weighted mode −0.99 0.73 0.173 0.37(0.08,1.51) -
[1] Sung H,Ferlay J,Siegel R L,et al. Global cancer statistics 2020: Globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J]. CA Cancer J Clin,2021,71(3):209-249. doi: 10.3322/caac.21660 [2] Xia C,Dong X,Li H,et al. Cancer statistics in China and United States,2022: Profiles,trends,anddeterminants[J]. Chin Med J (Engl),2022,135(5):584-590. doi: 10.1097/CM9.0000000000002108 [3] Chen S,Cao Z,Prettner K,et al. Estimates and projections of the global economic cost of 29 cancers in 204 countries and territories from 2020 to 2050[J]. JAMA Oncol,2023,9(4):465-472. doi: 10.1001/jamaoncol.2022.7826 [4] 开比努尔·再比布力,姑丽切克然·艾斯克,韩加. 坚果营养成分及保健功效的研究[J]. 粮食与食品工业,2022,29(3):32-36. doi: 10.3969/j.issn.1672-5026.2022.03.009 [5] Chang C,Xinyan G,Yan H,et al. Association between nut consumption and cancer risk: A meta-analysis[J]. Nutrition and cancer,2022,75(1):11-13. [6] Nieuwenhuis L,van den Brandt P A. Nut and peanut butter consumption and the risk of lung cancer and its subtypes: A prospective cohort study[J]. Lung Cancer,2019,128(1):57-66. [7] Hemani G,Zheng J,Elsworth B,et al. The MR-Base platform supports systematic causal inferenceacross the human phenome[J]. Elife,2018,7(1):e34408. [8] Rasooly D,Patel C J. Conducting a reproducible mendelian randomization analysis using the R analytic statistical environment.[J]. Current Protocols in Human Genetics,2019,101(1):82. doi: 10.1002/cphg.82 [9] Paraskevas A,Bellou E ,Keramida E,et al. Bronchial asthma and lung cancer[J]. Current Respira tory Medicine Reviews,2021,17(2): 78-82. [10] Forder A,Zhuang R,Souza V,et al. Mechanisms contributing to the comorbidity of COPD and lung cancer[J]. Int J Mol Sci,2023,24(3):2859-2861. doi: 10.3390/ijms24032859 [11] Levin M G, Judy R, Gill D, et al. Genetics of height and risk of atrial fibrillation: A Mendelian randomization study[J]. PLoS Med,2020,17(10):e1003288. [12] Fu R,Kim S J. Inferring causality from observational studies: the role of instrumental variable ana lysis[J]. Kidney Int,2021,99(6):1303-1308. doi: 10.1016/j.kint.2021.03.018 [13] Wang X,Jiang Y,Zhang N R,et al. Sensitivity analysis and power for instrumental variable studi es[J]. Biometrics,2018,74(4):1150-1160. doi: 10.1111/biom.12873 [14] 周启武,唐兴莹,李志宏,等. 5种市售坚果类产品部分营养成分及重金属元素分析[J]. 安徽农业科学,2022,50(19):193-197. [15] Luu H N,Cai H,Murff H J,et al. A prospective study of dietary polyunsaturated fatty acids int ke and lung cancer risk[J]. Int J Cancer,2018,143(9):2225-2237. doi: 10.1002/ijc.31608 [16] Yang J J,Yu D,Takata Y,et al. Dietary Fat intake and lung cancer risk: A pooled analysis[J]. J Clin Oncol,2017,35(26):3055-3064. doi: 10.1200/JCO.2017.73.3329 [17] Bai Y,Wang G,Fu W,et al. Circulating essential metals and lung cancer: risk assessment and p tential molecular effects[J]. Environ Int,2019,127(6):685-693. doi: 10.1016/j.envint.2019.04.021 [18] 李晓阳,吴志平,王梦馨,等. 表没食子儿茶素没食子酸酯抗癌分子机制及其应用的研究进展[J]. 草药,2019,50(13):3217-3229. [19] Zanoaga O,Braicu C,Jurj A,et al. Progress in research on the role of flavonoids in lung can cer[J]. Int J Mol Sci,2019,20(17):4291-4293. doi: 10.3390/ijms20174291 [20] Baer Benjamin R,Strawderman Robert L,Ertefaie Ashkan,et al. Discussion on Instrumental variable estimation of the causal hazard ratio[J]. Biometrics,2023,79(2):554-558. [21] Fei W,Dylan S,Nandita M. A general approach to evaluating the bias of 2-stage instrumental var iable estimators[J]. Statistics in Medicine,2018,37(12):1997-2015. doi: 10.1002/sim.7636