徐鋼春，魏廣蓮，李建林，徐跑，張呈祥，顧若波

(1.中國水產科學研究院淡水漁業(yè)研究中心，農業(yè)部淡水漁業(yè)和種質資源利用重點實驗室，江蘇無錫214081;2.南京農業(yè)大學無錫漁業(yè)學院，江蘇無錫214081;3.江陰市水產指導站，江蘇江陰214431)

刀鱭Coilia nasus隸屬于鯡形目、鳀科、鱭屬，俗稱刀魚，為中國長江重要的經濟洄游性魚類［1-2］，其脂肪豐滿、肉質細嫩鮮美，享有“長江三鮮”之一的美譽;但近年來刀鱭天然資源量急劇減少，種質小型化明顯，遠遠不能滿足市場的需求［3-4］。隨著刀鱭灌江納苗［5］、幼魚采捕［6］的成功實施后，人工繁殖試驗也取得了成功［7］，刀鱭的池塘養(yǎng)殖正逐漸成為新興的產業(yè)。程起群等［8］以線粒體細胞色素b基因片段序列為分子標記，分析了長江野生刀鱭和湖鱭Coilia nasus taihuensis的遺傳關系，認為湖鱭是刀鱭的一個地理種群且遺傳多樣性水平增加，同時還認為，隔離和適應是形態(tài)差異形成的重要因素［9］。隨著刀鱭人工養(yǎng)殖規(guī)模的擴大，很有必要對刀鱭淡水養(yǎng)殖群體和湖鱭的遺傳多樣性和遺傳分化狀況進行深入研究，從而促進刀鱭人工養(yǎng)殖的可持續(xù)發(fā)展。

在mtDNA上，D-loop區(qū)的進化速度最快，用它來進行種內、種群間的系統(tǒng)進化分析，是一種在分子水平上研究魚類種群遺傳多樣性比較經濟有效的手段［10-14］。本研究中，作者以灌江納苗后完全在淡水池塘環(huán)境中繁育、養(yǎng)殖的刀鱭群體與湖泊定居的湖鱭為試驗材料，運用PCR技術對D-loop區(qū)段進行擴增，旨在從分子水平上分析其遺傳多樣性及相互間的差異水平，為保護刀鱭種質資源和有效地開展刀鱭遺傳育種研究提供參考資料。

1 材料與方法

1.1 材料

試驗用養(yǎng)殖刀鱭11尾，為2004年灌江納苗養(yǎng)殖并繁育的子三代，于2010年8月采自中國水產科學研究院長江三鮮養(yǎng)殖與示范基地，體質量為17.7～85.5 g，平均為36.2 g;湖鱭8尾，于2010年9月采自太湖梅梁灣段，體質量為25.7～40.4 g，平均為33.5 g。從每尾魚剪取少許背部肌肉分別放入Eppendorf管中，于-20℃下保存。

1.2 方法

1.2.1 DNA提取 將肌肉剪碎，分別放入1.5 mL Eppendorf管中，加入470 μL SET緩沖液，再依次加入12.5 μL 20% 的 SDS 和 10 μL(10 mg/mL)蛋白酶K，混合后在55℃水浴中消化過夜，再分別用飽和酚、酚氯仿、氯仿/異戊醇抽提去除蛋白，最后用預冷的無水乙醇沉淀DNA，烘干后用TE溶解。所得基因組DNA樣品用紫外分光光度計測定DNA樣品的濃度和純度，同時用瓊脂糖凝膠電泳檢測基因組DNA的完整性并估測分子量。將DNA原樣稀釋至50 ng/μL，于4℃下保存?zhèn)溆谩?/p>

1.2.2 PCR擴增及測序 使用兩個通用引物DF1和DR2擴增D-loop區(qū)段基因，引物序列DF1:5'-CTAACTCCCAAAGCTAGAATTCT-3'，DR2:5'-ATCTTAGCATCTTCAGTG-3'［15］。擴增反應混合物總體積25 μL，包括 DNA 模板 1.5 μL，10×Reaction buffer 2.5 μL，200 μmol/L dNTP 0.5 μL，引物(濃度 0.2 μmol/L)各 1 μL，Taq 酶 (1 U)0.2 μL，用無菌ddH2O補至25 μL。PCR擴增程序為:94℃下預變性3 min;94℃下變性30 s，50℃下退火1 min，72℃下延伸1 min，共進行35個循環(huán);最后在72℃下再延伸8 min。PCR產物用14 g/L瓊脂糖凝膠電泳檢測，并拍照記錄。目標擴增產物用試劑盒回收純化后，由上海博尚生物有限公司進行雙向測序，本研究中，除用引物DF1和DR2由兩端向中部測序外，還根據所得的序列再設計一對引物 (F:5'-GCGGATTCTCCGTAGACAAC-3'，R:5'-GTAAAATTGTCTGGGTCTCC-3')由序列中部向兩端進行測序，以增加測序的準確性。

1.3 數(shù)據處理

序列用Clustal W 1.83軟件［16］分析比對。利用DnaSP 4.0軟件［17］統(tǒng)計單倍型及變異位點(S)、計算單倍型多樣性 (H)、核苷酸多樣性 (π)及平均核苷酸差異數(shù) (K)。用Mega 3.0軟件［18］中的Kimura雙參數(shù)法計算各單倍型間的遺傳距離。采用NJ法 (Neighbor-joining)對基因序列數(shù)據進行分析，Bootstrap置信值估算重復次數(shù)為1 000次。

2 結果

2.1 線粒體D-loop區(qū)段基因擴增結果

養(yǎng)殖刀鱭和湖鱭的D-loop區(qū)段基因均能被清晰穩(wěn)定地擴增，PCR產物電泳圖譜見圖1。從圖1可見，D-loop區(qū)段基因長度為1 200 bp左右。PCR產物經回收后測序，測出的序列參照GenBank中已有的刀鱭D-loop區(qū)段序列 (登錄號為NC_009579)，經Clastal X比對分析，結果表明同源性都高達99%以上，確定所得序列為刀鱭mtDNA D-loop區(qū)序列。

養(yǎng)殖刀鱭和湖鱭兩個群體D-loop區(qū)段基因片段中，A、T、G、C堿基含量分別為 33.3%、33.3%、14.2%、19.2%。其中，A+T堿基含量(66.6%)高于C+G的含量 (33.4%)。

圖1 養(yǎng)殖刀鱭和湖鱭mtDNA D-loop區(qū)段PCR電泳圖Fig.1 The PCR eletrophoretogram of mtDNA D-loop in the populations of farmed Coilia nasus and Coilia nasus taihuensis

2.2 遺傳多樣性分析

通過序列分析，養(yǎng)殖刀鱭和湖鱭的D-loop序列全長都有較大的變異，養(yǎng)殖刀鱭的D-loop序列長度為1 210～1 252 bp，若以刀鱭 (NC_009579)D-loop序列全長1 252 bp為標準，則11尾養(yǎng)殖刀鱭中，只有2尾魚序列長度沒有變化，其他9個個體存在1、3、38 bp 3種片段的缺失;湖鱭的D-loop序列長度為1 252～1 290 bp，8尾湖鱭中有3個個體存在38 bp片段的插入。

在兩個群體19尾個體中，共檢測到變異位點(S)35個，占全部序列的2.79%，其中單一多態(tài)位點14個，簡約信息位點21個。繪制堿基變異位點分布圖如圖2所示。由圖2可以看出，除了堿基在100以內的多態(tài)位點外，其他多態(tài)位點十分一致，這與D-loop區(qū)段基因具有一定的保守性和穩(wěn)定性相符合。本研究中共檢測到12種不同的單倍型，單倍型多態(tài)性 (H)為0.924;核苷酸多樣性(π)為0.0099，平均核苷酸差異數(shù) (K)為4.154;遺傳距離為0～0.020，平均遺傳距離為0.015;養(yǎng)殖刀鱭、湖鱭與七絲鱭 (EF_419803)的平均遺傳距離分別為0.0528、0.0537，而養(yǎng)殖刀鱭與湖鱭的平均遺傳距離為0.0148，明顯小于它們與七絲鱭的平均遺傳距離。養(yǎng)殖刀鱭和湖鱭群體內的各遺傳多樣性參數(shù)如表1所示。由表1可以看出，養(yǎng)殖刀鱭的各遺傳多樣性參數(shù)稍高于湖鱭，說明養(yǎng)殖刀鱭比湖鱭的遺傳多樣性要豐富。

2.3 分子系統(tǒng)樹

以七絲鱭(EF_419803)和刀鱭的D-loop區(qū)段序列(NC_009579)為參考進行系統(tǒng)分析，結果顯示，整個群體聚類成兩大支，七絲鱭為一大支，養(yǎng)殖刀鱭和湖鱭聚在一起組成另一大分支。在養(yǎng)殖刀鱭和湖鱭聚在一起組成的分支中，養(yǎng)殖刀鱭和湖鱭又各自聚在一起，形成獨立的單系類群(圖3)。

圖2 刀鱭兩個群體多態(tài)位點序列分布圖Fig.2 The base sequence of polymorphic fragments between this two populations

表1 養(yǎng)殖刀鱭和湖鱭兩個群體的遺傳多樣性Tab.1 Populational genetic diversity in farmed Coilia nasus and Coilia nasus taihuensis

圖3 基于D-loop區(qū)段序列構建的幾種鱭屬魚類的鄰接樹Fig.3 Molecular phylogenetic tree in several anchovy species based on the D-loop sequences construced by NJ method

3 討論

3.1 刀鱭和湖鱭種群控制區(qū)序列的堿基組成

所測定的19尾刀鱭mtDNA D-loop區(qū)序列中，堿基 A、T、G、C的平均含量分別為 33.3%、33.3%、14.2%、19.2%，從堿基組成可以看出，D-loop區(qū)段表現(xiàn)出很強的堿基組成偏向性，即在A、T、G、C 4種堿基中，G的含量明顯低于其他3種堿基的含量;此外，刀鱭mtDNA D-loop區(qū)中A+T堿基的含量 (66.6%)高于G+C堿基(33.4%)，符合脊椎動物mtDNA D-loop區(qū)堿基組成的特點［19］，并與劉煥章［20］、韓虎峰等［21］關于魚類mtDNA控制區(qū)的研究結果一致，同時與鱭屬其他種類一樣［22］，具有高比例的 A+T可能是mtDNA D-loop序列變異較快的原因之一。

3.2 遺傳變異及遺傳多樣性

本試驗中測定的養(yǎng)殖刀鱭和湖鱭兩個群體的D-loop序列中共檢測出35個變異位點，占全部序列的2.79%，共定義了12個單倍型，顯示出較為豐富的遺傳多樣性，與程起群等［8］、葛家春等［23］和馬春艷等［24］的研究結果一致。

人工養(yǎng)殖環(huán)境中因基礎群體數(shù)量過小、養(yǎng)殖過程中出現(xiàn)的近交衰退等因素可能導致養(yǎng)殖群體遺傳多樣性降低，進而影響?zhàn)B殖魚類的經濟性狀，最終對養(yǎng)殖業(yè)造成較大的危害［25-26］。本試驗中養(yǎng)殖刀鱭各遺傳多樣性參數(shù)稍高于湖鱭，說明養(yǎng)殖刀鱭的遺傳多樣性比湖鱭更為豐富，這點從養(yǎng)殖刀鱭D-loop序列長度存在比湖鱭更多樣的變化也能得到印證。湖鱭是長時期陸封型刀鱭群體［8］，由于奠基者效應 (Founder effect)等因素的影響，湖鱭的遺傳多樣性水平可能會較低。養(yǎng)殖刀鱭雖然沒有和外界群體發(fā)生基因交流，但可能由于其親本是由灌江納苗而培育繁殖的子三代，其有著較強的環(huán)境適應能力、生存能力和進化潛力，仍還保持著較高的遺傳多樣性。然而，在經過數(shù)代或者更長時間的人工繁育后，刀鱭種群的遺傳多樣性水平仍然可能因受到影響而出現(xiàn)衰退的現(xiàn)象。此外，本試驗中還發(fā)現(xiàn)，養(yǎng)殖刀鱭和湖鱭兩個群體的個體間D-loop序列長度的多態(tài)性，序列長度差異主要源于以38 bp為基本單位的不同次數(shù)的片段重復，唐文喬等［15］的研究也出現(xiàn)過相似的結果。有學者認為，這是由于某些魚類增加遺傳變異的一種方式，對維持物種的生存有一定的作用［27］。

Billington等［28］報道魚類種內遺傳距離一般在10%以內，對其他一些動物的D-loop區(qū)段基因系列分析表明，種內個體間的序列差異一般在0～4.06%，差異超過6%的個體間已有明顯的亞種或者種的分化［29］，而本試驗中養(yǎng)殖刀鱭、湖鱭與七絲鱭的差異分別為5.28%和5.37%，超出了一般種內個體間的序列差異范圍，但養(yǎng)殖刀鱭與湖鱭的平均遺傳距離為1.48%，仍屬于種內個體間的序列差異，未達到種或亞種的分化，這與唐文喬等［15］的研究結果相似。雖然系統(tǒng)樹中養(yǎng)殖刀鱭和湖鱭各自聚在一起，形成獨立的單系類群，但養(yǎng)殖刀鱭和湖鱭共同構成1個單系群，說明它們?yōu)閮煞N生態(tài)型種群，并非兩種有效物種，這與程起群等［8］的研究結果一致。

當然，由于野生刀鱭資源量急劇減少使得人工繁育、養(yǎng)殖、放流成為一種新的保護刀鱭資源的有效方法，所以，應該保持一定數(shù)量的人工繁育刀鱭親本，使其適時與洄游刀鱭進行基因交流，以保障養(yǎng)殖刀鱭較高的遺傳多樣性，同時防止養(yǎng)殖群體出現(xiàn)種質退化等現(xiàn)象，科學地保護刀鱭種質資源，達到資源的可持續(xù)利用。

［1］Whitehead P J P.FAO species catalogue.Vol.7.Clupeoid fishes of the world(suborder Clupeioidei).An annotated and illustrated catalogue of the herrings，sardines，pilchards，sprats，shads，anchovies and wolf-herrings.Part 1-Chirocentridae，Clupeidae and Pristigasteridae［M］.Rome:United Nations Development Programme，1985:470-471.

［2］張敏瑩，徐東坡，劉凱，等.長江下游刀鱭生物學及最大持續(xù)產量研究［J］.長江流域資源與環(huán)境，2005，14(6):694-698.

［3］Shi D，Gong D.Strategies of resource conservation for Coilia ectenes in the Yangtze Estuary［J］.Marine Science，2003，25:96-97.

［4］Li Y X，Xie S G，Li Z J，et al.Gonad development of an anadromous fish Coilia ectenes(Engraulidae)in lower reach of Yangtze River，China［J］.Fisheries Science，2007，73(6):1224-1230.

［5］張呈祥，鄭金良.長江刀鱭灌江納苗與養(yǎng)殖［J］.科學養(yǎng)魚，2006(7):26.

［6］沈林宏，戴玉紅，顧樹信，等.長江刀鱭幼魚的采集與運輸技術研究［J］.水產養(yǎng)殖，2011(5):4-6.

［7］Xu G C，Tang X，Zhang C X，et al.First studies of embryonic and larval development of Coilia nasus(Engraulidae)under controlled conditions［J］.Aquaculture Research，2011，42:593-601.

［8］程起群，溫俊娥，王云龍.刀鱭與湖鱭線粒體細胞色素b基因片段多態(tài)性及遺傳關系［J］.湖泊科學，2006，18(4):425-430.

［9］程起群，韓金娣.鱭屬兩種群的形態(tài)變異及綜合判別［J］.湖泊科學，2004，16(4):356-364.

［10］肖武漢，張亞平.魚類線粒體DNA的遺傳與進化［J］.水生生物學報，2000，24(4):384-391.

［11］Sugaya T，Ikeda M，F(xiàn)ujio Y.Comparison for the genetic variabilities of natural and seed populations of Japanese flounder based on PCR-RFLP analysis of mtDNA D-loop［J］.Fish Genetics Breeding Science，1999，28:65-73.

［12］單淇，董仕，吳海防，等.3個群體鳙魚mtDNA D-1oop區(qū)段的限制性片斷長度多態(tài)性分析［J］.中國水產科學，2006，13(2):174-180.

［13］譚書貞，董仕，邊春媛，等.長江流域3個群體草魚mtDNA D-loop區(qū)段的PCR-RFLP分析［J］.南開大學學報:自然科學版，2007，40(3):106-112.

［14］楊正玲，趙世海，龔疏影，等.州河鯉、建鯉和框鱗鏡鯉的mtDNA D-loop的 RFLP 分析［J］.淡水漁業(yè)，2008，38(3):40-45.

［15］唐文喬，胡雪蓮，楊金權.從線粒體控制區(qū)全序列變異看短頜鱭和湖鱭的物種有效性［J］.生物多樣性，2007，15(3):224-231.

［16］Thompson J D，Gibson T J，Plewwniak F，et al.The Clustal X windows interface:flexible strategies for multiple sequence alignment aided by quality analysis tools［J］.Nucleic Acids Research，1997，25(24):4876-4882.

［17］Rozas J，Sánchez-delbarrio J C，Messeguer X，et al.DnaSP，DNA polymerism analyses by the coalescent and other methods［J］.Bioinformatics，2003，19(18):2496-2497.

［18］Kumar S，Tamura K，Nei M.MEGA3:Integrated software for molecular evolutionary genetics analysis and sequence alignment［J］.Brief Bioinform，2004，5(2):150-163.

［19］Broughton R E，Milam J E，Roe B A.The complete sequence of the zebrafish(Danio rerio)mitochondrial genome and evolutionary patterns in vertebrate mitochondrial DNA［J］.Genome Research，2001，11(11):1958-1967.

［20］劉煥章.魚類線粒體DNA控制區(qū)的結構與進化:以鳑鲏魚類為例［J］.自然科學進展，2002，12(3):266-270.

［21］韓虎峰，代應貴.珠江水系光唇裂腹魚可渡河種群mtDNA D-loop序列多態(tài)性分析［J］.中國水產科學，2010，17(1):143-149.

［22］諸廷俊，楊金權，唐文喬.長江口鱭屬魚類線粒體DNA控制區(qū)結構分析［J］.上海水產大學學報，2008，17(2):152-157.

［23］葛家春，曹延，陳嬋娟，等.利用擴增片段長度多態(tài)性技術分析長江刀鱭的遺傳多樣性［J］.南京大學學報，2008，44(3):333-337.

［24］馬春艷，劉敏，馬凌波.長江口刀鱭遺傳多樣性擴增多態(tài)DNA(RAPD)分析［J］.海洋水產研究，2004，25(5):20-24.

［25］謝振宇，杜繼曾，陳學群，等.線粒體控制區(qū)在魚類種內遺傳分化中的意義［J］.遺傳，2006，28(3):362-368.

［26］王軍，全成干，蘇永全，等.宮井洋大黃魚遺傳多樣性的RAPD分析［J］.海洋學報，2001，23(3):87-91.

［27］Zhang S M，Deng H，Wang D Q，et al.Mitochondrial DNA length variation and heteroplasmy in Chinese sturgeon Acipenser sinensis［J］.Acta Genetica Sinica，1999，26:489-496.

［28］Billington N，Hebert P D N.Mitochondrial DNA diversity in fishes and its implications for introductions［J］.Canadian Journal of Fisheries Aquatic Sciences，1991，48:80-94.

［29］Yang X G，Wang Y Q，Zhou K Y.The authentication of oviductus ranae and their original animals by using molecular marker［J］.Biological﹠ Pharmaceutical Bulletin，2002，25:1035-1039.