凡纳滨对虾6个国内群体的遗传背景分析*

李旭鹏 栾 生 曹宝祥 罗 坤谭 建 孔张伟 孟宪红 孔 杰

凡纳滨对虾6个国内群体的遗传背景分析*

李旭鹏1,2栾 生1,2曹宝祥1罗 坤1,2谭 建1孔张伟1孟宪红1,2①孔 杰1,2①

(1. 中国水产科学研究院黄海水产研究所 青岛 266071；2. 青岛海洋科学与技术试点国家实验室海洋渔业科学与食物产出过程功能实验室 青岛 266071)

本研究选取了来自国内6个不同育苗场的凡纳滨对虾()群体，对8个SSR位点的基因型信息进行分析。结果显示，观察等位基因数(a)平均值为3.67~9.33，有效等位基因数(e)平均值为2.20~5.67，观察杂合度(o)平均值为0.12~0.71，期望杂合度(e)平均值为0.41~0.81，多态性信息含量(PIC)平均值为0.36~0.76。聚类分析结果显示，来自于同一个育苗场的对虾聚在一个分支。不同品牌来源的对虾分别聚在几个不同分支。6个群体间的遗传一致度为0.4229~0.8265，遗传距离为0.1905~0.8607。遗传分化系数(ST)是0.1837，群体内近交系数(IS)是0.2514，o

凡纳滨对虾；遗传多态性；遗传距离；SSR

凡纳滨对虾()的野生群体分布在太平洋东岸(墨西哥、厄瓜多尔、秘鲁等)。从20世纪末被引入我国开展养殖后，凡纳滨对虾已经发展成为国内最主要的海水养殖对虾种类，产量占海水养殖对虾总量的80%。国内养殖业对于凡纳滨对虾亲虾的依赖巨大。由于国内亲虾不足以满足生产，每年需要从国外进口大量成本高的亲虾。但是，通常这些亲虾的遗传背景和系谱信息无法得知。使用这些亲虾配对生产的后代苗种可能存在遗传多态性低的情况，随着培育代数的增加，近交变的更易发生，经济性状衰退明显。因此，培育国内优质凡纳滨对虾新品种意义重大。目前，已经证明遗传多态性的降低和近交情况能够增加隐性有害等位基因的表达，对选育群体的生长、存活、抗逆等性状造成近交衰退现象(Spielman, 2004; Moss, 2007; Goyard, 2008; Luo, 2014; 喻驰方等, 2013)。因此，在凡纳滨对虾遗传选育工作中，为维持育种群体的遗传多样性，避免近交，需要认清不同群体甚至个体间的分子遗传背景差异，为制定合理配种方案提供参考。

微卫星(SSR)标记是一种广泛存在于生物基因组内的共显性标记，拥有遗传多态性高的特点。研究表明，SSR适合于遗传多样性分析(Valles-Jimenez, 2004; Perez-Enriquez,2009; Rezaee, 2016)。国内在21世纪初期已有利用SSR对黄渤海不同中国对虾()群体的遗传多样性进行调查研究，区分不同地理群体间的亲缘关系距离，为种质资源利用和新品种选育提供了很好的参考(刘萍等, 2004; 孟宪红等, 2008)。由于SSR位点信息量大、分型技术成本低的优点，近年来依然被应用在对虾遗传多样性分析中(王军等, 2018; 赵志英等, 2018)。本研究使用8对SSR引物，随机对采集于我国的6个凡纳滨对虾群体进行分析，为遗传育种工作中对不同遗传背景的凡纳滨对虾资源的利用提供一些参考。

1 材料与方法

1.1 实验动物

随机从6家育苗场共收集133尾凡纳滨对虾，对其遗传背景进行分析。来自于6家育苗场的对虾分别编号为CP-1、CP-2、CP-3、CP-4、PRI和SIS群体。其中，CP-1、CP-2、CP-3、CP-4群体来自于同一个品牌。PRI和SIS群体分别来自2个不同品牌。21尾凡纳滨对虾来自CP-1群体；20尾来自CP-2群体；11尾来自CP-3群体；10尾来自CP-4群体；61尾来自PRI群体；10尾来自SIS群体。所有群体的对虾都孵化于同一年。在养殖场用消毒的剪刀取一条对虾游泳足，保存于95%酒精，带回实验室用于DNA提取和分析。

1.2 DNA提取

使用海洋动物组织基因组DNA提取试剂盒(天根)从对虾游泳足的肌肉组织中提取DNA。分别使用分光光度计和琼脂糖电泳检测DNA的质量。合格的DNA用于后期SSR分型实验。

1.3 SSR分型

使用8对已公开发表的凡纳滨对虾SSR引物(Meehan, 2003; Alcivar-Warren, 2007; Zhang, 2007)进行实验。引物序列信息详见表1。PCR反应体系包括12.5 μl Premix™ (TaKaRaVersion 2.0) (TaKaRa0.625 U、dNTPs 0.4 mmol/L、Tris-HCl 20 mmol/L、KCl 100 mmol/L、MgCl23 mmol/L) (TaKaRa)，1 μl模板DNA (200 ng/μl)，1 μl正向引物(10 μmol/L)，1 μl反向引物(10 μmol/L)， 9.5 μl去离子水。反应程序如下：98℃ 3 min；98℃ 30 s，60℃ 10 s，72℃ 1 min，25个循环；72℃延伸 3 min。使用ABI 3730XL测序仪(Applied Biosystems, 美国)获得每个SSR位点的基因型信息。

表1 引物信息

Tab.1 Primer information

1.4 遗传多样性分析

观察等位基因数(a)、有效等位基因数(e)、观察杂合度(o)、期望杂合度(e)、HWE检验、遗传一致度、遗传距离和-统计量使用POPGene(Version 1.32)软件(https://sites.ualberta.ca/~fyeh/popgene.html)进行计算。使用NTSYSpc (Version 2.10)软件(http://www. exetersoftware.com/cat/ntsyspc/ntsyspc.html)进行UPGMA聚类分析。使用Cervus (Version 3.0.7)软件(http://www. fieldgenetics.com)计算多态性信息含量(PIC)。

2 结果

2.1 遗传多态性分析

6个凡纳滨对虾群体的遗传多态性信息见表2。其中，a平均值为3.67~9.33。e平均值为2.20~5.67。o平均值为0.12~0.71。e平均值为0.41~0.81。除了在CP-4群体中的Lv12位点只检测到1种等位基因，其他SSR位点在6个群体中的等位基因数都大于1。 8个SSR位点的PIC平均值为0.36~0.76。M1103位点在CP-1、CP-3和PRI群体中都显著偏离HWE (<0.05)。Lv12位点在CP-2和PRI群体中都极显著偏离HWE(<0.01)。TuMXLv10.33位点在CP-2、CP-3、CP-4和PRI群体中都极显著偏离HWE(<0.01)。TuMXLv8.256位点在CP-4群体中显著偏离HWE (<0.05)；在CP-1、CP-2、CP-3和PRI群体中极显著偏离HWE(<0.01)。TuMXLv7.56位点在CP-4群体中显著偏离HWE(<0.05)。TuMXLv9.90位点在CP-3群体中显著偏离HWE(<0.05)；在PRI群体中极显著偏离HWE(<0.01)。TuMXLv7.21和TuMXLv9.103位点在所有群体中都没有发生偏离HWE。

表2 遗传多样性分析

Tab.2 Genetic diversity analysis

2.2 分子遗传距离分析

图1是个体的聚类分析结果。所有来自于CP-1、CP-3和CP-4群体的凡纳滨对虾集中在第1分支。其余来自于CP-2、SIS和PRI群体的凡纳滨对虾分别集中于另外3个分支。除了个体水平，在群体水平也进行了聚类分析(图2)，在2个水平的聚类分析结果是相符的。6个群体间的遗传一致度和遗传距离结果见表3。从表3可以看出，遗传一致度为0.4229~0.8265。遗传距离为0.1905~0.8607。其中，CP-1和CP-4群体间的遗传一致度最高，遗传距离最近。CP-2和SIS群体间的遗传一致度最低，遗传距离最远。

图1 个体水平聚类分析

图2 群体水平聚类分析

表3 遗传一致度和遗传距离

Tab.3 Genetic identity and genetic distance

注：对角线上方为遗传一致度，下方为遗传距离

Note: The numbers above the diagonal line are genetic identities, and the numbers below the diagonal line are genetic distances

2.3 F-统计量

对6个凡纳滨对虾群体进行-统计量分析(表4)。结果表明，8个SSR位点的群体内近交系数(IS)值为‒0.0426~0.6943。多数位点的IS>0。IS平均值是0.2514。8个SSR位点的总群体近交系数(IT)值为0.0416~0.7665，IT平均值为0.3889。8个SSR位点的遗传分化系数(ST)值为0.0808~0.3003。ST平均值为0.1837。

表4-统计量分析

Tab.4 The F-statistics analysis

2.4 等位基因的差异

韦恩图(图3)展示的是3个品牌凡纳滨对虾8个SSR位点的等位基因的分布情况。8个SSR位点共计获得等位基因个数是126。其中，只有12个等位基因类型是在3个品牌的凡纳滨对虾中都有发现。CP、PRI、SIS品牌的凡纳滨对虾各有36、29、4个等位基因类型是各自特有的，分别占各自所有等位基因类型的比例是40.9%、38.7%和12.5%。单独对CP的4个群体(CP-1、CP-2、CP-3和CP-4)进行分析(图4)，4个群体共计有88个等位基因类型，其中，4个群体共有的有10种。CP-1、CP-2、CP-3和CP-4群体各自拥有14、20、10、4种特有的等位基因，分别占各自所有等位基因类型的比例是28.0%、44.4%、23.8%和16.0%。

图3 3个品牌间SSR等位基因类型分布

图4 CP品牌的4个群体SSR等位基因类型分布

3 讨论

SSR标记广泛应用于生物的遗传多样性调查。在前期研究中，Xu等(2001)利用6对SSR引物对菲律宾的野生斑节对虾()遗传多样性进行研究，结果显示，其中1个野生群体的遗传多样性要明显低于其他3个野生群体，研究结果为了解斑节对虾的野生资源情况提供了参考。Rumisha等(2017)利用SSR标记调查了斑节对虾的遗传多态性与当地水域金属污染水平的相关性，结果显示，遗传多态性的改变与金属污染存在显著的相关性。王军等(2018)利用SSR对中国对虾人工选育群体和野生群体进行比较分析，结果显示，人工选育群体与野生群体间发生弱遗传分化。

本研究中，a平均值为3.67~9.33。在Perez- Enriquez等(2009)研究中，凡纳滨对虾SSR的a平均值为8.98。Cruz等(2004)和Rezaee等(2016)的研究中，凡纳滨对虾SSR的a值为5~10。上述研究结果和本研究结果相近。但是有趣的是，凡纳滨对虾与其他几种对虾相比，SSR的a值显得偏低很多。例如，Wang等(2016)研究表明，中国对虾的a为6~63。Zhang等(2015)研究表明，中国对虾SSR的a平均值为36.5。Luan等(2006)研究表明，日本囊对虾()SSR的a平均值为7.5~ 13.5。Xu等(2001)关于斑节对虾的报道中，SSR的a值为6~54。上述对虾种类都比凡纳滨对虾SSR的a值高。考虑到凡纳滨对虾SSR位点普遍存在a值偏低的情况，在进行遗传多样性分析时，最好能选用多态性高的SSR位点，或者适当增加分析位点的个数。已有报道中凡纳滨对虾SSR位点a值偏低的现象是否与人工选择相关也值得探索。

聚类分析结果表明，来自6个群体的凡纳滨对虾主要分为4支。其中，来自CP-1、CP-3和CP-4群体的凡纳滨对虾组成第1分支，表明这3个群体来源的对虾遗传距离较近。来自CP-2、PRI和SIS群体的对虾分别聚在其他3个分支，提示这3个群体来源的对虾之间遗传距离较远。有趣的是，CP-1、CP-2、CP-3和CP-4群体都是来源于同一个品牌，但CP-2群体与CP-1、CP-3和CP-4群体的遗传距离较远，甚至比SIS群体与CP-1、CP-3和CP-4群体的距离还远，显示该品牌群体间遗传背景较为丰富，品牌内遗传距离甚至高于品牌间遗传距离。通过聚类分析结果可以看出，使用本研究的8对SSR引物，可以将不同品牌间的个体准确区分开。但同一个品牌的群体间个体，还不能被全部准确区分开，原因可能是这些群体的亲本亲缘关系太近。

本研究中，o值为0.12~0.71。8个位点的杂合度在不同群体中的差别较大。PIC值为0.36~0.76，SSR位点的多态性较高，达到中高度水平(Botstein, 1980)，适合遗传多样性分析。ST平均值为0.1837，表明群体间已经达到高度遗传分化水平(0.15~0.25) (Balloux,2002)，18.37%的遗传分化存在于各群体间。

对于等位基因类型在不同品牌凡纳滨对虾中的分布情况，3个品牌共有的等位基因类型并不多。尤其是CP和PRI 2个品牌的凡纳滨对虾特有的等位基因类型比例较高。提示这2个品牌的对虾拥有更独特的遗传信息。但是与CP、PRI 2个品牌的对虾相比，SIS品牌的对虾特有的等位基因类型比例较低。CP品牌的4个群体中，CP-2特有的等位基因类型最高，达到44.4%，说明CP-2与其他3个CP品牌群体的遗传背景差异较大，这个结果和聚类分析的结果是一致的。假如在遗传选育中要将这6个群体引入基础选育群体，那么参考SSR分型结果，可将CP-1、CP-3、CP-4群体合并成1个群体看待，CP-2、SIS、PRI各自是3个独立群体。制定配种方案可优先考虑在这 4个新定义的群体间进行个体间配对，尽量避免群体内部个体的配对。

已有研究表明，控制近交在对虾育种中的重要性。近交对于凡纳滨对虾生长、存活、抗病毒等性状具有显著影响(Moss, 2007)。在中国对虾的研究中同样也表明，近交对于生长和存活性状具有显著影响(Luo, 2014)。在没有系谱参考的情况下，经过几代人工选育群体内近交水平就可能达到极高水平。而遗传信息丰富的基础选育群体也是保证人工选育更好开展的前提。因此，在以未知遗传背景的群体作为选育对象时，借助分子标记技术对群体、个体间的分子遗传背景进行分析，可以为配种方案制定提供重要的参考。本研究中计算的IS>0，表明群体内存在杂合度降低现象；综合考虑o

4 结论

本研究表明，不同品牌来源的凡纳滨对虾间的分子遗传距离较远，利用SSR标记可以准确区分。但是相同品牌的对虾，即使来自于不同群体，个体间的分子遗传距离可能很近。本研究检测的群体内可能已经存在一定近交现象。在遗传选育或者生产扩繁苗种过程中，为了控制近交、提高遗传多样性，在缺少候选亲本遗传背景信息时，可以尝试利用SSR标记对群体或个体的分子遗传背景进行调查，并计算分子遗传距离用于配种方案参考。

Alcivar-Warren A, Meehan-Meola D, Park SW,. ShrimpMap: A low-density, microsatellite-based linkage map of the pacific whiteleg shrimp,: Identification of sex-linked markers in linkage group 4. Journal of Shellfish Research, 2007, 26(4): 1259–1277

Balloux F, Lugon-Moulin N. 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 lengthpolymorphisms. American Journal of Human Genetics, 1980, 32(3): 314–331

Cruz P, Ibarra AM, Mejia-Ruiz H,. Genetic variability assessed by microsatellites in a breeding program of pacific white shrimp (). Marine Biotechnology, 2004, 6(2): 157–164

Goyard E, Goarant C, Ansquer D,. Cross breeding of different domesticated lines as a simple way for genetic improvement in small aquaculture industries: heterosis and inbreeding effects on growth and survival rates of the Pacific blue shrimp(). Aquaculture, 2008, 278(1–4): 43–50

Liu P, Meng XH, He YY,. Genetic diversity in three wild populations of shrimpin Yellow and Bohai Seas as revealed by microsatellite DNA. Oceanologia et Limnologia Sinica, 2004, 35(3): 252–257 [刘萍, 孟宪红, 何玉英, 等. 中国对虾()黄、渤海3个野生地理群遗传多样性的微卫星DNA分析. 海洋与湖沼, 2004, 35(3): 252–257]

Luan S, Kong J, Wang QY. Genetic variation of wild and cultured populations of the Kuruma prawn(Bate 1888) using microsatellites. Aquaculture Research, 2006, 37(8): 785–792

Luo K, Kong J, Luan S,. Effect of inbreeding on survival, WSSV tolerance and growth at the postlarval stage of experimental full-sibling inbred populations of the Chinese shrimp. Aquaculture, 2014, 420– 421(3): 32–37

Meehan D, Xu Z, Zuniga G,. High frequency and large number of polymorphic microsatellites in cultured shrimp,()(Crustacea: Decapoda). MarineBiotechnology (New York), 2003, 5(4): 311–330

Meng XH, Kong J, Wang QY,. Study on seven geographic populations of prawnbased on microsatellite DNA. Marine Fisheries Research, 2008, 29(5): 1–10 [孟宪红, 孔杰, 王清印, 等. 微卫星技术对黄、渤海海域7个不同地理群体中国对虾的遗传结构和遗传分化研究. 海洋水产研究, 2008, 29(5): 1–10]

Moss DR, Arce SM, Otoshi CA,. Effects of inbreeding on survival and growth of Pacific white shrimp(). Aquaculture, 2007, 272(S1): 30–37

Perez-Enriquez R, Hernández-Martínez F, Cruz P. Genetic diversity status of White shrimp()broodstock in Mexico. Aquaculture, 2009, 297 (1–4): 44–50

Rezaee S, Farahmand H, Nematollahi MA. Genetic diversity status of Pacific white shrimp () using SSR markers in Iran. AquacultureInternational, 2016, 24(2): 479–489

Rumisha C, Leermakers M, Elskens M,. Genetic diversity of the giant tiger prawnin relation to trace metal pollution at the Tanzanian coast. Marine Pollution Bulletin, 2017,114(2): 759–767

Spielman D, Brook BW, Briscoe DA,. Does inbreeding and loss of genetic diversity decrease disease resistance? ConservationGenetics, 2004, 5(4): 439–448

Valles-Jimenez R, Cruz P, Perez-Enriquez R. Population genetic structure of Pacific white shrimp () from Mexico to Panama: Microsatellite DNA variation. Marine Biotechnology,2004, 6(5): 475–484

Wang M, Wang W, Xiao G,. Genetic diversity analysis of spawner and recaptured populations of Chinese shrimp () during stock enhancement in the Bohai Bay based on an SSR marker. Acta Oceanologica Sinica, 2016, 35(8): 51–56

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 [王军, 王清印, 孔杰, 等. 中国明对虾人工选育群体与野生群体遗传多样性的SSR分析. 渔业科学进展, 2018, 39(2): 104–111]

Xu Z, Primavera JH, de la Pena LD,. Genetic diversity of wild and cultured Black Tiger Shrimp () in the Philippines using microsatellites. Aquaculture, 2001, 199(1): 13–40

Yu CF, Luan S, Zhang ZW,. Growth performance comparison of different inbreeding levels of families in. Progress in Fishery Sciences, 2013, 34(6): 75–79 [喻驰方, 栾生, 张志伟, 等. 中华卤虫不同近交水平家系生长性能比较. 渔业科学进展, 2013, 34(6): 75–79]

Zhao ZY, Liang LY, Bai LR. Analysis of genetic diversity among three wild populations ofusing microsatellite marker. Journal of Tropical Oceanography, 2018, 37(3): 65–72 [赵志英, 梁丽运, 白丽蓉. 斑节对虾3个野生群体遗传多样性的微卫星标记分析. 热带海洋学报, 2018, 37(3): 65–72]

Zhang K, Xiao G, Wang W,. Genetic variation analysis across six life periods in a natural population of the Chinese shrimp “” in Bohai Bay using SSR markers. Russian Journal of Marine Biology, 2015, 41(1): 10–16

Zhang L, Yang C, Zhang Y,. A genetic linkage map of Pacific white shrimp (): Sex-linked microsatellite markers and high recombination rates. Genetica,2007, 131(1): 37–49

Genetic Background Analysis of Six Groups of White Shrimpin China Using SSR Markers

LI Xupeng1,2, LUAN Sheng1,2, CAO Baoxiang1, LUO Kun1,2,TAN Jian1, KONG Zhangwei1, MENG Xianhong1,2①, KONG Jie1,2①

(1. Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071; 2. Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266071)

Controlling the inbreeding level is an important task in animal breeding. Genetic analysis such as genetic diversity and genetic distance estimate could provide information for the artificial selection of animals. In the present study, the genetic information of Pacific white shrimpfrom six different populations in China were analyzed using genotyping data of eight SSR loci. The results showed that the meanaranged from 3.67 to 9.33, the meanerange from 2.20 to 5.67, the meanoranged from 0.12 to 0.71, the meaneranged from 0.41 to 0.81, and the meanPIC ranged from 0.36 to 0.76. The cluster analysis results showed that the individual shrimps from populations of different companies were clustered in different nodes. Theindividual shrimps from three hatcheries from the same company were clustered in one node. Among the six populations of., the genetic identity ranged from 0.4229 to 0.8265, and the genetic distance ranged from 0.1905 to 0.8607. The mean value ofSTwas 0.1837. TheISvalue was0.2514. And theovalue was smaller than theevalue. The results suggested that genetic diversity between.from different populations could be high, while, genetic relationship of.in populations of the same company might be quite close. There might be an inbreeding phenomenon.

; Genetic diversity; Genetic distance; SSR

S917.4

A

2095-9869(2020)03-0103-08

10.19663/j.issn2095-9869.20190123001

孟宪红, 研究员, E-mail: mengxianhong@ysfri.ac.cn；孔 杰，研究员，E-mail: kongjie@ysfri.ac.cn

2019-01-23,

2019-04-21

* 国家自然科学基金联合基金项目(U1706203)、泰山学者种业人才团队项目、国家自然科学基金面上项目(41676148; 31572616)、现代农业产业技术体系专项资金(CARS-48)和山东省农业良种工程(2017LZN011)共同资助[This work was supported by Joint Fund of National Natural Science Foundation of China (U1706203), Taishan Scholar Program for Seed Industry, National Natural Science Foundation of China (41676148; 31572616), China Agriculture Research System (CARS-48), and Shandong Province Agricultural Seed Improvement Project (2017LZN011)]. 李旭鹏，E-mail: patrickxp@163.com

http://www.yykxjz.cn/

李旭鹏, 栾生, 曹宝祥, 罗坤, 谭建, 孔张伟, 孟宪红, 孔杰. 凡纳滨对虾6个国内群体的遗传背景分析. 渔业科学进展, 2020, 41(3): 103–110

Li XP, Luan S, Cao BX, Luo K, Tan J, Kong ZW, Meng XH, Kong J. Genetic background analysis of six groups of white shrimpin China using SSR markers. Progress in Fishery Sciences, 2020, 41(3): 103–110

MENG Xianhong, E-mail: mengxianhong@ysfri.ac.cn; KONG Jie, E-mail: kongjie@ysfri.ac.cn

(编辑 马璀艳)