劉 潔 高風(fēng)英 陳 剛 曹建萌 劉志剛 王 淼 盧邁新

珠江流域9個(gè)大鱗副泥鰍群體遺傳多樣性的微衛(wèi)星分析*

劉 潔1,2高風(fēng)英1①陳 剛2曹建萌1劉志剛1王 淼1盧邁新1①

(1. 中國(guó)水產(chǎn)科學(xué)研究院珠江水產(chǎn)研究所 農(nóng)業(yè)農(nóng)村部熱帶亞熱帶水產(chǎn)資源利用與養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室 廣州 510380；2. 廣東海洋大學(xué)水產(chǎn)學(xué)院 湛江 524025)

為探索我國(guó)珠江流域大鱗副泥鰍()群體的遺傳分化及親緣關(guān)系，本研究采用8個(gè)微衛(wèi)星分子標(biāo)記，對(duì)我國(guó)珠江流域9個(gè)大鱗副泥鰍群體(佛山、高要、封開(kāi)、肇慶、乳源、樂(lè)昌、韶關(guān)、河源和惠州)進(jìn)行了遺傳多樣性分析。結(jié)果顯示，8個(gè)微衛(wèi)星位點(diǎn)共檢測(cè)到 69個(gè)等位基因，平均等位基因(a)和有效等位基因(e)分別為8.6和4.0個(gè)，平均觀測(cè)雜合度(o)和期望雜合度(e)分別為0.4426和0.7030。9個(gè)大鱗副泥鰍群體間的遺傳分化系數(shù)(st)和基因流(m)分別為0.2452和0.7697，表明群體間遺傳分化水平較高，遺傳交流水平低。采用UPGMA法，對(duì) 9個(gè)群體基于遺傳距離進(jìn)行聚類，可分為兩大支，韶關(guān)、佛山和乳源群體聚為一支；另一支包含了其余的6個(gè)大鱗副泥鰍群體，分為3個(gè)小分支，分別為樂(lè)昌與肇慶群體、河源與惠州群體及高要與封開(kāi)群體。研究表明，珠江流域的9個(gè)大鱗副泥鰍群體具有較高的遺傳多樣性，且群體間存在著一定的遺傳分化，具有進(jìn)一步選育的價(jià)值。

大鱗副泥鰍；微衛(wèi)星；遺傳多樣性

大鱗副泥鰍()主要分布于我國(guó)東部地區(qū)的河流、水田及湖泊水域，其外部形態(tài)與泥鰍()極為相似(杜啟艷等, 2007)。大鱗副泥鰍肉質(zhì)細(xì)嫩、味道鮮美、營(yíng)養(yǎng)豐富，素有“水中人參”之稱(趙睿等, 2014)。近年來(lái)，我國(guó)向日本、韓國(guó)出口大量大鱗副泥鰍，隨著市場(chǎng)需求的增加，大鱗副泥鰍的養(yǎng)殖規(guī)模也不斷擴(kuò)大，而良種選育工作的滯后使苗種的人工繁育量已不能滿足養(yǎng)殖規(guī)模不斷擴(kuò)大的需要(李彩娟等, 2014)。此外，環(huán)境污染和近親繁殖導(dǎo)致野生泥鰍的資源量銳減，種質(zhì)資源不斷退化(邱楚雯等, 2017)。親魚(yú)捕撈后的異地養(yǎng)殖和人工養(yǎng)殖個(gè)體的逃逸干擾了泥鰍的自然分布和群體結(jié)構(gòu)。因此，研究大鱗副泥鰍的遺傳多樣性，對(duì)合理開(kāi)發(fā)和保護(hù)我國(guó)現(xiàn)有大鱗副泥鰍資源具有重要意義。

遺傳多樣性為生物多樣性的重要組成部分。一個(gè)群體或物種的遺傳多樣性越高或遺傳變異越豐富，越能適應(yīng)多變的環(huán)境，其進(jìn)化潛力也越大(沈浩等, 2001)。對(duì)物種遺傳多樣性的數(shù)據(jù)分析，可為物種資源的開(kāi)發(fā)和生物多樣性的保護(hù)提供科學(xué)依據(jù)(羅靜等, 2000)。遺傳多樣性的研究方法隨著遺傳學(xué)和分子生物學(xué)的發(fā)展，已到DNA分子水平。其中，微衛(wèi)星又稱簡(jiǎn)單重復(fù)序列(Simple sequence repeats, SSR)，是指以少數(shù)幾個(gè)核苷酸(2~6)經(jīng)多次重復(fù)組成的簡(jiǎn)單串聯(lián)重復(fù)序列(劉萍等, 2008)。SSR具有高多態(tài)性、共顯性、穩(wěn)定性強(qiáng)等優(yōu)點(diǎn)，被廣泛應(yīng)用于水產(chǎn)動(dòng)物遺傳多樣性評(píng)估、遺傳結(jié)構(gòu)和親緣關(guān)系分析等研究中(王軍等，2018; 李可等, 2018; 葉竹青等, 2015; 蘇天鳳等, 2010)。目前，有研究用到微衛(wèi)星技術(shù)對(duì)我國(guó)各地區(qū)大鱗副泥鰍進(jìn)行遺傳多樣性分析。李彩娟等(2014)采用17個(gè)多態(tài)性微衛(wèi)星標(biāo)記，發(fā)現(xiàn)四大湖泊大鱗副泥鰍自然群體具有豐富的遺傳多樣性。游翠紅等(2012)篩選出6對(duì)微衛(wèi)星引物，對(duì)我國(guó)長(zhǎng)江水系及廣州地區(qū)的7個(gè)野生大鱗副泥鰍群體遺傳結(jié)構(gòu)分析顯示，大鱗副泥鰍種群遺傳多樣性較低。珠江流域大鱗副泥鰍群體遺傳多樣性分析還未見(jiàn)報(bào)道。

本研究采用8個(gè)微衛(wèi)星，對(duì)我國(guó)珠江流域共9個(gè)群體大鱗副泥鰍的遺傳多樣性與遺傳分化進(jìn)行研究，充分探索大鱗副泥鰍的遺傳背景，以期為其種質(zhì)資源保護(hù)策略的制定以及良種選育基礎(chǔ)群體的選擇提供理論依據(jù)。

1 材料與方法

1.1 實(shí)驗(yàn)材料

2017年8~10月于珠江水系采集泥鰍活體，采集后于中國(guó)水產(chǎn)科學(xué)研究院珠江水產(chǎn)研究所實(shí)驗(yàn)室暫養(yǎng)，取尾鰭并保存于無(wú)水乙醇中。9個(gè)大鱗副泥鰍群體野生樣品分別采自珠江流域的西江(佛山，24尾；高要，30尾；封開(kāi)，35尾；肇慶，27尾)、北江(乳源，25尾；樂(lè)昌，30尾；韶關(guān)，34尾)和東江(河源，29尾；惠州，20尾) (圖1)，共計(jì)254尾個(gè)體。

鑒別泥鰍和大鱗副泥鰍的特異性引物：上游引物(5′CCGCCCAGGAGTATTTTATG3′)；下游引物(5′TGGCATTTCATCAGGGGTG3′) (侯吉倫等, 2015)。8個(gè)微衛(wèi)星(Mac37、Mac45、Mac49、Mac60、Mac425、Mac456、Mac576和Mac627)摘自泥鰍微衛(wèi)星連鎖圖譜中等位基因較多的位點(diǎn)(Kagayaki, 2008)，擴(kuò)增片段長(zhǎng)度為100~300 bp，便于電泳分型。上述引物均由生工生物工程(上海)股份有限公司合成。

1.2 實(shí)驗(yàn)方法

1.2.1 基因組DNA的提取 使用美基DNA提取試劑盒提取泥鰍尾鰭基因組DNA。用1%瓊脂糖凝膠120 V電泳20 min檢測(cè)所提DNA質(zhì)量；用微量分光光度計(jì)檢測(cè)DNA濃度和質(zhì)量；加ddH2O稀釋DNA至50 ng/μl，–20℃保存?zhèn)溆谩?/p>

圖1 大鱗副泥鰍采集地點(diǎn)()

FS：佛山；GY：高要；FK：封開(kāi)；ZQ：肇慶；RY：乳源；LC：樂(lè)昌；SG：韶關(guān)；HY：河源；HZ：惠州

FS: Foshan; GY: Gaoyao; FK: Fengkai; ZQ: Zhaoqing; RY: Ruyuan; LC: Lechang; SG: Shaoguan; HY: Heyuan; HZ: Huizhou

1.2.2 泥鰍和大鱗副泥鰍鑒別 用鑒別引物進(jìn)行PCR擴(kuò)增，PCR反應(yīng)體系(25 μl)：50 ng/μl DNA模板2.5 μl，10 μmol/L上下游引物各2.5 μl，10×LA PCR Buffer(含Mg2+）2.5 μl，2.5 mmol/L dNTP Mixture 2.0 μl，5 U/μl LADNA聚合酶0.25 μl，ddH2O 12.75 μl。PCR反應(yīng)程序：94℃預(yù)變性5 min；94℃變性30 s，54℃退火90 s，72℃延伸1 min，35個(gè)循環(huán)；72℃延伸10 min。用1%瓊脂糖凝膠120 V電泳20 min，檢測(cè)PCR擴(kuò)增產(chǎn)物，電泳圖譜顯示，被測(cè)樣品條帶大小鑒別泥鰍(285 bp)和大鱗副泥鰍(214 bp)。

1.2.3 PCR擴(kuò)增 采用10對(duì)微衛(wèi)星引物進(jìn)行PCR擴(kuò)增，同時(shí)，篩選最佳的退火溫度。PCR反應(yīng)體系(20 μl)：50 ng/μl DNA模板1.0 μl，10 μmol/L上下游引物各0.2 μl 2×EasyTaq?PCR SuperMix for PAGE 10.0 μl，ddH2O 8.6 μl。PCR反應(yīng)程序：94℃預(yù)變性 3 min；94℃變性30 s，退火30 s，72℃延伸30 s， 23個(gè)循環(huán)；72℃延伸5 min。各對(duì)引物的特異退火溫度見(jiàn)表1。

表1 本研究使用的微衛(wèi)星位點(diǎn)及其相關(guān)信息

Tab.1 Microsatellites used in this study and their associated information

*：(X)代表核心序列重復(fù)次數(shù)

*: (X)means repeat number of core sequence

1.2.4 聚丙烯凝膠電泳及銀染 擴(kuò)增產(chǎn)物經(jīng)12%非變性聚丙烯酰胺凝膠電泳分離，上樣量為5 μl，內(nèi)外電泳槽中各加適量的1×TBE電泳緩沖液，50 V電泳30 min后再150 V電泳150 min。銀染參考徐紹斌等(2002)的方法并適當(dāng)調(diào)整，同一位點(diǎn)等位基因由大到小記為A、B、C…V。

1.2.5 數(shù)據(jù)分析 使用Popgene32 1.32軟件分析各微衛(wèi)星基因座的等位基因數(shù)(Number of allele,a)、有效等位基因數(shù)(Effective mumber of allele,e)、觀測(cè)雜合度(Observed heterozygosity,o)、期望雜合度(Expected heterozygosity,e)、群體內(nèi)近交系數(shù)(Fixation index,is)、遺傳分化系數(shù)(Genetic differentiation coefficient,st)、基因流(Number of migrants per generation,m)、遺傳距離(Genetic distance,A)、遺傳相似系數(shù)(Genetic similarity index,) (Kagayaki, 2008)，使用PIC-CALC 0.6軟件計(jì)算各位點(diǎn)的多態(tài)信息含量(Polymorphism information content, PIC)。

基于群體間Nei氏標(biāo)準(zhǔn)遺傳距離，使用MEGA6.0構(gòu)建7個(gè)群體UPGMA(Unweighted pair group method with arithmetic means)系統(tǒng)進(jìn)化樹(shù)。

2 結(jié)果

2.1 泥鰍和大鱗副泥鰍鑒別

鑒別引物在所有樣品中均能擴(kuò)增出清晰的條帶，且能根據(jù)PCR擴(kuò)增的產(chǎn)物片段大小來(lái)區(qū)分泥鰍(285 bp)和大鱗副泥鰍(214 bp)。部分大鱗副泥鰍鑒定的電泳結(jié)果見(jiàn)圖2。

圖2 大鱗副泥鰍鑒定的部分電泳結(jié)果

2.2 微衛(wèi)星引物的PCR擴(kuò)增

選取的10對(duì)泥鰍微衛(wèi)星引物在9個(gè)大鱗副泥鰍群體中，有8對(duì)能擴(kuò)增出清晰的條帶，且重復(fù)性好，在12%非變性聚丙烯酰胺凝膠上能分離等位基因。帶型統(tǒng)計(jì)分析后發(fā)現(xiàn)，8個(gè)微衛(wèi)星位點(diǎn)共獲得了69個(gè)等位基因，微衛(wèi)星基因座的等位基因范圍為5~15個(gè)，平均每個(gè)位點(diǎn)檢測(cè)出8.6個(gè)等位基因，其中，Mac49存在等位基因最少(5個(gè))，Mac456最多(15個(gè))。部分位點(diǎn)的電泳結(jié)果見(jiàn)圖3。

2.3 9個(gè)大鱗副泥鰍群體遺傳多樣性

從表2可以看出，8個(gè)微衛(wèi)星位點(diǎn)的平均觀測(cè)等位基因?yàn)?.6個(gè)，平均有效等位基因數(shù)為4.0個(gè)，平均觀測(cè)雜合度為0.4426，平均期望雜合度為0.7030。8個(gè)位點(diǎn)在9個(gè)大鱗副泥鰍群體中表現(xiàn)出不同水平的多態(tài)。參照Botstein等(1980)提出的衡量基因變異程度高低的PIC指標(biāo)，除位點(diǎn)Mac45和Mac49為中度多態(tài)(0.250.5)。按照Balloux等(2002)的標(biāo)準(zhǔn)，0~0.05、0.05~0.15和0.15~0.25分別是低、中和高遺傳分化，本實(shí)驗(yàn)中8個(gè)微衛(wèi)星位點(diǎn)的st平均值為0.2452，表明大鱗副泥鰍群體間遺傳分化水平較高。世代遷移數(shù)量m可反映基因流水平(喬利英等, 2010)。本研究m平均值為0.7697 (m<1)，表明9個(gè)群體間遺傳交流水平低。

圖3 位點(diǎn)Mac425和Mac576的部分電泳結(jié)果

表2 8個(gè)微衛(wèi)星位點(diǎn)總遺傳參數(shù)

Tab.2 Total genetic parameters of 8 microsatellite loci

從表3可見(jiàn)，9個(gè)大鱗副泥鰍群體的等位基因數(shù)差別較小，其中，高要群體的等位基因數(shù)(a)最少(2.25)，韶關(guān)群體最多(5.38)；有效等位基因數(shù)(e)高要群體最少(1.98)，韶關(guān)群體最多(3.15)；各群體平均觀測(cè)雜合度(o)為0.3307~0.5626；期望雜合度(e)為0.3499~0.6144；PIC為0.4629~0.6777，其中，樂(lè)昌和河源群體表現(xiàn)為中度多態(tài)，其余7個(gè)群體表現(xiàn)為高度多態(tài)。

2.4 9個(gè)大鱗副泥鰍群體間的遺傳相似度和遺傳距離

遺傳距離和遺傳相似性指數(shù)是衡量群體親緣關(guān)系遠(yuǎn)近的重要指標(biāo)(朱滔等, 2013)。根據(jù)Nei指數(shù)法計(jì)算9個(gè)泥鰍群體間遺傳相似性指數(shù)和遺傳距離(表4)。9個(gè)大鱗副泥鰍群體間遺傳相似度為0.1965~0.9467，遺傳距離為0.0548~1.6271。肇慶與樂(lè)昌群體的遺傳相似度最高(0.9467)，遺傳距離最近(0.0548)；高要與乳源群體的遺傳相似度最低(0.9467)，遺傳距離最遠(yuǎn)(1.6271)。

根據(jù)遺傳距離運(yùn)用UPGMA法對(duì)7個(gè)泥鰍群體進(jìn)行了聚類分析。系統(tǒng)樹(shù)結(jié)果表明(圖4)，9個(gè)大鱗副泥鰍群體分為兩大支。其中，韶關(guān)和佛山群體聚在一起，再與乳源群體聚為一支；另一支包含了其余的6個(gè)大鱗副泥鰍群體，其中，樂(lè)昌與肇慶群體聚為一支，河源與惠州群體聚為一支，這2個(gè)分支聚為一個(gè)分支后，再與高要和封開(kāi)群體的這個(gè)分支聚為一大分支。

表3 8個(gè)微衛(wèi)星位點(diǎn)在9個(gè)大鱗副泥鰍群體中的遺傳多樣性參數(shù)

Tab.3 Genetic diversity indices of 8 microsatellite loci in 9 P. dabryanus populations

表4 9個(gè)大鱗副泥鰍群體的遺傳相似性(對(duì)角線之上)與遺傳距離(對(duì)角線之下)

Tab.4 Genetic identity (above diagonal) and genetic distance (below diagonal) among 9 P. dabryanus populations

圖4 基于Nei′s無(wú)偏遺傳距離構(gòu)建的9個(gè)大鱗副泥鰍群體UPGMA樹(shù)

3 討論

3.1 微衛(wèi)星引物多樣性

本研究8個(gè)微衛(wèi)星位點(diǎn)多態(tài)性豐富，共獲得 69個(gè)等位基因，平均等位基因?yàn)?.6個(gè)，與游翠紅等(2012)利用6對(duì)泥鰍微衛(wèi)星引物跨種擴(kuò)增大鱗副泥鰍，每個(gè)位點(diǎn)檢測(cè)出的平均等位基因數(shù)相近。PIC可反映基因座位變異程度的高低，同時(shí)，也關(guān)系到該位點(diǎn)的可用性和使用效率，PIC越大，提供的遺傳信息越高，越具有可用性(邢晶晶, 2012)。本研究的8個(gè)微衛(wèi)星位點(diǎn)中，除Mac45和Mac49兩個(gè)位點(diǎn)為中度多態(tài)(0.250.5)。這些高度多態(tài)性的基因座位對(duì)于大鱗副泥鰍群體間遺傳信息的進(jìn)一步研究具有重要意義。

3.2 種內(nèi)遺傳多樣性

9個(gè)大鱗副泥鰍群體平均PIC為0.5602，其中，樂(lè)昌和河源群體表現(xiàn)為中度多態(tài)，其余7個(gè)群體PIC均高于0.5，表現(xiàn)為高度多態(tài)，具有進(jìn)一步選育的潛力。李雅娟等(2010)采用5個(gè)微衛(wèi)星分子標(biāo)記對(duì)我國(guó)3種野生泥鰍進(jìn)行分析顯示，平均觀測(cè)雜合度和期望雜合度分別為0.8696和0.7283，平均多態(tài)信息含量為0.6193，均高于本研究結(jié)果的平均觀測(cè)雜合度(0.4477)、期望雜合度(0.5409)和PIC(0.5602)，在一定程度反映出泥鰍的遺傳多樣性高于大鱗副泥鰍。Dewoody等(2000)采用微衛(wèi)星標(biāo)記對(duì)13種淡水魚(yú)進(jìn)行分析，得出平均期望雜合度為0.46，與本研究結(jié)果中9個(gè)大鱗副泥鰍的平均期望雜合度相近。李彩娟等(2014)采用 17個(gè)微衛(wèi)星分子標(biāo)記，分析四大湖泊大鱗副泥鰍的平均期望雜合度為0.577，高于本研究的結(jié)果，可能研究對(duì)象中有一部分個(gè)體為逃逸到野生環(huán)境中的養(yǎng)殖個(gè)體的原因。

3.3 遺傳分化與親緣關(guān)系

遺傳分化指數(shù)(st)是評(píng)價(jià)群體間遺傳分化程度的重要指標(biāo)，本研究中9個(gè)大鱗副泥鰍群體的遺傳分化指數(shù)均值為0.2452，表明群體間遺傳分化水平較高(劉志剛等, 2018)。丁海燕等(2015)利用6個(gè)微衛(wèi)星標(biāo)記分析我國(guó)泥鰍，遺傳分化指數(shù)均值為0.4990，高于本研究的遺傳分化指數(shù)值，多種原因可造成遺傳分化的嚴(yán)重，人類活動(dòng)造成的遺傳漂變可能是最主要的原因。Wright(1931)認(rèn)為，基因流指數(shù)m>1，基因流則能發(fā)揮均質(zhì)化作用，即能有效抑制種群間遺傳分化的產(chǎn)生；m<1時(shí)，則表明各亞群間基因流受阻，群體可能走向遺傳分化。本研究中8個(gè)座位的m均值為0.7697(m<1)，表明9個(gè)群體間基因交流水平低，進(jìn)而導(dǎo)致遺傳分化程度較高。這可能也是由于大鱗副泥鰍自身遷移能力有限、地理距離、自然屏障的隔離和自然選擇等導(dǎo)致了較大程度的遺傳分化(王坤, 2009)。

聚類分析結(jié)果顯示，9個(gè)大鱗副泥鰍群體大多數(shù)按地理位置聚類。乳源和韶關(guān)群體位于珠江水系的北江及其支流，與其他地區(qū)群體因地理距離遠(yuǎn)而缺乏基因交流，使得與其他地區(qū)間遺傳距離較大，在聚類圖上這2個(gè)大鱗副泥鰍群體聚為一大支與其他6個(gè)群體分開(kāi)。另外的一大支中，位于珠江東江流域的河源和惠州群體聚在一起，位于珠江西江流域的高要和封開(kāi)群體聚在一起。在游翠紅等(2012)對(duì)我國(guó)7個(gè)大鱗副泥鰍群體的研究中發(fā)現(xiàn)，多個(gè)地區(qū)群體也表現(xiàn)出相似的遺傳分化現(xiàn)象?？梢?jiàn)地理隔離是影響大鱗副泥鰍群體間遺傳分化的主要原因。但在9個(gè)群體中，樂(lè)昌群體屬于珠江水系的北江，卻與位于珠江西江流域的肇慶群體遺傳相似度最高(0.9467)，遺傳距離最近(0.0548)。此外，位于珠江水系西江的佛山群體與位于珠江北江的韶關(guān)群體的遺傳距離(0.2201)也較近。單磊等(2009)對(duì)長(zhǎng)江中下游泥鰍的遺傳多樣性的研究也發(fā)現(xiàn)，武漢和太湖群體之間地理位置較遠(yuǎn)而遺傳關(guān)系較近，表現(xiàn)出遺傳距離與地理距離不一致的現(xiàn)象。丁海燕等(2015)對(duì)我國(guó)7個(gè)泥鰍群體的遺傳多樣性研究也表明，四大流域東部盆地各泥鰍群體遺傳距離卻與地理距離沒(méi)有明確的關(guān)聯(lián)性。這說(shuō)明人類活動(dòng)或養(yǎng)殖行為可能在一定程度上擾動(dòng)了泥鰍資源的分布，也可能與實(shí)驗(yàn)樣品量不足有關(guān)，后續(xù)應(yīng)擴(kuò)大采樣范圍和增加實(shí)驗(yàn)群體。本研究中，9個(gè)珠江水系的大鱗副泥鰍群體具有較高的遺傳多樣性，且群體間存在一定的遺傳分化，可對(duì)其進(jìn)行種質(zhì)保護(hù)，防止與其他苗種混雜，造成種質(zhì)退化。

Balloux F, Lugon MN. The estimation of population differentiation with microsatellite markers. Molecular Ecology, 2002, 11(2): 155–165

Botstein D, White RL, Skolnick M,. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. American Journal of Human Genetics, 1980, 32(3): 314–331

DeWoody JA, Avise JC. Microsatellite variation in marine, freshwater and anadromous fishes compared with other animals. Journal of Fish Biology, 2000, 56(3): 461–473

Ding HY, Zhu M, Peng C,. Preliminary analysis of the microsatellite diversity of Chinese loach). Journal of Fisheries of China, 2015, 39(10): 1433–1442 [丁海燕, 朱明, 彭沖, 等. 我國(guó)泥鰍微衛(wèi)星遺傳多樣性初步分析. 水產(chǎn)學(xué)報(bào), 2015, 39(10): 1433–1442]

Du QY, Nan P, Yan SG,. Observation on the embryonic development of. Journal of Henan Normal University (Natural Science), 2007(2): 146–149 [杜啟艷, 南平, 燕帥國(guó), 等. 豫北地區(qū)大鱗副泥鰍胚胎發(fā)育的觀察. 河南師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2007(2): 146–149]

Hou JL, Liu YF, Liu HJ,. Identification primer and identification method forandPatent: CN104498604A, 2015- 04-08 [侯吉倫, 劉永富, 劉海金, 等. 一種泥鰍和大鱗副泥鰍種質(zhì)鑒別引物和鑒別方法. 專利: CN104498604A, 2015-04-08]

Kagayaki M, Ichiro N, Katsutoshi A. Genetic linkage map of the loach(Teleostei: Cobitidae). Genetica, 2008, 132(3): 227–241

Li CJ, Ling QF, Ge C,. Genetic diversity offrom four big lakes in China. Jiangsu Journal of Agricultural Sciences, 2014, 30(5): 1087–1094 [李彩娟, 凌去非, 葛辰, 等. 中國(guó)4大湖泊大鱗副泥鰍群體遺傳多樣性分析. 江蘇農(nóng)業(yè)學(xué)報(bào), 2014, 30(5): 1087–1094]

Li K, Ma XF, Xie P,. Genetic diversity and population genetic structure ofbased on microsatellite genetic markers. Journal of Huazhong Agricultural University, 2018, 37(4): 112–120 [李可, 馬徐發(fā), 謝鵬, 等. 基于微衛(wèi)星標(biāo)記的擬鯉遺傳多樣性及群體遺傳結(jié)構(gòu)分析. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào), 2018, 37(4): 112–120]

Li YJ, Yu Z, Zhang MZ,. Genetic diversity in Chinese loach species analyzed by microsatellite DNA markers. Journal of Northeast Agricultural University, 2010, 41(12): 75–79 [李雅娟, 于卓, 張明昭, 等. 我國(guó)三種野生泥鰍遺傳多樣性的微衛(wèi)星分析. 東北農(nóng)業(yè)大學(xué)學(xué)報(bào), 2010, 41(12): 75–79]

Liu P, Song LP, Li J,. Development of microsatellite DNA markers and their applications in genetic studies in crabs. Periodical of Ocean University of China, 2008, 38(5): 712–718 [劉萍, 宋來(lái)鵬, 李健, 等. 蟹類微衛(wèi)星DNA標(biāo)記的篩選及其在遺傳學(xué)研究中的應(yīng)用. 中國(guó)海洋大學(xué)學(xué)報(bào), 2008, 38(5): 712–718]

Liu ZG, Lu MX, Cao JM,. Genetic diversity and genetic relationship analysis of tilapia “Yuemin No.1” and its breeding populations. Progress in Fishery Sciences, 2018, 39(6): 31–41 [劉志剛, 盧邁新, 曹建萌, 等. 羅非魚(yú)“粵閩1號(hào)”及其繁育群體的遺傳多樣性和遺傳關(guān)系分析. 漁業(yè)科學(xué)進(jìn)展, 2018, 39(6): 31–41]

Luo J, Yang JX, Zhang YP,. Genetic basis of fish diversity. Zoological Research, 2000, 21(2): 158–164 [羅靜, 楊君興, 張亞平. 魚(yú)類多樣性的遺傳基礎(chǔ). 動(dòng)物學(xué)研究, 2000, 21(2): 158–164]

Qiao LY, Yuan YN. The measurement indices and influencing factors of microsatellite polymorphisms. Chinese Animal Husbandry and Veterinary Medicine, 2010, 37(1): 107–111 [喬利英, 袁亞男. 微衛(wèi)星標(biāo)記遺傳多樣性的度量指標(biāo)及影響因素. 中國(guó)畜牧獸醫(yī), 2010, 37(1): 107–111]

Qiu CW, Wang HX. Advances in research on genetic diversity of. Fisheries Science and Technology Information, 2017, 44(2): 78–86 [邱楚雯, 王韓信. 泥鰍屬魚(yú)類遺傳多樣性的研究進(jìn)展. 水產(chǎn)科技情報(bào), 2017, 44(2): 78–86]

Shan L, Wei KJ, Zhang GR,. Genetic diversity of four diploid loach) populations in the middle and lower Yangtze River Basin. Journal of Huazhong Agricultural University, 2009, 28(4): 453–458 [單磊, 魏開(kāi)建, 張桂蓉, 等. 長(zhǎng)江中下游二倍體泥鰍4個(gè)種群的遺傳多樣性. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào), 2009, 28(4): 453–458]

Shen H, Liu DY. Summary of genetic diversity. Journal of Biology, 2001(3): 5–7 [沈浩, 劉登義. 遺傳多樣性概述. 生物學(xué)雜志, 2001(3): 5–7]

Su TF, Xiong XF, Jiang SG,. Microsatellite analysis of genetic relationship among 7 full-sib families of. South China Fisheries Science, 2010, 6(6): 1–7 [蘇天鳳, 熊小飛, 江世貴, 等. 斑節(jié)對(duì)蝦7個(gè)全同胞家系間親緣關(guān)系的微衛(wèi)星分析. 南方水產(chǎn), 2010, 6(6): 1–7]

Wang J, Wang QY, Kong J,. SSR Analysis on genetic diversity in breeding and wild populations of. Progress in Fishery Sciences, 2018, 39(2): 104–111[王軍, 王清印, 孔杰, 等. 中國(guó)明對(duì)蝦人工選育群體與野生群體遺傳多樣性的SSR分析. 漁業(yè)科學(xué)進(jìn)展, 2018, 39(2): 104–111]

Wang K. A primary study on moleculargenetic differences and growth of different loach populations in the Yangtze River Basin. Master′s Thesis of Soochow University, 2009 [王坤. 長(zhǎng)江流域不同泥鰍群體的分子遺傳差異及生長(zhǎng)的初步研究. 蘇州大學(xué)碩士研究生學(xué)位論文, 2009]

Wright. Evolution in mendelian population. Genetics, 1931, 16: 97–159

Xing JJ. The types of molecular markers and the research and applications of molecular technology on the aquatic creature. Chinese Journal of Fisheries, 2002(1): 61–70 [邢晶晶. 分子遺傳標(biāo)記及其技術(shù)在水產(chǎn)生物中的研究與應(yīng)用. 水產(chǎn)學(xué)雜志, 2002(1): 61–70]

Xu SB, Tao YF, Yang ZQ. A simple and rapid methods used for silver staining and gel preservation. Hereditas, 2002, 24(3): 335–336 [許紹斌, 陶玉芬, 楊昭慶. 簡(jiǎn)單快速的DNA銀染和膠保存方法. 遺傳, 2002, 24(3): 335–336]

Ye ZQ, Ling QF, Li CJ,. Genetic structure analysis of five families of. Journal of Hydroecology, 2015, 36(1): 66–73 [葉竹青, 凌去非, 李彩娟, 等. 5個(gè)大鱗副泥鰍家系的遺傳結(jié)構(gòu)分析. 水生態(tài)學(xué)雜志, 2015, 36(1): 66–73]

You CH, Tong JG, Yu XM. Microsatellite DNA analysis on genetic diversity of seven populations of. Journal of Hydroecology, 2012, 33(1): 84–91 [游翠紅, 童金茍, 俞小牧. 大鱗副泥鰍7個(gè)群體遺傳變異的微衛(wèi)星分析. 水生態(tài)學(xué)雜志, 2012, 33(1): 84–91]

Zhao R, Shen F, Na LH. Artificial breeding technology of. China Fisheries, 2014(9): 59–60 [趙睿, 沈芾, 那立海. 大鱗副泥鰍人工繁育技術(shù). 中國(guó)水產(chǎn), 2014(9): 59–60]

Zhu T, Liang XF, Peng MH,. Isolation and characterization of polymorphic EST-SSRs markers in threepopulation. Journal of Jinan University, 2013, 34(3): 347–352 [朱滔, 梁旭方, 彭敏燕, 等. 翹嘴鱖EST-SSR標(biāo)記的開(kāi)發(fā)及3個(gè)群體遺傳多態(tài)性分析. 暨南大學(xué)學(xué)報(bào), 2013, 34(3): 347–352]

Genetic Diversity Analysis of Nine Loach () Populations in the Pearl River Basin Based on Microsatellite Markers

LIU Jie1,2, GAO Fengying1①, CHEN Gang2, CAO Jianmeng1, LIU Zhigang1, WANG Miao1, LU Maixin1①

(1. Pearl River Fisheries Research Institute, Chinese Academy of Fishery Science, Key Laboratory of Tropical and Subtropical Fishery Resource Application and Cultivation, Ministry of Agriculture and Rural Affairs, Guangzhou 510380;2. College of Fisheries, Guangdong Ocean University, Zhanjiang 524025)

Human activities and environmental pollution have caused varying degrees of damage toresources in China. To explore the genetic differentiation and phylogenetic relationship of.populations in China’s Pearl River Basin, eight microsatellite markers were used to analyze the genetic diversity of nine geographicalpopulations (Foshan, Gaoyao, Fengkai, Zhaoqing, Ruyuan, Lechang, Shaoguan, Heyuan, and Huizhou) of.in the Pearl River Basin. The analysis detected 69 alleles in eight microsatellite loci. The average number of alleles (a) and effective alleles (e) per locus were 8.6 and 4.0, respectively, and the average observed heterozygosity (o) and expected heterozygosity (e) were 0.4426 and 0.7030, respectively. The average genetic differentiation coefficient (st) and number of migrants per generation (m) of the nine.populations were 0.2452 and 0.7697, respectively, which indicates that genetic differentiation among populations was high and genetic communication was low. Cluster analysis was performed by the unweighted pair group method with arithmetic mean (UPGMA) method, based on the genetic distance among nine populations. The resulting UPGMA tree was divided into two main branches. The Shaoguan, Foshan, and Ruyuan populations were monophyletic and formed oneclade. The remaining six.populations formed the other branch, which contained three smaller branches, which were composed of the Lechang and Zhaoqing populations, the Heyuan and Huizhou populations, the Gaoyao and Fengkai populations. These results suggest that the genetic diversity of nine.populations in the Pearl River Basin was high, and there was a certain genetic differentiation among the populations, which could be used for further breeding.

; Microsatellite; Genetic diversity

GAO Fengying, E-mail: fengyinggao2011@163.com;LU Maixin, E-mail: mx-lu@163.com

S931.5

A

2095-9869(2020)04-0050-08

10.19663/j.issn2095-9869.20190731002

http://www.yykxjz.cn/

劉潔, 高風(fēng)英, 陳剛, 曹建萌, 劉志剛, 王淼, 盧邁新. 珠江流域9個(gè)大鱗副泥鰍群體遺傳多樣性的微衛(wèi)星分析. 漁業(yè)科學(xué)進(jìn)展, 2020, 41(4): 50–57

Liu J, Gao FY, Chen G, Cao JM, Liu ZG, Wang M, Lu MX. Genetic diversity analysis of nine loach () populations in the Pearl River Basin based on microsatellite markers. Progress in Fishery Sciences, 2020, 41(4): 50–57

* 現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系專項(xiàng)資金(CARS-46)資助[This work was supported by China Agriculture Research System (CARS-46)]. 劉 潔，E-mail: ljie0624@163.com

高風(fēng)英，副研究員，E-mail: fengyinggao2011@163.com；盧邁新，研究員，E-mail: mx-lu@163.com

2019-07-31,

2019-08-15

(編輯 馮小花)