山東省3個(gè)克氏原螯蝦地理群體線粒體COⅠ基因的序列差異分析

楊玲　李寧　朱樹(shù)人等

摘要 利用PCR技術(shù)擴(kuò)增山東境內(nèi)的3個(gè)克氏原鰲蝦地理群體（濟(jì)南小清河、東平湖和微山湖）CO Ⅰ基因片段，得到長(zhǎng)度約為874bp的片段，分析比較了3群體87尾克氏原螯蝦CO Ⅰ基因序列的變異和遺傳結(jié)構(gòu)。結(jié)果顯示，87條序列共檢測(cè)到8個(gè)變異位點(diǎn)，存在6種單倍型。單倍型多樣性（H）為0.174，核苷酸多樣性（Pi）為0.0036，單倍型多樣性（H）以東平湖群體最高（0.243 ），微山湖群體其次（0.204），小清河群體最低（0.000）；3群體的核苷酸多樣性（H）均較低（<0.001）。分子變異等級(jí)分析（AMOVA）表明，3群體間總的遺傳分化指數(shù)（Fst）為0.023 85（P>0.05），群體間的遺傳變異僅占總遺傳變異的2.38%，而97.62%的遺傳變異源于群體內(nèi)，群體間遺傳分化指數(shù)與遺傳距離均較低。表明山東克氏原鰲蝦野生群體的遺傳多樣性較低，具有較低的遺傳分化水平。

關(guān)鍵詞 克氏原鰲蝦； CO Ⅰ；遺傳多樣性

中圖分類號(hào) S917；Q958.1 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2015）20-041-04

Abstract By PCR amplification， partial CO Ⅰ sequences with an aligned length of 874 bp were obtained from 3 Procambarus clarkii geographical population （Xiaoqing River in Jinan， Dongping Lake and Weishan Lake） in Shandong Province. Variation and genetic structure of partial CO Ⅰ sequences were analyzed from 87 individuals of these populations. A total of 6 haplotypes and 8 variable sites were detected. The highest haplotype diversity was in Dongping Lake population （0.243）， while the lowest haplotype diversity was Xiaoqing River population （0.000）. AMOVA showed that the FST among 3 population were 0.023 85 （P>0.05）. Genetic variation among populations accounted for only 2.38% of the total genetic variation， while 97.62% of genetic variation was from groups. Index of genetic differentiation and genetic distance between groups were lower. Results indicated that the genetic diversity of these 3 wild populations was low， and the 3 wild populations had lower genetic differentiation.

Key words Procambarus clarkii； CO Ⅰ； Genetic diversity

克氏原鰲蝦（Procambarus clarkii） 俗稱小龍蝦，屬螯蝦科原螯蝦屬，其肉味鮮美， 高蛋白低脂肪， 富含鈣、磷、鐵等礦物質(zhì)，深受消費(fèi)者喜愛(ài)。現(xiàn)已廣泛分布于我國(guó)中東部地區(qū)十余個(gè)省市，成為一些湖泊和溝渠的優(yōu)勢(shì)種群[1]，是我國(guó)重要的水產(chǎn)經(jīng)濟(jì)物種。分子標(biāo)記是分析群體遺傳多樣性的一種有效手段，學(xué)者利用AFLP、SSR、ISSR、RAPD等方法對(duì)克氏原螯蝦群體的遺傳多樣性進(jìn)行研究，黃羽等[2]、宋亮等[3]利用AFLP方法進(jìn)行克氏原螯蝦遺傳多樣性分析，揭示了長(zhǎng)江中下游地區(qū)6個(gè)群體的遺傳變異主要存在于群體內(nèi)，隨州、濟(jì)南、寧波和常熟4群體具有豐富的遺傳多樣性，且群體間存在較明顯的遺傳分化；王長(zhǎng)忠等[4]、邢智珺等[5]利用SSR方法進(jìn)行了克氏原螯蝦遺傳多樣性分析，揭示了長(zhǎng)江4個(gè)群體遺傳多樣性處于中等水平，江蘇8個(gè)地理群體遺傳多樣性水平較高且群體間存在中度分化；韓曉磊等[6]對(duì)2個(gè)克氏原螯蝦群體遺傳多樣性進(jìn)行了ISSR分析，顯示2野生群體存在一定程度的遺傳分化，遺傳變異主要存在于群體內(nèi)個(gè)體間；張龍崗等[7]、張瑋等[8]分別對(duì)克氏原螯蝦不同野生群體進(jìn)行了遺傳差異的RAPD分析。利用線粒體細(xì)胞色素C氧化酶Ⅰ亞基（CO Ⅰ）基因?qū)耸显r群體遺傳多樣性的研究相對(duì)較少，筆者利用CO Ⅰ基因?qū)ι綎|省內(nèi)3個(gè)克氏原螯蝦野生地理群體的遺傳多樣性狀況、遺傳結(jié)構(gòu)及遺傳變異水平進(jìn)行本底調(diào)查，以期為其種群遺傳多樣性評(píng)估、種質(zhì)資源保護(hù)和遺傳育種提供分子生物學(xué)依據(jù)。

1 材料與方法

1.1 材料

克氏原螯蝦分別采自于山東濟(jì)南小清河（QH）、濟(jì)寧微山湖（WS）和泰安東平湖（DP）的野生群體，每組群體樣本數(shù)量為28～31個(gè)； PCR引物由上海生工生物技術(shù)有限公司合成；DNA提取試劑盒（SK1206）購(gòu)自上海生工生物技術(shù)有限公司；Taq DNA聚合酶、dNTP 、MgCl2、Marker等為大連寶生物工程有限公司生產(chǎn)；PCR儀為Takara TP650型。

1.2 方法

1.2.1

基因組DNA的提取。取克氏原螯蝦肉約20 mg，用DNA提取試劑盒提取樣品基因組DNA，經(jīng)OD260/OD280比值和瓊脂糖凝膠電泳檢測(cè)濃度和純度。

1.2.2

引物設(shè)計(jì)及PCR擴(kuò)增。根據(jù)GenBank（NC_016926.1）中的相應(yīng)序列，設(shè)計(jì)并合成一對(duì)引物，用于CO Ⅰ基因序列的擴(kuò)增，引物序列為

CO ⅠS： 5′ACGCAACGATGATTTTTTTCT3′；

CO ⅠA： 5′CATCCATCCCTACCGTAAATA 3′。

PCR擴(kuò)增體系50 μl，擴(kuò)增條件：94 ℃預(yù)變性5 min；94 ℃變性45 s， 51 ℃退火50 s，72 ℃ 延伸50 s，40個(gè)循環(huán)；72 ℃延伸10 min。PCR擴(kuò)增產(chǎn)物經(jīng)1%的瓊脂糖凝膠電泳檢測(cè)后，送上海生工生物技術(shù)有限公司，純化后采用上述引物進(jìn)行雙向測(cè)序。

1.2.3

序列分析。所測(cè)得的DNA序列經(jīng)DNASTAR軟件拼接并輔以人工校，與從GenBank中下載的克氏原螯蝦相應(yīng)序列用Clustal W方法進(jìn)行比對(duì)分析和確定長(zhǎng)度，采用MEGA 5.2 軟件統(tǒng)計(jì)序列的堿基組成和轉(zhuǎn)換/顛換比率（Ts/Tv ratios）、變異位點(diǎn)，采用Kimura 2parameter模型計(jì)算群體間的遺傳距離，采用Kimura 2parameter距離矩陣，鄰接法（NeighborJoining，NJ）構(gòu)建單倍型分子系統(tǒng)樹(shù)，自舉檢測(cè)（bootstrap）1000次計(jì)算各分支置信度；用DNASP V5計(jì)算單倍型多樣性指數(shù)（Hd）和核苷酸多樣性指數(shù)（π）；使用Arlequin3.5軟件中的分子變異分析（AMOVA）方法估算遺傳變異在群體內(nèi)和群體間的分布及遺傳分化指數(shù)（Fstatistics，F(xiàn)st）。

2 結(jié)果與分析

2.1 PCR擴(kuò)增結(jié)果

采用引物CO ⅠS和CO ⅠA對(duì)3個(gè)不同地理群體克氏原螯蝦的基因組DNA進(jìn)行PCR擴(kuò)增，結(jié)果都擴(kuò)增出一條約900 bp條帶（圖1）。

2.2 CO Ⅰ基因片段的序列分析

2.2.1

序列分析。對(duì)克氏原螯蝦3群體87個(gè)樣本線粒體CO Ⅰ基因進(jìn)行了測(cè)序，結(jié)果經(jīng)DNASTAR軟件拼接并輔以人工校對(duì)，去除部分端序列后，獲得874 bp的片段，經(jīng)BLAST分析，與GenBank（NC_016926.1）中的克氏原螯蝦CO Ⅰ基因片段同源性高達(dá)99%，確定該產(chǎn)物為克氏原螯蝦CO Ⅰ基因的部分序列。

用CLUSTALW軟件對(duì)測(cè)得序列及GenBank中（NC_016926.1）的同源序列進(jìn)行比對(duì)分析，共有8個(gè)變異位點(diǎn)（占位點(diǎn)總數(shù)的0.9%），其中有5個(gè)轉(zhuǎn)換（transition），3個(gè)顛換（transversion），轉(zhuǎn)換顛換比（Ts/Tv）為1.67∶1.00， 變異位點(diǎn)多發(fā)生在第3密碼子上，沒(méi)有堿基的插入與缺失。在編碼的296個(gè)氨基酸序列中，有1處氨基酸替代，由密碼子第一位點(diǎn)上的核苷酸替代引起。

87個(gè)樣本共檢測(cè)到6種單倍型（表1），其中3群體的共享單倍型1個(gè)（Hap1），單倍型Hap2、Hap3、Hap4為DP群體所特有，單倍型Hap5、Hap6為WS群體所特有，即DP 群體擁有4種單倍型，WS群體有3種單倍型，而QH群體只有1種單倍型。除共享單倍型外的5種單倍型出現(xiàn)頻率較低，只有1～2個(gè)樣本。

2.2.2

堿基組成和遺傳距離。

采用MGEA 5.2軟件統(tǒng)計(jì)3個(gè)群體克氏原螯蝦CO Ⅰ基因部分序列的堿基組成和G+G含量（表2）。G+C含量（32.0%）顯著低于A+T含量（68%），其中堿基C的含量最低，只有12.1%，且在1、2、3位的含量變化很大，分別為12.0%、23.4%和1.0%。DP群體內(nèi)的相對(duì)遺傳距離（Pdistance）為0～0.007，WS群體內(nèi)的遺傳距離為0～0.002，QH群體內(nèi)的遺傳距離為0。

克氏原螯蝦CO Ⅰ基因6種單倍型序列間的遺傳距離見(jiàn)表3，利用Kimura雙參數(shù)模型計(jì)算，置信度估算采用Kimura的雙參數(shù)法，重復(fù)數(shù)1 000，結(jié)果顯示各單倍型間的遺傳距離在0.001～0.007。

2.2.3

不同群體遺傳多樣性及遺傳分化分析。利用DnaSP 5.1軟件進(jìn)行3群體克氏原螯蝦CO Ⅰ基因遺傳多樣性分析（表4），結(jié)果表明，DP和WS群體具有較高的單倍型多樣性（Hd），分別為0.243和0.204，而小清河群體單倍型多樣性非常低，為0。3群體核苷酸多樣性（Pi）水平與平均核苷酸差異數(shù)（K）都較低。

分子變異等級(jí)分析（AMOVA）結(jié)果表明，群體間遺傳分化指數(shù)Fst=0.023 85（ P>0.05），表明在整個(gè)遺傳變異中不同群體間的遺傳分化只占2.38%，其余97.62%的遺傳分化來(lái)自群體內(nèi)的個(gè)體間。群體間的遺傳分化指數(shù)Fst及3群體間的遺傳距離見(jiàn)表5，顯示3群體間DP和WS間的遺傳分化指數(shù)最低（Fst=-0.004 82），DP和QH群體間的遺傳分化指數(shù)最高（Fst=0.052 59），表明根據(jù)線粒體CO Ⅰ基因的序列，3個(gè)群體中沒(méi)有顯著的遺傳分化（P>0.05）。

從3群體間的遺傳距離來(lái)看，DP和WS群體的遺傳距離最大（0.000 42），DP和QH群體的遺傳距離最?。?.000 12），表明遺傳距離與地理距離沒(méi)有相關(guān)性。

2.2.4

系統(tǒng)發(fā)育和聚類分析。根據(jù)6個(gè)單倍型序列和GenBank中檢索的同源序列（NC_016926.1），采用雙參數(shù)模型（Kimura 2parameter）作為距離參數(shù)，鄰接法（neighberjoining，NJ）構(gòu)建系統(tǒng)進(jìn)化樹(shù)，自舉檢測(cè)（bootstrap）1000次計(jì)算各分支置信度（圖2）。由圖2可知，Hap1、Hap2和GenBank匯成一支，另外4個(gè)單倍型匯成一支，可以看出DP群體擁有的4種單倍型在各分支上都有分布，各群體間存在較為廣泛的共享單倍型，未表現(xiàn)出明顯的地理位置與單倍型之間的對(duì)應(yīng)關(guān)系。

利用MGEA 5.2軟件構(gòu)建了87個(gè)克氏原螯蝦個(gè)體的UPGMA系統(tǒng)進(jìn)化樹(shù)（圖3），也揭示了克氏原螯蝦山東3群體不是按照地理位置形成對(duì)應(yīng)族群的，絕大多數(shù)個(gè)體聚為一大類后和WS群體的部分個(gè)體群聚為一支， DP群體的部分個(gè)體自成一支。

3 討論

動(dòng)物線粒體 CO Ⅰ基因進(jìn)化速度較快，適合種群水平差異的檢測(cè)，也可用于種間分析，彭士明等[9]、馬凌波等[10]、姜虎成等[11] 、朱立靜等[12]通過(guò)分析銀鯧、青蟹、日本沼蝦、四角蛤喇等物種群體的線粒體CO Ⅰ基因序列差異研究其遺傳結(jié)構(gòu)和系統(tǒng)演化。該研究通過(guò)分析山東3個(gè)不同地理群體克氏原螯蝦CO Ⅰ基因序列的差異來(lái)探討其遺傳結(jié)構(gòu)、遺傳分化和進(jìn)化關(guān)系。

3.1 遺傳多樣性

種群內(nèi)和種群間的遺傳多樣性水平是反映物種遺傳背景的科學(xué)資料，該研究表明，在克氏原螯蝦3群體87個(gè)樣本擴(kuò)增的 CO Ⅰ基因874 bp片段中，存在8個(gè)突變位點(diǎn)，有6種單倍型。3群體核苷酸多樣性（Pi）為0.000 29±0.000 12，平均核苷酸差異數(shù)（K）0.251±0.082。單倍型多樣性分析表明，DP和 WS群體較高，QH群體最低，與張龍崗等[7]研究的群體遺傳多樣性為WS >DP >QH的RAPD分析結(jié)果相似?？傮w而言，3群體的遺傳多樣性處于較低的水平，各群體間存在較為廣泛的共享單倍型，未表現(xiàn)出明顯的地理位置與單倍型之間的對(duì)應(yīng)關(guān)系。這是因?yàn)榭耸显r為外來(lái)物種，因奠基者效應(yīng)或瓶頸效應(yīng)的影響，由于克氏原螯蝦適應(yīng)環(huán)境的能力非常強(qiáng)，環(huán)境的選擇壓力對(duì)其影響不大，形成基因突變和重組的概率很小，也未見(jiàn)由國(guó)外多次引入的報(bào)道，故其遺傳多樣性增加的可能性微乎其微，故而呈現(xiàn)出遺傳多樣性降低的現(xiàn)象。相對(duì)而言，DP和WS群體的遺傳多樣性水平比QH群體高，這是因?yàn)闁|平湖和微山湖一帶小龍蝦養(yǎng)殖業(yè)發(fā)展迅猛，和長(zhǎng)江流域存在基因交流，人為因素使得其遺傳多樣性水平相對(duì)較高；而QH群體因生存環(huán)境相對(duì)封閉，克氏原螯蝦定居性較強(qiáng)，游泳能力弱，自身遷徙能力差，與其他群體缺乏交流，加之近些年過(guò)度捕撈，造成自然資源銳減和種群內(nèi)近親交配，導(dǎo)致遺傳多樣性水平逐漸降低。王長(zhǎng)忠等[4]利用AFLP技術(shù)進(jìn)行研究，結(jié)果表明，克氏原螯蝦群體遺傳多樣性為寧波群體>隨州群體>濟(jì)南群體>常熟群體，也證明濟(jì)南群體處于較低的遺傳多樣性水平。

3.2 遺傳分化

AMOVA分析的Fst值是用來(lái)測(cè)量群體間遺傳分化的指標(biāo)，F(xiàn)st值越小說(shuō)明群體間發(fā)生遺傳分化越低，一般認(rèn)為Fst<0.05認(rèn)為遺傳分化程度很??；0.050.25表示有顯著差異。該研究中群體間Fst=0.023 85（P>0.05），表明群體間的遺傳分化程度較低，在整個(gè)遺傳變異中群體間占2.38%，其余的遺傳變異來(lái)自于群體內(nèi)，3個(gè)群體中DP和QH間的Fst值最大（0.052 59），兩者間的遺傳分化最大， DP和WS間的Fst最小（-0.004 82），即兩者間的遺傳分化最小，與張龍崗等[7]的研究結(jié)果類似。DP和WS群體間的遺傳距離最大，DP和QH群體間的遺傳距離最?。?.000 12），可能與DP和WS群體CO Ⅰ基因具有獨(dú)有單倍型有關(guān)，導(dǎo)致遺傳分化大的遺傳距離反而小。很多情況下，遺傳分化大的2個(gè)群體間遺傳距離不一定最大[13]，群體間Fst值雖然大小與其趨勢(shì)不一致，但該研究Fst值都不顯著，說(shuō)明3個(gè)群體間不存在遺傳分化，這個(gè)Fst和遺傳距離不是對(duì)應(yīng)的。

UPGMA系統(tǒng)關(guān)系樹(shù)也反映出它們的系統(tǒng)進(jìn)化情況，絕大多數(shù)個(gè)體聚為一大類后和WS群體的部分個(gè)體群聚為一支， DP群體的部分個(gè)體自成一支，由此推測(cè)，可能東平湖群體更接近于外來(lái)原始種群，微山湖群體和小清河群體由東平湖群體演化而來(lái)。

該研究山東3個(gè)克氏原螯蝦野生群體間的遺傳分化較低，群體間還沒(méi)有形成互相隔離的遺傳結(jié)構(gòu)，鑒于此，在克氏原螯蝦選育過(guò)程中，應(yīng)以東平湖群體作為基礎(chǔ)群體，并適當(dāng)引進(jìn)外地良種，以獲得更多遺傳多樣性水平，促進(jìn)克氏原螯蝦養(yǎng)殖業(yè)的健康發(fā)展。

參考文獻(xiàn)

[1]王衛(wèi)民.軟殼克氏原螯蝦在我國(guó)開(kāi)發(fā)利用前景[J].水生生物學(xué)報(bào)，1999，23（4）：375-381.

[2] 黃羽，戴銀根，畢成武，等.長(zhǎng)江中下游地區(qū)6個(gè)克氏原螯蝦群體遺傳多樣性分析[J].南昌大學(xué)學(xué)報(bào)，2011，33（3）：243-247.

[3] 宋亮，韓曉磊，夏蟬，等.不同地區(qū)克氏原螯蝦遺傳多樣性的AFLP分析[J].江蘇農(nóng)業(yè)學(xué)報(bào)，2012，28（2）：359-364.

[4] 王長(zhǎng)忠，李忠，梁宏偉，等.長(zhǎng)江下游地區(qū)4個(gè)克氏原螯蝦群體的遺傳多樣性分析[J].生物多樣性，2009，17（5）：518-523.

[5] 邢智珺，江虎成，陸偉，等.江蘇8個(gè)克氏原螯蝦群體遺傳多樣性微衛(wèi)星分析[J].上海海洋大學(xué)學(xué)報(bào)，2014，23（5）：656-661.

[6] 韓曉磊，徐松維，李小蕊，等.克氏原螯蝦群體遺傳多樣性的ISSR分析[J].江蘇農(nóng)業(yè)科學(xué)，2011，39（4）：26-29.

[7] 張龍崗，楊玲，劉羽清，等.山東克氏原螯蝦3個(gè)地理群體遺傳差異的RAPD 分析[J].長(zhǎng)江大學(xué)學(xué)報(bào)，2014，11（17）：37-40.

[8] 張瑋，阮龍，陳義紅，等.克氏原螯蝦4個(gè)地理群體遺傳差異的RAPD分析[J].安徽農(nóng)業(yè)科學(xué)，2010，38（22）：11861-11862，11876.

[9] 彭士明，施兆鴻，侯俊利，等.銀鯧3個(gè)野生群體線粒體CO Ⅰ基因的序列差異分析[J].上海海洋大學(xué)學(xué)報(bào)，2009，18（4）：398-402.

[10] 馬凌波，張鳳英，喬振國(guó)，等.中國(guó)東南沿海青蟹線粒體CO Ⅰ基因部分序列分析[J].水產(chǎn)學(xué)報(bào)，2006，30（4）：463-467.

[11] 姜虎成，馮建彬，丁懷宇，等.淮河水系日本沼蝦群體遺傳結(jié)構(gòu)和系統(tǒng)演化的線粒體CO Ⅰ序列分析[J].動(dòng)物學(xué)雜志，2012，47（2）：73-84.

[12] 朱立靜，陳淑吟，許曉風(fēng)，等.四角蛤喇江蘇群體線粒體CO Ⅰ基因片段序列研究[J].江蘇農(nóng)業(yè)科學(xué)，2010，38（4）：33-35，97.

[13] 張東亞，汪登強(qiáng)，劉紹平，等.怒江頻危魚類缺須盆唇魚基于線粒體Cyt b 序列的群體遺傳結(jié)構(gòu)分析[J].中國(guó)水產(chǎn)科學(xué)，2009，16（4）：477-485.