高天翔畢瀟瀟趙林林李創(chuàng)舉

(1.中國海洋大學(xué)海洋生物多樣性與進(jìn)化研究所, 青島 266003; 2.中國水產(chǎn)科學(xué)院長江水產(chǎn)研究所, 荊州 434000)

基于線粒體Cytb基因全序列的松江鱸群體遺傳結(jié)構(gòu)分析

高天翔1畢瀟瀟1趙林林1李創(chuàng)舉2

(1.中國海洋大學(xué)海洋生物多樣性與進(jìn)化研究所, 青島 266003; 2.中國水產(chǎn)科學(xué)院長江水產(chǎn)研究所, 荊州 434000)

對松江鱸(Trachidermus fasciatusHeckel)中國沿海7個群體和日本有明海群體的線粒體Cytb基因全序列進(jìn)行了測定、分析。結(jié)果顯示: 47個個體共檢測到31個單倍型, 8個群體均呈現(xiàn)出較高的單倍型多樣性(0.60—1.00)和較低的核苷酸多樣性(0.0005—0.0041)的特點(diǎn)。AMOVA分析結(jié)果及單倍型鄰接關(guān)系樹和單倍型網(wǎng)絡(luò)關(guān)系圖均顯示松江鱸分為中國和日本兩個世系, 而中國世系的7個群體未呈現(xiàn)出明顯的地理遺傳結(jié)構(gòu)?；诤塑账酜imura雙參數(shù)替代模型計(jì)算得出的中國和日本兩個世系的凈遺傳距離, 再參照其他硬骨魚類線粒體Cytb基因2%/Ma(百萬年)的分歧速率, 推測松江鱸中日兩個世系間分化時間約為41萬年前。對中國世系進(jìn)行群體歷史動態(tài)分析, 中性檢驗(yàn)結(jié)果均為負(fù)值且顯著, 核苷酸不配對分布呈單峰型, 表明松江鱸中國世系曾發(fā)生過群體擴(kuò)張, 其擴(kuò)張時間大約為12萬年前。

松江鱸;Cytb基因; 全序列; 遺傳結(jié)構(gòu); 群體擴(kuò)張

松江鱸(Trachidermus fasciatus Heckel)隸屬于鲉形目(Scorpaeniformes)、杜父魚科(Cottidae)、松江鱸屬(Trachidermus), 俗稱四鰓鱸、花鼓魚、媳婦魚等。松江鱸是一年生肉食性降河洄游魚類, 其營養(yǎng)價值和藥用價值高[1], 與黃河鯉魚、松花江鮭魚、黑龍江興凱湖鲌并列為“中國四大淡水名魚”。松江鱸為西北太平洋所特有, 僅分布于中國、朝鮮半島和日本南部, 我國分布于黃渤海和東海沿岸以及沿岸諸河下游及河口[2,3]。近年來, 由于工業(yè)有害污水及農(nóng)藥化肥等對江河口和沿海的嚴(yán)重污染, 松江鱸數(shù)量急劇下降; 與此同時, 攔河建壩、水閘等水利設(shè)施的大量興建, 截?cái)嗔怂山|的洄游通道, 導(dǎo)致其數(shù)量急劇下降, 進(jìn)入內(nèi)陸水域的松江鱸魚苗數(shù)量大為減少, 嚴(yán)重影響了其資源的恢復(fù), 在一些原有分布記錄的水域已難覓其蹤跡。松江鱸現(xiàn)為我國野生動物重點(diǎn)保護(hù)二級水生動物[4], 2004年被列入中國物種紅色名錄。

線粒體DNA(Mitochondrial DNA, mtDNA)結(jié)構(gòu)簡單, 呈母性遺傳, 進(jìn)化速度快且不發(fā)生重組, 是一種應(yīng)用較廣的分子標(biāo)記, 其中一些基因已被廣泛用于魚類群體遺傳學(xué)和系統(tǒng)發(fā)育關(guān)系研究[5—7]。細(xì)胞色素b (Cytochrome b, Cytb)基因進(jìn)化速度適中, 不僅容易使用通用引物擴(kuò)增和測序, 而且它還是mtDNA上唯一的結(jié)構(gòu)和功能被了解的較為清楚的蛋白編碼基因[8,9]。Cytb被認(rèn)為是了解種質(zhì)資源狀況和群體遺傳結(jié)構(gòu)的理想工具, 近年來已成為魚類遺傳多樣性研究中常用的標(biāo)記之一[7]。

目前對松江鱸魚的報(bào)道多集中在資源調(diào)查與保護(hù)、繁育及生理生態(tài)等方面[3,10—12], 關(guān)于遺傳多樣性的研究相對較少。王金秋等[13]、徐建榮等[14,15]、劉海林等[16]分別用同工酶、AFLP、SSR和線粒體控制區(qū)標(biāo)記對我國黃渤海松江鱸群體的遺傳多樣性進(jìn)行了初步研究, 但樣品分布點(diǎn)均較少。本研究基于線粒體Cytb基因全序列分析了日本有明海和中國渤海、黃海及東海近海松江鱸群體的遺傳結(jié)構(gòu), 旨在為松江鱸種質(zhì)資源的保護(hù)和管理提供參考依據(jù)。

1 材料與方法

1.1 實(shí)驗(yàn)材料

松江鱸樣品于2007至2010年采自于日本有明海(Ariake Sea)、大連(DL)、丹東(DD)、秦皇島(QHD)、東營(DY)、榮成(RC)、文登(WD)、杭州灣(HZB)8個地點(diǎn)(圖1)。除日本有明海松江鱸樣本數(shù)為12尾 外, 其他群體均為5尾,共計(jì)47尾。樣品經(jīng)形態(tài)學(xué)鑒定后,均用95%乙醇固定, 于?20℃冰箱中保存。

圖1 松江鱸采樣地點(diǎn)Fig.1 Sampling sites of T.fasciatus

1.2 實(shí)驗(yàn)方法

基因組DNA的提取 取松江鱸肌肉組織約100 mg, 采用標(biāo)準(zhǔn)的酚-氯仿方法提取基因組DNA,將乙醇沉淀后的基因組DNA于100 μL超純水溶解, 4℃保存?zhèn)溆谩?/p>

目的片段的PCR擴(kuò)增 根據(jù)GenBank中報(bào)道的相近種Cytb全序列及現(xiàn)有的通用引物[17], 共采用了3對引物對松江鱸線粒體Cytb基因全序列進(jìn)行擴(kuò)增, 引物分別為:

PCR 反應(yīng)體系總體積為50μL, 其中: 10×PCR緩沖液(200 mmol/L Tris-HCl, pH 8.4; 200 mmol/L KCl; 100 mmol/L (NH4)2SO4; 15 mmol/L MgCl2) 5 μL, dNTP 200 μmol/L, 引物各0.2 μmol/L, Taq酶1.25 U (大連寶生物公司), 模板DNA 20 ng, 加水至50 μL。應(yīng)用Eppendorf熱循環(huán)儀進(jìn)行PCR反應(yīng), 反應(yīng)條件如下: 首先94℃預(yù)變性5min, 之后進(jìn)行35個循環(huán),每個循環(huán)包括94℃變性45s、50℃退火45s、72℃延伸45s, 最后72℃延伸10 min。以上反應(yīng)均設(shè)陰性對照以排除DNA污染的情況。用1.5%瓊脂糖凝膠電泳檢測PCR擴(kuò)增產(chǎn)物, 凝膠成像系統(tǒng)拍照。

純化回收及測序 PCR產(chǎn)物經(jīng)1%TAE瓊脂糖凝膠電泳檢測, 產(chǎn)物電泳后于紫外燈下切割目的條帶, 用DNA膠回收試劑盒進(jìn)行產(chǎn)物的純化回收。純化產(chǎn)物送上海桑尼生物有限公司進(jìn)行雙向測序,測序反應(yīng)采用與PCR反應(yīng)一致的引物。

1.3 數(shù)據(jù)分析

序列分析 所有序列均由Dnastar軟件包進(jìn)行編輯、校對和排序, 并對排序結(jié)果進(jìn)行分析和手工校正。利用Arlequin3.1[18]軟件計(jì)算堿基組成、變異位點(diǎn)數(shù)、單倍型數(shù)、平均核苷酸差異數(shù)、單倍型多樣性和核苷酸多樣性。

群體結(jié)構(gòu)分析 應(yīng)用Modeltest3.7[19]對線粒體Cytb基因全序列進(jìn)行最佳替換模型篩選, 并估計(jì)相關(guān)參數(shù)。根據(jù)篩選的模型和相關(guān)參數(shù), 利用Paup4.0[20]構(gòu)建單倍型鄰接關(guān)系樹來探討單倍型的遺傳結(jié)構(gòu)。利用Mega4.0[22]軟件基于Kimura雙參數(shù)模型計(jì)算中國世系和日本世系之間的凈遺傳距離以估算中日松江鱸的分化時間。使用分子變異分析(AMOVA)來評估群體間遺傳變異, 通過1000次重抽樣來檢驗(yàn)不同遺傳結(jié)構(gòu)水平上協(xié)方差的顯著性。采用分化固定指數(shù)(Fst)來評價兩兩群體間的遺傳差異, 通過1000次重抽樣來檢驗(yàn)兩兩群體間Fst的顯著性。通過設(shè)定兩種AMOVA分析來檢驗(yàn)松江鱸的群體遺傳結(jié)構(gòu): 一是將8個群體都劃分為一個組群以驗(yàn)證群體間是否具有顯著的遺傳分化; 二是按樣品的地理來源將松江鱸8個群體劃分為2個組群,分別對應(yīng)不同的海區(qū)以驗(yàn)證是否存在顯著的地理結(jié)構(gòu), 日本有明海群體劃分為日本組群, 其余7個群體劃分為中國組群。

群體歷史動態(tài)分析 在Arlequin3.1[18]軟件中, 采用中性檢驗(yàn)和核苷酸不配對分布兩種方法來檢測松江鱸的群體歷史動態(tài)。首先由Tajima’s D[23]檢驗(yàn)和Fu’s Fs[24]檢驗(yàn)來進(jìn)行中性檢驗(yàn), 其次進(jìn)行核苷酸不配對分布檢測。采用最小方差法來檢驗(yàn)核苷酸不配對的觀測值和群體擴(kuò)張模型下的預(yù)期分布之間是否一致, 采用廣義非線性最小方差法來估算擴(kuò)張參數(shù)τ。根據(jù)公式τ= 2ut計(jì)算實(shí)際的擴(kuò)張時間, 其中u為所研究的整個序列長度的突變速率, t是自群體擴(kuò)張開始到現(xiàn)在的時間。參照其他硬骨魚類的線粒體Cytb基因2%/Ma(百萬年)的分歧速率[25,26]來估算松江鱸群體擴(kuò)張時間。

2 結(jié)果

2.1 序列變異及遺傳多樣性

所分析的8個群體、47尾松江鱸個體的線粒體Cytb基因全序列的長度均為1141 bp, 其中C、T、A和G含量分別為30.0%、28.2%、23.0%和16.8%, A+T的含量(51.2%)略高于C+T的含量(48.8%), Cytb基因組成具有較大的偏向性, G的含量較低。在所有47個個體的1141 bp的序列中, 共檢測到38個變異位點(diǎn), 約占總序列長度的3.33%, 簡約信息位點(diǎn)30個, 單一變異位點(diǎn)8個。這些變異位點(diǎn)共定義39個核苷酸替代(37個轉(zhuǎn)換和2個顛換), 未發(fā)現(xiàn)插入、缺失現(xiàn)象(圖2)。47個個體中共檢測到31單倍型, 所有單倍型已經(jīng)提交到GenBank, 登錄號為JX079997-JX080027。在所有單倍型中, 2個單倍型為共享單倍型, 其他29個單倍型均為群體特有單倍型; 共享單倍型中, Hap1為東營、大連和杭州灣三個群體共享, Hap16為文登和丹東兩個群體共享。

單倍型多樣性指數(shù)整體較高, 而核苷酸多樣性指數(shù)相對較低(表1)。其中, HZW、WD、QHD和DD群體的單倍型多樣性指數(shù)最高, 均為1.00, RC和DY群體最低, 均為0.6; WD群體的核苷酸多樣性指最高, 為0.0041, DY群體最低, 為0.0005。所有樣本作為一個群體進(jìn)行數(shù)據(jù)分析, 單倍型多樣性指數(shù)0.97, 核苷酸多樣性指數(shù)為0.0059。

2.2 群體遺傳結(jié)構(gòu)

從松江鱸8個地理群體的遺傳距離(表2)來看,群體內(nèi)的遺傳距離為0.0005—0.0041, 中國組群各群體間遺傳距離為0.0014—0.0058, 日本群體和中國各群體間的遺傳距離0.0095—0.0120(表2)。此結(jié)果表明, 松江鱸日本群體和中國各群體間的遺傳距離明顯大于中國各群體間的遺傳距離, 中國各群體間的遺傳距離與群體內(nèi)的遺傳距離處于同一水平,無明顯遺傳分化。兩兩群體之間的Fst值(表2)顯示,日本群體和中國各群體之間的Fst值為0.8008—0.8922, 統(tǒng)計(jì)檢驗(yàn)均顯著(P<0.05); 而中國各群體之間的Fst值為0.0736—0.3984, 統(tǒng)計(jì)檢驗(yàn)均不顯著(P>0.05); 中國組群和日本組群之間的Fst值為0.7499, 統(tǒng)計(jì)檢驗(yàn)顯著(P<0.05)。AMOVA 分析結(jié)果顯示松江鱸中國組群和日本組群之間已經(jīng)產(chǎn)生一定程度的遺傳分化, 而中國各群體內(nèi)不存在顯著的遺傳結(jié)構(gòu)?；贙imura雙參數(shù)模型計(jì)算得到的松江鱸中國組群和日本組群間凈遺傳距離為0.082, 以線粒體Cytb基因序列2%/百萬年的核苷酸分化速率計(jì)算[25,26], 中國組群和日本組群的分化年代大約在41萬年前。

圖2 松江鱸線粒體Cytb 基因序列變異位點(diǎn)Fig.2 Variable sites of mitochondrial Cytb gene in T.fasciatus

以圖們江杜父魚(Cottus hangiongensis)和賴氏杜父魚(Cottus reinii)(GenBank登錄號分別為NC_014851、NC_004404)為外群, 基于最適模型進(jìn)行1000次自展檢驗(yàn), 利用PAUP4.0構(gòu)建鄰接系統(tǒng)樹(圖3)。結(jié)果顯示松江鱸群體明顯分為中國和日本兩個世系, 而中國各群體之間沒有明顯的地理遺傳結(jié)構(gòu)。基于簡化的中介網(wǎng)絡(luò)法構(gòu)建的松江鱸單倍型網(wǎng)絡(luò)關(guān)系圖(圖4)亦明顯分為中國和日本兩個世系, 其中中國世系內(nèi)單倍型Hap1位于網(wǎng)絡(luò)圖中心, 其他單倍型均經(jīng)過一步或幾步突變連接到Hap1, 由此推測Hap1為中國群體的原始單倍型。

表1 松江鱸各群體的樣品信息和遺傳多樣性參數(shù)Tab.1 Sample information and genetic diversity parameters in different populations of T.fasciatus

表2 松江鱸群體內(nèi)遺傳距離(對角線)及兩兩群體間遺傳距離(上對角線)和遺傳分化系數(shù)(Fst)(下對角線)Tab.2 Pairwise genetic distances within (diagonal) population, and genetic distance (below diagonal), fixation index (Fst) (above diagonal) between populations of T.fasciatus

2.3 群體歷史動態(tài)

松江鱸中國世系內(nèi)不存在明顯的遺傳分化, 將其看作一個群體進(jìn)行歷史動態(tài)分析; 中性檢驗(yàn)表明(表2), Tajima’s D 和Fu’s Fs值均為負(fù)值, 并且統(tǒng)計(jì)檢驗(yàn)均達(dá)顯著水平(P<0.05), 單倍型核苷酸不配對分布成單峰形, 與群體擴(kuò)張下的預(yù)期相符合(圖5),提示松江鱸經(jīng)歷了群體歷史擴(kuò)張事件。采用最小方差法對核苷酸不配對的觀測值和群體擴(kuò)張模型下的預(yù)期值之間的一致性進(jìn)行擬合優(yōu)度檢驗(yàn), 結(jié)果顯示SSD值和Raggedness指數(shù)均較小, 統(tǒng)計(jì)檢驗(yàn)均不顯著(P>0.05)(表3), 表明所觀測的核苷酸不配對分布沒有顯著偏離群體擴(kuò)張模型。采用廣義非線性最小方差法估算松江鱸中國世系擴(kuò)張參數(shù)τ值為2.793 (95% CI: 1.947—5.178), 據(jù)此推算松江鱸的群體擴(kuò)張時間約為12萬年前(95%CI: 7.9萬年—22.7萬年)。

3 討論

圖 3 基于Cytb基因全序列構(gòu)建的松江鱸鄰接系統(tǒng)樹(外群為圖們江杜父魚和賴氏杜父魚)Fig.3 NJ phylogenetic trees based on Cytb gene complete sequences of T.fasciatus (C.Hangiongensis and C.Reinii as outgroup)

松江鱸8個群體的線粒體Cytb基因全序列中,共有38個變異位點(diǎn)、31單倍型, 單倍型多樣度為0.97, 呈現(xiàn)出較高的遺傳多樣性, 但其核苷酸多樣性僅為0.0059, 處于一個較低水平, 與劉海林等[16]對黃渤海松江鱸控制區(qū)序列研究的結(jié)果相一致, 這也符合Grant和Bowen[27]提出的海水魚類不同單倍型多樣性和核苷酸多樣性間的第二種類型, 即較高的單倍型多樣性和較低的核苷酸多樣性。西北太平洋區(qū)域有些魚類如黃姑魚[28]、花鱸[29]和銀鯧[30]也具有這種類型的遺傳多樣性模式, 顯示該類型魚類群體可能經(jīng)歷過歷史擴(kuò)張事件, 即由一個較小的有效群體在短時間內(nèi)快速成長為一個大的群體。在群體擴(kuò)張的過程中, 隨著群體數(shù)量的增加, 單倍型多樣度會有所提高, 但沒有足夠的時間來積累核苷酸產(chǎn)生的變異, 因此會造成單倍型多樣度較高而核苷酸多樣度較低的遺傳多樣性模式[26]。

生物地理群之間基因流的阻隔和環(huán)境因子之間的差異可導(dǎo)致不同群體之間遺傳結(jié)構(gòu)的顯著差異,當(dāng)群體間存在較大的基因流時, 其遺傳分化不明顯。然而不同地區(qū)的生境差異是產(chǎn)生和保持地理群體間遺傳結(jié)構(gòu)差異的重要因素[31]。在本研究中, AMOVA分析結(jié)果、單倍型鄰接樹和單倍型網(wǎng)絡(luò)關(guān)系圖均表明松江鱸中國沿海各群體和日本有明海群體間的遺傳差異顯著, 可明顯分為中國和日本兩大世系, 其結(jié)果與白姑魚[32]研究結(jié)果相似?；诤塑账酜imura雙參數(shù)替代模型計(jì)算的松江鱸中國和日本兩個世系的凈遺傳距離達(dá)0.082, 表明中日兩個世系分歧時間約為41萬年以前, 其分化可能是受到更新世冰期的影響。很多研究者發(fā)現(xiàn)更新世冰期對于海洋魚類的群體遺傳結(jié)構(gòu)產(chǎn)生了巨大的影響[29,33]。第四紀(jì)末期, 全球氣候經(jīng)歷了一系列冰期-間冰期的變化, 在過去的約80萬年里, 氣候波動主要以一個約為10萬年的周期進(jìn)行[34]。第五次冰期結(jié)束時間大約為42萬年前[35], 與松江鱸中日世系的分化時間相吻合。伴隨著第五次冰盛期的到來, 海平面下降約120—140 m, 松江鱸中國和日本群體可能會被隔離在相應(yīng)邊緣海的避難所里; 冰期過后, 隨海平面上升, 被隔離的群體重新擴(kuò)張, 但中日之間長距離的地理隔離和松江鱸幼體有限的擴(kuò)散能力等原因?qū)е缕浠蚪涣髦袛? 隨著時間的推移, 中國松江鱸群體和日本群體逐漸分化成兩個不同的世系。然而,中國世系的松江鱸各群體之間沒有明顯的地理遺傳結(jié)構(gòu), 即各群體不存在顯著的遺傳差異, 與之前關(guān)于松江鱸ISSR、AFLP和線粒體控制區(qū)的研究結(jié)果[13—16]相一致。分析可能有以下原因: 松江鱸在中國沿海群體擴(kuò)張的時間尚較晚, 群體缺乏足夠的時間在遷移和漂變之間取得平衡, 因此其各群體之間不存在顯著的遺傳分化[36]; 松江鱸繁殖和幼魚成長階段均在河口近海完成, 從2月到6月持續(xù)將近五個月[1], 此階段有可能隨海洋沿岸流被動擴(kuò)散, 使不同群體基因交流的機(jī)會增加, 遺傳分化程度降低;其次, 松江鱸大都生活在河流至河口近海區(qū)域, 其生活的地理構(gòu)架和環(huán)境因子相似, 且松江鱸屬廣溫廣鹽性種類, 不同地區(qū)的溫度和鹽度對其構(gòu)成的選擇壓力不大, 因此其群體間的遺傳差異較小; 此外,松江鱸為一年生魚種, 近年來的增養(yǎng)殖也可能造成群體間的基因交流。將中國世系當(dāng)作一個大的群體進(jìn)行群體歷史動態(tài)分析, 中性檢驗(yàn)的結(jié)果均為負(fù)值且顯著, 核苷酸不配對分布呈單峰型, 表明此群體歷史上曾經(jīng)發(fā)生過群體擴(kuò)張, 計(jì)算得到松江鱸中國世系的擴(kuò)張時間大約為12萬年前?？赡苁悄┐伪诮Y(jié)束以后, 海平面上升后, 松江鱸的棲息地發(fā)生了大面積擴(kuò)張, 導(dǎo)致其群體發(fā)生了擴(kuò)張。

圖 4 基于簡約法構(gòu)建的松江鱸單倍型網(wǎng)絡(luò)關(guān)系圖Fig.4 Median-networks showing genetic relationship among haplotypes of T.fasciatus based on parsimony method

遺傳多樣性不僅是生物多樣性的基礎(chǔ)和核心,也是生物物種進(jìn)化潛能的保證[37], 遺傳多樣性的降低將導(dǎo)致物種對環(huán)境適應(yīng)能力的下降, 復(fù)雜多變的生存環(huán)境將對其群體的生存產(chǎn)生重大的威脅。松江鱸的遺傳多樣性總體水平較高, 說明其資源量雖然銳減可能并未影響到遺傳多樣性, 若及時采取有效的保護(hù)策略, 其資源量恢復(fù)的潛在能力較大。研究結(jié)果顯示松江鱸中國群體和日本群體已經(jīng)產(chǎn)生了一定的遺傳分化, 在不同研究結(jié)果難以獲得參比的情況下, 中國沿海的松江鱸是否作為一個隨機(jī)交配的群體進(jìn)行管理需要謹(jǐn)慎對待。今后, 有必要采用其他分子標(biāo)記(如SSR、AFLP、SNP等)對松江鱸群體結(jié)構(gòu)及其遺傳多樣性狀況進(jìn)行全面研究, 促進(jìn)其資源的保護(hù)和可持續(xù)發(fā)展。

圖 5 松江鱸中國世系的核苷酸不配對分布圖(柱狀圖顯示觀測值, 曲線為群體擴(kuò)張模型下的預(yù)期分布)Fig.5 Nucleotide mismatch distribution of T.fasciatus in Chinese lineage (Bars are the observed pairwise differences and solid line represents the expected mismatch distributions under the sudden expansion model)

表3 松江鱸中國群體中性檢驗(yàn)、核苷酸不配對分布及擬合優(yōu)度檢驗(yàn)Tab.3 Neutral test, mismatch distribution of Nucleotide and test of goodness of fit of T.fasciatus populations in China

致謝:

吳建新、國寶富、郭棟、王茂林、宋林及日本九州大學(xué)鬼谷徳雄博士等幫助采集松江鱸樣品, 在此表示感謝。

[1] Wang J Q.Advance in studies on the ecology and reproductive biology of Trachidermus fasciatus Heckel [J].Acta Hydrobiologica Sinica, 1999, 23(6): 729—734 [王金秋.松江鱸的生態(tài)學(xué)和繁殖生物學(xué)的研究進(jìn)展.水生生物學(xué)報(bào), 1999, 23(6): 729—734]

[2] Wang Y H.Notes on the scientific name, type-locality and geographical distribution of roughskin sculpin (Trachidermis fasciatus Heckel 1840) [J].Marine Fisheries, 2006, 28(4): 209—303 [王幼槐.關(guān)于松江鱸學(xué)名和模式產(chǎn)地以及地理分布之探討.海洋漁業(yè), 2006, 28(4): 209—303]

[3] Wang J Q, Cheng G, Tang Z P.The distribution of roughskin sculpin (Trachidermis fasciatus Heckel) in Yalu River Basin, China [J].Journal of Fudan University (Natural Science), 2001, 4(5): 471—475 [王金秋, 成功, 唐作鵬.鴨綠江流域中國境內(nèi)松江鱸的分布.復(fù)旦學(xué)報(bào)(自然科學(xué)版), 2001, 4(5): 471—475]

[4] Yue P Q, Chen Y Y.China Red Data Book of Endangered Animals—Fish [M].Beijing: Academic Press.1998, 240—243 [樂佩琦, 陳宜瑜.中國瀕危動物紅皮書—魚類.北京:科學(xué)出版社.1998, 240—243]

[5] Guo L, Li J, Wang Z S, et al.Phylogenetic relationships of noodle-fishes (Osmeriformes: Salangidae) based on four mitochondrial genes [J].Acta Hydrobiologica Sinica, 2011, 35(3): 449-459 [郭立, 李 雋, 王忠鎖, 等.基于四個線粒體基因片段的銀魚科魚類系統(tǒng)發(fā)育.水生生物學(xué)報(bào), 2011, 35(3): 449-459]

[6] Liu Z J, Cordes J F.DNA marker technologies and their applications in aquaculture genetics [J].Aquaculture, 2004, 238(1—4): 1—37

[7] Xiao W H, Zhang Y P.Genetics and evolution of mitochondrial DNA in fish [J].Acta Hydrobiologica Sinica, 2000, 24(4): 384—391 [肖武漢, 張亞平.魚類線粒體DNA 的遺傳與進(jìn)化.水生生物學(xué)報(bào), 2000, 24(4): 384—391]

[8] Zardoyar R, Meyer A.Phylogenetic performance mitochondrial protein coding genes in resolving relationship among vertebrates [J].Molecular Biology and Evolution, 1996, 13(7): 933—942

[9] Lwin D W, Kocher T D, Wilson A C.Evolution of the cytoehrome b gene of mammals [J].Journal of Molecular Evo1ution, 1991, 32(2): 128—144

[10] Wang J Q, Pan L D, Liang T H, et al.Preliminary study on embryonic development of the roughskin sculpin, Trachidermus fasciatus (Scorpaeniformes: Cottidae) [J].Journal of Fudan University (Natural Science), 2004, 43(2): 250—254 [王金秋, 潘連德, 梁天紅, 等.松江鱸魚(Trachidermus fasciatus)胚胎發(fā)育的初步觀察.復(fù)旦學(xué)報(bào)(自然科學(xué)版), 2004, 43(2): 250—254]

[11] Takeshita N, Onikura N, Matsui S, et al.Embryonic, larval and juvenile development of the roughskin sculpin, Trachidermus fasciatus (Scorpaeniformes, Cottidae) [J].Ichthyological Research, 1997, 44(3): 257—265

[12] Norio O, Naohiko T, Seiichi M, et al.Estimation of the frequency of maturity and sexual differences in the maturation period of the roughskin sculpin Trachidermus fasciatus [J].Fisheries Science, 2007, 73(3): 735—737

[13] Wang J Q, Shi C, Cheng G.Isozyme Analysis for a wild population of roughskin sculpin, Trachidermus fasciatus, from the Yalu River Basin in China [J].Journal of Fudan University (Natural Science), 2002, 41(6): 688—693 [王金秋, 石椿, 成功.鴨綠江流域中國境內(nèi)松江鱸魚(Trachidermus fasciatus)天然群體的同工酶分析.復(fù)旦學(xué)報(bào)(自然科學(xué)版), 2002, 41(6): 688—693]

[14] Xu J R, Han X L, Li N, et al.Analysis of genetic diversity in roughskin sculpin Trachidermus fasciatus by AFLP markers [J].Journal of Dalian Fisheries University, 2008, 23(6): 437—441 [徐建榮, 韓曉磊, 李寧, 等.松江鱸群體遺傳多樣性的AFLP 分析.大連水產(chǎn)學(xué)院學(xué)報(bào), 2008, 23(6): 437—441]

[15] Xu J R, Han X L, Yu J F, et al.Analysis of genetic diversity of Trachidermus fasciatus by ISSR [J].Freshwater Fisheries, 2009, 139(11): 21—24 [徐建榮, 韓曉磊, 郁建鋒, 等.松江鱸群體遺傳多樣性的ISSR分析.淡水漁業(yè), 2009, 139(11): 21—24]

[16] Liu H L, Zhang Q, Tang Y L, et al.Structure and genetic diversity of mtDNA D-Loop sequences among Trachidermus fasciatus stocks in Yellow Sea and Bohai Sea of China [J].Marine Science Bulletin, 2010, 29(3): 283—288 [劉海林,章群, 唐優(yōu)良, 等.黃渤海松江鱸魚線粒體控制區(qū)結(jié)構(gòu)與序列多態(tài)性分析.海洋通報(bào), 2010, 29(3): 283—288]

[17] Xiao W H, Zhang Y P, Liu H Z.Molecular systematics of Xenocyprinae (Teleostei: Cyprinidae): taxonomy, biogeography and coevolution of a special group restricted in East Asia [J].Molecular Phylogenetics and Evolution, 2001, 18(2): 163—173

[18] Excoffier L, Laval G, Schneider S.Arlequin (version 3.0): An integrated software package for population genetics data analysis [J].Evolutionary Bioinformatics, 2005, 1: 47—50

[19] Posada D, Crandall K A.Modeltest: testing the model of DNA substitution [J].Bioinformatics, 1998, 14(9): 817—818

[20] Swofford D L.Phylogenetic Analysis Using Parsimony (and other methods), Version 4 [Z].Sunderland: Sinauer Associates, MA.2002

[21] Bandelt H J, Macaulay V A, Richards M, et al.Median networks: speedy construction and gredy reduction, one simulation, and two case studies from human mtDNA [J].Molecular Phylogenetics and Evolution, 2000, 16(1): 8—28

[22] Tamura K, Dudley J, Nei M, et al.MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0 [J].Molecular Biology and Evolution, 2007, 24(8): 1596—1599

[23] Tajima F.Statistical method fortesting the neutral muation hypothesis by DNA polymorphism [J].Genetics, 1989, 123(3): 585—595

[24] Fu Y X.Statistical tests of neutrality of mutations against population growth, hitch hiking and background selection [J].Genetics, 1997, 147(2): 915—925

[25] Bermingham E S, McCafferty A, Martin P.Molecular Systematics of Fishes [M].New York: Academic Press.1997, 113—128

[26] Bowen B W, Bass A L, Rocha L A, et al.Phylogeography of the trumpetfishes (Aulostomus): ring species complex on a global scale [J].Evolution, 2001, 55(5), 1029—1039

[27] Grant W S, Bowen B W.Shallow population histories in deep evolutionary linages of marine fishes: Insights from sardines and anchovies and lessons for conservation [J].Journal of Heredity, 1998, 89(5): 415—426

[28] Han ZQ, Gao T X, Yanagimoto T, et al.Genetic population structure of Nibea albiflora in the Yellow and East China seas [J].Fisheries Science, 2008, 74(3): 544—552

[29] Liu J X, Gao T X, Yokogawa K, et al.Differential population structuring and demographic history of two closely related fish species, Japanese sea bass (Lateolabrax japonicus) and spotted sea bass (Lateolabrax maculatus) in Northwestern Pacific [J].Molecular Phylogenetics and Evolution, 2006, 39(3): 799—811

[30] Zhao F, Zhuang P, Zhang L Z, et al.Population genetic structure of Pampus argenteus in the South Yellow Sea and East China Sea based on the mitochondrial Cytb sequence [J].Acta Hydrobiologica Sinica, 2011, 35(5): 746—752 [趙峰,莊平, 章龍珍, 等.基于線粒體Cytb基因的黃海南部和東海銀鯧群體遺傳結(jié)構(gòu)分析.水生生物學(xué)報(bào), 2011, 35(5): 746—752]

[31] Levinton J S, Suchanek T H.Geographic variation, niche breadth and genetic differentiation at different geographic scales in the mussels Mytilus californianus and M.edulis [J].Marine Biology, 1978, 49(4): 363—375

[32] Han Z Q, Gao T X, Yanagimoto T, et al.Deep phylogeographic break among Pennahia argentata (Sciaenidae, Perciformes) populations in the Northwestern Pacific [J].Fisheries Science, 2008, 74(4): 770—780

[33] Liu J X, Gao T X, Zhuang Z M, et al.Late Pleistocene divergence and subsequent population expansion of two closely related fish species, Japanese anchovy (Engraulis japonicus) and Australian anchovy (Engraulis australis) [J].Molecular Phylogenetic and Evolution, 2006, 40(3): 712—723

[34] Lambeck K, Esat T M, Potter E K.Links between climate and sea levels for the past three million years [J].Nature, 2002, 419: 199—206

[35] Petit J R, Jouzel J, Raynaud D, et al.Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica [J].Nature, 1999, 399: 429—436

[36] Slatkin M.Isolation by distance in equilibrium and non equilibrium populations [J].Evolution, 1993, 47(1): 264—279

[37] Loeschcke V, Tomiuk J, Jian S K, et al.Conservation Genetics [M].Basel: Birkhauser Veriag.1994, 37—53

POPULATION GENETIC STRUCTURE OF ROUGHSKIN SCULPIN TRACHIDERMUS FASCIATUS BASED ON THE MITOCHONDRIAL CYTb SEQUENCE

GAO Tian-Xiang1, BI Xiao-Xiao1, ZHAO Lin-Lin1and LI Chuang-Ju2
(1.Institute of Evolution & Marine Biodiversity, Ocean University of China, Qingdao 266003, China; 2.Yangtze River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Jingzhou 434000, China)

Roughskin sculpin, Trachidermus fasciatus is a catadromous fish, and it is distributed along the coast of Northwest Pacific.The fish population has been greatly decreased in recent years due to degradation of habitats, water pollution and dam construction.In order to protect the fish population, it is necessary to study genetic structure of the fish population, which is very important for fishery management and conservation.Forty seven individuals used in this study were collected from 8 sampling locations in the coast of China and Japan, i.e., Hangzhou Bay (HZB), Wendeng (WD), Rongcheng (RC), Dongying (DY), Qinhuangdao (QHD), Dalian (DL), Dandong (DD) and Ariake Sea.The mitochondrial cytochromeb(Cytb) gene sequence ofT.fasciatuswas amplified by polymerase chain reaction (PCR) technology and sequenced on an automatic sequencer with both forward and reverse primers.The sequences were analyzed by using some software, which included Dnastar, Arlequin3.1, Mega4.0, Modeltest3.7 and Paup4.0.The results were as follows: the Cytb sequence of each specimen was 1141 bp and used for genetic diversity analysis.There was neither insertion nor deletion found among 38 mutations of nucleotide acids and 31 haplotypes were found among 47 individuals from 8 sampling sites.All 8 populations were characterized by high haplotype diversity (0.60—1.00) and low nucleotide diversity (0.0005—0.0041).The analysis result of molecular variance (AMOVA), median-networks and NJ phylogenetic trees showed that there was large genetic differentiation between those two regions, while there was no genetic variance among populations of the Chinese lineage.Net average genetic distance between two lineages was 0.82%.Applying divergence rate of Cytb gene sequence, the divergence of lineage China and Japan occurred about 410000 years ago.The mismatch distribution of pairwise nucleotide and the negatively selective neutrality test suggest that a recent population expansion had occurred in Chinese lineages, and the time was estimated to be about 120000 years before present (during the late Pleistocene).Two distinct lineages found inT.fasciatuswere probably related to the Pleistocene coastal glaciations and the long geographical distance.Dispersal ability, coastal currents and the population expansion may be responsible for the homegeneity among the Chinese populations.

Trachidermus fasciatus; Cytb gene; Complete sequence; Genetic structure; Population expansion

Q347

A

1000-3207(2013)02-0199-09

10.7541/2013.5

2011-12-28;

2012-10-29

海洋公益項(xiàng)目(201105005); 中國水產(chǎn)科學(xué)院長江水產(chǎn)研究所開放課題(LFBCU0713)資助

高天翔(1962—), 遼寧遼陽人; 博士; 主要從事漁業(yè)資源生物學(xué)及群體遺傳學(xué)研究。E-mail: gaozhang@ouc.edu.cn

李創(chuàng)舉, E-mail: lcj@yfi.ac.cn