舟山帶魚線粒體COI基因SNP位點分析

卞光明，王娜泠，胡則輝，王躍斌，胡成碩，柴學軍

（浙江海洋大學海洋與漁業研究所，浙江省海洋水產研究所，浙江舟山 316021)

舟山帶魚線粒體COI基因SNP位點分析

卞光明，王娜泠，胡則輝，王躍斌，胡成碩，柴學軍

（浙江海洋大學海洋與漁業研究所，浙江省海洋水產研究所，浙江舟山 316021)

為探討舟山帶魚種質狀況，本研究利用直接測序法對舟山帶魚線粒體COI基因進行SNP挖掘分析，獲得長度一致的985 bp基因片段，檢測到22個SNP位點，對SNP位點在樣品中分布情況分析，發現部分SNP位點分布情況類似甚至完全相同，推測基于線粒體COI基因片段不同位點突變可能存在關聯性。

帶魚；線粒體COI基因；SNP；關聯性

舟山帶魚Trichiurus lepturus是舟山漁場最重要的海洋捕撈經濟魚類，是全國首批海鮮類地理標志證明商標海鮮，其眼睛呈黑色、骨小體肥、背脊無凸骨、鱗片易脫落、肉質鮮而不腥，與其他帶魚相比，具有DHA成分、無大骨粒異物感，被稱為“世界上最好吃的帶魚”[1]。由于舟山帶魚捕撈強度過度、人工繁育技術難度較大以及生存環境變化，舟山漁場甚至出現無帶魚可撈現象，隨著休漁、種質資源保護區、產出控制等措施實施，帶魚資源狀況雖有所好轉，但群體小型化、年齡結構失衡、產卵期延長以及單位漁獲量下降等現象仍較為嚴峻[2]。

線粒體DNA（mitochondrial DNA，mtDNA）嚴格遵守母系遺傳、分子結構簡單穩定，具有可溯源性等，常應用于進化生物學及群體遺傳學研究。細胞色素C氧化酶I亞基（cytochrome c oxidase subunit I,COI）攜帶的遺傳信息量大、進化速率適中、中間解析能力強，已普遍應用于銀鯧Pampus argenteus[3]、黃姑魚屬Nibea[4]、四指馬鲅屬 Eleutheronoma[5]、太湖新銀魚 Neosalanx taihuensis[6]、脊尾白蝦 Exopalaemon carinicauda[7]、中華絨螯蟹Eriocheir sinensis[8]、大眼蟹屬Macrophthalmus[9]、紫貽貝Mytilus galloprovincialis[10]以及厚殼貽貝Mytilus coruscus[11]等水生動物種群遺傳結構分析、分子進化及系統發育。SNP分布廣、覆蓋率高、遺傳穩定、多態性豐富，直接測序法操作簡單、技術成熟、無需通過多個位點構建遺傳圖譜，特別適合魚類等較小目的基因片段研究[12]。目前所報道的研究主要采用 16S rRNA[13]、18S rDNA[14]、D-loop[15-16]、COI[17]、SSR[18]分析不同海域帶魚群體，未見具體分析舟山帶魚線粒體COI基因SNP位點特性研究。本研究利用直接測序法對舟山帶魚的線粒體COI基因區域進行SNP位點挖掘和分析，以期為舟山帶魚遺傳學研究、種質資源的研究提供基礎資料。

1 材料與方法

1.1 實驗材料

實驗所用60個帶魚樣品由雙拖船取自舟山朱家尖附近海域（122°47′E，29°97′N），取樣時間為2016年4月中旬，現場取樣并剪下背部肌肉于95%酒精中保存，帶回實驗室置于-20℃保存待用。

1.2 實驗方法

1.2.1 DNA提取

按照天根海洋動物基因組提取試劑盒操作說明提取總基因組DNA，經0.8%瓊脂糖凝膠電泳檢測DNA完整性后，用Nanodrop 2000紫外分光光度計檢測其質量與濃度，稀釋為50 ng/μL，-20℃保存待用。

1.2.2 引物設計

從GenBank中下載帶魚線粒體全基因組序列[19]，用Oligo7軟件[20]設計引物擴增目的片段，目的片段大小在1 200 bp左右，引物序列為COI-F：CTTATCCGAGCAGAACTAAG，COI-R：ACGGCTCCTATTGATAGAAC，退火溫度55℃左右，由上海生工生物工程股份有限公司合成。

1.2.3 PCR擴增及測序

PCR 反應總體系為 50 μL，模板 DNA 100 ng，0.5 μL 上、下游引物（10 μmol/L），1 μL dNTP（10 mmol/μL），5 μL 10×Buffer（含Mg2+），0.5 μL Taq DNA聚合酶（5 U/μL），去離子水補充至50 μL。PCR反應條件：95℃預變性3 min，然后進行35個循環，每個循環包括94℃ 30 s，55℃ 1 min，72℃ 1 min，最后72℃延伸7 min。PCR產物經1.4%瓊脂糖凝膠電泳檢測后，委托上海生工生物工程股份有限公司進行雙向測序。

1.3 數據分析

運用ContigExpress Project[21]對測序結果進行拼接，再使用BioEdit[22]對所得序列進行比對篩選，利用Chromas 2.4.1分析序列峰圖，以測序波峰整齊，沒有雜峰，突變堿基位點樣品數不少于2個的位點記作有效突變位點，以排除測序機器誤讀或錯讀產生的變異位點。

2 結果與分析

2.1 PCR瓊脂糖凝膠電泳圖

PCR擴增產物凝膠電泳圖如圖1所示，Marker為2 000 bp，從下至上依次為（100 bp、250 bp、500 bp、750 bp、1 000 bp、2 000 bp），從電泳圖上可觀察到，產物片段的大小與設計的目的片段的長度一致（1 200 bp），且條帶明亮、清晰，沒有雜帶，滿足直接測序的要求。

2.2 SNP位點峰圖

圖1 部分PCR產物電泳圖譜Fig.1 Electrophoresis results of part PCR product

本文對三個突變頻率大于25%的SNP位點進行峰圖分析（23位點、500位點、773位點），位點具體突變形式如圖2～4所示，三個SNP位點波峰整齊、無雜峰。

SNP23突變形式：C/T轉換

SNP500突變形式：T/C轉換

SNP 773突變形式：G/A轉換

圖2 COI基因23位點突變形式Fig.2 Mutation forms of COI gene at 23rdlocus

圖3 COI基因500位點突變形式Fig.3 Mutation forms of COI gene at 500thlocus

圖4 COI基因733位點突變形式Fig.4 Mutation forms of COI gene at 733rdlocus

2.3 SNP位點統計分析

使用BioEdit軟件對測序序列比對分析，刪除兩端不穩定序列，獲得985 bp基因片段，共檢測到22個變異位點，檢測到SNP位點具體情況見表1，773位點、23位點、500位點以及953位點突變率較高，分別為43.67%、25.00%、25.00% 以及23.33%。

表1 COI基因序列部分區域SNP位點Tab.1 SNPs identified in COI gene subregio

表2為SNP位點在60個帶魚樣品中的具體分布情況，從表中可以發現，815位點、857位點、968位點以及983位點發生突變可能具有一定的關聯性，位點突變率均為18.33%，在11個樣品中均檢測到且分布情況完全一致；位點119、位點404以及位點977的位點突變率為6.67%；位點353與位點413的位點突變率為6.67%，在樣品中的分布情況也呈現一致，雖然突變數量較少（位點突變率僅為6.67%），但不排除相關位點之間的突變存在關聯性。除此以外，440位點、500位點以及953位點在樣品中的分布情況相似度也較高。

表2 SNP位點在樣品中的分布情況Tab.2 The distribution of SNP locus in samples

續表

3 討論

本文采用直接測序法對帶魚COI基因進行測序，有效避免傳統基因克隆方法對變異位點分型產生誤讀、錯讀等影響。在獲得的985 bp片段中檢測到22個變異位點，變異率為0.021 3，與帶魚其他相關研究結果相近[17]。根據Chromas 2.4.1分析序列峰圖所示，變異位點峰圖清晰、無雜帶，能夠對SNP位點進行精準的統計。對線粒體COI基因片段SNP位點在60個舟山帶魚樣品中的分布情況統計的過程中，發現一些SNP位點變異頻率以及分布情況存在相似性甚至完全相同，關于不同位點變異是否存在關聯性，還有待于進一步研究。

[1]舟山帶魚[DB/OL].(2012-12-27)[2016-08-15].http://www.wangdanian.cn/wenda/5324.html.

[2]浙江省海洋水產研究所資源室.2015年上半年浙江海洋漁業資源動態評析[DB/OL].(2015-09-11)[2016-08-15].http://www.zjhys.cn/news/hydt/20159/1111331020951901.html.

[3]彭士明,施兆鴻,侯俊利.基于線粒體D-loop區與COI基因序列比較分析養殖與野生銀鯧群體遺傳多樣性[J].水產學報,2010,34(1):19-25.

[4]CHAI X J,HU Z H,WANG Y B,et al.Genetic diversity and the relationship of Nibea japonica and Nibea miichthiodies using mitochondrial DNA gene[J].Biochemical Systematics and Ecology,2013,48(2):274-280.

[5]鄧春興.基于COI基因序列的中國東南沿海四指馬鲅屬魚類的遺傳多樣性分析[D].廣州:暨南大學,2014:23-45.

[6]張 迪,雷光春,龔 成,等.基于COI基因序列的太湖新銀魚遺傳多樣性[J].湖泊科學,2012,24(2):299-306.

[7]馬 朋,劉 萍,李 健,等.脊尾白蝦3個野生群體線粒體COI基因的遺傳多樣性及其系統發育分析[J].漁業科學進展,2011,32(6):50-56.

[8]葛家春,許志強,李曉暉,等.利用線粒體COI序列分析4水系中華絨螯蟹群體遺傳學特征[J].中國水產科學,2011,18(1):16-22.

[9]徐敬明.2種大眼蟹線粒體COI基因序列的比較研究[J].水生態學雜志,2010,3(2):42-47.

[10]沈玉幫,張俊彬,馮冰冰,等.基于線粒體COI序列分析對紫貽貝群體遺傳多樣性的研究分析[J].海洋通報,2011,30(4):435-440.

[11]楊振雄,毛陽麗,宋 娜,等.浙江和福建沿海厚殼貽貝Mytilus coruscus群體的COI序列比較分析[J].海洋湖沼通報,2014(2):82-87.

[12]柏盈盈,賈智英,葛會爭,等.TLR9基因SNP位點檢測及其在7種鯉中的遺傳多樣性研究[J].淡水漁業,2013,43(6):19-28.

[13]蒙子寧.中國近海小黃魚、帶魚和小帶魚的遺傳多樣性以及八種石首魚類的分子系統進化[D].廈門:廈門大學,2003:91-108.

[14]張繼民,劉 霜,趙建民,等.帶魚（Trichiurus haumela）魚卵DNA的提取及其18S rDNA初步分析[J].海洋通報,2009,28(6):62-65.

[15]肖永雙.西北太平洋五種海洋魚類的分子系統地理學研究[D].青島:中國海洋大學,2011:101-116.

[16]鄭文娟,杜一超,林 潔,等.基于線粒體DNA D-loop區部分序列分析舟山海域帶魚種群遺傳結構[J].水生生物學報,2015,39(2):408-423.

[17]李鵬飛,朱文斌,賀舟挺,等.東海帶魚DNA條形碼的建立及COI序列變異分析[J].浙江海洋學院學報:自然科學版,2013,32(1):6-9.

[18]畢金貞.帶魚微衛星標記的開發和微衛星標記在牙鲆選擇育種中的應用[D].青島:中國海洋大學,2010:6-40.

[19]LIU Y,CUI Z.The complete mitochondrial genome sequence of the cutlassfish Trichiurus japonicus(Perciformes:Trichiuridae):Genome characterization and phylogenetic considerations[J].Marine Genomics,2009,2(2):133-142.

[20]RYCHLIK W.OLIGO 7 primer analysis software[J].Methods in Molecular Biology,2007,402:35-60.

[21]CLARK M J,CHEN R,LAM H Y,et al.Performance comparison of exome DNA sequencing technologies[J].Nature Biotechnol,2011,29(10):908-914.

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

Analysis on Mitochondrial COI Gene SNP Loci of Trichiurus lepturus from Zhoushan Sea Area

BIAN Guang-ming,WANG Na-ling,HU Ze-hui,et al
(Marine and Fisheries Institute of Zhejiang Ocean University,Zhejiang Marine Fisheries Institute,Zhoushan 316021,China)

To investigate the germplasm resources of Trichiurus lepturus collected from Zhoushan coastal sea area,SNP loci and analysis based on mtDNA COI gene were carried out using direct sequencing technique.The fragments with the length of 985 bp were obtained,detecting 22 SNP locus.According to the similar or even identical distribution of SNP locus in samples,there may be a correlation between mutations of mtDNA COI gene fragments at different sites.

Trichiurus lepturus;mitochondrial COI gene;SNP;correlation

Q953+.3

：A

2016-11-10

浙江省科技廳合作與轉化項目（2013F50003）；舟山市公益項目（2015C31013）

卞光明（1992-），男，江蘇連云港人，碩士研究生，研究方向：種質與種苗工程.E-mail:mingbg@163.com

柴學軍（1972-），男，教授級高工，研究方向：魚類繁育與海洋生物技術.E-mail:chaixj6530@sohu.com

1008－830X(2017)01－0019－06