双酚A暴露与胰腺癌的关联：来自人群流行病学、多组学筛选和细胞验证的证据

《Ecotoxicology and Environmental Safety》：Linking bisphenol a exposure to pancreatic cancer: Evidence from population epidemiology, multi-omics screening, and cellular validation

　　

双酚A（BPA）是一种普遍存在的环境内分泌干扰物；其与胰腺癌的关联及潜在作用机制尚不清楚。本研究建立了一个整合人群流行病学、计算毒理学、机器学习、分子动力学模拟和实验验证的综合框架，进行了系统性研究。对NHANES 2003–2016数据（N=11,894名成年人）的分析显示，高尿BPA暴露（四分位数4对比四分位数1–3）与癌症死亡率显著增加相关（OR=1.43，95% CI：1.04–1.96，P=0.027；PAF=9.51%）。对597个BPA靶基因和3933个胰腺癌相关基因的交集分析确定了246个重叠候选基因，通过15种算法和175种组合的机器学习筛选，从中优先筛选出14个核心基因。基于这14个基因构建的诊断列线图在GSE15471和TCGA-PAAD+GTEx数据集中分别达到AUC为0.984和0.997。分子动力学模拟（200 ns × 3次重复）结合MM/GBSA结合自由能计算表明，AHR对BPA的结合亲和力（?57.55±7.20 kcal/mol）约为GPRC5A（?28.87±12.96）的两倍，与已知AHR配体相当。在功能上，BPA表现出细胞类型依赖的双相效应，抑制正常胰腺导管上皮细胞（hTERT-HPNE）活力，同时增强AsPC-1和MiaPaCa-2胰腺癌细胞的集落形成，而AHR敲低减弱了这种增殖优势。总之，这些发现提供了人群水平、计算、生物物理和功能证据，支持BPA与胰腺癌的关联，并确定AHR为潜在的关键分子中介。

**论文解读：双酚A暴露与胰腺癌的关联——来自人群流行病学、多组学筛选和细胞验证的证据**

**研究背景、存在问题及研究目的**
双酚A（BPA）是一种广泛使用的工业化合物，作为聚碳酸酯塑料和环氧树脂的主要单体，年产量超过150亿磅，存在于食品容器、水管、婴儿奶瓶等数百种消费品中，并通过水源、空气和皮肤接触等多种途径导致人类持续暴露。BPA已被确认为环境内分泌干扰物，其与多种恶性肿瘤的关联备受关注。胰腺作为关键内分泌器官，被认为是BPA作用的重要靶器官，但BPA与胰腺癌的关联及其潜在分子机制尚不明确。胰腺癌临床预后极差，总体5年生存率约5%，预计2030年将成为癌症相关死亡的第二大原因，多因早期症状隐匿导致多数患者确诊时已处于晚期，错失最佳治疗时机。因此，阐明BPA暴露在胰腺癌中的作用、评估其作为诊断因素的潜力并识别有意义的分子靶点，对改善治疗策略和早期检测具有重要意义。

**研究内容与主要结论**
研究人员构建了整合人群流行病学、计算毒理学、机器学习、分子动力学模拟和实验验证的多维研究框架，系统探究BPA与胰腺癌的关联及潜在机制。基于NHANES 2003–2016数据（11,894名成年人）的分析显示，高尿BPA暴露（四分位数4对比1–3）与癌症死亡率显著增加相关（OR=1.43，95% CI：1.04–1.96，P=0.027），归因危险度百分比（PAF）为9.51%，估计每年美国约409例超额胰腺癌死亡。通过BPA靶基因与胰腺癌相关基因的交集分析，识别出246个重叠基因，经15种机器学习算法和175种组合的共识筛选，确定14个核心基因（AR、NTRK3、ESR2、RAC1、CYP2E1、AHR、INHBA、CCR2、EPHA1、GPRC5A、CRYAB、ACACB、EPHX2、MGLL）。基于这些基因构建的诊断列线图在GSE15471和TCGA-PAAD+GTEx数据集中分别达到AUC 0.984和0.997。分子动力学模拟（200 ns×3次重复）结合MM/GBSA计算表明，芳烃受体（AHR）对BPA的结合亲和力（?57.55±7.20 kcal/mol）约为GPRC5A的两倍，与已知AHR配体（如TCDD）相当。功能实验显示，BPA抑制正常胰腺导管上皮细胞（hTERT-HPNE）活力，但促进胰腺癌细胞（AsPC-1和MiaPaCa-2）集落形成，且AHR敲低可减弱此增殖优势。该研究发表在《Ecotoxicology and Environmental Safety》。

**主要技术方法**
研究人员利用NHANES 2003–2016队列（N=11,894名成年人，来自美国国家健康与营养调查）进行流行病学分析；通过SwissTargetPrediction、ChEMBL、SEA、DrugBank、TargetNet数据库预测BPA靶基因，并从GeneCards、OMIM、CTD、GEO数据库获取胰腺癌相关基因；采用15种机器学习算法（包括Lasso、Ridge、弹性网络、逐步GLM、SVM、随机森林、XGBoost等）及175种组合进行核心基因筛选；运用分子动力学模拟（200 ns×3次重复）结合MM/GBSA计算评估结合自由能；使用CETSA、CCK-8、集落形成、活/死细胞染色、免疫荧光、RT-qPCR及慢病毒shRNA/siRNA敲低等技术进行细胞功能验证。

**研究结果**
**3.1 BPA暴露与癌症死亡率增加相关**
分析NHANES数据发现，尿BPA检出率91.8%，几何均值1.61 ng/mL。最高四分位数（Q4）参与者相比Q1–Q3，癌症死亡风险增加43%（OR=1.43，95% CI：1.04–1.96，P=0.027），并呈现显著剂量-反应趋势（P-trend=0.046）。亚组分析显示，男性（OR=1.22，P=0.047）和非西班牙裔黑人（OR=1.32，P=0.021）风险更高。

**3.2 BPA靶点与胰腺癌相关基因的交集**
从597个BPA潜在靶基因与3933个胰腺癌相关基因中鉴定出246个重叠基因。功能富集分析显示，这些基因显著富集于PI3K-Akt、MAPK、HIF-1、NF-κB等癌症相关通路，以及核受体活性、配体激活转录因子活性等分子功能。

**3.3 基于机器学习的核心靶点筛选**
在GSE16515训练集和GSE119224、GSE15471验证集中，通过175种算法组合筛选，确定14个核心基因（在≥10种算法中保留）。最优模型Lasso+plsRglm在三个数据集的AUC分别为0.983、1.000、0.936。

**3.4 基于核心靶点的诊断预测模型构建**
基于14个基因构建列线图，在GSE15471中AUC为0.984，校准曲线良好；在TCGA-PAAD+GTEx外部验证队列中AUC为0.997。单基因ROC分析显示AHR（AUC=0.915）、GPRC5A（AUC=0.930）、INHBA（AUC=0.963）诊断效能最高，且三者均在胰腺癌组织中显著高表达。

**3.5 免疫浸润、GSEA、单细胞转录组及蛋白表达分析**
CIBERSORT免疫浸润分析显示，胰腺癌组织中B细胞 naive、CD8 T细胞、NK细胞等减少，而B细胞记忆、浆细胞、M0巨噬细胞等增加。AHR阳性调节Treg细胞；GPRC5A与M1巨噬细胞负相关；INHBA与CD274、CD86等正相关。单细胞分析表明，AHR与应激反应相关（cor=0.326），GPRC5A与分化相关（cor=0.598），INHBA与侵袭相关（cor=0.464）。HPA免疫组化证实AHR和GPRC5A在胰腺癌组织中蛋白表达强阳性。

**3.6 分子对接与分子动力学模拟**
分子对接显示BPA与AHR结合能为?7.55 kcal/mol，优于GPRC5A（?4.91）和INHBA（?4.50）。三组独立200 ns分子动力学模拟结合MM/GBSA计算表明，AHR-BPA结合自由能（?57.55±7.20 kcal/mol）约为GPRC5A（?28.87±12.96）和INHBA（?30.81±14.51）的两倍，且与已知AHR配体（TCDD：?65.3±5.1；苯并[a]芘：?48.2±6.8 kcal/mol）相当。能量分解显示，范德华力和脂溶性相互作用是主要驱动力。CETSA实验证实BPA可增强AHR蛋白热稳定性。

**3.7 BPA通过AHR抑制正常胰管上皮细胞生长并促进胰腺癌细胞生长**
BPA浓度依赖性抑制hTERT-HPNE细胞活力及集落形成，并促进细胞死亡；但对AsPC-1和MiaPaCa-2细胞短期活力无显著影响，长期暴露则促进集落形成。BPA暴露上调AHR mRNA和蛋白表达。AHR敲低后，BPA对AsPC-1和MiaPaCa-2细胞的增殖促进作用显著减弱，提示AHR介导了该效应。

**讨论与结论**
讨论部分指出，本研究通过人群流行病学、计算毒理学、生物信息学、机器学习、分子动力学模拟和实验验证的多维框架，系统揭示了BPA与胰腺癌的关联。NHANES分析提供了人群水平的环境健康相关性证据，14个核心基因的诊断模型显示出卓越的区分能力，其中AHR、GPRC5A、INHBA为最具诊断潜力的靶点。功能富集表明BPA可能通过PI3K-Akt、MAPK、HIF-1等通路影响胰腺癌，免疫浸润分析提示BPA可能通过重塑肿瘤免疫微环境发挥作用。分子动力学模拟和CETSA实验证实AHR是BPA的直接靶点，且BPA通过AHR发挥细胞类型依赖的双相效应。结论部分翻译如下：总之，这些多层面发现提示BPA是胰腺癌的潜在环境贡献因素，并强调AHR作为候选分子中介，但需进一步的前瞻性和机制验证。