古基因組解密末次盛冰期前后的東亞北部人群動態

毛曉偉，平婉菁，付巧妹

古基因組解密末次盛冰期前后的東亞北部人群動態

毛曉偉，平婉菁，付巧妹

中國科學院古脊椎動物與古人類研究所，北京 100044

現代人的起源與演化一直是地球與生命科學領域的研究焦點。過去10多年，圍繞世界各地考古遺址人群開展的古基因組研究，在舊石器時代晚期以來歐洲、北亞等區域人群的演化問題上取得重大進展。然而，東亞人群古基因組研究則相對滯后，尤其是對中國人類古基因組的研究極度匱乏，使東亞人群遺傳演化和遷徙融合的歷史成為現代人演化進程里的重要“缺環”。但是，近幾年中國科學家在這方面的研究取得了突破性進展。

2017年，中國科學院古脊椎動物與古人類研究所(簡稱“中科院古脊椎所”)付巧妹團隊從我國北京4萬年前的田園洞人化石中獲取中國以至東亞最古老現代人的基因組序列，數據分析表明以其為代表的古東亞人群與同時期某一古歐洲人群存在特殊的遺傳聯系，首次為探索整個歐亞大陸人群的遺傳歷史開啟東亞篇章[1]。2020年，其團隊首次系統性研究距今5000年青藏高原人群的線粒體基因組揭示青藏高原人群由低海拔向高海拔地區少量的人群流動和之后一段時間高海拔地區內部的人群擴張[2]。同年，其團隊又相繼發表我國南方1.1萬年前人類古線粒體基因組及東亞自9500年來南北方人群的大規模、系統性古基因組研究，表明萬年來東亞南北方人群在內頻繁遷徙互動、南島語族人群和南方內陸古人群的聯系、向外通過沿海地帶不斷擴散交流的復雜歷史[3,4]。2021年，其團隊再度針對我國北方沿海的山東人群及河南仰韶遺址進行規模性古線粒體基因組研究，揭示其人群自新石器時代以來的遺傳連續性[5,6]。這些研究不僅從時間和空間上填補了東亞人類古基因組的巨大空缺，揭開東亞人群與歐亞其他人群此前未知的遺傳聯系以及遷徙互動在其演化過程中的重要作用，更為人們解析了東亞北部自4萬年前田園洞人至近萬年前古北方沿海人群之間長達數萬年的遺傳歷史斷層。

冰河時期(距今約7～1.15萬年)的氣候變化對歐洲人群的遷徙和種群大小變化構成影響[7]，考古證據表明其對亞洲高緯度和高海拔地區的人群歷史也有類似的影響[8]，但缺少相關的遺傳學證據與針對性研究。在這一時期，尤其是末次盛冰期前后的東亞北部人群是否因氣候的劇烈變化而發生相應人群的更替或遷徙交流？東亞人群是否攜有與環境適應性相關的特有基因？這些都是古基因組研究亟待探索而頗具挑戰的關鍵問題。

2021年5月27日，雜志在線發表了中科院古脊椎所付巧妹團隊與云南大學張虎才研究組(第一作者為中科院古脊椎所毛曉偉副研究員，通訊作者為中科院古脊椎所付巧妹研究員和云南大學張虎才教授)合作開展的東亞北部人群的古基因組研究[9]。在該研究中，團隊通過最前沿的古基因組實驗技術獲取我國東北地區黑龍江省距今33,000～3400年的25個早期人類樣本的古基因組數據(其中11例個體超過1萬年，圖1)，填補東亞北部人群遺傳歷史的相關時間斷層，并結合北京田園洞人基因組數據展開遺傳演化與群體遺傳學分析，系統繪制出4萬年來東亞北部人群的群體動態演化圖譜。此外，研究發現東亞人群特有的環境適應性基因至少在1.9萬年前出現，對此前利用當今現代人基因組所估算出的突變時間予以精確修正，為探究東亞北部人群的適應性遺傳變異情況及相關基因的選擇機制提供了新的思路和證據。

圖1 新樣本和已發表的東亞北方個體信息地理和時間分布圖[9]

該研究指出生活在末次盛冰期前、距今約3.3萬年的AR33K個體與距今4萬年田園洞人的遺傳成分相似。此前有研究表明同時期來自蒙古東北部、距今3.4萬年的索爾吉特人(Salkhit)同時攜有中國田園洞人和古西伯利亞北方人群“亞納人”(Yana)相關的祖源成分[10]，而本研究顯示AR33K個體與同時期索爾吉特人(Salkhit)不一樣，他不含有“亞納人”(Yana)相關的遺傳成分，與田園洞人是聚類在一起。這一結果說明在末次盛冰期之前，田園洞人相關人群在東亞北部廣泛存在；同時，田園洞人相關人群在蒙古地區與亞那人(Yana)相關人群混合，而此相關人群在黑龍江地區存在并保持著連續性(圖2)。

末次盛冰期約在2.65萬年前開始，一直持續到約1.9萬年前，這是距現在最近的極寒冷時期。為進一步探究末次盛冰期氣候變化如何影響東亞北部人群，該研究針對生存在末次盛冰期末期、距今1.9萬年的AR19K個體的古基因組展開分析。這是東亞首例末次盛冰期前后的人類基因組，分析表明在末次盛冰期前后(1.9萬年)北部發生了人群特殊變化，即之前廣泛存在的田園洞人相關人群可能消失，東亞古北方人群此時已在黑龍江出現。一方面，其數據顯示AR19K個體相較于4萬年前田園洞人，他并沒有與同一地區的3.3萬年前的個體更近，間接反映了人群的不連續。另一方面，研究顯示AR19K個體相比于東亞古南方沿海人群(距今8400～4200年的奇和洞、亮島、溪頭村、曇石山遺址個體)來說，與東亞古北方沿海人群(距今9500～7700年的扁扁洞、小荊山、小高遺址、博山遺址個體)的遺傳關系更近。由此表明，在末次盛冰期即將結束的時候，東亞南北方人群的差異就已經形成，并且最早的東亞古北方人群(以AR19K個體為代表的人群)在黑龍江地區出現。這一結果將此前研究所發現的東亞南北方人群遺傳分化的時間(9500年前[3])向前推進了近萬年。

末次盛冰期之后，東亞北部人群的流動情況又是如何？此前有關新石器時代早期黑龍江流域鬼門洞的狩獵采集人群(來自俄羅斯遠東地區、距今7700年的DevilsCave_N個體)的古基因組研究表明，黑龍江流域人群在距今約8000年以來存在遺傳連續性[11,12]，然而距今8000年之前的人群情況尚不可知。本研究通過獲取末次盛冰期之后、距今8000年之前黑龍江地區的16例個體基因組數據(距今約1.4～0.9萬年)，對該地區這段空白的人群遺傳歷史展開研究并獲得有價值的發現。研究顯示，末次盛冰期后的黑龍江地區人群相比于其他東亞人群而言，與新石器時代早期鬼門洞人群的遺傳成分最為接近。這一結果表明，黑龍江流域人群的遺傳連續性可能早在1.4萬年前就已經開始，這比此前假說所提出的時間向前推進了6000年。該結論與黑龍江地區陶器和玉器最早出現在1.5萬年前的考古記錄相符[13]。另外，本研究發現生活在末次盛冰期后、距今1.4萬年前的AR14K個體與古西伯利亞人群(Ancient Paleo- Siberians)[12]有特殊的遺傳聯系。古西伯利亞人群是在美洲之外與美洲原住民最相關的人群，而數據顯示AR14K個體比鬼門洞人群更符合古西伯利亞人群所混有東亞祖源成分的來源。

圖2 末次盛冰期前后的東亞北部人群動態示意圖[9]

此外，研究對東亞北部人群的適應性遺傳變異情況展開探索。此前有全基因組關聯分析表明基因的V370A突變與東亞典型體征如更粗的毛發、更多的汗腺和鏟形門齒有關[14,15]，因此理解該基因突變的自然選擇機制，對研究東亞人群的適應性變異至關重要。目前對該基因的選擇機制有兩種假說：一是在約2萬年前的低紫外線環境中，該選擇增加了母乳中的維生素D含量[16]；二是在約3萬年前的溫暖潮濕氣候中，該選擇可調控體溫調節性出汗[17]。此外，相關研究利用當今現代人基因組估算出基因V370A的突變時間在距今1.14萬年(置信區間在距今4300～43,700年)。然而，本研究通過對大規模的古東亞人群基因組研究，更直接、深入地揭示該基因的出現時間和可能的選擇機制。研究發現，

除田園洞人及末次盛冰期之前的AR33K個體以外，基因的V370A突變出現在包括AR19K在內的所有已知的古東亞人群之中。這一結果表明，基因的V370A突變率在末次盛冰期階段或之后不久的時間里可能升高，此外也證明了原假說二的觀點：3萬年前基因的V370A突變在溫暖潮濕環境中被選擇的可能性較小。該結論使得所估算突變基因出現時間的上限提早并更為精確(圖2)。

隨著古基因組技術與遺傳學分析手段的不斷進步及國內專業團隊與考古機構多方、跨學科合作的大力開展，針對東亞人群的古基因組研究已然迎來飛速發展的時刻。本研究首次在東亞地區開展跨度為4萬年的大規模人類古基因組研究，且首次利用古基因組在適應性方面探究東亞人群重要表性特征的演化來源，揭示東亞北部距今40,000～3400年的的人群動態遺傳歷史，并為進一步探索東亞人群與環境的關系提供了重要遺傳學證據。此外，該合作研究團隊在通過古基因組技術與遺傳學分析手段的深度融合揭秘東亞人群演化方面樹立新的研究范式，為推進跨學科學術生態的建立提供借鑒和范本。

前沿技術的探索運用，是本研究得以取得重要進展與發現的前提，而在未來的東亞人群演化研究中亦將繼續發揮重要作用。在古基因組獲取方面，付巧妹團隊通過延展應用或共同創新的古基因組關鍵技術，極大擴展了古基因組研究樣本的時間軸線與區域范圍：如共同開發古核DNA捕獲技術，實現從大量環境微生物DNA里“釣取”出極其微量的內源DNA，使古基因組研究進入全基因組時代；如參與開發的新一代古DNA片段提取技術，將非冰凍層人類古DNA破譯的時間推進到40萬年前；如延展應用沉積物古DNA富集與種屬鑒定方法，成功從青藏高原“土”中獲取東亞首例已滅絕古人類DNA等等，成為世界上可實施操作這些古基因組關鍵技術的僅有的幾個國際頂尖實驗室之一，為國內科研團隊獨立、自主開展中國人群遺傳演化研究提供了技術保障。在遺傳學分析手段方面，該研究所運用的D統計量、人群混合結構、古老型人遺傳成分比例等分析，都為挖掘、準確解讀與深度分析不斷增長的、更大規模、海量而復雜的東亞人類古基因組數據提供了方法基礎。東亞有著極其豐富的考古遺存與復雜的人群歷史，古基因組研究與考古學、病理學、微生物學、計算機科學等更多學科的結合，將為人們從更多維度、更精確、細節性地了解東亞人類神秘的過去——他們的飲食、環境、遷徙、健康、功能性基因的選擇機制與適應性結果等打開更多新的窗口。

[1] Yang MA, Gao X, Theunert C, Tong H, Aximu-Petri A, Nickel B, Slatkin M, Meyer M, P??bo S, Kelso J, Fu Q. 40,000-year-old individual from Asia provides insight into early population structure in Eurasia., 2017, 27(20): 3202–3208.

[2] Ding M, Wang T, Min-Shan Ko A, Chen H, Wang H, Dong G, Lu H, He W, Wangdue S, Yuan H, He Y,Cai L, Chen Z, Hou G, Zhang D, Zhang Z, Cao P, Dai Q, Feng X, Zhang M, Wang H, Yang MA, Fu Q. Ancient mitogenomes show plateau populations from last 5200 years partially contributed to present-day Tibetans., 2020, 287(1923):20192968.

[3] Yang MA, Fan XC, Sun B, Chen CY, Lang JF, Ko YC, Tsang CH, Chiu Hunglin, Wang TY, Bao QC, Wu XH, Hajdinjak M, Ko AMS, Ding MY, Cao P, Yang RW, Liu F, Nickel B, Dai QY, Feng XT, Zhang LZ, Sun CK, Ning C, Zeng W, Zhao YS, Zhang M, Gao X, Cui YQ, Reich D, Stoneking M, Fu QM. Ancient DNA indicates human population shifts and admixture in northern and southern China., 2020, 369(6501): 282–288.

[4] Bai F, Zhang X, Ji X, Cao P, Feng X, Yang R, Peng M, Pei S, Fu Q. Paleolithic genetic link between Southern China and Mainland Southeast Asia revealed by ancient mitochondrial genomes., 2020, 65(12): 1125–1128.

[5] Liu JC, Zeng W, Sun B, Mao XW, Zhao YS, Wang F, Li ZG, Luan FS, Guo JF, Zhu C, Wang ZM, Wei CM, Zhang M, Cao P, Liu F, Dai QY, Feng XT, Yang RW, Hou WH, Ping WJ, Wu XH, Bennett EA, Liu YC, Fu QM. Maternal genetic structure in ancient Shandong between 9500 and 1800 years ago., 2021, doi: 10.1016/j.scib.2021. 01.029.

[6] Miao B, Liu Y, Gu W, Wei Q, Wu Q, Wang W, Zhang M, Ding M, Wang T, Liu J, Liu F, Cao P, Dai Q, Yang R, Feng X, Ping W, Hou W, Yuan H, Fu Q. Maternal genetic structure of a Neolithic population of the Yangshao culture., doi:10.1016/j.jgg.2021.04.005.

[7] Fu Q, Posth C, Hajdinjak M, Petr M, Mallick S, Fernandes D, Furtw?ngler A, Haak W, Meyer M, Mittnik A, Nickel B, Peltzer A, Rohland N, Slon V, Talamo S, Lazaridis I, Lipson M, Mathieson I, Schiffels S, Skoglund P, Derevianko AP, Drozdov N, Slavinsky V, Tsybankov A, Cremonesi RG, Mallegni F, Gély B, Vacca E, Morales MR, Straus LG, Neugebauer-Maresch C, Teschler-Nicola M, Constantin S, Moldovan OT, Benazzi S, Peresani M, Coppola D, Lari M, Ricci S, Ronchitelli A, Valentin F, Thevenet C, Wehrberger K, Grigorescu D, Rougier H, Crevecoeur I, Flas D, Semal P, Mannino MA, Cupillard C, Bocherens H, Conard NJ, Harvati K, Moiseyev V, Drucker DG, Svoboda J, Richards MP, Caramelli D, Pinhasi R, Kelso J, Patterson N, Krause J, P??bo S, Reich D. The genetic history of Ice Age Europe., 2016, 534(7606): 200–205.

[8] Rybin EP, Khatsenovich AM, Gunchinsuren B, Olsen JW, Zwyns N. The impact of the LGM on the development of the Upper Paleolithic in Mongolia., 2016, 425: 69–87.

[9] Mao X, Zhang H, Qiao S, Liu Y, Chang F, Xie P, Zhang M, Wang T, Li M, Cao P, Yang R, Liu F, Dai Q, Feng X, Ping W, Lei C, Olsen JW, Bennett EA, Fu Q. The deep population history of northern East Asia from the Late Pleistocene to the Holocene., 2021, doi: 10.1016/ j.cell.2021.04.040.

[10] Massilani D, Skov L, Hajdinjak M, Gunchinsuren B, Tseveendorj D, Yi S, Lee J, Nagel S, Nickel B, Devièse T, Tom H, Meyer M, Kelso J, Peter BM, P??bo S. Denisovan ancestry and population history of early East Asians., 2020, 370(6516): 579–583.

[11] Siska V, Jones ER, Jeon S, Bhak Y, Kim HM, Cho YS, Kim H, Lee K, Veselovskaya E, Balueva T, Gallego-Llorente M, Hofreiter M, Bradley DG, Eriksson A, Pinhasi R, Bhak J, Manica A. Genome-wide data from two early Neolithic East Asian individuals dating to 7700 years ago., 2017, 3(2): e1601877.

[12] Sikora M, Pitulko VV, Sousa VC, Allentoft ME, Vinner L, Rasmussen S, Margaryan A, de Barros Damgaard P, de la Fuente C, Renaud G, Yang MA, Fu Q, Dupanloup I, Giampoudakis K, Nogués-Bravo D, Rahbek C, Kroonen G, Peyrot M, McColl H, Vasilyev SV, Veselovskaya E, Gerasimova M, Pavlova EY, Chasnyk VG, Nikolskiy PA, Gromov AV, Khartanovich VI, Moiseyev V, Grebenyuk PS, Fedorchenko AY, Lebedintsev AI, Slobodin SB, Malyarchuk BA, Martiniano R, Meldgaard M, Arppe L, Palo JU, Sundell T, Mannermaa K, Putkonen M, Alexandersen V, Primeau C, Baimukhanov N, Malhi RS, Sj?gren KG, Kristiansen K, Wessman A, Sajantila A, Lahr MM, Durbin R, Nielsen R, Meltzer DJ, Excoffier L, Willerslev E. The population history of northeastern Siberia since the Pleistocene., 2019, 570(7760): 182–188.

[13] Yue JP, Li YQ, Yang SX. Neolithisation in the southern Lesser Khingan Mountains: lithic technologies and ecological adaptation., 2019, 93(371): 1144–1160.

[14] Tan J, Yang Y, Tang K, Sabeti PC, Jin L, Wang S. The adaptive variant EDARV370A is associated with straight hair in East Asians., 2013, 132(10): 1187–1191.

[15] Fujimoto A, Ohashi J, Nishida N, Miyagawa T, Morishita Y, Tsunoda T, Kimura R, Tokunaga K. A replication study confirmed the EDAR gene to be a major contributor to population differentiation regarding head hair thickness in Asia., 2008, 124(2): 179–185.

[16] Hlusko LJ, Carlson JP, Chaplin G, Elias SA, Hoffecker JF, Huffman M, Jablonski NG, Monson TA, O'Rourke DH, Pilloud MA, Scott GR. Environmental selection during the last ice age on the mother-to-infant transmission of vitamin D and fatty acids through breast milk., 2018, 115(19): E4426–E4432.

[17] Kamberov YG, Wang S, Tan J, Gerbault P, Wark A, Tan L, Yang Y, Li S, Tang K, Chen H, Powell A, Itan Y, Fuller D, Lohmueller J, Mao J, Schachar A, Paymer M, Hostetter E, Byrne E, Burnett M, McMahon AP, Thomas MG, Lieberman DE, Jin L, Tabin CJ, Morgan BA, Sabeti PC. Modeling recent human evolution in mice by expression of a selected EDAR variant., 2013, 152(4): 691–702.

URI: https://kns.cnki.net/kcms/detail/11.1913.R.20210528.1013.002.html

2021-05-27

毛曉偉，博士，副研究員，研究方向：群體遺傳、數量遺傳。E-mail: maoxiaowei@ivpp.ac.cn

平婉菁，碩士，中級工程師，研究方向：實驗技術。E-mail: pingwanjing@ivpp.ac.cn

毛曉偉和平婉菁并列第一作者。

付巧妹，博士，研究員，研究方向: 演化遺傳、群體遺傳。E-mail: fuqiaomei@ivpp.ac.cn

10.16288/j.yczz.21-187

2021-05-28 16:32:40