• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    真核生物起源研究進(jìn)展

    2020-10-29 05:37:18高志偉王龍
    遺傳 2020年10期
    關(guān)鍵詞:真核宿主線粒體

    高志偉,王龍

    真核生物起源研究進(jìn)展

    高志偉1,王龍2

    1. 南京大學(xué)匡亞明學(xué)院,南京 210023 2. 南京大學(xué)生命科學(xué)學(xué)院,醫(yī)藥生物技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,南京 210023

    作為重大進(jìn)化謎題,真核生物起源的研究對于解碼真核基因組、闡釋真核細(xì)胞內(nèi)部結(jié)構(gòu)之間的關(guān)系有重要啟示作用。在1977年美國微生物學(xué)家Carl Woese發(fā)現(xiàn)古細(xì)菌并提出三域生命之樹之后,大量研究顯示古細(xì)菌與真核生物在進(jìn)化上存在著密切聯(lián)系。21世紀(jì)以來,系統(tǒng)發(fā)育分析方法不斷改進(jìn),泉古菌門(Crenarchaeota)、廣古菌門(Euryarchaeota)之外與真核生物更加相似的新古細(xì)菌門類也相繼被發(fā)現(xiàn),這些證據(jù)更加支持將真核生物與古細(xì)菌合并為一域,形成二域生命之樹。目前,通過宏基因組技術(shù)發(fā)現(xiàn)的Asgard古細(xì)菌是與真核生物進(jìn)化距離最近的原核生物。然而,真核生物祖先的身份以及線粒體起源的時(shí)間等核心問題仍是學(xué)術(shù)界爭論的焦點(diǎn)。本文結(jié)合近年來國內(nèi)外研究成果,從生命之樹的形態(tài)變化與真核生物演變的具體機(jī)制兩個(gè)角度梳理了目前對真核生物起源的認(rèn)知過程、現(xiàn)有水平和研究前景,以期為揭示真核生物起源進(jìn)程的后續(xù)研究提供參考與指引。

    真核生物起源;生命之樹;古細(xì)菌;三域;二域;內(nèi)共生;線粒體;Asgard超門

    真核生物起源被稱為生物界僅次于生命起源的第2大謎團(tuán)[1]。加拿大著名微生物學(xué)家Roger Stanier曾指出:“原核生物與真核生物細(xì)胞之間的差異是當(dāng)今世界中已發(fā)現(xiàn)的最大的單一不連續(xù)性”[2]。雖然許多原核生物已經(jīng)具有帶有膜的內(nèi)部組織結(jié)構(gòu)[3],但真核生物的復(fù)雜程度卻遠(yuǎn)高于此。真核生物具有核膜包被的細(xì)胞核,將轉(zhuǎn)錄與翻譯過程分隔開,為基因調(diào)控提供更多選擇[4];眾多細(xì)胞器組成的內(nèi)膜系統(tǒng)使細(xì)胞具有精密高效的運(yùn)輸和代謝功能[5]。這些結(jié)構(gòu)增強(qiáng)了真核生物的環(huán)境適應(yīng)能力,使其演化為現(xiàn)存物種中數(shù)目占大多數(shù)的生命形式。

    從原核到真核這一巨大轉(zhuǎn)變通常在教科書中被過度簡化為一個(gè)迅速的單一事件[6]。然而從時(shí)間上看,化石記錄顯示目前最早的真細(xì)菌[7,8]和古細(xì)菌化石距今約35億年,相比之下,真核生物化石距今約18億年[9],與利用分子鐘計(jì)算的17~19億年時(shí)間相一致[10],由此推算出的原核生物和真核生物的出現(xiàn)時(shí)間相差17億年之久。闡釋這一漫長歷史過程中真核生物的演變歷程對人們認(rèn)知真核生物乃至人類自身的細(xì)胞和基因組結(jié)構(gòu)至關(guān)重要。

    真核生物的祖先是什么?真核生物演變的具體機(jī)制是什么?這兩大核心問題也成為學(xué)界長久以來爭論的焦點(diǎn)。其中,尤其是真核生物在生命之樹上所處的關(guān)鍵節(jié)點(diǎn)位置,以及線粒體的起源問題上,依然存在有很大的分歧。近年來,與真核生物進(jìn)化密切相關(guān)的Asgard古細(xì)菌超門的發(fā)現(xiàn)[11,12]以及其在實(shí)驗(yàn)室條件下培養(yǎng)研究的實(shí)現(xiàn)[13],為解決這一僵局帶來了光明。因此,本文在回顧歷史的基礎(chǔ)上,重點(diǎn)聚焦近年來在該領(lǐng)域取得的重大研究進(jìn)展,從生命之樹和演化機(jī)制兩個(gè)角度深入闡述真核生物的起源問題。

    1 生命之樹中的真核生物

    1.1 古細(xì)菌的發(fā)現(xiàn)和三域生命之樹的提出

    1977年,美國微生物學(xué)家Carl Woese通過16S/18S rRNA序列系統(tǒng)發(fā)育分析偶然發(fā)現(xiàn)了古細(xì)菌。他們據(jù)此將原核生物分為兩大類——真細(xì)菌和古細(xì)菌,并將真細(xì)菌、古細(xì)菌和真核生物各分為一個(gè)“域”(Domain),作為比“界”(Kingdom)高一級的分類系統(tǒng),分別命名為細(xì)菌域(Bacteria)、古菌域(Archaea)和真核域(Eukarya)[14]。自此,真核生物在分類學(xué)中的位置迎來巨變,利用保守基因?yàn)樗屑?xì)胞生命建立分類系統(tǒng)的方法開始流行。20世紀(jì)80年代末,基于基因重復(fù)(gene duplication)[15]和SSU rRNA[16]的研究發(fā)現(xiàn)真核生物雖與古細(xì)菌形態(tài)差異大,但屬于姐妹群(sister groups),具有一個(gè)獨(dú)立于真細(xì)菌之外的共同祖先。據(jù)此,Woese等[16]提出了著名的三域(three primary domains, 3D)生命之樹(tree of life/ universal tree),包含細(xì)菌域、古菌域(分為泉古菌門和廣古菌門)和真核域(圖1)[17],長期以來得到各類教科書和文獻(xiàn)的廣泛認(rèn)可[18]。

    20世紀(jì)80年代到21世紀(jì)初,許多研究發(fā)現(xiàn)古細(xì)菌和真核生物在分子層面上高度相關(guān),直接或間接地為Carl Woese的三域生命之樹提供證明。在遺傳信息復(fù)制、轉(zhuǎn)錄和翻譯水平上,真核生物與古細(xì)菌具有許多真細(xì)菌不具有的相似性。例如,免疫學(xué)實(shí)驗(yàn)[19]和進(jìn)化計(jì)算[20]都證明古細(xì)菌依賴于DNA的RNA聚合酶與真核生物的結(jié)構(gòu)同源性。真核生物與古細(xì)菌的啟動(dòng)子和轉(zhuǎn)錄起始因子、30種共有的核糖體蛋白、翻譯因子、DNA復(fù)制相關(guān)的部分酶、組蛋白等都體現(xiàn)出真細(xì)菌所缺乏的相似性[21]。此外,真核生物和古細(xì)菌具有多種相似的、維持細(xì)胞基本生存的酶和代謝通路,如液泡ATPase、分泌通路、輔酶A生物合成的部分酶以及嘧啶和精氨酸生物合成的部分酶等[21]。古細(xì)菌細(xì)胞分裂系統(tǒng)內(nèi)的部分蛋白甚至與真核生物的內(nèi)體分選復(fù)合體同源[22,23]。如此多的進(jìn)化關(guān)聯(lián)性使古細(xì)菌成為揭示真核生物起源的重要研究對象[24]。

    當(dāng)然,古細(xì)菌和真核生物之間也存在許多差別。最為顯著的是膜脂,古細(xì)菌采用醚鍵相連,而真核生物采用酯鍵相連;古細(xì)菌的甘油磷酸骨架使用甘油-1-磷酸,而真核生物和真細(xì)菌使用甘油-3-磷酸,兩者由不相關(guān)的酶合成。古細(xì)菌還擁有獨(dú)特的產(chǎn)甲烷代謝途徑、DNA拓?fù)洚悩?gòu)酶和DNA聚合酶等[21]。此外,真核生物與真細(xì)菌也有代謝[25,26]和膜化學(xué)性質(zhì)[27]的部分相似之處。

    圖1 Carl Woese的生命之樹

    Carl Woese依據(jù)16S/18S rRNA系統(tǒng)發(fā)育分析發(fā)現(xiàn)古細(xì)菌后,將自然界中的生物分為真核域、細(xì)菌域和古菌域三域。根據(jù)文獻(xiàn)[17]修改繪制。

    1.2 二域與三域的爭端

    雖然Carl Woese提出的三域系統(tǒng)(圖2A)受到廣泛認(rèn)可,但同時(shí)期,Lake等[28]依據(jù)核糖體結(jié)構(gòu)證據(jù)、Rivera和Lake[29]依據(jù)延伸因子EF-1α同系物保守區(qū)有共同氨基酸插入的證據(jù),提出真核生物是泉古菌門的一個(gè)姐妹群,即真核生物是從古細(xì)菌域之內(nèi)產(chǎn)生的,是二級域,而一級域僅為真細(xì)菌和古細(xì)菌,形成二域(two primary domains, 2D)生命之樹(圖2B)。人們對應(yīng)地提出了許多演化模型支持二域系統(tǒng),實(shí)質(zhì)均為古細(xì)菌與細(xì)菌之間的相互作用[30],如支持泉古菌門(Crenarchaeota)與細(xì)菌相互作用的泉古菌假說(eocyte hypothesis)[28,29],廣古菌門(Euryarchaeota)與細(xì)菌相互作用的氫假說(hydrogen hypothesis)[31]和互養(yǎng)共棲假說(syntrophic hypo-thesis)[32]等。

    21世紀(jì)的第一個(gè)10年中,生物信息學(xué)蓬勃發(fā)展。從系統(tǒng)發(fā)育重建的角度看,三域與二域的爭端依然激烈。在2003~2009年發(fā)表的7篇文獻(xiàn)中,通過真核生物、真細(xì)菌、古細(xì)菌物種取樣,測定基因組中的系統(tǒng)發(fā)生信號(phylogenetic signal)標(biāo)記,進(jìn)行大規(guī)模系統(tǒng)發(fā)育學(xué)分析,3篇文獻(xiàn)得出三域系統(tǒng)結(jié)果[33~35],4篇文獻(xiàn)得出二域系統(tǒng)結(jié)果[36~39]。7篇文獻(xiàn)采用的數(shù)據(jù)和方法比較如表1所示。

    值得注意的是,上述7篇文獻(xiàn)使用數(shù)據(jù)在很大程度上有重疊,大多基于蛋白相鄰類的聚簇(Cluster of Orthologous Groups of proteins, COG)數(shù)據(jù)庫,最終卻得到了截然不同的結(jié)論,集中體現(xiàn)了系統(tǒng)發(fā)育學(xué)方法在解決真核生物起源時(shí)面臨的典型難題:

    圖2 三域與二域生命之樹的對比

    A:三域系統(tǒng);B:二域系統(tǒng)。AB圖中的1、2、3分別對應(yīng)真核域、古菌域和細(xì)菌域的最近共同祖先;真核域與古菌域是兩個(gè)距離遙遠(yuǎn)的姐妹群,具有一個(gè)最近共同祖先4。B圖中真核域的最近共同祖先1直接產(chǎn)生自古菌域內(nèi),所以古菌域的最近共同祖先2同時(shí)也是古菌域和真核域的最近共同祖先,這與A圖的情況不同。根據(jù)文獻(xiàn)[21]修改繪制。

    表1 2003~2009年間發(fā)表的7篇闡釋生命之樹中古細(xì)菌與真核生物關(guān)系的文章比較

    3D:三域系統(tǒng);2D:二域系統(tǒng);HGT:水平基因轉(zhuǎn)移。表格根據(jù)Gribaldo等[21]的綜述拓展。

    (1)長期的堿基替換導(dǎo)致基因組中遺留的系統(tǒng)發(fā)育信號很微弱,若使用單基因位點(diǎn)或較短的16S rRNA、18S rRNA分析跨域的物種,背景噪聲大,可以得到的有效信息很少,如前述Harris等[33]與Yutin等[35]兩文的方法正如此。同時(shí),單個(gè)基因可能受基因特異性事件影響(如基因重復(fù)、丟失和轉(zhuǎn)移),不能反映真實(shí)的物種進(jìn)化情況[40,41]。

    (2)為解決單基因信號微弱的問題,許多研究采用了分散基因串聯(lián)(concatenation)的方法[34,38,39],即將多個(gè)不連續(xù)的基因直接串聯(lián)組合為“普遍基因”(universal gene),序列對齊后再統(tǒng)一分析。但此舉的局限性在于真核生物起源分析的范圍廣,真核生物、古細(xì)菌和真細(xì)菌3個(gè)類群跨域分析時(shí)很可能將非同源蛋白強(qiáng)制對齊,在組合基因時(shí)人為引入偏差[21]。此外,可供大范圍跨域分析的“普遍基因”本身數(shù)量也十分有限。

    (3)同樣為解決單基因信號微弱的問題,有研究采用全基因組比對,但局限性在于可用的全基因組數(shù)據(jù)比單基因數(shù)據(jù)少,且計(jì)算更加復(fù)雜。選取物種數(shù)少時(shí),結(jié)果受抽樣本身影響大,例如Rivera和Lake[36]僅僅使用了7個(gè)全基因組,偶然誤差較大。但該方法隨著未來基因組數(shù)據(jù)的快速增加和組學(xué)分析方法的不斷進(jìn)步,有望提供更多可靠的信息。

    (4)不同物種間的水平基因轉(zhuǎn)移(horizontal gene transfer, HGT)和基因丟失都可能抹去進(jìn)化印記。HGT和基因丟失可能發(fā)生在真核生物進(jìn)化的各個(gè)階段(圖3),機(jī)制包括轉(zhuǎn)化、轉(zhuǎn)導(dǎo)和接合等[42~44]。在Cox等[38]的研究中,檢測到選取的80個(gè)基因有30個(gè)受到HGT影響,數(shù)據(jù)舍棄不用。Akanni等[45]利用超樹方法分析了392個(gè)真細(xì)菌和51個(gè)古細(xì)菌基因組,在古細(xì)菌進(jìn)化樹基部發(fā)現(xiàn)了兩個(gè)來自真細(xì)菌的大規(guī)模HGT,基因來源涉及δ-變形桿菌、梭菌和放線菌。對此,Cotton等[46]和Cohen等[47]認(rèn)為,與轉(zhuǎn)錄、翻譯相關(guān)的核心基因的負(fù)選擇壓大,受HGT的影響應(yīng)比一般代謝通路的基因小,是更加合適的分析對象。

    (5)長支吸引(long branch attraction, LBA),由于不同域間的基因差別可能很大,在同源性高的物種之外可能形成進(jìn)化距離遙遠(yuǎn)但仍然相互靠近的長支,造成偏差[48,49]。早在1999年,Tourasse等[50]便認(rèn)為三域樹結(jié)果是受LBA影響產(chǎn)生的。

    (6)建立的模型本身設(shè)計(jì)過于簡單,模擬能力差[48]。近10年來,隨著模型的復(fù)雜化和基因數(shù)據(jù)的積累,人們盡量規(guī)避上述6點(diǎn)對結(jié)果的影響。2011~ 2013年發(fā)表的3篇文獻(xiàn)結(jié)果都支持二域系統(tǒng)[48,51,52]。2014年,Rochette等[53]的系統(tǒng)發(fā)育分析顯示,真核生物與古細(xì)菌間的聯(lián)系比真細(xì)菌緊密得多,但與泉古菌門和廣古菌門沒有密切聯(lián)系,真核生物分支接近古細(xì)菌并很可能從中發(fā)出。與此同時(shí),泉古菌和廣古菌以外古細(xì)菌門類的發(fā)現(xiàn)為二域系統(tǒng)提供了更多的支持,相關(guān)內(nèi)容將在下一部分闡述。

    圖3 真核生物進(jìn)化過程中的水平基因轉(zhuǎn)移(HGT)與基因丟失

    HGT可能發(fā)生在古細(xì)菌祖先→FECA→LECA→現(xiàn)存真核生物演化進(jìn)程中的各個(gè)階段。線粒體和葉綠體的內(nèi)共生過程也為HGT的發(fā)生提供條件。HGT和基因丟失均能削弱現(xiàn)存生物基因組中遺留的系統(tǒng)發(fā)育信號。根據(jù)文獻(xiàn)[17]修改繪制。

    1.3 新古細(xì)菌類群的發(fā)現(xiàn)對生命之樹的影響

    2004年,Venter等[54]首次采用不經(jīng)實(shí)驗(yàn)室培養(yǎng)直接測定環(huán)境中微小生物基因組的“全基因組鳥槍測序”方法,即通常所說的宏基因組(metagenome)技術(shù),研究海水中的微生物樣本,發(fā)現(xiàn)了148個(gè)新細(xì)菌種類和120萬個(gè)此前未知的基因。此后,人們不斷利用宏基因組技術(shù)發(fā)現(xiàn)多個(gè)泉古菌門、廣古菌門之外的古細(xì)菌新門,為生命之樹帶來巨大變化。

    2008年后,奇古菌門(Thaumarchaeota)[55]、初古菌門(Korarchaeota)[56]、曙古菌門(Aigarchaeota)[57]相繼被發(fā)現(xiàn)。其中奇古菌門是被發(fā)現(xiàn)的第一個(gè)新古細(xì)菌門類,而曙古菌門是第一個(gè)完全由宏基因組學(xué)數(shù)據(jù)重建的基因組,并且發(fā)現(xiàn)了與真核生物同源的泛素系統(tǒng)。依據(jù)16S rRNA測序結(jié)果,這3個(gè)新門與泉古菌門一起歸為一個(gè)單系群(monophyletic group),稱為“TACK”超門(TACK superphylum)或變形古菌(Proteoarchaeota)[58,59]。通過改進(jìn)的最大似然估計(jì)或貝葉斯網(wǎng)絡(luò)模型,人們發(fā)現(xiàn)真核生物從TACK超門內(nèi)分支出來或是TACK超門的一個(gè)姐妹群(圖4A)[17,60]。這表明真核生物分支是從古細(xì)菌域內(nèi)發(fā)出的,支持上一部分所敘述的二域系統(tǒng),否定三域系統(tǒng)。新發(fā)現(xiàn)的TACK超家族成員含有大量一般只在真核生物中表達(dá)或與其相關(guān)的真核生物標(biāo)志蛋白(eukaryotic signature proteins, ESPs),如肌動(dòng)蛋白、微管蛋白的古細(xì)菌同源物、轉(zhuǎn)錄翻譯因子類似的蛋白等[11]。TACK超門成為追溯真核細(xì)胞起源的重要研究材料。

    2014年以來,通過改進(jìn)的宏基因組技術(shù)——解碼宏基因組(genome-resolved metagenomics),人們又在TACK超門中發(fā)現(xiàn)新的古細(xì)菌門:深古菌門(Bathyarchaeota)、地古菌門(Geoarchaeota)和韋斯特古菌門(Verstraetearchaeota) (圖4B)[17]。廣古菌門外還存在著一個(gè)深枝DPANN超門(DPANN super-phylum),包含丙鹽古菌門(Diapherotrites)、小古菌門(Parvarchaeota)、謎古菌門(Aenigmarchaeota)、納古菌門(Nanoarchaeota)和納鹽古菌門(Nanohaloar-chaeota) (圖4B)[17]。這些發(fā)現(xiàn)極大地豐富了生命之樹中古細(xì)菌的分支。

    近5年來的一個(gè)重要突破是發(fā)現(xiàn)了與真核生物遺傳關(guān)系更近的Asgard古細(xì)菌。Asgard古細(xì)菌大多生活在海底,分類上原先屬于海洋底棲類群B (Marine Benthic Group B, MBG-B)、古代古細(xì)菌群(Ancient Archaeal Group, AAG)、深海古細(xì)菌群(Deep-sea Archaeal Group, DSAG)和海洋熱液噴口群(Marine Hydrothermal Vent Group)[61]。2015年,瑞典烏普薩拉大學(xué)Ettema課題組在北冰洋Loki’s Castle附近的海底沉積物中利用宏基因組技術(shù)重建出一種新的古細(xì)菌基因組,命名為洛基古菌門(Lokiarchaea),蛋白質(zhì)組中的重要部分(175個(gè)蛋白質(zhì),占比3.3%)都與真核蛋白相似,為ESPs,參與膜的變形和細(xì)胞形狀的形成過程[11]。此后,又發(fā)現(xiàn)了與Lokiarchaea門相近的3個(gè)新門:奧丁古菌門(Odinarchaeota)、索爾古菌門(Thorarchaeota)和海姆達(dá)爾古菌門(Heimdallarchaeota)[12,62]。這4個(gè)新門被一并命名為Asgard超門(Asgard superphylum)[11,12](圖4B),是真核生物的姐妹群。其中,海姆達(dá)爾古菌門與真核生物距離最近[63]。

    圖4 發(fā)現(xiàn)新門類古細(xì)菌后的生命之樹

    A:新發(fā)現(xiàn)TACK超門古細(xì)菌。真核生物是TACK超門的一個(gè)姐妹群,或者直接從TACK超門中產(chǎn)生。B:新發(fā)現(xiàn)DPANN超門和Asgard超門古細(xì)菌。DPANN超門形成古細(xì)菌的一個(gè)深枝,Asgard超門與真核生物進(jìn)化距離最近,真核生物是Asgard超門的一個(gè)姐妹群,或者直接從Asgard超門中產(chǎn)生。根據(jù)文獻(xiàn)[17]修改繪制。

    不過,新的發(fā)現(xiàn)也引發(fā)了學(xué)術(shù)界的爭議。一方面,Da Cunha等[64]指出,Ettema課題組的取樣存在污染,最終結(jié)果也因?yàn)槿〉窖莼俣瓤斓男蛄卸斐蒐BA效應(yīng),并且去除數(shù)據(jù)集中的蛋白延伸因子EF2便足夠破壞真核生物與Lokiarchaea的密切關(guān)系。他們認(rèn)為Lokiarchaea和Asgard超門是廣古菌門的姐妹群,而非真核生物的姐妹群,且真核生物為單系群,根源通往TACK超門中Thaumarchaeota的分支,最終結(jié)果是三域系統(tǒng)。Ettema課題組的Spang等對此進(jìn)行了反駁,認(rèn)為沒有證據(jù)表明取樣污染,去除EF2依然能夠證明Asgard與真核生物屬于一個(gè)單系群,而Da Cunha等使用了不合適的方法導(dǎo)致三域系統(tǒng)結(jié)果[65]。此后雙方仍有爭論,但未有結(jié)果[66,67]。最近,Williams等[49]針對Da Cunha等[64]的結(jié)果,使用相同軟件的不同最大似然方法包分析相同數(shù)據(jù),得到與原文不同的二域系統(tǒng)結(jié)果;他們進(jìn)一步使用多個(gè)新模型分析35個(gè)基因的對齊結(jié)果,真核生物均與Asgard古細(xì)菌形成姐妹群;基于“超樹”和基因合并對齊方法分析3000多個(gè)古細(xì)菌和真核生物基因的結(jié)果顯示真核生物是從古細(xì)菌域內(nèi)產(chǎn)生的,因此認(rèn)為“生命之樹”應(yīng)當(dāng)為二域系統(tǒng)。另一方面,Harish等[68~70]以蛋白質(zhì)折疊作為分析對象,將前人證明三域系統(tǒng)的模型[71,72]擴(kuò)展為非平穩(wěn)模型(non- stationary model),推算出的生命之樹樹根將細(xì)胞生命分為無核生物(Akaryote,即古細(xì)菌與真細(xì)菌)和真核生物二域,形成了一個(gè)與此前二域、三域生命之樹均不同的二域系統(tǒng)(圖5)。分析結(jié)果還指出,Asgard古細(xì)菌與所有其他的古細(xì)菌是姐妹群,難以作為一個(gè)整體處理。這一結(jié)論說明細(xì)胞生命內(nèi)最深層次的分裂位于原核與真核之間。相比常用的蛋白質(zhì)氨基酸序列[49],Harish等[73]認(rèn)為蛋白質(zhì)結(jié)構(gòu)域在生化角度上非冗余和更保守,更適合作為進(jìn)化標(biāo)記進(jìn)行系統(tǒng)發(fā)育學(xué)分析。然而,這個(gè)新的二域生命之樹是選取不同進(jìn)化分子標(biāo)記的結(jié)果,違背了此前眾多研究所證明的真核生物與古細(xì)菌間的密切關(guān)系,因此爭論依然存在。

    對于Asgard古細(xì)菌與真核生物的密切關(guān)系也存在不同解讀。Fournier和Poole[74]便試圖以現(xiàn)代哺乳動(dòng)物的祖先合弓綱(Synapsida)與現(xiàn)代爬行動(dòng)物的進(jìn)化關(guān)系作為平行案例,說明僅從真核生物與Asgard古細(xì)菌的譜系關(guān)系無法解決二域/三域爭論,此外還應(yīng)確定古細(xì)菌和真細(xì)菌之間的根源。

    雖有上述爭論,但發(fā)現(xiàn)新的、進(jìn)化上更接近于真核生物的Asgard種類依然是重要的研究方向。Seitz等[75]發(fā)現(xiàn)了新的Asgard古細(xì)菌海拉古菌門(Helarchaeota)。深圳大學(xué)李猛課題組利用開源數(shù)據(jù)庫中的現(xiàn)有數(shù)據(jù)分析Asgard古細(xì)菌16S rRNA序列,發(fā)現(xiàn)了5個(gè)此前未知的Asgard亞組(subgroup),并且為Asgard古細(xì)菌的轉(zhuǎn)錄活性混合營養(yǎng)型生活方式提供了證據(jù)[76]。他們利用幾個(gè)沿海沉積物樣本重建出15個(gè)基因組,其中包括一個(gè)新的Asgard門類“Gerdarchaeota”[77]。最近,他們設(shè)計(jì)的新SSU rRNA引物完成了Asgard中一種不能與通用引物匹配的新類群的擴(kuò)增工作[78]。Ettema課題組未發(fā)表的結(jié)果也顯示發(fā)現(xiàn)了5個(gè)新Asgard類群,其中Idunnarchaeota是目前與真核生物最接近的類群[79]。

    總之,隨著模型的復(fù)雜化、基因序列和蛋白質(zhì)序列數(shù)據(jù)的積累以及新門類古細(xì)菌的發(fā)現(xiàn),學(xué)術(shù)界對于生命之樹的爭論已逐漸由三域系統(tǒng)向二域系統(tǒng)學(xué)說傾斜。主流觀點(diǎn)認(rèn)為真核生物產(chǎn)生自古細(xì)菌域內(nèi),Asgard超門是目前與其進(jìn)化距離最接近的古細(xì)菌。在未來,生命之樹的形態(tài)以及真核生物的定位仍可能不斷改變,但Carl Woese的三域生命之樹仍將長期存在于教科書與人們的腦海中。

    2 真核生物的起源機(jī)制

    2.1 真核生物起源機(jī)制的研究現(xiàn)狀

    關(guān)于真核生物的起源,首先需要明確的兩個(gè)概念是FECA和LECA。FECA (first eukaryotic common ancestor)指所有真核生物的最早共同祖先;而LECA (last eukaryotic common ancestor)指所有現(xiàn)存真核生物距今最近的共同祖先[17]。FECA與LECA間的階段常被稱為進(jìn)化樹上的莖期(stem phase)[80]。大量原生生物和多細(xì)胞生物的基因組學(xué)和細(xì)胞生物學(xué)研究已證明,LECA擁有現(xiàn)代真核生物的大部分遺傳和細(xì)胞結(jié)構(gòu)特征,如細(xì)胞核、線粒體、細(xì)胞骨架、囊泡和運(yùn)輸系統(tǒng)、細(xì)胞周期等[6]。保守估計(jì)LECA基因組有至少4000~5000個(gè)基因[80,81]。

    圖5 有爭議的新二域生命之樹

    依據(jù)蛋白質(zhì)折疊分析結(jié)果構(gòu)建的生命之樹將生物分為真核生物和無核生物兩域。根據(jù)文獻(xiàn)[73]修改繪制。

    從理論上分析,原核生物向真核生物的轉(zhuǎn)變大致涉及以下方面:由厭氧向兼性厭氧、好氧代謝方式的轉(zhuǎn)變,共生與內(nèi)共生,吞噬作用的產(chǎn)生,基因組嵌合體的變化,細(xì)胞核的出現(xiàn),信號通路的復(fù)雜化,被膜細(xì)胞器的出現(xiàn),細(xì)胞周期的出現(xiàn),細(xì)胞體積的增大等等[1]。研究真核生物起源機(jī)制的關(guān)鍵在于以FECA和LECA為坐標(biāo),理清莖期中這些進(jìn)化事件的先后關(guān)系和具體過程。

    結(jié)合實(shí)驗(yàn)數(shù)據(jù)與理論分析,一些局部的問題已經(jīng)得到充分闡釋,例如GTPase、longins、COPs和ESCRT-III復(fù)合物等真核生物標(biāo)志性蛋白的進(jìn)化歷程,內(nèi)膜系統(tǒng)和核孔復(fù)合物等真核細(xì)胞結(jié)構(gòu)的組建過程已在大量綜述中有詳盡描述[6,82,83]。

    然而,從整體上解決真核生物的起源機(jī)制始終繞不開線粒體的起源問題(圖6)[6,84]。作為細(xì)胞的能量工廠,線粒體為真核生物在“大氧化事件”后在有氧環(huán)境中的輻射進(jìn)化提供了必要條件;真核生物內(nèi)膜系統(tǒng)與細(xì)胞核的產(chǎn)生也與線粒體的內(nèi)共生密切相關(guān)。

    2.2 線粒體與真核生物起源

    線粒體普遍存在于現(xiàn)代真核生物體內(nèi),且隨不同物種展現(xiàn)出豐富的多樣性[85]。許多真核生物(多營寄生生活)為適應(yīng)厭氧環(huán)境,將線粒體簡化為線粒體相關(guān)細(xì)胞器(mitochondrion-related organelles, MROs)[86,87]。極端的例子如sp.的線粒體已完全喪失[88]。但學(xué)界普遍認(rèn)為線粒體的簡化喪失是發(fā)生在LECA之后比較新的進(jìn)化行為,LECA擁有線粒體結(jié)構(gòu)[30,81,89]。

    圖6 線粒體內(nèi)共生是引發(fā)真核生物起源的重要事件

    現(xiàn)有的主流觀點(diǎn)認(rèn)為一種類似于現(xiàn)存Asgard超門的古細(xì)菌宿主內(nèi)共生一種α-變形菌后形成胞內(nèi)的線粒體,原核細(xì)胞逐漸演化為真核細(xì)胞。根據(jù)文獻(xiàn)[6]修改繪制。

    真核生物祖先最初通過內(nèi)共生途徑獲得線粒體。1905年,俄國植物學(xué)家Konstantin Mereschkowsky便提出葉綠體由內(nèi)共生形成的構(gòu)想;20世紀(jì)20年代,Wallin將此推廣至線粒體[90]。40多年后,美國生物學(xué)家Lynn Sagan (后來的Lynn Margulis)重拾內(nèi)共生理論,在著名的一文中提出真核生物的細(xì)胞器包括線粒體和葉綠體在內(nèi)都是由內(nèi)共生細(xì)菌演化而來的[91],引發(fā)巨大爭議。不過,隨后對于線粒體、葉綠體的基因和蛋白的分析顯示它們都擁有與真核譜系相距甚遠(yuǎn)的原核特征,分別對應(yīng)α-變形菌和藍(lán)細(xì)菌,為內(nèi)共生學(xué)說提供了佐證[92~94]。后來在一書中,Margulis正式提出線粒體是原始真核生物吞噬革蘭氏陰性好氧細(xì)菌后,通過長期內(nèi)共生作用形成的[90],逐漸得到學(xué)術(shù)界的廣泛認(rèn)可。自此,線粒體的內(nèi)共生起源成為真核生物起源研究的核心點(diǎn)。

    2.2.1 線粒體內(nèi)共生的早期與晚期之爭

    長期以來,學(xué)界在線粒體內(nèi)共生發(fā)生的時(shí)間問題上存在爭論。主要觀點(diǎn)分為兩派:一派認(rèn)為宿主吞噬好氧的α-變形菌前,已經(jīng)先通過不同途徑形成細(xì)胞核,成為具有真核細(xì)胞特征(含細(xì)胞核、動(dòng)態(tài)細(xì)胞骨架、內(nèi)膜系統(tǒng))的過渡體,具有原始吞噬功能等,被稱為晚期線粒體模型(“Mito-late” models);另一派認(rèn)為,宿主首先與好氧的α-變形菌內(nèi)共生形成帶有線粒體的原核細(xì)胞,隨后促進(jìn)原核細(xì)胞演化出細(xì)胞核、內(nèi)膜系統(tǒng)等真核生物特征,被稱為早期線粒體模型(“Mito-early” models)。

    20世紀(jì)80年代,研究者在現(xiàn)存真核生物中發(fā)現(xiàn)了不含線粒體的物種,如賈第蟲屬()、毛滴蟲屬()等,大多是寄生蟲。英國演化生物學(xué)家Cavalier-Smith將它們歸類為源真核生物(archezoans)[95],祖先是細(xì)胞內(nèi)部已具有復(fù)雜結(jié)構(gòu)而尚未內(nèi)共生線粒體的初級真核生物,屬于真核生物樹中的早期分支譜系[96]。傳統(tǒng)的晚期線粒體假說認(rèn)為,細(xì)胞內(nèi)部的復(fù)雜結(jié)構(gòu)以及吞噬營養(yǎng)能力是線粒體起源的前提,細(xì)胞核的出現(xiàn)是為了避免細(xì)胞骨架牽拉對DNA的損傷[97,98]。基于此的真核生物起源機(jī)制有很多,幾個(gè)影響深遠(yuǎn)的模型如(圖7,A~D)所示[32,97,99~101]。內(nèi)共生學(xué)說提出者M(jìn)argulis本人認(rèn)為吞噬宿主是一種類似現(xiàn)代spp. (屬于廣古菌門)的嗜酸無壁古細(xì)菌,細(xì)胞核在嵌合細(xì)胞形成后自發(fā)進(jìn)化形成(圖7A)[100,102];著名的互養(yǎng)共棲假說(圖7B)認(rèn)為真核生物來自產(chǎn)甲烷菌與δ-變形桿菌的共生作用[32];Cavalier-Smith本人則認(rèn)為真核生物與古細(xì)菌演化自一個(gè)共同祖先,真核域與古菌域歸屬于新壁總域(Neomura) (圖7D)[98]。

    但是,針對Cavalier-Smith提出的源真核生物概念,21世紀(jì)以來的研究表明,此前發(fā)現(xiàn)的所有不含線粒體的真核生物均含有產(chǎn)氫體(hydrogenosome)或mitosome等細(xì)胞器,與線粒體同源,是上文所說的MROs[1,30,105]。早期支持源真核生物概念的研究也被認(rèn)為受到分析偏差影響,導(dǎo)致源真核生物概念不再流行[30]。

    此后,另一派早期線粒體假說獲得了許多擁護(hù)者。幾個(gè)經(jīng)典的Mito-early模型如(圖7,E~G)所示[31,103,104],相似之處在于都是一對古細(xì)菌和真細(xì)菌首先形成共生關(guān)系,真細(xì)菌是古細(xì)菌代謝產(chǎn)物的受體,隨后內(nèi)共生,形成細(xì)胞核等其他結(jié)構(gòu),其中包含著名的氫假說(圖7E)。

    早期線粒體假說最有力的論據(jù)在于現(xiàn)存生物中尚未發(fā)現(xiàn)不含線粒體或其同源物的真核細(xì)胞過渡形態(tài),這似乎暗示線粒體的內(nèi)共生應(yīng)發(fā)生在一個(gè)較早的時(shí)期,很可能早于細(xì)胞核。Von Dohlen等[106]在非吞噬作用的原核細(xì)胞體內(nèi)發(fā)現(xiàn)存在其他原核細(xì)胞的現(xiàn)象也表明,內(nèi)共生或許不需要真核生物祖先擁有復(fù)雜的吞噬營養(yǎng)系統(tǒng)。

    從線粒體與真核細(xì)胞其他結(jié)構(gòu)的關(guān)系角度也支持早期線粒體假說。Martin等[4,107~110]提出了對線粒體、細(xì)胞核和真核細(xì)胞分裂形式出現(xiàn)順序的推測。在線粒體始祖進(jìn)入宿主胞質(zhì)形成內(nèi)共生體后,兩者間的基因重組和啟動(dòng)子的相互關(guān)聯(lián)導(dǎo)致內(nèi)共生體(原始線粒體)基因在宿主細(xì)胞質(zhì)內(nèi)表達(dá)[107]。這種嵌合的原核基因組模糊了宿主與內(nèi)共生體間的區(qū)別,彼此出現(xiàn)基因片段交換,如從線粒體基因組向宿主基因組轉(zhuǎn)移II型內(nèi)含子等。然而,原核生物的共轉(zhuǎn)錄翻譯機(jī)制無法解決內(nèi)含子剪接速度比蛋白翻譯速度慢從而產(chǎn)生缺陷肽鏈的矛盾,細(xì)胞核的出現(xiàn)則為轉(zhuǎn)錄和翻譯提供了物理屏障[4]。隨后染色體被局限在細(xì)胞核內(nèi),失去與細(xì)胞膜的接觸,在細(xì)胞分裂時(shí)無法自動(dòng)分離,使細(xì)胞體積較原核生物顯著增加[108]。同時(shí),胞質(zhì)內(nèi)由于富含線粒體產(chǎn)生的ATP,蛋白表達(dá)水平顯著升高。為解決染色體分離問題,可能出現(xiàn)了依賴于微管的染色體分離方式,最終形成減數(shù)分裂、真核細(xì)胞周期和有絲分裂過程[109,110]。此外,線粒體早期內(nèi)共生釋放外膜囊泡(outer membrane vesicles, OMVs)的過程也有利于解釋真核生物內(nèi)膜系統(tǒng)的快速形成[111],詳細(xì)機(jī)制可參閱McBride[112]的綜述。

    圖7 有關(guān)線粒體起源的模型

    A~D為晚期線粒體模型,E~G為早期線粒體模型。A:Margulis認(rèn)為古細(xì)菌宿主熱原體菌()首先與螺旋菌(Spirochaete)內(nèi)共生,逐漸演化出細(xì)胞核,后形成線粒體[100]。B:互養(yǎng)共棲假說,產(chǎn)甲烷古細(xì)菌與δ-變形菌在代謝上存在氫的交換。一個(gè)產(chǎn)甲烷菌被多個(gè)δ-變形菌包圍后,δ-變形菌間的相互作用促進(jìn)胞質(zhì)融合,形成內(nèi)質(zhì)網(wǎng)、細(xì)胞核等結(jié)構(gòu),δ-變形菌基因組隨著基因的轉(zhuǎn)移逐漸消失,線粒體的內(nèi)共生最后發(fā)生[32]。C:內(nèi)生核假說,泉古菌與G-真細(xì)菌互相作用,泉古菌形成細(xì)胞核,真細(xì)菌形成細(xì)胞質(zhì),再與α-變形菌內(nèi)共生形成線粒體[101]。D:Cavalier-Smith認(rèn)為真核生物與古細(xì)菌演化自一個(gè)共同祖先,同屬新壁總域[97,99]。E:氫假說,宿主是嚴(yán)格依賴氫且嚴(yán)格自養(yǎng)的厭氧古細(xì)菌,共生體能進(jìn)行有氧呼吸,但通過厭氧異養(yǎng)代謝為宿主提供氫氣[31]。F:古細(xì)菌宿主在形成真核細(xì)胞的其他復(fù)雜結(jié)構(gòu)前首先與耗氧的α-變形菌內(nèi)共生形成線粒體[103]。G:基于硫元素的代謝循環(huán),宿主是產(chǎn)H2S的古細(xì)菌,線粒體祖先是耗H2S的α-變形菌[104]。根據(jù)文獻(xiàn)[30]修改繪制。

    近10年來,有關(guān)Mito-late和Mito-early模型的問題仍有爭論。延續(xù)Mito-early假說風(fēng)潮,主流觀念認(rèn)為線粒體祖先先于宿主吞噬作用出現(xiàn),直接侵入古細(xì)菌細(xì)胞內(nèi)發(fā)生內(nèi)共生,導(dǎo)致真核細(xì)胞的產(chǎn)生[53,113]。然而,發(fā)現(xiàn)Asgard古細(xì)菌Lokiarchaeota后,基因組重建出大量細(xì)胞骨架相關(guān)的ESPs,開啟了具有吞噬作用的原核生物譜系從古細(xì)菌域內(nèi)部進(jìn)化的可能性,支持Mito-late模型[11,114]。但后續(xù)基于宏基因組重建的研究顯示,Asgard古細(xì)菌不具備吞噬營養(yǎng)功能,吞噬作用的形成涉及來自古細(xì)菌、真細(xì)菌以及真核生物新基因的組合[115],宿主是否具有吞噬功能依然未知。不過,限于主流觀念大多是依據(jù)對當(dāng)今真核生物組織和功能方面的觀察所作的進(jìn)化推斷,始終缺乏確切的實(shí)驗(yàn)分析證據(jù),Hampl等[116]便認(rèn)為沒有充分的理由肯定線粒體產(chǎn)生于細(xì)胞吞噬作用之前。

    總的來說,目前Mito-early假說更為流行,但也面臨著其他一些無法解釋的現(xiàn)象:真核生物的嵌合性基因組中,大多數(shù)來自細(xì)菌的基因無法回溯到線粒體的祖先α-變形菌。2016年,Pittis和Ga-baldón[2,117]首次通過比較真核生物與距離最近的原核生物蛋白間的親緣距離,發(fā)現(xiàn)α-變形菌與真細(xì)菌來源蛋白之間的進(jìn)化距離短于真細(xì)菌與古細(xì)菌來源蛋白之間的進(jìn)化距離,推測線粒體內(nèi)共生前,真核細(xì)胞祖先可能已與其他真細(xì)菌內(nèi)共生形成嵌合基因組,胞內(nèi)有內(nèi)質(zhì)網(wǎng)和高爾基體,但還不是進(jìn)化成熟的真核細(xì)胞。Ettema[118]將其稱為介于Mito-early和Mito-late間的中期線粒體(Mito-intermediate)假說。但此結(jié)果引發(fā)許多質(zhì)疑。William和Nelson-Sathi[119]認(rèn)為作者采用了錯(cuò)誤的統(tǒng)計(jì)和分析方法;Esposti也認(rèn)為Pittis和Gabaldón的方法早已被證明失準(zhǔn)[53],低估了α-變形菌祖先對真核生物譜系的影響[120]。此外,目前認(rèn)為與真核生物祖先關(guān)系最密切的Asgard古細(xì)菌本身也存在大量其他真細(xì)菌基因[17]。又鑒于位于現(xiàn)存α-變形菌綱深枝的磁球菌()具有一個(gè)α-變形菌成分僅占1/3的嵌合基因組,Esposti[120]推測久遠(yuǎn)的線粒體祖先本身很可能就是一個(gè)嵌合基因組,而非僅有α-變形菌的血統(tǒng)[121],所以Pittis與Gabaldón的最終結(jié)論顯得有些草率。

    至今,線粒體內(nèi)共生的早期與晚期之爭依舊無法定論,不過可以確定的是,線粒體的內(nèi)共生發(fā)生在TACK和Asgard古細(xì)菌基因組中所編碼的ESPs基因出現(xiàn)之后。若真如Esposti所說,線粒體祖先是一個(gè)嵌合基因組的話,這段長達(dá)30余年的爭論便可以畫上句號。

    2.2.2 線粒體祖先的確定

    探索線粒體始祖細(xì)菌的性質(zhì)有助于人們了解具體的內(nèi)共生過程。長久以來,線粒體祖先與α-變形菌綱相關(guān)的結(jié)論為學(xué)術(shù)界廣泛接受[93],但與其中哪一譜系更加接近的問題始終沒有解決,給重建線粒體與現(xiàn)存α-變形菌譜系的最近祖先(last common ancestor)帶來困難。

    大量通過單基因或全基因組基因串聯(lián)建樹的結(jié)果常常將線粒體與立克次氏體目(Rickettsiales)相關(guān)聯(lián)[122,123]。而當(dāng)引入Pelagibacterales目(SAR11)基因組數(shù)據(jù)作系統(tǒng)發(fā)育分析時(shí),則發(fā)現(xiàn)Pelagibacterales目是線粒體-立克次氏體目一枝的姐妹群[124],或直接是線粒體的姐妹群[125]。后來的分析表明,這些結(jié)果是受到線粒體和Pelagibacterales目之間趨同相似性影響帶來的偏差[126,127]。隨后,依賴于核基因編輯和線粒體基因編輯的線粒體蛋白數(shù)據(jù)庫的分析結(jié)果支持線粒體發(fā)源自立克次氏體目內(nèi)[128,129]。

    但是,上述將線粒體與立克次氏體目關(guān)聯(lián)的研究也存在缺陷。第一,線粒體、立克次氏體目以及上文提到的Pelagibacterales目(SAR11)與其他α-變形菌綱的物種相比進(jìn)化更加迅速且已經(jīng)失去了許多基因,在基因組層面A+T占比頗高,很容易招致LBA偏差[89]。對此,一些研究的應(yīng)對措施是選用較早由線粒體基因組轉(zhuǎn)移到核基因組的基因進(jìn)行分析[128],但不可忽視的一點(diǎn)是前文所提到的核基因組中存在的堿基替換、基因丟失和HGT問題。第二,推測的立克次氏體目在地球上的存在時(shí)間不夠長,目內(nèi)處于主干地位的屬僅有大約5億年歷史[120]。第三,立克次氏體目的代謝方式也與線粒體差別很大,普遍不存在現(xiàn)代真核生物的厭氧性狀,不編碼調(diào)節(jié)線粒體嵴的MICOS蛋白和功能等同于線粒體ATP合酶抑制亞基的相關(guān)結(jié)構(gòu)[120]。

    另一種通過單基因數(shù)據(jù)集和多基因串聯(lián)數(shù)據(jù)集推斷網(wǎng)絡(luò)的方法將線粒體與其他α-變形菌綱分支相連,如紅螺菌目(Rhodospirillales)的紅螺菌屬()和根瘤菌目(Rhizobiales)的蒼白桿菌屬()[89]。從代謝相似性來看,有分析認(rèn)為線粒體可能更加接近于紅螺菌屬[120]。還有研究認(rèn)為線粒體來自α-變形菌綱內(nèi)的一個(gè)獨(dú)立深枝[127]。

    2018年,Martijn等[130]通過宏基因組重建得到12個(gè)α-變形菌綱譜系(包含幾個(gè)新發(fā)現(xiàn)的譜系),將線粒體放在了這些譜系的進(jìn)化分支之外,是所有測序的α-變形菌綱序列的姐妹群,與此前的眾多結(jié)果存在矛盾。然而,2019年Fan等[131]的結(jié)果則又顯著支持線粒體來源于α-變形菌綱內(nèi)的傳統(tǒng)觀點(diǎn),并指出Martijn等在規(guī)避異源基因的處理時(shí)可能存在模型過度擬合的缺陷。

    形成這一僵局的原因來自多方面。首先,許多現(xiàn)存的α-變形菌綱生物受到與其他細(xì)菌之間的雙向水平基因轉(zhuǎn)移(HGT)影響[132],線粒體始祖本身甚至可能就有一個(gè)嵌合體基因組[120];其次,線粒體基因在進(jìn)化上比一般自由生活的細(xì)菌更加迅速[128],且線粒體基因組向核基因的轉(zhuǎn)移與時(shí)間呈現(xiàn)非線性關(guān)系[133],為模型的建立帶來困難;最后,α-變形菌綱作為一個(gè)成員變異性極大的綱,現(xiàn)存成員間的進(jìn)化關(guān)系仍在調(diào)整,例如Mu?oz-Gómez等[134]的研究已將Holosporales目從“廣義的立克次氏體目”(Rickett-siales sensu lato)中分離出來。環(huán)境微生物基因組的測定也在不斷發(fā)現(xiàn)新的α-變形菌綱物種。

    綜上,線粒體始祖的判定仍是探索真核生物起源道路上亟需解決的問題。

    2.2.3 古細(xì)菌宿主的確定

    長期以來的系統(tǒng)發(fā)育基因組學(xué)研究已經(jīng)基本確定線粒體內(nèi)共生的宿主與古細(xì)菌相關(guān)[84]。宏基因組技術(shù)出現(xiàn)之前,古細(xì)菌分為廣古菌門和泉古菌門兩門。宿主身份有頗多說法。泉古菌假說支持宿主屬于泉古菌門[28,29];氫假說(圖7E)[31]與互養(yǎng)共棲假說(圖7B)[32]支持宿主屬于廣古菌門中的產(chǎn)甲烷菌(methanogens);Searcy等[135]與Margulis等[100]支持宿主屬于廣古菌門中的熱原體菌() (圖7A)。這些模型雖彼此對立,但共同點(diǎn)是線粒體內(nèi)共生僅涉及一個(gè)古細(xì)菌宿主和一個(gè)真細(xì)菌內(nèi)共生體。然而,與此不同的其他一些模型則認(rèn)為宿主是兩種原核生物的嵌合體,這些觀點(diǎn)被統(tǒng)一稱作三元假說(ternary hypothesis)。其中較為流行的是內(nèi)生核假說(endokaryotic hypothesis) (圖7C),推測宿主的細(xì)胞核起源于一種古細(xì)菌,而細(xì)胞質(zhì)起源于一種真細(xì)菌,以此解釋真核生物基因組中無法追溯到α-變形菌綱的細(xì)菌來源基因[101,136]。然而,后期的系統(tǒng)發(fā)育檢驗(yàn)并沒有發(fā)現(xiàn)這樣一個(gè)細(xì)菌譜系,這表明內(nèi)生核假說缺乏證據(jù)支持[53]。面對眾多模型長期不可調(diào)和的僵局,2012年,Godd[137]提出了另一種三元假說。在融合納古菌屬() (當(dāng)時(shí)未獨(dú)立為門,屬于廣古菌門)和燃球菌屬() (屬于泉古菌門)兩種古細(xì)菌的基因組后,Godd發(fā)現(xiàn)可以解釋當(dāng)時(shí)的大部分?jǐn)?shù)據(jù),提出宿主可能是兩種古細(xì)菌的嵌合體(chimera)。這一設(shè)想具有很強(qiáng)的創(chuàng)新性,但仍然需要大量的系統(tǒng)發(fā)育分析和細(xì)胞學(xué)研究檢驗(yàn)才能令人信服。

    在結(jié)合宏基因組和單細(xì)胞基因組分析技術(shù)發(fā)現(xiàn)TACK超門和Asgard超門以后,內(nèi)共生宿主的身份逐漸明朗,也基本消除了第三個(gè)原始真核生物譜系存在的可能性[11,17]。Asgard古細(xì)菌中大量ESPs的出現(xiàn)同時(shí)否定了其他一些有關(guān)病毒影響真核生物起源的理論。這些理論認(rèn)為真核生物的大量基因,甚至細(xì)胞核本身都來自核質(zhì)大DNA病毒(nucleo-cytoplasmic large DNA viruses, NCLDVs)[138,139],其基于的證據(jù)主要是在NCLDVs中發(fā)現(xiàn)了許多與真核生物相關(guān)的基因。而如今Asgard古細(xì)菌中大量ESPs的發(fā)現(xiàn)以及隨后糾正的系統(tǒng)發(fā)育分析則表明,NCLDVs中的真核樣基因?qū)嶋H來自于真核生物病毒的HGT。病毒的高進(jìn)化速率為此前的系統(tǒng)發(fā)育重建帶來偏差,導(dǎo)致了NCLDVs是真核生物基因甚至是細(xì)胞核源頭的判斷。事實(shí)上,NCLDVs可能是Asgard古細(xì)菌或相關(guān)譜系的古老成員[90]。

    如今,對Asgard古細(xì)菌的分析研究已經(jīng)成為了解真核生物祖先的重要手段。Asgard超門呈現(xiàn)出豐富的代謝方式。宏基因組重建分析推測Asgard超門普遍存在氮和硫代謝循環(huán)[61],不能進(jìn)行吞噬作用[115],Lokiarchaeota存在依賴氫的代謝途徑[140],Thorar-chaeota存在依賴無機(jī)碳和有機(jī)碳的混合代謝以及產(chǎn)乙酸途徑[62,141],Lokiarchaeota和Thorarchaeota存在鹵代有機(jī)物代謝途徑[142],Heimdallarchaeota存在光能營養(yǎng)途徑[143,144],Helarchaeota存在厭氧碳?xì)浠衔锎x途徑[75],部分Lokiarchaeota編碼鎳鐵氫化酶[142]?;谶@些代謝研究提出的真核生物起源新模型也不斷涌現(xiàn),如電子或氫從有機(jī)異養(yǎng)古細(xì)菌宿主流向細(xì)菌共生體的“逆流模型”(reverse flow model) (圖8A)[145]。此外,對水底生物墊內(nèi)Asgard古細(xì)菌與周圍其他細(xì)菌的基因組綜合研究也有助于闡釋線粒體祖先與宿主內(nèi)共生前的共生關(guān)系,例如Sagha?等[146]便利用16S/18S rRNA宏條形碼和多元統(tǒng)計(jì)分析技術(shù),在鹽沼池水中發(fā)現(xiàn)了Lokiarchaeota與此前知之甚少的細(xì)菌門類TA06間的強(qiáng)正相關(guān)關(guān)系,可能存在共生作用。

    由于大量Asgard古細(xì)菌的研究結(jié)論都只是基于宏基因組分析的結(jié)果,細(xì)胞生物學(xué)層面的研究驗(yàn)證便顯得尤為重要。然而,Asgard古細(xì)菌的生長環(huán)境大多處于海底沉積物中,且可能依賴其他共生細(xì)菌生存,這為Asgard古細(xì)菌的分離培養(yǎng)和蛋白的提純帶來巨大困難,相關(guān)工作開展緩慢。Ak?l和Robinson[147]取得了標(biāo)志性成果。他們利用質(zhì)粒載體表達(dá)Asgard古細(xì)菌的前纖維蛋白編碼基因,發(fā)現(xiàn)表達(dá)產(chǎn)物能夠成功與哺乳動(dòng)物的肌動(dòng)蛋白結(jié)合,即使相距2、30億年的進(jìn)化距離依舊與哺乳動(dòng)物的相關(guān)蛋白呈現(xiàn)結(jié)構(gòu)和功能的相似性,這為Asgard超門與真核生物的相似性提供了一個(gè)細(xì)胞學(xué)證據(jù)。

    圖8 真核生物起源的逆流模型與E3模型的對比

    A:逆流模型;B:E3模型。兩個(gè)模型的共同點(diǎn)是異養(yǎng)古細(xì)菌宿主與α-變形菌內(nèi)共生形成線粒體;不同點(diǎn)在于E3模型比逆流模型多引入了一個(gè)共生菌——硫酸鹽還原細(xì)菌。

    近年來,對Asgard古細(xì)菌的研究終于迎來里程碑式的成果。日本科學(xué)家Imachi等[13]花費(fèi)12年時(shí)間利用連續(xù)流動(dòng)生物反應(yīng)器,成功從深海沉積物中分離,并首次在實(shí)驗(yàn)室環(huán)境下提純培養(yǎng)出一個(gè)Asgard古細(xì)菌菌株MK-D1,屬于Lokiarchaeota。該研究的原始目的是尋找甲烷降解菌[148]。分離到的Lokiarchaeota門古細(xì)菌生長緩慢,滯后期長達(dá)3~6個(gè)月,二倍擴(kuò)增需14~25天。進(jìn)一步觀察顯示,這是一種小型厭氧球菌,細(xì)胞可以形成凸起結(jié)構(gòu)、氣泡和胞外小泡,沒有明顯的真核樣細(xì)胞器,具備厭氧代謝氨基酸的能力,依賴共生菌一起生存。Imachi等據(jù)此提出了新的真核生物起源模型——糾纏–吞噬–內(nèi)生(entangle–engulf–endogenize, E3)模型:厭氧古細(xì)菌宿主與一種硫酸鹽還原細(xì)菌共生,共生細(xì)菌是宿主氨基酸代謝產(chǎn)生的脂肪酸和氫氣的受體,隨后宿主與另一好氧細(xì)菌內(nèi)共生形成線粒體,共涉及3種原核生物。這與此前僅涉及兩個(gè)原核生物的“逆流模型”[145]存在不同(圖8)。逆流模型認(rèn)為古細(xì)菌代謝產(chǎn)生的氫氣直接被線粒體祖先消耗,無需引入假設(shè)的另一個(gè)共生伴侶(圖8A),而E3模型中的線粒體起源是處于共生關(guān)系下的厭氧宿主古細(xì)菌對環(huán)境中氧氣濃度增加的適應(yīng)行為(圖8B)。

    在未來,有關(guān)古細(xì)菌宿主的研究無疑將借鑒Imachi等的經(jīng)驗(yàn),從環(huán)境中發(fā)現(xiàn)并在實(shí)驗(yàn)室條件下培養(yǎng)研究與真核生物進(jìn)化距離更近的微生物。

    3 結(jié)語與展望

    近幾十年來對真核生物起源的研究積累了大量數(shù)據(jù)、模型和觀點(diǎn),這是一個(gè)涉及進(jìn)化生物學(xué)、微生物學(xué)、生物化學(xué)、分子生物學(xué)、細(xì)胞生物學(xué)和古生物學(xué)等多領(lǐng)域的宏大問題。新技術(shù)的發(fā)明(如16S/18S rRNA測序、基因組測序、宏基因組測序和單細(xì)胞全基因組測序)、新物種的發(fā)現(xiàn)(如古細(xì)菌域、TACK超門和Asgard超門)以及分析方法的改進(jìn)不斷推動(dòng)我們對于真核生物起源的認(rèn)知。

    然而爭論依然存在,主流觀點(diǎn)將生命之樹上的真核生物分支放在古細(xì)菌域內(nèi),從Asgard超門內(nèi)發(fā)出或是Asgard超門的姐妹群,對應(yīng)二域系統(tǒng)生命之樹;真核生物的始祖是一個(gè)類似Asgard超門的古細(xì)菌,在與一種α-變形菌綱細(xì)菌內(nèi)共生形成線粒體后促進(jìn)其他細(xì)胞內(nèi)部結(jié)構(gòu)的進(jìn)化。

    未來的研究重點(diǎn),一方面在于改進(jìn)系統(tǒng)發(fā)育分析方法,并延續(xù)環(huán)境微生物的宏基因組重建分析工作。在確定古細(xì)菌宿主的角度上,發(fā)現(xiàn)新的Asgard古細(xì)菌,甚至是其他與真核生物更接近的古細(xì)菌新門類。在確定線粒體祖先的角度上,亟需理清α-變形菌綱的譜系關(guān)系。目前可獲取到基因組數(shù)據(jù)的α-變形菌綱物種大多與醫(yī)學(xué)與農(nóng)業(yè)相關(guān),或者是部分可在實(shí)驗(yàn)室條件下分離培養(yǎng)的菌種[120,130],可能存在一定的取樣偏差,不一定能完全代表現(xiàn)存α-變形菌的多樣性。因此,在海水、地下水和土壤等環(huán)境中發(fā)現(xiàn)未培養(yǎng)微生物,擴(kuò)展α-變形菌綱物種數(shù)的工作可能仍是重中之重。此外,發(fā)展基因樹–物種樹(gene tree-species tree)綜合分析[149]等更加復(fù)雜的計(jì)算方法也有望為線粒體的祖先提供更加精確的定位。另一方面,實(shí)驗(yàn)室中的細(xì)胞生物學(xué)研究亟需跟進(jìn)。在確定古細(xì)菌宿主的角度上,借鑒Imachi等[13]的工作,實(shí)驗(yàn)室條件下直接培養(yǎng)Asgard古細(xì)菌進(jìn)行基因組測序可以避免宏基因組重建時(shí)的誤差,直接觀察細(xì)胞形態(tài)和代謝過程。但這面臨兩個(gè)問題,一是需要改進(jìn)生物反應(yīng)器系統(tǒng),縮短Asgard古細(xì)菌的分離提純時(shí)間;二是需要研究自然條件下Asgard古細(xì)菌的生活環(huán)境和共生細(xì)菌等,摸索其在實(shí)驗(yàn)室中的最適培養(yǎng)條件,加快繁殖速度以供實(shí)驗(yàn)需求。在確定線粒體祖先的角度上,Hu[150]通過冷凍電鏡發(fā)現(xiàn)HeLa細(xì)胞線粒體存在與許多現(xiàn)存α-變形菌綱生物相似的出芽增殖方式,不同于傳統(tǒng)主流觀點(diǎn)所描述的二分裂增殖方式,而這一現(xiàn)象難以在常用的熒光顯微鏡下觀察到。此結(jié)果說明或許可以將新一代顯微成像技術(shù)應(yīng)用于線粒體祖先的研究中,從形態(tài)學(xué)層面比較線粒體與立克次氏體等現(xiàn)存α-變形菌的異同,為系統(tǒng)發(fā)育學(xué)分析提供必要補(bǔ)充。此外,實(shí)驗(yàn)室內(nèi)對真核細(xì)胞進(jìn)化某些步驟的直接模擬實(shí)驗(yàn)也有利于我們還原真核生物始祖的內(nèi)共生效應(yīng)和細(xì)胞器的進(jìn)化歷程。相關(guān)研究工作已在開展之中,例如2018年Mehta等[151]報(bào)道的在酵母胞質(zhì)內(nèi)培養(yǎng)大腸桿菌多達(dá)120天的工作。

    最后,值得提出的一點(diǎn)是,現(xiàn)在采用的Asgard古細(xì)菌等研究對象并不古老,也同真核生物一樣走過了十幾億年的輻射進(jìn)化歷程?,F(xiàn)存代謝方式豐富多樣的Asgard古細(xì)菌不可能完全代表原始的真核生物。未來的目標(biāo)將是通過分析大量現(xiàn)存與即將涌現(xiàn)的數(shù)據(jù),提出一個(gè)適用于所有真核生物特征的真核生物起源模型。

    致謝

    本文在撰稿及寫作過程中得到南京大學(xué)生命科學(xué)學(xué)院楊四海教授的悉心指導(dǎo),特此表示由衷的感謝。

    [1] Blackstone NW. An evolutionary framework for under-standing the origin of eukaryotes., 2016, 5(2): 18.

    [2] Gabaldón T. Relative timing of mitochondrial endosy-mbiosis and the “pre-mitochondrial symbioses” hypo-thesis., 2018, 70(12): 1188–1196.

    [3] Fuerst JA. Intracellular compartmentation in plancto-mycetes., 2005, 59: 299–328.

    [4] Martin W, Koonin EV. Introns and the origin of nucleus-cytosol compartmentalization., 2006, 440(7080): 41–45.

    [5] Lin L, Shi AB. Endocytic recycling pathways and the regulatory mechanisms., 2019, 41(6): 451–468.林瓏, 史岸冰. 細(xì)胞內(nèi)吞循環(huán)運(yùn)輸通路及其分子調(diào)控機(jī)制. 遺傳, 2019, 41(6): 451–468.

    [6] Dacks JB, Field MC, Buick R, Eme L, Gribaldo S, Roger AJ, Brochier-Armanet C, Devos DP. The changing view of eukaryogenesis–fossils, cells, lineages and how they all come together., 2016, 129(20): 3695–3703.

    [7] Shen Y, Buick R, Canfield DE. Isotopic evidence for microbial sulphate reduction in the early Archaean era., 2001, 410(6824): 77–81.

    [8] Ueno Y, Ono S, Rumble D, Maruyama S. Quadruple sulfur isotope analysis of ca. 3.5 Ga Dresser Formation: New evidence for microbial sulfate reduction in the early Archean., 2008, 72(23): 5675–5691.

    [9] Yan YZ, Liu ZL. Significance of eucaryotic organisms in the microfossil flora of Changcheng system., 1993, 10(2): 167–180.

    [10] Parfrey LW, Lahr DJG, Knoll AH, Katz LA. Estimating the timing of early eukaryotic diversification with multigene molecular clocks., 2011, 108(33): 13624–13629.

    [11] Spang A, Saw JH, J?rgensen SL, Zaremba-Niedzwiedzka K, Martijn J, Lind AE, van Eijk R, Schleper C, Guy L, Ettema TJG. Complex archaea that bridge the gap between prokaryotes and eukaryotes., 2015, 521(7551): 173–179.

    [12] Zaremba-Niedzwiedzka K, Caceres EF, Saw JH, B?ckstr?m D, Juzokaite L, Vancaester E, Seitz KW, Anantharaman K, Starnawski P, Kjeldsen KU, Stott MB, Nunoura T, Banfield JF, Schramm A, Baker BJ, Spang A, Ettema TJ. Asgard archaea illuminate the origin of eukaryotic cellular complexity., 2017, 541(7637): 353–358.

    [13] Imachi H, Nobu MK, Nakahara N, Morono Y, Ogawara M, Takaki Y, Takano Y, Uematsu K, Ikuta T, Ito M, Matsui Y, Miyazaki M, Murata K, Saito K, Sakai S, Song CH, Tasumi E, Yamanaka Y, Yamaguchi T, Kamagata Y, Tamaki H, Takai K. Isolation of an archaeon at the prokaryote–eukaryote interface., 2020, 577(7791): 519–525.

    [14] Woese CR, Fox GE. Phylogenetic structure of the prokaryotic domain: the primary kingdoms., 1977, 74(11): 5088–5090.

    [15] Iwabe N, Kuma K, Hasegawa M, Osawa S, Miyata T. Evolutionary relationship of archaebacteria, eubacteria, and eukaryotes inferred from phylogenetic trees of duplicated genes., 1989, 86(23): 9355–9359.

    [16] Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya., 1990, 87(12): 4576–4579.

    [17] Eme L, Spang A, Lombard J, Stairs CW, Ettema TJG. Archaea and the origin of eukaryotes., 2017, 15(12): 711–723.

    [18] Pace NR. Time for a change., 2006, 441(7091): 289–289.

    [19] Huet J, Schnabel R, Sentenac A, Zillig W. Archae-bacteria and eukaryotes possess DNA‐dependent RNA polymerases of a common type., 1983, 2(8): 1291–1294.

    [20] Zillig W, Klenk HP, Palm P, Pühler G, Gropp F, Garrett RA, Leffers H. The phylogenetic relations of DNA- dependent RNA polymerases of archaebacteria, eukaryotes, and eubacteria., 1989, 35(1): 73–80.

    [21] Gribaldo S, Poole AM, Daubin V, Forterre P, Brochier-Armanet C. The origin of eukaryotes and their relationship with the Archaea: are we at a phylogenomic impasse?, 2010, 8(10): 743–752.

    [22] Lind?s AC, Karlsson EA, Lindgren MT, Ettema TJ, Bernander R. A unique cell division machinery in the Archaea., 2008, 105(48): 18942–18946.

    [23] Samson RY, Obita T, Freund SM, Williams RL, Bell SD. A role for the ESCRT system in cell division in archaea., 2008, 322(5908): 1710–1713.

    [24] Li J, Zhang KQ. Genetic diversity of microorganisms., 2012, 34(11): 1399–1408.李娟,張克勤. 微生物的遺傳多樣性. 遺傳, 2012, 34(11): 1399–1408.

    [25] Rivera MC, Jain R, Moore JE, Lake JA. Genomic evidence for two functionally distinct gene classes., 1998, 95(11): 6239–6244.

    [26] Canback B, Andersson SG, Kurland CG. The global phylogeny of glycolytic enzymes., 2002, 99(9): 6097–6102.

    [27] Peretó J, López-García P, Moreira D. Ancestral lipid biosynthesis and early membrane evolution., 2004, 29(9): 469–477.

    [28] Lake JA, Henderson E, Oakes M, Clark MW. Eocytes: a new ribosome structure indicates a kingdom with a close relationship to eukaryotes., 1984, 81(12): 3786–3790.

    [29] Rivera MC, Lake JA. Evidence that eukaryotes and eocyte prokaryotes are immediate relatives., 1992, 257(5066): 74–76.

    [30] Embley TM, Martin W. Eukaryotic evolution, changes and challenges., 2006, 440(7084): 623–630.

    [31] Martin W, Müller M. The hydrogen hypothesis for the first eukaryote., 1998, 392(6671): 37–41.

    [32] Moreira D, López-García P. Symbiosis between methanogenic archaea and δ-proteobacteria as the origin of eukaryotes: the syntrophic hypothesis., 1998, 47(5): 517–530.

    [33] Harris JK, Kelley ST, Spiegelman GB, Pace NR. The genetic core of the universal ancestor., 2003, 13(3): 407–412.

    [34] Ciccarelli FD, Doerks T, Von Mering C, Creevey CJ, Snel B, Bork P. Toward automatic reconstruction of a highly resolved tree of life., 2006, 311(5765): 1283–1287.

    [35] Yutin N, Makarova KS, Mekhedov SL, Wolf YI, Koonin EV. The deep archaeal roots of eukaryotes., 2008, 25(8): 1619–1630.

    [36] Rivera MC, Lake JA. The ring of life provides evidence for a genome fusion origin of eukaryotes., 2004, 431(7005): 152–155.

    [37] Pisani D, Cotton JA, McInerney JO. Supertrees disentangle the chimerical origin of eukaryotic genomes., 2007, 24(8): 1752–1760.

    [38] Cox CJ, Foster PG, Hirt RP, Harris SR, Embley TM. The archaebacterial origin of eukaryotes., 2008, 105(51): 20356–20361.

    [39] Foster PG, Cox CJ, Embley TM. The primary divisions of life: a phylogenomic approach employing composition- heterogeneous methods., 2009, 364(1527): 2197–2207.

    [40] Boussau B, Daubin V. Genomes as documents of evolutionary history., 2010, 25(4): 224–232.

    [41] Som A. Causes, consequences and solutions of phylogenetic incongruence., 2015, 16(3): 536–548.

    [42] Lerat E, Daubin V, Ochman H, Moran NA. Evolutionary origins of genomic repertoires in bacteria., 2005, 3(5): e130.

    [43] Soucy SM, Huang J, Gogarten JP. Horizontal gene transfer: building the web of life., 2015, 16(8): 472–482.

    [44] Wagner A, Whitaker RJ, Krause DJ, Heilers JH, Van Wolferen M, Van Der Does C, Albers SV. Mechanisms of gene flow in archaea., 2017, 15(8): 492–501.

    [45] Akanni WA, Siu-Ting K, Creevey CJ, McInerney JO, Wilkinson M, Foster PG, Pisani D. Horizontal gene flow from Eubacteria to Archaebacteria and what it means for our understanding of eukaryogenesis., 2015, 370(1678): 20140337.

    [46] Cotton JA, McInerney JO. Eukaryotic genes of archaebacterial origin are more important than the more numerous eubacterial genes, irrespective of function., 2010, 107(40): 17252–17255.

    [47] Cohen O, Gophna U, Pupko T. The complexity hypothesis revisited: connectivity rather than function constitutes a barrier to horizontal gene transfer., 2011, 28(4): 1481–1489.

    [48] Williams TA, Foster PG, Cox CJ, Embley TM. An archaeal origin of eukaryotes supports only two primary domains of life., 2013, 504(7479): 231–236.

    [49] Williams TA, Cox CJ, Foster PG, Sz?ll?si GJ, Embley TM. Phylogenomics provides robust support for a two-domains tree of life., 2020, 4(1): 138–147.

    [50] Tourasse NJ, Gouy M. Accounting for evolutionary rate variation among sequence sites consistently changes universal phylogenies deduced from rRNA and protein-coding genes., 1999, 13(1): 159–168.

    [51] Guy L, Ettema TJ. The archaeal ‘TACK’ superphylum and the origin of eukaryotes., 2011, 19(12): 580–587.

    [52] Williams TA, Foster PG, Nye TM, Cox CJ, Embley TM. A congruent phylogenomic signal places eukaryotes within the Archaea., 2012, 279(1749): 4870–4879.

    [53] Rochette NC, Brochier-Armanet C, Gouy M. Phylogenomic test of the hypotheses for the evolutionary origin of eukaryotes., 2014, 31(4): 832–845.

    [54] Venter JC, Remington K, Heidelberg JF, Halpern AL, Rusch D, Eisen JA, Wu D, Paulsen I, Nelson KE, Nelson W, Fouts DE, Levy S, Knap AH, Lomas MW, Nealson K, White O, Peterson J, Hoffman J, Parsons R, Baden-Tillson H, Pfannkoch C, Rogers YH, Smith HO. Environmental genome shotgun sequencing of the Sargasso Sea., 2004, 304(5667): 66–74.

    [55] Brochier-Armanet C, Boussau B, Gribaldo S, Forterre P. Mesophilic Crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota., 2008, 6(3): 245–252.

    [56] Elkins JG, Podar M, Graham DE, Makarova KS, Wolf Y, Randau L, Hedlund BP, Brochier-Armanet C, Kunin V, Anderson I, apidus A, Goltsman E, Barry K, Koonin EV, Hugenholtz P, Kyrpides N, Wanner G, Richardson P, Keller M, Stettera K. A korarchaeal genome reveals insights into the evolution of the Archaea., 2008, 105(23): 8102–8107.

    [57] Nunoura T, Takaki Y, Kakuta J, Nishi S, Sugahara J, Kazama H, Chee GJ, Hattori M, Kanai A, Atomi H, Takai K, Takami H. Insights into the evolution of Archaea and eukaryotic protein modifier systems revealed by the genome of a novel archaeal group., 2011, 39(8): 3204–3223.

    [58] Guy L, Saw JH, Ettema TJ. The archaeal legacy of eukaryotes: a phylogenomic perspective., 2014, 6(10): a016022.

    [59] Petitjean C, Deschamps P, López-García P, Moreira D. Rooting the domain archaea by phylogenomic analysis supports the foundation of the new kingdom Proteoarchaeota., 2015, 7(1): 191– 204.

    [60] Raymann K, Brochier-Armanet C, Gribaldo S. The two-domain tree of life is linked to a new root for the Archaea., 2015, 112(21): 6670–6675.

    [61] MacLeod F, Kindler GS, Wong HL, Chen R, Burns BP. Asgard archaea: diversity, function, and evolutionary implications in a range of microbiomes., 2019, 5(1): 48–61.

    [62] Seitz KW, Lazar CS, Hinrichs KU, Teske AP, Baker BJ. Genomic reconstruction of a novel, deeply branched sediment archaeal phylum with pathways for aceto-genesis and sulfur reduction., 2016, 10(7): 1696–1705.

    [63] López-García P, Moreira D. Eukaryogenesis, a syntrophy affair., 2019, 4(7): 1068–1070.

    [64] Da Cunha V, Gaia M, Gadelle D, Nasir A, Forterre P. Lokiarchaea are close relatives of Euryarchaeota, not bridging the gap between prokaryotes and eukaryotes., 2017, 13(6): e1006810.

    [65] Spang A, Eme L, Saw JH, Caceres EF, Zaremba- Niedzwiedzka K, Lombard J, Guy L, Ettema TJG. Asgard archaea are the closest prokaryotic relatives of eukaryotes., 2018, 14(3): e1007080.

    [66] Da Cunha V, Gaia M, Nasir A, Forterre P. Asgard archaea do not close the debate about the universal tree of life topology., 2018, 14(3): e1007215.

    [67] Narrowe AB, Spang A, Stairs CW, Caceres EF, Baker BJ, Miller CS, Ettema TJG. Complex evolutionary history of translation Elongation Factor 2 and diphthamide biosynthesis in Archaea and parabasalids., 2018, 10(9): 2380–2393.

    [68] Harish A, Tunlid A, Kurland CG. Rooted phylogeny of the three superkingdoms., 2013, 95(8): 1593–1604.

    [69] Harish A, Kurland CG. Akaryotes and Eukaryotes are independent descendants of a universal common ancestor., 2017, 138: 168–183.

    [70] Harish A. What is an archaeon and are the Archaea really unique?, 2018, 6: e5770.

    [71] Caetano-Anollés G, Caetano-Anollés D. An evolu-tionarily structured universe of protein architecture., 2003, 13(7): 1563–1571.

    [72] Yang S, Doolittle RF, Bourne PE. Phylogeny determined by protein domain content., 2005, 102(2): 373–378.

    [73] Harish A, Morrison DA. The deep(er) roots of Eukaryotes and Akaryotes., 2020, 9(112): 112.

    [74] Fournier GP, Poole AM. A briefly argued case that Asgard archaea are part of the eukaryote tree., 2018, 9: 1896.

    [75] Seitz KW, Dombrowski N, Eme L, Spang A, Lombard J, Sieber JR, Teske AP, Ettema TJ, Baker BJ. Asgard archaea capable of anaerobic hydrocarbon cycling., 2019, 10(1): 1–11.

    [76] Cai MW, Liu Y, Zhou ZC, Yang YC, Pan J, Gu JD, Li M. Asgard archaea are diverse, ubiquitous, and transcrip-tionally active microbes., 2018: 374165.

    [77] Cai MW, Liu Y, Yin XR, Zhou ZC, Friedrich MW, Richter-Heitmann T, Nimzyk R, Kulkarni A, Wang XW, Li WJ, Pan J, Yang YC, Gu JD, Li M. Highly diverse Asgard archaea participate in organic matter degradation in coastal sediments., 2019: 858530.

    [78] Zhang RY, Zou B, Yan YW, Jeon CO, Li M, Cai M, Quan ZX. Design of targeted primers based on 16S rRNA sequences in meta-transcriptomic datasets and identification of a novel taxonomic group in the Asgard archaea., 2020, 20(1): 25.

    [79] Caceres EF. Genomic and evolutionary exploration of Asgard archaea., 2019, 60–61.

    [80] Koonin EV. The origin and early evolution of eukaryotes in the light of phylogenomics., 2010, 11(5): 209.

    [81] Koumandou VL, Wickstead B, Ginger ML, Van Der Giezen M, Dacks JB, Field MC. Molecular paleontology and complexity in the last eukaryotic common ancestor., 2013, 48(4): 373–396.

    [82] Rout MP, Field MC. The evolution of organellar coat complexes and organization of the eukaryotic cell., 2017, 86: 637–657.

    [83] Gould SB. Membranes and evolution., 2018, 28(8): R381–R385.

    [84] López-García P, Moreira D. Open questions on the origin of eukaryotes., 2015, 30(11): 697–708.

    [85] Gray MW. Mosaic nature of the mitochondrial proteome: Implications for the origin and evolution of mitochondria., 2015, 112(33): 10133–10138.

    [86] Santos HJ, Makiuchi T, Nozaki T. Reinventing an organelle: the reduced mitochondrion in parasitic protists., 2018, 34(12): 1038–1055.

    [87] Tachezy J. Hydrogenosomes and mitosomes: mitochondria of anaerobic eukaryotes., 2019, 1-20.

    [88] Karnkowska A, Vacek V, Zubá?ová Z, Treitli SC, Petr?elková R, Eme L, Novák L, ?ársky V, Barlow LD, Herman EK, Soukal P, Hroudová M, Dole?al P, Stairs CW, Roger AJ, Eliá? M, Dacks JB, Vl?ek ?, Hampl V. A eukaryote without a mitochondrial organelle., 2016, 26(10): 1274–1284.

    [89] Roger AJ, Mu?oz-Gómez SA, Kamikawa R. The origin and diversification of mitochondria., 2017, 27(21): R1177–R1192.

    [90] López-García P, Eme L, Moreira D. Symbiosis in eukaryotic evolution., 2017, 434: 20–33.

    [91] Sagan L. On the origin of mitosing cells., 1967, 14(3): 255–274.

    [92] Bonen L, Cunningham RS, Gray MW, Doolittle WF. Wheat embryo mitochondrial 18S ribosomal RNA: evidence for its prokaryotic nature., 1977, 4(3): 663–671.

    [93] Schwartz RM, Dayhoff MO. Origins of prokaryotes, eukaryotes, mitochondria, and chloroplasts., 1978, 199(4327): 395–403.

    [94] Yang D, Oyaizu Y, Oyaizu H, Olsen GJ, Woese CR. Mitochondrial origins., 1985, 82(13): 4443–4447.

    [95] Cavalier-Smith T. Eukaryotes with no mitochondria., 1987, 326(6111): 332–333.

    [96] Cavalier-Smith T. Archaebacteria and archezoa., 1989, 339(6220): 100–101.

    [97] Cavalier-Smith T. The phagotrophic origin of eukaryotes and phylogenetic classification of Protozoa., 2002, 52(2): 297–354.

    [98] Cavalier-Smith T. Origin of the cell nucleus, mitosis and sex: roles of intracellular coevolution., 2010, 5(1): 7.

    [99] Cavalier-Smith T. Only six kingdoms of life., 2004, 271(1545): 1251–1262.

    [100] Margulis L, Dolan MF, Whiteside JH. “Imperfections and oddities” in the origin of the nucleus., 2005, 31(2): 175–191.

    [101] Sapp J. Microbial Phylogeny and Evolution: Concepts and Controversies. Oxford: Oxford University Press on Demand, 2005, 184–206.

    [102] Margulis L, Dolan MF, Guerrero R. The chimeric eukaryote: origin of the nucleus from the karyomastigont in amitochondriate protists., 2000, 97(13): 6954–6959.

    [103] Vellai T, Takács K, Vida G. A new aspect to the origin and evolution of eukaryotes., 1998, 46(5): 499–507.

    [104] Searcy DG. Origins of mitochondria and chloroplasts from sulfur-based symbioses. In: Hartman H, Matsuno K eds. The Origin and Evolution of the Cell. Singapore: World Scientific, 1992, 47–78.

    [105] van der Giezen M, Tovar J, Clark CG. Mitochondrion- derived organelles in protists and fungi., 2005, 244: 175–225.

    [106] Von Dohlen CD, Kohler S, Alsop ST, McManus WR. Mealybug β-proteobacterial endosymbionts contain γ-proteobacterial symbionts., 2001, 412(6845): 433–436.

    [107] Henze K, Martin W. How do mitochondrial genes get into the nucleus?, 2001, 17(7): 383–387.

    [108] Martin WF, Garg S, Zimorski V. Endosymbiotic theories for eukaryote origin., 2015, 370(1678): 20140330.

    [109] Lane N, Martin W. The energetics of genome complexity., 2010, 467(7318): 929–934.

    [110] Garg SG, Martin WF. Mitochondria, the cell cycle, and the origin of sex via a syncytial eukaryote common ancestor., 2016, 8(6): 1950–1970.

    [111] Gould SB, Garg SG, Martin WF. Bacterial vesicle secretion and the evolutionary origin of the eukaryotic endomembrane system., 2016, 24(7): 525–534.

    [112] McBride HM. Mitochondria and endomembrane origins., 2018, 28(8): R367–R372.

    [113] Yutin N, Wolf MY, Wolf YI, Koonin EV. The origins of phagocytosis and eukaryogenesis., 2009, 4(1): 9.

    [114] Embley TM, Williams TA. Steps on the road to eukaryotes., 2015, 521(7551): 169–170.

    [115] Burns JA, Pittis AA, Kim E. Gene-based predictive models of trophic modes suggest Asgard archaea are not phagocytotic., 2018, 2(4): 697–704.

    [116] Hampl V, ?epi?ka I, Eliá? M. Was the mitochondrion necessary to start eukaryogenesis?, 2019, 27(2): 96–104.

    [117] Pittis AA, Gabaldón T. Late acquisition of mitochondria by a host with chimaeric prokaryotic ancestry., 2016, 531(7592): 101–104.

    [118] Ettema TJ. Evolution: Mitochondria in the second act., 2016, 531(7592): 39–40.

    [119] Martin WF, Roettger M, Ku C, Garg SG, Nelson-Sathi S, Landan G. Late Mitochondrial origin is an artifact., 2017, 9(2): 373–379.

    [120] Degli Esposti M. Phylogeny and Evolution of Bacteria and Mitochondria. CRC Press, 2018.

    [121] Ku C, Nelson-Sathi S, Roettger M, Sousa FL, Lockhart PJ, Bryant D, Hazkani-Covo E, McInerney JO, Landan G, Martin WF. Endosymbiotic origin and differential loss of eukaryotic genes., 2015, 524(7566): 427– 432.

    [122] Andersson SG, Zomorodipour A, Andersson JO, Sicheritz-Pontén T, Alsmark UCM, Podowski RM, N?slund AK, Eriksson AS, Winkler HH, Kurland CG. The genome sequence of Rickettsia prowazekii and the origin of mitochondria., 1998, 396(6707): 133– 140.

    [123] Fitzpatrick DA, Creevey CJ, McInerney JO. Genome phylogenies indicate a meaningful α-proteobacterial phylogeny and support a grouping of the mitochondria with the Rickettsiales., 2006, 23(1): 74–85.

    [124] Williams KP, Sobral BW, Dickerman AW. A robust species tree for the alphaproteobacteria., 2007, 189(13): 4578–4586.

    [125] Thrash JC, Boyd A, Huggett MJ, Grote J, Carini P, Yoder RJ, Robbertse B, Spatafora JW, Rappé MS, Giovannoni SJ. Phylogenomic evidence for a common ancestor of mitochondria and the SAR11 clade., 2011, 1: 13.

    [126] Brindefalk B, Ettema TJ, Viklund J, Thollesson M, Andersson SG. A phylometagenomic exploration of oceanic alphaproteobacteria reveals mitochondrial relatives unrelated to the SAR11 clade., 2011, 6(9): e24457.

    [127] Rodríguez-Ezpeleta N, Embley TM. The SAR11 group of alpha-proteobacteria is not related to the origin of mitochondria., 2012, 7(1): e30520.

    [128] Wang Z, Wu M. An integrated phylogenomic approach toward pinpointing the origin of mitochondria., 2015, 5: 7949.

    [129] Wang Z, Wu M. Phylogenomic reconstruction indicates mitochondrial ancestor was an energy parasite., 2014, 9(10): e110685.

    [130] Martijn J, Vosseberg J, Guy L, Offre P, Ettema TJ. Deep mitochondrial origin outside the sampled alphaproteo-bacteria., 2018, 557(7703): 101–105.

    [131] Fan L, Wu DF, Goremykin V, Xiao J, Xu YB, Garg S, Zhang CL, Martin WF, Zhu RX. Mitochondria branch within Alphaproteobacteria., 2019, 715870.

    [132] Le PT, Pontarotti P, Raoult D. Alphaproteobacteria species as a source and target of lateral sequence transfers., 2014, 22(3): 147–156.

    [133] Janou?kovec J, Tikhonenkov DV, Burki F, Howe AT, Rohwer FL, Mylnikov AP, Keeling PJ. A new lineage of eukaryotes illuminates early mitochondrial genome reduction., 2017, 27(23): 3717–3724.e5.

    [134] Mu?oz-Gómez SA, Hess S, Burger G, Lang BF, Susko E, Slamovits CH, Roger AJ. An updated phylogeny of the Alphaproteobacteria reveals that the parasitic Rickettsiales and Holosporales have independent origins., 2019, 8: e42535.

    [135] Searcy DG. Metabolic integration during the evolutionary origin of mitochondria., 2003, 13(4): 229–238.

    [136] Lake JA, Rivera MC. Was the nucleus the first endosymbiont?, 1994, 91(8): 2880.

    [137] Godde JS. Breaking through a phylogenetic impasse: a pair of associated archaea might have played host in the endosymbiotic origin of eukaryotes., 2012, 2(1): 29.

    [138] Forterre P, Ga?a M. Giant viruses and the origin of modern eukaryotes., 2016, 31: 44–49.

    [139] Bell PJL. The viral eukaryogenesis hypothesis: a key role for viruses in the emergence of eukaryotes from a prokaryotic world environment., 2009, 1178(1): 91–105.

    [140] Sousa FL, Neukirchen S, Allen JF, Lane N, Martin WF. Lokiarchaeon is hydrogen dependent., 2016, 1(5): 1–3.

    [141] Liu Y, Zhou ZC, Pan J, Baker BJ, Gu JD, Li M. Comparative genomic inference suggests mixotrophic lifestyle for Thorarchaeota., 2018, 12(4): 1021–1031.

    [142] Manoharan L, Kozlowski JA, Murdoch RW, L?ffler FE, Sousa FL, Schleper C. Metagenomes from coastal marine sediments give insights into the ecological role and cellular features of Loki-and Thorarchaeota., 2019, 10(5): e02039–02019.

    [143] Pushkarev A, Inoue K, Larom S, Flores-Uribe J, Singh M, Konno M, Tomida S, Ito S, Nakamura R, Tsunoda SP, Philosof A, Sharon I, Yutin N, Koonin EV, Kandori H, Béjà O. A distinct abundant group of microbial rhodopsins discovered using functional metagenomics., 2018, 558(7711): 595–599.

    [144] Bulzu PA, Andrei A?, Salcher MM, Mehrshad M, Inoue K, Kandori H, Beja O, Ghai R, Banciu HL. Casting light on Asgardarchaeota metabolism in a sunlit microoxic niche., 2019, 4(7): 1129–1137.

    [145] Spang A, Stairs CW, Dombrowski N, Eme L, Lombard J, Caceres EF, Greening C, Baker BJ, Ettema TJG. Proposal of the reverse flow model for the origin of the eukaryotic cell based on comparative analyses of Asgard archaeal metabolism., 2019, 4(7): 1138–1148.

    [146] Sagha? A, Gutiérrez-Preciado A, Deschamps P, Moreira D, Bertolino P, Ragon M, López-García P. Unveiling microbial interactions in stratified mat communities from a warm saline shallow pond., 2017, 19(6): 2405–2421.

    [147] Ak?l C, Robinson RC. Genomes of Asgard archaea encode profilins that regulate actin., 2018, 562(7727): 439–443.

    [148] Imachi H, Aoi K, Tasumi E, Saito Y, Yamanaka Y, Saito Y, Yamaguchi T, Tomaru H, Takeuchi R, Morono Y, Inagaki F, Takai K. Cultivation of methanogenic com-munity from subseafloor sediments using a continuous- flow bioreactor., 2011, 5(12): 1913–1925.

    [149] Sz?ll?si GJ, Tannier E, Daubin V, Boussau B. The inference of gene trees with species trees., 2015, 64(1): e42–e62.

    [150] Hu GB. Whole cell cryo-electron tomography suggests mitochondria divide by budding., 2014, 20(4): 1180–1187.

    [151] Mehta AP, Supekova L, Chen JH, Pestonjamasp K, Webster P, Ko Y, Henderson SC, McDermott G, Supek F, Schultz PG. Engineering yeast endosymbionts as a step toward the evolution of mitochondria., 2018, 115(46): 11796–11801.

    Progress in elucidating the origin of eukaryotes

    Zhiwei Gao1, Long Wang2

    Knowledge of the origin of eukaryotes is key to broadening our understanding of the eukaryotic genome and the relationship among internal structures within a eukaryotic cell. Since the discovery of archaea in 1977 and the proposal of three-domain tree of life by the American microbiologist Carl Woese, the intimate relationship in evolution between eukaryotes and archaea has been demonstrated by considerable experiments and analyses. From the beginning of the 21st century, with the development of phylogenetic methods and the discovery of new archaeal phyla more related to eukaryotes, increasing evidence has shown that Eukarya and Archaea should be merged into one domain, leading to a two-domain tree of life. Nowadays, the Asgard superphylum discovered via metagenomic analysis is regarded as the closest prokaryotes to eukaryotes. Nevertheless, several key questions are still under debate, such as what the ancestors of the eukaryotes were and when mitochondria emerged. Here, we review the current research progress regarding the changes of the tree of life and the detailed eukaryotic evolutionary mechanism. We show that the recent findings have greatly improved our knowledge on the origin of eukaryotes, which will pave the way for future studies.

    eukaryogenesis; tree of life; archaea; three domains; two domains; endosymbiosis; mitochondria; Asgard superphylum

    2020-04-18;

    2020-05-06

    國家自然科學(xué)基金項(xiàng)目(編號:31970236)資助[Supported by the National Natural Science Foundation of China (No. 31970236)]

    高志偉,在讀本科生,專業(yè)方向:生物化學(xué)與分子生物學(xué)。E-mail: 171240503@smail.nju.edu.cn

    王龍,博士,副教授,研究方向:分子遺傳與進(jìn)化。E-mail: wanglong@nju.edu.cn

    10.16288/j.yczz.20-107

    2020/6/5 10:08:53

    URI: http://kns.cnki.net/kcms/detail/11.1913.R.20200604.1524.001.html

    (責(zé)任編委: 劉鋼)

    猜你喜歡
    真核宿主線粒體
    棘皮動(dòng)物線粒體基因組研究進(jìn)展
    真核翻譯起始因子-5A2在肝內(nèi)膽管癌中的表達(dá)及意義
    線粒體自噬與帕金森病的研究進(jìn)展
    病原體與自然宿主和人的生態(tài)關(guān)系
    科學(xué)(2020年3期)2020-11-26 08:18:22
    龜鱉類不可能是新冠病毒的中間宿主
    表現(xiàn)為扁平苔蘚樣的慢性移植物抗宿主病一例
    人乳頭瘤病毒感染與宿主免疫機(jī)制
    NF-κB介導(dǎo)線粒體依賴的神經(jīng)細(xì)胞凋亡途徑
    人醛縮酶A干擾RNA真核表達(dá)載體的構(gòu)建
    真核翻譯起始因子5A2在胰腺癌中的表達(dá)及與預(yù)后的相關(guān)性
    秋霞在线观看毛片| 少妇被粗大的猛进出69影院 | 三级国产精品片| 乱码一卡2卡4卡精品| 欧美人与善性xxx| 久久99一区二区三区| 精品一区二区免费观看| 在线观看美女被高潮喷水网站| 亚洲精品日本国产第一区| 亚洲美女搞黄在线观看| 国产成人精品福利久久| 久久久久久久久久久久大奶| 亚洲欧美成人综合另类久久久| 国产视频首页在线观看| 黑人欧美特级aaaaaa片| xxxhd国产人妻xxx| 国产综合精华液| 51国产日韩欧美| 美女中出高潮动态图| videosex国产| 亚洲欧洲国产日韩| 性色avwww在线观看| 欧美丝袜亚洲另类| 午夜福利,免费看| 国产成人精品久久久久久| 成人国语在线视频| 色婷婷久久久亚洲欧美| 国产日韩欧美在线精品| 高清黄色对白视频在线免费看| av免费在线看不卡| 男男h啪啪无遮挡| 99热这里只有精品一区| 国产成人freesex在线| 91国产中文字幕| 亚洲av电影在线观看一区二区三区| 亚洲精品国产av成人精品| 国产在线一区二区三区精| 麻豆成人av视频| 2022亚洲国产成人精品| 日韩大片免费观看网站| 国产精品嫩草影院av在线观看| 精品亚洲成a人片在线观看| 亚洲av二区三区四区| 国产欧美日韩综合在线一区二区| 精品国产国语对白av| 亚洲精品国产av蜜桃| 亚洲在久久综合| 亚洲精品,欧美精品| 99九九线精品视频在线观看视频| 久热久热在线精品观看| 岛国毛片在线播放| 日本vs欧美在线观看视频| 男男h啪啪无遮挡| 亚洲精品乱码久久久v下载方式| 日韩,欧美,国产一区二区三区| 黄片无遮挡物在线观看| 亚洲色图综合在线观看| 乱人伦中国视频| 色94色欧美一区二区| 精品久久久久久久久亚洲| 青青草视频在线视频观看| 久久精品国产亚洲网站| 九草在线视频观看| 少妇人妻久久综合中文| 99久久精品国产国产毛片| 午夜福利视频在线观看免费| 欧美三级亚洲精品| 高清av免费在线| av女优亚洲男人天堂| 日韩视频在线欧美| 亚洲精品视频女| 国产男女内射视频| 下体分泌物呈黄色| 欧美老熟妇乱子伦牲交| 欧美日韩精品成人综合77777| 日本av手机在线免费观看| 婷婷成人精品国产| 欧美日韩综合久久久久久| 搡女人真爽免费视频火全软件| 男男h啪啪无遮挡| 亚洲人成网站在线播| 久久久久国产网址| 日日撸夜夜添| 老司机影院成人| 日韩视频在线欧美| 交换朋友夫妻互换小说| 久久午夜综合久久蜜桃| 少妇人妻精品综合一区二区| 熟妇人妻不卡中文字幕| 一级爰片在线观看| 国产色婷婷99| 岛国毛片在线播放| 亚洲精品日韩在线中文字幕| 亚洲怡红院男人天堂| 少妇的逼好多水| 91精品国产国语对白视频| 麻豆乱淫一区二区| 久久精品熟女亚洲av麻豆精品| 大陆偷拍与自拍| 日产精品乱码卡一卡2卡三| 黑人高潮一二区| 欧美日韩视频高清一区二区三区二| 日本欧美国产在线视频| 婷婷成人精品国产| 久久精品熟女亚洲av麻豆精品| 欧美老熟妇乱子伦牲交| 国产免费又黄又爽又色| 乱人伦中国视频| 亚洲国产日韩一区二区| 天天躁夜夜躁狠狠久久av| 十分钟在线观看高清视频www| 在线观看免费高清a一片| 秋霞在线观看毛片| 国产av一区二区精品久久| 寂寞人妻少妇视频99o| 男人操女人黄网站| 七月丁香在线播放| 国产国语露脸激情在线看| 久久99精品国语久久久| 交换朋友夫妻互换小说| 寂寞人妻少妇视频99o| 最后的刺客免费高清国语| 丝袜脚勾引网站| 久久精品国产亚洲av天美| 亚洲欧美色中文字幕在线| 乱码一卡2卡4卡精品| 能在线免费看毛片的网站| 日韩精品有码人妻一区| 亚洲国产精品999| 亚洲四区av| 97超视频在线观看视频| 亚洲不卡免费看| 国产国拍精品亚洲av在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲综合色惰| 大香蕉久久网| 久热久热在线精品观看| 国产综合精华液| 大香蕉久久成人网| 国产精品国产av在线观看| 免费日韩欧美在线观看| 99久久中文字幕三级久久日本| 九草在线视频观看| 国产一级毛片在线| 欧美性感艳星| 国产极品天堂在线| 久热这里只有精品99| 国产精品嫩草影院av在线观看| 亚洲精品乱久久久久久| 国产有黄有色有爽视频| 久久久久久久久久久丰满| 美女福利国产在线| av福利片在线| a级毛片在线看网站| 亚洲精品色激情综合| 日日摸夜夜添夜夜爱| 69精品国产乱码久久久| 亚洲精品美女久久av网站| 国产精品久久久久久久久免| 精品少妇久久久久久888优播| 日本与韩国留学比较| kizo精华| 三级国产精品片| 9色porny在线观看| 蜜臀久久99精品久久宅男| 国产av一区二区精品久久| 日韩 亚洲 欧美在线| 国产精品欧美亚洲77777| 日韩av在线免费看完整版不卡| 99久久精品一区二区三区| 日韩 亚洲 欧美在线| 亚洲欧洲精品一区二区精品久久久 | 麻豆成人av视频| 亚洲av综合色区一区| 天堂8中文在线网| 91aial.com中文字幕在线观看| 欧美人与善性xxx| 亚洲国产成人一精品久久久| 免费观看无遮挡的男女| 久久精品久久久久久噜噜老黄| 国产精品一区二区三区四区免费观看| 中文字幕免费在线视频6| 成人无遮挡网站| 久久久久久久久大av| 高清午夜精品一区二区三区| 91成人精品电影| 欧美 亚洲 国产 日韩一| 亚洲av福利一区| 欧美精品人与动牲交sv欧美| 欧美97在线视频| 欧美xxxx性猛交bbbb| 国产精品欧美亚洲77777| 这个男人来自地球电影免费观看 | 免费观看无遮挡的男女| 日本vs欧美在线观看视频| 亚洲国产精品一区二区三区在线| 亚洲欧美色中文字幕在线| 精品少妇黑人巨大在线播放| 国产欧美日韩一区二区三区在线 | 国产午夜精品一二区理论片| 亚洲精品第二区| √禁漫天堂资源中文www| 日本与韩国留学比较| www.av在线官网国产| 亚洲精品一区蜜桃| 五月天丁香电影| 日韩强制内射视频| 欧美日韩视频精品一区| 波野结衣二区三区在线| 色婷婷av一区二区三区视频| 大片免费播放器 马上看| 国内精品宾馆在线| 99热6这里只有精品| 啦啦啦啦在线视频资源| 久久久久久久大尺度免费视频| 亚洲图色成人| 欧美亚洲日本最大视频资源| 春色校园在线视频观看| 黄片播放在线免费| 五月玫瑰六月丁香| 国产 一区精品| 欧美日韩一区二区视频在线观看视频在线| 一区二区日韩欧美中文字幕 | 一区二区av电影网| 男女免费视频国产| 国产精品99久久久久久久久| 亚洲少妇的诱惑av| 色94色欧美一区二区| 日韩精品有码人妻一区| 高清av免费在线| 搡老乐熟女国产| 国产成人免费无遮挡视频| 麻豆精品久久久久久蜜桃| 久久久久久人妻| 黑人巨大精品欧美一区二区蜜桃 | 秋霞在线观看毛片| 中文字幕精品免费在线观看视频 | 免费播放大片免费观看视频在线观看| 亚洲av日韩在线播放| 亚洲精品美女久久av网站| 青春草视频在线免费观看| 国产一级毛片在线| 美女国产高潮福利片在线看| 国产淫语在线视频| 黑人巨大精品欧美一区二区蜜桃 | 精品人妻熟女av久视频| 高清不卡的av网站| 色哟哟·www| 精品一区二区免费观看| 免费人妻精品一区二区三区视频| 久久国产精品大桥未久av| a级片在线免费高清观看视频| 国产在线一区二区三区精| 日韩一区二区三区影片| 高清欧美精品videossex| 国产综合精华液| 国产成人精品婷婷| av免费观看日本| 美女大奶头黄色视频| .国产精品久久| kizo精华| 国产精品国产三级国产专区5o| 国产成人av激情在线播放 | 国产 一区精品| 日韩一本色道免费dvd| 在线天堂最新版资源| 精品少妇黑人巨大在线播放| 在线观看免费视频网站a站| 精品人妻一区二区三区麻豆| 一个人免费看片子| 天堂俺去俺来也www色官网| 女性被躁到高潮视频| 免费观看无遮挡的男女| 亚洲av电影在线观看一区二区三区| 五月玫瑰六月丁香| 亚洲国产av影院在线观看| 91aial.com中文字幕在线观看| 天天操日日干夜夜撸| 女人久久www免费人成看片| 97超碰精品成人国产| 久久久a久久爽久久v久久| 狠狠精品人妻久久久久久综合| 久久精品夜色国产| 一级爰片在线观看| 老熟女久久久| 国产一级毛片在线| 人妻系列 视频| 精品人妻熟女av久视频| 一本大道久久a久久精品| 男男h啪啪无遮挡| 91在线精品国自产拍蜜月| 国产精品一区二区在线不卡| 国产国拍精品亚洲av在线观看| 亚洲av.av天堂| 日本wwww免费看| 制服诱惑二区| 大香蕉久久成人网| 亚洲精品视频女| 国产精品 国内视频| 大话2 男鬼变身卡| 人体艺术视频欧美日本| 亚洲av二区三区四区| 国产欧美亚洲国产| 99re6热这里在线精品视频| 亚洲精华国产精华液的使用体验| 精品视频人人做人人爽| 亚洲第一av免费看| 一个人免费看片子| 丁香六月天网| 日本91视频免费播放| 久久精品国产鲁丝片午夜精品| 99热这里只有精品一区| 中文字幕亚洲精品专区| 18禁动态无遮挡网站| 久久久久久久大尺度免费视频| 国产一级毛片在线| 亚洲精品456在线播放app| 国产视频内射| 99热全是精品| 99热这里只有是精品在线观看| 韩国av在线不卡| 国产精品国产三级国产专区5o| 亚洲精品成人av观看孕妇| 51国产日韩欧美| 99热6这里只有精品| 一级黄片播放器| 秋霞在线观看毛片| 99热这里只有是精品在线观看| 三级国产精品片| 另类精品久久| 91精品国产九色| 国产欧美日韩一区二区三区在线 | 国产女主播在线喷水免费视频网站| 三级国产精品欧美在线观看| 国产成人精品无人区| 大片电影免费在线观看免费| 久久精品国产亚洲网站| 人妻 亚洲 视频| 国产免费福利视频在线观看| 日日撸夜夜添| 99久久综合免费| 久久久久久久大尺度免费视频| 波野结衣二区三区在线| 国产成人精品久久久久久| 亚洲精品自拍成人| 男女啪啪激烈高潮av片| 2022亚洲国产成人精品| av免费观看日本| 色婷婷av一区二区三区视频| 亚州av有码| 久久久精品94久久精品| 视频中文字幕在线观看| 黑人猛操日本美女一级片| 久久人人爽av亚洲精品天堂| 国精品久久久久久国模美| 能在线免费看毛片的网站| av女优亚洲男人天堂| 日韩一区二区视频免费看| freevideosex欧美| 91久久精品国产一区二区三区| 久久久久精品性色| 青春草视频在线免费观看| 满18在线观看网站| 亚洲国产欧美在线一区| 高清欧美精品videossex| 80岁老熟妇乱子伦牲交| 成年人免费黄色播放视频| 久久人妻熟女aⅴ| freevideosex欧美| 一区二区av电影网| 亚洲欧洲日产国产| 最近中文字幕高清免费大全6| 精品久久久精品久久久| 黄色视频在线播放观看不卡| 美女视频免费永久观看网站| 精品亚洲乱码少妇综合久久| 亚洲精品aⅴ在线观看| 99热网站在线观看| 国产男人的电影天堂91| 久久久久人妻精品一区果冻| 啦啦啦视频在线资源免费观看| 高清欧美精品videossex| 国产精品国产三级国产专区5o| 自线自在国产av| 久热这里只有精品99| 少妇高潮的动态图| 亚洲精品aⅴ在线观看| 观看av在线不卡| 一级片'在线观看视频| xxxhd国产人妻xxx| 毛片一级片免费看久久久久| 一本久久精品| 一边摸一边做爽爽视频免费| 国产深夜福利视频在线观看| 在线观看免费视频网站a站| 国产成人精品福利久久| 日韩三级伦理在线观看| 久久精品国产亚洲网站| 亚洲人成网站在线观看播放| 国产精品.久久久| 日韩人妻高清精品专区| 视频区图区小说| 91久久精品国产一区二区三区| 国产免费一区二区三区四区乱码| 99热6这里只有精品| 乱码一卡2卡4卡精品| 亚洲色图 男人天堂 中文字幕 | 99re6热这里在线精品视频| 亚洲综合色网址| 大片免费播放器 马上看| 秋霞在线观看毛片| 在线观看免费日韩欧美大片 | 亚洲不卡免费看| 久久人妻熟女aⅴ| 看非洲黑人一级黄片| 国产精品免费大片| 久久韩国三级中文字幕| 少妇熟女欧美另类| 欧美精品高潮呻吟av久久| 久久久久国产网址| 亚洲激情五月婷婷啪啪| 我的女老师完整版在线观看| 18禁在线播放成人免费| 国产精品嫩草影院av在线观看| 晚上一个人看的免费电影| 久久这里有精品视频免费| 老司机影院毛片| 日本av免费视频播放| 午夜久久久在线观看| 99九九在线精品视频| 丝袜在线中文字幕| 精品人妻在线不人妻| 九九在线视频观看精品| 精品久久久久久久久亚洲| 国产亚洲午夜精品一区二区久久| 一区二区三区免费毛片| 国产无遮挡羞羞视频在线观看| 夜夜爽夜夜爽视频| 免费av不卡在线播放| 九色成人免费人妻av| 国产乱人偷精品视频| 精品久久久久久久久亚洲| 天堂俺去俺来也www色官网| 亚洲av二区三区四区| 制服诱惑二区| 青春草亚洲视频在线观看| 国产精品嫩草影院av在线观看| 亚洲精品一区蜜桃| 国产在线免费精品| 制服诱惑二区| 中国国产av一级| 最近2019中文字幕mv第一页| 免费高清在线观看日韩| 国产男人的电影天堂91| 欧美日韩国产mv在线观看视频| 久久精品国产亚洲网站| 日韩欧美一区视频在线观看| 成年美女黄网站色视频大全免费 | 丝袜脚勾引网站| 久久国产精品男人的天堂亚洲 | 大片免费播放器 马上看| 久久久精品区二区三区| 日韩一区二区视频免费看| 我的老师免费观看完整版| 久久久久久久久大av| 99久久综合免费| 国产深夜福利视频在线观看| 成年人免费黄色播放视频| 亚洲熟女精品中文字幕| 黄片播放在线免费| 麻豆精品久久久久久蜜桃| 一区二区日韩欧美中文字幕 | 99热6这里只有精品| 赤兔流量卡办理| 亚洲精品自拍成人| av一本久久久久| 人妻少妇偷人精品九色| h视频一区二区三区| 欧美精品人与动牲交sv欧美| 精品久久久久久久久av| 国产精品一区www在线观看| tube8黄色片| 制服诱惑二区| 人妻少妇偷人精品九色| 99国产综合亚洲精品| 亚洲av福利一区| 中文字幕制服av| 男的添女的下面高潮视频| 国产成人免费无遮挡视频| 秋霞伦理黄片| 亚洲国产av新网站| 国产精品秋霞免费鲁丝片| 日本vs欧美在线观看视频| 中文字幕av电影在线播放| 99久久人妻综合| 一级毛片黄色毛片免费观看视频| 精品一区二区免费观看| 一本大道久久a久久精品| 美女中出高潮动态图| 欧美精品人与动牲交sv欧美| 国产精品欧美亚洲77777| 国产精品国产三级国产专区5o| 99国产精品免费福利视频| 最黄视频免费看| 欧美精品亚洲一区二区| 纵有疾风起免费观看全集完整版| 在线观看美女被高潮喷水网站| 在线亚洲精品国产二区图片欧美 | 久久精品国产自在天天线| 久久韩国三级中文字幕| 国产亚洲精品第一综合不卡 | 母亲3免费完整高清在线观看 | 精品99又大又爽又粗少妇毛片| 亚洲av欧美aⅴ国产| 国产免费一区二区三区四区乱码| 亚洲无线观看免费| 国产精品国产三级国产av玫瑰| 国产精品蜜桃在线观看| 国产精品一区二区三区四区免费观看| 极品人妻少妇av视频| 亚洲国产精品成人久久小说| 欧美激情 高清一区二区三区| 我的女老师完整版在线观看| 黑人巨大精品欧美一区二区蜜桃 | 国产午夜精品一二区理论片| 一区在线观看完整版| 精品久久久精品久久久| 秋霞伦理黄片| 观看av在线不卡| 99久久精品国产国产毛片| 午夜免费男女啪啪视频观看| 男人操女人黄网站| 日韩免费高清中文字幕av| 午夜视频国产福利| 看非洲黑人一级黄片| 亚洲av免费高清在线观看| 亚洲欧美成人精品一区二区| 国产国语露脸激情在线看| 成年人午夜在线观看视频| 99re6热这里在线精品视频| 我的女老师完整版在线观看| 在线免费观看不下载黄p国产| 精品人妻熟女毛片av久久网站| 国产成人免费观看mmmm| 亚洲欧美日韩另类电影网站| 亚洲国产日韩一区二区| 中文精品一卡2卡3卡4更新| 成人二区视频| 亚洲第一区二区三区不卡| 9色porny在线观看| 久久久国产欧美日韩av| 亚洲精品视频女| 欧美日韩国产mv在线观看视频| 国产成人午夜福利电影在线观看| 一区二区av电影网| 午夜激情福利司机影院| 成人亚洲精品一区在线观看| 一边亲一边摸免费视频| 18禁动态无遮挡网站| 男的添女的下面高潮视频| 久久国产亚洲av麻豆专区| 秋霞伦理黄片| 国产一区二区在线观看av| 韩国av在线不卡| 插逼视频在线观看| 美女xxoo啪啪120秒动态图| av在线老鸭窝| 一区二区三区精品91| 一本大道久久a久久精品| 99re6热这里在线精品视频| 精品午夜福利在线看| 女的被弄到高潮叫床怎么办| 亚洲第一区二区三区不卡| 久久青草综合色| 日韩欧美一区视频在线观看| 精品久久久久久电影网| 青春草国产在线视频| 亚洲欧美精品自产自拍| 久久青草综合色| 久久毛片免费看一区二区三区| 成人国产av品久久久| 激情五月婷婷亚洲| 乱人伦中国视频| 亚洲色图 男人天堂 中文字幕 | 美女大奶头黄色视频| 国语对白做爰xxxⅹ性视频网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品偷伦视频观看了| 国产片特级美女逼逼视频| 夫妻性生交免费视频一级片| 美女xxoo啪啪120秒动态图| 欧美97在线视频| 少妇人妻久久综合中文| 欧美成人精品欧美一级黄| freevideosex欧美| 国产精品嫩草影院av在线观看| 亚洲中文av在线| 久久99热这里只频精品6学生| 国产高清不卡午夜福利| av卡一久久| 免费看光身美女| 青春草亚洲视频在线观看| 一级毛片aaaaaa免费看小| 人人妻人人澡人人看| 热re99久久精品国产66热6| 亚洲精品乱码久久久久久按摩| 2021少妇久久久久久久久久久| 九草在线视频观看| 青青草视频在线视频观看| 亚洲人成网站在线观看播放| 亚洲综合色网址| 又黄又爽又刺激的免费视频.| 又大又黄又爽视频免费| 婷婷成人精品国产| 久久99蜜桃精品久久|