紫紅笛鯛線粒體全序列與笛鯛屬魚(yú)類分化年代的貝葉斯分析

廖 杰, 王中鐸, 郭昱嵩, 劉楚吾

(廣東海洋大學(xué) 水產(chǎn)學(xué)院 南海水產(chǎn)經(jīng)濟(jì)動(dòng)物增養(yǎng)殖廣東普通高校重點(diǎn)實(shí)驗(yàn)室, 廣東 湛江 524025)

紫紅笛鯛線粒體全序列與笛鯛屬魚(yú)類分化年代的貝葉斯分析

廖 杰, 王中鐸, 郭昱嵩, 劉楚吾

(廣東海洋大學(xué) 水產(chǎn)學(xué)院 南海水產(chǎn)經(jīng)濟(jì)動(dòng)物增養(yǎng)殖廣東普通高校重點(diǎn)實(shí)驗(yàn)室, 廣東 湛江 524025)

利用常規(guī) PCR未純化產(chǎn)物直接測(cè)序的方法, 獲得紫紅笛鯛(Lutjanus argentimaculatus)全長(zhǎng)為16 543 bp線粒體DNA全序列, GenBank登錄號(hào)為JN182927。結(jié)合GenBank中的已有數(shù)據(jù)比較分析表明: 紫紅笛鯛與笛鯛屬(Lutjanus)線粒體基因組成基本一致, 符合脊椎動(dòng)物的特征, 各基因間進(jìn)化速率差異顯著; 13個(gè)蛋白編碼基因拼接序列貝葉斯建樹(shù)的后驗(yàn)率顯著高于單基因或全序列, 更適用于笛鯛屬的貝葉斯法系統(tǒng)分化分析和年代推斷; 13個(gè)蛋白編碼基因及其拼接序列進(jìn)行的系統(tǒng)進(jìn)化分析均支持紫紅笛鯛與勒氏笛鯛(Lutjanus russellii)及海雞母笛鯛(Lutjanus rivulatus)親緣關(guān)系較近; 基于13個(gè)蛋白編碼基因拼接序列的貝葉斯松散分子鐘推算出笛鯛屬魚(yú)類兩個(gè)主要支系大約在 61.5Ma前古新世的達(dá)寧階與蒙蒂階交界時(shí)期發(fā)生分化。

紫紅笛鯛(Lutjanus argentimaculatus); 線粒體全序列; 貝葉斯分析; 松散分子鐘; 笛鯛屬魚(yú)類

紫紅笛鯛(Lutjanus argentimaculatus)隸屬硬骨魚(yú)綱(Osteichthyes), 鱸形目(Perciformes), 笛鯛科(Lutjanidae), 笛鯛屬(Lutjanus)。主要分布于印度洋和太平洋中部及西部, 包括中國(guó)的南海和東海南部, 是一種優(yōu)質(zhì)的食用魚(yú)類。有關(guān)紫紅笛鯛的研究很多, 涉及人工養(yǎng)殖、種群遺傳多樣性以及它在笛鯛屬中的關(guān)系等多個(gè)方面。邢玉娜等[1]、區(qū)又君等[2]分別研究了紫紅笛鯛育苗和幼魚(yú)養(yǎng)殖密度與攝食關(guān)系; 張俊彬等[3]研究了中國(guó)南海部分地區(qū)紫紅笛鯛養(yǎng)殖和野生群體的遺傳多樣性。Russell等[4]綜述了澳大利亞紫紅笛鯛生態(tài)特征、群體遺傳結(jié)構(gòu)以及養(yǎng)殖技術(shù), 并記錄紫紅笛鯛存在獨(dú)特的生活史:幼魚(yú)和早期成魚(yú)都可以生活在淡水溪流中。Guo等[5]基于線粒體基因(COX2和Cytb)的系統(tǒng)進(jìn)化表明紫紅笛鯛與紅鰭笛鯛(Lutjanus erythopterus)關(guān)系密切;王中鐸等[6]綜合利用3個(gè)線粒體基因(COX1、COX2和Cytb)和2個(gè)核基因(RAG1和RAG2)構(gòu)建了笛鯛屬魚(yú)類的系統(tǒng)進(jìn)化樹(shù), 揭示出紫紅笛鯛在進(jìn)化上最早與大體型、紅體色、近岸底棲的紅鰭笛鯛和千年笛鯛(Lutjanus sebae)一支分化, 獨(dú)立于其他笛鯛屬魚(yú)類, 自成一支。

近年來(lái)發(fā)展起來(lái)的貝葉斯模型參數(shù)可直接進(jìn)行量化,并以后驗(yàn)概率來(lái)表示各分支的可信度而不需用自引導(dǎo)法(Bootstrap)進(jìn)行檢驗(yàn), 而且在計(jì)算速度和可靠性上明顯優(yōu)于其他算法, 十分利于線粒體全序列等較大數(shù)據(jù)集的分析[7-8]。 另外, 貝葉斯松散分子鐘利用化石標(biāo)定點(diǎn)計(jì)算分化時(shí)間, 推測(cè)出實(shí)驗(yàn)對(duì)象出現(xiàn)或發(fā)生年代已經(jīng)引起魚(yú)類進(jìn)化研究者的廣泛重視[9-11]。隨著測(cè)序技術(shù)的進(jìn)步使線粒體全序列的獲得變得越來(lái)越容易, 用線粒體全序列作為遺傳進(jìn)化分析的材料, 將會(huì)獲得更加全面的信息進(jìn)而得到更合理的結(jié)論, 這一技術(shù)在魚(yú)類進(jìn)化研究中的應(yīng)用越來(lái)越廣泛[12-14]。一直以來(lái)對(duì)笛鯛屬魚(yú)類的分類學(xué)地位和形態(tài)分類學(xué)標(biāo)準(zhǔn)都存在爭(zhēng)議[6,11,13]。Guo等[5-6]利用線粒體基因(COX2和Cytb)和核基因重建了笛鯛屬的系統(tǒng)進(jìn)化歷程, 并且先后獲取了多種笛鯛屬魚(yú)類的線粒體DNA(mtDNA)全序列。

本研究通過(guò)測(cè)定紫紅笛鯛全序列進(jìn)而確定該物種在笛鯛屬的分類地位, 并通過(guò)貝葉斯分析推斷笛鯛屬各支系的分化年代同時(shí)分析比較線粒體全序列及其各編碼區(qū)段貝葉斯建樹(shù)的優(yōu)劣性, 為揭示魚(yú)類演化研究積累基礎(chǔ)資料。

1 材料與方法

1.1 實(shí)驗(yàn)材料和總DNA的提取

紫紅笛鯛0.8 kg, 性腺發(fā)育未成熟, 購(gòu)自湛江市霞山區(qū)水產(chǎn)品交易市場(chǎng)。采用形態(tài)與DNA 條碼結(jié)合的方法鑒定實(shí)驗(yàn)對(duì)象[14], 參照盧圣棟等[15]的方法從肌肉或鰭條中提取樣本DNA, 4℃保存?zhèn)溆谩?/p>

1.2 LA-PCR純化線粒體DNA

引物參考王中鐸等[10](表 1), 均由上海 Sangon生物技術(shù)公司合成。采用長(zhǎng)距 PCR(LA-PC R)引物L(fēng)A1和LA2擴(kuò)增得到覆蓋整個(gè)線粒體DNA序列的大片段, PCR反應(yīng)體系為: 總體積25 μL, 含1×PCR LA Buffer, dNTPs 15 nmol, 引物各 20 pmol, LATaqDNA聚合酶(TaKaRa) 2.5U, 總DNA 約20 ng。擴(kuò)增條件: 94 ℃預(yù)變性1 min; 98℃變性10 s, 68℃退火16 min, 循環(huán)30次; 最后72℃延伸10 min。1%的瓊脂糖電泳對(duì)擴(kuò)增產(chǎn)物進(jìn)行檢測(cè)。

1.3 常規(guī)PCR擴(kuò)增獲取線粒體DNA片段

以稀釋20倍的LA-PCR 產(chǎn)物為模板, 使用表1所列常規(guī)PCR引物擴(kuò)增獲取長(zhǎng)度600～1 000 bp左右的相互重疊的小片段。反應(yīng)體系為: 總體積50 μL, 含1×PCR Buffer, dNTPs 10 nmol, 引物各20 pmol, rTaqDNA聚合酶(TaKaRa) 2.5 U, DNA約20 ng。PCR 反應(yīng)的循環(huán)參數(shù)為: 94 ℃預(yù)變性3 min, 94 ℃變性42 s, 50～55 ℃退火30 s, 72 ℃延伸55 s, 共30個(gè)循環(huán), 最終72 ℃延伸10 min, 用1.0%瓊脂糖凝膠電泳檢測(cè)擴(kuò)增產(chǎn)物, 凝膠成像儀拍照。PCR擴(kuò)增產(chǎn)物委托廣州華大基因公司進(jìn)行序列雙向測(cè)定。

1.4 全序列拼接與分析

用Invitrogen軟件包的ContigExpress進(jìn)行序列拼接, 利用Editseq程序(http://www.dnastar. com/)查找序列進(jìn)行開(kāi)放閱讀框, 獲得序列互補(bǔ)鏈以及統(tǒng)計(jì)各密碼子所占比例; 利用 MEGA 4.0 軟件完成蛋白編碼序列的翻譯和氨基酸編碼情況的調(diào)查[16]; 在線查找tRNA (http://lowelab.ucsc.edu/tRNAscan-SE)[17];用 Sequin 軟件(http://www.ncbi.nlm.nih.gov/Sequin/ index.html)遞交序列至 GenBank。從 GenBank下載的其余7種笛鯛的線粒體全序列(6種為本研究室測(cè)定)。分別利用ATPase6、Cytb、COX1、COX2、ND1、ND2、ND4、ND5等基因以及D-loop區(qū)進(jìn)行遺傳距離分析。

1.5 基于線粒體全序列和編碼基因及其拼接序列的貝葉斯分析

參考Drummond等[18]的方法利用Beast 1.6.1軟件包(http://beast.bio.ed.ac.uk/)構(gòu)建貝葉斯樹(shù), 并利用貝葉斯松散分子鐘(uncorrected relaxed lognormal clock)推斷分化年代[19]。軟件運(yùn)行設(shè)置如下: (GTR+G+I+relaxed, MCMC=1000000, burn in=1000)利用 FigTree v1.3.1(http://tree.bio.ed.ac.uk/)生成樹(shù)圖。所選的外群物種為斑馬魚(yú)(Danio rerio)和青鳉(Oryzias latipes), 化 石 標(biāo) 定 點(diǎn) 則 為 鲀 科(Tetraodontidae)魚(yú)類紅鰭東方豚(Takifugu rubripes)和凹鼻鲀(Tetraodon nigroviridis), 它們的分化上下限為56.0 ～32.25 Ma[20]。

1.6 不同基因貝葉斯后驗(yàn)概率分析

分別建立全序列、13個(gè)蛋白編碼基因及其拼接序列的貝葉斯分析樹(shù), 以笛鯛主分支, 紫紅笛鯛所在支以及所有分支平均貝葉斯后驗(yàn)概率最為判斷各基因和拼接序列貝葉斯建樹(shù)的優(yōu)劣性。

2 實(shí)驗(yàn)結(jié)果

2.1 片段擴(kuò)增、序列拼接及分析

LA-PCR獲得的長(zhǎng)距產(chǎn)物(圖1)為模板進(jìn)行常規(guī)PCR(圖2), 擴(kuò)增產(chǎn)物直接測(cè)序后拼接獲得16 543 bp的紫紅笛鯛線粒體全序列。序列中(A、T、C、G)4種堿基分別占 28.25%, 24.68%, 30.95%, 16.12%, A+T＞G+C。用Sequin上傳至Gene Bank, 登錄號(hào)為JN182927。紫紅笛鯛mtDNA基因組包括13個(gè)蛋白編碼基因, 22個(gè)tRNA基因, 2個(gè)rRNA基因和1個(gè)控制區(qū)。基因分布呈不均一性, 其中輕鏈(L-strand)僅編碼8個(gè)tRNA基因(Pro、Glu、Ser、Tyr、Cys、Asn、Ala、Gln)及ND6基因, 其余基因皆由重鏈(H-strand)編碼。

2.2 與已有笛鯛線粒體全序列的比較分析

分析得到紫紅笛鯛與7種笛鯛的ATPase6、Cyt b、COX 1、COX 2、ND 1、ND 2、ND 4、ND 5、D-loop和線粒體全序列的遺傳距離(表 2)。從表中可見(jiàn), 單基因的物種間遺傳距離平均值為線粒體全序列間的1.4～1.9倍。基于基因片段所得的K2P遺傳距離平均值比較可以看出: 紫紅笛鯛與馬拉巴笛鯛(Lutjanus malabarius)、紅鰭笛鯛和千年笛鯛的距離較遠(yuǎn)(≥0.239), 而與其他笛鯛相對(duì)較近(≤0.195); 其中在

表1 線粒體DNA基因組序列測(cè)定引物Tab. 1 PCR primers for mtDNA complete genome

ATPase6間的遺傳距離比較中紫紅笛鯛與其他 4種笛鯛的遺傳距離值比其他基因片段的大, 而線粒體全序列間遺傳距離比較的結(jié)果與 12個(gè)基因片段的相似,均為馬拉巴笛鯛、紅鰭笛鯛和千年笛鯛的遺傳距離較遠(yuǎn)。所有基因中海雞母笛鯛(L. rivulatus)的距離都最小,它相對(duì)其他笛鯛與紫紅笛鯛的關(guān)系更為密切。

圖1 LA-PCR電泳檢測(cè)圖Fig. 1 Electrophoresis of LA-PCR

圖2 部分常規(guī)PCR電泳檢測(cè)圖Fig. 2 Electrophoresis of common PCR

2.3 基于貝葉斯松散分子鐘的笛鯛分化年代推斷分析

分別建立 8種笛鯛線粒體全序列、線粒體各個(gè)編碼基因以及線粒體COX1基因的貝葉斯分析樹(shù),得到多個(gè)貝葉斯推斷圖(圖3～圖5)。3圖中標(biāo)定化石鲀科(Tetraodontidae)種都落在化石記錄范圍的中心?；诰€粒體全序列的推斷結(jié)果顯示笛鯛屬的分化年代處于40.86Ma, 而基于13個(gè)蛋白編碼基因拼接序列的分析則顯示分化發(fā)生在 61.49Ma, 兩者推斷的主分支年代差異很大。圖 3中個(gè)小分支發(fā)生都位于32.7～8.2Ma這個(gè)相對(duì)較短的時(shí)間, 而圖 4中各小分支的出現(xiàn)年代跨度則小得多(57.8～48.7Ma)。兩個(gè)系統(tǒng)進(jìn)化樹(shù)中8種笛鯛都可以明確的劃分成兩個(gè)支系。進(jìn)化樹(shù)和遺傳距離各自的比較結(jié)果相吻合。遺傳距離較遠(yuǎn)的紅鰭笛鯛, 馬拉巴笛鯛和千年笛鯛聚為一支, 紫紅笛鯛所在一支的另五種笛鯛聚為另一支。

表2 K2P紫紅笛鯛與7種笛鯛各基因區(qū)段遺傳距離比較Tab. 2 K2P Genetic distances between Lutjanus argentimaculatus and other snappers

在 3個(gè)圖中紫紅笛鯛所處的位置存在差異, 基于線粒體全序列的貝葉斯進(jìn)化樹(shù)笛鯛分支不夠細(xì)致,而另外兩個(gè)圖中各個(gè)小支明顯分開(kāi)。圖3和圖5中紫紅笛鯛和海雞母笛鯛(Lutjanus rivulatus)關(guān)系較緊密, 圖4中紫紅笛鯛和勒氏笛鯛(Lutjanus russellii)先聚為一支再和海雞母笛鯛并列。

線粒體DNA全序列、13個(gè)基因拼接序列以及單獨(dú)的基因序列建立貝葉斯分析樹(shù), 笛鯛主分支、紫紅笛鯛分支以及平均后驗(yàn)概率統(tǒng)計(jì)如表 3所示?；趍tDNA蛋白編碼基因拼接序列的貝葉斯分析, 無(wú)論是各支還是總后驗(yàn)概率上都是最高點(diǎn), 而基于線粒體全序列貝葉斯分析樹(shù)的各支后驗(yàn)概率都很低。不同 mtDNA蛋白編碼基因得出的后驗(yàn)概率差別很大,COX1、ND2、ATPase6總后驗(yàn)概率最高,ND3、ND4、ND4L、COX3中等,ND1、ND5、ND6、COX2、Cytb、ATPase8較差,COX1、ND2、ATPase63個(gè)基因的拼接序列后驗(yàn)概率也相對(duì)較高。紫紅笛鯛所在支的后驗(yàn)概率在0.9以上的有ND2、ND3、ND4L樹(shù)圖均顯示紫紅笛鯛和勒氏笛鯛聚為一小支, 這和13個(gè)mtDNA蛋白編碼基因拼接序列聚類一致。

圖3 笛鯛線粒體全序列貝葉斯樹(shù)圖Fig. 3 Bayesian evolutionary tree of Mitochondrial complete sequence

圖4 線粒體蛋白編碼序列基因貝葉斯樹(shù)Fig. 4 Bayesian evolutionary tree of Mitochondrial protein coding gene

圖5 笛鯛屬魚(yú)類COX 1基因全序列貝葉斯樹(shù)Fig. 5 Bayesian evolutionary tree of COX 1 complete sequences of snappers

表3 線粒體DNA部分區(qū)段貝葉斯后驗(yàn)概率Tab. 3 Bayesian posterior probability of mtDNA partial genes region

3 分析與討論

本研究基于獲取的紫紅笛鯛 mtDNA全序列和已存在的笛鯛mtDNA全序列, 探討了適用于貝葉斯系統(tǒng)進(jìn)化分析的最佳序列, 明確了紫紅笛鯛與其他7種笛鯛的進(jìn)化關(guān)系, 并推斷了笛鯛兩大支系的分化年代。

3.1 線粒體DNA全序列的高效獲取與分析

相比以直接提取獲得的總DNA為模板, 運(yùn)用長(zhǎng)距PCR的方法獲得紫紅笛鯛線粒體DNA長(zhǎng)距片段,保證了總DNA很少量時(shí)模板量足夠使用, 同時(shí)提高了模板質(zhì)量, 采用常規(guī) PCR產(chǎn)物直接送測(cè), 測(cè)序結(jié)果分析簡(jiǎn)便, 整個(gè)周期短[10,12]。此方法可方便、快速獲得高質(zhì)量的線粒體全序列。之前對(duì)笛鯛魚(yú)類研究一般常用 MEGA進(jìn)行建樹(shù)分析[6,24], 這種方式只能獲得物種間的進(jìn)化關(guān)系, 不能獲得分化年代方面的結(jié)果。本研究嘗試運(yùn)用BEAST 軟件對(duì)笛鯛全序列、mtDNA編碼區(qū)基因、mtDNA編碼區(qū)基因拼接數(shù)據(jù)集進(jìn)行了貝葉斯分析, 并采用松散分子鐘模型推斷分歧的發(fā)生時(shí)間[21], 與MrBayes相比[22], BEAST操作直觀, 而且配套有 Figtree程序, 數(shù)據(jù)分析參數(shù)可進(jìn)行及時(shí)調(diào)整, 生成的進(jìn)化樹(shù)和進(jìn)化年代可進(jìn)行方便地編輯和標(biāo)記。

3.2 基于貝葉斯法的最佳 mtDNA區(qū)段選擇

對(duì)進(jìn)化關(guān)系研究時(shí)線粒體基因組 DNA及其 13個(gè)編碼基因的適用性, 以往研究者已經(jīng)進(jìn)行了一系列基于NJ等傳統(tǒng)方法的分析。Zardoya等[23]比較了不同基因區(qū)分脊椎動(dòng)物門(mén)下的各綱的能力, 認(rèn)為COX1、Cytb、ND2、ND4、ND5都非常適合綱之間進(jìn)化研究; 陳姝君等[24]比較了硬骨魚(yú)類目間不同基因的適應(yīng)性, 認(rèn)為COX1、COX3、ND4、ND5最為出色; 郭昱嵩[25]則認(rèn)為COX 2、ND 6最為適合鱸形目各屬間進(jìn)化分析。王中鐸[6]對(duì)笛鯛屬的研究揭示出Cytb、COX1、COX2是最適合系統(tǒng)進(jìn)化分析的單基因, 13個(gè)蛋白編碼基因拼接得到組合片段建立的NJ樹(shù)的置信度比單個(gè)基因有顯著提高, 但低于線粒體全序列?？梢?jiàn), 基于分類階元不同, 區(qū)段的適用性存在差異。

本研究運(yùn)用貝葉斯法分析了笛鯛屬系統(tǒng)發(fā)生,平均后驗(yàn)概率的高低次序?yàn)? 13個(gè)蛋白編碼基因拼接序列(0.97)＞COX1＋ND2＋ATPase6(0.92)＞COX1、ND2、ATPase6、ND3和COX3(0.91～0.86)＞全序列(0.83)＞余下 8個(gè)單基因(0.82～0.59)。這與王中鐸等NJ法進(jìn)行系統(tǒng)進(jìn)化分析時(shí)的最佳線粒體區(qū)段存在差異, 同時(shí)他的研究的結(jié)果表明全序列會(huì)優(yōu)于任一編碼區(qū)。全序列的貝葉斯后驗(yàn)率明顯低于13個(gè)蛋白編碼基因拼接序列, 甚至低于COX1、ND2、ATPase6等基因, 以 13個(gè)蛋白編碼基因拼接序列進(jìn)行貝葉斯法年代推斷比全序列更可信。

3.3 笛鯛屬分化年代的貝葉斯法推斷

本研究通過(guò)紫紅笛鯛和各種笛鯛遺傳距離比較分析及貝葉斯進(jìn)化樹(shù)分析表明, 雖然紫紅笛鯛在形態(tài)上屬于大型笛鯛種, 但是其生活史表現(xiàn)出特殊性[4]。表2中K2P遺傳距離統(tǒng)計(jì)顯示紫紅笛鯛和另外3種大型笛鯛(馬拉巴笛鯛、千年笛鯛、紅鰭笛鯛)遺傳距離較遠(yuǎn), 而與勒氏笛鯛、海雞母笛鯛較近。貝葉斯樹(shù)圖更直觀地表明這種關(guān)系。

Bruno等[26]研究發(fā)現(xiàn)貝葉斯法推斷系統(tǒng)發(fā)育往往會(huì)高估節(jié)點(diǎn)的后驗(yàn)概率, 認(rèn)為只有后驗(yàn)概率大于95%的節(jié)點(diǎn)才是可信的?；谌蛄泻筒煌虻姆治鲋懈髦У暮篁?yàn)概率都未達(dá)到這一標(biāo)準(zhǔn), 基于 13個(gè)蛋白編碼基因拼接的貝葉斯分析中各個(gè)節(jié)點(diǎn)的后驗(yàn)率均大于95%, 支持紫紅笛鯛與勒氏笛鯛聚類?；趍tDNA的13個(gè)編碼區(qū)拼接序列, 貝葉斯松散分子鐘推測(cè)笛鯛屬兩大類群的分化年代距今約在61.5Ma年左右, 根據(jù)世界地質(zhì)年代表記錄這一時(shí)期處于古新世期間(65～53 Ma)的達(dá)寧階與蒙蒂階交界時(shí)期。

[1] 邢玉娜, 尹紹武, 陳國(guó)華,等. 紫紅笛鯛Lutjanus argentimaculatus(Forsska)繁殖和池塘育苗研究[J]. 現(xiàn)代漁業(yè)信息, 2005, 8: 25-27.

[2] 區(qū)又君, 勾效偉, 廖銳, 等. 紫紅笛鯛幼魚(yú)對(duì)集群密度和投餌定額的攝食反應(yīng)[J]. 動(dòng)物學(xué)雜志, 2008, 43(3): 34-38.

[3] 張俊彬, 黃良民, 陳真然. 紫紅笛鯛遺傳多樣性的AFLP分析[J]. 海洋學(xué)報(bào). 2005, 3, 27(2): 165-171.

[4] Russell D J, McDougall A J, Fletcher A S, et al. Biology, management and genetic stock structure of mangrove jack (Lutjanus argentimaculatus) in Australia[J]. FRDC, 2003, 18(9): 1456-1563.

[5] Guo Y S, Wang Z D, Liu C W, et al. Sequencing and analysis of the complete mitochondrial DNA of Russell’s snapper (L. russellii) [J]. Natural Science, 2008, 18: 1233-1238.

[6] 王中鐸, 譚圍, 郭昱嵩, 等.笛鯛屬線粒體基因組編碼序列的系統(tǒng)進(jìn)化分析能力評(píng)估[J]. 水產(chǎn)學(xué)報(bào), 2010, 34(6): 656-663.

[7] 張艷春. 大口鳒Psettodes erumei線粒體全序列的研究和鰈形目魚(yú)類系統(tǒng)進(jìn)化分析[D]. 青島: 中國(guó)海洋大學(xué), 2009.

[8] 鄭楠, 張?jiān)? 敖濱, 等. 貝葉斯推論重建乳頭瘤病毒的進(jìn)化系統(tǒng)[J]. 病毒學(xué)報(bào), 2002, 19(4): 293-300.

[9] Kawahara R, Miya M, Mabuchi K, et al. Interrelationships of the 11 gasterosteiform families (sticklebacks, pipefishes, and their relatives): A new perspective based on whole mitogenome sequences from 75 higher teleosts [J]. Molecular Phylogenetics and Evolution, 2008, 46: 224-236.

[10] 王中鐸, 郭昱嵩, 劉楚吾, 等. 軍曹魚(yú)線粒體 DNA全序列與 鲹魚(yú)宗系的系統(tǒng)進(jìn)化[J]. 水生物學(xué)報(bào), 2011, 35(2): 1-9.

[11] 李志強(qiáng), 郭寶成, 李俊兵, 等. 貝葉斯聯(lián)合模型與中國(guó)洞穴魚(yú)類分化時(shí)間的估算[J]. 科學(xué)通報(bào), 2008, 53(13): 1560-1569.

[12] 譚圍, 王中鐸, 郭昱嵩, 等. 孟加拉笛鯛線粒體基因組序列結(jié)構(gòu)及其進(jìn)化[J]. 中國(guó)生物化學(xué)與分子生物學(xué)報(bào), 2009, 25(3): 287-291.

[13] Miya M, Nishida M. Use of mitogenomic information in teleostean molecular phylogentics: a tree-based exploration under the maximum-parsimony optimality criterion[J]. Molecula Phylogenetics and Evolution, 2000, 17(3): 437-455.

[14] 王中鐸, 郭昱嵩, 劉楚吾, 等. 笛鯛屬魚(yú)類DNA分子條碼、系統(tǒng)進(jìn)化和成種機(jī)制[J]. 中國(guó)科學(xué), 2010, 6: 516-521.

[15] 盧圣棟. 現(xiàn)代分子生物學(xué)實(shí)驗(yàn)技術(shù)[M]. 北京: 中國(guó)協(xié)和醫(yī)科大學(xué)出版社, 1999: 6.

[16] Kumar S, Nei M, Dudley J, et al. MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences[J]. Brief Bioinform, 2008, 9(4): 299-306.

[17] Lowe T M, Eddy S R. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence[J]. Nucleic Acids Research, 1997, 25(5): 955-964.

[18] Drummond A J, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees[J]. BMC Evolutionary Biology, 2007, 7: 214.

[19] Drummond A J, Ho S Y, Phillips M J, et al. Relaxed phylogenetics and dating with confidence[J]. PLoS Biology, 2006, 214(7):699-710.

[20] Benton M J, Donoghue P C J. Paleontological evidence to date the tree of life [J]. Molecular Biology and Evolution, 2007, 24: 26-53.

[21] Leache A D, Reeder T W. Molecular systematics of the Eastern Fence Lizard (Sceloporus undulatus): a comparison of parsimony, likelihood, and bayesian approaches[J]. Systematic Biology, 2002, 51: 44-68.

[22] Ronquist F, Huelsenbeck J P. MrBayes 3: Bayesian phylogenetic inference under mixed models[J]. Bioinformatics, 2003, 19: 1572-1574.

[23] Zardoya R, Meyer A. Phylogenetic performance of mitochondrial protein-coding genes in resolving genes in resolving relationgships among vertebrates[J]. Molecular Phylogenetic and Evolution, 2000, 17(3): 437-455.

[24] 陳姝君, 赫崇波, 木云雷, 等.硬骨魚(yú)類線粒體基因系統(tǒng)發(fā)育信息效率分析[J].中國(guó)水產(chǎn)科學(xué), 2008, 15(1): 12-21.

[25] 郭昱嵩, 王中鐸, 焦陽(yáng), 等. 鱸形目線粒體DNA蛋白編碼基因的適應(yīng)性分析[J]. 自然科學(xué)進(jìn)展, 2008, 18(9): 7355-7358.

[26] Bruno W J. Weighted neighbor joining, a likelihood-based approach to distancebased phylogeny reconstruction[J]. Mol Biol, 2000, 17(4): 189-197.

(本文編輯: 譚雪靜)

Complete mitochondria sequence ofLutjanus argentimaculatusand Bayesian analysis

LIAO Jie, WANG Zhong-duo, GUO Yu-song, LIU Chu-wu

(Fisheries College, Key Laboratory of Aquaculture in South China Sea for Aquatic Economic Animal of Guangdong Higher Education Institutes, Guangdong Ocean University, Zhanjiang 524025, China)

Dec.,30,2011

Lutjanus argentimaculatus; mitochondrion genome; bayesian analysis; uncorrected relaxed lognormal clock; snapper fish

We got 16 543 bp Red snapper (Lutjanus argentimaculatus) complete mitochondrial DNA sequence through common PCR, with the GenBank registration number being JN182927. Our sequence and the data in Gen-Bank, red snapper’s mitochondrion gene composition were consistent with other snappers, according to the characteristics of vertebrates. Bayesian tree of complete sequence and 13 protein-coding genes and their splicing sequence indicated that: red snapper was close toL. russelliiandL. rivulatus. The evolution rate differences of different mitochondrial sections were obvious.COX1,ND2 andATPase6 Bayesian posterior rates were significantly higher than the other ten in 13 protein-coding mitochondrial genes of. The 13 protein-coding gene splicing sequences Bayesian posterior rates were significantly higher than completed sequence, and were more suitable for Bayesian uncorrected relaxed lognormal clock analysis. Thirteen protein coding gene splicing sequence’s Bayesian uncorrected relaxed lognormal clock work out snappers’ two main branches broke out in 61.5 Ma before, which was the junction point of Danian and Selandian in Paleocene.

Q349

A

1000-3096(2013)01-0062-08

2011-12-30;

2012-04-06

國(guó)家自然科學(xué)基金項(xiàng)目(31101904, 31201996); 廣東省自然科學(xué)基金項(xiàng)目(10452408801004224)

廖杰(1986-), 男, 湖南湘潭人, 碩士, 主要從事海洋生物研究工作, 電話: 13809887920, E-mail: liaojie05102@163.com;王中鐸,

, 電話: 13659795480, E-mail: aduofa@gmail.com