编辑推荐：
本文聚焦于Rab GTPases这一关键膜转运调节蛋白家族，研究人员通过整合哺乳动物进化分析与人类突变数据，旨在揭示其在家系水平上的进化与突变模式。研究发现，在更多哺乳动物谱系中发生正选择的Rab蛋白，其人类同源物的Switch I结构域倾向于积累更多有害突变，表明该区域在物种间具有变异灵活性，同时可能容纳与人类疾病相关的变异。这项研究为理解保守蛋白家族的进化如何影响人类疾病易感性提供了新视角。
在细胞内部，物质如何被精准地运输到指定位置，犹如一个高度组织化的物流系统，而Rab GTPases（鸟苷三磷酸酶）正是这个系统中的核心“交通指挥官”。作为Ras超家族的重要成员，Rab蛋白通过在其活性（GTP结合）与非活性（GDP结合）状态间切换，精确调控着细胞内吞、囊泡运输、自噬等关键膜转运过程。然而，一旦这些“指挥官”功能失调，例如因编码突变或表达异常，就可能引发一系列严重的人类疾病，包括多种癌症和帕金森病等神经退行性疾病。尽管针对单个Rab蛋白在疾病和进化背景下的研究已非常深入，但一个根本性问题仍未得到充分解答：作为一整个蛋白家族，Rab GTPases在漫长的哺乳动物进化历程中，其序列演变遵循着怎样的整体模式？这些跨物种的进化动态，又如何影响了其在现代人类群体中与疾病相关的突变积累（即“突变负担”）？解答这个问题，对于从更深层次理解疾病易感性的进化根源、识别关键的致病性变异区域具有重要意义。
为此，来自都柏林大学学院的研究团队Unmani Sidor、Graham M. Hughes和Jeremy C. Simpson开展了一项开创性的整合研究。他们巧妙地将比较基因组学与人类群体遗传学数据相结合，构建了一个能够同时审视物种间进化与物种内变异的研究框架。研究人员首先对62种胎盘哺乳动物中的54个Rab蛋白进行了大规模的进化分析，利用分支位点模型（aBSREL, MEME）检测了它们在各个哺乳动物谱系中经历“正选择”（即有利于物种适应的突变被保留）的程度。随后，他们从UniProt数据库中获取了对应人类Rab蛋白的突变数据，并创新性地定义了三个结构域水平的量化指标——约束分数、损伤耐受性和总体突变负担——来系统评估人类Rab蛋白及其各个功能结构域（包括GTP结合区、Switch I/II结构域和高变区）对群体水平变异的耐受情况。
研究团队运用了几项关键的技术方法：1）利用TOGA工具从基因组比对中推断直系同源基因，并结合Clustal Omega进行多序列比对，为54个Rab基因在62种哺乳动物中构建了序列数据集；2）使用HyPhy软件套件中的aBSREL和MEME方法进行正选择检测，以识别在特定进化分支或位点上经历适应性变化的Rab基因；3）基于UniProt注释和PolyPhen-2预测分数，对人类49个Rab蛋白的错义变异进行致病性分类和注释，并计算上述三个结构域水平的突变负担指标；4）通过Jukes-Cantor距离计算，量化了各哺乳动物物种Switch I结构域序列与人类参考序列之间的差异程度，以评估该区域的进化分歧。
2.1. Rab在哺乳动物中的进化程度存在差异
对54个Rab蛋白的分析显示，其中53个在至少一个哺乳动物谱系中显示出显著的正选择信号。这表明Rab家族在哺乳动物进化过程中广泛经历了适应性序列变化。研究发现，进化最为广泛的Rab成员包括RAB36、RAB33B、RAB34和RAB17，它们在大量哺乳动物物种中都表现出正选择迹象。然而，Rab的进化程度与物种水平的特征（如癌症死亡率风险、寿命和体重）之间并未发现显著相关性，提示Rab的适应性变化可能并非直接驱动这些宏观物种差异的主要因素。
2.2. 人类Rab的突变负担反映了哺乳动物进化动态
研究人员将哺乳动物的进化数据与人类Rab的突变数据联系起来分析。他们发现，那些在更多哺乳动物物种中经历正选择的Rab蛋白，其对应的人类同源物往往在某些关键结构域表现出更高的突变负担。这一关联在三个量化指标（约束分数、损伤耐受性、总体突变负担）中均得到体现。值得注意的是，结构域特异性分析显示，突变负担的升高主要集中于特定的功能区域，而非均匀分布在整个蛋白序列上。
2.3. Switch I区域是突变负担的焦点
进一步的相关性分析揭示，Rab的进化程度与人类突变负担之间的关联，在Switch I结构域和GTP结合基序中最为强烈。尤其是Switch I结构域，在所有指标和标准化方法下都显示出与Rab在哺乳动物中进化范围最显著的正相关。这表明Switch I是一个进化上灵活、同时又能在人类群体中容纳潜在有害变异的“热点”区域。Switch I结构域在Rab蛋白中负责GTP结合并介导与效应蛋白的相互作用，其较高的变异容忍度可能有助于在维持核心GTPase功能的同时，实现物种特异性的调控微调。
2.4. Switch I区域的差异揭示了谱系特异性模式
除了突变负担，研究还评估了Switch I结构域在物种间的序列分歧程度。通过计算各哺乳动物与人类之间Switch I区域的Jukes-Cantor距离，研究人员发现某些物种（如大狐蝠、浣熊、红腿白臀叶猴）在多个Rab蛋白的Switch I区域都表现出与人类较高的序列差异。这种模式暗示了Switch I的变异具有谱系特异性，而非孤立事件。然而，Switch I的分歧度与物种的癌症死亡率风险或体重同样没有普遍关联，仅与少数特定Rab（如RAB36、RAB9A、RAB34）的物种寿命存在或正或负的相关性。
归纳研究结论与讨论
这项研究的主要结论在于，它首次在家系水平上系统描绘了Rab GTPases在哺乳动物中的进化全景图，并成功将其与人类疾病相关的遗传变异联系起来。研究发现，RAB36、RAB33B、RAB34和RAB17是哺乳动物中进化最为广泛的Rab成员，它们频繁经历正选择，暗示这些蛋白可能在适应不同的膜转运需求中扮演特殊角色。更为重要的是，研究揭示了一个关键机制：那些在跨物种进化中表现出高度灵活性的Rab蛋白，其人类版本在关键的Switch I调控结构域倾向于积累更多有害突变。这意味着Switch I结构域在进化上扮演了一个“双重角色”：一方面，其序列的灵活性允许其在不同哺乳动物谱系中发生适应性改变；另一方面，这种对变异的容忍度也可能使得该区域在现代人类群体中成为疾病相关突变的潜在“蓄水池”。
在讨论中，作者指出，尽管Rab蛋白的失调与癌症进展密切相关，但本研究并未发现Rab的进化与物种水平的癌症死亡率风险或寿命等特征有直接关联。这表明Rab的序列变化可能更多地反映了对膜转运通路本身的微调适应，而非直接贡献于宏观的癌症防御机制或寿命延长。同时，作者也谨慎地解释了进化与突变负担的关联：这种关联可能源于某些蛋白区域对序列变化的内在容忍性，这种容忍性既允许其在物种间适应，也允许其在人群内积累变异，但并不一定意味着直接的因果关系。
该研究的整合分析框架——将跨物种的正选择分析与人类群体遗传学的结构域水平突变负担评估相结合——具有很强的普适性和启发性。它不仅深化了我们对Rab这一重要蛋白家族进化与功能关系的理解，更重要的是，提供了一种可推广的研究范式。未来，这一框架可被应用于研究其他与疾病相关的保守蛋白家族（如整个Ras超家族或其他细胞内运输调节因子），从而更系统地揭示进化历史如何塑造当代人类的疾病易感性，为发现新的疾病治疗靶点和理解疾病的进化基础开辟了新路径。
