鄭真真，許強(qiáng)華、2，戴小杰、2，朱江峰、2

(1.上海海洋大學(xué) 海洋科學(xué)學(xué)院，上海201306;2.大洋漁業(yè)資源可持續(xù)開發(fā)省部共建教育部重點(diǎn)實驗室，農(nóng)業(yè)部大洋漁業(yè)資源環(huán)境科學(xué)觀測實驗站，國家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心，上海201306)

大青鯊Prionace glauca是一種大型中上層鯊魚，隸屬于真鯊目Carcharhiniformes、真鯊科Carcharhinidae、大青鯊屬Prionace Cantor，是3種最豐富的大洋性鯊魚之一［1］。作為一種大型的軟骨魚類，大青鯊廣泛分布于太平洋、大西洋、印度洋(簡稱三大洋)海域，其分布范圍最北可達(dá)挪威，最南可達(dá)智利。大青鯊在太平洋主要分布于20°N ～20°S，在大西洋主要分布于從紐芬蘭到阿根廷的東部海域，以及從挪威到南非的西部海域［2-3］。根據(jù)大西洋金槍魚養(yǎng)護(hù)委員會(ICCAT)1990—2010年的統(tǒng)計數(shù)據(jù)，大青鯊在金槍魚的兼捕種類中占的比重最大，大青鯊數(shù)量正在急劇下降，已引起了ICCAT 的重點(diǎn)關(guān)注［4］。結(jié)合當(dāng)前的捕撈情況和資源現(xiàn)狀，保護(hù)大青鯊種群已成為漁業(yè)管理中一項重要任務(wù)。

線粒體DNA(mitochondrial DNA，mtDNA)呈共價閉合的環(huán)狀雙鏈，具有分子小、結(jié)構(gòu)簡單、嚴(yán)格遵守母系遺傳的特點(diǎn)，能有效地揭示群體內(nèi)外的遺傳變異，現(xiàn)已成為魚類分子生態(tài)學(xué)、遺傳多樣性研究的重要標(biāo)記［5-6］。線粒體控制區(qū)(D - loop)具有較高的突變和累積，是研究種內(nèi)遺傳分化的重要工具，被廣泛應(yīng)用于各物種的遺傳分化和遺傳結(jié)構(gòu)研究［7］。陳驍?shù)龋?］利用線粒體控制區(qū)對中國南部海域條紋斑竹鯊Chiloscyllium plagiosum 的遺傳多態(tài)性進(jìn)行了研究，結(jié)果表明，條紋斑竹鯊的基因流主要在淺海近岸水域擴(kuò)散，遺傳變異程度受地理結(jié)構(gòu)和距離隔離的影響;Mendonca等［9］采用COⅠ基因部分序列和D-loop 部分序列評估大西洋西部斜鋸牙鯊Rhizoprionodon acutus種群核苷酸的差異，結(jié)果表明，大西洋斜鋸牙鯊和加勒比斜鋸牙鯊兩個不同物種之間存在較高水平的遺傳差異;Keeney等［10］利用線粒體控制區(qū)對大西洋墨西哥灣169 個黑鰭真鯊Carcharhinus limbatus 幼體樣本的遺傳異質(zhì)性進(jìn)行了相關(guān)研究，結(jié)果表明，該海域的黑鰭真鯊并沒有出現(xiàn)明顯的種群遺傳結(jié)構(gòu)。

近年來，不少學(xué)者對大青鯊的繁殖生物學(xué)特征［11］、資源狀況［12-15］、行為特性［16-18］、棲息環(huán)境［19］和生長［4］等進(jìn)行了研究，為大青鯊群體的資源保護(hù)和可持續(xù)發(fā)展提供了科學(xué)資料。戴小杰等［20］對東太平洋熱帶海域大青鯊繁殖生物學(xué)特征的研究表明，大青鯊是繁殖能力很強(qiáng)的大型鯊魚之一。大青鯊為季節(jié)性的高度洄游性種類，對其的標(biāo)記回捕記錄和漁獲數(shù)據(jù)表明，南北大西洋大青鯊為同一種群。利用線粒體DNA 非編碼區(qū)序列和微衛(wèi)星標(biāo)記分析顯示，北大西洋東西部的大青鯊種群也無顯著性差異［21-23］。但迄今為止，全球范圍內(nèi)對大青鯊的遺傳多樣性以及群體遺傳結(jié)構(gòu)等方面的研究尚未見報道。本研究中，以三大洋中大青鯊的5個采樣群體為研究對象，通過分析其mtDNA 的D-loop 區(qū)部分序列，分析其遺傳多樣性和種群遺傳結(jié)構(gòu)，并對其地理分布進(jìn)行初步探討，以期為全球大青鯊資源的合理開發(fā)、管理和保護(hù)提供科學(xué)依據(jù)。

1 材料與方法

1.1 材料

試驗用165 尾大青鯊成魚的肌肉組織樣本取自2010年9月—2012年2月中國科學(xué)觀察員在三大洋公海海域金槍魚延繩釣漁船的漁獲物。該船作業(yè)期間在赤道附近海域有5 個采樣位點(diǎn)，分別為中東太平洋(Central East Pacific，CEP)、中西太平洋(Central West Pacific，CWP)、西南大西洋(Southwest Atlantic，SWA)、中東大西洋(Central East Atlantic，CEA)和印度洋(Indian Ocean，INO)，采樣位點(diǎn)及采樣時間見表1。剪取大青鯊肌肉組織(約20 g)浸入無水乙醇中，于冰箱(-20 ℃)中保存，帶回實驗室后保存于冰箱(-80 ℃)中備用。

1.2 方法

1.2.1 基因組DNA 的提取 取每尾大青鯊50 ～100 mg 肌肉組織剪碎，使用組織基因組DNA 快速提取試劑盒(北京艾德萊生物科技有限公司)提取大青鯊基因組DNA。用10 g/L 瓊脂糖凝膠電泳檢測基因組DNA 的完整性，用紫外分光光度計檢測提取的DNA 濃度和純度。

表1 大青鯊的采樣位點(diǎn)、采樣時間Tab.1 The sampling locations and date for blue shark Prionace glauca

1.2.2 D-loop 區(qū)基因的PCR 擴(kuò)增和測序 以大青鯊基因組DNA 為模板，隨機(jī)抽取3 個大青鯊樣本直接采用引物5' TTGGCTCCCAAAGCCAARATTCTG 3'(正向)、5'GTGGCTTAGCAAGARGTCTTG 3'(反向)(上海杰李生物技術(shù)有限公司)，擴(kuò)增出的D-loop 區(qū)基因片段長約1100 bp，PCR 擴(kuò)增體系(共20 μL):Taq 聚合酶0.2 μL，10 ×buffer 2.0 μL，dNTP 1.6 μL，正反向引物各1 μL(10 pmol/μL)，模板1.0 μL，用雙蒸水ddH2O 補(bǔ)足至20 μL。PCR 反應(yīng)條件:95 ℃下預(yù)變性7 min;95℃下變性45 s，56 ℃下退火60 s，72 ℃下延伸90 s，共進(jìn)行35 個循環(huán);最后在72 ℃下再延伸7 min，4 ℃下保存。PCR 擴(kuò)增產(chǎn)物用10 g/L 瓊脂糖凝膠電泳分離，得到目的片段，送交生工生物工程(上海)股份有限公司測序。根據(jù)3 個樣本的測序結(jié)果，自行設(shè)計另外一對引物:PGCR -F，5' TGCACGTTAGTACCACGCTAA 3';PGCR - R，5'AACTGGCCCTACATAACCTGC 3'。擴(kuò)增出的D -loop 基因片段長約700 bp，PCR 擴(kuò)增體系同上。PCR 反應(yīng)條件:95 ℃下預(yù)變性3 min;95 ℃下變性30 s，55 ℃下退火30 s，72 ℃下延伸60 s，共進(jìn)行30 個循環(huán);最后在72 ℃下再延伸10 min，4℃下保存。送交生工生物工程(上海)股份有限公司測序，得到165 個大青鯊樣本的序列片段。

1.2.3 數(shù)據(jù)分析 測序結(jié)果先經(jīng)Blast 程序進(jìn)行序列比對檢索，結(jié)合DNA Star 軟件檢測其同源性，使用BioEdit 軟件進(jìn)行控制區(qū)序列的對位排列(alignment)［24］，并對序列輔以人工校對、編輯和比對，確保序列無誤。采用Mega 4.0 軟件［25］進(jìn)行序列分析，包括序列長度、堿基組成和親緣關(guān)系樹的構(gòu)建，大青鯊單倍型系統(tǒng)發(fā)育樹的構(gòu)建采用鄰接法(Neighbor -joining，NJ)［26-27］，序列間的遺傳距離以Kimura 雙參數(shù)法(Kimura 2 -parameter)［28］進(jìn)行估計。系統(tǒng)發(fā)育樹分支的置信度采用自引導(dǎo)法(bootstrap analysis，BP)［29］重復(fù)檢測(1000次重復(fù))。應(yīng)用DnaSP 4.0 軟件［30］計算單倍型多樣性指數(shù)(Hd)與核苷酸多樣性指數(shù)(π)。應(yīng)用Arlequin 3.0 軟件［31］進(jìn)行分子方差分析(AMOVA)，其中，采用單倍型頻率計算Fst值，根據(jù)公式M=1/4(1/Fst-1)計算基因流Nm。

2 結(jié)果與分析

2.1 mtDNA D－loop 區(qū)基因序列特征

經(jīng)PCR 擴(kuò)增后，獲得165 個長度為694 bp 的線粒體控制區(qū)(D-loop 區(qū))的堿基序列，序列中各堿基在DNA 雙鏈中平均所占的比例:A(腺嘌呤)為33.46%，T(胸腺嘧啶)為36.40%，C(胞嘧啶)為18.61%，G(鳥嘌呤)為11.53%，G+C 含量為30.14%，僅占A+T 含量的二分之一，大青鯊mtDNA D -loop 區(qū)基因片段中G 含量都比較低，平均含量為11.53%(表2)，明顯低于其他3種堿基的含量，表現(xiàn)出顯著的堿基組成偏向性。

表2 大青鯊mtDNA D－loop 區(qū)基因序列的堿基組成Tab.2 Nucleotide compositions of control region sequence of mitochondrial DNA in blue shark

2.2 遺傳多樣性和單倍型分布

在165 尾大青鯊mtDNA D-loop 區(qū)域的694 bp序列中共檢測到110 個變異位點(diǎn)，變異率為15.85%，定義145 個單倍型，分別為H1 ～H145。5 個海域的大青鯊群體mtDNA D -loop 統(tǒng)計結(jié)果顯示，165 尾個體的單倍型多樣性指數(shù)為0.997 3 ±0.001 4，核苷酸多樣性指數(shù)為0.015 10 ±0.000 91(表3)。165 尾大青鯊個體中共檢測到145種單倍型。除了CEA 群體擁有自己特有的單倍型，其他群體都擁有共享單倍型。SWA和CWP 群體同時擁有單倍型H24，而單倍型H66、H71、H76、H78、H79 同時出現(xiàn)在CWP和CEP 群體，單倍型H79、H92 同時出現(xiàn)在CEP和INO 群體(表4)。這些大量存在的共享單倍型表明，全球大青鯊群體間存在較強(qiáng)的基因交流。

表3 大青鯊群體線粒體DNA 控制區(qū)序列的遺傳多樣性參數(shù)Tab.3 Genetic diversity parameters of blue shark based on control region sequence data of mitochondrial DNA

表4 各采樣位點(diǎn)的大青鯊單倍型分布Tab.4 Haplotype distribution of blue shark sampled at various localities

2.3 系統(tǒng)發(fā)育樹

對大青鯊145 個單倍型采用鄰接法構(gòu)建系統(tǒng)發(fā)育樹。結(jié)果顯示:各單倍型廣泛分布在單倍型鄰接樹上，未形成明顯的分支，不存在顯著的譜系結(jié)構(gòu)?；诰€粒體DNA D-loop區(qū)構(gòu)建的單倍型間系譜關(guān)系顯示:145種單倍型只形成單一的網(wǎng)絡(luò)，各單倍型間呈現(xiàn)高水平的平行演化(由于145 個單倍型未形成明顯分支，不存在顯著譜系結(jié)構(gòu)，故圖略)。

2.4 群體遺傳結(jié)構(gòu)分析

用分子方差分析(AMOVA)法對群體間的遺傳結(jié)構(gòu)進(jìn)行分析，結(jié)果顯示，99.28%的遺傳變異出現(xiàn)在種群內(nèi)，而0.72%的遺傳變異出現(xiàn)在種群間;群體間的遺傳分化指數(shù)Fst為0.007 21(表5)，說明群體間不存在顯著的遺傳分化。

表5 大青鯊種群的分子方差分析Tab.5 AMOVA analysis of blue shark populations

利用Fst進(jìn)一步檢驗大青鯊群體間的遺傳結(jié)構(gòu)，根據(jù)兩兩群體間的Fst值可以發(fā)現(xiàn)，所有的Fst值均較小，除SWA 群體與CEP 群體間統(tǒng)計檢驗極顯著(P=0.009＜0.01)外，其他任何兩個群體間的統(tǒng)計檢驗均不顯著(P ＞0.05)。本研究結(jié)果表明，除西南大西洋與中東太平洋的大青鯊群體間存在著顯著的遺傳分化外，其他采樣區(qū)域內(nèi)的大青鯊群體均不存在顯著的遺傳結(jié)構(gòu)。除此之外，5 個采樣點(diǎn)間大青鯊的基因交流(Nm)也比較頻繁，CWP 群體與CEA 群體間Nm最小，為11.83(表6)。

表6 采樣區(qū)域大青鯊種群間的遺傳分化分析及基因流NmTab.6 Fst analysis and Nm values of blue shark populations in various sampling regions

3 討論

大青鯊是分布最為廣泛、種群數(shù)量最多的一類大洋性鯊魚［32］，盡管大洋交界處會有不同地理種群的個體混合，但標(biāo)志重捕研究表明，大青鯊的運(yùn)動范圍局限在其本身所在的海域，各海域間大青鯊的基因交流很少［33］。一般認(rèn)為:群體之間的Nm＜1 時，說明群體可能由于遺傳漂變而發(fā)生了分化;而Nm＞1 時，則表明群體間的基因流水平較高，群體間遺傳分化較小;當(dāng)Nm＞4 時，表明種群間的基因交流更為充分，遺傳分化更?。?4］。本研究中得到的Nm值均比4 大得多，這說明大青鯊作為一種高度洄游的物種，它們之間的基因交流非常頻繁。根據(jù)ICCAT 的研究結(jié)果［35］，大西洋的大青鯊至少存在3 個不同的群體，分別為北大西洋、南大西洋和地中海群體。但Kohler等［21］、Shivji等［22］、Lessa等［23］根據(jù)線粒體DNA 控制區(qū)序列及微衛(wèi)星標(biāo)記分析得出，采樣區(qū)域內(nèi)北大西洋的東部和西部大青鯊種群不存在遺傳分化，這與本研究結(jié)果相符。對北大西洋和南大西洋大青鯊種群進(jìn)行的微衛(wèi)星分析，也證實了北大西洋和南大西洋種群之間存在的遺傳差異［36］。由于受到北美洲大陸的地理性隔離，北大西洋種群和北太平洋種群被認(rèn)為是兩種不同的群體。對北大西洋大青鯊的標(biāo)志重捕研究，顯示了大青鯊的季節(jié)性緯度運(yùn)動特點(diǎn)［37-38］，釋放點(diǎn)的大青鯊向東、西兩個方向上的遷移都很常見，但跨越赤道的運(yùn)動卻極少，結(jié)果為“北大西洋大青鯊為單一種群”這一假說提供了證據(jù)。此外，Karl等［39］采用線粒體控制區(qū)測定技術(shù)和微衛(wèi)星技術(shù)對大西洋西部低鰭真鯊種群的遺傳關(guān)系進(jìn)行了相關(guān)研究，表明相毗鄰的種群之間沒有顯著的遺傳差異，且遺傳多樣性水平較低;Hoelzel等［40］檢測了采集自大西洋西北部、大西洋東北部、地中海、印度洋和西太平洋姥鯊Cetorhinus maximus 的mtDNA控制區(qū)的多樣性，發(fā)現(xiàn)其遺傳變異度相當(dāng)?shù)?，不同海域間沒有太大的區(qū)別。本研究結(jié)果也顯示，大青鯊種群的遺傳差異不顯著，其原因可能是由于軟骨魚類的DNA 進(jìn)化速率比其他魚類慢造成的，同時大青鯊具有生長速度慢、性成熟晚、繁殖周期長、繁殖率低等生物特征，而且大青鯊在大西洋中又是金槍魚延繩釣的主要兼捕對象，所以捕撈壓力對于大青鯊種群的遺傳分化造成了較大的影響。

遺傳多樣性是生物多樣性的重要組成部分，是保護(hù)生物學(xué)研究的核心之一，對遺傳多樣性的研究具有重要的理論和實踐意義［41］。單倍型多樣性和核苷酸多樣性是衡量一個物種或種群線粒體DNA變異程度的兩個重要指標(biāo)［42］。Martin［43］在比較了鯊魚和哺乳動物線粒體序列后發(fā)現(xiàn)，鯊魚的線粒體進(jìn)化速率遠(yuǎn)低于哺乳動物和許多高等脊椎動物，即使是公認(rèn)進(jìn)化速率最快的控制區(qū)，其進(jìn)化速率也比其他脊椎動物低。Duncan等［44］對采自全球各大洋20 個海域的27l 尾路氏雙髻鯊Sphyrnalewini 線粒體控制區(qū)序列進(jìn)行分析時，僅發(fā)現(xiàn)22 個多態(tài)位點(diǎn)，24 個單倍型，表明路氏雙髻鯊線粒體DNA 變異程度較低。但本研究中對三大洋5 個海域的165 尾大青鯊線粒體D - loop 區(qū)序列進(jìn)行分析時，發(fā)現(xiàn)了110 個多態(tài)位點(diǎn)，145 個單倍型，分子方差分析得到Fst值為0.007 21。當(dāng)0＜Fst＜0.05 時，表明群體遺傳分化水平較弱，因此，作者認(rèn)為，大青鯊的線粒體DNA(D -loop 區(qū))具有較高的變異水平，只是由于不同群體之間強(qiáng)烈的基因交流，導(dǎo)致其遺傳分化水平較低。

迄今為止，對大青鯊種群遺傳結(jié)構(gòu)的研究多數(shù)是基于標(biāo)記重捕技術(shù)開展的。利用線粒體DNA 標(biāo)記等分子手段對大青鯊種群進(jìn)行的遺傳結(jié)構(gòu)研究很少，并且僅集中在大西洋。國內(nèi)外利用線粒體DNA 標(biāo)記對印度洋和太平洋大青鯊種群遺傳結(jié)構(gòu)的研究尚處于空白，這在很大程度上限制了對全球大青鯊種群遺傳結(jié)構(gòu)的全面認(rèn)識，從而影響到對大青鯊種群的有效保護(hù)和管理。本研究是首次利用線粒體D-loop 分子標(biāo)記對三大洋的大青鯊群體遺傳結(jié)構(gòu)進(jìn)行的研究，為大青鯊資源的管理和保護(hù)提供了一定的理論依據(jù)。但本研究中5 個海域的采樣樣本數(shù)量相對不足，需要在今后的研究中進(jìn)一步增加采樣點(diǎn)和樣本數(shù)，以期獲得更準(zhǔn)確的數(shù)據(jù)。

［1］Mejuto J，García-Cortés B.Reproductive and distribution parameters of the blue shark Prionace glauca，on the basis of on -board observations at sea in the Atlantic，Indian and Pacific oceans［J］.Collect Vol Sci Pap ICCAT，2005，58(3):974 -1000.

［2］Nakano H.Age，reproduction and migration of blue shark in the North Pacific Ocean［J］.Bulletin - National Research Institute of Far Seas Fisheries，1994，31:141 -190.

［3］Tavares R，Ortiz M，Arocha F.Population structure，distribution and relative abundance of the blue shark(Prionace glauca)in the Caribbean Sea and adjacent waters of the North Atlantic［J］.Fisheries Research，2012，129:137 -152.

［4］高春霞，田思泉，戴小杰，等.熱帶中東大西洋擬錐齒鯊生物學(xué)的初步研究［J］.上海海洋大學(xué)學(xué)報，2013，22(2):289 -294.

［5］郭新紅，劉少軍，劉巧，等.魚類線粒體DNA 研究新進(jìn)展［J］.遺傳學(xué)報，2004，31(9):983 -1000.

［6］閆杰，許強(qiáng)華，陳新軍，等.東太平洋公海莖柔魚種群遺傳結(jié)構(gòu)初步研究［J］.水產(chǎn)學(xué)報，2012，35(11):1617 -1623.

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

［8］陳驍，楊圣云，潘聰.我國南部海域條紋斑竹鯊線粒體DNA 控制區(qū)遺傳多態(tài)性研究［J］.海洋學(xué)報，2009，30(6):115 -121.

［9］Mendon?a F F，Ussami L H F，Hashimoto D T，et al.Identification and characterization of polymorphic microsatellite loci in the blue shark Prionace glauca，and cross - amplification in other shark species［J］.Journal of Fish Biology，2012，80(7):2643 -2646.

［10］Keeney D B，Heupel M R，Hueter R E，et al.Microsatellite and mitochondrial DNA analyses of the genetic structure of blacktip shark(Carcharhinus limbatus)nurseries in the northwestern Atlantic，Gulf of Mexico，and Caribbean Sea［J］.Molecular Ecology，2005，14(7):1911 -1923.

［11］吳峰，戴小杰，姜潤林.熱帶中東大西洋海域大青鯊繁殖生物學(xué)研究［J］.海洋湖沼通報，2012(3):30 -36.

［12］Lessa R P.Distribution and relative abundance of the blue shark，Prionace glauca，in the southwestern equatorial Atlantic Ocean［J］.Fishery Bulletin，1994，92:474 -480.

［13］Montealegre-Quijano S，Vooren C M.Distribution and abundance of the life stages of the blue shark(Prionace glauca)in the Southwest Atlantic［J］.Fisheries Research，2010，101(3):168-179.

［14］Vas P.The abundance of the blue shark，Prionace glauca，in the western English Channel［J］.Environmental Biology of Fishes，1990，29(3):209 -225.

［15］Zhu Jiangfeng，Dai Xiaojie，Xu Liuxiong，et al.Reproductive biology of female blue shark Prionace glauca in the southeastern Pacific Ocean［J］.Environmental Biology of Fishes，2011，91(1):95 -102.

［16］Sciarrotta T，Nelson D.Diel behavior of the blue shark，Prionace glauca，near Santa Catalina Island［J］.California Fish Bull，1977，75(3):519 -528.

［17］Queiroz N，Lima F P，Maia A，et al.Movement of blue shark，Prionace glauca，in the north-east Atlantic based on mark -recapture data［J］.Journal of the Marine Biological Association of the United Kingdom，2005，85(5):1107 -1112.

［18］Stevens，J D，Bradford R W，West G L.Satellite tagging of blue sharks(Prionace glauca)and other pelagic sharks off eastern Australia:depth behaviour，temperature experience and movements［J］.Marine Biology，2010，157(3):575 -591.

［19］宋利明，胡振新.馬紹爾群島海域大青鯊棲息地綜合指數(shù)［J］.水產(chǎn)學(xué)報，2011，35(8):1208 -1216.

［20］戴小杰，許柳雄.東太平洋熱帶海域大青鯊繁殖生物學(xué)特征［J］.水產(chǎn)學(xué)報，2005，29(4):565 -569.

［21］Kohler N E，Turner P A，Hoey J J，et al.Tag and recapture data for three pelagic shark species:blue shark(Prionace glauca)，shortfinmako(Isurus oxyrinchus)，and porbeagle(Lamna nasus)in the north Atlantic ocean［J］.Collect Vol Sci Pap ICCAT，2002，54(4):1231 -1260.

［22］Shivji M，Clarke S，Pank M，et al.Genetic identification of pelagic shark body parts for conservation and trade monitoring［J］.Conservation Biology，2002，16(4):1036 -1047.

［23］Lessa R，Santana F M，Hazin F H.Age and growth of the blue shark(Prionace glauca)(Linnaeus，1758)off northeastern Brazil［J］.Fisheries Research，2004，66(1):19 -30.

［24］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.

［25］Kumar S，Tamura K，Nei M.MEGA3:integrated software for molecular evolutionary genetics analysis and sequence alignment［J］.Briefings in Bioinformatics，2004，5(2):150 -163.

［26］Ma B，Xin L，Zhang K Z.A new quartet approach for reconstructing phylogenetic trees:quartet joining method［J］.Journal of Combinatorial Optimization，2008，16(3):293 -306.

［27］Zhang W，Sun Z R.Random local neighbor joining:a new method for reconstructing phylogenetic trees［J］.Molecular Phylogenetics and Evolution，2008，47(1):117 -128.

［28］Kimura M A.Simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences［J］.Journal of Molecular Evolution，1980，16(2):111 -120.

［29］Felsenstein J.Confidence limits on phylogenies:an approach using the bootstrap［J］.Evolution，1985:783 -791.

［30］Rozas J，Juan C S，Messeguer X，et al.DnaSP，DNA polymorphism analyses by the coalescent and other methods［J］.Bioinformatics，2003，19(18):2496 -2497.

［31］Excoffier L，Smouse P E，Quattro J M.Analysis of molecular variance inferred from metric distances among DNA haplotypes:application to human mitochondrial DNA restriction data［J］.Genetics，1992，131(2):479 -491.

［32］Kubodera T，Watanabe H，Ichii T.Feeding habits of the blue shark，Prionace glauca，and salmon shark，Lamna ditropis，in the transition region of the Western North Pacific［J］.Reviews in Fish Biology and Fisheries，2007，17(2/3):111 -124.

［33］Mandelman J W，Cooper P W，Werner T B，et al.Shark bycatch and depredation in the U.S.Atlantic pelagic longline fishery［J］.Reviews in Fish Biology and Fisheries，2008，18(4):427 -442.

［34］樂小亮，趙爽，劉海林，等.基于線粒體細(xì)胞色素b 基因全序列的海南長臀鮠南渡江種群遺傳變異分析［J］.生態(tài)科學(xué)，2010，29(3):247 -250.

［35］Campana S E，Joyce W，Manning M J.Bycatch and discard mortality in commercially caught blue sharks Prionace glauca assessed using archival satellite pop -up tags［J］.Marine Ecology Progress Series，2009，387:241 -253.

［36］Fitzpatrick S，Shivji M S，Chapman D D，et al.Development and characterization of 10 polymorphic microsatellite loci for the blue shark，Prionace glauca，and their cross shark-species amplification［J］.Conservation Genetics Resources，2011，3(3):523 -527.

［37］Stevens J.First results of shark tagging in the north -east Atlantic，1972 -1975［J］.Mar Biol Assoc UK，1976，56(4):929 -937.

［38］Kohler N E，Turner P A.Shark tagging:a review of conventional methods and studies［J］.Environmental Biology of Fishes，2001，60(1/3):191 -224.

［39］Karl S A，Castro A L F，Lopez J A，et al.Phylogeography and conservation of the bull shark(Carcharhinus leucas)inferred from mitochondrial and microsatellite DNA［J］.Conservation Genetics，2011，12(2):371 -382.

［40］Hoelzel A R，Shivji M S，Magnussen J，et al.Low worldwide genetic diversity in the basking shark(Cetorhinus maximus)［J］.Biology Letters，2006，2(4):639 -642.

［41］尚占環(huán)，姚愛興.生物遺傳多樣性研究方法及其保護(hù)措施［J］.寧夏農(nóng)學(xué)院學(xué)報，2002，23(1):66 -69.

［42］Nei M.Molecular evolutionary genetics［J］.Journal of Human Evolution，1989，18(8):808 -810.

［43］Martin A P，Naylor G J，Palumbi S R.Rates of mitochondrial DNA evolution in sharks are slow compared with mammals［J］.Nature，1992，357:153 -155.

［44］Duncan K M，Martin A P，Bowen B W，et al.Global phylogeography of the scalloped hammerhead shark(Sphyrna lewini)［J］.Molecular Ecology，2006，15(8):2239 -2251.