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

    分子絕對進(jìn)化速率與物種分歧時(shí)間之間的定量關(guān)系

    2015-03-22 02:19:29李可群
    生物學(xué)雜志 2015年2期
    關(guān)鍵詞:核苷酸定量位點(diǎn)

    李可群

    (同濟(jì)大學(xué) 化學(xué)系, 上海 200092)

    分子絕對進(jìn)化速率與物種分歧時(shí)間之間的定量關(guān)系

    李可群

    (同濟(jì)大學(xué) 化學(xué)系, 上海 200092)

    分子進(jìn)化;絕對進(jìn)化速率;分子鐘;定量關(guān)系;系統(tǒng)發(fā)育分析

    系統(tǒng)發(fā)育研究不僅有助重建地球上所有生物體的進(jìn)化歷史,而且還可以揭示生物學(xué)領(lǐng)域的一些基本問題。清晰了解各生物物種進(jìn)化歷程及不同物種之間的進(jìn)化關(guān)系,是進(jìn)一步研究和探索生物學(xué)其他學(xué)科的基礎(chǔ)[1]。

    分子系統(tǒng)發(fā)育分析大多基于這樣一種事實(shí)基礎(chǔ): 在各種不同的發(fā)育譜系(lineages)及足夠大的進(jìn)化時(shí)間尺度中,許多氨基酸或核苷酸序列位點(diǎn)的進(jìn)化速率是幾乎恒定不變的,即基于所謂“分子鐘”概念[2]。但在實(shí)際工作中,人們發(fā)現(xiàn)多數(shù)氨基酸和核苷酸序列的“分子鐘”都有一定的缺陷, 早在“分子鐘”概念提出不久的1971年,Ohta和Kimura就指出分子進(jìn)化速率并不恒定。相反在不同支系之間存在差異,甚至在同一支系不同進(jìn)化時(shí)間分子進(jìn)化速率也有變化[3-6]。而對于分子絕對進(jìn)化速率與物種分歧時(shí)間之間的變化關(guān)系,不同文獻(xiàn)給出了不同甚至是相反的結(jié)論。賀福初[7]指出物種分子絕對進(jìn)化速率隨物種分歧時(shí)間的減小而下降,即存在“減速進(jìn)化”現(xiàn)象,文獻(xiàn)[8]也指出高等分類元存在進(jìn)化速率減慢現(xiàn)象。同時(shí)實(shí)際工作中也發(fā)現(xiàn)從魚類、兩棲類、鳥類再到哺乳類,線粒體的分子進(jìn)化速率是增速的[9]。

    準(zhǔn)確地估計(jì)物種分歧時(shí)間和推斷其進(jìn)化歷史,是分子進(jìn)化遺傳學(xué)和系統(tǒng)發(fā)育學(xué)的一個(gè)重要研究課題[10]。因此了解生物學(xué)體系中分子絕對進(jìn)化速率與物種分歧時(shí)間之間的定量關(guān)系顯得十分重要。本文將利用一些常見“分子鐘”分子在這一領(lǐng)域進(jìn)行初步的探討。

    1 血紅蛋白α鏈體系的研究

    1.1 血紅蛋白α鏈絕對進(jìn)化速率的計(jì)算

    若假定分子進(jìn)化速率不變,本文作者給出了分子絕對進(jìn)化速率的計(jì)算公式[11]:

    (1)

    式(1)中,k1、k2分別為兩個(gè)同源序列的分子絕對進(jìn)化速率,t為從共同祖先序列進(jìn)化的時(shí)間,或稱兩個(gè)物種的分歧時(shí)間,p為序列比對所得的種間差異率,而k1t、k2t分別為兩個(gè)同源序列的平均突變概率。d1、d2分別是兩個(gè)同源序列與共同祖先序列相比的差異位點(diǎn)數(shù),d為兩個(gè)同源序列中相同位點(diǎn)氨基酸(或核苷酸)不相同的位點(diǎn)數(shù)。n0為被比較同源序列氨基酸或核苷酸位點(diǎn)數(shù)。為了表述和計(jì)算的方便, 本文中分子絕對進(jìn)化速率的單位采用氨基酸(或核苷酸)位點(diǎn)數(shù)/10億年,進(jìn)化時(shí)間(或物種分歧時(shí)間)的單位采用10億年。

    圖1 8種脊椎動(dòng)物的系統(tǒng)進(jìn)化示意圖

    我們首先計(jì)算了Kimura在其所著的《Theneutraltheoryofmolecularevolution》一書[12]中所舉的用來說明分子進(jìn)化速率恒定的經(jīng)典體系,即8種脊椎動(dòng)物(鯊魚、鯉魚、蟾螈、雞、針鼴鼠、袋鼠、狗和人)中血紅蛋白α鏈的絕對進(jìn)化速率。為了使用最新數(shù)據(jù),我們從美國國家生物技術(shù)信息中心(NCBI)GenBank數(shù)據(jù)庫中檢索了該8種脊椎動(dòng)物的血紅蛋白α鏈序列數(shù)據(jù),它們的Accessionnumber分別為鯊魚(Callorhinchusmilii,AFM89055.1)、鯉魚(Cyrinuscarpio,AGJ70090.1)、蟾螈(Tarchagranulose, 0402193A)、雞(GallusGallus,NP_001004376.1)、針鼴鼠(Tachyglossusaculeatus,P01977.1)、袋鼠(Macropusgiganteus,P01975.1)、狗(Canislupusfamiliaris,P60529.1)和人(Homosapiens,NP_000549.1)。這8種脊椎動(dòng)物的物種分歧時(shí)間來自Kimura的書,它們的系統(tǒng)進(jìn)化圖請參見圖1,圖中t1至t7分別取0.468、0.423、0.368、0.291、0.232、0.139和0.081,單位為10億年。

    圖中8種脊椎動(dòng)物的血紅蛋白α鏈序列數(shù)據(jù)使用NCBI提供的ProteinBlast軟件進(jìn)行對齊(align)并計(jì)算出物種之間的種間差異率,所得數(shù)據(jù)見表1。

    表1 8種脊椎動(dòng)物血紅蛋白α鏈的比較

    由于公式(1)直接求解較為困難,可采用非線性最優(yōu)化方法求解。若有n種同源序列兩兩進(jìn)行比較,則其目標(biāo)函數(shù)為

    (2)

    式(2)中ki、kj分別為被比較的第i種和第j種同源序列的分子絕對進(jìn)化速率,pij為它們的種間差異率,tij為它們的物種分歧時(shí)間。我們使用優(yōu)化軟件Lingo 11.0計(jì)算出8種脊椎動(dòng)物血紅蛋白α鏈的分子絕對進(jìn)化速率見表2。

    表2 8種脊椎動(dòng)物血紅蛋白α鏈的分子絕對進(jìn)化速率

    1.2 血紅蛋白α鏈分子絕對進(jìn)化速率與物種分歧時(shí)間的關(guān)系

    我們在研究分子絕對進(jìn)化速率與物種分歧時(shí)間之間的定量關(guān)系時(shí),發(fā)現(xiàn)如果將血紅蛋白α鏈分子絕對進(jìn)化速率與物種分歧時(shí)間用下列公式進(jìn)行擬合

    (3)

    式(3)中k、t分別為前述的分子絕對進(jìn)化速率和物種分歧時(shí)間,A、B為系數(shù),可發(fā)現(xiàn)明顯的線性關(guān)系,參見圖2。

    圖2 血紅蛋白α鏈分子絕對進(jìn)化速率與物種分歧時(shí)間之間的定量關(guān)系Fig 2 Quantitative relationship between hemoglobin α chain molecular

    注:圖中沿物種分歧時(shí)間軸從左至右各點(diǎn)依次為鯊魚、鯉魚、蟾螈、雞、針鼴鼠、袋鼠和人。

    1.3 使用分子絕對進(jìn)化速率與物種分歧時(shí)間之間的定量關(guān)系公式擬合數(shù)據(jù)

    為了進(jìn)一步驗(yàn)證上述定量關(guān)系并探究其產(chǎn)生的原因,我們使用式(3)重新擬合其中前7種脊椎動(dòng)物血紅蛋白α鏈的種間差異率數(shù)據(jù)。

    1.3.1 進(jìn)化路徑各節(jié)點(diǎn)區(qū)間絕對進(jìn)化速率之間的關(guān)系

    若某序列由共同祖先序列進(jìn)化的路徑包含n+1個(gè)不同節(jié)點(diǎn),各節(jié)點(diǎn)區(qū)間的分子絕對進(jìn)化速率和進(jìn)化時(shí)間分別為k1、k2、……、kn以及t1、t2、……、tn。由于序列分子絕對進(jìn)化速率計(jì)算公式由泊松分布及未突變概率e-kt推導(dǎo)而來[11-12],不難理解對于包含n段不同節(jié)點(diǎn)區(qū)間的序列進(jìn)化路徑,存在下列關(guān)系式:

    (4)

    1.3.2 7種脊椎動(dòng)物血紅蛋白α鏈的種間差異率數(shù)據(jù)的擬合

    由于血紅蛋白α鏈分子為經(jīng)典的“分子鐘”分子,其許多位于不同譜系的分子在較大時(shí)間尺度內(nèi)絕對進(jìn)化速率數(shù)值近似恒定(參見表2),因此不難理解在7個(gè)脊椎動(dòng)物體系中主干枝和各進(jìn)化分枝分別可近似視為一個(gè)“分子鐘”在起作用。以圖3示意的兩姊妹物種為例說明擬合過程。

    圖3 主干枝遺傳距離計(jì)算示意圖

    設(shè)物種i和j在它們各自進(jìn)化分枝中分子絕對進(jìn)化速率分別為ki、kj,而它們的物種分歧時(shí)間分別為ti、tj,則它們在各自進(jìn)化分枝時(shí)間段的平均突變概率為kiti和kjtj,將式(3)代入有

    (5)

    (6)

    式(5)和式(6)中,Ai、Aj以及Bi、Bj均為系數(shù),由于它們對應(yīng)著兩個(gè)不同的姊妹物種,因此系數(shù)數(shù)值上可能存在差異。而計(jì)算物種j時(shí)還需包含主干枝ti至tj時(shí)間段平均突變概率,該平均突變概率可用下列公式來計(jì)算

    m=k0j(t0-tj)-k0i(t0-ti)

    (7)

    式(7)中,t0為主干枝“分子鐘”的分歧時(shí)間。k0i和k0j分別為主干枝分子在時(shí)間t0至ti或tj時(shí)間段平均絕對進(jìn)化速率。將式(3)代入式(7)有

    (8)

    由于主干枝上可視為同一“分子鐘”,故不難理解式(8)中系數(shù)A0和B0使用相同數(shù)值。將式(8)中eB0用k0替代并移到自然指數(shù)項(xiàng)外有

    (9)

    實(shí)際計(jì)算表明,式(9)中兩個(gè)自然指數(shù)項(xiàng)上指數(shù)絕對值大多較小,因此可近似將兩個(gè)自然指數(shù)項(xiàng)分別用泰勒級數(shù)展開并取前兩項(xiàng)有

    (10)

    則物種j在被比較的進(jìn)化時(shí)間段內(nèi)的總平均突變概率為

    (11)

    將式(5)和式(11)所得的平均突變概率代入式(2)目標(biāo)函數(shù)中,同樣地我們使用優(yōu)化軟件Lingo11.0 重新擬合了圖1中前7種脊椎動(dòng)物的種間差異率數(shù)據(jù)??梢缘玫?種脊椎動(dòng)物在其各自的進(jìn)化分枝內(nèi)的分子絕對進(jìn)化速率等結(jié)果如表3所示,各物種的k0也為它們各自的系數(shù)B的自然指數(shù)eB。

    表3 7種脊椎動(dòng)物體系重新優(yōu)化的結(jié)果

    需注意的是, 由于優(yōu)化過程需滿足方程數(shù)大于變量數(shù)的條件,因此參與計(jì)算的物種需大于或等于6種。

    1.4 分子絕對進(jìn)化速率與物種分歧時(shí)間之間定量關(guān)系公式的物理意義

    我們知道化學(xué)動(dòng)力學(xué)中化學(xué)反應(yīng)速率常數(shù)k與熱力學(xué)溫度T的定量關(guān)系可用阿侖尼烏斯公式來表達(dá):

    (12)

    式(12)中,Ea為化學(xué)反應(yīng)活化能(即活化分子最低能量與反應(yīng)物平均能量的差值),R為氣體常數(shù),B為與反應(yīng)體系有關(guān)的系數(shù)。對比本文公式(3), 可以發(fā)現(xiàn)兩者公式形式非常相似,而事實(shí)上分子進(jìn)化理論中所稱的絕對進(jìn)化速率嚴(yán)格地應(yīng)稱為絕對進(jìn)化速率常數(shù),因此我們把公式(3)改寫為

    (13)

    式(13)中Ea為序列位點(diǎn)突變活化能,t為進(jìn)化時(shí)間或物種分歧時(shí)間,B為與物種有關(guān)的系數(shù),R為與物種無關(guān)的常數(shù)。Kimura認(rèn)為氨基酸和核苷酸位點(diǎn)突變?yōu)殡S機(jī)過程,即所研究的序列位點(diǎn)上氨基酸或核苷酸可視為獨(dú)立質(zhì)點(diǎn)而被其他外來氨基酸或核苷酸質(zhì)點(diǎn)隨機(jī)取代,本文作者利用這一思路提出了計(jì)算分子絕對進(jìn)化速率的公式,即本文公式(1), 上述血紅蛋白α鏈以及本文后面所舉的例子都表明其計(jì)算誤差很小,平均殘差僅為0.02~0.03左右甚至更小, 這說明將序列位點(diǎn)突變視為質(zhì)點(diǎn)隨機(jī)取代過程這種假設(shè)是可行的。而取代反應(yīng)也為化學(xué)反應(yīng),因此將本文公式(3)改寫為式(13)的形式是合理的。

    不過,序列位點(diǎn)的突變與化學(xué)反應(yīng)也存在區(qū)別,化學(xué)反應(yīng)中,系統(tǒng)熱力學(xué)溫度上升,反應(yīng)物平均能量增加,可發(fā)生化學(xué)反應(yīng)的活化分子百分?jǐn)?shù)增加,從而反應(yīng)速率常數(shù)變大;而序列位點(diǎn)突變中,隨著進(jìn)化時(shí)間的增加,突變位點(diǎn)數(shù)增加,序列會(huì)出現(xiàn)結(jié)構(gòu)張力發(fā)生改變等因素使體系能量變化,若體系能量上升,可發(fā)生突變的位點(diǎn)百分?jǐn)?shù)增加,使位點(diǎn)絕對突變速率亦即分子絕對進(jìn)化速率增加,其與進(jìn)化時(shí)間的定量關(guān)系公式形式為式(13);若突變位點(diǎn)數(shù)增加導(dǎo)致體系能量下降,可發(fā)生突變的位點(diǎn)百分?jǐn)?shù)減小,分子絕對進(jìn)化速率與進(jìn)化時(shí)間的定量關(guān)系公式形式為式(14)

    (14)

    將式(13)和式(14)改寫成指數(shù)形式,并用k0替代其中的eB項(xiàng),則有

    (15)

    (16)

    式(15)和式(16)中,當(dāng)進(jìn)化時(shí)間或物種分歧時(shí)間分別趨向無窮大時(shí),分子絕對進(jìn)化速率k將等于k0,即兩式中的k0分別為它們的分子極限絕對進(jìn)化速率。

    由于突變活化能一般為正值,我們不難由式(15)和式(16)來理解本文序言部分提及的隨著物種分歧時(shí)間的減小而出現(xiàn)分子絕對進(jìn)化速率增加和減小的兩種情形。對于同一進(jìn)化分枝而言,若屬于突變位點(diǎn)數(shù)增加時(shí)體系能量也增加,即屬于式(15)描述的情形,進(jìn)化時(shí)間的減小會(huì)導(dǎo)致分子絕對進(jìn)化速率減小;而若屬于突變位點(diǎn)數(shù)增加時(shí)而體系能量下減,即式(16)描述的情形,進(jìn)化時(shí)間的減小會(huì)導(dǎo)致分子絕對進(jìn)化速率的增加。對于從同一主干枝進(jìn)化出的姊妹物種,若它們分子的位點(diǎn)突變活化能相差不大,也會(huì)出現(xiàn)類似現(xiàn)象,即賀福初提到的所謂“減速進(jìn)化現(xiàn)象”以及所謂線粒體分子進(jìn)化的“加速進(jìn)化現(xiàn)象”。

    另外,物種進(jìn)化時(shí)可能會(huì)有序列的活化能數(shù)值的較大變化,從而出現(xiàn)分子絕對進(jìn)化速率的較大變化,即出現(xiàn)所謂物種進(jìn)化時(shí)“分子鐘”的突然變異[13-14]和物種進(jìn)化過程可能存在多個(gè)串聯(lián)存在的“分子鐘”的現(xiàn)象[15]。

    1.5 分子進(jìn)化過程中“雙重分子鐘”現(xiàn)象的揭示

    如果我們將表3中7種血紅蛋白α鏈的分子極限絕對進(jìn)化速率與物種分歧時(shí)間按式(3)相同形式的公式(17)作圖,也可發(fā)現(xiàn)明顯的線性關(guān)系,參見圖4。

    (17)

    圖4 血紅蛋白α鏈分子極限絕對進(jìn)化速率與物種分歧時(shí)間之間的定量關(guān)系

    由于物種分子極限絕對進(jìn)化速率與進(jìn)化時(shí)間或物種分歧時(shí)間之間存在上述定量關(guān)系,它們也應(yīng)存在活化能概念和所謂物種分子極限絕對進(jìn)化速率的極限值。用公式可表達(dá)為

    (18)

    式(18)中Eb為控制物種分子極限絕對進(jìn)化速率“分子鐘”的活化能,k00為分子極限絕對進(jìn)化速率的極限值。將式(18)代入式(15)有

    (19)

    式(19)中E0為兩個(gè)“分子鐘”的活化能之和。由該式可以看出,盡管生物分子存在序列位點(diǎn)突變和控制物種分子極限進(jìn)化速率進(jìn)化的兩個(gè)“分子鐘”,但表觀上看僅表現(xiàn)出一個(gè)“分子鐘”在起作用,且其極限分子絕對進(jìn)化速率為k00。不過,當(dāng)兩個(gè)活化能數(shù)值相差較大時(shí),數(shù)值較大的活化能因其對應(yīng)的指數(shù)項(xiàng)數(shù)值隨時(shí)間變化較小而可與數(shù)值較小的活化能在二次優(yōu)化過程中分別求解。

    由圖4的血紅蛋白α鏈的例子可以看出,7種脊椎動(dòng)物中存在多組相同的分子極限絕對進(jìn)化速率的進(jìn)化活化能,它們的分子極限絕對進(jìn)化速率的進(jìn)化極限值k00也各自分別相等,它們具體的生物學(xué)意義有待進(jìn)一步的研究。

    2 其他分子體系的研究

    為了不失一般性,我們還對其他分子體系進(jìn)行了研究,取得了相似的結(jié)果。

    2.1 8種脊椎動(dòng)物的血紅蛋白β鏈

    我們還計(jì)算了Kimura在書中為了說明分子進(jìn)化速率恒定性而舉的另一個(gè)經(jīng)典例子,即血紅蛋白β鏈。8種脊椎動(dòng)物血紅蛋白β鏈的序列數(shù)據(jù)同樣來自于NCBI的GenBank數(shù)據(jù)庫,它們的Accessionnumber分別為: 鯊魚(Callorhinchusmilii,AFM89007.1)、鯉魚(Cyrinuscarpio,P02139.1)、蟾螈(Triturusaristatus,P10785.1)、雞(GallusGallus,NP_990820.1)、針鼴鼠(Tachyglossusaculeatus,P02110.2)、袋鼠(Macropusgiganteus, 690945A)、狗(Canislupusfamiliaris,P60524.1)和人(Homosapiens,NP_000509.1)。同樣地使用NCBI提供的ProteinBlast軟件進(jìn)行對齊(align)并計(jì)算出物種之間的種間差異率,所得數(shù)據(jù)見表4。

    表4 8種脊椎動(dòng)物血紅蛋白β鏈的比較

    我們分別直接使用本文公式(1)和使用公式(3)對表4數(shù)據(jù)進(jìn)行了擬合,結(jié)果如下。

    表5 8種脊椎動(dòng)物血紅蛋白β鏈的優(yōu)化結(jié)果

    圖5 血紅蛋白β鏈分子極限絕對進(jìn)化速率與物種分歧時(shí)間之間的定量關(guān)系

    Fig5Quantitativerelationshipbetweenhemoglobinβchainmolecularextremeabsoluteevolutionaryratesandtheirtaxadivergencetimes

    血紅蛋白β鏈體系的計(jì)算結(jié)果表明生物分子進(jìn)化確實(shí)存在控制序列位點(diǎn)突變和物種分子極限絕對進(jìn)化速率進(jìn)化的“雙重分子鐘”現(xiàn)象。

    2.2 其他蛋白質(zhì)體系

    2.3 線粒體核苷酸序列

    我們同樣也研究了一些線粒體核苷酸序列,具體計(jì)算結(jié)果請見表6。

    表6 一些線粒體核苷酸序列計(jì)算結(jié)果

    計(jì)算過程可能使用的物種及其 accession number 分別為:鯊魚(Carcharchinus leucas, NC_023522.1)、鯉魚(Cyprinus carpio, KJ511883.1)、蟾螈(Triturus cristatus, NC_015790.1)、雞(Gallus Gallus, KF826490.1)、針鼴鼠(Tachyglossus aculeatus, NC_003321.1)、袋鼠(Lagostrophus fasciatus, NC_008447.1)和狗(Canis lupus familiaris, NC_008092.1)。

    由表6計(jì)算結(jié)果可以看出,所得平均殘差很小,說明擬合效果較好。同時(shí)還可以看出,各體系中均存在一組或多組很好的線性關(guān)系,說明“雙重分子鐘”也存在于線粒體核苷酸序列中。

    我們使用本文提出的公式擬合了一些常見“分子鐘”分子在7個(gè)脊椎動(dòng)物體系中的蛋白質(zhì)和核苷酸序列比對數(shù)據(jù),它們的平均偏差僅為0.02~0.03左右甚至更小,因此我們認(rèn)為本文提出的生物分子絕對進(jìn)化速率與物種分歧時(shí)間之間的定量計(jì)算公式是符合實(shí)驗(yàn)事實(shí)的。并且“雙重分子鐘”現(xiàn)象可能廣泛存在于各種蛋白質(zhì)和核苷酸分子進(jìn)化體系。

    本文給出的公式和結(jié)論可為理解多數(shù)分子的進(jìn)化速率在不同譜系甚至同一譜系不同進(jìn)化時(shí)間中存在差異的現(xiàn)象提供幫助,同時(shí)也可為解決包括令人困擾的早期生物進(jìn)化在內(nèi)的生物系統(tǒng)發(fā)育問題提供新的思路和方法。另外,本文還提出了分子位點(diǎn)突變活化能的概念,為研究生物物種進(jìn)化過程中的能量變化提供了可能。

    [1]黎一葦,于 黎,張亞平. 系統(tǒng)發(fā)育研究中的“長枝吸收”假象概述[J]. 遺傳, 2007, 29(6): 659-667.

    [2]潘星華,傅繼梁. 基因的分子進(jìn)化:原理與方法[J]. 自然雜志, 1995, 17(4):189-193.

    [3]柯葉艷,齊文同. 前寒武紀(jì)生物起源時(shí)間的化石和分子鐘研究[J]. 地質(zhì)論評, 2002, 48(5): 457-462.

    [4]Wray G A, Levinton J S, Shapiro L H. Molecular evidence for deep precambrian divergences among Metazoan phyla[J]. Science, 1996, 274: 568-573.

    [5]Li W H, Wu C I. Rates of nucleotide substitution are higher in rodent than in man[J]. Molecular Biology and Evolution, 1987, 4(1):74-83.

    [6]Ayala F J, Rzhetsky A, Ayala F J. Origin of the metazoan phyla: molecular clocks confirm paleontological estimates[J]. Proc Natl Acad Sci, 1998, 95:606-611.

    [7]賀福初. 分子減速進(jìn)化的普遍性研究[J]. 科學(xué)通報(bào),1996, 41(24): 2264-2268.

    [8]羅 靜. 分子鐘及其存在的問題[J]. 人類學(xué)學(xué)報(bào),2000, 19(2): 151-159.

    [9]Adachi A P, Cao Y, Hasegawa M. Tempo and mode of mitochondrial DNA evolution in vertebrates at the amino acid sequence level: Rapid evolution in warm-blooded vertebrates[J]. Journal of molecular evolution, 1993, 36: 270-281.

    [10]徐宏發(fā). 分子系統(tǒng)學(xué)研究進(jìn)展[J]. 生態(tài)學(xué)雜志,2001, 20(3): 41-46.

    [11]李可群. 分子絕對進(jìn)化速率計(jì)算公式的推導(dǎo)及應(yīng)用方法[J]. 高師理科學(xué)刊, 2015, 35(1):19-21.

    [12]Kimura M. The neutral theory of molecular evolution[M]. Cambridge University Press, Cambridge, 1983.

    [13]宮田隆著, 劉績生譯. 什么是分子鐘[J]. 世界科學(xué). 1987,11:22-23.

    [14]張英培. 分子分類的若干問題[J]. 動(dòng)物學(xué)研究, 1994, 15(1): 1-10.

    [15]呂寶忠. 躍入分子水平的群體研究學(xué)和進(jìn)化遺傳學(xué)[J]. 生物學(xué)雜志, 1988, 24(4): 5-8,17.

    Quantitative analysis of relationship between absolute evolutionary rates and taxa divergence times

    LI Ke-qun

    (Department of Chemistry, Tongji University, Shanghai 200092, China)

    molecular evolution; absolute evolutionary rate; molecular o′clock; quantitative analysis; phylogenetic analysis

    2014-08-18;

    2014-09-24

    李可群,講師,研究方向?yàn)榉肿舆z傳進(jìn)化。

    Q

    A

    2095-1736(2015)02-0070-06

    doi∶10.3969/j.issn.2095-1736.2015.02.070

    猜你喜歡
    核苷酸定量位點(diǎn)
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進(jìn)展
    徐長風(fēng):核苷酸類似物的副作用
    肝博士(2022年3期)2022-06-30 02:48:28
    鎳基單晶高溫合金多組元置換的第一性原理研究
    上海金屬(2021年6期)2021-12-02 10:47:20
    CLOCK基因rs4580704多態(tài)性位點(diǎn)與2型糖尿病和睡眠質(zhì)量的相關(guān)性
    顯微定量法鑒別林下山參和園參
    Acknowledgment to reviewers—November 2018 to September 2019
    二項(xiàng)式通項(xiàng)公式在遺傳學(xué)計(jì)算中的運(yùn)用*
    當(dāng)歸和歐當(dāng)歸的定性與定量鑒別
    中成藥(2018年12期)2018-12-29 12:25:44
    10 種中藥制劑中柴胡的定量測定
    中成藥(2017年6期)2017-06-13 07:30:35
    慢性HBV感染不同狀態(tài)下HBsAg定量的臨床意義
    国产亚洲av高清不卡| 日韩欧美在线二视频| 色尼玛亚洲综合影院| 久久久久国内视频| 午夜福利免费观看在线| av国产免费在线观看| 老汉色∧v一级毛片| 天天一区二区日本电影三级| 男女那种视频在线观看| 99久久国产精品久久久| 最近视频中文字幕2019在线8| av中文乱码字幕在线| 女同久久另类99精品国产91| 日韩欧美在线乱码| 色综合欧美亚洲国产小说| 好男人在线观看高清免费视频| 黄频高清免费视频| 97超视频在线观看视频| 黄色 视频免费看| 不卡av一区二区三区| 亚洲欧美日韩高清在线视频| 欧美3d第一页| 真人一进一出gif抽搐免费| 欧美成人性av电影在线观看| 哪里可以看免费的av片| 一a级毛片在线观看| 好男人电影高清在线观看| 国产精品一区二区三区四区免费观看 | 天堂影院成人在线观看| 日本撒尿小便嘘嘘汇集6| 毛片女人毛片| 一个人观看的视频www高清免费观看 | 美女黄网站色视频| 巨乳人妻的诱惑在线观看| 久久久久久九九精品二区国产| 白带黄色成豆腐渣| 午夜日韩欧美国产| 最新美女视频免费是黄的| 成人鲁丝片一二三区免费| 好男人在线观看高清免费视频| 亚洲国产欧美一区二区综合| 精品久久久久久久毛片微露脸| 久久中文看片网| 五月伊人婷婷丁香| 成人国产一区最新在线观看| 国产视频一区二区在线看| 国产免费男女视频| 人妻丰满熟妇av一区二区三区| 成熟少妇高潮喷水视频| 久久久久性生活片| 国产亚洲精品久久久com| 国产成人精品无人区| 亚洲无线在线观看| 在线免费观看的www视频| 啦啦啦免费观看视频1| 伊人久久大香线蕉亚洲五| 黑人巨大精品欧美一区二区mp4| 亚洲va日本ⅴa欧美va伊人久久| 亚洲在线观看片| 国产av一区在线观看免费| 狂野欧美激情性xxxx| 国产麻豆成人av免费视频| a在线观看视频网站| av视频在线观看入口| 极品教师在线免费播放| 欧美一区二区国产精品久久精品| 五月玫瑰六月丁香| 99久久精品热视频| 美女大奶头视频| 神马国产精品三级电影在线观看| 国产一级毛片七仙女欲春2| 久久国产精品影院| 亚洲最大成人中文| 国产成+人综合+亚洲专区| 国产1区2区3区精品| 欧美一级a爱片免费观看看| 伦理电影免费视频| 色av中文字幕| 女人被狂操c到高潮| 国产日本99.免费观看| 别揉我奶头~嗯~啊~动态视频| 18禁裸乳无遮挡免费网站照片| 成人午夜高清在线视频| 久久中文字幕一级| 久久久国产精品麻豆| 一二三四社区在线视频社区8| 久久天堂一区二区三区四区| 国产精品香港三级国产av潘金莲| 成人国产一区最新在线观看| 一区二区三区高清视频在线| 一区二区三区激情视频| 欧美日韩瑟瑟在线播放| 亚洲成av人片在线播放无| 成人精品一区二区免费| 看免费av毛片| 啦啦啦观看免费观看视频高清| 最新美女视频免费是黄的| 中文字幕最新亚洲高清| 99热这里只有精品一区 | 国产精品亚洲美女久久久| 免费观看精品视频网站| 亚洲第一电影网av| 熟女少妇亚洲综合色aaa.| 97超级碰碰碰精品色视频在线观看| 欧美成狂野欧美在线观看| 中国美女看黄片| а√天堂www在线а√下载| 亚洲国产欧美人成| 九色成人免费人妻av| 1024香蕉在线观看| 国产精品av视频在线免费观看| 久久精品国产亚洲av香蕉五月| 国产精品自产拍在线观看55亚洲| 成人性生交大片免费视频hd| 怎么达到女性高潮| 日本成人三级电影网站| 巨乳人妻的诱惑在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 午夜成年电影在线免费观看| 欧美黄色淫秽网站| 欧美日韩福利视频一区二区| 男插女下体视频免费在线播放| 岛国在线免费视频观看| 国产午夜精品久久久久久| 一进一出好大好爽视频| 亚洲自拍偷在线| 99久久久亚洲精品蜜臀av| 久久久久免费精品人妻一区二区| 国产精品久久电影中文字幕| 精品久久蜜臀av无| 国内久久婷婷六月综合欲色啪| 亚洲av中文字字幕乱码综合| 老司机深夜福利视频在线观看| 亚洲熟女毛片儿| 999精品在线视频| 国产高清视频在线播放一区| 久久久国产欧美日韩av| 色综合亚洲欧美另类图片| 淫秽高清视频在线观看| 男女视频在线观看网站免费| 成年人黄色毛片网站| 国产麻豆成人av免费视频| 久久久久性生活片| 国产人伦9x9x在线观看| 久久精品国产综合久久久| 岛国在线观看网站| 国产高清视频在线播放一区| 亚洲精品一区av在线观看| 国产成人系列免费观看| 国产一区在线观看成人免费| 久久久成人免费电影| 久久精品影院6| а√天堂www在线а√下载| 在线观看午夜福利视频| 无人区码免费观看不卡| 亚洲天堂国产精品一区在线| 欧美高清成人免费视频www| 亚洲熟女毛片儿| 久久国产乱子伦精品免费另类| 五月伊人婷婷丁香| www.999成人在线观看| 欧美一区二区国产精品久久精品| 色播亚洲综合网| 国产精品 国内视频| av视频在线观看入口| 日日干狠狠操夜夜爽| 叶爱在线成人免费视频播放| 噜噜噜噜噜久久久久久91| 99在线人妻在线中文字幕| 丰满的人妻完整版| 久久久久亚洲av毛片大全| 午夜免费成人在线视频| 日本免费一区二区三区高清不卡| 黄色成人免费大全| 久久久久免费精品人妻一区二区| 亚洲人成网站高清观看| 欧美极品一区二区三区四区| 国产单亲对白刺激| 热99re8久久精品国产| 两人在一起打扑克的视频| 老司机午夜十八禁免费视频| 热99re8久久精品国产| 女警被强在线播放| 一本一本综合久久| 亚洲aⅴ乱码一区二区在线播放| or卡值多少钱| 黄色视频,在线免费观看| 国产一区二区三区在线臀色熟女| 男女床上黄色一级片免费看| 午夜a级毛片| 91av网一区二区| 久久久水蜜桃国产精品网| 国产精品久久久久久精品电影| 在线视频色国产色| 精品人妻1区二区| 午夜福利高清视频| 国内少妇人妻偷人精品xxx网站 | 午夜福利免费观看在线| 欧美日韩中文字幕国产精品一区二区三区| 国产av不卡久久| 18禁国产床啪视频网站| 国产91精品成人一区二区三区| 国产91精品成人一区二区三区| 宅男免费午夜| 久久中文看片网| 欧美日韩中文字幕国产精品一区二区三区| 露出奶头的视频| 制服丝袜大香蕉在线| 精品熟女少妇八av免费久了| 成人鲁丝片一二三区免费| 久久久久国产一级毛片高清牌| 精品久久蜜臀av无| 欧美3d第一页| 久久久久精品国产欧美久久久| 可以在线观看毛片的网站| 久久久久免费精品人妻一区二区| h日本视频在线播放| 国产高潮美女av| 亚洲成人中文字幕在线播放| 国产精品,欧美在线| 国产成人福利小说| 舔av片在线| 搡老妇女老女人老熟妇| 亚洲av熟女| 一本久久中文字幕| 国产成人一区二区三区免费视频网站| 中出人妻视频一区二区| 亚洲天堂国产精品一区在线| 欧美最黄视频在线播放免费| 亚洲国产欧美网| 亚洲欧美精品综合一区二区三区| 国产精品综合久久久久久久免费| 一级a爱片免费观看的视频| 亚洲美女黄片视频| 中文亚洲av片在线观看爽| 亚洲电影在线观看av| 亚洲色图av天堂| 国内精品美女久久久久久| 女人高潮潮喷娇喘18禁视频| 亚洲片人在线观看| 天堂√8在线中文| 我要搜黄色片| a在线观看视频网站| 亚洲成a人片在线一区二区| www.精华液| 亚洲 欧美一区二区三区| 99国产精品一区二区三区| 我的老师免费观看完整版| 天天躁狠狠躁夜夜躁狠狠躁| 久久天堂一区二区三区四区| 久久香蕉国产精品| 母亲3免费完整高清在线观看| 在线观看日韩欧美| 欧美极品一区二区三区四区| 日日摸夜夜添夜夜添小说| 亚洲最大成人中文| 1000部很黄的大片| 中文字幕人成人乱码亚洲影| 亚洲av成人av| 51午夜福利影视在线观看| 给我免费播放毛片高清在线观看| xxxwww97欧美| 真实男女啪啪啪动态图| 欧美乱色亚洲激情| 全区人妻精品视频| 久久久国产欧美日韩av| 中文字幕人妻丝袜一区二区| 观看美女的网站| 亚洲成av人片免费观看| 国产v大片淫在线免费观看| 男女那种视频在线观看| 亚洲自偷自拍图片 自拍| 欧美日韩黄片免| 国产又色又爽无遮挡免费看| 国产三级黄色录像| 亚洲精华国产精华精| 首页视频小说图片口味搜索| 黑人欧美特级aaaaaa片| 日韩国内少妇激情av| 国产一区二区三区视频了| 欧美在线黄色| 久久久久久国产a免费观看| 天天躁日日操中文字幕| 欧美3d第一页| 亚洲 欧美一区二区三区| 婷婷亚洲欧美| 亚洲精品一区av在线观看| 亚洲18禁久久av| 啪啪无遮挡十八禁网站| x7x7x7水蜜桃| 人妻丰满熟妇av一区二区三区| h日本视频在线播放| 午夜免费成人在线视频| 一二三四社区在线视频社区8| 中亚洲国语对白在线视频| 两人在一起打扑克的视频| 国产欧美日韩精品一区二区| 婷婷精品国产亚洲av| 日韩欧美 国产精品| 在线免费观看的www视频| 99久久精品一区二区三区| 亚洲专区国产一区二区| 搡老岳熟女国产| 噜噜噜噜噜久久久久久91| 欧美色视频一区免费| 欧美大码av| 成人国产综合亚洲| 欧美日韩黄片免| 国产精品一及| 日韩av在线大香蕉| 国产单亲对白刺激| 99久久99久久久精品蜜桃| 999精品在线视频| 在线免费观看的www视频| 日本 欧美在线| 淫妇啪啪啪对白视频| av天堂中文字幕网| 国产一区在线观看成人免费| 叶爱在线成人免费视频播放| 后天国语完整版免费观看| 午夜福利成人在线免费观看| 欧美绝顶高潮抽搐喷水| av中文乱码字幕在线| 欧美性猛交黑人性爽| 精品熟女少妇八av免费久了| 人妻丰满熟妇av一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲美女视频黄频| 制服丝袜大香蕉在线| aaaaa片日本免费| 亚洲av中文字字幕乱码综合| 三级国产精品欧美在线观看 | bbb黄色大片| 亚洲在线自拍视频| 国产视频一区二区在线看| 巨乳人妻的诱惑在线观看| 亚洲九九香蕉| 99热6这里只有精品| 国产成人av激情在线播放| av福利片在线观看| av欧美777| 一夜夜www| 成人鲁丝片一二三区免费| 中文资源天堂在线| 国产精品久久久久久亚洲av鲁大| 亚洲精品国产精品久久久不卡| 一本一本综合久久| 久久国产精品人妻蜜桃| 变态另类成人亚洲欧美熟女| 18禁黄网站禁片免费观看直播| or卡值多少钱| 国产精品美女特级片免费视频播放器 | 亚洲va日本ⅴa欧美va伊人久久| 午夜免费成人在线视频| 视频区欧美日本亚洲| 国产成人aa在线观看| 国产欧美日韩精品一区二区| 我的老师免费观看完整版| 天堂√8在线中文| 国产欧美日韩一区二区精品| 久久国产精品人妻蜜桃| 九九在线视频观看精品| 久久中文看片网| 日本一二三区视频观看| 九色国产91popny在线| 成人特级黄色片久久久久久久| 香蕉av资源在线| 又大又爽又粗| 制服人妻中文乱码| 欧美三级亚洲精品| 久久久久久久久久黄片| 免费观看的影片在线观看| 伦理电影免费视频| 在线观看免费视频日本深夜| 91麻豆精品激情在线观看国产| 99riav亚洲国产免费| 高清在线国产一区| 99视频精品全部免费 在线 | 成人18禁在线播放| av天堂中文字幕网| 18禁国产床啪视频网站| 日本一二三区视频观看| 精品一区二区三区av网在线观看| 亚洲中文字幕一区二区三区有码在线看 | 狠狠狠狠99中文字幕| 成人三级做爰电影| 色老头精品视频在线观看| 99re在线观看精品视频| 男女视频在线观看网站免费| 色尼玛亚洲综合影院| 精品不卡国产一区二区三区| 88av欧美| 精品熟女少妇八av免费久了| 亚洲成av人片在线播放无| 亚洲国产精品999在线| 欧美黑人巨大hd| 国产高清三级在线| 国产爱豆传媒在线观看| 九九热线精品视视频播放| 女人被狂操c到高潮| 网址你懂的国产日韩在线| 丰满人妻一区二区三区视频av | 琪琪午夜伦伦电影理论片6080| av片东京热男人的天堂| 日韩精品中文字幕看吧| 国产精品久久久人人做人人爽| 久久九九热精品免费| 一级毛片高清免费大全| 亚洲国产精品999在线| 最近最新免费中文字幕在线| 黑人欧美特级aaaaaa片| 精品电影一区二区在线| 观看美女的网站| 国产av一区在线观看免费| 人妻夜夜爽99麻豆av| 听说在线观看完整版免费高清| 99久久综合精品五月天人人| 精品不卡国产一区二区三区| 变态另类丝袜制服| 欧美丝袜亚洲另类 | 免费在线观看影片大全网站| 最近视频中文字幕2019在线8| 精品日产1卡2卡| 久久精品91蜜桃| 噜噜噜噜噜久久久久久91| 夜夜躁狠狠躁天天躁| 国内精品美女久久久久久| 听说在线观看完整版免费高清| 91麻豆av在线| 国产精品亚洲一级av第二区| 99在线人妻在线中文字幕| 国产69精品久久久久777片 | 日韩欧美免费精品| 国产激情久久老熟女| 老汉色∧v一级毛片| 偷拍熟女少妇极品色| 国内久久婷婷六月综合欲色啪| 1000部很黄的大片| 黄频高清免费视频| 最近最新中文字幕大全电影3| 18禁观看日本| 国产成人精品久久二区二区91| 淫妇啪啪啪对白视频| 亚洲乱码一区二区免费版| 这个男人来自地球电影免费观看| 少妇裸体淫交视频免费看高清| 熟女少妇亚洲综合色aaa.| 亚洲精品在线美女| 中文字幕熟女人妻在线| 欧美中文综合在线视频| 老汉色∧v一级毛片| 欧美日韩乱码在线| 欧美大码av| 波多野结衣高清无吗| 国产高清有码在线观看视频| 黑人欧美特级aaaaaa片| av天堂在线播放| 亚洲精华国产精华精| 成人欧美大片| 国产真人三级小视频在线观看| 久久久久九九精品影院| 日韩欧美在线乱码| av欧美777| 日韩有码中文字幕| 色老头精品视频在线观看| 天堂动漫精品| 欧美色欧美亚洲另类二区| 日本成人三级电影网站| 久久久久久大精品| svipshipincom国产片| 男人舔奶头视频| 最近最新免费中文字幕在线| 国产1区2区3区精品| 在线观看一区二区三区| 亚洲人成电影免费在线| 欧美激情久久久久久爽电影| 亚洲av成人av| 老司机福利观看| 国产乱人视频| 最近最新免费中文字幕在线| 国产亚洲精品综合一区在线观看| 香蕉av资源在线| 亚洲一区高清亚洲精品| 麻豆久久精品国产亚洲av| av视频在线观看入口| 丁香六月欧美| 亚洲美女视频黄频| 两人在一起打扑克的视频| 精品一区二区三区视频在线观看免费| 18禁美女被吸乳视频| 丁香六月欧美| 1024香蕉在线观看| 成年女人毛片免费观看观看9| 一进一出抽搐gif免费好疼| 欧美黑人巨大hd| 91九色精品人成在线观看| 五月玫瑰六月丁香| 亚洲欧美精品综合久久99| 国产乱人视频| 男人的好看免费观看在线视频| 日本熟妇午夜| 亚洲av日韩精品久久久久久密| 制服人妻中文乱码| 久久国产乱子伦精品免费另类| 午夜精品在线福利| 国产毛片a区久久久久| 噜噜噜噜噜久久久久久91| 久久久精品大字幕| 久久婷婷人人爽人人干人人爱| 最新在线观看一区二区三区| 国内精品美女久久久久久| 男女视频在线观看网站免费| 网址你懂的国产日韩在线| 欧美在线一区亚洲| 757午夜福利合集在线观看| 岛国在线免费视频观看| 久久久久国产精品人妻aⅴ院| 久久中文字幕人妻熟女| 中文资源天堂在线| 午夜福利免费观看在线| 真人做人爱边吃奶动态| 国产真人三级小视频在线观看| 性色av乱码一区二区三区2| 亚洲欧美日韩东京热| 欧美在线一区亚洲| h日本视频在线播放| 久久亚洲真实| 欧美日韩综合久久久久久 | 欧美+亚洲+日韩+国产| 亚洲成av人片免费观看| 欧美三级亚洲精品| 成人性生交大片免费视频hd| 日韩精品青青久久久久久| 久久这里只有精品19| 一级a爱片免费观看的视频| 非洲黑人性xxxx精品又粗又长| 亚洲乱码一区二区免费版| 99热这里只有精品一区 | 男人和女人高潮做爰伦理| 麻豆久久精品国产亚洲av| 午夜激情福利司机影院| 亚洲国产精品sss在线观看| 欧美另类亚洲清纯唯美| 亚洲,欧美精品.| 国产精品野战在线观看| 嫩草影院精品99| 窝窝影院91人妻| 舔av片在线| 超碰成人久久| 最近最新中文字幕大全免费视频| 无限看片的www在线观看| 99国产精品一区二区蜜桃av| h日本视频在线播放| 最近最新免费中文字幕在线| 亚洲电影在线观看av| 精品一区二区三区四区五区乱码| 美女免费视频网站| www.熟女人妻精品国产| 亚洲精品美女久久久久99蜜臀| 人人妻人人看人人澡| 亚洲欧美一区二区三区黑人| 亚洲av片天天在线观看| 亚洲av熟女| 精品一区二区三区视频在线观看免费| 成人性生交大片免费视频hd| 国产三级在线视频| 日韩欧美在线二视频| 国产精品,欧美在线| 中文字幕高清在线视频| 99国产精品一区二区蜜桃av| 三级男女做爰猛烈吃奶摸视频| 嫩草影院入口| 网址你懂的国产日韩在线| 欧美成人一区二区免费高清观看 | 男人的好看免费观看在线视频| 国内少妇人妻偷人精品xxx网站 | 国内揄拍国产精品人妻在线| 91老司机精品| 90打野战视频偷拍视频| 亚洲国产精品合色在线| 香蕉丝袜av| 黄色丝袜av网址大全| 国内揄拍国产精品人妻在线| 久久热在线av| 国产欧美日韩一区二区三| 97超级碰碰碰精品色视频在线观看| 久久这里只有精品19| 国产精品久久久久久亚洲av鲁大| 国内精品久久久久久久电影| 手机成人av网站| 成年人黄色毛片网站| 亚洲av成人不卡在线观看播放网| 最近在线观看免费完整版| 亚洲真实伦在线观看| x7x7x7水蜜桃| 久久国产乱子伦精品免费另类| 久久伊人香网站| 少妇人妻一区二区三区视频| 99热这里只有精品一区 | 热99在线观看视频| 人妻久久中文字幕网| 国产成人系列免费观看| 色视频www国产| 午夜福利在线在线| 91麻豆av在线| 1024香蕉在线观看| 欧美中文综合在线视频| АⅤ资源中文在线天堂| 国产午夜精品论理片| 日日干狠狠操夜夜爽| 国产午夜福利久久久久久| av片东京热男人的天堂| 少妇人妻一区二区三区视频| 久久人妻av系列| 久久99热这里只有精品18|