《Biochemistry and Biophysics Reports》：Mammalian evolution and human mutation burden in Rab GTPases
編輯推薦：
本文聚焦于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超家族或其他細胞內運輸調節因子），從而更系統地揭示進化歷史如何塑造當代人類的疾病易感性，為發現新的疾病治療靶點和理解疾病的進化基礎開辟了新路徑。
