小清河河口鄰近海域口蝦蛄(Oratosquilla oratoria)COI序列特征及遺傳多樣性分析

劉 英, 王全超, 李海匯, 紀(jì)瑩璐, 陳琳琳, 李寶泉

小清河河口鄰近海域口蝦蛄()序列特征及遺傳多樣性分析

劉 英1, 2, 王全超1, 李海匯3, 紀(jì)瑩璐4, 陳琳琳1, 李寶泉1, 2

(1. 中國(guó)科學(xué)院煙臺(tái)海岸帶研究所, 山東 煙臺(tái) 264003; 2. 中國(guó)科學(xué)院大學(xué), 北京 100049; 3. 壽光市漁業(yè)發(fā)展中心, 山東 壽光 262700; 4. 國(guó)家海洋局北海預(yù)報(bào)中心, 山東 青島 266061)

受河流沿岸工農(nóng)業(yè)發(fā)展以及人口急劇增加影響, 小清河河口及鄰近海域受到嚴(yán)重污染?？谖r蛄()作為小清河河口及鄰近海域重要漁業(yè)資源之一, 承受著巨大的捕撈和生存壓力。為了掌握口蝦蛄種群在該區(qū)域的種群遺傳多樣性, 本研究于2020年6—10月在小清河河口及鄰近海域采集口蝦蛄生物樣本, 進(jìn)行了基于線粒體基因的種群遺傳多樣性分析。研究結(jié)果表明, 基于216條長(zhǎng)度為655 bp的基因片段, 共定義90個(gè)單倍型, 63個(gè)多態(tài)位點(diǎn), 其中包括59個(gè)轉(zhuǎn)換位點(diǎn), 1個(gè)顛換位點(diǎn)和3個(gè)轉(zhuǎn)換與顛換同時(shí)存在的位點(diǎn), 核苷酸多樣性指數(shù)和單倍型多樣性分別為0.970 0和0.005 1。與其他海域相比, 小清河河口種群遺傳多樣性處于相對(duì)較高水平; 同時(shí)遺傳結(jié)構(gòu)分析發(fā)現(xiàn)黃渤海與東海及南海分布的口蝦蛄種群遺傳變異較大, 但黃渤海海域內(nèi)遺傳變異較小。通過(guò)本研究基本掌握了小清河河口及鄰近海域口蝦蛄種群遺傳多樣性背景, 也可為口蝦蛄資源管理以及種質(zhì)資源庫(kù)建立等提供科學(xué)的基礎(chǔ)理論依據(jù)。

小清河河口及鄰近海域; 口蝦蛄; 線粒體; 遺傳多樣性

口蝦蛄(), 俗稱爬蝦、蝦爬子, 隸屬于節(jié)肢動(dòng)物門(Arthropoda)、甲殼綱(Crustacea)、口足目(Stomatopoda)、蝦蛄科(Squillidae), 自然分布于西太平洋溫?zé)釒\水海區(qū), 在中國(guó)沿海均廣泛分布, 其中以黃渤海產(chǎn)量最大[1]?？谖r蛄棲息于潮下帶泥沙質(zhì)海底, 其食性主要以貝類及甲殼類為主[2]?？谖r蛄具有個(gè)體大、產(chǎn)量高、繁殖周期短、易捕撈等優(yōu)點(diǎn), 使其具有很高的經(jīng)濟(jì)價(jià)值, 是中國(guó)沿海重要的海產(chǎn)經(jīng)濟(jì)種類, 同時(shí)其肉質(zhì)鮮美、營(yíng)養(yǎng)價(jià)值高, 沿海民眾喜食。因此, 口蝦蛄多年來(lái)均作為近海漁業(yè)主要捕撈對(duì)象之一[3]。特別在渤海, 口蝦蛄已成為支柱漁業(yè)之一[4]。但近年來(lái), 由于全球氣候變化、環(huán)境污染、養(yǎng)殖疾病危害以及過(guò)度捕撈等導(dǎo)致口蝦蛄野生資源急劇減少[5], 并進(jìn)而影響或降低口蝦蛄野生種群的遺傳多樣性。

鑒于口蝦蛄重要的經(jīng)濟(jì)價(jià)值, 國(guó)內(nèi)對(duì)該物種已開展了較多研究, 涉及到生物學(xué)特征[6-8]、生理生態(tài)[9]、漁業(yè)資源[10-11]和人工繁育[12-13]等方面。隨著分子生物學(xué)技術(shù)發(fā)展以及口蝦蛄資源管理和種質(zhì)資源保護(hù)的需要, 對(duì)口蝦蛄種群遺傳多樣性的研究也日漸增多。劉海映等[14]基于隨機(jī)擴(kuò)增多態(tài)性DNA標(biāo)記(random amplified polymorphic DNA, RAPD)分析了大連海域口蝦蛄遺傳多樣性; Liu等[15]報(bào)道了口蝦蛄線粒體全長(zhǎng)序列, 自此關(guān)于利用線粒體基因作為遺傳標(biāo)記研究其種群遺傳多樣性的報(bào)道也日漸增多; Zhang等[16]使用線粒體基因作為遺傳標(biāo)記分析了口蝦蛄在中國(guó)海域的地理分布、擴(kuò)散與遺傳差異; Du等[17]以作為遺傳標(biāo)記分析口蝦蛄在渤海、黃海和東海的遺傳多樣性; 隋宥珍等[18]基于細(xì)胞色素(cytochrome,)基因分析了口蝦蛄黃渤海種群、東海種群、南海種群的遺傳多樣性; 隋宥珍等[19]基于細(xì)胞色素氧化酶亞基I(cytochrome oxidase subunit I,)基因分析了東海海域口蝦蛄遺傳多樣性。上述研究基本摸清了中國(guó)不同海域口蝦蛄種群的遺傳多樣性現(xiàn)狀和不同地理種群的遺傳分化。但對(duì)個(gè)別區(qū)域, 尤其是具有重要生態(tài)意義的河口區(qū)如小清河河口, 口蝦蛄種群的遺傳多樣性研究迄今還尚未詳細(xì)開展。

小清河河口及鄰近海域曾經(jīng)是渤海萊州灣重要的漁業(yè)區(qū), 是渤海生態(tài)系統(tǒng)中漁業(yè)生物的產(chǎn)卵場(chǎng)與索餌場(chǎng)。小清河是典型受污染河流之一, 沿河污染物隨河流入海, 使河口及鄰近海域污染加重[20-22], 且導(dǎo)致河口內(nèi)出現(xiàn)低氧區(qū)[23]。小清河河口及鄰近海域的污染加重以及不斷增大的捕撈強(qiáng)度, 必然導(dǎo)致該區(qū)口蝦蛄種群數(shù)量的維持和生存受到影響。因此有必要對(duì)口蝦蛄開展遺傳多樣性研究, 摸清該海域口蝦蛄的遺傳分化潛力, 以維持資源生產(chǎn)力。線粒體DNA(mitochondrial DNA, mtDNA)為母系遺傳, 基因結(jié)構(gòu)簡(jiǎn)單, 且容易提取擴(kuò)增, 是研究海洋生物遺傳多樣性的重要手段?；蚴蔷€粒體基因的重要組成部分, 具有結(jié)構(gòu)簡(jiǎn)單、進(jìn)化速率適中、具有較高的多態(tài)性以及易被通用引物擴(kuò)增的優(yōu)勢(shì), 已被廣泛用于種群遺傳多樣性研究。本研究采集了小清河河口及鄰近海域口蝦蛄生物樣本, 基于口蝦蛄種群基因標(biāo)記研究其遺傳多樣性, 評(píng)估小清河河口及鄰近海域口蝦蛄種質(zhì)資源現(xiàn)狀, 為該種群健康維持、資源合理開發(fā)利用以及相關(guān)部門管理決策的制定提供理論性參考。

1 材料與方法

1.1 材料

2020年6—10月于小清河河口及其鄰近海域(表1), 采用漁業(yè)拖網(wǎng)方式, 采集口蝦蛄生物樣品。隨機(jī)選取部分樣品放于體積分?jǐn)?shù)95%的乙醇中暫存并轉(zhuǎn)移至實(shí)驗(yàn)室備用。

1.2 DNA提取、擴(kuò)增及測(cè)序

取口蝦蛄背部肌肉, 使用海洋動(dòng)物組織基因組DNA提取試劑盒[天根生化科技(北京)有限公司]按照操作說(shuō)明書提取每個(gè)口蝦蛄個(gè)體總DNA。用NanoDrop 1000分光光度計(jì)(NanoDrop, Wilmington, DE, USA)測(cè)定提取DNA的純度和完整性, 并進(jìn)行1%瓊脂糖凝膠電泳。合格DNA放于–20 ℃暫存?zhèn)溆谩?/p>

COI基因片段擴(kuò)增引物為通用引物, 引物序列: LCO1490: 5′-GGTCAACAAATCATAAAGATATTGG-3′;HCO2198: 5′-TAAACTTCAGGGTGACCAAAAAATCA-3′[24]。PCR擴(kuò)增在50 μL反應(yīng)體系中進(jìn)行。50 μL反應(yīng)體系包括: 25 μL含DNA聚合酶反應(yīng)混合物[天根生化科技(北京)有限公司], 17 μL ddH2O, 上游引物和下游引物各2 μL(100 μmol/L)及模板DNA 4 μL。PCR擴(kuò)增按照以下熱循環(huán)程序進(jìn)行: 94 ℃初步變性5 min; 隨后94 ℃變性45 s, 40 ℃退火45 s, 72 ℃延伸30 s, 進(jìn)行35個(gè)循環(huán); 最后在72 ℃下延伸10 min, 10 ℃保存。PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢驗(yàn)合格后, 送測(cè)序公司[生工生物工程(上海)股份有限公司]進(jìn)行雙向測(cè)序。

表1 口蝦蛄樣本采集信息

1.3 數(shù)據(jù)分析

通過(guò)BioEdit軟件對(duì)獲得的序列進(jìn)行可視化、修剪、編輯、拼接以及人工校對(duì), 獲得完整無(wú)誤的序列, 并在NCBI(https://www.ncbi.nlm.nih.gov/)上進(jìn)行復(fù)檢。在MEGA 7.0.26進(jìn)行CLUSTAL比對(duì)和序列堿基組成及變異分析。采用DAMEB 7軟件分析基因序列的堿基替換飽和度。使用DNAsp v6.12.03計(jì)算多態(tài)位點(diǎn)()、總變異位點(diǎn)()、單倍型數(shù)量()、單倍型多樣性指數(shù)(H)、遺傳多樣性指數(shù)()、平均核苷酸差異數(shù)()等遺傳多樣性參數(shù), 分析種群遺傳多樣性。

為明確小清河河口及鄰近海域口蝦蛄種群遺傳結(jié)構(gòu)現(xiàn)狀, 采用從GenBank上額外下載的唐山、青島、舟山和廣州口蝦蛄種群的序列, 進(jìn)行系統(tǒng)發(fā)育樹構(gòu)建和遺傳結(jié)構(gòu)分析。序列包括唐山種群(GenBank登錄號(hào): MF173560-MF173566)、青島種群(GenBank登錄號(hào): MF173526-MF173534)、舟山種群(GenBank登錄號(hào): MF173576-MF173598)[25]以及廣州種群(GenBank登錄號(hào): HQ621811-HQ621829)[26]。分子系統(tǒng)發(fā)育樹的重建是在MEGA 7.0.26中使用基于Kimura雙參數(shù)法模型(Kimura 2-parameter)遺傳距離的非加權(quán)分組平均法(unweighted pair group method using arithmetic average, UPGMA)進(jìn)行, 并選用bootstrap法進(jìn)行1 000次重復(fù)自舉檢驗(yàn), 評(píng)估系統(tǒng)發(fā)育樹。遺傳結(jié)構(gòu)分析是利用Arlequin 3.5.2.2軟件開展分子方差分析(AMOVA)和遺傳分化系數(shù)ST計(jì)算。

2 結(jié)果

2.1 小清河河口及鄰近海域口蝦蛄COI序列

獲得序列后經(jīng)剪切、人工校正及NCBI復(fù)檢后, 共獲得216條長(zhǎng)度為655 bp的口蝦蛄基因部分序列片段。序列片段無(wú)插入和缺失。所有序列A、G、C、T的平均含量分別為36.6%、17.1%、18.4%和27.9%, A+T含量明顯高于G+C含量。4種堿基在密碼子第1位上分布較為均勻(A=26.48%, G=32.96%, C=16.90%, T=24.65%), 第2位上T的含量較高(A=13.30%, G=18.81%, C=25.23%, T=42.66% ), 第3位上呈現(xiàn)G和C的含量較低的偏倚現(xiàn)象(A=43.83%, G=4.32%, C=9.22%, T=42.63%)。

所有序列共發(fā)現(xiàn)66個(gè)變異位點(diǎn), 其中63個(gè)多態(tài)位點(diǎn), 由30個(gè)單一突變位點(diǎn)和33個(gè)簡(jiǎn)約信息位點(diǎn)組成。63個(gè)多態(tài)位點(diǎn)中包括59個(gè)轉(zhuǎn)換位點(diǎn), 1個(gè)顛換位點(diǎn)和3個(gè)轉(zhuǎn)換與顛換同時(shí)存在的位點(diǎn)。4個(gè)變異位點(diǎn)發(fā)生在密碼子的第1位, 59個(gè)變異位點(diǎn)發(fā)生在第3位, 無(wú)變異位點(diǎn)發(fā)生在第2位。655 bp的基因序列共定義218個(gè)氨基酸, 且堿基的替換并未導(dǎo)致氨基酸序列發(fā)生變化。

堿基的替換飽和分析見圖1, 成對(duì)堿基的轉(zhuǎn)換率、顛換率都與遺傳距離成線性相關(guān), 其中, 94.18%的堿基轉(zhuǎn)換和5.84%的堿基顛換能分別得到相關(guān)線性分析。序列堿基替換飽和度分析結(jié)果也表明, 標(biāo)準(zhǔn)堿基替換飽和值(ss.c)顯著大于實(shí)際檢測(cè)到的基因堿基替換值(ss)(<0.05)。

2.2 小清河河口口蝦蛄種群遺傳多樣性以及與其他海域的比較

該海域口蝦蛄種群H為0.970 0,為0.005 1,為3.330 0。所有序列共定義90個(gè)單倍型, 其中, 單倍型Hap-1、Hap-16、和Hap-27分別為20、19和15個(gè)個(gè)體所有, 是優(yōu)勢(shì)單倍型。

圖1 小清河河口口蝦蛄種群COI基因堿基替換飽和分析

本研究口蝦蛄種群?jiǎn)伪缎投鄻有耘c粵東海種群、東海種群相當(dāng), 略高于渤海灣種群、皮口種群、綏中種群和青島種群; 與Zhu等[25]研究的唐山種群、龍口種群、青島種群、連云港種群以及七星列島種群?jiǎn)伪缎投鄻有韵喈?dāng), 但略高于東營(yíng)種群、小洋山種群、舟山種群和三門種群; 核苷酸多樣性低于粵東海種群, 東海種群, 高于渤海灣種群、遼東灣種群和青島種群, 低于Zhu等[25]研究中的青島種群、連云港種群和七星列島種群, 高于唐山種群、小洋山種群和三門種群(表2)。

2.3 口蝦蛄種群遺傳結(jié)構(gòu)

小清河河口口蝦蛄種群與唐山種群、青島種群、舟山種群以及廣州種群之間的親緣關(guān)系和遺傳結(jié)構(gòu)見圖2。分子系統(tǒng)發(fā)育樹顯示所有單倍型聚為兩支: 小清河河口種群首先與青島種群聚集在一起, 然后再與唐山種群匯合, 組成一支, 廣州種群和舟山種群組成則組成另一支。

AMOVA分析結(jié)果顯示, 5個(gè)口蝦蛄種群遺傳變異主要來(lái)自于種群間, 且5個(gè)口蝦蛄種群間具有顯著的遺傳分化水平(ST=0.857 38,<0.05)(表3)。不同種群間成對(duì)ST范圍在–0.020 62～0.900 91之間(表4)。結(jié)果表明, 唐山種群、小清河河口種群及青島種群與舟山種群和廣州種群間存在顯著遺傳分化(<0.05), 但唐山種群、小清河河口種群和青島種群間不存在顯著遺傳分化(>0.05), 廣州種群和舟山種群間不存在顯著遺傳分化(>0.05)。

表2 其他海域口蝦蛄種群遺傳多樣性

圖2 基于5個(gè)口蝦蛄種群COI單倍型序列遺傳距離構(gòu)建的UPGMA系統(tǒng)發(fā)育樹

3 討論

遺傳多樣性是生物多樣性的基礎(chǔ), 開展遺傳多樣性的研究, 對(duì)于漁業(yè)資源管理的可持續(xù)利用以及漁業(yè)資源的管理和保護(hù)都具有重要意義[29]。一般而言, 遺傳多樣性水平高的物種, 具有較強(qiáng)的環(huán)境適應(yīng)能力、生存能力和進(jìn)化潛力。因此, 物種遺傳多樣性一旦降低, 可導(dǎo)致物種適應(yīng)及生存能力降低、物種退化甚至滅絕[30]。

物種在基因水平上的遺傳變異最直接的表現(xiàn)就是堿基組成的差異。目前已觀測(cè)到的無(wú)脊椎動(dòng)物線粒體基因組中堿基組成基本規(guī)律均為A+T含量顯著高于G+C含量[31]。本研究中口蝦蛄種群的序列堿基組成A+T含量高于G+C含量, 該結(jié)果與三疣梭子蟹()[32]、脊尾白蝦()[33]、芋螺[34]的基因片段研究結(jié)果相似。本研究中, 口蝦蛄的序列堿基轉(zhuǎn)換與顛換之比為9, 遠(yuǎn)大于基因變異達(dá)到飽和狀態(tài)的臨界值2[35], 說(shuō)明小清河河口的序列變異還未達(dá)到飽和狀態(tài)。堿基轉(zhuǎn)換數(shù)大于顛換數(shù), 符合在動(dòng)物線粒體進(jìn)化過(guò)程中堿基的替換主要以轉(zhuǎn)換為主的結(jié)論[36]。并且一般在線粒體DNA進(jìn)化過(guò)程中, 親緣關(guān)系較近物種間, 轉(zhuǎn)換更容易頻繁發(fā)生, 而顛換則在親緣關(guān)系較遠(yuǎn)的物種間更容易發(fā)生。同一物種內(nèi)部, 轉(zhuǎn)換數(shù)量高于顛換數(shù)量[37]。同時(shí),基因片段的堿基飽和性分析顯示小清河河口口蝦蛄種群實(shí)際堿基替換值顯著低于期望堿基替換飽和值, 表明基因堿基替換未達(dá)到飽和狀態(tài)?；诨蛭葱Ｕ霓D(zhuǎn)換、顛換與遺傳距離的散點(diǎn)分布圖得到的線性關(guān)系也表明小清河河口口蝦蛄種群的基因序列未達(dá)到飽和狀態(tài), 因此, 線粒體基因適合用作研究口蝦蛄群體遺傳多樣性的分子標(biāo)記。

表3 基于COI序列的口蝦蛄5個(gè)地理種群間的AMOVA分析

注:<0.05差異顯著

表4 基于COI序列的口蝦蛄5個(gè)地理種群間的FST值

注: *表示種群間差異顯著,<0.05差異顯著

在所有堿基變異中, 有59個(gè)變異位點(diǎn)發(fā)生在密碼子的第3位上, 僅有4個(gè)變異位點(diǎn)發(fā)生在密碼子的第1位上, 無(wú)變異位點(diǎn)發(fā)生在密碼子的第2位上。可見, 密碼子第3位的變異率明顯高于第1位和第2位, 這與密碼子第1位相對(duì)保守、第2位非常保守以及第3位具有搖擺性的特點(diǎn)相符, 這與某些魚類研究結(jié)果一致[38-39]。同時(shí)本結(jié)果檢測(cè)到的所有核苷酸變異均為同義突變。這進(jìn)一步說(shuō)明密碼子簡(jiǎn)并性大多發(fā)生在密碼子第3位上且在物種穩(wěn)定上起著重要作用。

Grant和Bowen[40]根據(jù)線粒體基因序列的單倍型多樣性和核苷酸多樣性, 將海洋生物的遺傳多樣性分為4種類型: 第1類, 低單倍型多樣性(H>0.5), 低核苷酸多樣性(<0.005)型, 近期種群瓶頸或是由單一、少數(shù)種群所產(chǎn)生的建立者效應(yīng); 第2類, 高單倍型多樣性(H>0.5), 低核苷酸多樣性(<0.005)型, 種群瓶頸, 隨后是種群快速增長(zhǎng)和突變積累; 第3類,低單倍型多樣性(H<0.5), 高核苷酸多樣性(>0.005)型, 種群差異, 這種差異可能來(lái)自孤立種群間的二次接觸或者是一個(gè)大而穩(wěn)定的種群中發(fā)生過(guò)瓶頸效應(yīng); 第4類, 高單倍型多樣性(H>0.5), 高核苷酸多樣性(>0.005)型, 種群由一個(gè)大而穩(wěn)定的種群經(jīng)過(guò)長(zhǎng)時(shí)間演化所產(chǎn)生, 或是兩個(gè)不同系群的種群二次接觸產(chǎn)生。本研究口蝦蛄種群總的單倍型多樣性H為0.970 0, 核苷酸多樣性為0.005 1, 屬于第4類種群狀態(tài), 即種群屬于高單倍型多樣性和高核苷酸多樣性模式, 種群遺傳多樣性較高, 處于長(zhǎng)期穩(wěn)定狀態(tài)。本研究口蝦蛄種群?jiǎn)伪缎投鄻有院秃塑账岫鄻有栽谥袊?guó)沿海多數(shù)海域中處于相對(duì)較高水平, 這說(shuō)明小清河河口口蝦蛄種群的遺傳多樣性在我國(guó)沿海處于相對(duì)較高的水平, 可考慮作為種質(zhì)資源選育及建立種質(zhì)資源庫(kù)的選擇之一。

種群遺傳結(jié)構(gòu)受地理隔離、洋流、生境、人類活動(dòng)、遷移和歷史因素等多種因素的影響[41-42]。唐山種群、小清河河口種群和青島種群位于黃渤海, 這3個(gè)種群間兩兩相互聚為姐妹群, AMOVA和ST分析也顯示3個(gè)種群間沒(méi)有顯著遺傳分化, 表明3個(gè)種群間極可能存在基因交流, 而且小清河河口種群與青島種群交流更為頻繁。一方面在黃渤海水交換過(guò)程中, 黃海沿岸流流動(dòng)方向基本終年不變, 沿山東半島北岸從渤海進(jìn)入北黃海, 然后南下到達(dá)南黃海[43], 另一方面口蝦蛄具有大約2個(gè)月的幼蟲浮游期。因此, 口蝦蛄幼體可在海洋環(huán)流的驅(qū)動(dòng)下發(fā)生擴(kuò)散, 從而維持口蝦蛄種群連通性[44]。此外, 本研究中AOMVA、ST和分子系統(tǒng)發(fā)育樹也均顯示以長(zhǎng)江為界, 長(zhǎng)江口以北口蝦蛄種群和長(zhǎng)江口以南的種群間存在顯著的遺傳分化。但目前, 關(guān)于中國(guó)沿?？谖r蛄種群遺傳結(jié)構(gòu)尚未得到統(tǒng)一。Cheng和Sha[45]通過(guò)線粒體基因標(biāo)記和核基因標(biāo)記分析認(rèn)為, 在長(zhǎng)江入?？趦蓚?cè)存在兩個(gè)口蝦蛄亞種。一些學(xué)者分別通過(guò)對(duì)和分析, 認(rèn)為東海和黃海的口蝦蛄種群的遺傳分化的原因可能是由長(zhǎng)江淡水的注入形成的物理屏障和局部適應(yīng)造成[17, 25]。也有學(xué)者通過(guò)以線粒體基因和核基因?yàn)檫z傳標(biāo)記, 認(rèn)為口蝦蛄在冰川時(shí)期受臺(tái)灣海峽分隔, 長(zhǎng)期隔離形成了南海種群孤立于其他海域種群的一個(gè)顯著的南北種群結(jié)構(gòu)[16]。因此未來(lái)對(duì)中國(guó)沿?？谖r蛄種群的遺傳研究中, 在選擇遺傳標(biāo)記時(shí), 應(yīng)綜合考慮遺傳標(biāo)記的進(jìn)化速度、多態(tài)性、遺傳信息的全面性等。

4 結(jié)論與展望

口蝦蛄作為小清河河口及鄰近海域常年優(yōu)勢(shì)種之一, 面臨著巨大的捕撈和生存壓力。研究結(jié)果表明小清河河口及鄰近海域口蝦蛄種群遺傳多樣性目前還處于較高水平, 為以后口蝦蛄的苗種選育、野生種質(zhì)資源保護(hù)和合理開發(fā)利用提供了新的依據(jù)。為了維持該海域口蝦蛄種群的遺傳多樣性水平, 建議加強(qiáng)口蝦蛄資源保護(hù), 基于現(xiàn)有資源量的科學(xué)評(píng)估, 進(jìn)行合理采捕, 如控制網(wǎng)具網(wǎng)目大小, 減小網(wǎng)目尺寸, 以利于口蝦蛄幼體的存活和自然種群的延續(xù), 保障口蝦蛄野生資源可持續(xù)利用。同時(shí), 由于口蝦蛄在4月下旬便進(jìn)入繁殖期, 而目前休漁期在5月1日至9月1日, 有必要考慮適當(dāng)延長(zhǎng)休漁時(shí)間至4月下旬。

[1] 吳強(qiáng), 陳瑞盛, 黃經(jīng)獻(xiàn), 等. 萊州灣口蝦蛄的生物學(xué)特征與時(shí)空分布[J]. 水產(chǎn)學(xué)報(bào), 2015, 39(8): 1166- 1177.

WU Qiang, CHEN Ruisheng, HUANG Jingxian, et al. Fishery biology characteristics, temporal and spatial distribution ofin Laizhou Bay, Bohai Sea[J]. Journal of Fisheries of China, 2015, 39(8): 1166-1177.

[2] 寧加佳, 杜飛雁, 王雪輝, 等. 基于穩(wěn)定同位素的口蝦蛄食性分析[J]. 水產(chǎn)學(xué)報(bào), 2016, 40(6): 903-910.

NING Jiajia, DU Feiyan, WANG Xuehui, et al. Feeding habits of mantis shrimp based on stable isotope analysis[J]. Journal of Fisheries of China, 2016, 40(6): 903- 910.

[3] 劉修澤, 郭棟, 王愛勇, 等. 遼東灣海域口蝦蛄的資源特征及變化[J]. 水生生物學(xué)報(bào), 2014, 38(3): 602- 608.

LIU Xiuze, GUO Dong, WANG Aiyong, et al. The re-source characteristics and their variation ofin Liaodong Bay[J]. Acta Hydrobiolo-gica Sincia, 2014, 38(3): 602-608.

[4] 李顯森, 吳亞飛, 尤宗博, 等. 渤?？谖r蛄三重刺網(wǎng)漁獲組成及其捕撈性能分析[J]. 海洋漁業(yè), 2016, 38(5): 516-524.

LI Xiansen, WU Yafei, YOU Zongbo, et al. Analysis on the catch composition by trammel net and its fishing performance forin the Bohai Sea[J]. Marine Fisheries, 2016, 38(5): 516-524.

[5] 谷德賢, 劉茂利. 天津海域口蝦蛄群體結(jié)構(gòu)及資源量分析[J]. 河北漁業(yè), 2011, 212(8): 24-26.

GU Dexian, LIU Maoli. Analysis on the Population structure and abundance of in Tianjin Sea Area[J]. Hebei Fisheries, 2011, 212(8): 24-26.

[6] 王春琳, 徐善良, 梅文驤, 等. 口蝦蛄的生物學(xué)基本特征[J]. 浙江水產(chǎn)學(xué)院學(xué)報(bào), 1996, 15(1): 60-62.

WANG Chunling, XU Shanliang, MEI Wenxiang, et al. A biological basic character of[J]. Journal of Zhenjian College of Fisheries, 1996, 15(1): 60-62.

[7] 林月嬌, 劉海映, 徐海龍, 等. 大連近海口蝦蛄形態(tài)參數(shù)關(guān)系的研究[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào), 2008, 23(3): 215-217.

LIN Yuejiao, LIU Haiying, XU Hailong, et al. Morpho-metry of mantis shrimpin Dalian coast[J]. Journal of Dalian Ocean University, 2008, 23(3): 215-217.

[8] 孫東昱, 王蕾, 邱盛堯, 等. 海陽(yáng)和煙臺(tái)近岸口蝦蛄群體的形態(tài)比較[J]. 漁業(yè)科學(xué)進(jìn)展, 2021, 42(1): 154-164.

SUN Dongyu, WANG Lei, QIU Shengyao, et al. Morphological comparison ofstocks along the coast of Haiyang and Yantai[J].Progress in Fishery Sciences, 2021, 42(1): 154-164.

[9] 王春琳, 母昌考, 丁愛俠, 等. 口蝦蛄()同工酶的組織特異性及生化遺傳分析[J]. 海洋與湖沼, 2004, 35(3): 258-263.

WANG Chunlin, MU Changkao, DING Aixia, et al. Tissue specificity and biochemical genetic analysis of isozyme on crustacean[J]. Ocea-nologia et Limnologia Sinica, 2004, 35(3): 258-263.

[10] 潘國(guó)良, 張龍, 朱增軍, 等. 浙江南部近岸海域春季口蝦蛄生物量的時(shí)空分布[J]. 海洋與湖沼, 2013, 44(2): 366-370.

PAN Guoliang, ZHANG Long, ZHU Zengjun, et al. Spatial-temporal distribution of the biomass ofin the coastal waters of south Zhejiang during spring[J]. Oceanologia et Limnologia Sinica, 2013, 44(2): 366-370.

[11] 李鵬程, 張崇良, 任一平, 等. 山東近?？谖r蛄空間分布特征及其與環(huán)境因子的關(guān)系[J]. 中國(guó)水產(chǎn)科學(xué), 2020, 27(12): 1515-1523.

LI Pengcheng, ZHANG Chongliang, REN Yiping, et al. Relationship betweenspatial distribution and environmental factors in coastal Shangdong[J]. Journal of Fishery Sciences of China, 2020, 27(12): 1515-1523.

[12] 王春琳, 葉選怡, 丁愛俠, 等. 蝦蛄繁殖生物學(xué)與繁育技術(shù)研究[J]. 海洋湖沼通報(bào), 2002(3): 58-64.

WANG Chunlin, YE Xuanyi, DING Aixia, et al. Studies on reproductive biology and breeding technolo-gy of squilla[J]. Transactions of Oceanology and Lim-nology, 2002(3): 58-64.

[13] 張年國(guó), 潘桂平, 周文玉, 等. 不同餌料和水溫對(duì)口蝦蛄成活、生長(zhǎng)及育肥性能的影響[J]. 動(dòng)物學(xué)雜志, 2019, 54(3): 425-435.

ZHANG Nianguo, PAN Guiping, ZHOU Wenyu, et al. Effects of different diets and water temperatures on survival, growth, fattening performance of[J]. Chinese Journal of Zoology, 2019, 54(3): 425-435.

[14] 劉海映, 王桂娥, 王秀利. 大連海域口蝦蛄資源遺傳多樣性的分析[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào), 2009, 24(4): 350-353.

LIU Haiying, WANG Guie, WANG Xiuli. Genetic di-versity analysis of mantis shrimpfrom Dalian coast[J]. Journal of Dalian Ocean Uni-versity, 2009, 24(4): 350-353.

[15] LIU Y, CUI Z X. The complete mitochondrial genome of the mantid shrimp(Crustacea: Malacostraca: Stomatopoda): Novel non-coding regions features and phylogenetic implications of the Stomatopoda[J]. Comparative Biochemistry and Physiology Part D: Genomics and Proteomics, 2010, 5(3): 190-198.

[16] ZHANG D Z, DING G, GE B M, et al. Geographical distribution, dispersal and genetic divergence of the mantis shrimp(Stomatopoda: Squillidae) in China Sea[J]. Biochemical Systematics and Ecology, 2016, 65: 1-8.

[17] DU X W, CAI S S, YU C G, et al. Population genetic structure of mantis shrimps Oratosquilla oratoria: Testing the barrier effect of the Yangtze River outflow[J]. Biochemical Systematics and Ecology, 2016, 66: 12-18.

[18] 隋宥珍, 劉連為, 徐開達(dá), 等. 基于線粒體基因的口蝦蛄種群遺傳結(jié)構(gòu)研究[J]. 大連海洋大學(xué)學(xué)報(bào), 2019, 34(3): 355-361.

SUI Youzhen, LIU Lianwei, XU Kaida, et al. Population genetic structure of mantis shrimpbased on the partial mitochondrial DNA cytochromegene[J]. Journal of Dalian Ocean University, 2019, 34(3): 355-361.

[19] 隋宥珍, 周永東, 盧占暉, 等. 東海海域口蝦蛄種群遺傳多樣性[J]. 動(dòng)物學(xué)雜志, 2016, 51(2): 291-300.

SUI Youzhen, ZHOU Yongdong, LU Zhanhui, et al. Genetic diversity of mantis shrimp () populations in the East China Sea[J]. Chinese Journal of Zoology, 2019, 51(2): 291-300.

[20] 劉建強(qiáng). 萊州灣海洋工程建設(shè)對(duì)小清河口環(huán)境影響數(shù)值研究[D]. 青島: 中國(guó)海洋大學(xué), 2012.

LIU Jianqiang. Numerical study on the environment effects of marine engineering on the mouth of Xiaoqing River in Laizhou Bay[D]. Qingdao: Ocean University of China, 2012.

[21] LUO X X, ZHANG S S, YANG J Q, et al. Macrobenthic community in the Xiaoqing River Estuary in Laizhou Bay, China[J]. Journal of Ocean University of China, 2013, 12(3): 366-372.

[22] LUO X X, LI S, YANG J Q, et al. Benthic habitat quality assessment based on biological indices in Xiaoqing River estuary and its adjacent sea of Laizhou Bay, China[J]. Journal of Ocean University of China, 2017, 16(3): 537-546.

[23] 孟春霞, 鄧春梅, 姚鵬, 等. 小清河口及鄰近海域的溶解氧[J]. 海洋環(huán)境科學(xué), 2005, 24(3): 25-28.

MENG Chunxia, DNEG Chunmei, YAO Peng, et al. Dissolved oxygen in the Xiaoqinghe Estuary and Adja-cent Waters[J]. Marine Environmental Science, 2005, 24(3): 25-28.

[24] FOLMER O, BLACK M, HOEH W, et al. DNA primers for amplification of mitochondrial cytochrome coxidase subunit I from diverse metazoan invertebrates[J]. Molecular Marine Biology and Biotechnology, 1994, 3(5): 294-299.

[25] ZHU W B, DU X W, HAN Z Q, et al. The genetic divergence ofbetween the East China Sea and Yellow Sea: physical barrier and possible local adaptation[J]. Journal of the Marine Biological Association of the United Kingdom, 2019, 99(3): 631-638.

[26] 黃映萍, 王瑩, 苗素英. 粵東海域口蝦蛄遺傳多樣性[J]. 動(dòng)物學(xué)雜志, 2011, 46(2): 82-89.

HUANG Yingping, WANG Ying, MIAO Suying. Genetic diversity analysis of mantis shrimp () from the eastern coast of Guangzhou provin-ce[J]. Chinese Journal of Zoology, 2011, 46(2): 82-89.

[27] 董鑫, 邢坤, 隋宥珍, 等. 基于線粒體基因序列的4個(gè)海域口蝦蛄群體的遺傳多樣性研究[J]. 海洋科學(xué), 2015, 39(7): 29-36.

DONG Xin, XING Kun, SUI Youzhen, et al. Genetic diversity offrom four dea waters based on the mitochondrialgene sequences analysis[J]. Marine Sciences, 2015, 39(7): 29-36.

[28] 張代臻, 丁鴿, 張華彬, 等. 渤海灣口蝦蛄線粒體基因的遺傳多樣性研究[J]. 南京師大學(xué)報(bào) (自然科學(xué)版), 2010, 33(4): 80-83.

ZHANG Daizhen, DING Ge, ZHANG Huabin, et al. Genetic diversity offrom different sea area of Bohai Bay by mitchondrialgene[J]. Journal of Nanjing Normal University (Natural Science Edition), 2010, 33(4): 80-83.

[29] RYMAN N, UTTER F, LAIKRE L. Protection of intraspecific biodiversity of exploited fishes[J]. Reviews in Fish Biology and Fisheries, 1995, 5(4): 417-446.

[30] 季維智, 宿兵. 遺傳多樣性研究的原理和方法[M]. 杭州: 浙江科學(xué)技術(shù)出版社, 1999.

JI Weizhi, SU Bing.Principles and methodologies of genetic diversity studies[M]. Hangzhou: Zhejiang Science and Technology Publishing House, 1999.

[31] 楊建敏, 鄭小東, 王如才, 等. 3種鮑基因片段序列的初步研究[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2003, 33(1): 36-40.

YANG Jianmin, ZHENG Xiaodong, WANG Rucai, et al. Comparative study on mtDNAgene sequence of three species of abalone[J]. Periodical of Ocean University of China, 2003, 33(1): 36-40.

[32] 王景, 張鳳英, 蔣科技, 等. 基于線粒體基因序列的三疣梭子蟹東海區(qū)群體遺傳多樣性分析[J]. 海洋漁業(yè), 2015, 37(2): 114-121.

WANG Jing, ZHANG Fengying, JIANG Keji, et al. Genetic diversity ofbased on the mitochondrial Cytochrome Oxidase Subunit I sequence from the East China Sea[J], Marine Fisheries, 2015, 37(2): 114-121.

[33] 馬朋, 劉萍, 李健, 等. 脊尾白蝦3個(gè)野生群體線粒體基因的遺傳多樣性及其系統(tǒng)發(fā)育分析[J]. 漁業(yè)科學(xué)進(jìn)展, 2011, 32(6): 50-56.

MA Peng, LIU Ping, LI Jian, et al. The genetic diversity and phylogenetic analysis ofgene in mitochondrial DNA of three populations of[J]. Progress in Fishery Sciences, 2011, 32(6): 50-26.

[34] 陳琴, 彭超, 趙峰, 等. 基于序列研究南海芋螺遺傳多樣性[J]. 基因組學(xué)與應(yīng)用生物學(xué), 2020, 39(9): 3955-3960.

CHEN Qin, PENG Chao, ZHAO Feng, et al. Study on genetic diversity of Cone snail in South China Sea based onsequences[J]. Genomics and Applied Biology, 2020, 39(9): 3955-3960.

[35] KNIGHT A, MINDELL D P. Substitution bias, weigh-ting of DNA sequence evolution, and the phylogenetic position of Fea's Viper[J]. Systematic Biology, 1993, 42(1): 18-31.

[36] IRWIN D M, KOEHER T D, WILSON A C. Evolution of the cytochromegene of mammals[J]. Jounal of Molecular Evolution, 1991, 32(2): 128-144.

[37] KOCHER T D, THOMAS W K, MEYER A, et al. Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved prime-rs[J]. Proceedings of the National Academy of Sciences, 1989, 86(16): 6196-6200.

[38] 李鵬飛, 賀舟挺, 徐開達(dá), 等. 基于基因和序列的東海鰳魚種質(zhì)資源遺傳變異研究[J]. 浙江海洋學(xué)院學(xué)報(bào) (自然科學(xué)版), 2013, 32(2): 93-98.

LI Pengfei, HE Zhouting, XU Kaida, et al. Genetic diversity ofin the East China Sea revealed byandsequence[J]. Journal of Zhejiang Ocean University (Natural Science), 2013, 32(2): 93-98.

[39] 李鵬飛, 周永東, 徐漢祥. 大黃魚、鮸魚及美國(guó)紅魚線粒體DNA的基因序列比較[J]. 南方水產(chǎn), 2008, 4(3): 43-47.

LI Pengfei, ZHOU Yongdong, XU Hanxiang. Sequence analysis on mitochondrial (mtDNA)genes inand[J]. South China Fisheries Science, 2008, 4(3): 43-47.

[40] GRANT W S, BOWEN B W. Shallow population his-tories in deep evolutionary lineages of marine fishes: Insights from sardines and anchovies and lessons for conservation[J]. Journal of Heredity, 1998, 89(5): 415- 426.

[41] GROSBERG R K, CUNNINGHAM C W. Genetic stru-cture in the sea from population to community[M]. Sunderland: Marine Community Ecology, 2001.

[42] LEHMANN T, HAWLEY W A, GREBERT H, et al. The effective population size ofin Kenya: implications for population structure[J]. Molecular Biology and Evolution, 1998, 15(3): 264-276.

[43] GUAN B X. Patterns and structures of the currents in Bohai, Huanghai and East China Seas[M]. Dordrecht: Oceanology of China Seas, 1994.

[44] YANG M, LI X Z. Population genetic structure of the mantis shrimp(Crustacea: Squillidae) in the Yellow Sea and East China Sea[J]. Journal of Oceanology and Limnology, 2018, 36(3): 905-912.

[45] CHENG J, SHA Z L. Cryptic diversity in the Japanese mantis shrimp(Crustacea: Squil-lidae): Allopatric diversification, secondary contact and hybridization[J]. Scientific Reports, 2017, 7: 1-13.

Analysis of thesequence characteristics and genetic diversity ofin the Xiaoqing River Estuary and its adjacent sea

LIU Ying1, 2, WANG Quan-chao1, LI Hai-hui3, JI Ying-lu4, CHEN Lin-lin3, LI Bao-quan1, 2

(1. Yantai Institute of Coastal Zone Research, Chinese Academy of Sciences, Laboratory of Protection of Biological Resources Biology in Coastal Waters, Yantai 264003, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. Shouguang Marine Fishery Development Center, Shouguang 262700, China; 4. North China Sea Marine Forecasting Center of State Oceanic Administration, Qingdao 266061, China)

The Xiaoqing River Estuary has been severely polluted due to the quick development of industry and agriculture as well as the rapid increase in the population in the region. Mantis shrimp () is an important fishery resource in this area that is under great threat from overfishing and pollution. In this study, samples were collected in the Xiaoqing River Estuary and its adjacent sea during June–August 2020 to understand the population genetic diversity based ongene. A total of 90 haplotypes and 63 polymorphic loci, including 59 transition loci, 1 transversion loci, and 3 loci where both a transition and transversion existed, were defined based on 216gene fragments with a length of 655 bp. The nucleotide diversity and haplotype diversity were 0.970 0 and 0.005 1, respectively.Compared with other sea areas, we found that the genetic diversity of the Xiaoqing River Estuary population was relatively high. Meanwhile, genetic structure analysis showed an obvious difference in the genetic variation of the mantis shrimp population on a large spatial scale, e.g., the Yellow Sea and the Bohai Sea, the East China Sea and the South China Sea. Whereas, it did not show any difference on a relatively smaller spatial scale, e.g., the Yellow Sea and the Bohai Sea. The genetic diversity background of the mantis shrimp population in the Xiaoqing River Estuary and its adjacent sea could provide a theoretical basis for resource management and germplasm bank establishment.

XiaoqingRiver Estuary and its adjacent sea areas;;gene; genetic diversity

Jun. 18, 2021

S917.4

A

1000-3096(2022)11-0038-09

10.11759/hykx20210618001

2021-06-18;

2021-07-06

美麗中國(guó)生態(tài)文明建設(shè)科技工程專項(xiàng)(XDA23050304); 山東省自然基金面上項(xiàng)目(ZR2021MD084)

[The Strategic Priority Research Program of the Chinese Acade-my of Sciences, No. XDA23050304; the Shandong Provincial Natural Scie-n-ce Foundation, No. ZR2021MD084]

劉英(1996—), 女, 四川內(nèi)江人, 碩士, 主要從事大型底棲動(dòng)物生理生態(tài)研究, E-mail: 3104931657@qq.com; 李寶泉(1972—),通信作者, 男, 研究員, 主要從事海洋生物多樣性、底棲動(dòng)物生理生態(tài)和生態(tài)健康評(píng)價(jià)研究, E-mail: bqli@yic.ac.cn

(本文編輯: 楊 悅)