整合免疫细胞eQTL与孟德尔随机化分析鉴定乳腺癌易感性相关遗传支持基因

《Cancer Genetics》：Integrative analysis of immune-cell eQTLs and Mendelian randomization identifies genetically supported genes associated with breast cancer susceptibility

　　

研究人员通过整合免疫细胞表达数量性状位点（expression quantitative trait loci, eQTL）与全基因组关联研究（genome-wide association study, GWAS）数据，优先识别出与乳腺癌风险相关的基因。免疫细胞中APOBEC3A/APOBEC3B的表达与较低的乳腺癌风险相关。多资源孟德尔随机化（Mendelian randomization, MR）、共定位（colocalization）分析及基于汇总数据的孟德尔随机化（summary-data-based Mendelian randomization, SMR）分析为候选因果基因提供了支持。优先识别的基因涉及抗原呈递和雌激素信号通路。药物-基因分析则凸显了优先基因中潜在的治疗靶点。

乳腺癌是全球女性最常见的恶性肿瘤和癌症相关死亡的主要原因，尽管诊疗不断进步，其分子和临床异质性仍使遗传易感机制解析面临巨大挑战。大规模全基因组关联研究（genome-wide association study, GWAS）已发现数百个乳腺癌易感位点，但多数位于非编码区，靶基因不明。免疫调控在乳腺癌发生发展中具有重要作用，而既往孟德尔随机化（Mendelian randomization, MR）研究多基于循环免疫表型或组织水平分子特征，难以捕捉细胞类型特异的调控效应。该研究发表在《Cancer Genetics》上，旨在整合免疫细胞表达数量性状位点（expression quantitative trait loci, eQTL）资源与乳腺癌GWAS数据，通过两样本MR框架系统识别具有遗传支持的乳腺癌易感相关基因。

研究人员采用两样本MR设计，以TenK10K单细胞eQTL资源为主要暴露数据集，并纳入OneK1K、DICE、eQTLGen和Soskic四个独立eQTL资源进行验证；结局数据来自FinnGen R12乳腺癌GWAS（24,270例病例，222,078例对照，欧洲血统）。方法上采用逆方差加权（inverse-variance weighted, IVW）和Wald比率方法估计因果效应，结合贝叶斯共定位（coloc）、SuSiE精细定位、LD检查及基于汇总数据的孟德尔随机化（summary-data-based Mendelian randomization, SMR）进行遗传验证；利用GTEx v10组织eQTL进行跨组织验证，并在BCAC队列（76,192例病例，63,082例对照）中独立重复；同时使用单细胞RNA测序和空间转录组数据验证基因表达，并进行功能富集、蛋白质相互作用网络和药物-基因相互作用分析。

研究结果分述如下。免疫细胞eQTL资源与乳腺癌GWAS数据概览：TenK10K包含28种免疫细胞类型、12,860个基因的97,833个独立工具变量；其余四个资源亦提供不同规模的eQTL记录。基于单细胞eQTL的MR分析及遗传验证揭示免疫细胞特异性的乳腺癌关联：MR分析共获得516个FDR<0.05的显著关联，涉及131个基因和28种细胞类型；最强信号位于染色体22q的APOBEC3A/APOBEC3B基因座，所有6个相关关联均呈保护效应，提示免疫细胞中APOBEC3表达升高与乳腺癌风险降低相关；结合共定位、LD检查和SuSiE精细定位，433个关联（106个基因）具有收敛遗传证据。SMR分析验证MR优先的免疫细胞特异性关联：在五种eQTL资源中，42个基因至少在1个数据集中达到SMR显著且通过HEIDI检验，17个在≥2个数据集中验证，YBEY、KCNN4和ATG10在4个数据集中获得验证；TenK10K中MR与SMR效应方向一致率为95.4%。基于GTEx的SMR分析揭示MR优先基因的组织特异性调控效应：乳腺组织中10个基因达到FDR显著，其中APOBEC3B信号最强（p_SMR=1.09×10^-8，OR=0.85），且所有FDR显著基因的效应方向均与MR一致。BCAC独立复制和ER亚型分析：506个可测试关联中346个（68.4%）达到名义显著，方向一致率88.3%；FGFR2和SLC22A5的效应局限于雌激素受体阳性（ER+）乳腺癌，而ESR1在雌激素受体阴性（ER-）乳腺癌中效应更强。乳腺癌相对于其他妇科癌症的特异性：在子宫内膜癌、卵巢癌和宫颈癌GWAS中，106个具有收敛证据的基因中87%无FDR显著关联，支持这些关联的乳腺癌特异性。乳腺肿瘤中的单细胞和空间转录组验证：在26例乳腺癌患者的100,064个单细胞和6个空间转录组切片中，12个优先基因呈区室特异性表达；HLA-DRB1和HLA-DQA1在免疫丰富区域高表达，并与髓系/树突状细胞丰度呈正相关；APOBEC3A主要表达于髓系细胞。功能富集突出免疫和抗原呈递通路：基因本体（GO）分析显示“单核细胞分化正调控”等4个术语达到FDR显著；KEGG通路未达FDR显著，但名义显著富集于自身免疫及抗原处理呈递通路。细胞类型特异性富集仅显示调节性T细胞的标称富集：28种免疫细胞类型中，仅调节性T细胞（Treg）达到名义富集显著（12.45%，P=0.0497），多重校正后无细胞类型显著。蛋白相互作用分析识别有限但可解释的连接：STRING高置信度网络中106个基因仅形成6条边，最大连通组件为XBP1–ESR1–FGFR2–CSF1，ESR1处于中心位置。药物-基因相互作用分析识别潜在治疗靶点：24个基因在DGIdb中有药物相互作用记录，ESR1（229条已批准药物记录）和FGFR2（67条）拥有最多已批准药物，包括内分泌治疗药物和FGFR抑制剂，提示其潜在治疗价值。

讨论部分指出，该研究通过整合多种免疫eQTL资源与乳腺癌GWAS，结合MR、共定位、SMR和GTEx验证，提供了免疫细胞调控及APOBEC介导的DNA编辑与乳腺癌风险相关联的收敛遗传证据。APOBEC3A/APOBEC3B位点的保护效应在独立队列中复制，且与肿瘤细胞中APOBEC3的致突变作用形成对比，提示细胞背景依赖的生物学功能。功能网络和药物注释聚焦于ESR1等激素相关基因，支持其治疗靶点潜力。研究局限性包括欧洲血统样本可能限制普适性、多数暴露仅由单一SNP工具变量、共定位区域定义可能降低功效等。研究结论可翻译为：本研究整合多种免疫eQTL资源与乳腺癌GWAS数据，通过两样本MR框架识别了具有遗传支持的免疫细胞特异性乳腺癌易感基因；MR、共定位、SMR及GTEx验证共同提供了收敛证据，将免疫细胞调控和APOBEC介导的DNA编辑与乳腺癌风险联系起来，网络分析进一步聚焦于激素相关基因。这些发现加深了对乳腺癌免疫介导遗传易感性的理解，并为未来功能研究和治疗开发提供了候选靶点。