眼黑色素瘤的进化起源：与rs12913832 G等位基因频率和纬度的关联

《Cancer Medicine》：Evolutionary Origin of Ocular Melanoma: Associations With rs12913832 G Allele Frequency and Latitude

　　

背景：眼黑色素瘤（OM）的发病率在世界范围内存在差异，但地理纬度与色素沉着遗传学的相对影响仍存在争议。研究人员调查了国家层面rs12913832衍生（G）蓝眼等位基因频率（G等位基因频率）或纬度是否能更好地解释OM发病率的国际差异。方法：研究人员利用国际癌症研究机构（IARC）的年龄标准化OM发病率数据、公开数据库中的G等位基因频率以及基于夜间灯光数据得出的群体加权重心，对29个国家进行了一项观察性生态学研究。研究人员进行了相关分析、线性回归和中介分析，以评估纬度、G等位基因频率和OM发病率之间的关联。结果：在单变量分析中，纬度和G等位基因频率均与OM发病率正相关。在多变量回归中，纬度不再显著，而G等位基因频率解释了显著更多的方差，并且其解释的方差（偏r²）比纬度多61倍。中介分析表明，纬度对发病率的影响中约92%是通过G等位基因频率介导的。结论：在群体水平上，G等位基因频率是比地理纬度更强的OM发病率预测因子。作为一项生态学分析，本研究确定了群体层面的关联，而非个体层面的遗传风险，并提示色素沉着遗传学是OM发病率国际差异的重要原因。

眼黑色素瘤（ocular melanoma, OM）的发病率在全球范围内呈现显著差异，从印度的每年每百万人0.1例到丹麦的11.3例不等。长期以来，地理纬度被认为是潜在影响因素，因为发病率随赤道距离增加而升高，尤其在欧洲裔人群中更为明显。然而，这一纬度梯度是否直接由紫外线（UV）暴露驱动，抑或是由于人口中色素沉着相关遗传变异（如决定蓝眼颜色的rs12913832 G等位基因）的纬度分布所介导，一直是争论的焦点。已有候选基因研究和全基因组关联研究（GWAS）将rs12913832 G等位基因与葡萄膜黑色素瘤（uveal melanoma, UM）风险直接关联，但缺乏系统性的跨国比较。为了阐明这两种因素的解释力，研究人员利用29个国家的年龄标准化OM发病率数据、rs12913832 G等位基因频率（G allele frequency）以及基于夜间灯光数据计算的人口加权重心（population-weighted centroid），进行了生态学分析。该研究发表于《Cancer Medicine》。

研究人员开展了一项观察性生态学研究，主要技术方法包括：从国际癌症研究机构（IARC）获取29个国家的年龄标准化OM发病率数据（2008–2012年为主，部分国家时间范围有差异）；从耶鲁大学维护的等位基因频率数据库（ALFRED）提取各人群的rs12913832 G等位基因频率，并根据各国人口构成计算加权平均值；利用Hall等（2019）发布的基于夜间灯光数据的人口加权重心坐标（2013年估算），作为各国人口分布的中心点，以绝对纬度（度）表示地理纬度。统计方法包括Pearson相关分析、简单线性回归、多元线性回归、偏相关分析，以及中介分析（bootstrap 1000次模拟），评估纬度效应中由G等位基因频率介导的比例。此外，进行了留一法敏感性分析和随机改变G等位基因频率±25个百分点的敏感性检验。

研究结果如下：

**描述性统计（Descriptive Statistics）**：29个国家的年发病率从印度的0.1例/百万/年到丹麦的11.3例/百万/年；G等位基因频率从撒哈拉以南非洲、日本、韩国的小于1%到芬兰的88%不等。人口加权重心纬度范围从肯尼亚的-1.2°到挪威的63.4°。

**相关与回归分析（Correlation and Regression Analyses）**：单变量分析显示，纬度与OM发病率呈正相关（r=0.76, p<0.001），G等位基因频率也与发病率正相关（r=0.92, p<0.001）。简单线性回归中，纬度每增加1°与发病率增加0.16例/百万/年相关；G等位基因频率每增加1个百分点与发病率增加0.10例/百万/年相关。在纳入两者的多元回归模型中（n=29），G等位基因频率仍显著（β=0.09, p<0.001），而纬度不再显著（β=0.01, p=0.46）。标准化回归系数显示G等位基因频率的效果（β_std=0.90）远大于纬度（β_std=0.06）。

**偏相关分析（Partial Correlations）**：G等位基因频率的偏相关系数r_partial≈0.78，解释了纬度未解释方差中约61%（偏r²=0.61）；纬度的偏相关系数r_partial≈0.09，解释方差小于1%（偏r²<0.01）。即G等位基因频率在联合模型中解释的方差是纬度的61倍。

**中介分析（Mediation Analysis）**：纬度对发病率的总效应为0.16，其中平均因果中介效应（ACME）为0.15（95%CI: 0.12–0.22, p<0.001），平均直接效应（ADE）为0.01（95%CI: –0.01–0.07, p=0.31）。约92%的纬度效应通过G等位基因频率中介。

**敏感性分析（Sensitivity Analyses）**：留一法排除任一国家后，G等位基因频率的回归系数均保持在0.089–0.095之间（均p<0.001），纬度系数始终不显著（p>0.54）。随机扰动G等位基因频率±25个百分点（10,000次迭代）后，9980次中G等位基因频率显著（p<0.05），而纬度仅5777次显著。结果稳健。

在讨论部分，研究人员指出，尽管UV暴露是皮肤黑色素瘤的明确病因，但对OM的直接作用有限，因为大部分UV被角膜和晶状体吸收。然而，rs12913832 G等位基因通过下调OCA2表达减少真黑素（eumelanin）合成，使葡萄膜黑素细胞更易受氧化应激损伤，且纯合GG基因型患者常伴有更具侵袭性的细胞遗传学特征。进化分析表明，该突变约在6000–10000年前起源于黑海地区，而该区域更早的选择信号（超过42000年）可能涉及更广的色素沉着相关变异。研究局限性包括：生态学设计无法推断个体水平因果关系；国家层面数据可能掩盖域内变异；G等位基因频率基于特定人群加权平均，存在估算误差；OM作为整体分析无法分离UV敏感的结膜黑色素瘤和虹膜黑色素瘤。结论部分翻译如下：

**结论**：本研究的分析表明，在群体水平上，rs12913832的衍生（G）等位基因是眼黑色素瘤发病率的显著强于地理纬度的预测因子。纬度和发病率之间的几乎所有关联都是由G等位基因频率的差异介导的，考虑遗传变异后纬度贡献极小。由于该研究是生态学设计，这些结果识别了国家层面的关联，并未建立个体水平的遗传风险；然而，后者已由单独的候选基因研究和全基因组关联研究支持。综上所述，这些发现表明色素沉着相关遗传因素是与眼黑色素瘤发病率国际差异相关的主要因素。需进一步研究以阐明该关系背后的生物学机制。