朱威霖，肖 珊，楊春玲，陳慧芳，3，何震晗，3，曾地剛，陳秀荔*，彭 敏*

(1.廣西水產(chǎn)科學(xué)研究院/廣西水產(chǎn)遺傳育種與健康養(yǎng)殖重點實驗室，廣西 南寧 530021；2.廣西水產(chǎn)引育種中心，廣西 南寧 530031；3.廣西大學(xué)動物科學(xué)技術(shù)學(xué)院，廣西 南寧 530004)

【研究意義】黃鰭棘鯛(Acanthopagruslatus)又名闊黒鯛，隸屬于鱸形目(Percoiformes)鯛科(Sparide)棘鯛屬(Acanthopagrus)，為淺海暖水性底層魚類，原產(chǎn)于太平洋和北印度洋，廣泛分布在我國、日本、朝鮮、菲律賓和科威特等國家，在我國主要分布于東南沿海地區(qū)，在廣東俗稱黃腳魚立、黃絲魚立、黃立魚或黃墻，在福建俗稱黃翅，在臺灣俗稱烏鯮或赤鰭仔[1-2]。黃鰭棘鯛肉質(zhì)鮮美、營養(yǎng)價值較高，口感極佳，是一種名貴的海洋經(jīng)濟(jì)魚類，也是我國海水魚類養(yǎng)殖的重要對象[3]。近年來，由于環(huán)境污染、過度捕撈、近親繁殖、人工放流和養(yǎng)殖逃逸等引起黃鰭棘鯛種質(zhì)資源混亂，自然種群嚴(yán)重衰退，導(dǎo)致種質(zhì)資源嚴(yán)重退化、抗病力減弱及養(yǎng)殖性能下降[4]。線粒體DNA(mtDNA)是細(xì)胞核外唯一具有自我復(fù)制功能的DNA，其進(jìn)化效率遠(yuǎn)高于單拷貝核基因，較少的mtDNA樣本即可代表1個群體[5]。mtDNA由于具有進(jìn)化速度快、母系遺傳、基因組較小、易于測序分析等特點，已成為分子進(jìn)化研究中的首選材料[6-9]，其中線粒體細(xì)胞色素b基因(Cytb)序列是mtDNA中最常用的分子標(biāo)記之一，已在魚類種群遺傳學(xué)和系統(tǒng)發(fā)育學(xué)中廣泛應(yīng)用[10-14]，但未見在華南沿海黃鰭棘鯛遺傳多樣性方面的研究報道。因此，基于線粒體Cytb基因序列分析華南沿海黃鰭棘鯛種群遺傳結(jié)構(gòu)，對促進(jìn)華南沿海黃鰭棘鯛資源的可持續(xù)利用及遺傳多樣性保護(hù)具有重要意義。【前人研究進(jìn)展】胡靜等[15]通過mtDNACOⅠ和Cytb基因序列分析發(fā)現(xiàn)，我國南海波紋唇魚的遺傳多樣性處于較低水平，存在遺傳分化現(xiàn)象但分化不顯著。李敏等[16]基于線粒體Cytb基因序列分析發(fā)現(xiàn)，我國近海和陸架的花斑蛇鯔遺傳分化不顯著，在漁業(yè)上可作為一個單元來管理。朱傳坤等[15]基于線粒體Cytb基因和D-loop序列分析發(fā)現(xiàn)，洪澤湖、千島湖和西湖3個華鳈群體存在明顯的遺傳分化現(xiàn)象。沈朕等[18]研究顯示，線粒體Cytb基因和D-loop區(qū)均可作為檢測大瀧六線魚遺傳多樣性的有效分子標(biāo)記。李大命等[19]分析太湖新銀魚Cytb和COⅠ基因序列的多態(tài)性發(fā)現(xiàn)，太湖新銀魚的遺傳多樣性處于較低水平，種群遺傳結(jié)構(gòu)不穩(wěn)定。近年來，國內(nèi)外關(guān)于黃鰭棘鯛線粒體遺傳學(xué)的研究主要有：劉紅艷等[20]對海南?？谑小⒏＝◤B門市和廣東珠海市3個地理群體海捕黃鰭鯛線粒體D-loop基因進(jìn)行多態(tài)性分析：Xia等[21]對我國水域黃鰭棘鯛mtDNA控制區(qū)序列變異遺傳進(jìn)行分析；江世貴等[22]對4種鯛科魚類Cytb基因序列及分子系統(tǒng)學(xué)進(jìn)行分析；徐田軍等[23]對舟山海域4種鯛科魚類線粒體Cytb基因全序列進(jìn)行克隆分析；黃元佳等[24]基于線粒體COⅠ基因的海陵島鯛科魚類系統(tǒng)發(fā)育分析；陳詠霞等[25]基于線粒體COⅠ基因序列的中國鯛科魚類系統(tǒng)進(jìn)化關(guān)系分析；Hsu等[26]基于線粒體和擴(kuò)增片段長度多態(tài)性的黃鰭棘鯛、黑棘鯛(Acanthopagrusschlegelii)、桔鰭棘鯛(A.sivicolus)、臺灣棘鯛(A.taiwanensis)、琉球黃鰭棘鯛(A.chinshira)和太平洋棘鯛(A.pacificus)區(qū)分。【本研究切入點】為了保護(hù)和恢復(fù)黃鰭棘鯛海洋資源，我國華南沿海各地常年大量放流黃鰭棘鯛魚苗，但目前未見針對黃鰭棘鯛遺傳結(jié)構(gòu)變化長期跟蹤及進(jìn)行遺傳多樣性分析的相關(guān)研究報道?！緮M解決的關(guān)鍵問題】比較我國華南沿海黃鰭棘鯛8個地理群體的線粒體Cytb基因序列，分析其特征和遺傳多樣性背景，了解我國黃鰭棘鯛群體的遺傳差異和分化情況，并對其遺傳多樣性進(jìn)行評價，為黃鰭棘鯛種質(zhì)資源的保護(hù)和利用提供參考依據(jù)。

1 材料與方法

1.1 試驗材料

供試黃鰭棘鯛樣本于2017年12月至2018年12月分別采自我國華南沿海的福建廈門(以下簡稱XM)、廣東汕尾(以下簡稱SW)、廣東陽江(以下簡稱YJ)、海南?？?以下簡稱HK)、海南三亞(以下簡稱SY)、廣西北海(以下簡稱BH)、廣西欽州(以下簡稱QZ)和廣西防城港東興(以下簡稱FC)，共320尾。每尾約取5.0 g肌肉組織存放于95 %乙醇中，-20 ℃保存?zhèn)溆?。以瓊州海峽為界，又將黃鰭棘鯛群體劃分為東海域(EAST)組群和西海域(WEST)組群。

1.2 DNA提取、PCR擴(kuò)增與測序

每尾黃鰭棘鯛分別剪取肌肉組織50.0 mg，采用醋酸銨法[27]提取總DNA，用微量核酸蛋白分析儀測定樣品DNA濃度和純度，然后以1 %瓊脂糖凝膠電泳檢測樣品DNA質(zhì)量，檢測合格的樣品DNA置于-20 ℃保存?zhèn)溆?。以合格的樣品DNA為模板進(jìn)行PCR擴(kuò)增，PCR引物F(5′-GACTTGAAAAACCACCGTTG-3′)和R(5′-GGCGTTCGGCTTACAAAA-3′)委托生工生物工程(上海)股份有限公司合成。PCR反應(yīng)體系50.0 μl：2×ExTaqMastermix Bufffer 25.0 μl，10 μmol/L上、下游引物各2.0 μl，100 ng/μl DNA模板2.0 μl，ddH2O補(bǔ)足至50.0 μl；擴(kuò)增程序：95 ℃預(yù)變性5 min；95 ℃ 30 s，62 ℃ 30 s，72 ℃ 1 min，進(jìn)行40個循環(huán)；72 ℃延伸10 min，4 ℃保存。擴(kuò)增產(chǎn)物經(jīng)1 %瓊脂糖凝膠電泳檢測合格后送至深圳華大基因科技有限公司測序。

1.3 統(tǒng)計分析

利用Chromas 2.33(http://www.technelysium com.au/chromas.html)對測序結(jié)果進(jìn)行編輯、校對和剪切；利用SeqMan 7.1.0(https://www.dnastar.com/)進(jìn)行拼接，以獲得Cytb基因全長序列；采用MEGA X(https://www.megasoftware.net/)分析序列的堿基組成及變異位點。運(yùn)用DnaSP 5.10(http://www.ub.edu/dnasp/)分析黃鰭棘鯛的單倍型(h)、單倍型多樣性指數(shù)(Hd)、核苷酸遺傳多樣性指數(shù)(π)；利用MEGA X分析序列計算Kimura 2-parameter遺傳距離和分析遺傳結(jié)構(gòu)，然后采用遺傳距離、基因流(Nm)、遺傳分化指數(shù)(Fst)評價兩兩群體間的遺傳差異，重復(fù)抽樣10000次以檢驗Fst結(jié)果的差異顯著性，并構(gòu)建非加權(quán)組平均法(UPGMA)系統(tǒng)發(fā)育進(jìn)化樹，以及進(jìn)行群體間聚類分析。通過種群分化測試(Exact test of population differentiation)進(jìn)行隨機(jī)交配群(Population panmixia)假設(shè)即單倍型在群體間隨機(jī)分布的檢驗，以AMOVA分析估算遺傳變異的分布，10 000次隨機(jī)重復(fù)抽樣后進(jìn)行差異顯著性分析。采用Tajima′sD檢驗[28]和Fu′sFs檢驗[29]檢測種群進(jìn)化是否嚴(yán)格遵循中性理論。根據(jù)公式t=Tτ/2u[30]計算群體自擴(kuò)張開始至今的時間(t)，其中T為世代時間，τ為擴(kuò)張時間參數(shù)(在Arlequin軟件中根據(jù)核苷酸錯配分布計算獲得)，u為序列的突變速率[u=2μk，式中，k為序列長度，μ為每個核苷酸位點的突變速率(每百萬年突變1.0 %～2.5 %[31])]。

2 結(jié)果與分析

2.1 線粒體Cytb基因序列分析結(jié)果

經(jīng)PCR擴(kuò)增、測序和校對，獲得8個黃鰭棘鯛群體的320條線粒體Cytb基因全長序列(1141 bp)，Cytb基因序列中堿基T、C、A和G的平均含量分別為30.3 %、29.4 %、24.8 %和15.5 %，且A+T平均含量(55.1 %)高于C+G平均含量(44.9 %)，表現(xiàn)出較強(qiáng)的少G偏倚性，符合脊椎動物線粒體DNA序列特征。在Cytb基因序列片段使用的密碼子中，T、C、A和G在第1、2和3位點的平均含量存在明顯差異，其中G平均含量變化最明顯，在第1、2和3位點分別為26.4 %、13.9 %和6.0 %，說明密碼子的堿基使用頻率存在明顯偏向性。

在黃鰭棘鯛Cytb基因全長序列中未檢測到插入或缺失位點，但檢測到63個變異位點，占總位點數(shù)的5.52 %，其中簡約信息位點25個、單一變異位點38個；在XM群體中發(fā)現(xiàn)變異位點最多，為23個，在FC群體中發(fā)現(xiàn)變異位點最少，僅10個。在63個變異位點中，59個為轉(zhuǎn)換、4個為顛換，轉(zhuǎn)換與顛換比(R值)為14.75，說明線粒體Cytb基因序列突變未達(dá)飽和狀態(tài)，受進(jìn)化噪音影響的可能性較小。其中，52個變異位點發(fā)生在密碼子第3位上，10個變異位點發(fā)生在密碼子第1位上，1個變異位點發(fā)生在密碼子第2位上。

黃鰭棘鯛Cytb基因部分序列編碼380個氨基酸殘基，其中亮氨酸平均含量最高，半胱氨酸平均含量最低。4種堿基(T、C、A和G)在氨基酸密碼子3個位點出現(xiàn)頻率差異的趨勢一致，均為在密碼子第1位點無明顯的堿基偏倚性，在第2位點存在明顯的T偏倚和反G偏倚，在第3位點存在反G偏倚。將黃鰭棘鯛群體Cytb基因的核苷酸序列翻譯成氨基酸序列，檢測到8個氨基酸變異位點，其中7個由密碼子第1位突變引起，另1個由密碼子第2位突變引起，其他核苷酸突變未引起氨基酸變異，屬同義突變，據(jù)此定義9個氨基酸為單倍型，說明黃鰭棘鯛的Cytb氨基酸序列具有較低的多態(tài)性。

2.2 單倍型及遺傳多樣性分析結(jié)果

8個黃鰭棘鯛群體的320條線粒體Cytb基因序列定義62個單倍型，其中，42個單倍型僅在一個個體中檢測到，在其余20個單倍型中，有18個是群體間共享單倍型，另外2個單倍型在不止1個個體中發(fā)現(xiàn)，但這些個體僅分布在同1個群體。Hap1為優(yōu)勢單倍型，在8個群體中均有出現(xiàn)；單倍型Hap4和Hap2分別由8和7個群體共享；單倍型Hap5、Hap6、Hap8、Hap11、Hap16和Hap37均由3個以上群體共享；單倍型Hap3、Hap7、Hap9、Hap15、Hap23、Hap26、Hap34、Hap40和Hap54均由2個以上群體共享；而Hap10為BH群體獨(dú)有，Hap12～Hap14為FC群體獨(dú)有，Hap17～Hap22、Hap24、Hap25和Hap27為HK群體獨(dú)有，Hap28～Hap32為QZ群體獨(dú)有，Hap33、Hap35和Hap36為SW群體獨(dú)有，Hap38、Hap39和Hap41～Hap44為SY群體獨(dú)有，Hap46～Hap53和Hap55～Hap59為XM群體獨(dú)有，Hap60～Hap62為YJ群體獨(dú)有。在62個單倍型中，7個單倍型為2個組群共享，13個單倍型為WEST組群獨(dú)享，42個單倍型為EAST組群獨(dú)享。62個單倍型的出現(xiàn)頻率存在較明顯差異，其中以Hap1出現(xiàn)的頻率最高，占總個體數(shù)的38.75 %，其次為Hap4和Hap2，分別占總個體數(shù)的23.75 %和8.44 %。說明2個組群共享、WEST組群獨(dú)享和EAST組群獨(dú)享的單倍型是華南沿海黃鰭棘鯛在長期進(jìn)化過程中形成的較穩(wěn)定的優(yōu)勢基因型。

利用DnaSP 5.10分析8個黃鰭棘鯛群體的單倍型多樣性和核苷酸多樣性，結(jié)果(表1)表明，8個黃鰭棘鯛地理群體的Hd在0.57179～0.83718，其中以XM群體最高，為0.83718，其次為BH群體和HK群體，分別為0.74872和0.72821，QZ群體最低，為0.57179；π在0.00089～0.00162，其中BH群體最高，為0.00162，其次為XM群體和HK群體，分別為0.00152和0.00148，SW群體最低，為0.00089。WEST組群和EAST組群的Hd分別為0.66471和0.69714，π分別為0.00128和0.00120，二者差異不明顯。將所有樣本作為一個整體進(jìn)行分析，發(fā)現(xiàn)黃鰭棘鯛群體的Hd為0.78638，π為0.00164。

2.3 群體系統(tǒng)進(jìn)化關(guān)系分析結(jié)果

根據(jù)黃鰭棘鯛Cytb基因全長序列，以黑棘鯛線粒體Cytb基因序列(GenBank登錄號NC_018553.1)作為外群，利用MEGA X分析序列計算Kimura 2-parameter的遺傳距離并構(gòu)建單倍型非加權(quán)組平均法(UPGMA)系統(tǒng)發(fā)育進(jìn)化樹及進(jìn)行群體間聚類分析，結(jié)果(圖1)顯示，62個單倍型系統(tǒng)發(fā)育進(jìn)化樹可分為三大分支：Group A、Group B和Group C，各節(jié)點的支持率均較低(<30 %)，只有極少數(shù)兩兩單倍型分支間的支持率超過50 %；不同群體的單倍型未按取樣地點有規(guī)律地分布在系統(tǒng)發(fā)育進(jìn)化樹上，但大部分單倍型按瓊州海峽東、西海域分布在系統(tǒng)發(fā)育進(jìn)化樹上，即分為EAST和WEST兩個組群。聚類分析結(jié)果(圖2)也表明，8個黃鰭棘鯛地理群體以瓊州海峽分隔為瓊州海峽東和瓊州海峽西兩大部分。Hap1和Hap4是高頻單倍型，二者在所有群體中均有出現(xiàn)，以Hap4為主體的Group B包含了WEST組群(FC、QZ和BH群體)絕大部分個體(87.74 %)，以Hap1為中心的Group A則涵蓋了EAST組群(HK、SY、YJ、SW和XM群體)絕大部分個體(87.38 %)。

2.4 群體遺傳結(jié)構(gòu)分析結(jié)果

利用MEGA X分析測序獲得320尾黃鰭棘鯛樣品的線粒體Cytb基因序列，得到8個黃鰭棘鯛地理群體間的遺傳距離(表2)。由表2可知，群體內(nèi)遺傳距離為0.00089～0.00162，各地理群體內(nèi)的遺傳距離排序為BH>XM>HK>FC>SY>YJ>QZ>SW；從群體間遺傳距離來看，QZ群體與XM群體的遺傳距離最遠(yuǎn)，為0.00232，而遺傳距離最近的為SW群體與YJ群體(0.00093)；8個黃鰭棘鯛地理群體間的Nm為0.39049～無窮大(inf)，WEST組群內(nèi)兩兩間的Nm為8.77670～inf，EAST組群內(nèi)兩兩群體間的Nm均為inf，WEST組群與EAST組群兩兩群體間的Nm為0.39049～1.36673，BH群體與EAST組群間的Nm為1.00378～1.36673，WEST組群與EAST組群的Nm以BH群體為過渡。說明EAST組群和WEST組群內(nèi)基因交流頻繁，組群間由于地理隔離而限制了群體的基因交流。

由表3可知，8個黃鰭棘鯛地理群體間的Fst在-0.01300～0.56149，其中，SW群體與QZ群體的Fst最大(0.56149)，且差異極顯著(P<0.01，下同)，說明其遺傳分化程度最大，而SY群體與YJ群體間

表1 黃鰭棘鯛的單倍型及遺傳多樣性信息

的Fst(最小-0.01300)，說明其遺傳分化程度最小。在EAST組群內(nèi)，F(xiàn)st均為負(fù)值，但相互間差異不顯著(P>0.05，下同)；而在WEST組群內(nèi)，最大分化出現(xiàn)在BH群體與QZ群體間(Fst=0.05390，P=0.00901)，且差異顯著(P<0.05，下同)；而WEST組群與EAST組群兩兩群體間的Fst為0.26785～0.56149，且差異極顯著。相對于組群內(nèi)，WEST組群與EAST組群的兩兩群體間遺傳距離和遺傳分化指數(shù)較大，基因流較小，說明WEST組群與EAST組群存在明顯的遺傳分化。綜合8個黃鰭棘鯛地理群體的聚類分析結(jié)果及組群間的Fst分析結(jié)果可知，瓊州海峽與海南島形成的天然地理屏障限制了8個黃鰭棘鯛地理群體的基因交流，且WEST組群與EAST組群存在明顯的遺傳分化，在漁業(yè)管理上可將瓊州海峽以東和瓊州海峽以西的黃鰭棘鯛群體作為2個管理單位(MUs)進(jìn)行管理。

對以樣品來源劃分的EAST組群(HK、SY、YJ、SW和XM群體)和WEST組群(FC、QZ和BH群體)黃鰭棘鯛群體，通過AMOVA計算其遺傳變異分布，結(jié)果(表4)顯示組群間、組群內(nèi)群體間和群體內(nèi)個體間的變異分別為41.26 %、0.09 %和58.65 %，經(jīng)檢驗均達(dá)顯著或極顯著差異水平，進(jìn)一步說明華南沿海黃鰭棘鯛EAST組群和WEST組群內(nèi)群體間分化較小，組群間由于地理隔離限制了組群間的基因交流，而WEST組群和EAST組群內(nèi)的群體基因交流因不受限制而無群體分化。將所有黃鰭棘鯛個體按單倍型類群可分成三部分(Group A、Group B和Group C)進(jìn)行AMOVA分析，結(jié)果(表4)顯示70.98 %變異來自單倍型類群間。

Exact檢驗結(jié)果(表3)表明，單倍型在群體間的隨機(jī)交配假設(shè)檢驗水平不均勻，其在EAST組群與WEST組群間表現(xiàn)為顯著(P=0.00000)，不符合單倍型在群體間群體是隨機(jī)分布的假設(shè)，而在WEST組群內(nèi)(P=0.49038)和EAST組群內(nèi)(P=0.70135)表現(xiàn)不顯著，符合單倍型分別在WEST組群和EAST組群內(nèi)群體是隨機(jī)分布的假設(shè)。

2.5 群體動態(tài)分析結(jié)果

由表5可知，單倍型類群Group C的Tajima’sD和Fu’sFs檢驗結(jié)果均為負(fù)值，但未達(dá)顯著差異水平，而個體數(shù)占絕大部分的單倍型類群Group A和Group B的Tajima’sD和Fu’sFs檢驗結(jié)果為負(fù)值，但均存在顯著差異，單倍型類群總體的Tajima’sD和Fu’sFs檢驗結(jié)果也存在顯著差異，說明序列中含有比中性進(jìn)化模型更多的核苷酸位點變化，預(yù)示被研究群體可能經(jīng)歷了一個擴(kuò)張歷史。核苷酸錯配分析結(jié)果(圖3)表明，單倍型類群Group A和Group B均為單峰分布，Group C為雙峰分布，但單倍型類群總體上呈雙峰分布，可能是多個單倍型類群疊加的結(jié)果?？梢?，Tajima’sD和Fu’sFs檢驗及核苷酸錯配分析結(jié)果均表明華南沿海黃鰭棘鯛群體可能經(jīng)歷了種群擴(kuò)張事件。采用Arlequin分析獲得擴(kuò)張時間參數(shù)τ=2.42773，以每百萬年突變2.0 %為線粒體Cytb基因的核苷酸進(jìn)化速率，根據(jù)黃鰭棘鯛存在1～2齡雄性性腺發(fā)育成熟后2～3齡轉(zhuǎn)變?yōu)榇菩孕阅孓D(zhuǎn)現(xiàn)象，本研究參考施曉峰等[32]、洪萬樹等[33]的方法以生物學(xué)最小齡2年齡作為一個世代的時間跨度，并根據(jù)公式τ=2ut，計算得到華南沿海黃鰭棘鯛群體擴(kuò)張時間約為5.319萬年前的更新世晚期。綜上所述，采用無限突變位點模型Tajima’sD和Fu’sFs、線粒體DNA單倍型錯配分布等方法對8個黃鰭棘鯛地理種群的歷史變動情況進(jìn)行調(diào)查，可判斷華南沿海黃鰭棘鯛群體可能經(jīng)歷了種群擴(kuò)張事件，擴(kuò)張時間約在5.319萬年前。

表2 8個黃鰭棘鯛地理群體間的遺傳距離(對角線下方)、群體內(nèi)遺傳距離(對角線上)和Nm(對角線上方)比較

表3 8個黃鰭棘鯛地理群體間的Fst(對角線下方)及其Exact 檢驗(對角線上方)結(jié)果

表5 Tajima’s D和Fu’s Fs檢驗結(jié)果

3 討 論

遺傳多樣性是物種長期進(jìn)化或群體持續(xù)生存并適應(yīng)不斷變化環(huán)境而進(jìn)化的結(jié)果，物種遺傳多樣性或變異性越豐富，表明物種的進(jìn)化潛力越大，對環(huán)境改變響應(yīng)的進(jìn)化能力越強(qiáng)[34-35]。核苷酸多樣性指數(shù)和線粒體DNA的單倍型多樣性指數(shù)是衡量DNA多態(tài)程度的兩個重要指標(biāo)[36]。本研究中，8個黃鰭棘鯛地理種群基于線粒體Cytb基因序列分析的單倍型多樣性指數(shù)在0.57179～0.83718，群體的單倍型多樣性指數(shù)為0.78638，核苷酸多樣性指數(shù)在0.00089～0.00162，總體為0.00164。相對于Grant和Bowen[37]的分類結(jié)果，本研究中8個黃鰭棘鯛群體基于線粒體Cytb基因全長序列的遺傳多樣性符合高Hd低π類型，具有較高的Hd(>0.5)和較低的π(<0.005)可能是種群歷史上經(jīng)歷瓶頸效應(yīng)后快速擴(kuò)張的表現(xiàn)；與同樣基于線粒體Cytb基因序列標(biāo)記分析的同海域其他魚類相比，8個黃鰭棘鯛群體的多樣性低于花斑蛇鯔[16]、短尾大眼鯛[38]、大黃魚[39]和黃鰭馬面鲀[40]，高于龍頭魚[41]，與波紋唇魚[15]和藍(lán)圓鲹[42]接近?？梢?，基于線粒體Cytb基因的黃鰭棘鯛多樣性處于中等偏低水平。魚類的遺傳多樣性受捕撈壓力影響，過度捕撈會顯著降低其遺傳多樣性[43]。因此，在漁業(yè)資源開發(fā)利用中。正確認(rèn)識和評價魚類的遺傳多樣性水平，是可持續(xù)開發(fā)利用該漁業(yè)資源的基礎(chǔ)。

Tajima[28]研究認(rèn)為，應(yīng)用Tajima’sD和Fu’sFs檢驗推測種群歷史時，若Tajima’sD與Fu’sFs檢驗結(jié)果呈負(fù)值，且在統(tǒng)計學(xué)上達(dá)顯著水平，說明序列中含有比中性進(jìn)化模型更多的核苷酸位點變化，預(yù)示著種群可能曾經(jīng)經(jīng)歷一個擴(kuò)張歷史。本研究中，單倍型類群總體的Tajima’sD和Fu’sFs檢驗結(jié)果均為顯著負(fù)值，且個體數(shù)占絕大部分的單倍型類群Group A和Group B的核苷酸錯配分布圖為單峰，說明華南沿海黃鰭棘鯛群體發(fā)生了群體擴(kuò)張，并依據(jù)Arlequin分析所得擴(kuò)張時間參數(shù)推算擴(kuò)張時間約在5.319萬年前的更新世晚期。

單倍型間的遺傳距離是衡量一個物種或群體線粒體變異程度的重要指標(biāo)，也是衡量群體間遺傳變異程度的可靠參數(shù)，遺傳距離越大表明群體間親緣關(guān)系越遠(yuǎn)。本研究利用MEGA X根據(jù)線粒體Cytb基因序列計算8個黃鰭棘鯛群體間的遺傳距離為0.00093～0.00232，群體內(nèi)的遺傳距離為0.00089～0.00162，群體間和群體內(nèi)的遺傳距離較小，均遠(yuǎn)低于0.05000，且群體間遺傳距離與群體內(nèi)遺傳距離無明顯差異，說明8個黃鰭棘鯛群體間存在較近的親緣關(guān)系，與Shaklee[44]對魚類種群、種和屬水平上的分類依據(jù)吻合；UPGMA聚類分析結(jié)果表明，大部分單倍型按瓊州海峽東、西海域分布在系統(tǒng)發(fā)育進(jìn)化樹上，Group A涵蓋了EAST組群絕大部分個體(87.38 %)，Group B包含了WEST組群絕大部分個體(87.74 %)；群體聚類分析結(jié)果也表明，8個黃鰭棘鯛地理群體以瓊州海峽分隔為兩部分。一般認(rèn)為，群體間的Nm<1.00000表示群體可能由于遺傳漂變而發(fā)生分化，而Nm>1.00000表示群體間的Nm水平較高，群體間遺傳分化較??；當(dāng)Nm>4.00000時，種群間的基因交流更充分，遺傳分化更小[45]。本研究中EAST組群和WEST組群內(nèi)的Nm均>4.00000，表明組群內(nèi)的基因交流充分，遺傳分化更小；而EAST組群與WEST組群間，以BH群體為過渡，其余Nm均小于1.00000，表明組群間可能發(fā)生遺傳分化。Wright[46]提出的遺傳分化標(biāo)準(zhǔn)為：Fst≤0.05，無遺傳分化；0.05

4 結(jié) 論

以瓊州海峽為界分隔的東海域和西海域兩個華南沿海黃鰭棘鯛組群(EAST和WEST)間的Fst為0.26785～0.56149，遺傳分化明顯。因此，可將我國華南沿海黃鰭棘鯛作為2個管理單位(MUs)進(jìn)行漁業(yè)管理。