《Advanced Science》：Mapping Genetic Regulation of Transcription to Identify Functional Variants and Genes Associated with Pancreatic Cancer Risk
编辑推荐：
本期推荐一篇在胰腺癌遗传学与功能基因组学领域的重要研究。该工作通过对482例胰腺组织进行表达数量性状位点（eQTL）荟萃分析，构建了胰腺组织的全基因组调控图谱，并与全基因组关联研究（GWAS）数据整合，系统性地鉴定了与胰腺癌风险相关的功能性单核苷酸多态性（SNP）及其靶基因。研究重点验证了保护性变异rs11102484，揭示了其通过减弱ZNF263结合，削弱远端沉默子-启动子远程互作，从而上调肿瘤抑制基因ST7L表达的分子机制，并证明ST7L通过抑制AKT/β-catenin信号通路发挥抑癌功能。这项工作为解读胰腺癌GWAS发现提供了功能机制见解，并鉴定了新的潜在干预靶点。
胰腺导管腺癌是全球范围内最具致死性的恶性肿瘤之一，其五年生存率极低。尽管全基因组关联研究（GWAS）已鉴定出多个胰腺癌风险位点，但绝大多数位于非编码区，其生物学功能及调控的靶基因仍不明确，这限制了我们从遗传关联到致病机制的理解。为了系统解析胰腺癌风险变异的调控机制，本研究整合了大规模的表达数量性状位点（eQTL）数据和跨人群的GWAS荟萃分析结果。
胰腺组织的全基因组调控图谱绘制
研究人员首先对来自癌症基因组图谱（TCGA）的177例胰腺癌组织和基因型-组织表达（GTEx）项目的305例正常胰腺组织，共计482个样本，进行了顺式eQTL（cis-eQTL）荟萃分析。通过固定效应逆方差荟萃，显著提升了统计效能，最终在错误发现率（FDR）小于0.05的水平上，鉴定出1,123,483个显著的SNP-基因对，对应13,758个eGenes（受遗传调控的基因）。这些eGenes在功能上显著富集于DNA修复、脂肪酸代谢、mTORC1信号传导和p53通路等癌症相关通路。进一步分析显示，eGenes与胰腺癌中高频体细胞突变基因（如KRAS、TP53、SMAD4）高度重合，并且其表达水平与肿瘤免疫微环境中的免疫细胞浸润（如巨噬细胞、CD8+T细胞）显著相关，提示这些遗传调控基因在胰腺癌发生发展中具有广泛的生物学和临床相关性。
整合分析优先确定胰腺癌风险相关功能变异与基因
为了从GWAS信号中筛选出可能具有调控功能的风险变异，研究者将上述eQTL结果与一项包含5130例病例和5776例对照的跨人群GWAS荟萃分析数据（结合了中国人群和PanC4联盟数据）进行整合。通过筛选在GWAS中具有提示性关联（Pmeta< 1 × 10-5，且在两个人群中均P < 0.05）且与基因表达显著相关（eQTL FDR < 0.05）的SNP-基因对，最终锁定了82个潜在的功能性SNP，它们映射到15个候选靶基因。这些基因包括ETAA1、GSDMB、PDX1等已知与癌症相关进程（如DNA损伤应答、免疫调控）有关的基因，也包含了ST7L等有待深入研究的基因。这些变异位于1p13.2、13q12.2等多个基因组区域，为理解胰腺癌的遗传基础提供了新的候选分子。
保护性变异rs11102484的功能验证与机制解析
在众多候选变异中，位于1p13.2区域的rs11102484展现出最强的功能潜力证据。独立队列（569例病例和2691例对照）的验证及合并分析（总计5699例病例和8467例对照）均确认，其G等位基因是胰腺癌的保护性因素，可显著降低患病风险（合并比值比OR = 0.85，95%置信区间CI = 0.80–0.90，P = 4.83 × 10-8）。eQTL分析显示，rs11102484的G等位基因与抑癌基因ST7L的表达上调显著相关。
为阐明其调控机制，研究人员开展了一系列功能实验：
- 1.
表观遗传特征与调控活性：rs11102484位于染色质开放区域，并带有多种组蛋白修饰标记，提示其具有调控元件特性。双荧光素酶报告基因实验证明，该区域具有沉默子活性，且G等位基因比风险C等位基因的沉默作用更弱。
- 2.
远程染色质互作：染色体构象捕获（3C）实验和公共Hi-C数据分析均证实，rs11102484所在区域与112 kb外的ST7L启动子存在物理相互作用。这种互作具有等位基因特异性，在携带保护性GG基因型的细胞中互作频率最低，而在CC基因型细胞中最高，表明保护性等位基因削弱了沉默子对启动子的远程抑制。
- 3.
关键转录因子ZNF263的介导作用：生物信息学预测和实验验证（电泳迁移率变动分析EMSA、染色质免疫共沉淀ChIP-qPCR）发现，rs11102484位于转录抑制因子ZNF263的结合基序内。保护性G等位基因会削弱ZNF263与该位点的结合，进而减少ZNF263在ST7L启动子区的富集。功能上，敲低ZNF263可导致ST7L表达上调并抑制细胞增殖，模拟了保护性等位基因的效应。
靶基因ST7L的抑癌功能与下游信号通路
ST7L（抑制肿瘤发生性7样蛋白）在多种癌症中被报道为抑癌基因。本研究证实，ST7L在胰腺癌组织中表达显著下调。功能实验表明，在胰腺癌细胞系（BxPC-3和PANC-1）中敲低ST7L可促进细胞增殖和克隆形成，而过表达ST7L则产生相反的抑癌效果。
在机制层面，Western blot实验显示，敲低ST7L会激活AKT/β-catenin信号通路，表现为磷酸化AKT（p-AKT (Ser473)）和β-catenin蛋白水平的升高；而过表达ST7L则抑制该通路。这与既往研究中ST7L抑制AKT/GSK3β/β-catenin轴的结果一致。此外，蛋白质互作网络分析预测ST7L可能与AKT1、TAOK2、MAP2K3、RINT1等蛋白存在相互作用，暗示其可能参与更广泛的信号网络调控。
研究意义与总结
本研究通过整合大规模eQTL与GWAS数据，构建了胰腺组织的遗传调控图谱，并系统性地鉴定了与胰腺癌风险相关的功能性变异及靶基因。研究者不仅通过计算整合优先出候选名单，更以rs11102484-ST7L调控轴为范例，完成了从群体遗传关联、独立验证到分子机制阐释的完整功能研究闭环。研究揭示了一个完整的致病通路：保护性等位基因rs11102484(G)通过削弱转录抑制因子ZNF263的结合，减弱远端沉默子对ST7L启动子的抑制作用，从而上调ST7L表达；高表达的ST7L蛋白进而抑制促癌的AKT/β-catenin信号通路，最终抑制胰腺癌细胞生长，降低个体患病风险。
这项工作将非编码风险变异与特定基因功能和信号通路直接联系起来，显著推进了对胰腺癌GWAS发现的生物学解读，为未来开发基于遗传学的风险预测、早期检测策略乃至针对ST7L等相关通路的新型疗法提供了重要的科学依据。同时，研究所构建的eQTL资源及优先的候选基因列表，也为胰腺癌遗传学和生物学领域的后续研究提供了宝贵的数据基础和方向指引。
