Yan Li, He Qianchuan, Verma Shiv P, Zhang Xu, Giel Ann-Sophie, Maj Carlo, Graz Kathryn, Naderi Elnaz, Chen Jianhong, Ali Mourad Wagdy, Gharahkhani Puya, Shu Xiang, Offit Kenneth, Shah Pari M, Gerdes Hans, Molena Daniela, Srivastava Amitabh, MacGregor Stuart, Palles Claire, Thieme René, Vieth Michael, Gockel Ines, Vaughan Thomas L, Schumacher Johannes, Buas Matthew F
Department of Biostatistics and Bioinformatics, Roswell Park Comprehensive Cancer Center, Buffalo, NY, USA.
Public Health Sciences Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA.
HGG Adv. 2025 Apr 10;6(2):100399. doi: 10.1016/j.xhgg.2025.100399. Epub 2025 Jan 3.
Inherited genetics represents an important contributor to risk of esophageal adenocarcinoma (EAC), and its precursor Barrett's esophagus (BE). Genome-wide association studies have identified ∼30 susceptibility variants for BE/EAC, yet genetic interactions remain unexamined. To address challenges in large-scale G×G scans, we combined knowledge-guided filtering and machine learning approaches, focusing on genes with (1) known/plausible links to BE/EAC pathogenesis (n = 493) or (2) prior evidence of biological interactions (n = 4,196). Approximately 75 × 10 SNP×SNP interactions were screened via hierarchical group lasso (glinternet) using BEACON GWAS data. The top ∼2,000 interactions retained in each scan were prioritized using p values from single logistic models. Identical scans were repeated among males only (78%), with two independent GWAS datasets used for replication. In overall and male-specific primary replications, 11 of 187 and 20 of 191 interactions satisfied p < 0.05, respectively. The strongest evidence for secondary replication was for rs17744726×rs3217992 among males, with consistent directionality across all cohorts (P = 2.19 × 10); rs3217992 "T" was associated with reduced risk only in individuals homozygous for rs17744726 "G." Rs3217992 maps to the CDKN2B 3' UTR and reportedly disrupts microRNA-mediated repression. Rs17744726 maps to an intronic enhancer region in BLK. Through in silico prioritization and experimental validation, we identified a nearby proxy variant (rs4841556) as a functional modulator of enhancer activity. Enhancer-gene mapping and eQTLs implicated BLK and FAM167A as targets. The first systematic G×G investigation in BE/EAC, this study uncovers differential risk associations for CDKN2B variation by BLK genotype, suggesting novel biological dependency between two risk loci encoding key mediators of tumor suppression and inflammation.
遗传因素是食管腺癌(EAC)及其癌前病变巴雷特食管(BE)发病风险的重要影响因素。全基因组关联研究已经确定了约30个与BE/EAC相关的易感变异,但基因间的相互作用仍未得到研究。为应对大规模基因-基因(G×G)扫描中的挑战,我们结合了知识引导的筛选和机器学习方法,重点关注以下两类基因:(1)与BE/EAC发病机制有已知/合理联系的基因(n = 493),或(2)有生物相互作用先前证据的基因(n = 4196)。使用BEACON全基因组关联研究(GWAS)数据,通过分层组套索(glinternet)筛选了约75×10个单核苷酸多态性(SNP)×SNP相互作用。每次扫描中保留的约前2000个相互作用,使用单逻辑模型的p值进行优先级排序。仅在男性中重复相同的扫描(78%),使用两个独立的GWAS数据集进行验证。在总体和男性特异性的初次验证中,187个相互作用中有11个、191个相互作用中有20个分别满足p < 0.05。二次验证的最有力证据是男性中的rs17744726×rs3217992,在所有队列中具有一致的方向性(P = 2.19×10);仅在rs17744726为“G”纯合子的个体中,rs3217992的“T”与风险降低相关。rs3217992定位于CDKN2B的3'非翻译区(UTR),据报道会破坏微小RNA介导的抑制作用。rs17744726定位于BLK的内含子增强子区域。通过计算机模拟优先级排序和实验验证,我们确定了一个附近的替代变异(rs4841556)作为增强子活性的功能调节剂。增强子-基因图谱和表达数量性状基因座(eQTL)表明BLK和FAM167A是靶点。本研究是BE/EAC领域的首次系统性G×G研究,揭示了CDKN2B变异因BLK基因型而异的风险关联,提示了两个编码肿瘤抑制和炎症关键介质的风险位点之间存在新的生物学依赖性。