《Advanced Science》：Genomics Insights Into High-Latitude Adaptation of Tibetan Macaques
編輯推薦：
本文系統解析了藏獼猴（Macaca thibetana）適應高緯度寒冷環境的遺傳機制。研究通過構建高質量參考基因組（BUSCO完整性達99.41%），結合基因編輯（CRISPR-Cas9）和定量CT（QCT）等技術，首次揭示TBX6基因Pro71Thr突變通過調控尾椎發育（RIPPLY1/FGF8表達上調）導致短尾表型；發現脂代謝相關基因（DGAT2/DYSF/CAV1/PRKD1/LEPR/CPE）的正向選擇和基因家族擴張（如氧化磷酸化通路），驅動腹腔脂肪積累提升9.3倍。群體基因組學進一步揭示東西部種群在脂代謝基因（如ELOVL6）上的差異化選擇，為靈長類環境適應提供了創新性見解。
藏獼猴高緯度適應性的基因組學機制
1 引言
非人靈長類在高緯度地區的分布極為有限，而藏獼猴（Macaca thibetana）作為完全分布于北回歸線以北的稀有物種，為研究寒冷適應機制提供了獨特模型。其典型短尾（相對尾長僅12%）、厚密體毛及大型體型（雄性平均體重>15 kg）符合伯格曼法則（Bergmann's rule）和艾倫法則（Allen's rule），暗示其形態與生理特征可能為高緯度寒冷環境適應的結果。然而，此前基因組注釋不完整（BUSCO評分僅81.8%）限制了對其遺傳基礎的深入解析。
2 結果
2.1 基因組優化與系統發育
本研究通過NCBI真核基因組注釋流程注釋29,108個基因，平均基因長度提升至69,527 bp，BUSCO完整性達99.41%。系統發育分析顯示藏獼猴與阿薩姆猴（M. assamensis）親緣最近，分化時間約2.06百萬年。
2.2 基因組變異特征
藏獼猴基因組中檢測到448個顯著擴張的基因家族（涉及1,975個基因），富集于氧化磷酸化（mcc00190）、糖酵解/糖異生（mcc00010）等能量代謝通路。同時鑒定出345個譜系特異性正選擇基因（PSGs），功能集中于脂質儲存（GO:0019915）等過程。結構變異（SV）分析發現1,175個藏獼猴特異性SV（TMSSVs），其中661個與蛋白編碼基因重疊。
2.3 TBX6突變導致短尾表型
在尾發育關鍵基因TBX6中發現藏獼猴特異性純合錯義突變c.211C>A（p.Pro71Thr），CADD評分18.55提示有害效應。通過CRISPR-Cas9構建TBX6Pro71Thr點突變小鼠模型，發現突變體尾椎數量顯著減少（24–25節 vs. 野生型28–29節），且下游基因RIPPLY1（E10.5時期表達上調6倍）和FGF8（E11.5時期上調4倍）表達異常，證實該突變通過干擾體節發生通路調控尾長。
2.4 短尾表型的多基因調控基礎
除TBX6外，其他短尾獼猴物種（如日本獼猴M. fuscata和蘇拉威西冠毛獼猴M. nigra）亦存在TBX6非同義突變，提示該基因在多物種尾長演化中的保守作用。同時發現FGFR1基因p.Val102Ile突變（CADD評分24.7）及WNT3A基因第二內含子區300 bp插入可能通過影響FGF和Wnt信號通路協同調控尾發育。
2.5 脂代謝適應機制
QCT分析顯示藏獼猴腹腔脂肪體積（1,190.68 cm3）為恒河猴（127.99 cm3）的9.3倍。基因組學證據包括：（1）脂儲存基因DGAT2（催化三酰甘油合成）、DYSF（脂滴形成）和CAV1（膽固醇轉運）受正選擇；（2）脂肪生成基因PRKD1兩個錯義突變（p.Ala381Thr/p.Val558Met）促進脂質積累；（3）食欲調控基因LEPR（p.Pro1011Ser）表達下調0.19倍，而CPE基因390 bp缺失導致其表達上調1.73倍。qRT-PCR驗證上述基因在脂肪組織中表達顯著改變。此外，乳酸脫氫酶（LDH）和NADH脫氫酶等基因家族擴張，增強能量代謝靈活性。
2.6 群體遺傳學分析
26只野生藏獼猴的全基因組重測序揭示東西部種群（以梵凈山FJS為界）顯著分化。東部種群（如黃山HS群體）遺傳多樣性最低（π=3.49×10?4），近交系數最高；西部種群（如FJS）多樣性最高（π=9.05×10?4）。選擇性清除分析顯示東部種群在脂肪酸延伸基因ELOVL6等脂合成通路基因上受到強烈選擇，而西部種群偏好酮體代謝基因（如BDH1/BDH2）。
3 討論
3.1 短尾的遺傳基礎
藏獼猴短尾表型由TBX6、FGFR1和WNT3A等多基因變異協同調控，其中TBX6Pro71Thr突變通過擾亂體節發生邊界建立導致尾椎減少。這一適應特征可能由高緯度寒冷環境驅動，以減少熱量散失并促進陸地活動轉化。
3.2 脂代謝的高緯度適應
藏獼猴通過基因突變（如CPE缺失）、正選擇（如DGAT2）及基因家族擴張（如氧化磷酸化通路）等多層次遺傳變異，實現脂質高效儲存與動員，與北極熊、因紐特人等寒冷適應物種的代謝策略趨同。與北非獼猴（M. sylvanus）的平行進化分析進一步驗證脂代謝通路在靈長類高緯度適應中的保守性。
3.3 種群演化與保護啟示
藏獼猴種群約24萬年前分化為東西支系，東部種群近期經歷瓶頸效應，遺傳多樣性嚴重衰減。全球變暖背景下，已特化出短尾、厚脂特征的藏獼猴（尤其是東部種群）面臨嚴峻適應挑戰，需優先制定保護策略。
4 材料與方法
涵蓋基因組注釋、基因家族分析、正選擇檢測、基因編輯模型構建、群體遺傳學分析等標準化流程，確保研究結論的可重復性。
