基于線粒體COI基因的江蘇湖泊鳙群體遺傳多樣性分析

摘 要：為掌握江蘇湖泊鳙群體的遺傳資源狀況，利用COⅠ基因序列研究了6個(gè)湖泊（太湖、滆湖、長(zhǎng)蕩湖、高郵湖、洪澤湖和駱馬湖）鳙群體的遺傳多樣性及遺傳結(jié)構(gòu)。結(jié)果顯示，在分析的630 bp COⅠ基因片段中，堿基A+T含量（55.3%）高于C+G含量（44.7%）。223條COⅠ序列檢測(cè)出14個(gè)變異位點(diǎn)，定義9個(gè)單倍型，6個(gè)群體的單倍型多樣性為0.494～0.665，核苷酸多樣性為0.000 9～0.001 9，總體單倍型多樣性和核苷酸多樣性分別為0.576和0.001 2，具有高單倍型多樣性和低核苷酸多樣性的特點(diǎn)，表明6個(gè)湖泊鳙群體遺傳多樣性較低。分子方差分析結(jié)果顯示，遺傳變異主要來(lái)自于群體內(nèi)（99.39%），總遺傳分化指數(shù)為0.006 1，群體間遺傳分化指數(shù)為-0.022 8～0.044 5，遺傳距離為0.001 2～0.002 7，表明群體間沒(méi)有明顯的遺傳分化。單倍型系統(tǒng)發(fā)育樹(shù)和網(wǎng)絡(luò)結(jié)構(gòu)圖顯示，單倍型在6個(gè)群體中混雜分布，未形成特定的湖泊分支。中性檢驗(yàn)和歧點(diǎn)分布分析表明，整個(gè)鳙群體經(jīng)歷過(guò)顯著的種群擴(kuò)張。結(jié)果表明，江蘇湖泊鳙群體的遺傳多樣性較低，遺傳結(jié)構(gòu)趨于同質(zhì)化，需要采取措施提高鳙的遺傳多樣性。

關(guān)鍵詞：鳙；COⅠ基因；遺傳多樣性；遺傳結(jié)構(gòu)；增殖放流

鳙（Aristichthys nobilis）俗稱花鰱、胖頭魚(yú)，隸屬于鯉形目（Cypriniformes）、鯉科（Cyprinidae）、鳙屬（Aristichthys），廣泛分布于全國(guó)各大江河及其附屬湖泊［1］。鳙是我國(guó)“四大家魚(yú)”之一，具有生長(zhǎng)速度快、營(yíng)養(yǎng)豐富、經(jīng)濟(jì)價(jià)值高、養(yǎng)殖技術(shù)成熟等優(yōu)勢(shì)，已成為我國(guó)重要的大宗淡水養(yǎng)殖魚(yú)類。鳙為濾食性魚(yú)類，主要攝食浮游動(dòng)物，在修復(fù)水體環(huán)境及維持水域生態(tài)系統(tǒng)平衡穩(wěn)定中具有重要作用［2-4］。鳙屬于江湖洄游型魚(yú)類，不能在湖泊中繁殖，加之江湖阻隔、環(huán)境污染、過(guò)度捕撈等不利因素的影響，湖泊鳙天然資源嚴(yán)重衰退，個(gè)體低齡化、小型化現(xiàn)象嚴(yán)重［1，5-8］。為增加湖泊鳙資源量及修復(fù)湖泊生態(tài)環(huán)境，多年來(lái)我國(guó)開(kāi)展了大規(guī)模的增殖放流活動(dòng)，這對(duì)于恢復(fù)湖泊鳙資源量起到了重要作用，同時(shí)也可能會(huì)對(duì)湖泊鳙的種質(zhì)資源遺傳多樣性及遺傳結(jié)構(gòu)產(chǎn)生影響［9-11］。為了避免人工放流群體對(duì)野生群體的負(fù)面影響，應(yīng)在放流之前評(píng)估野生群體的遺傳特征［12］。

遺傳多樣性是生物多樣性的重要組成部分，也是生物適應(yīng)環(huán)境的基礎(chǔ)和物種進(jìn)化的動(dòng)力，掌握生物的遺傳多樣性對(duì)于科學(xué)保護(hù)和合理利用生物資源具有重要指導(dǎo)作用。線粒體DNA（mitochondrial DNA，mtDNA）是研究物種群體遺傳多樣性常用的分子標(biāo)記，具有結(jié)構(gòu)簡(jiǎn)單、母系遺傳、不易發(fā)生基因重組且進(jìn)化速度快等諸多特點(diǎn)［13］。其中，細(xì)胞色素氧化酶Ⅰ亞基（cytochrome oxidase subunit Ⅰ，COⅠ）基因的結(jié)構(gòu)和功能較為清楚，進(jìn)化速率較快，且有通用引物便于擴(kuò)增和測(cè)序，因此被廣泛應(yīng)用于水生動(dòng)物尤其是魚(yú)類的遺傳多樣性和遺傳結(jié)構(gòu)研究［14-16］。

長(zhǎng)江中下游區(qū)域是我國(guó)湖泊集中分布區(qū)，也是我國(guó)規(guī)劃的內(nèi)陸6個(gè)主要水生生物增殖放流區(qū)域之一。鳙在長(zhǎng)江中下游區(qū)域的放流數(shù)量居第2位，其功能定位主要為漁民增收及生物凈化［9］。本研究擬采用COⅠ基因序列作為分子標(biāo)記，探究增殖放流對(duì)江蘇省6個(gè)湖泊（太湖、滆湖、長(zhǎng)蕩湖、高郵湖、洪澤湖和駱馬湖）鳙群體遺傳多樣性和遺傳結(jié)構(gòu)的影響，以期為科學(xué)保護(hù)和合理利用鳙種質(zhì)資源及規(guī)范人工放流提供參考依據(jù)。

1 材料和方法

1.1 樣本采集和處理

2020—2022年在長(zhǎng)江流域的太湖、滆湖和長(zhǎng)蕩湖及淮河流域的高郵湖、洪澤湖和駱馬湖開(kāi)展?jié)O業(yè)資源監(jiān)測(cè)，共采集鳙樣本223尾，現(xiàn)場(chǎng)進(jìn)行種類鑒定及生物學(xué)測(cè)量，采樣水域及各群體樣本數(shù)量見(jiàn)表1。采集樣本后，剪取鳙的尾鰭組織，放入2 mL離心管，加入無(wú)水乙醇保存于4 ℃冰箱。

1.2 基因組提取、PCR擴(kuò)增和測(cè)序

使用DNA試劑盒（TaKaRa公司）提取鳙的基因組DNA，將DNA溶于ddH2O，經(jīng)1%瓊脂糖凝膠電泳檢測(cè)后，保存于-20 ℃冰箱備用。

通過(guò)引物FishF1（5′-TCAACCAACCACAAAGACATTGGCAC-3′）和FishR1（5′-TAGACTTCTGGGTGGCCAAAGAATCA-3′）擴(kuò)增COⅠ基因序列［17］。PCR反應(yīng)體系共50 μL，其中包括：2×Premix TaqTM 25 μL（TaKaRa Taq 1.25 U/25 μL、0.2 mmol/L dNTP、1.5 mmol/L Mg2+）、上下游引物各2 μL（10 μmol/L）、DNA模板2 μL和ddH2O 19 μL。PCR擴(kuò)增程序?yàn)椋?4 ℃預(yù)變性4 min；94 ℃變性40 s，55 ℃復(fù)性40 s，72 ℃延伸50 s，35個(gè)循環(huán)；72 ℃延伸5 min，4 ℃保存。PCR產(chǎn)物經(jīng)1.5%瓊脂糖凝膠電泳檢測(cè)，合格的PCR產(chǎn)物送生工生物工程（上海）股份有限公司進(jìn)行雙向測(cè)序。

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

通過(guò)BioEdit 7.0.9.0軟件［18］讀取序列，并對(duì)序列進(jìn)行編輯和校對(duì)。使用CLUSTAL X2.0軟件［19］進(jìn)行多重比對(duì)和排序，獲取長(zhǎng)度一致性序列。通過(guò)DanSP 5.0軟件［20］計(jì)算群體的單倍型數(shù)量、單倍型多樣性（haplotype diversity，Hd）、核苷酸多樣性（nucleotide diversity，Pi）、變異位點(diǎn)（variable sites）、單一信息位點(diǎn)（singleton sites）、簡(jiǎn)約信息位點(diǎn)（parsimony-informative sites）等。采用MEGA 7.0.26軟件［21］統(tǒng)計(jì)序列的堿基組成，計(jì)算群體間的遺傳距離。以鰱為外類群，采用鄰接法（neighbor-joining，NJ）構(gòu)建單倍型系統(tǒng)發(fā)育樹(shù)，系統(tǒng)樹(shù)中節(jié)點(diǎn)的自舉置信水平應(yīng)用自引導(dǎo)估計(jì)（循環(huán)次數(shù)為1 000次）。采用Network 4.5.1軟件［22］構(gòu)建單倍型網(wǎng)絡(luò)結(jié)構(gòu)圖，檢測(cè)各個(gè)單倍型間的進(jìn)化關(guān)系。使用Arlequin 3.5軟件［23］進(jìn)行分子方差分析（analysis of molecular variance，AMOVA），檢測(cè)遺傳變異來(lái)源，計(jì)算群體間遺傳分化指數(shù)（Fst）。通過(guò)Tajima’s D、Fu’s Fs中性檢驗(yàn)和錯(cuò)配分布分析推測(cè)鳙群體的歷史動(dòng)態(tài)。

2 結(jié)果

2.1 COⅠ序列變異及遺傳多樣性

本研究獲得了223尾鳙的線粒體COⅠ基因序列，其長(zhǎng)度為630 bp，序列中無(wú)堿基的插入和缺失。223條COⅠ序列的平均堿基組成為A（26.6%），T（28.7%），C（27.2%），G（17.5%）。A+T的含量（55.3%）高于C+G的含量（44.7%），堿基組成表現(xiàn)出明顯的偏倚。在223條序列中，共有14個(gè)變異位點(diǎn)，占總位點(diǎn)數(shù)的2.22%，其中單突變位點(diǎn)有9個(gè)，簡(jiǎn)約信息位點(diǎn)有5個(gè)。

6個(gè)鳙群體的單倍型多樣性（Hd）為0.494～0.665，核苷酸多樣性（Pi）為0.000 9～0.001 9，其中太湖群體的遺傳多樣性最高，長(zhǎng)蕩湖群體的遺傳多樣性最低。6個(gè)群體的總體單倍型多樣性和核苷酸多樣性分別為0.576和0.001 2，呈現(xiàn)出高Hd和低Pi模式（見(jiàn)表1）。

2.2 群體單倍型組成及系統(tǒng)發(fā)育

6個(gè)鳙群體定義9個(gè)單倍型（Hap 1～Hap 9），其中單倍型Hap 1、Hap 2和Hap 7為6個(gè)群體的共享單倍型，個(gè)體數(shù)量有206個(gè)；單倍型Hap 3、Hap 4和Hap 6為部分群體所共享，個(gè)體數(shù)量有14個(gè)；單倍型Hap 5、Hap 8和Hap 9為太湖群體所獨(dú)享，個(gè)體數(shù)量有3個(gè)（見(jiàn)表2）。

單倍型系統(tǒng)發(fā)育樹(shù)顯示，單倍型Hap 1～Hap 8聚為一支，單倍型Hap 9聚為一支，單倍型隨機(jī)分布在6個(gè)群體中，并未形成相應(yīng)的地理分支（見(jiàn)圖1）。單倍型網(wǎng)絡(luò)結(jié)構(gòu)圖呈星型（見(jiàn)圖2），也沒(méi)有形成特定的地理譜系結(jié)構(gòu)，進(jìn)一步支持了單倍型系統(tǒng)發(fā)育樹(shù)的結(jié)果。單倍型Hap 2位于中心位置，且個(gè)體數(shù)量最多，推測(cè)為祖先單倍型。

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

分子方差結(jié)果顯示，6個(gè)群體中，群體內(nèi)的遺傳變異占總體變異的99.39%，群體間的遺傳變異為0.61%，6個(gè)群體間遺傳分化指數(shù)為0.006 1（Fstlt;0.05），表明群體間的遺傳變異主要來(lái)自群體內(nèi)個(gè)體間（見(jiàn)表3）。

兩兩群體間的遺傳分化指數(shù)為-0.022 8～0.044 5（Fstlt;0.05），遺傳距離為0.001 2～0.002 7（見(jiàn)表4），表明群體間親緣關(guān)系較近，沒(méi)有產(chǎn)生明顯的遺傳分化。

2.4 群體歷史動(dòng)態(tài)

中性檢驗(yàn)結(jié)果顯示，太湖、長(zhǎng)蕩湖、洪澤湖和駱馬湖的Tajima’s D值和Fu’s Fs值均為負(fù)值，只有太湖群體的統(tǒng)計(jì)檢驗(yàn)具有顯著性差異，其他群體的統(tǒng)計(jì)檢驗(yàn)沒(méi)有顯著性差異（見(jiàn)表5）。滆湖和高郵湖群體的Tajima’s D值和Fu’s Fs值為正值，統(tǒng)計(jì)檢驗(yàn)也沒(méi)有顯著性差異。整個(gè)鳙群體的Tajima’s D值和Fu’s Fs值為負(fù)值，統(tǒng)計(jì)檢驗(yàn)具有顯著性差異，且歧點(diǎn)分布曲線呈明顯的單峰型（見(jiàn)圖3），表明整個(gè)鳙群體經(jīng)歷過(guò)種群擴(kuò)張。

3 討論

3.1 鳙群體的遺傳多樣性

單倍型多樣性（Hd）和核苷酸多樣性（Pi）是評(píng)價(jià)物種遺傳多樣性的兩個(gè)重要指標(biāo)，其數(shù)值越大，表明群體遺傳變異越豐富［24］。核苷酸多樣性考慮了各種單倍型在群體中所占的比例，更能精確反映群體的多態(tài)程度［25］。本研究基于COⅠ序列分析了江蘇6個(gè)湖泊鳙群體的遺傳多樣性，得出單倍型多樣性為0.494～0.665，核苷酸多樣性為0.000 9～0.001 9，總體單倍型和核苷酸多樣性分別為0.576和0.001 2。根據(jù)Grant等［26］提出的單倍型多樣性和核苷酸多樣性大小標(biāo)準(zhǔn)，太湖、滆湖、高郵湖和洪澤湖群體屬于高Hd和低Pi模式，長(zhǎng)蕩湖和駱馬湖群體屬于低Hd和低Pi模式，6個(gè)群體的核苷酸多樣性均明顯小于0.005，表明鳙群體的遺傳多樣性較低。與鳙相似，湖泊水庫(kù)鰱的遺傳多樣性也較低［12，27-28］，這是因?yàn)楹?kù)中的鰱、鳙不能自然繁殖，其資源量主要來(lái)源于人工放流群體，由于放流群體的親本數(shù)量有限及人工選擇作用，更容易發(fā)生瓶頸效應(yīng)和近交衰退，造成一些稀有等位基因的丟失，進(jìn)而導(dǎo)致放流群體的遺傳多樣性普遍低于野生群體［29-30］。比較而言，長(zhǎng)江鳙群體的遺傳多樣性比較豐富［10，31-32］，可通過(guò)灌江納苗及江湖連通，將長(zhǎng)江鳙種質(zhì)資源補(bǔ)充到湖泊中，提高湖泊鳙群體的遺傳多樣性。

3.2 鳙群體的遺傳結(jié)構(gòu)

通常采用遺傳分化指數(shù)判定群體間的遺傳差異大小。Wright［33］認(rèn)為，F(xiàn)st值在0～0.05表示群體分化較弱，在0.05～0.15表示群體中等分化，在0.15～0.25表示群體分化較大，F(xiàn)st值大于0.25表示群體分化極大。本研究顯示，鳙群體間的Fst為-0.022 8～0.044 5，均小于0.05，表明群體間沒(méi)有產(chǎn)生明顯的遺傳分化。分子方差分析表明，大部分的遺傳變異來(lái)自于群體內(nèi)個(gè)體間，僅有0.61%的遺傳變異源于群體間。群體間的遺傳距離為0.001 2～0.002 7，明顯低于地理群體間分化閾值［34］。群體擁有多個(gè)共享單倍型，且共享單倍型的個(gè)體數(shù)量占據(jù)絕對(duì)優(yōu)勢(shì)（92.4%），均說(shuō)明6個(gè)群體間親緣關(guān)系較近，遺傳結(jié)構(gòu)趨于同質(zhì)化。這種遺傳結(jié)構(gòu)模式可能受以下因素的影響：（1）增殖放流。湖泊鳙增殖放流的苗種主要來(lái)源于省內(nèi)的原良種場(chǎng)，通過(guò)微衛(wèi)星分子標(biāo)記檢測(cè)發(fā)現(xiàn)，鳙親本沒(méi)有明顯的遺傳分化［35］，所以放流群體的遺傳變異較小。（2）江淮下游河網(wǎng)密布，鳙群體間直接或間接相連通，有利于群體間基因交流［36］。（3）鳙的運(yùn)動(dòng)能力和適應(yīng)能力強(qiáng)，活動(dòng)范圍廣，群體間地理位置較近，沒(méi)有地理隔離。根據(jù)鳙群體的遺傳結(jié)構(gòu)特征，可以將6個(gè)湖泊鳙群體視為一個(gè)遺傳單元進(jìn)行管理。

3.3 湖泊鳙的資源管理和保護(hù)措施

本研究通過(guò)線粒體COⅠ基因序列，獲得了6個(gè)湖泊鳙群體的遺傳多樣性水平和遺傳結(jié)構(gòu)特征。結(jié)果顯示，鳙群體的遺傳多樣性較低，其中太湖、滆湖、高郵湖和洪澤湖群體的遺傳多樣性相對(duì)較高，而長(zhǎng)蕩湖和駱馬湖群體的遺傳多樣性相對(duì)較低。6個(gè)鳙群體間遺傳變異較小，沒(méi)有產(chǎn)生遺傳分化，應(yīng)將6個(gè)群體視為一個(gè)進(jìn)化單元進(jìn)行管理。增殖放流是影響湖泊鳙群體遺傳多樣性及遺傳結(jié)構(gòu)的重要因素。

根據(jù)研究結(jié)果，建議采取以下措施保護(hù)和恢復(fù)湖泊鳙的種質(zhì)資源遺傳多樣性：（1）強(qiáng)化漁政管理，禁止非法捕撈，落實(shí)涉漁工程生態(tài)補(bǔ)償措施，開(kāi)展生態(tài)環(huán)境治理、修復(fù)。（2）加強(qiáng)增殖放流監(jiān)督管理，科學(xué)規(guī)范開(kāi)展增殖放流，提高放流苗種的成活率。開(kāi)展親本和苗種遺傳多樣性檢測(cè)，增加親本數(shù)量和種質(zhì)更新速度，提升放流苗種的遺傳多樣性。（3）通過(guò)灌江納苗或江湖連通等方式，將遺傳多樣性較為豐富的長(zhǎng)江鳙種質(zhì)資源補(bǔ)充到湖泊中，增加湖泊鳙群體的資源量及遺傳多樣性。

參考文獻(xiàn)

［1］倪勇，伍漢霖.江蘇魚(yú)類志［M］.北京：中國(guó)農(nóng)業(yè)出版社，2006.

［2］劉建康，謝平.用鰱鳙直接控制微囊藻水華的圍隔試驗(yàn)和湖泊實(shí)踐［J］.生態(tài)科學(xué)，2003，22（3）：193-198.

［3］楊姣姣，過(guò)龍根，尹成杰，等.富營(yíng)養(yǎng)化初期湖泊放養(yǎng)鰱（Hypophthalmichthys molitrix）、鳙（Aristichthys nobilis）控藻生態(tài)效果的初步評(píng)估［J］.湖泊科學(xué)，2019，31（2）：386-396.

［4］趙旭昊，徐東坡，任瀧，等.基于Ecopath模型的太湖鰱鳙生態(tài)容量評(píng)估［J］.中國(guó)水產(chǎn)科學(xué)，2021，28（6）：785-795.

［5］馮照軍，王光標(biāo)，趙彥禹，等.江蘇駱馬湖濕地魚(yú)類資源及其保護(hù)［J］.四川動(dòng)物，2007，26（1）：126-129.

［6］唐晟凱，張彤晴，孔優(yōu)佳，等.滆湖魚(yú)類學(xué)調(diào)查及漁獲物分析［J］.水生態(tài)學(xué)雜志，2009，30（6）：20-24.

［7］毛志剛，谷孝鴻，曾慶飛，等.太湖魚(yú)類群落結(jié)構(gòu)及多樣性［J］.生態(tài)學(xué)雜志，2011，30（12）：2836-2842.

［8］林明利，張?zhí)昧郑~少文，等.洪澤湖魚(yú)類資源現(xiàn)狀、歷史變動(dòng)和漁業(yè)管理策略［J］.水生生物學(xué)報(bào)，2013，37（6）：1118-1127.

［9］張照鵬，董芳，杜浩，等.長(zhǎng)江中下游區(qū)增殖放流現(xiàn)狀與對(duì)策研究［J］.淡水漁業(yè)，2021，51（6）：19-28.

［10］張敏瑩，劉凱，徐東坡，等.長(zhǎng)江下游鳙放流群體和天然捕撈群體遺傳多樣性的微衛(wèi)星分析［J］.江西農(nóng)業(yè)大學(xué)學(xué)報(bào)，2013，35（3）：579-586.

［11］陳會(huì)娟.長(zhǎng)江中游四大家魚(yú)放流親本對(duì)早期資源和遺傳多樣性的影響研究［D］.重慶：西南大學(xué)，2019.

［12］汪鄂洲，李全宏，徐念，等.基于COnbsp; Ⅰ和Cytb序列的丹江口水庫(kù)鰱群體遺傳結(jié)構(gòu)分析［J］.水生態(tài)學(xué)雜志，2022，43（4）：103-110.

［13］郭新紅，劉少軍，劉巧，等.魚(yú)類線粒體DNA研究新進(jìn)展［J］.遺傳學(xué)報(bào)，2004，31（9）：983-1000.

［14］李大命，唐晟凱，劉燕山，等.基于線粒體COⅠ基因序列的江蘇省4個(gè)太湖新銀魚(yú)群體遺傳多樣性分析［J］.海洋漁業(yè)，2020，42（3）：277-286.

［15］歐琳，張余，陳曉芳，等.基于線粒體Cyt b和COⅠ基因的洞庭湖區(qū)養(yǎng)殖克氏原螯蝦遺傳多樣性分析［J］.水產(chǎn)科技情報(bào)，2022，49（3）：137-142.

［16］李福貴，程林慧，何坤，等.基于COⅠ基因序列的泥鰍、大鱗副泥鰍選育群體雜交子代遺傳多樣性分析［J］.水產(chǎn)科技情報(bào)，2023，50（1）：20-28.

［17］WARD R D，ZEMLAK T S，INNES B H，et al.DNA barcoding Australia’s fish species［J］.Philosophical Transactions of the Royal Society of London，2005，360（1462）：1847-1857.

［18］HALL T A.BioEdit：a user-friendly biological sequence alignment editor and analysis program for Windows 95/98 NT［J］.Nucleic Acids Symposium Series，1999，41：95-98.

［19］LARKIN M A，BLACKSHIELDS G，BROWN N P，et al.Clustal W and clustal X version 2.0［J］.Bioinformatics，2007（21）：2947-2948.

［20］LIBRADO P，ROZAS J.DnaSP v5：a software for comprehensive analysis of DNA polymorphism data［J］.Bioinformatics，2009，25（11）：1451-1452.

［21］KUMAR S，STECHER G，TAMURA K.MEGA7：Molecular evolutionary genetics analysis version 7.0 for bigger datasets［J］.Molecular Biology and Evolution，2016，33（7）：1870-1874.

［22］BANDELT H J，F(xiàn)ORSTER P，RHL A.Median-joining networks for inferring intraspecific phylogenies［J］.Molecular Biology and Evolution，1999，16（1）：37-48.

［23］EXCOFFIER L，LISCHER H E L.Arlequin suite ver 3.5：a new series of programs to perform population genetics analyses under Linux and Windows［J］.Molecular Ecology Resources，2010，10（3）：564-567.

［24］范啟，何舜平.長(zhǎng)江流域歺 又魚(yú)種群遺傳多樣性和遺傳結(jié)構(gòu)分析［J］.水生生物學(xué)報(bào)，2014，38（4）：627-635.

［25］王丹，程慶武，楊鎮(zhèn)宇，等.三峽庫(kù)區(qū)鲌屬魚(yú)類線粒體COⅠ基因遺傳多樣性的初步分析［J］.水生生物學(xué)報(bào)，2015，39（5）：1054-1058.

［26］GRANT W，BOWEN B.Shallow population histories in deep evolutionary lineages of marine fishes：insights from sardines and anchovies and lessons for conservation［J］.Journal of Heredity，1998，89（5）：415-426.

［27］李大命，楊子萍，劉燕山，等.基于線粒體COⅠ序列的江淮下游湖泊鰱群體遺傳多樣性和遺傳結(jié)構(gòu)分析［J］.淡水漁業(yè)，2023，53（4）：3-11.

［28］李大命，劉燕山，唐晟凱，等.基于Cyt b基因序列的江淮下游湖泊鰱群體遺傳多樣性分析［J］.海洋漁業(yè)，2023，45（4）：423-433.

［29］吳旭，嚴(yán)美姣，李鐘杰.肖四海湖野生和人工放流鱖群體遺傳結(jié)構(gòu)分析［J］.水生生物學(xué)報(bào)，2010，34（3）：562-568.

［30］張燕萍，陳文靜，汪登強(qiáng)，等.鄱陽(yáng)湖水系草魚(yú)野生及增殖放流群體遺傳多樣性分析［J］.江蘇農(nóng)業(yè)科學(xué)，2013，41（9）：207-211.

［31］張德春，張錫元，楊代淑，等.長(zhǎng)江鳙遺傳多樣性的研究［J］.武漢大學(xué)學(xué)報(bào)（自然科學(xué)版），1999，45（6）：857-860.

［32］沙航，羅相忠，鄒桂偉，等.長(zhǎng)江中游鳙群體的微衛(wèi)星遺傳多樣性分析［J］.淡水漁業(yè)，2020，50（4）：12-17.

［33］WRIGHT S.Evolution and the genetics of populations.4：Variability within and among natural populations［M］.Chicago，Illinois：The University of Chicago Press，1978.

［34］YANG X G，WANG Y Q，ZHOU K Y，et al.Authentication of oviductus ranae and its original animals using molecular marker［J］.Biological amp; Pharmaceutical Bulletin，2002，25（8）：1035-1039.

［35］馮曉婷，張桂寧，薛向平，等.基于SSR標(biāo)記的長(zhǎng)江下游原良種場(chǎng)鳙親本和后備親本種質(zhì)資源現(xiàn)狀分析［J］.中國(guó)水產(chǎn)科學(xué)，2020，27（5）：589-597.

［36］李秀啟，劉峰，冷春梅，等.基于線粒體DNA控制區(qū)的南四湖湖鱭（Coilia ectenes taihuensis）群體遺傳結(jié)構(gòu)和種群擴(kuò)散［J］.湖泊科學(xué)，2015，27（4）：686-692.

Genetic diversity analysis of Aristichthys nobilis populations from

freshwater lakes in Jiangsu province based on mitochondrial DNA" COⅠ gene

LI Daming1， LIU Yanshan1， TANG Shengkai1， ZHANG Zeng1，

WANG Chaoqun1， MU Huan2， WANG Hua3

（1. Freshwater Fisheries Research Institute of Jiangsu Province，

Key Laboratory of Fisheries Resources in Inland Waters of Jiangsu Province，

Nanjing 210017，China;

2. Hongze Lake Fisheries Administration Committee Office of Jiangsu Province，Huaian 223300，China;

3. Fisheries Management Commission of Ge Lake of Jiangsu Province，Changzhou 213100，China）

Abstract： To understand the genetic status of Aristichthys nobilis populations from freshwater lakes in Jiangsu province of China，the genetic diversity and genetic structure of A. nobilis from six freshwater lakes（Tai Lake，Ge Lake，Changdang Lake，Gaoyou Lake，Hongze Lake and Luoma Lake） were analyzed based on" COⅠ gene sequences.The results showed that in the analyzed 630 bp" COⅠ gene fragment，the A+T base composition（55.3%） was higher than the G+C combination（44.7%）.A total of 14 variable sites were detected and 9 haplotypes were defined from 223" COⅠ sequences.The haplotype diversity of the six populations ranged from 0.494 to 0.665，and nucleotide diversity varied 0.000 9 to 0.001 9. The total haplotype diversity（Hd） and nucleotide diversity（Pi） was 0. 576 and 0. 001 2，respectively，indicating high haplotype diversity and low nucleotide diversity pattern，demonstrating low genetic diversity. Analysis of molecular variance（AMOVA） showed that genetic variations mainly occurred within the populations（99.39%），the total genetic differentiation index（Fst） was 0.006 1.The pairwise Fst ranged from -0.022 8 to 0.044 5，and the genetic distance among populations varied from 0.001 2 to 0.002 7，indicating no significant genetic differentiation among populations.The haplotype neighbor-joining phylogenetic tree and network structure diagram both showed that the haplotypes were randomly distributed in six populations，without forming a specific geographic clade.Both neutral tests and mismatch distribution analysis indicate that A. nobilis had experienced significant population expansion.The study showed that the genetic diversity of A. nobilis populations from freshwater lakes in Jiangsu province was relatively low，and their genetic structure tended to be homogenized.So，more measures should be taken to improve the genetic diversity of A. nobilis.

Key words： Aristichthys nobilis; COⅠ gene; genetic diversity; genetic structure; stock enhancement