留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

大口黑鲈3个养殖群体的遗传多样性分析

孙成飞 谢汶峰 胡婕 董浚键 田园园 吴灶和 叶星

引用本文:
Citation:

大口黑鲈3个养殖群体的遗传多样性分析

    作者简介: 孙成飞 (1985—),男,博士,助理研究员,从事水产遗传育种研究。E-mail: scfsec@163.com;
    通讯作者: 叶星, gzyexing@163.com
  • 中图分类号: S917.4

Genetic diversity analysis of three cultured populations of Micropterus salmoides

    Corresponding author: Xing YE, gzyexing@163.com
  • CLC number: S917.4

  • 摘要: 该研究通过毛细管电泳法筛选微卫星引物,从GenBank的51对微卫星引物中共筛选出12对具有特异性和多态性的引物,并合成荧光标记引物,对广东大口黑鲈 (Micropterus salmoides) 主养区3个养殖群体进行STR分型,以分析其遗传多样性。结果显示所检测的12个微卫星位点中,3个位点呈现高多态性,4个具有中度多态性,其余5个位点为低多态性。3个群体的遗传多样性水平偏低,期望杂合度 (He) 分别为0.312 5、0.360 6和0.328 4。群体间遗传分化指数及AMOVA分析显示,群体间遗传分化属于低水平,遗传变异有85.83%是由群体内不同个体间的差异造成的。群体间遗传距离分析显示,3个养殖群体间的遗传距离和遗传相似度比较接近。遗传结构分析表明3个养殖群体来自于同一个亚群。研究结果提示现阶段中国大口黑鲈养殖群体的遗传多样性已显著下降,因此在选育种过程应高度重视提高与维持种群遗传多样性的问题。
  • 图 1  毛细管电泳检测引物MiSaTPW011部分扩增产物

    Figure 1.  Detection for partial PCR products of MiSaTPW011 with capillary electrophoresis

    图 2  引物MiSaTPW011的3个个体的STR分型结果

    Figure 2.  STR typing results of three individuals of primer MiSaTPW011

    图 3  基于Nei's遗传距离构建的3个大口黑鲈群体的UPGMA聚类树

    Figure 3.  An UPGMA dendrograrn of three populations of M. salmoides based on Nei's genetic distance

    图 4  lnP(D)值随K的变化趋势图

    Figure 4.  Tendency chart of lnP(D) with change of Κ

    图 5  3个大口黑鲈群体的遗传结构图 (K=2)

    Figure 5.  Genetic structure map of three populations of M. salmoides (K=2)

    表 1  大口黑鲈微卫星引物序列及退火温度

    Table 1.  Primer sequences and annealing temperatures of microsatellite marker

    位点
    loci
    引物序列 (5'−3')
    primer sequence
    片段长度/bp
    size
    重复序列
    repeat sequence
    退火温度/℃
    annealingtemperature
    参考文献
    Reference
    mdo3AGGTGCTTTGCGCTACAAGT
    CTGCATGGCTGTTATGTTGG
    119~143(CA)2046.2[11]
    mdo4TCTGAACAACTGCATTTAGACTG
    CTAATCCCAGGGCAAGACTG
    151~155(CA)1148.6[11]
    mdo7TCAAACGCACCTTCACTGAC
    GTCACTCCCATCATGCTCCT
    179~205(CA)1253[11]
    lma120TGTCCACCCAAACTTAAGCC
    TAAGCCCATTCCCAATTCTCC
    204~224(GT)2860[10, 13]
    MiSaTPW011CAACATGGACGCTACTAT
    CAACCATCACATGCTTCT
    170~194(AGAT)1360[14]
    MiSaTPW038AGTCAACATTCAAACCACTTTCCCAC
    TTGTATGTTAGGATCGGCTGATTGTG
    227~312(AGAT)6+5+9+560[14]
    MiSaTPW025CCAAGGTCAGGTTTAAC
    ACCTTTGTGCTGTTCTGTC
    277~298(AGAT)1155[14]
    MiSaTPW060TATAGTTTGGTCCAGCAGGTGGCGT
    TGTGGAATGACATTTAGCCGAGGCC
    294~561(AGAT)14+10+3+2+4+9+1060[14]
    MiSaTPW116CCAAACAGCCTACAGAATGTGCCT
    AGTCCCTCGACTGAATGTGTGCAA
    191~219(AC)2155[14]
    MiSaTPW001AGTAAAGGACCACCCTTGTCCA
    GCCTGGTCATTAGGTTTCGGAG
    293~299(AC)1660[14]
    MiSaTPW058CATTCCTGAGAGACGTGGTTGCTG
    TGTGACTTGCTCTGTGAAAGGTGC
    143~169(AAGT)660[14]
    MiSaTPW154TGCGGGCTAATAGGCCTGTCTG
    GCACCACACAATCACACCTGGTAT
    182~228(AC)1855[14]
    下载: 导出CSV

    表 2  大口黑鲈12个微卫星位点的等位基因数、杂合度及多态信息含量

    Table 2.  Number of alleles, heterozygosity and polymorphic information content of 12 microsatellite loci of M. salmoides

    位点
    locus
    等位基因数
    Na
    有效等位基因数
    Ne
    Shannon指数
    I
    观测杂合度
    Ho
    期望杂合度
    He
    多态信息含量
    PIC
    mdo3 1 1.00 0.000 0.000 0.000 0.040
    mdo4 4 1.89 0.797 0.086 0.473 0.439
    mdo7 2 1.13 0.225 0.118 0.112 0.158
    lma120 4 2.89 1.211 0.552 0.658 0.604
    MISATPW011 4 2.14 0.838 0.854 0.536 0.425
    MISATPW038 4 1.72 0.658 0.523 0.420 0.442
    MISATPW025 3 2.49 0.980 0.484 0.601 0.526
    MISATPW060 6 3.52 1.427 0.667 0.720 0.693
    MISATPW116 1 1.00 0.000 0.000 0.000 0.020
    MISATPW001 2 1.12 0.219 0.115 0.109 0.102
    MISATPW058 2 1.04 0.101 0.000 0.041 0.040
    MISATPW154 4 1.65 0.749 0.208 0.397 0.363
    平均 mean 3.08 1.80 0.600 0.301 0.339 0.321
    下载: 导出CSV

    表 3  3个大口黑鲈养殖群体遗传多样性

    Table 3.  Genetic diversity of three populations of M. salmoides

    群体
    population
    等位基因数
    Na
    有效等位基因数
    Ne
    Shannon指数
    I
    观测杂合度
    Ho
    期望杂合度
    He
    多态信息含量
    PIC
    西樵 XJ 2.583 3 1.659 2 0.533 7 0.294 7 0.312 5 0.294 2
    沙头 ST 2.666 7 1.911 0 0.620 6 0.304 9 0.360 6 0.335 6
    烟南 YN 2.583 3 1.736 1 0.554 1 0.301 5 0.328 4 0.296 4
    下载: 导出CSV

    表 4  3个大口黑鲈养殖群体的遗传分化指数

    Table 4.  Fixation index (FST) of three populations of M. salmoides

    西樵 XJ 沙头 ST 烟南 YN
    西樵 XJ
    沙头 ST 0.017 9
    烟南 YN 0.012 1 0.028 9
    下载: 导出CSV

    表 5  3个大口黑鲈养殖群体AMOVA分析

    Table 5.  AMOVA analysis among three populations of M. salmoides

    变异来源
    source of variation
    自由度
    df
    平方和
    sum of squares
    方差组分
    variance component
    %
    群体间 among populations 95 192.859 0.296 76 14.17
    群体内 within populations 96 146.500 1.526 04 85.83**
    总变异 total variation 191 339.359 1.778 07
     注:**. 1 023次模拟检验后为极显著 (P<0.01)  Note: **. very significant difference after 1 023 permutation tests (P<0.01)
    下载: 导出CSV

    表 6  3个大口黑鲈养殖群体的遗传相似度S (对角线上) 和遗传距离Da (对角线下)

    Table 6.  Genetic similarity (above diagonal) and genetic distance (below diagonal) of three populations of M. salmoides

    西樵 XJ 沙头 ST 烟南 YN
    西樵 XJ 0.982 5 0.986 4
    沙头 ST 0.017 7 0.975 4
    烟南 YN 0.013 7 0.024 9
    下载: 导出CSV
  • [1] BAI J J, LUTZ-CARRILLO D J, QUAN Y C, et al. Taxonomic status and genetic diversity of cultured largemouth bass Micropterus salmoides in China[J]. Aquaculture, 2008, 278(1/2/3/4): 27-30
    [2] 农业部渔业局. 中国渔业统计年鉴[M]. 北京: 中国农业出版社, 2018: 25-33.
    [3] PRESTI F T, WASKO A P. A review of microsatellite markers and their application on genetic diversity studies in parrots[J]. Open J Genet, 2014, 4(2): 69-77 doi: 10.4236/ojgen.2014.42010
    [4] 董在杰, 刘念, 傅建军, 等. 6个野生与选育鲤群体的微卫星遗传分析[J]. 南方水产科学, 2018, 14(4): 46-55
    [5] 张小谷, 童金苟, 熊邦喜. 微卫星标记在鱼类遗传及育种研究中的应用[J]. 农业生物技术学报, 2006, 14(1): 117-121 doi: 10.3969/j.issn.1674-7968.2006.01.023
    [6] 梁素娴, 孙效文, 白俊杰, 等. 微卫星标记对中国引进加州鲈养殖群体遗传多样性的分析[J]. 水生生物学报, 2008, 32(5): 694-700
    [7] 李镕, 白俊杰, 李胜杰, 等. 大口黑鲈选育群体遗传结构的微卫星分析[J]. 广东海洋大学学报, 2010, 30(3): 11-15 doi: 10.3969/j.issn.1673-9159.2010.03.003
    [8] 樊佳佳, 白俊杰, 李胜杰, 等. 大口黑鲈微卫星DNA指纹图谱的构建和遗传结构分析[J]. 水生生物学报, 2012, 36(4): 600-609
    [9] 黄太寿, 白俊杰, 李胜杰, 等. 我国加州鲈养殖现状和绿色发展问题探讨[J]. 中国水产, 2017(12): 44-47
    [10] NEFF B D, FU P, GROSS M R. Microsatellite evolution in sunfish (Centrarchidae)[J]. Can J Fish Aquat Sci, 1999, 56(7): 1198-1205 doi: 10.1139/f99-068
    [11] MALLOY T P, van den BUSSCHE R A, COUGHLIN W D, et al. Isolation and characterization of microsatellite loci in smallmouth bass, Micropterus dolomieu (Teleostei: Centrarchidae), and cross-species amplification in spotted bass, M. punctulatus[J]. Mol Ecol, 2000, 9(11): 1946-1948 doi: 10.1046/j.1365-294x.2000.01096-16.x
    [12] RICHARD M, THORPE R S. Highly polymorphic microsatellites in the lacertid Gallotia galloti from the western Canary Islands[J]. Mol Ecol, 2000, 9(11): 1919-1920 doi: 10.1046/j.1365-294x.2000.01052.x
    [13] LUTZ-CARRILLO D J, NICE C C, BONNER T H, et al. Admixture analysis of Florida largemouth bass and northern largemouth bass using microsatellite loci[J]. Trans Am Fish Soc, 2006, 135(3): 779-791 doi: 10.1577/T04-221.1
    [14] LUTZ-CARRILLO D J, HAGEN C, DUECK L A, et al. Isolation and characterization of microsatellite loci for Florida largemouth bass, Micropterus salmoides floridanus, and other micropterids[J]. Mol Ecol Resour, 2008, 8(1): 178-184 doi: 10.1111/j.1471-8286.2007.01917.x
    [15] TAMURA K, PETERSON D, PETERSON N, et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods[J]. Mol Biol Evol, 2011, 28(10): 2731-2739 doi: 10.1093/molbev/msr121
    [16] KALINOWSKI S T, TAPER M L, MARSHALL T C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment[J]. Mol Ecol, 2007, 16(5): 1099-1106 doi: 10.1111/j.1365-294X.2007.03089.x
    [17] EXCOFFIER L, LISCHER H E. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows[J]. Mol Ecol Resour, 2010, 10(3): 564-567 doi: 10.1111/men.2010.10.issue-3
    [18] EVANNO G, REGNAUT S, GOUDET J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study[J]. Mol Ecol, 2005, 14(8): 2611-2620 doi: 10.1111/mec.2005.14.issue-8
    [19] 曾庆凯, 孙成飞, 董浚键, 等. 翘嘴鳜3个不同群体的遗传多样性分析[J]. 基因组学与应用生物学, 2017, 36(8): 3241-3250
    [20] 孙成飞, 叶星, 董浚键, 等. 罗氏沼虾6个养殖群体遗传多样性的微卫星分析[J]. 南方水产科学, 2015, 11(2): 20-26 doi: 10.3969/j.issn.2095-0780.2015.02.003
  • [1] 杜博龚世园童馨黄桂菊喻达辉 . 皱纹盘鲍和盘鲍南方养殖群体遗传变异的微卫星分析. 南方水产科学, 2007, 3(6): 22-29.
    [2] 孙成飞叶星董浚键田园园梁健辉 . 罗氏沼虾6个养殖群体遗传多样性的微卫星分析. 南方水产科学, 2015, 11(2): 20-26. doi: 10.3969/j.issn.2095-0780.2015.02.003
    [3] 朱冰樊佳佳白俊杰姜鹏 . 金草鱼与中国4个草鱼群体的微卫星多态性比较分析. 南方水产科学, 2017, 13(2): 51-58. doi: 10.3969/j.issn.2095-0780.2017.02.007
    [4] 孙立元郭华阳朱彩艳马振华江世贵张殿昌 . 卵形鲳鲹育种群体遗传多样性分析. 南方水产科学, 2014, 10(2): 67-71. doi: doi:10.3969/j.issn.2095-0780.2014.02.010
    [5] 范嗣刚王婧璇黄桂菊刘宝锁郭奕惠喻达辉1 . 合浦珠母贝选育家系的遗传多样性分析. 南方水产科学, 2016, 12(5): 90-96. doi: 10.3969/j.issn.2095-0780.2016.05.011
    [6] 李敏李玉芳张鹏陈作志 . 基于线粒体控制区序列的南海圆舵鲣种群遗传结构分析. 南方水产科学, 2016, 12(4): 88-95. doi: 10.3969/j.issn.2095-0780.2016.04.011
    [7] 柳明喻达辉黄桂菊卢传亮 . 中国海南三亚大珠母贝不同年代种群的遗传变异研究. 南方水产科学, 2011, 7(1): 26-31. doi: 10.3969/j.issn.2095-0780.2011.01.004
    [8] 赵彦花区又君温久福李加儿周慧 . 基于微卫星标记的黄唇鱼遗传多样性研究. 南方水产科学, 2019, 15(4): 1-6. doi: 10.12131/20180261
    [9] 唐首杰毕详王成辉张飞明张友良谢志强 . 团头鲂3个选育群体遗传潜力的微卫星分析. 南方水产科学, 2017, 13(2): 59-68. doi: 10.3969/j.issn.2095-0780.2017.02.008
    [10] 董在杰刘念傅建军朱文彬王兰梅苏胜彦 . 6个野生与选育鲤群体的微卫星遗传分析. 南方水产科学, 2018, 14(4): 46-55. doi: 10.3969/j.issn.2095-0780.2018.04.006
    [11] 吴文婷张磊宋志民赵翠汤晓荣刘涛陈伟洲 . 5个紫菜属物种丝状体的微卫星初步遗传分析. 南方水产科学, 2012, 8(4): 29-36. doi: 10.3969/j.issn.2095-0780.2012.04.005
    [12] 黄小帅徐煜胡晓娟徐武杰苏浩昌文国樑杨铿曹煜成 . 利用微卫星标记分析7个凡纳滨对虾引进群体子一代的遗传多样性. 南方水产科学, 2019, 15(1): 54-62. doi: 10.12131/20180135
    [13] 杨喜书章群余帆洋吕金磊底晓丹邵伟军黄镇宇卢丽锋 . 华南6水系与澜沧江-湄公河攀鲈线粒体ND2基因的遗传多样性分析. 南方水产科学, 2017, 13(3): 43-50. doi: 10.3969/j.issn.2095-0780.2017.03.006
    [14] 李胜杰樊佳佳姜鹏白俊杰孙建国吴建开费志平 . 大口黑鲈HSC70-1基因多态性及其双倍型与生长性状的关联分析. 南方水产科学, 2018, 14(6): 74-80. doi: 10.12131/20180086
    [15] 熊丹李敏陈作志李永振李玉芳黄梓荣 . 南海短尾大眼鲷的种群遗传结构分析. 南方水产科学, 2015, 11(2): 27-34. doi: 10.3969/j.issn.2095-0780.2015.02.004
    [16] 匡天旭帅方敏陈蔚涛李新辉 . 西江鲫的遗传多样性与群体结构. 南方水产科学, 2018, 14(5): 29-35. doi: 10.3969/j.issn.2095-0780.2018.05.004
    [17] 李敏张鹏李玉芳陈森张魁孔啸兰陈作志 . 南海扁舵鲣种群遗传结构和遗传多样性评价. 南方水产科学, 2015, 11(5): 82-89. doi: 10.3969/j.issn.2095-0780.2015.05.010
    [18] 莫艳秀王晓清莫永亮 . 长吻鮠遗传多样性的RAPD分析. 南方水产科学, 2010, 6(6): 77-80. doi: 10.3969/j.issn.1673-2227.2010.06.014
    [19] 刘伟苏胜彦董在杰张建桥马良骁李灵玲曲疆奇袁新华 . 3个鲤群体的微卫星标记与生长性状相关性分析. 南方水产科学, 2012, 8(3): 17-24. doi: 10.3969/j.issn.2095-0780.2012.03.003
    [20] 朱彩艳叶卫夏军红符云周发林江世贵 . 广东1个鲮原种群体的种质特征及遗传多样性分析. 南方水产科学, 2005, 1(4): 1-5.
  • 加载中
图(5)表(6)
计量
  • 文章访问数:  1016
  • HTML全文浏览量:  439
  • PDF下载量:  20
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-11
  • 录用日期:  2018-11-01
  • 网络出版日期:  2018-12-20
  • 刊出日期:  2019-04-01

大口黑鲈3个养殖群体的遗传多样性分析

    作者简介:孙成飞 (1985—),男,博士,助理研究员,从事水产遗传育种研究。E-mail: scfsec@163.com
    通讯作者: 叶星, gzyexing@163.com
  • 1. 中国水产科学研究院珠江水产研究所,农业农村部热带亚热带水产资源利用与养殖重点实验室,广东 广州 510380
  • 2. 上海海洋大学,上海 201306
  • 3. 仲恺农业工程学院,广东 广州 510380

摘要: 该研究通过毛细管电泳法筛选微卫星引物,从GenBank的51对微卫星引物中共筛选出12对具有特异性和多态性的引物,并合成荧光标记引物,对广东大口黑鲈 (Micropterus salmoides) 主养区3个养殖群体进行STR分型,以分析其遗传多样性。结果显示所检测的12个微卫星位点中,3个位点呈现高多态性,4个具有中度多态性,其余5个位点为低多态性。3个群体的遗传多样性水平偏低,期望杂合度 (He) 分别为0.312 5、0.360 6和0.328 4。群体间遗传分化指数及AMOVA分析显示,群体间遗传分化属于低水平,遗传变异有85.83%是由群体内不同个体间的差异造成的。群体间遗传距离分析显示,3个养殖群体间的遗传距离和遗传相似度比较接近。遗传结构分析表明3个养殖群体来自于同一个亚群。研究结果提示现阶段中国大口黑鲈养殖群体的遗传多样性已显著下降,因此在选育种过程应高度重视提高与维持种群遗传多样性的问题。

English Abstract

  • 大口黑鲈 (Micropterus salmoides),俗称加州鲈,隶属鲈形目、鲈亚目、太阳鱼科、黑鲈属。原自然分布于北美洲,于1983年从中国台湾引入广东省[1]。大口黑鲈适温范围较广,生长快且易起捕,同时其肉质鲜美、无肌间刺,因而深受养殖者和消费者欢迎,在我国多个省市与地区得以迅速推广养殖。近几年大口黑鲈人工配合饲料的突破与应用,使大口黑鲈养殖区域与规模迅速扩大。2017年,我国养殖大口黑鲈年产量已达45.7万吨[2]。但近两三年来养殖大口黑鲈抗病性能下降,病害问题频发。大口黑鲈引进中国已三十多年,虽然这期间也从原产地进行过为数不多的小规模引进,但相对于养殖规模而言,亲本的种质更新极为有限。因此抗病性能下降是否与种质质量与遗传多样性下降有关,进一步开展选择育种的种质遗传基础如何,都亟需对中国现有大口黑鲈养殖群体的遗传多样性进行评估。

    微卫星 (microsatellite) 标记,又称简单重复序列 (simple sequence repeats, SSR) 或短串联重复序列 (simple tandem repeats, STR),是均匀分布于真核生物基因组中的简单重复序列,由2~6个核苷酸的串联重复片段构成,具有分布广泛、共显性遗传、多态性高等特点,广泛应用于多个生命学科领域[3-4]。20世纪90年代初,微卫星标记开始应用于鱼类遗传及育种研究中[5]。利用微卫星研究大口黑鲈的遗传多样性在近年有所报道。梁素娴等[6]用21对微卫星引物对广东的3个大口黑鲈养殖群体进行遗传多样性分析发现,3个不同群体平均杂合度处于中等偏低水平 (期望杂合度He分别为0.375、0.403和0.368),遗传多样性不高。李镕等[7]利用11对微卫星引物对大口黑鲈第2~第4代选育群体的遗传结构进行分析发现,与F2相比,F4的遗传多样性减少8.76%,表明选育群体的遗传多样性在逐步下降。樊佳佳等[8]用43个微卫星标记对目前国内养殖的大口黑鲈群体、2009年引进的佛罗里达亚种、2010年引进的佛罗里达亚种和北方亚种这4个群体进行遗传多样性检测,发现国内养殖的大口黑鲈群体遗传多样性水平明显低于国外新引进的大口黑鲈群体。以上研究表明,我国大口黑鲈养殖群体出现遗传多样性显著下降、种质退化及种群结构单一等问题。

    广东省是中国大口黑鲈成鱼的主养区和主要苗种培育基地。2016年广东省大口黑鲈养殖面积约5 333.4 hm2,产量为23.5万吨,占全国总产量的62.8%。广东省大口黑鲈成鱼养殖主要集中在珠江三角洲的佛山市,其养殖产量近20万吨。佛山市的顺德区和南海区有大口黑鲈鱼苗培育专业村,仅南海区九江镇的大口黑鲈鱼苗生产量就接近40亿尾,约占广东全省生产量的75%,占全国生产量的60%[9]。因此佛山地区的大口黑鲈养殖群体极具代表性,研究其养殖群体的遗传多样性和遗传结构,可以反映出中国大口黑鲈养殖群体的遗传多样性现状。本研究选取广东省佛山市3个不同地方的大口黑鲈养殖群体作为研究对象,利用微卫星标记技术对其进行遗传多样性分析,旨在为当前大口黑鲈养殖群体的种质资源调查和进一步的选育种工作提供支撑。

    • 实验所需大口黑鲈样品分别为佛山市南海区西樵镇的西樵群体 (编号为XJ)、九江镇沙头片区的沙头群体 (编号为ST)、九江镇九江片区烟南村的烟南群体 (编号为YN),这3个养殖群体的苗种来自不同苗种场。从3个养殖群体中各随机选取32尾个体,分别剪取尾鳍,–20 ℃保存于95%的乙醇中,用于DNA的提取。

    • 利用天根生化科技 (北京)有限公司的海洋动物组织基因组DNA提取试剂盒,提取大口黑鲈尾鳍的基因组DNA。用1%的琼脂糖凝胶进行电泳,并用酶标仪检测其质量和浓度,−20 ℃冰箱保存。

    • 根据GenBank及文献上登录的大口黑鲈微卫星引物序列[10-14],选择并设计合成51对引物,进行预扩增,并用Qsep100全自动核酸蛋白分析系统 (中国,台湾) 检测引物的特异性和多态性,最终确定其中的12对重复性好、特异性强和多态性高的微卫星引物 (表1) 用于研究。引物委托生工生物工程 (上海) 股份有限公司进行合成。各对引物分别带有TAMRA、HEX和FAM的特异性荧光标记,用于扩增和测算目的条带的大小。

      位点
      loci
      引物序列 (5'−3')
      primer sequence
      片段长度/bp
      size
      重复序列
      repeat sequence
      退火温度/℃
      annealingtemperature
      参考文献
      Reference
      mdo3AGGTGCTTTGCGCTACAAGT
      CTGCATGGCTGTTATGTTGG
      119~143(CA)2046.2[11]
      mdo4TCTGAACAACTGCATTTAGACTG
      CTAATCCCAGGGCAAGACTG
      151~155(CA)1148.6[11]
      mdo7TCAAACGCACCTTCACTGAC
      GTCACTCCCATCATGCTCCT
      179~205(CA)1253[11]
      lma120TGTCCACCCAAACTTAAGCC
      TAAGCCCATTCCCAATTCTCC
      204~224(GT)2860[10, 13]
      MiSaTPW011CAACATGGACGCTACTAT
      CAACCATCACATGCTTCT
      170~194(AGAT)1360[14]
      MiSaTPW038AGTCAACATTCAAACCACTTTCCCAC
      TTGTATGTTAGGATCGGCTGATTGTG
      227~312(AGAT)6+5+9+560[14]
      MiSaTPW025CCAAGGTCAGGTTTAAC
      ACCTTTGTGCTGTTCTGTC
      277~298(AGAT)1155[14]
      MiSaTPW060TATAGTTTGGTCCAGCAGGTGGCGT
      TGTGGAATGACATTTAGCCGAGGCC
      294~561(AGAT)14+10+3+2+4+9+1060[14]
      MiSaTPW116CCAAACAGCCTACAGAATGTGCCT
      AGTCCCTCGACTGAATGTGTGCAA
      191~219(AC)2155[14]
      MiSaTPW001AGTAAAGGACCACCCTTGTCCA
      GCCTGGTCATTAGGTTTCGGAG
      293~299(AC)1660[14]
      MiSaTPW058CATTCCTGAGAGACGTGGTTGCTG
      TGTGACTTGCTCTGTGAAAGGTGC
      143~169(AAGT)660[14]
      MiSaTPW154TGCGGGCTAATAGGCCTGTCTG
      GCACCACACAATCACACCTGGTAT
      182~228(AC)1855[14]

      表 1  大口黑鲈微卫星引物序列及退火温度

      Table 1.  Primer sequences and annealing temperatures of microsatellite marker

    • PCR反应总体积20 μL,含有模板DNA 1 μL、上下游引物 (10 μmol·L–1) 各1 μL、ddH2O 7 μL、Premix rTaq (TaKaRa,中国大连) 10 μL,在96孔板中进行PCR反应。扩增程序为94 ℃预变性3 min;94 ℃变性30 s,退火 (各引物退火温度参照表1中的退火温度) 30 s,72 ℃延伸30 s,35个循环;72 ℃延伸10 min,4 ℃保存。

    • 引物筛选的预扩增产物通过Qsep100全自动核酸蛋白分析系统进行毛细管电泳检测。将PCR扩增产物室温放置2~3 min,开盖后按要求放在样品托上进行毛细管电泳。电泳参数为6 kV下进样2 s再电泳300 s。电泳完毕后,用Q-ANALYZER软件对结果进行矫正与分析,对扩增产物分子量的差异大小进行个体各微卫星基因座的分型。

      荧光引物的扩增产物通过STR基因分型检测,由生工生物工程 (上海) 股份有限公司完成。使用毛细管电泳分离STR片段,激光诱导的荧光检测,由设备检测并数据化存储,自动数据采集并分型。荧光颜色分离,计算片段大小,确定各个等位座的基因分型。整个过程通过ABI 310遗传分析仪 (美国) 完成。

    • 通过STR分型确定标记在个体中的基因型,运用POPGENE 32软件进行统计分析,计算各个群体微卫星位点的等位基因数 (Na)、有效等位基因数 (Ne)、观测杂合度 (Ho)、He、Shannon指数 (I)、遗传相似度 (S) 及群体间遗传距离 (Da)。运用FSTAT软件算出遗传分化指数 (FST)。根据Nei氏遗传距离的大小,运用MEGA 5软件构建系统树[15]。运用CERVUS 3.0软件统计出多态信息含量 (polymorphism information content, PIC)[16]。运用ARLEQUIN 3.1软件进行群体间FST计算、群体分子方差分析 (AMOVA)[17]。利用STRUCTURE 2.3软件分析群体的遗传结构[18]

    • 从不同地区采集到的大口黑鲈养殖群体中,随机选取12个大口黑鲈样本,利用51对微卫星引物进行PCR扩增,通过Qsep100全自动核酸蛋白分析系统进行毛细管电泳检测产物(图1),筛选出12对多态性好、特异性强的微卫星引物用于大群体的扩增分析 (表1)。

      图  1  毛细管电泳检测引物MiSaTPW011部分扩增产物

      Figure 1.  Detection for partial PCR products of MiSaTPW011 with capillary electrophoresis

    • 用筛选出的12对微卫星引物扩增3个大口黑鲈群体共96个个体的DNA,将产物进行STR分型 (图2)。结果表明,各对引物在3个群体中均能扩增到目的条带,并体现了程度不一的多态性(表2)。其中,Na为1~6,Ne为1~3.52,I为0~1.427,Ho为0~0.854,He为0~0.720,PIC为0.02~0.693。

      图  2  引物MiSaTPW011的3个个体的STR分型结果

      Figure 2.  STR typing results of three individuals of primer MiSaTPW011

      位点
      locus
      等位基因数
      Na
      有效等位基因数
      Ne
      Shannon指数
      I
      观测杂合度
      Ho
      期望杂合度
      He
      多态信息含量
      PIC
      mdo3 1 1.00 0.000 0.000 0.000 0.040
      mdo4 4 1.89 0.797 0.086 0.473 0.439
      mdo7 2 1.13 0.225 0.118 0.112 0.158
      lma120 4 2.89 1.211 0.552 0.658 0.604
      MISATPW011 4 2.14 0.838 0.854 0.536 0.425
      MISATPW038 4 1.72 0.658 0.523 0.420 0.442
      MISATPW025 3 2.49 0.980 0.484 0.601 0.526
      MISATPW060 6 3.52 1.427 0.667 0.720 0.693
      MISATPW116 1 1.00 0.000 0.000 0.000 0.020
      MISATPW001 2 1.12 0.219 0.115 0.109 0.102
      MISATPW058 2 1.04 0.101 0.000 0.041 0.040
      MISATPW154 4 1.65 0.749 0.208 0.397 0.363
      平均 mean 3.08 1.80 0.600 0.301 0.339 0.321

      表 2  大口黑鲈12个微卫星位点的等位基因数、杂合度及多态信息含量

      Table 2.  Number of alleles, heterozygosity and polymorphic information content of 12 microsatellite loci of M. salmoides

      3个大口黑鲈群体的遗传多样性见表3。其中3个群体的Na为2.583 3~2.666 7,Ne为1.659 2~1.911 0,I为0.533 7~0.620 6,Ho为0.294 7~0.304 9,He为0.312 5~0.360 6,PIC为0.294 2~0.335 6。沙头群体的6个遗传多样性参数略高于西樵和烟南群体,西樵和烟南群体的比较接近,且西樵群体的遗传多样性参数最低。

      群体
      population
      等位基因数
      Na
      有效等位基因数
      Ne
      Shannon指数
      I
      观测杂合度
      Ho
      期望杂合度
      He
      多态信息含量
      PIC
      西樵 XJ 2.583 3 1.659 2 0.533 7 0.294 7 0.312 5 0.294 2
      沙头 ST 2.666 7 1.911 0 0.620 6 0.304 9 0.360 6 0.335 6
      烟南 YN 2.583 3 1.736 1 0.554 1 0.301 5 0.328 4 0.296 4

      表 3  3个大口黑鲈养殖群体遗传多样性

      Table 3.  Genetic diversity of three populations of M. salmoides

    • 3个大口黑鲈养殖群体间的遗传分化指数 (FST)为0.012 1~0.028 9 (表4),属低水平遗传分化 (FST<0.05)。其中西樵和烟南群体间的遗传分化水平最低 (FST=0.012 1),沙头群体和烟南群体间的遗传分化水平最高 (FST=0.028 9,表4)。

      西樵 XJ 沙头 ST 烟南 YN
      西樵 XJ
      沙头 ST 0.017 9
      烟南 YN 0.012 1 0.028 9

      表 4  3个大口黑鲈养殖群体的遗传分化指数

      Table 4.  Fixation index (FST) of three populations of M. salmoides

      AMOVA分析显示,14.17%的遗传变异来自群体间,85.83%的遗传变异来自群体内 (P<0.01)。结果表明,遗传变异主要存在个体间,个体间的遗传变异程度远大于群体间 (表5)。

      变异来源
      source of variation
      自由度
      df
      平方和
      sum of squares
      方差组分
      variance component
      %
      群体间 among populations 95 192.859 0.296 76 14.17
      群体内 within populations 96 146.500 1.526 04 85.83**
      总变异 total variation 191 339.359 1.778 07
       注:**. 1 023次模拟检验后为极显著 (P<0.01)  Note: **. very significant difference after 1 023 permutation tests (P<0.01)

      表 5  3个大口黑鲈养殖群体AMOVA分析

      Table 5.  AMOVA analysis among three populations of M. salmoides

      3个群体间Nei's遗传距离 (Da)介于0.013 7~0.024 9,遗传相似度 (S)介于0.975 4~0.986 4 (表6),其中西樵群体和沙头群体的遗传距离最近 (Da=0.013 7)、遗传相似度最高 (S=0.986 4),而烟南群体和沙头群体的遗传距离最远 (Da=0.024 9)、遗传相似度最低 (S=0.975 4)。由遗传距离大小构建了UPGMA聚类树,聚类分析结果显示,3个大口黑鲈群体归为两大类,西樵与烟南群体为一类,沙头群体单独为一类 (图3)。

      西樵 XJ 沙头 ST 烟南 YN
      西樵 XJ 0.982 5 0.986 4
      沙头 ST 0.017 7 0.975 4
      烟南 YN 0.013 7 0.024 9

      表 6  3个大口黑鲈养殖群体的遗传相似度S (对角线上) 和遗传距离Da (对角线下)

      Table 6.  Genetic similarity (above diagonal) and genetic distance (below diagonal) of three populations of M. salmoides

      图  3  基于Nei's遗传距离构建的3个大口黑鲈群体的UPGMA聚类树

      Figure 3.  An UPGMA dendrograrn of three populations of M. salmoides based on Nei's genetic distance

    • 运用STRUCTURE软件对3个大口黑鲈群体的遗传结构进行分析,执行K=2~7的假设,设定10次重复,结果显示对数似然函数值lnP(D)随亚群数K的增加而降低,K=1时,lnP(D)最大,从而确定出最适宜的亚群数为1 (图4)。K=2时,3个群体的遗传结构图几乎没有分别 (图5),也说明3个养殖群体都来自一个亚群。

      图  4  lnP(D)值随K的变化趋势图

      Figure 4.  Tendency chart of lnP(D) with change of Κ

      图  5  3个大口黑鲈群体的遗传结构图 (K=2)

      Figure 5.  Genetic structure map of three populations of M. salmoides (K=2)

    • 本研究采用Qsep100全自动核酸蛋白分析系统,通过毛细管电泳筛选微卫星引物,在51对微卫星引物中挑选出12对多态性相对较好的引物进行大群体遗传多样性的检测。在检测阶段则是采用ABI310遗传分析仪通过荧光标记的引物进行PCR扩增,扩增产物经毛细管电泳后进行STR分析。Qsep100全自动核酸蛋白分析系统和ABI310遗传分析仪均采用毛细管电泳系统,2个系统都具有操作简单、速度快、准确性高的优点。不同之处在于,Qsep100系统在电泳介质中加入能结合核酸的荧光物质,利用激光诱导其发光来检测核酸[19],ABI310系统则是先对引物进行荧光标记,使得PCR产物带有荧光,毛细管电泳时直接检测产物荧光来确定产物大小[20]。Qsep100系统成本较高的部分是带有荧光的电泳介质 (即卡夹),ABI310系统成本较高的部分是荧光引物的合成,Qsep100系统的总成本略低。但ABI310系统毛细管的峰图则更清晰,更容易排除杂带的干扰,检测准确性更高。因此综合2套系统的优势,本研究在引物筛选阶段采用Qsep100系统,在检测分型阶段则采用ABI310系统,可以在保证准确性的前提下提高检测效率,节省成本。

      对本研究筛选的12对微卫星引物的遗传参数分析发现,12对引物中有5对引物属于低度多态 (PIC<0.25),4对引物属于中度多态 (0.25<PIC<0.5),只有3对引物属于高度多态 (PIC>0.5)。其中有2对引物 (mdo3和MISATPW116)只有1个等位基因 (Na=1),其He均为0,说明这2对引物没有多态性。Malloy等[11]和Lutz-Carrillo等[13-14]通过磁珠富集法筛选到大量的大口黑鲈微卫星引物,本研究选择其中的51对引物用于扩增,并从中筛选出12对多态性高且特异性好的微卫星引物。在Malloy等[11]和Lutz-Carrillo等[13-14]的研究中,这12对引物大都属于高度多态性,其中有10对引物He为0.55~0.94,其遗传多样性参数均显著高于本研究所获得的参数值。另外,在他们的研究中报道有高度多态性的引物,大部分在本研究的3个群体中已不再具有多态性[10-14]。因此,从本研究12对引物在大口黑鲈3个养殖群体扩增所获得的多样性参数来看,本研究中的3个养殖群体的遗传多样性较Malloy等报道过的群体已显著下降。推测当前养殖的大口黑鲈群体大都是由最初引进的大口黑鲈群体繁育的后代,经过多年的近交繁育,遗传多样性已显著下降。

      对本研究选取的广东佛山3个大口黑鲈养殖群体的遗传多样性分析发现,3个群体的He为0.312 5~0.360 6。而樊佳佳等[8]用微卫星标记对国内养殖群体和3个引进群体进行遗传多样性检测,其He为0.454 9~0.613 8,高于本研究的分析结果。本研究所选用的MISATPW001和mdo72位点,PIC值已下降为0.102和0.158,而在樊佳佳等[8]的结果中,国内养殖群体中这2个位点的PIC分别为0.141 1和0.239 2,这均说明近几年我国大口黑鲈养殖群体遗传多样性下降明显。另外,3个群体之间,沙头群体 (ST)的遗传多样性最高 (He=0.360 6),各遗传参数均高于另外2个群体,可能与苗种来源有关。

      本研究的3个大口黑鲈养殖群体遗传分化指数 (FST)为0.012 1~0.028 9,遗传分化水平较低,说明3个群体遗传差异不明显。3个群体Da介于0.013 7~0.024 9,非常接近,另外,运用STRUCTURE软件进行遗传结构分析也表明3个养殖群体应来自于一个亚群。这佐证了上述观点,即当前大口黑鲈养殖群体极可能均由最初引进的大口黑鲈群体繁育而来。当前大口黑鲈养殖群体结构单一,遗传多样性低,抗病性能下降的问题除了与养殖模式及管理有关外,也极可能与种质退化、遗传多样性下降有关。由于遗传多样性下降,导致其对环境的适应能力变差,进而更容易被病原等侵袭。因此应高度重视提高与维持大口黑鲈种群遗传多样性的问题,并通过分子选育种技术进行优良品种的培育,以促进大口黑鲈养殖产业的健康可持续发展。

参考文献 (20)

目录

    /

    返回文章
    返回