在東亞與東南亞的交匯處，華南地區(qū)一直扮演著連接南北的關(guān)鍵走廊角色，這里不僅是現(xiàn)代解剖學(xué)意義上人類早期擴(kuò)散的樞紐，更是水稻馴化與南方農(nóng)耕文明的搖籃。然而，熱帶、亞熱帶地區(qū)炎熱潮濕的氣候嚴(yán)重限制了古代DNA的保存與獲取，使得我們對(duì)這一關(guān)鍵區(qū)域長(zhǎng)期、深層次的人口歷史過(guò)程認(rèn)識(shí)存在顯著缺口。特別是，從父系遺傳角度，中國(guó)古代粟作與稻作農(nóng)民對(duì)東南亞人群的遺傳貢獻(xiàn)規(guī)模、時(shí)間線以及具體的人口動(dòng)力學(xué)模式，仍有許多未解之謎。為了填補(bǔ)這一空白，一項(xiàng)發(fā)表在《SCIENCE ADVANCES》上的研究應(yīng)運(yùn)而生，它通過(guò)迄今最大規(guī)模的中國(guó)Y染色體變異數(shù)據(jù)集，為我們揭開了新石器時(shí)代農(nóng)民多次南下的神秘面紗。

研究人員為開展此項(xiàng)研究，主要運(yùn)用了以下幾個(gè)關(guān)鍵技術(shù)方法：首先，他們從中國(guó)53個(gè)民族語(yǔ)言和地理多樣性人群中采集樣本，利用DNBSEQ-T7平臺(tái)進(jìn)行了Y染色體靶向測(cè)序，獲得了463個(gè)高覆蓋度Y染色體序列，并整合了已發(fā)表的1044個(gè)序列、1748個(gè)東南亞低覆蓋度序列以及780個(gè)來(lái)自歐亞大陸東部的古代基因組數(shù)據(jù)，構(gòu)建了包含1507個(gè)個(gè)體的核心數(shù)據(jù)集。其次，他們采用了嚴(yán)格的生物信息學(xué)流程進(jìn)行數(shù)據(jù)質(zhì)控、比對(duì)和變異檢測(cè)。再者，研究運(yùn)用了最大似然法(RAxML-NG)和貝葉斯方法(BEAST)構(gòu)建高分辨率、時(shí)間標(biāo)定的系統(tǒng)發(fā)育樹，以解析譜系關(guān)系和推斷分化時(shí)間。此外，還通過(guò)貝葉斯天際線圖(BSP)分析了關(guān)鍵單倍群的有效種群大小變化歷史，并通過(guò)空間統(tǒng)計(jì)分析、單倍型共享分析和網(wǎng)絡(luò)分析等手段，探究了譜系的地理分布與人群間的遺傳聯(lián)系。

為探究由中國(guó)向東南亞多次南遷塑造的父系遺傳景觀，研究團(tuán)隊(duì)生成了來(lái)自華南53個(gè)群體的463個(gè)高覆蓋度Y染色體序列。通過(guò)整合多個(gè)數(shù)據(jù)集，分析顯示東亞人群的父系遺傳多樣性高于東南亞人群。 neutrality 檢驗(yàn)表明，除東亞南亞語(yǔ)系人群外，大多數(shù)語(yǔ)言多樣性人群經(jīng)歷了近期由遷移驅(qū)動(dòng)的擴(kuò)張。值得注意的是，東南亞部分侗臺(tái)語(yǔ)人群顯示出比東亞同類人群更高的多樣性，這可能反映了來(lái)自東亞的多波遷移以及與當(dāng)?shù)厝巳旱幕旌�。在新測(cè)序個(gè)體中，O2a、O1b和O1a單倍群占主導(dǎo)，這些譜系在苗瑤和侗臺(tái)語(yǔ)人群中尤為突出，凸顯了它們?cè)诟赶底V系南向擴(kuò)散中的核心作用。

通過(guò)對(duì)古今個(gè)體Y染色體系統(tǒng)發(fā)育樹的高分辨率分析，研究揭示了東亞和東南亞詳細(xì)的父系遺傳圖譜。O1和O2是兩個(gè)主導(dǎo)譜系。O1a譜系廣泛分布于中國(guó)西南和沿海地區(qū)，與百越、良渚文化緊密相關(guān)，并在侗臺(tái)、苗瑤和藏緬人群中均有出現(xiàn)。古代證據(jù)將O1a與粟作和稻作農(nóng)業(yè)早期發(fā)展相關(guān)的人群聯(lián)系起來(lái)。O1b1a譜系地理和民族語(yǔ)言分布更廣，主要見于華南和東南亞的侗臺(tái)、南亞和藏緬人群中。O2譜系中的O2a2a1a2亞支在漢族和苗瑤人群中頻率較高，其中O2a2a1a2-M7在廣西獨(dú)山的古代個(gè)體中可追溯到約8707年前，是該父系分支最早記錄之一。O2a2b在侗臺(tái)和藏緬人群中常見，其衍生支系O2a2b1a1a1a4a在現(xiàn)代藏緬人群中頻率很高。

通過(guò)整合古代和現(xiàn)代基因組數(shù)據(jù)，研究證明了多個(gè)語(yǔ)言各異的東亞人群共同貢獻(xiàn)了東南亞人群的父系血統(tǒng)，具有復(fù)雜的多源性。特別是，在長(zhǎng)江中游大溪文化中發(fā)現(xiàn)的O2a2a1a2-M7單倍群，在現(xiàn)代苗瑤、南方漢族以及東亞和東南亞的侗臺(tái)語(yǔ)人群中占主導(dǎo)。研究還發(fā)現(xiàn)，通常與古代東北亞人和現(xiàn)代藏緬或漢語(yǔ)人群相關(guān)的N和D單倍群的上游分支也出現(xiàn)在東南亞古代基因組中，支持了北部漢藏語(yǔ)人群對(duì)東南亞父系基因庫(kù)的影響。與南島語(yǔ)相關(guān)的單倍群O2a2b2-N6在大汶口和曇石山古代個(gè)體中均有發(fā)現(xiàn)，在現(xiàn)代南島語(yǔ)人群中普遍存在，表明中國(guó)新石器時(shí)代南北農(nóng)民與當(dāng)代東南亞南島語(yǔ)群體之間存在直接遺傳聯(lián)系。O1a相關(guān)譜系在侗臺(tái)語(yǔ)人群中廣泛分布，反映了農(nóng)業(yè)時(shí)期侗臺(tái)語(yǔ)人群向東南亞大陸的擴(kuò)張。

為了更好地闡明東亞祖先人群的擴(kuò)張及其對(duì)華南和東南亞稻作社區(qū)形成的時(shí)間貢獻(xiàn)，研究人員基于衍生變異的積累速率重建了五個(gè)時(shí)間標(biāo)定的Y染色體系統(tǒng)發(fā)育樹。這些包括迄今最全面的東亞系統(tǒng)發(fā)育樹。BEAST推導(dǎo)的系統(tǒng)發(fā)育拓?fù)浣Y(jié)構(gòu)和推斷的擴(kuò)張事件與之前基于最大似然法的重建結(jié)果高度一致。值得注意的是，在新報(bào)告的數(shù)據(jù)集中，研究檢測(cè)到單倍群O1a-M119、O1b1a1a-M95、O2a2a1a2a1-F1276和O2a2b1a1a1-CTS4111的擴(kuò)張信號(hào)。所有這些主要擴(kuò)張都可以追溯到新石器時(shí)代轉(zhuǎn)型的開始，并在新石器時(shí)代中晚期加劇，對(duì)侗臺(tái)、苗瑤和漢藏等語(yǔ)言多樣性群體的父系譜系做出了實(shí)質(zhì)性貢獻(xiàn)。

為了闡明關(guān)鍵主導(dǎo)Y染色體單倍群的種群歷史，研究構(gòu)建了貝葉斯天際線圖。單倍群O2a2a1a2-M7在大約4500年前開始顯著的人口擴(kuò)張，在3500年前左右達(dá)到頂峰，這一時(shí)期與從新石器時(shí)代晚期到青銅時(shí)代早期的過(guò)渡期吻合。之后該譜系經(jīng)歷了一個(gè)明顯的瓶頸期。相比之下，單倍群O2a2b1a1a1a4-Z25921和O1b1a1a1b-F789分別在大約10000和13000年前開始持續(xù)擴(kuò)張。單倍群O1b1a1a1a-F1252在3000至2000年前之間顯著擴(kuò)張，而主要沿海分布的O1a-M119在4100至3000年前之間經(jīng)歷了初期人口擴(kuò)張。

研究進(jìn)一步通過(guò)來(lái)自華西生物樣本庫(kù)的36871個(gè)微陣列基因型，分析了關(guān)鍵Y染色體單倍群在中國(guó)和東南亞的地理分布。皮爾遜相關(guān)性分析顯示，這些關(guān)鍵單倍群與緯度呈顯著負(fù)相關(guān)，表明存在明顯的北向南梯度。單倍群O2a2a1a2-M7、O2a2b1a1a1a4-Z25921、O1a-M119、O1b1a1a1a-F1252和O1b1a1a1b-F789主要聚集在華南和東南亞。這種模式揭示了強(qiáng)烈的空間關(guān)聯(lián)性，表明這些父系譜系極大地塑造了該地區(qū)早期人群的遺傳結(jié)構(gòu)。研究還統(tǒng)計(jì)估計(jì)了高分辨率Y染色體系統(tǒng)發(fā)育樹中節(jié)點(diǎn)出現(xiàn)和譜系積累的速率。譜系積累速率在末次盛冰期前后較低，但在大約4000年前的青銅時(shí)代急劇上升。研究發(fā)現(xiàn)在3500至4000年前之間有約230個(gè)譜系出現(xiàn)，表明早期中華文明的崛起及相關(guān)文化創(chuàng)新導(dǎo)致了顯著的人口增長(zhǎng)。

為了闡明東南亞人群與不同民族語(yǔ)言的中國(guó)人群之間的遺傳關(guān)系，包括他們與不同地理來(lái)源的古代人群的聯(lián)系，研究分析了群體水平的單倍群共享模式，并重建了關(guān)鍵奠基譜系的網(wǎng)絡(luò)分析。單倍群共享模式揭示了東亞和東南亞人群之間顯著的人群分層。在第二和第三級(jí)單倍群分辨率上出現(xiàn)了精細(xì)的亞結(jié)構(gòu)，反映了民族語(yǔ)言和地理分化。值得注意的是，來(lái)自東亞和東南亞的侗臺(tái)、苗瑤和某些藏緬語(yǔ)人群形成了獨(dú)特的遺傳簇，凸顯了它們更緊密的遺傳關(guān)系和各自語(yǔ)系內(nèi)可能的共享血統(tǒng)。網(wǎng)絡(luò)分析進(jìn)一步識(shí)別了與O1a、O1b1a1a1a1a1a、O1b1a1a1a1b、O1b1a1a1b、O2a2a和O2a2b1a1等多個(gè)譜系相關(guān)的擴(kuò)張事件。譜系O1a1a2a1a3呈現(xiàn)出層次化的網(wǎng)絡(luò)結(jié)構(gòu)：來(lái)自東亞的侗臺(tái)語(yǔ)人群聚集在上游譜系O1a1a2a附近，而來(lái)自東南亞的侗臺(tái)語(yǔ)人群則位于更衍生的分支上，暗示了從華南向東南亞的南遷軌跡。類似地，單倍群O2a2a1a2a1a1a2a1a1a1呈現(xiàn)出幾乎完全由苗瑤語(yǔ)人群組成的“年輪”模式，反映了從東亞到東南亞的遺傳連續(xù)性。

通過(guò)整合深時(shí)父系遺傳重建與最新的考古學(xué)、語(yǔ)言學(xué)和古DNA證據(jù)，研究揭示了東亞與東南亞交匯處精細(xì)且復(fù)雜的人口歷史。研究結(jié)果支持了至少五次不同的南遷波次的存在，對(duì)應(yīng)于東亞南部和東南亞人類形成的“兩層模型”中的第二層，主要由中國(guó)稻作和粟作人群的擴(kuò)張驅(qū)動(dòng)。其中，最早的一波與南亞語(yǔ)系的擴(kuò)散一致，隨后的兩波與南島語(yǔ)和侗臺(tái)語(yǔ)群體的擴(kuò)散相關(guān)，而最近的遷移則與苗瑤語(yǔ)和漢藏語(yǔ)人群的擴(kuò)張有關(guān)。

這項(xiàng)研究通過(guò)整合來(lái)自38378名無(wú)關(guān)中國(guó)個(gè)體的大規(guī)模Y染色體測(cè)序和基因分型數(shù)據(jù)、1748個(gè)東南亞低覆蓋度基因組以及780個(gè)來(lái)自歐亞大陸東部的古代Y染色體序列，提供了華南和東南亞父系遺傳結(jié)構(gòu)和人口歷史的詳細(xì)視圖。綜合系統(tǒng)發(fā)育分析揭示了中國(guó)新石器農(nóng)業(yè)擴(kuò)張背景下，人類向東南亞遷移的多波次浪潮。這些發(fā)現(xiàn)與先前的常染色體和線粒體DNA研究一致，但提供了父系視角的補(bǔ)充證據(jù)，并在分化與擴(kuò)張時(shí)間、方向性以及文化驅(qū)動(dòng)力（特別是農(nóng)業(yè)創(chuàng)新和語(yǔ)言多樣化）方面提供了更高的分辨率。

高分辨率Y染色體系統(tǒng)發(fā)育樹精煉了“兩層模型”，明確了關(guān)鍵的父系譜系O1a、O1b和O2a2a1a2，它們將東亞古代粟作和稻作人群與當(dāng)今東南亞人群聯(lián)系起來(lái)。雖然先前的常染色體分析揭示了古代人群的雙重血統(tǒng)成分，但Y染色體數(shù)據(jù)揭示了這些古代農(nóng)民與東南亞現(xiàn)代侗臺(tái)、苗瑤和南亞語(yǔ)人群之間直接的父系連續(xù)性，表明來(lái)自中國(guó)的顯著父系基因流動(dòng)。相比之下，來(lái)自云南和長(zhǎng)江流域的古代線粒體DNA顯示出母系連續(xù)性和有限的女性介導(dǎo)的基因流動(dòng)，支持了性別偏向性遷移模型。

關(guān)鍵父系單倍群與民族語(yǔ)言歸屬之間的強(qiáng)對(duì)應(yīng)關(guān)系，為南島語(yǔ)、南亞語(yǔ)、侗臺(tái)語(yǔ)、漢藏語(yǔ)和苗瑤語(yǔ)等語(yǔ)系擴(kuò)散的人口機(jī)制提供了令人信服的遺傳證據(jù)。研究結(jié)果表明，侗臺(tái)人群主要攜帶O1a和O1b譜系，苗瑤人群富含O2a2a1a2及其亞支，南亞語(yǔ)人群則常攜帶O1b1a1a1a和O1b1a1a1b譜系，南島語(yǔ)相關(guān)譜系則支持了從華南沿海出發(fā)的新石器時(shí)代航海擴(kuò)張模型。此外，在古代華南個(gè)體中檢測(cè)到的古代東北亞譜系，支持了與藏緬語(yǔ)相關(guān)的南向遷移假說(shuō)。

本研究的一個(gè)關(guān)鍵發(fā)現(xiàn)是，主要Y染色體單倍群的擴(kuò)張時(shí)間與粟作和稻作農(nóng)業(yè)的出現(xiàn)相吻合。時(shí)間標(biāo)定的Y染色體系統(tǒng)發(fā)育樹和古今基因組的綜合分析表明，主要擴(kuò)張，特別是O2a2a1a2-M7、O1a-M119和O1b1a1a1b，發(fā)生在大約4000至3500年前。這一時(shí)期與長(zhǎng)江中下游地區(qū)農(nóng)業(yè)活動(dòng)的強(qiáng)化、社會(huì)分層以及早期國(guó)家的出現(xiàn)相吻合。地理分析顯示單倍群頻率隨緯度升高而降低，支持了人口從東北亞向南遷移的定向模型。東南亞人群，特別是侗臺(tái)和苗瑤語(yǔ)人群基因組多樣性的增加，可歸因于農(nóng)業(yè)擴(kuò)張過(guò)程中的系列奠基者效應(yīng)、本地混合和生態(tài)適應(yīng)。

總之，本研究通過(guò)整合來(lái)自不同民族語(yǔ)言的現(xiàn)代群體的高分辨率Y染色體數(shù)據(jù)和可追溯到全新世的古代基因組，重建了華南和東南亞的父系人口歷史。由此產(chǎn)生的系統(tǒng)發(fā)育樹和人口模型表明，父系譜系的擴(kuò)張與農(nóng)業(yè)的發(fā)展、社會(huì)復(fù)雜性的興起以及歐亞大陸東部早期國(guó)家的出現(xiàn)相吻合。來(lái)自東亞的多波次男性主導(dǎo)的遷移，主要與新石器時(shí)代轉(zhuǎn)型相關(guān)，塑造了東南亞的父系遺傳景觀。O1a、O1b和O2a2a1a2等譜系在當(dāng)今侗臺(tái)、苗瑤、南亞和南島語(yǔ)人群的形成中扮演了核心角色，其祖先根源可追溯至長(zhǎng)江和黃河流域。這些遷移是結(jié)構(gòu)化且反復(fù)發(fā)生的，而非均一的。譜系地理學(xué)和基于網(wǎng)絡(luò)的分析揭示了早期東亞農(nóng)民與現(xiàn)代東南亞群體之間持久的父系聯(lián)系，強(qiáng)調(diào)了遺傳連續(xù)性而非替代。這些發(fā)現(xiàn)共同挑戰(zhàn)了過(guò)度簡(jiǎn)化的人口替代模型，轉(zhuǎn)而支持一個(gè)以分層遷移、混合和文化轉(zhuǎn)化為特征的復(fù)雜人口歷史。