郝君,楊薔,鮑迪,梁愛軍,張興華,董仕

(天津師范大學(xué)生命科學(xué)學(xué)院天津市動(dòng)植物抗性重點(diǎn)實(shí)驗(yàn)室,天津300387)

6種魚mtDNA D-loop及其鄰近區(qū)段的序列比較分析

郝君,楊薔,鮑迪,梁愛軍,張興華,董仕

(天津師范大學(xué)生命科學(xué)學(xué)院天津市動(dòng)植物抗性重點(diǎn)實(shí)驗(yàn)室,天津300387)

采用PCR技術(shù)對(duì)烏克蘭鱗鯉Cyprinus carpio、鯽Carassis auratus、鰱Hypophthalmichthys molitrix、鳙Aristichthys nobilis、草魚Ctenopharyngodon idellus和烏蘇里擬鲿Pseudobagrus ussuriensis共28個(gè)個(gè)體(分屬于天津市動(dòng)植物抗性重點(diǎn)實(shí)驗(yàn)室應(yīng)用PCR-RFLP方法已檢測(cè)出的不同單倍型)的mtDNA D-loop及其鄰近區(qū)段進(jìn)行擴(kuò)增并測(cè)序,獲得了大小為1 500～1 800 bp的擴(kuò)增產(chǎn)物。用Clustal 1.83軟件與GenBank中的鯉、金魚、鰱、草魚、長(zhǎng)吻鮠5種魚的mtDNA全序列進(jìn)行排序比較,確定擴(kuò)增產(chǎn)物包含完整的tRNAPro、D-loop、tRNAPhe序列以及12S rRNA長(zhǎng)度約為400 bp的部分序列。將6種魚共28個(gè)樣本的序列提交GenBank后獲得序列號(hào)為KC292921～KC292948。運(yùn)用Mega 4.0軟件計(jì)算出6種魚的堿基組成和堿基差異,基于Kimura雙參數(shù)模型計(jì)算種間的遺傳距離,并與GenBank中5種魚的mtDNA序列一起構(gòu)建NJ系統(tǒng)樹。序列結(jié)構(gòu)分析表明,在6種魚序列的4個(gè)區(qū)段中,D-loop區(qū)段在種內(nèi)、種間的差異性均高于另外3個(gè)區(qū)段。6種魚NJ系統(tǒng)樹的結(jié)果與傳統(tǒng)分類方法一致,可為魚類分類以及種類鑒定提供科學(xué)依據(jù)。另外,6種魚各種單倍型均在所測(cè)得序列中找到相應(yīng)的酶切位點(diǎn),印證了RFLP實(shí)驗(yàn)結(jié)果,在酶切位點(diǎn)以外也發(fā)現(xiàn)了同種魚的不同單倍型間存在堿基差異。

烏克蘭鱗鯉;鯽;鰱;鳙;草魚;烏蘇里擬鲿;D-loop;NJ系統(tǒng)樹

魚類線粒體DNA(Mitochondrial DNA,mtDNA)是一種環(huán)狀共價(jià)閉合的雙鏈超螺旋分子,為核外遺傳物質(zhì)。具有嚴(yán)格的母系遺傳、結(jié)構(gòu)簡(jiǎn)單、分子量小、無(wú)組織特異性和進(jìn)化速率快等特點(diǎn),被廣泛應(yīng)用于物種起源、遺傳分化、種內(nèi)種間系統(tǒng)發(fā)生關(guān)系及原種鑒定等方面[1-9]。

D-loop區(qū)是mtDNA上最重要的非編碼區(qū),受進(jìn)化壓力小,變異程度高,進(jìn)化速率是mtDNA其他區(qū)段的2～5倍。D-loop區(qū)位于 tRNAPro和 tRNAPhe之間,包含與DNA復(fù)制和轉(zhuǎn)錄相關(guān)的序列,是分析相近物種種間和種內(nèi)遺傳多樣性的有效分子遺傳標(biāo)記,也是近年來(lái) mtDNA 的研究熱點(diǎn)[6-7,10-12]。在測(cè)序技術(shù)未被廣泛應(yīng)用于mtDNA研究之前,對(duì)魚類的種群結(jié)構(gòu)、物種多樣性方面的研究大多基于RFLP分析[13-14]。隨著PCR技術(shù)以及自動(dòng)化測(cè)序技術(shù)的發(fā)展,已有許多研究人員將D-loop區(qū)作為分子標(biāo)記應(yīng)用于魚類遺傳分化的研究[1-7]。

鯉Cyprinus carpio與鯽Carassis auratus是中國(guó)最重要的淡水養(yǎng)殖魚類,在世界范圍內(nèi)廣泛養(yǎng)殖,歷史悠久,并形成許多亞種和品系,鯉品種之一的烏克蘭鱗鯉的養(yǎng)殖已產(chǎn)生了巨大的經(jīng)濟(jì)效益[15];鰱Hypophthalmichthys molitrix、鳙Aristichthys nobilis和草魚Ctenopharyngodon idellus是中國(guó)特產(chǎn)魚類,與青魚 Mylopharyngodon piceus一起被譽(yù)為中國(guó)“四大家魚”,是世界水產(chǎn)養(yǎng)殖業(yè)發(fā)展的重要養(yǎng)殖種類[16-19];烏蘇里擬鲿 Pseudobagrus ussuriensis主要分布于中國(guó)黑龍江水系,其營(yíng)養(yǎng)豐富,肉質(zhì)鮮美,具有很高的經(jīng)濟(jì)價(jià)值[20-21]。本研究中,作者選用天津市動(dòng)植物抗性重點(diǎn)實(shí)驗(yàn)室進(jìn)行過(guò)mtDNA D-loop區(qū)PCR-RFLP分析的6種魚的部分個(gè)體為試驗(yàn)材料,對(duì)其D-loop及其鄰近區(qū)段進(jìn)行PCR擴(kuò)增和序列測(cè)定,對(duì)測(cè)序結(jié)果進(jìn)行種內(nèi)、種間遺傳結(jié)構(gòu)差異分析,構(gòu)建系統(tǒng)樹,探討D-loop及其鄰近區(qū)段作為分子標(biāo)記對(duì)系統(tǒng)分化研究的意義。同時(shí)對(duì)6種魚PCR擴(kuò)增區(qū)段的序列結(jié)構(gòu)與PCR-RFLP檢測(cè)結(jié)果進(jìn)行比較,旨在為魚類分類以及種類鑒定提供科學(xué)依據(jù)。

1 材料與方法

1.1 材料

試驗(yàn)用烏克蘭鱗鯉2尾(每種單倍型各1尾)[15],鯽3尾(采自子牙河天津段)[22],鰱12尾(每種單倍型各1尾)[23],鳙7尾(每種單倍型各1尾)[24],草魚2尾(每種單倍型各1尾)[25],烏蘇里擬鲿2尾(采自天津市天祥水產(chǎn)養(yǎng)殖場(chǎng))。

1.2 方法

1.2.1 基因組DNA的提取 從每尾鮮活試驗(yàn)魚體上剪下大小約1.0 cm×0.5 cm的尾鰭鰭條,用苯酚-氯仿抽提法提取,于4℃下保存?zhèn)溆肹26]。

1.2.2 mtDNA的D-loop及其鄰近區(qū)段的擴(kuò)增及凝膠電泳 以提取的 DNA為模板,使用位于mtDNA的D-1oop區(qū)段兩側(cè)的引物對(duì)每尾試驗(yàn)魚的D-loop及其鄰近區(qū)段進(jìn)行擴(kuò)增。所用引物為L(zhǎng)15923:5′TTA AAG CAT CGG TCT TGT AA 3′[27]和H1067:5′ATA GTG GGG TAT CTA ATC CCA GTT 3′[28],兩個(gè)引物的結(jié)合位點(diǎn)分別位于tRNAThr和12S rRNA編碼區(qū)上。PCR方法與擴(kuò)增片段的確認(rèn)參照梁愛軍等[15]的方法。將PCR產(chǎn)物委托深圳華大基因科技有限公司進(jìn)行測(cè)序。測(cè)序得到的序列依據(jù)峰圖進(jìn)行確認(rèn)。

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

以從GenBank中獲得的鯉序列(X61010.1)、金魚序列(AB111951.1)、鰱序列(EU315941.1)、草 魚 序 列 (EU391390.1) 和 長(zhǎng) 吻 鮠 序 列(GU596454.1)為參考序列,用Clustal 1.83軟件對(duì)測(cè)序結(jié)果進(jìn)行同源序列比對(duì)分析[29],確定PCR擴(kuò)增產(chǎn)物所包含的各個(gè)區(qū)段。將比對(duì)好的6種魚28個(gè)個(gè)體的序列提交到GenBank庫(kù)中獲得序列號(hào)。用Mega 4.0軟件統(tǒng)計(jì)6種魚D-loop及其鄰近區(qū)段的堿基組成和堿基差異[30],基于Kimura-雙參數(shù)計(jì)算6種魚群體間的遺傳距離[31],并依此與Gen-Bank中5種魚的序列一起構(gòu)建NJ系統(tǒng)樹。

2 結(jié)果

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

6種魚mtDNA D-loop及其鄰近區(qū)段的PCR擴(kuò)增產(chǎn)物經(jīng)瓊脂糖凝膠電泳后,得到大小不同的擴(kuò)增片段(圖1)。PCR產(chǎn)物電泳結(jié)果顯示條帶清晰單一,未發(fā)現(xiàn)非特異性條帶。其中1號(hào)草魚個(gè)體的片段長(zhǎng)度約為1 800 bp,12、13號(hào)烏蘇里擬鲿兩個(gè)個(gè)體的片段長(zhǎng)度約為1 500 bp,其他個(gè)體的片段長(zhǎng)度約為1 600 bp。

圖1 6種魚mtDNA D-loop及其鄰近區(qū)段的PCR電泳圖Fig.1 The PCR electrophoretogram of mtDNA D-loop and adjacent regions in six fish species

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

將測(cè)序結(jié)果與GenBank中的序列進(jìn)行比對(duì)后發(fā)現(xiàn),6種魚使用引物L(fēng)15923和H1067擴(kuò)增出的片段包括完整的tRNAPro、D-loop、tRNAPhe序列以及12S rRNA長(zhǎng)度約為400 bp的部分序列。將6種魚共28個(gè)個(gè)體的序列提交到GenBank中,獲得的序列號(hào)為KC292921～KC292948。確定了每種魚D-loop的位置以及片段長(zhǎng)度。其中,鯉D-loop區(qū)段的長(zhǎng)度為926～929 bp,鯽為925 bp,草魚為928 bp(GenBank序列號(hào)為KC292921的1號(hào)草魚個(gè)體的69號(hào)堿基插入了一段長(zhǎng)度為231 bp的重復(fù)序列,其中包括草魚D-loop區(qū)段前部的203個(gè)堿基和tRNAPro尾部的28個(gè)堿基),鰱為933～938 bp,鳙為935～939 bp,烏蘇里擬鲿為893 bp。

2.3 6種魚序列的堿基組成與差異

2.3.1 堿基組成 用Mega 4.0計(jì)算出6種魚的D-loop區(qū)段及擴(kuò)增所得全部序列的堿基頻率(表1)。每種魚的A+T含量都明顯高于G+C含量,符合mtDNA的特點(diǎn),表明6種魚的堿基組成具有明顯偏向性[32-33]。

表1 6種魚D-loop區(qū)段及擴(kuò)增所得序列的堿基頻率Tab.1 Nucleotide frequencies in D-loop region and amplified fragment in six fish species %

2.3.2 種內(nèi)和種間的堿基差異 將6種魚進(jìn)行排序比較后發(fā)現(xiàn),6種魚序列結(jié)構(gòu)的差異性表現(xiàn)在種間和種內(nèi)兩方面,種間的差異性遠(yuǎn)大于種內(nèi)。而且6種魚的tRNAPro、tRNAPhe和12S rRNA三個(gè)區(qū)段的兩種差異性均遠(yuǎn)小于D-loop區(qū)段(表2)。不同種間在D-loop區(qū)段的起始端就表現(xiàn)出明顯的堿基差異,與其他5種魚親緣關(guān)系較遠(yuǎn)的烏蘇里擬鲿的差別性更加明顯。

表2 6種魚的擴(kuò)增片段中種內(nèi)和種間的堿基差異Tab.2 The base difference in amplified fragment between intraspecies and interspecies

2.3.3 PCR-RFLP酶切結(jié)果與序列分析的比較

使用相同引物對(duì)6種魚的mtDNAD-loop及其鄰近區(qū)段進(jìn)行了RFLP分析[15,22-25](烏蘇里擬鲿的結(jié)果未發(fā)表),在序列分析結(jié)果中均找到了RFLP的相應(yīng)酶切位點(diǎn),印證了RFLP實(shí)驗(yàn)結(jié)果。同時(shí)在酶切位點(diǎn)以外也發(fā)現(xiàn)了同種魚的不同單倍型間存在堿基差異(GenBank號(hào):KC292921～KC292948)。

2.4 6種魚的遺傳距離與系統(tǒng)發(fā)育樹

2.4.1 遺傳距離 基于 Kimura雙參數(shù)法(Kimura's 2-Parameter)計(jì)算6種魚種間D-loop區(qū)段的遺傳距離(表3)和擴(kuò)增所得全部序列的遺傳距離(表4)。其中,鳙與烏蘇里擬鲿的遺傳距離最大,為0.534668;鰱與鳙的遺傳距離最小,為0.106186。

表3 6種魚D-loop區(qū)段種間的平均Kimura's 2-Parameter遺傳距離Tab.3 The genetic distance of D-loop in six species based on Kimura's 2-Parameter model

表4 6種魚擴(kuò)增所得序列中種間的平均Kimura's 2-Parameter遺傳距離Tab.4 The genetic distance of amplified fragment in six species based on Kimura's-Parameter model

2.4.2 系統(tǒng)發(fā)育樹 運(yùn)用Mega 4.0軟件,基于Kimura雙參數(shù)模型,進(jìn)行1 000次抽樣重復(fù)檢測(cè),將對(duì)6種魚測(cè)序所得片段與GenBank中獲得的鯉序列(X61010.1)、金魚序列(AB111951.1)、鰱序列(EU315941.1)、草魚序列(EU391390.1)和長(zhǎng)吻鮠序列(GU596454.1)構(gòu)建NJ系統(tǒng)樹,如圖2所示。結(jié)果顯示,除鰱分支的支持率為85%外,其余均在90%以上,并可明顯分出6個(gè)分支。

圖2 基于Kimura's 2-Parameter模型構(gòu)建的PCR擴(kuò)增片段的NJ系統(tǒng)樹Fig.2 NJ phylogenetic tree of amplified fragment based on the Kimura's 2-Parameter model

3 討論

近年來(lái),對(duì)魚類mtDNA的研究逐步向序列分析的方向發(fā)展,序列遺傳多樣性分析及結(jié)構(gòu)研究已成為D-loop區(qū)的研究重點(diǎn)。本試驗(yàn)中運(yùn)用軟件對(duì)6種魚的mtDNAtRNAPro、D-loop、tRNAPhe和部分12S rRNA序列進(jìn)行分析,得到6種魚堿基含量、堿基差異、遺傳距離和NJ系統(tǒng)樹等相關(guān)數(shù)據(jù)。

與GenBank中5種魚的同源序列比對(duì)后,確定了6種魚的PCR產(chǎn)物中所包含的各個(gè)區(qū)段以及mtDNA D-loop區(qū)的位置和片段長(zhǎng)度。發(fā)現(xiàn)6種魚D-loop區(qū)的長(zhǎng)度在種內(nèi)和種間都存在差異。吳海防等[25]報(bào)道的草魚單倍型Ⅱ比一般草魚長(zhǎng)出一段DNA為231 bp的重復(fù)序列,基于mtDNA D-loop區(qū)段用PCR-RFLP法檢測(cè)出的草魚兩種單倍型的差異源于此重復(fù)序列。

本研究中選取了除鯽之外的5種魚mtDNA D-loop及其鄰近區(qū)段經(jīng)PCR-RFLP法檢測(cè)的單倍型,經(jīng)過(guò)序列分析,印證了每種魚的不同單倍型在限制酶酶切位點(diǎn)上的差異性,并且除酶切位點(diǎn)以外也存在堿基差異。經(jīng)過(guò)酶切位點(diǎn)比對(duì),發(fā)現(xiàn)單酶切之后會(huì)出現(xiàn)一些較小的酶切片段,這些片段在使用瓊脂糖凝膠電泳分辨時(shí)受到限制,因而可以認(rèn)為,在很多用PCR-RFLP檢測(cè)分析的報(bào)道中酶切片段長(zhǎng)度之和會(huì)小于 PCR 擴(kuò)增的片段長(zhǎng)度[13-14,23-25]。與PCR-RFLP法相比,序列分析的方法精確度更高。

6種魚序列差異性的研究證明了mtDNA的D-loop區(qū)段的變異程度比其他序列高。在6種魚的NJ系統(tǒng)樹中,親緣關(guān)系較近的鰱和鳙,鯉和鯽都可以明顯的進(jìn)行區(qū)分,并且具有較高的支持率。本研究結(jié)果表明,D-loop區(qū)段或者D-loop區(qū)段結(jié)合其鄰近區(qū)段適合系統(tǒng)發(fā)育研究,并且可以為魚類種類識(shí)別以及遺傳多樣性分析等方面提供科學(xué)依據(jù)。

[1]肖武漢,張亞平.魚類線粒體DNA的遺傳與進(jìn)化[J].水生生物學(xué)報(bào),2000,24(4):384-391.

[2]楊明生,熊邦喜.動(dòng)物mtDNA的研究及在魚類生態(tài)學(xué)中的應(yīng)用[J].孝感學(xué)院學(xué)報(bào),2005,25(3):23-27.

[3]汪泰初,劉朝良,肖林珍.線粒體基因組(mtDNA)的研究進(jìn)展[J].安徽農(nóng)業(yè)科學(xué),2006,34(10):2068,2071.

[4]潘寶平,卜文俊.線粒體基因組的遺傳與進(jìn)化研究進(jìn)展[J].生物學(xué)通報(bào),2005,40(8):1-3.

[5]鄒習(xí)俊,韓雪,韓虎峰.魚類mtDNA及其非編碼區(qū)的研究概況[J].貴州畜牧獸醫(yī),2009,33(3):23-25.

[6]徐鋼春,魏廣蓮,李建林,等.基于線粒體DNA D-loop序列分析養(yǎng)殖刀鱭與湖鱭的遺傳多樣性[J].大連海洋大學(xué)學(xué)報(bào), 2012,27(5):448-452.

[7]朱雪蓮,王志勇,韓志強(qiáng).借助mtDNA控制區(qū)序列分析金魚與不同地域鯽的親緣關(guān)系[J].大連水產(chǎn)學(xué)院學(xué)報(bào),2010,25(2): 153-157.

[8]黃立群,儀慧蘭,崔松林,等.歷山中國(guó)大鯢線粒體片段序列的測(cè)定及其遺傳差異研究[J].大連海洋大學(xué)學(xué)報(bào),2012,27(6): 513-517.

[9]程起群,馬春艷,繆炯,等.基于線粒體12Sr RNA序列研究鳳鱭兩個(gè)野生群體的遺傳多樣性[J].大連水產(chǎn)學(xué)院學(xué)報(bào),2007, 22(5):387-391.

[10]Rosel P E,Dizon A E,Haygood M G.Variability of the mitochondrial control region in population of the harbour porpoise,Phocoena phocoena,on interoceanic and regional scales[J].Can J Fish Aquat Sci,1995,52:1421-1429.

[11]Rosel P E,Haygood M G,Perrin W F.Phylogenetic relationship among the true porpoise(Cetacean:Phocoenidae)[J].Molecular Phylogenetics and Evolution,1995,4(4):463-474.

[12]Zhu D,Jamieson B G,Hugall A,et al.Sequence evolution and phylogenetic signal in control region and cytochrome b sequences of rainbow fishes(Melanotaeniidae)[J].Mol Biol Evol,1994,11: 672-683.

[13]Sugaya T,Ikeda M,Taniguchi N.Typing of mtDNA in the kuruma prawn based on the PCR-RFLP analysis of mtDNA control region [J].Fish Genetics and Breeding Science,2001,31:11-15.

[14]Ikeda M,Taniguchi N.Genetic variation and divergence in populations of ayu Plecoglossus altivelis,including endangered subspecies,inferred from PCR-RFLP analysis of the mitochondrial DNA D-loop region[J].Fisheries Science,2002,68:18-26.

[15]梁愛軍,王淞,鮑迪,等.烏克蘭鱗鯉同工酶及mtDNA D-loop區(qū)段的RFLP分析[J].安徽農(nóng)業(yè)科學(xué),2011,39(27):17109-17113.

[16]李思發(fā),吳力釗,王強(qiáng),等.長(zhǎng)江、珠江、黑龍江鰱、鳙、草魚種質(zhì)資源研究[M].上海:上?？茖W(xué)技術(shù)出版社,1990:28-39.

[17]宋平,李小迎,熊全沫.鰱鳙線粒體DNA的九種限制性內(nèi)切酶酶切圖譜的比較[J].水產(chǎn)學(xué)報(bào),1994,18(3):221-230.

[18]張錫元,張德春,楊代淑,等.長(zhǎng)江鰱遺傳多樣性的隨機(jī)擴(kuò)增多態(tài)DNA分析[J].水產(chǎn)學(xué)報(bào),1999,23(增刊):7-14.

[19]沈玉幫,張俊彬,李家樂(lè).草魚種質(zhì)資源研究進(jìn)展[J].中國(guó)農(nóng)學(xué)通報(bào),2011,27(7):369-373.

[20]潘偉志,陳軍,王鵬,等.烏蘇里擬鲿人工繁育技術(shù)研究[J].水產(chǎn)學(xué)雜志,2006,19(1):31-35.

[21]陳軍,潘偉志,趙春剛,等.烏蘇里擬鲿苗種培養(yǎng)技術(shù)[J].黑龍江水產(chǎn),2006(3):13-17.

[22]劉琳.應(yīng)用DNA及同工酶遺傳標(biāo)記檢測(cè)分析二、三倍體鯽魚的遺傳結(jié)構(gòu)[D].天津:天津師范大學(xué),2008.

[23]王淞,曹曉霞,谷口順彥,等.4個(gè)群體鰱mtDNA D-loop的PCR-RFLP分析[J].淡水漁業(yè),2010,40(4):3-9.

[24]單淇,董仕,吳海防,等.3個(gè)群體鳙魚mtDNA D-loop區(qū)段的限制性片段長(zhǎng)度多態(tài)性分析[J].中國(guó)水產(chǎn)科學(xué),2006,13(2): 174-180.

[25]吳海防,董仕,單淇,等.3個(gè)群體草魚mtDNA D-Loop的PCR-RFLP分析[J].水產(chǎn)科學(xué),2006,25(4):184-188.

[26]Perez-Enriquze R,Taniguchi N.Genetic structure of red sea bream(Pagrus major)population of Japan and the Southwest Pacific,using microsatellite DNA markers[J].Fisheries Science, 1999,65(1):23-30.

[27]Iguchi K,Tanimura Y,Nishida M.Sequence divergence in the mtDNA control region of amphidromous and landlocked forms of ayu[J].Fisheries Science,1997,63:901-905.

[28]Takahashi Y,Kobayashi T,Ueda T.Genetical variation of the D-loop region in mitochondrial DNA using PCR method in metropolitan bitterling[J].Fish Genetics and Breeding Science,1994,20: 39-45.

[29]Thompson J D,Gibson T J,Plewniak 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.

[30]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.

[31]李彬彬,黃培春,鐘復(fù)光.生物學(xué)軟件在線粒體DNA序列多態(tài)性分析中的應(yīng)用[J].生物信息學(xué),2010,8(2):153-155.

[32]張方,米志勇,毛鐘榮,等.草魚(Ctenopharyngodon idellus)線粒體DNA控制區(qū)及其兩側(cè)tRNA基因的克隆與結(jié)構(gòu)分析[J].中國(guó)生物化學(xué)與分子生物學(xué)報(bào),1999,15(3):417-422.

[33]Broughton R E,Roe B A.The complete sequence of the zebra fish(Danio rerio)mitochondrial genome and evolutionary patterns in vertebrate mitochondrial DNA[J].Genome Res,2001,11(11): 1958-1967.

The sequence comparison of mtDNA D-loop and adjacent regions in six fish species

HAO Jun,YANG Qiang,BAO Di,LIANG Ai-jun,ZHANG Xing-hua,DONG Shi

(College of Life Science,Tianjin Normal University,Tianjin Key Laboratory of Animal and Plant Resistance,Tianjin 300387,China)

The mtDNA D-loop and its adjacent regions were amplified and sequenced in 28 individuals from different haplotypes detected using PCR-RFLP in Ukraine scale carp(Cyprinus carpio),crucian carp(Carassius auratus),silver carp(Hypophthalmichthys molitrix),big head carp(Aristichthys nobilis),grass carp(Ctenopharyngodon idellus)and Ussuri bullhead(Pseudobgrus ussuriensis)by PCR.The amplified fragment length of mtDNA D-loop and its adjacent regions in the fish were varied from 1 500 bp to 1 800 bp,and are compared with the mtDNA complete sequences of five fish species from GenBank by Clustal 1.83 software.It was found that the amplified fragment was comprised of complete tRNAPro,D-loop,tRNAPhesequences and about 400 bp part of the 12S rRNA sequence.Sequences of 28 individuals were submitted to GenBank for serial numbers(KC292921-KC292948), and the base composition and difference of the six species,and the genetic distance were calculated based on Kimura's 2-Parameter model and NJ phylogenetic tree with five sequences from GenBank was constructed by Mega 4.0 software.The sequence structure analysis showed that the intraspecific and interspecific differences in D-loop region were higher in the four region sequences of six species than those in other species.NJ phylogenetic trees of six species were consistent with the traditional classification methods,showing useful for fish species classification and identification.In addition,restriction sites of all kinds of haplotypes were found in sequences,which is consistent with PCR-RFLP experimental results.More variation out of the restriction sites were observed compared with RFLP analysis.

Ukraine scale carp;crucian carp;silver carp;big head carp;grass carp;Ussuri bullhead;D-loop; NJ phylogenetic tree

S965.113;Q953

A

2012-12-13

天津市動(dòng)植物抗性重點(diǎn)實(shí)驗(yàn)室開放研究基金資助項(xiàng)目

郝君(1987-),女,碩士研究生。E-mail:haohaojunjun2645@163.com

董仕(1959-),男,博士,研究員。E-mail:skyds@mail.tjnu.edu.cn

2095-1388(2013)02-0160-06