白 云，王 妍，胡 瑩，張興華，金萬(wàn)昆，高永平，楊建新，董 仕

( 1.天津師范大學(xué) 生命科學(xué)學(xué)院，天津市動(dòng)植物抗性重點(diǎn)實(shí)驗(yàn)室，天津 300387；2.國(guó)家級(jí)天津市換新水產(chǎn)良種場(chǎng)，天津 301500 )

鯉魚(yú)(Cyprinuscarpio)是我國(guó)重要的淡水經(jīng)濟(jì)魚(yú)類。我國(guó)的鯉魚(yú)可分為4個(gè)亞種：西鯉(C.carpiocarpio)，分布于我國(guó)的新疆以及歐洲；鯉(C.carpiohaematopterus)，分布于北到黑龍江，南到長(zhǎng)江、閩江以及臺(tái)灣；華南鯉(C.carpiorubrofuscus)，也叫海鯉、元江鯉、團(tuán)鯉，分布于西江、云南元江以及海南島；杞麓鯉(C.carpiochilia)，分布于云南湖泊中[1]。自20世紀(jì)70年代開(kāi)始，我國(guó)開(kāi)展了多方面的鯉魚(yú)遺傳育種研究工作，培育出很多新品種用于養(yǎng)殖生產(chǎn)或作為雜交用親本，產(chǎn)生了巨大的經(jīng)濟(jì)效益。如興國(guó)紅鯉(C.carpiosinguonensis)、荷包紅鯉(C.carpiowuyuanensis)、豐鯉(C.carpiosinguonensis♀×C.carpiovar.mirrorsplittered♂)、建鯉(C.carpiovar.jian)、德國(guó)鏡鯉(C.carpioL mirror)、散鱗鏡鯉(C.carpio)、烏克蘭鱗鯉(C.carpio)、津新鯉(C.carpiojinxin)等[2]。

線粒體DNA是核外遺傳物質(zhì)，為共價(jià)閉合的環(huán)狀超螺旋結(jié)構(gòu)，具有分子量小、遵循母性遺傳、進(jìn)化速度快等特點(diǎn)，成為研究群體間親緣關(guān)系的重要遺傳標(biāo)記，線粒體DNA的D-loop區(qū)段也稱控制區(qū)，具有高突變、進(jìn)化快等特點(diǎn)，是近年來(lái)線粒體DNA的研究熱點(diǎn)[3-7]。鄭文娟等[8]對(duì)舟山小黃魚(yú)(Pseudosciaenapolyactis)進(jìn)行了線粒體DNA D-loop區(qū)段序列變異的遺傳多樣性分析，袁振興等[9]對(duì)都柳江鯉魚(yú)、鯽魚(yú)(Carassiusauratus)和草魚(yú)(Ctenopharyngodonidellus)線粒體DNA的控制區(qū)進(jìn)行了遺傳多樣性分析，董新培等[10]對(duì)不同地理群體烏鱧(Channaargus)進(jìn)行了線粒體DNA控制區(qū)結(jié)構(gòu)及遺傳多樣性分析。

烏克蘭鱗鯉于1998年由天津市水產(chǎn)技術(shù)推廣站通過(guò)農(nóng)業(yè)部引種中心引進(jìn)，2002年天津市技術(shù)推廣站將親本轉(zhuǎn)至天津市換新水產(chǎn)良種場(chǎng)，2005年通過(guò)全國(guó)水產(chǎn)原種和良種審定委員會(huì)審定[11]。德國(guó)鏡鯉原產(chǎn)于德國(guó)巴伐利亞，由農(nóng)業(yè)部于1984年自德國(guó)引進(jìn)[12-14]?？蝼[鏡鯉(C.carpiovar.specularis)簡(jiǎn)稱框鯉，由德國(guó)鏡鯉經(jīng)選育而來(lái)[15-16]。紅鏡鯉(C.carpioL mirror)是國(guó)家級(jí)天津市換新水產(chǎn)良種場(chǎng)在德國(guó)鏡鯉中發(fā)現(xiàn)的紅色個(gè)體。州河鯉(C.carpioheamatopterus)是天津市薊州區(qū)于橋水庫(kù)中的一種地方野鯉[17-19]。本研究以烏克蘭鱗鯉、德國(guó)鏡鯉、框鱗鏡鯉、紅鏡鯉和州河鯉為材料，進(jìn)行了線粒體DNA D-loop區(qū)段的遺傳多樣性檢測(cè)分析，以期為鯉魚(yú)的遺傳育種提供資料。

1 材料與方法

1.1 材料

5種類型鯉魚(yú)的采樣時(shí)間、地點(diǎn)、樣本含量以及體長(zhǎng)、體質(zhì)量見(jiàn)表1，其中德國(guó)鏡鯉和紅鏡鯉剪取鰭條作為研究材料。

表1 5種類型鯉魚(yú)的采樣時(shí)間、地點(diǎn)、樣本含量及體長(zhǎng)、體質(zhì)量

注：*，德國(guó)鏡鯉及紅鏡鯉采集了鰭條.

1.2 總DNA的提取

剪取約1.0 cm×0.5 cm的鰭條，采用苯酚—氯仿抽提法提取總DNA，4 ℃保存?zhèn)溆肹20]。

1.3 線粒體DNA D-loop及其鄰近區(qū)段的擴(kuò)增

使用位于線粒體DNA D-loop區(qū)段兩側(cè)的引物進(jìn)行PCR擴(kuò)增，引物序列為，L15923： 5′-TTA AAG CAT CGG TCT TGT AA-3′[21]，H1067： 5′-ATA GTG GGG TAT CTA ATC CCA GTT-3′[22]，結(jié)合位點(diǎn)分別位于tRNAPro和12S rRNA兩個(gè)編碼區(qū)上，PCR反應(yīng)體系參照文獻(xiàn)[23]的方法。獲得的PCR產(chǎn)物經(jīng)電泳檢測(cè)后，委托華大基因公司進(jìn)行測(cè)序，得到測(cè)序結(jié)果后，再將序列與峰圖進(jìn)行比對(duì)。

1.4 數(shù)據(jù)處理

使用Bioxm軟件將序列轉(zhuǎn)為Fast格式，經(jīng)Seqverter進(jìn)行整合，用ClustalX 1.83將序列進(jìn)行分析，與峰圖比對(duì)后獲得排列清晰的序列，再用Mega 6.0軟件得出堿基差異、堿基頻率，基于Kimura′s 2-Parameter計(jì)算組內(nèi)和組間遺傳距離，并構(gòu)建鄰接系統(tǒng)樹(shù)。

2 結(jié) 果

2.1 線粒體DNA D-loop及其鄰近區(qū)段的PCR擴(kuò)增

5種類型鯉魚(yú)線粒體DNA D-loop及鄰近區(qū)段PCR擴(kuò)增均獲得了約1.6 kb的DNA片段。部分

1.3.6 PPCI術(shù)后3月，患者復(fù)查超聲心動(dòng)圖，了解心功能情況，記錄左室射血分?jǐn)?shù)（left ventricle ejection fration，LVEF）；觀察LVEF＞50%的發(fā)生率。

個(gè)體的PCR產(chǎn)物電泳結(jié)果見(jiàn)圖1。

2.2 PCR產(chǎn)物序列分析

PCR擴(kuò)增產(chǎn)物測(cè)序后，與峰圖比對(duì)，獲得堿基排列順序清晰的DNA片段長(zhǎng)度為1342～1344 bp，經(jīng)與GenBank中的KU050703比對(duì)，序列中包括tRNAPro區(qū)段13 bp，D-loop區(qū)段928、930 bp或931 bp，tRNAPhe區(qū)段69 bp以及12S rRNA區(qū)段329、331 bp或332 bp。97尾個(gè)體中共有8種單倍型，其中6種單倍型的序列提交GenBank，獲得序列號(hào)為MG786481～MG786483，MG786485～MG786487，另外兩種與郝君等[5]得到的序列一致，序列號(hào)為KC292935、KC292936。

圖1 5種類型鯉魚(yú)部分個(gè)體的線粒體DNA D-loop及其鄰近區(qū)段的PCR擴(kuò)增結(jié)果1、2：烏克蘭鱗鯉；3、4：德國(guó)鏡鯉；5、6：框鱗鏡鯉；7、8：紅鏡鯉；9、10：州河鯉.

2.3 5種類型鯉魚(yú)的堿基差異和堿基組成

2.3.1 堿基差異

5種類型鯉魚(yú)8種單倍型間共有28個(gè)位點(diǎn)存在堿基差異(表2)。每種單倍型中包含不同類型鯉魚(yú)的個(gè)體數(shù)(表3)。每種類型鯉魚(yú)在單倍型中有一定傾向性，可以將8種單倍型分為2組。單倍型Ⅰ～Ⅳ為第1組，包含烏克蘭鱗鯉、德國(guó)鏡鯉、紅鏡鯉的全部個(gè)體以及框鱗鏡鯉中近70%的個(gè)體，其中單倍型Ⅰ和Ⅱ包含了烏克蘭鱗鯉的94.1%、德國(guó)鏡鯉的93.75%、框鱗鏡鯉的69.57%、紅鏡鯉的68.42%，占這4種類型鯉魚(yú)的80%。單倍型Ⅵ～Ⅷ為第2組，包含州河鯉的90.91%。

表2 鯉魚(yú)8種單倍型的堿基差異

表3 5種類型鯉魚(yú)在不同單倍型中的分布

2.3.2 堿基組成

使用Mega 6.0軟件分別得出5種類型鯉魚(yú)的線粒體DNA D-loop區(qū)段、線粒體DNA D-loop及鄰近區(qū)段的堿基頻率(表4)。A+T的含量明顯高于G+C的含量，這與魚(yú)類線粒體DNA的蛋白質(zhì)編碼基因的特點(diǎn)[24-25]一致。

2.4 5種類型鯉魚(yú)的遺傳距離和系統(tǒng)發(fā)育樹(shù)

2.4.1 單倍型間的遺傳距離和系統(tǒng)發(fā)育樹(shù)

2.4.1.1 D-loop及鄰近區(qū)段的遺傳距離和系統(tǒng)發(fā)育樹(shù)

表4 5種類型鯉魚(yú)D-loop及鄰近區(qū)段的堿基頻率 %

表5 鯉魚(yú)8種單倍型的D-loop及鄰近區(qū)段間的遺傳距離

圖2 8種單倍型D-loop及鄰近區(qū)段的鄰接系統(tǒng)樹(shù)

2.4.1.2 D-loop區(qū)段的遺傳距離和系統(tǒng)發(fā)育樹(shù)

將8種單倍型與以上3個(gè)GenBank序列以及鯉魚(yú)AY345336、AB307063比對(duì)，經(jīng)ClustalX軟件比對(duì)后獲得926～929 bp的D-loop區(qū)段，單倍型間遺傳距離見(jiàn)表6，構(gòu)建的系統(tǒng)樹(shù)見(jiàn)圖3。由表6可見(jiàn)，鯉魚(yú)AP017364與單倍型Ⅰ遺傳距離為0，松浦鏡鯉KU050703、鯉魚(yú)AY345336與單倍型Ⅱ的遺傳距離為0，鯉魚(yú)KJ511883、鯉魚(yú)AB307063與單倍型Ⅴ的遺傳距離為0，5種類型鯉魚(yú)8種單倍型間的遺傳距離與表5中數(shù)據(jù)相比均有不同程度的增大。由圖3可見(jiàn)，第1組與鯉魚(yú)AP017364、AY345336(具有俄羅斯血統(tǒng)的鏡鯉[28-29])、松浦鏡鯉KU050703聚在一起；第2組與鯉魚(yú)KJ511883、 AB307063聚在一起，兩組差異明顯。

圖3 8種單倍型D-loop區(qū)段的鄰接系統(tǒng)樹(shù)

2.4.2 不同類型鯉魚(yú)間的遺傳距離及系統(tǒng)樹(shù)

2.4.2.1 D-loop及鄰近區(qū)段的遺傳距離和系統(tǒng)發(fā)育樹(shù)

用Mega 6.0軟件得出D-loop及鄰近區(qū)段間的組間及組內(nèi)遺傳距離(表7)并構(gòu)建系統(tǒng)樹(shù)(圖4)，由表7可見(jiàn)，德國(guó)鏡鯉和框鱗鏡鯉的組內(nèi)遺傳距離分別為0.00462和0.00451，遺傳多樣性最高，烏克蘭鱗鯉和紅鏡鯉的遺傳多樣性稍低于德國(guó)鏡鯉和框鱗鏡鯉，州河鯉組內(nèi)遺傳距離為0.00065，遺傳多樣性最低。德國(guó)鏡鯉、框鱗鏡鯉與紅鏡鯉距離較小，為0.00458～0.00544；烏克蘭鱗鯉與框鱗鏡鯉、紅鏡鯉距離稍大，為0.00624～0.00684。州河鯉與烏克蘭鱗鯉、德國(guó)鏡鯉、紅鏡鯉的遺傳距離均約為0.01，與框鱗鏡鯉的遺傳距離稍低，為0.00676；由圖4可見(jiàn)，烏克蘭鱗鯉與鯉魚(yú)AP017364，紅鏡鯉與松浦鏡鯉，州河鯉與鯉魚(yú)KJ511883聚在一起。

表6 鯉魚(yú)8種單倍型的D-loop區(qū)段間的遺傳距離

表7 5種類型鯉魚(yú)的D-loop及鄰近區(qū)段間的遺傳距離

圖4 5種類型鯉魚(yú)D-loop及臨近區(qū)段的鄰接系統(tǒng)樹(shù)

2.4.2.2 D-loop區(qū)段的遺傳距離和系統(tǒng)發(fā)育樹(shù)

D-loop區(qū)段的遺傳距離見(jiàn)表8，構(gòu)建的系統(tǒng)樹(shù)見(jiàn)圖5。由圖5可見(jiàn)，鯉魚(yú)AP017364與烏克蘭鱗鯉、德國(guó)鏡鯉聚在一起，松浦鏡鯉KU050703、鯉魚(yú)AY345336與紅鏡鯉聚在一起。州河鯉與其他4種類型鯉魚(yú)的遺傳距離均約為0.01(表8)，鯉魚(yú)KJ511883、AB307063與州河鯉聚在一起。

2.5 5種類型鯉魚(yú)的遺傳多樣性指數(shù)

用DnaSp_v5軟件計(jì)算出5種類型鯉魚(yú)的單倍型多樣性指數(shù)、核苷酸多樣性指數(shù)及平均核苷酸差異性，結(jié)果見(jiàn)表9。由表9可見(jiàn)，5種類型鯉魚(yú)的單倍型多樣性指數(shù)差異較?。坏聡?guó)鏡鯉和框鱗鏡鯉的核苷酸多樣性指數(shù)和平均核苷酸差異性稍大于烏克蘭鱗鯉和紅鏡鯉，州河鯉明顯小于其他4種類型鯉魚(yú)。

圖5 5種類型鯉魚(yú)D-loop區(qū)段的鄰接系統(tǒng)樹(shù)

表8 5種類型鯉魚(yú)的D-loop區(qū)段間的遺傳距離

表9 5種類型鯉魚(yú)的遺傳多樣性指數(shù)

3 討 論

3.1 5種類型鯉魚(yú)中單倍型可分為歐洲血統(tǒng)和亞洲血統(tǒng)

經(jīng)對(duì)烏克蘭鱗鯉、德國(guó)鏡鯉、框鱗鏡鯉、紅鏡鯉和州河鯉共5種類型鯉魚(yú)的線粒體DNA D-loop及鄰近區(qū)段和僅D-loop區(qū)段的序列分析，共檢測(cè)出8種單倍型。單倍型Ⅰ～Ⅳ包含了烏克蘭鱗鯉、德國(guó)鏡鯉、框鱗鏡鯉和紅鏡鯉中90.7%的個(gè)體，GenBank上的鯉魚(yú)AP017364與單倍型ⅠD-loop區(qū)段，具有俄羅斯血統(tǒng)的鏡鯉AY345336、松浦鏡鯉KU050703與單倍型ⅡD-loop區(qū)段的堿基序列相同，遺傳距離為0。單倍型Ⅰ和單倍型Ⅱ分別占到烏克蘭鱗鯉、德國(guó)鏡鯉、框鱗鏡鯉和紅鏡鯉4種類型鯉魚(yú)的24.0%和56.0%，兩者達(dá)到4種類型鯉魚(yú)的80%。單倍型Ⅲ、Ⅳ與單倍型Ⅰ、Ⅱ遺傳關(guān)系較近。

單倍型Ⅴ～Ⅷ包含了州河鯉的全部個(gè)體、框鱗鏡鯉的7個(gè)個(gè)體。州河鯉為天津市薊州區(qū)于橋水庫(kù)中的1種地方野鯉[17-19]，并且單倍型Ⅴ～Ⅷ的組內(nèi)平均遺傳距離僅為0.00099。

因此，筆者認(rèn)為單倍型Ⅰ～Ⅳ歐洲血統(tǒng)，單倍型Ⅴ～Ⅷ為亞洲血統(tǒng)。至于單倍型Ⅴ中含有框鱗鏡鯉的原因有待進(jìn)一步探討。

3.2 單倍型間的遺傳差異

歐洲血統(tǒng)的Ⅰ～Ⅳ單倍型間堿基差異位點(diǎn)數(shù)較多，為21個(gè)，其中4個(gè)缺失或插入，1個(gè)轉(zhuǎn)換，16個(gè)顛換，單倍型間遺傳距離為0.00224～0.00902，單倍型間遺傳差異較大。亞洲血統(tǒng)的Ⅴ～Ⅷ單倍型間堿基差異位點(diǎn)數(shù)較少，為5個(gè)，其中3個(gè)缺失或插入，2個(gè)顛換，單倍型間遺傳距離為0.00075～0.00149，單倍型間遺傳差異較小。

歐洲血統(tǒng)的單倍型Ⅱ比較特殊，其與其他7種單倍型間的遺傳距離比較相近。單倍型Ⅱ中含有烏克蘭鱗鯉、德國(guó)鏡鯉、框鱗鏡鯉和紅鏡鯉中56.0%的個(gè)體，與具有俄羅斯血統(tǒng)的鏡鯉AY345336[24-25]、松浦鏡鯉KU050703的遺傳距離為0，不含有州河鯉。據(jù)此，筆者認(rèn)為單倍型Ⅱ?qū)儆跉W洲血統(tǒng)。

3.3 5種類型鯉魚(yú)遺傳多樣性

烏克蘭鱗鯉、德國(guó)鏡鯉、框鱗鏡鯉、紅鏡鯉和州河鯉共5種類型鯉魚(yú)中，德國(guó)鏡鯉的組內(nèi)遺傳距離、遺傳多樣性指數(shù)最大，遺傳多樣性最高，與胡雪松等[13]對(duì)德國(guó)鏡鯉的3個(gè)群體進(jìn)行微衛(wèi)星分析，得出德國(guó)鏡鯉選育群體仍保持著較高遺傳多樣性的結(jié)論一致。紅鏡鯉與烏克蘭鱗鯉、德國(guó)鏡鯉和框鱗鏡鯉相比，組內(nèi)遺傳距離、遺傳多樣性指數(shù)最小，遺傳多樣性最低。

州河鯉組內(nèi)遺傳距離、核苷酸多樣性指數(shù)和平均核苷酸差異性明顯小于其他4種類型鯉魚(yú)，遺傳多樣性較低。州河鯉屬于天津地方野鯉，為天津市薊州區(qū)于橋水庫(kù)所特有，其形態(tài)及品質(zhì)和其他鯉魚(yú)有很多不同之處[13,30]?？赡苡捎陂L(zhǎng)期生殖隔離，使得州河鯉在小群體的基礎(chǔ)上繁殖，其遺傳多樣性相對(duì)較低。這與楊正玲等[30]應(yīng)用線粒體DNA D-loop的RFLP方法分析得出的結(jié)論相同。

作為育種材料，不同類型鯉魚(yú)在我國(guó)魚(yú)類遺傳育種研究以及生產(chǎn)應(yīng)用中占有舉足輕重的地位，至今已培育出很多鯉魚(yú)或鯉鯽雜交新品種[2]。所以保護(hù)不同類型鯉魚(yú)的種質(zhì)資源、促進(jìn)鯉魚(yú)養(yǎng)殖業(yè)健康發(fā)展顯得非常重要。