留言板

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

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

基于线粒体ND2基因序列的少鳞遗传多样性研究

郑德育 郭易佳 杨天燕 高天翔 郑瑶 袁冬皓 斯舒谨

引用本文:
Citation:

基于线粒体ND2基因序列的少鳞遗传多样性研究

    作者简介: 郑德育(1998—),男,本科生,从事海洋生物资源与环境研究。E-mail: 723259288@qq.com;
    通讯作者: 杨天燕, hellojelly1130@163.com
  • 中图分类号: S 913.4

Genetic diversity analysis of Sillago japonica based on mitochondrial DNA ND2 gene

    Corresponding author: Tianyan YANG, hellojelly1130@163.com ;
  • CLC number: S 913.4

  • 摘要: 以莱州、胶南、舟山、厦门、汕头和北海6个群体119尾少鳞 (Sillago japonica)为研究对象,采用PCR扩增测序获得长度为450 bp的线粒体DNA NADH脱氢酶亚基2 (ND2)基因片段,共检测到77个变异位点,其中简约信息位点30个,单变异位点28个,无碱基缺失。119条序列定义了61个单倍型,平均单倍型多样性(Hd)和核苷酸多样性(π)分别为0.945 3±0.015 5和0.009 718±0.005 445。6个群体间的平均遗传距离为0.008 3,遗传分化指数FST均小于0.05,各群体间无显著遗传分化。AMOVA分析得出少鳞的遗传变异主要来自于种群内个体间(99.96%)。中性检验的Tajima's D和Fu's Fs统计值均为负值且显著偏离中性,核苷酸不配对分布图呈现明显的单峰分布,表明少鳞历史上经历了群体扩张事件,估算扩张时间大约在 (0.12~0.29)百万年前的第四纪更新世晚期。
  • 图 1  少鳞采样地点图

    Figure 1.  Sampling sites of S. japonica

    图 2  基于ND2序列构建的少鳞NJ树

    Figure 2.  NJ phylogenetic tree of haplotypes based on ND2 sequences of S. japonica

    图 3  基于少鳞ND2基因单倍型的核苷酸不配对分布图

    Figure 3.  Mismatch distribution based on ND2 gene haplotypes of S. japonica

    表 1  基于线粒体ND2基因的少鳞6个种群遗传多样性

    Table 1.  Genetic diversity in seven S. japonica populations based on mtDNA ND2 gene

    种群
    population
    样本数
    number of individuals
    单倍型数目
    number of haplotypes
    多态位点数目
    number of polymorphic sites
    单倍型多样性
    Hd
    核苷酸多样性
    π
    莱州 LZ1912160.924 0±0.045 80.005 788±0.003 723
    胶南 JN2013170.957 9±0.025 50.007 800±0.004 739
    舟山 ZS2013340.947 4±0.032 30.013 019±0.007 349
    厦门 XM2017210.984 2±0.020 50.010 282±0.005 993
    汕头 ST2014480.915 8±0.054 60.014 225±0.007 953
    北海 BH2013190.910 5±0.053 80.007 678±0.004 677
    总计 total11951770.945 3±0.015 50.009 718±0.005 445
    下载: 导出CSV

    表 2  少鳞不同群体间平均遗传距离

    Table 2.  Average genetic distances among different S. japonica populations

    种群
    population
    北海
    BH
    胶南
    JN
    莱州
    LZ
    汕头
    ST
    厦门
    XM
    胶南 JN0.007 74
    莱州 LZ0.006 850.006 88
    汕头 ST0.008 640.008 510.007 66
    厦门 XM0.009 110.009 020.008 500.010 02
    舟山 ZS0.008 040.007 950.007 420.008 980.009 17
    下载: 导出CSV

    表 3  基于ND2序列单倍型频率的少鳞群体FST (对角线下) 和相应P (对角线上)

    Table 3.  F-Statistics (below diagonal) and FST P value (above diagonal) from haplotype frequencies of S. japonica

    种群
    population
    北海
    BH
    胶南
    JN
    莱州
    LZ
    汕头
    ST
    厦门
    XM
    舟山
    ZS
    BH0.625 000.201 170.898 440.535 160.564 45
    JN−0.008 440.434 570.693 360.457 030.517 58
    LZ 0.000 930.012 840.638 670.018 550.021 48
    ST−0.005 06−0.011 29−0.002 920.355 470.654 30
    XM−0.001 54−0.004 820.044 380.003 940.826 17
    ZS−0.002 41−0.005 850.035 13−0.013 66−0.014 16
    下载: 导出CSV

    表 4  少鳞群体线粒体ND2基因的AMOVA分析

    Table 4.  AMOVA analysis of S. japonica populations based on mtDNA ND2 gene

    变异来源
    source of variation
    自由度
    df
    平方和
    SS
    方差分量
    variance components
    变异百分比
    percentage of variation
    种群间 among populations 5 9.720 0.000 78Va 0.04
    种群内 within populations 113 217.911 1.928 41Vb 99.96
    总计 total 118 227.630 1.929 19
    固定值 fixation index (FST) 0.000 41
     注:Va. 组群间方差、Vb. 组群内群体间方差
     Note: Va. variance between groups; Vb. variance between groups within the group
    下载: 导出CSV

    表 5  少鳞群体中性检验及分化时间

    Table 5.  Neutral test and differentiation time of fit for S. japonica

    中性检验
    neutral test
    莱州
    LZ
    胶南
    JN
    舟山
    ZS
    厦门
    XM
    汕头
    ST
    北海
    BH
    Tajima's D−1.918 06−1.678 23−0.721 48−1.244 16−2.098 64−1.378 17
    P0.012 000.033 000.272 000.090 000.008 000.072 00
    Fu's Fs−7.368 12−6.792 29−3.610 54−12.466 49−4.380 38−6.687 65
    P0.000 000.001 000.050 000.000 000.030 000.003 00
    Tau2.23.73.54.35.33.4
    T (Mya)0.120.210.190.240.290.18
    下载: 导出CSV

    表 6  少鳞种群中性检验

    Table 6.  Neutral test of population for S. japonica

    数量
    number
    Tajima's D Fu's Fs偏离方差
    SSD
    Rg
    DP FsP
    119 −2.223 54 0 −25.916 43 0 0.001 89 0.020 24
    下载: 导出CSV
  • [1] NELSON J S. Fishes of the world[M]. New Jersey: John Wiley & Sons, 2016: 503.
    [2] MCKAY R J. An annotated and illustrated catalogue of the sillago, smelt or Indo-Pacific whiting species known to date[R]. Rome: FAO, 1992: 1-83.
    [3] SANO J. Fisheries management by spawning per recruit analysis and yield per recruit analysis for Sillago japonica around the coastal waters of Itoshima [Japan] area[R]. Bulletin of Fukuoka Fisheries & Marine Technology Research Center, Fukuoka, 2004: 46-47.
    [4] SHIMASAKI Y, OSHIMA Y, INOUE S, et al. Effect of tributyltin on reproduction in Japanese whiting, Sillago japonica[J]. Mar Environ Res, 2006, 62(S): S245-S248.
    [5] OOZEKI Y, HWANG P P, HIRANO R. Larval development of the Japanese whiting, Sillago japonica[J]. Jpn J Ichthyol, 1992, 39(1): 59-66.
    [6] KASHIWAGI M, KONDO S, YOSHIDA W, et al. Effects of temperature and salinity on hatching success of Japanese whiting Sillago japonica eggs[J]. Suisan Zoshoku, 2000, 48(4): 637-642.
    [7] SULISTIONO S, WATANABE S, YOKOTA M. Reproduction of the Japanese whiting, Sillago japonica, in Tateyama Bay[J]. Aquacult Sci, 1999, 47(2): 209-214.
    [8] RAHMAN S M, MAJHI S K, SUZUKI T A, et al. Suitability of cryoprotectants and impregnation protocols for embryos of Japanese whiting Sillago japonica[J]. Cryobiology, 2008, 57(2): 170-174. doi: 10.1016/j.cryobiol.2008.08.002
    [9] RAHMAN S M, STRUESSMANN C A, SUZUKI T, et al. Electroporation enhances permeation of cryoprotectant (dimethyl sulfoxide) into Japanese whiting (Sillago japonica) embryos[J]. Theriogenology, 2013, 79(5): 853-858. doi: 10.1016/j.theriogenology.2013.01.002
    [10] SULISTIONO S, YOKOTA M, KITADA S, et al. Age and growth of Japanese whiting Sillago japonica in Tateyama Bay[J]. Fish Sci, 1999, 65(1): 117-122.
    [11] ARAYAMA K, IMAI H, KOHNO H, et al. Early life story of Japanese whiting Sillago japonica occurring in the surf zone of sandy beaches Tateyama Bay, central Japan[J]. Nippon Suisan Gakkaishi, 2003, 69(3): 359-367. doi: 10.2331/suisan.69.359
    [12] 杨亚峰, 宋娜, 肖家光, 等. 莱州湾少鳞的形态特征描述[J]. 齐鲁渔业, 2016, 33(10): 8-10.
    [13] 潘晓哲, 高天翔. 基于耳石形态的属鱼类鉴别[J]. 动物分类学报, 2010, 35(4): 799-805.
    [14] 薛泰强, 杜宁, 高天翔. 基于线粒体COI及Cytb基因的4种科鱼类系统发育研究[J]. 中国海洋大学学报(自然科学版), 2010, 40(S1): 91-98.
    [15] 肖家光. 基于线粒体基因组全序列的属鱼类系统发育研究[D]. 青岛: 中国海洋大学, 2015: 44-53.
    [16] GAO T X, YANG T Y, YANAGIMOTO T, et al. Levels and patterns of genetic variation in Japanese whiting (Sillago japonica) based on mitochondrial DNA control region[J]. Mitochondrial DNA Pt A, 2019, 30(1): 172-183. doi: 10.1080/24701394.2018.1467411
    [17] 王林燕. 基于微卫星标记的中国和少鳞群体遗传学研究[D]. 青岛: 中国海洋大学, 2014: 36-64.
    [18] VELLEND M, GEBER M A. Connections between species diversity and genetic diversity[J]. Ecol Lett, 2005, 8(7): 767-781. doi: 10.1111/ele.2005.8.issue-7
    [19] JUAN Y, ZHONG Z Q, FEN L. Mitochondrial DNA and its application to the molecular population genetics of fish[J]. Ecol Sci, 2008, 27(4): 272-276.
    [20] WILSON A C, CANN R L, CARR S M, et al. Mitochondrial DNA and two perspectives on evolutionary genetics[J]. Biol J Linn Soc, 2010, 26(4): 375-400.
    [21] 杨喜书, 章群, 余帆洋, 等. 华南6水系与澜沧江-湄公河攀鲈线粒体ND2基因的遗传多样性分析[J]. 南方水产科学, 2017, 13(3): 43-50. doi: 10.3969/j.issn.2095-0780.2017.03.006
    [22] 阮燕如. 基于线粒体ND2基因序列的华南地区斑鳢遗传多样性研究[D]. 广州: 暨南大学, 2014: 52-53.
    [23] 伊西庆. 中国东部6个大型湖泊翘嘴鲌(Culter alburnus)遗传多样性的线粒体ND2基因序列分析[D]. 广州: 暨南大学, 2009: 31-35.
    [24] GEORGE A L, CALDIERARO J B, CHARTRAND K M. Population genetics of the blue shiner, Cyprinella caerulea[J]. Southeast Nat, 2008, 7(4): 637-650. doi: 10.1656/1528-7092-7.4.637
    [25] VERISSIMO A, MCDOWELL J R, GRAVES J E. Genetic population structure and connectivity in a commercially exploited and wide-ranging deepwater shark, the leafscale gulper (Centrophorus squamosus)[J]. Mar Freshw Res, 2012, 63(6): 505-512. doi: 10.1071/MF11237
    [26] SAMBROOK J, FRITSCH E F, MANIATIS T. Molecular cloning: a laboratory manual[M]. New York: Cold Spring Harbor Laboratory Press, 1982: 76-82.
    [27] CLEWLEY J P. Macintosh sequence analysis software. DNAStar's LaserGene[J]. Mol Biotechnol, 1995, 3(3): 221-224. doi: 10.1007/BF02789332
    [28] ROZAS J, FERRERMATA A, SÁNCHEZDELBARRIO J C, et al. DnaSP 6: DNA sequence polymorphism analysis of large datasets[J]. Mol Biol Evol, 2017, 34(12): 3299-3302. doi: 10.1093/molbev/msx248
    [29] 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
    [30] KUMAR S, STECHER G, TAMURA K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Mol Biol Evol, 2016, 33(7): 1870-1874. doi: 10.1093/molbev/msw054
    [31] ROGERS A R, HARPENDING H. Population growth makes waves in the distribution of pairwise genetic differences[J]. Mol Biol Evol, 1992, 9(3): 552-569.
    [32] BERMINGHAM E S, MCCAFFERTY A. Molecular systematics of fishes[M]. New York: Academic Press, 1997: 113-126.
    [33] FERGUSON J H. On the use of genetic divergence for identifying species[J]. Biol J Linn Soc, 2015, 75(4): 509-516.
    [34] SKIBINSKI D F. DNA tests of neutral theory: applications in marine genetics[M]. Berlin: Springer Netherlands, 2000: 137-152.
    [35] FU Y X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection[J]. Genetics, 1997, 147(2): 915-925.
    [36] BONIN A, NICOLE F, POMPANON F, et al. Population adaptive index: a new method to help measure intraspecific genetic diversity and prioritize populations for conservation[J]. Conserv Biol, 2007, 21(3): 697-708. doi: 10.1111/cbi.2007.21.issue-3
    [37] NEI M. Molecular evolutionary genetics[M]. New York: Columbia University Press,1987: 92-145.
    [38] 王秀亮. 玉筋鱼群体遗传多样性及其适应进化研究[D]. 舟山: 浙江海洋大学, 2017: 24-25.
    [39] XU S Y, SUN D R, SONG N, et al. Local adaptation shapes pattern of mitochondrial population structure in Sebastiscus marmoratus[J]. Environ Biol Fish, 2017, 100(7): 763-774. doi: 10.1007/s10641-017-0602-5
    [40] WRIGHT S. Evolution and the genetics of populations[M]. Chicago: University of Chicago Press, 1968: 76-79.
    [41] HEWITT G M. Genetic consequences of climatic oscillations in the quaternary[J]. Philos T R Soc B, 2004, 359(1442): 183-195. doi: 10.1098/rstb.2003.1388
    [42] 刘海松. 地貌学及第四纪地质学[M]. 北京: 地质出版社, 2013: 10-11.
    [43] 沈浪, 陈小勇, 李媛媛. 生物冰期避难所与冰期后的重新扩散[J]. 生态学报, 2002, 22(11): 1983-1990. doi: 10.3321/j.issn:1000-0933.2002.11.026
  • [1] 王春晓高风英卢迈新刘志刚朱华平叶星 . 2个尼罗罗非鱼群体GHSR基因5侧翼序列的多态性及其遗传多样性分析. 南方水产科学, 2015, 11(1): 18-25. doi: 10.3969/j.issn.2095-0780.2015.01.003
    [2] 匡天旭帅方敏陈蔚涛李新辉 . 西江鲫的遗传多样性与群体结构. 南方水产科学, 2018, 14(5): 29-35. doi: 10.3969/j.issn.2095-0780.2018.05.004
    [3] 杨喜书章群余帆洋吕金磊底晓丹邵伟军黄镇宇卢丽锋 . 华南6水系与澜沧江-湄公河攀鲈线粒体ND2基因的遗传多样性分析. 南方水产科学, 2017, 13(3): 43-50. doi: 10.3969/j.issn.2095-0780.2017.03.006
    [4] 孙立元郭华阳朱彩艳马振华江世贵张殿昌 . 卵形鲳鲹育种群体遗传多样性分析. 南方水产科学, 2014, 10(2): 67-71. doi: doi:10.3969/j.issn.2095-0780.2014.02.010
    [5] 孙成飞叶星董浚键田园园梁健辉 . 罗氏沼虾6个养殖群体遗传多样性的微卫星分析. 南方水产科学, 2015, 11(2): 20-26. doi: 10.3969/j.issn.2095-0780.2015.02.003
    [6] 朱彩艳叶卫夏军红符云周发林江世贵 . 广东1个鲮原种群体的种质特征及遗传多样性分析. 南方水产科学, 2005, 1(4): 1-5.
    [7] 李莉好喻达辉黄桂菊杜博符云童馨郭奕惠叶卫 . 吉富罗非鱼不同选育群体的遗传多样性. 南方水产科学, 2007, 3(5): 40-48.
    [8] 孙成飞谢汶峰胡婕董浚键田园园吴灶和叶星 . 大口黑鲈3个养殖群体的遗传多样性分析. 南方水产科学, 2019, 15(2): 64-71. doi: 10.12131/20180203
    [9] 原居林朱俊杰王高学 . 秦岭细鳞鲑黑河种群和湑水河种群的遗传多样性分析. 南方水产科学, 2009, 5(3): 63-66. doi: 10.3969/j.issn.1673-2227.2009.03.011
    [10] 孙奉玉宋忠魁赵鹏聂振平苏琼王芳宇 . 广西沿海及其邻近海区拟穴青蟹群体遗传多样性的RAPD分析. 南方水产科学, 2012, 8(2): 30-35. doi: 10.3969/j.issn.2095-0780.2012.02.005
    [11] 符云钟金香颉晓勇叶卫林碧海陈辉崇张汉华 . 罗非鱼3个养殖群体的遗传多样性及特异性AFLP标记研究. 南方水产科学, 2008, 4(6): 50-55.
    [12] 颉晓勇苏天凤陈文张志颜远义江世贵 . 凡纳滨对虾6个养殖群体遗传多样性的比较分析. 南方水产科学, 2008, 4(6): 42-49.
    [13] 黄小帅徐煜胡晓娟徐武杰苏浩昌文国樑杨铿曹煜成 . 利用微卫星标记分析7个凡纳滨对虾引进群体子一代的遗传多样性. 南方水产科学, 2019, 15(1): 54-62. doi: 10.12131/20180135
    [14] 杜博龚世园童馨黄桂菊喻达辉 . 皱纹盘鲍和盘鲍南方养殖群体遗传变异的微卫星分析. 南方水产科学, 2007, 3(6): 22-29.
    [15] 周发林江世贵姜永杰黄建华马之明 . 海南三亚斑节对虾野生种群线粒体16SrRNA基因和控制区序列的多态性研究. 南方水产科学, 2006, 2(6): 13-18.
    [16] 李敏张鹏李玉芳陈森张魁孔啸兰陈作志 . 南海扁舵鲣种群遗传结构和遗传多样性评价. 南方水产科学, 2015, 11(5): 82-89. doi: 10.3969/j.issn.2095-0780.2015.05.010
    [17] 莫艳秀王晓清莫永亮 . 长吻鮠遗传多样性的RAPD分析. 南方水产科学, 2010, 6(6): 77-80. doi: 10.3969/j.issn.1673-2227.2010.06.014
    [18] 范嗣刚王婧璇黄桂菊刘宝锁郭奕惠喻达辉1 . 合浦珠母贝选育家系的遗传多样性分析. 南方水产科学, 2016, 12(5): 90-96. doi: 10.3969/j.issn.2095-0780.2016.05.011
    [19] 荣朝振祖国掌胡建华孙守旗孙棠丽 . 泥鳅线粒体DNA控制区结构分析及遗传多样性研究. 南方水产科学, 2011, 7(5): 55-62. doi: 10.3969/j.issn.2095-0780.2011.05.009
    [20] 赵彦花区又君温久福李加儿周慧 . 基于微卫星标记的黄唇鱼遗传多样性研究. 南方水产科学, 2019, 15(4): 127-132. doi: 10.12131/20180261
  • 加载中
图(3)表(6)
计量
  • 文章访问数:  2025
  • HTML全文浏览量:  750
  • PDF下载量:  14
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-02-25
  • 录用日期:  2019-04-24
  • 网络出版日期:  2019-05-31
  • 刊出日期:  2019-10-01

基于线粒体ND2基因序列的少鳞遗传多样性研究

    作者简介:郑德育(1998—),男,本科生,从事海洋生物资源与环境研究。E-mail: 723259288@qq.com
    通讯作者: 杨天燕, hellojelly1130@163.com
  • 浙江海洋大学水产学院,浙江 舟山 316022

摘要: 以莱州、胶南、舟山、厦门、汕头和北海6个群体119尾少鳞 (Sillago japonica)为研究对象,采用PCR扩增测序获得长度为450 bp的线粒体DNA NADH脱氢酶亚基2 (ND2)基因片段,共检测到77个变异位点,其中简约信息位点30个,单变异位点28个,无碱基缺失。119条序列定义了61个单倍型,平均单倍型多样性(Hd)和核苷酸多样性(π)分别为0.945 3±0.015 5和0.009 718±0.005 445。6个群体间的平均遗传距离为0.008 3,遗传分化指数FST均小于0.05,各群体间无显著遗传分化。AMOVA分析得出少鳞的遗传变异主要来自于种群内个体间(99.96%)。中性检验的Tajima's D和Fu's Fs统计值均为负值且显著偏离中性,核苷酸不配对分布图呈现明显的单峰分布,表明少鳞历史上经历了群体扩张事件,估算扩张时间大约在 (0.12~0.29)百万年前的第四纪更新世晚期。

English Abstract

  • 少鳞 (Sillago japonica)是鲈形目、科、属的一种小型海洋经济鱼类,主要分布在红海、印度洋、太平洋西部和南部,盛产于我国南海和东海的浅水区或河口区[1-2]。少鳞肉质鲜美、营养丰富,深受人们喜爱,也是近岸渔业捕捞对象和游钓鱼种。20世纪末开始,日本学者围绕少鳞渔业管理[3]、繁殖发育学[4-7]、组织胚胎学[8-9]和资源生物学[10-11]领域开展了大量基础性研究工作,近年来国内有关其形态学和遗传学方面的研究主要集中在外部形态特征描述[12]、耳石形态多元统计分析[13]、分子系统学研究[14-15]、线粒体控制区[16]和微卫星分子标记开发[17]。遗传多样性反映了生物对环境变化的适应能力,是物种多样性的基础和决定性因素,体现在表型、染色体、蛋白质和DNA多个层次水平的变异上,开展遗传多样性研究对于探讨物种的系统进化、遗传分化和资源保护都具有十分重要的意义[18]

    线粒体DNA (Mitochondrial DNA,mtDNA)具有分子量小、结构简单、排列紧凑,几乎不发生重组且遵循母系遗传的特性,目前已广泛应用在鱼类分子系统学、生物地理学和种群遗传学等诸多领域[19-20]。在mtDNA包含的13个蛋白质编码基因中,NADH脱氢酶亚基2 (NADH dehydrogenase subunit 2,ND2)基因由于进化速率相对较快,能较好地反映种内与群体间的遗传变异及系统分类地位,正逐渐应用于鱼类系统发育和种群划分研究中[21-25]。本研究采集了中国沿海6个地理种群的少鳞样本,基于线粒体ND2基因序列开展其群体遗传多样性的研究,以期为少鳞种质资源的保护提供理论依据。

    • 本研究所用少鳞样本分别于2017—2018年间采自莱州(LZ,19尾)、胶南(JN,20尾)、舟山(ZS,20尾)、厦门(XM,20尾)、汕头(ST,20尾)、北海(BH,20尾),共计119尾样品(图1)。所有样品取背部新鲜肌肉组织存放于5 mL离心管,加入无水乙醇固定,4 ℃保存备用。

      图  1  少鳞采样地点图

      Figure 1.  Sampling sites of S. japonica

    • 剪取适量已浸泡在无水乙醇中的背部肌肉组织,采用传统的酚-氯仿法[26]提取少鳞基因组DNA,用无水乙醇沉淀和75%乙醇润洗后经自然风干的DNA溶解于100 μL灭菌水中,置于4 ℃保存备用。

      参考少鳞线粒体基因组全序列(GenBank登录号KR363149),采用Primer primer 6.0软件设计扩增ND2基因的引物,分别为上游引物ND2-1F (5'-CACG AACGCCCCTATACTCA-3')和下游引物ND2-1R (5'-CCTGGTAGGTTGTTAGGGGA-3')。PCR总反应体系为25 μL,包括0.25 μL浓度为5 U·μL−1Taq酶,各1 μL正反向引物(10 μmol·L−1),2 μL的dNTPs (2.5 mmol·L−1),2.5 μL的10×buffer缓冲液(含Mg2+),1 μL的模板DNA (50~100 ng),以及17.25 μL的灭菌双蒸水。PCR反应程序为95 ℃预变性5 min;95 ℃变性0.5 min,52 ℃退火0.5 min,72 ℃延伸1 min,总计循环35次;72 ℃延伸5 min。用1.5%琼脂糖凝胶电泳检测PCR扩增产物,凝胶成像系统拍照后,挑选条带亮度较高的PCR产物送至上海美吉生物有限公司进行测序。

    • 所获序列均由DNAstar[27]软件包进行编辑、校对和排序;采用DnaSP 6.0[28]和Arlequin[29]软件进行多态位点数目(number of polymorphic sites)、单倍型数目(number of haplotype)、核苷酸多样性(nucleotide diversity)和单倍型多样性(haplotype diversity)等计算;使用MEGA 7.0[30]软件基于Kimura 2-parameter model参数计算少鳞群体间平均遗传距离、碱基组成和碱基转换/颠换,并采用邻接法(neighbor-joining,NJ)基于1 000次Bootstrap重抽样构建单倍型系统发育树;运用Network软件基于中间连接法(median-joining)绘制少鳞单倍型网络图。使用Arlequin[29]软件统计少鳞群体间的遗传分化水平,分子变异水平分析(analysis of molecular variance analysis,AMOVA)种群间遗传变异情况;使用Tajima's D和Fu's Fs中性检验以及核苷酸不配对分布(mismatch distribution)来估算种群的历史动态变化。使用公式τ=2ut[31]估算种群扩张时间,其中τ为扩张时间参数(Tau),u为进化速率,u=μk,其中μ为每个碱基的变异速率,k为所分析序列的长度,t表示自扩张以来所经历的代数,扩张时间T=t×代时,μ参考Bermingham和Mccafferty[32]的研究,采用每百万年2%的突变速率。

    • 使用DNAstar软件包中SeqMan程序进行序列比对,得到长度为450 bp的ND2基因片段。119尾少鳞样本平均碱基组成为A=20.96%,T=30.38%,G=14.05%,C=34.6%,A+T含量(51.34%)略高于G+C (48.66%),其中C含量最高,G含量最低,碱基组成偏向于嘧啶碱基。所获得序列中多态位点数目77个,占所获得序列长度的19.84%,没有发现碱基的插入和缺失。119条序列共检测到61个单倍型,从各群体所含单倍型分布来看,厦门群体检测到的单倍型数量最多,而莱州群体检测到的单倍型数量最少。遗传多样性参数显示,厦门群体的单倍型多样性最高(0.984 2±0.020 5),北海群体单倍型多样性最低(0.910 5±0.053 8);汕头群体核苷酸多样性最高(0.014 225±0.007 953),莱州群体核苷酸多样性最低(0.005 788±0.003 723,表1)。

      种群
      population
      样本数
      number of individuals
      单倍型数目
      number of haplotypes
      多态位点数目
      number of polymorphic sites
      单倍型多样性
      Hd
      核苷酸多样性
      π
      莱州 LZ1912160.924 0±0.045 80.005 788±0.003 723
      胶南 JN2013170.957 9±0.025 50.007 800±0.004 739
      舟山 ZS2013340.947 4±0.032 30.013 019±0.007 349
      厦门 XM2017210.984 2±0.020 50.010 282±0.005 993
      汕头 ST2014480.915 8±0.054 60.014 225±0.007 953
      北海 BH2013190.910 5±0.053 80.007 678±0.004 677
      总计 total11951770.945 3±0.015 50.009 718±0.005 445

      表 1  基于线粒体ND2基因的少鳞6个种群遗传多样性

      Table 1.  Genetic diversity in seven S. japonica populations based on mtDNA ND2 gene

    • 在MEGA软件中采用Kimura 2-parameter模型计算少鳞6个群体间的遗传距离(表2),种群间的遗传距离在0.006 85~0.010 02,其中厦门和汕头的遗传距离最大 (0.010 02);莱州和北海的遗传距离最小 (0.006 85),遗传分化尚未达到种级水平[33]

      种群
      population
      北海
      BH
      胶南
      JN
      莱州
      LZ
      汕头
      ST
      厦门
      XM
      胶南 JN0.007 74
      莱州 LZ0.006 850.006 88
      汕头 ST0.008 640.008 510.007 66
      厦门 XM0.009 110.009 020.008 500.010 02
      舟山 ZS0.008 040.007 950.007 420.008 980.009 17

      表 2  少鳞不同群体间平均遗传距离

      Table 2.  Average genetic distances among different S. japonica populations

      少鳞不同群体间遗传多样性指数(FST)分析发现两两群体间的FST介于−0.014 16 ~0.044 38,除莱州和厦门、舟山群体之间差异显著外,其余群体间P均大于0.05 (表3)。采用分子方差分析(AMOVA)对6个少鳞群体遗传变异来源和结构进行统计(表4)。结果表明,少鳞遗传变异主要来自群体内个体间(99.96%),其P检验也不显著(P=0.426 20±0.015 26),表明少鳞群体间存在广泛的基因交流,使得遗传分化程度较低、遗传差异不显著。

      种群
      population
      北海
      BH
      胶南
      JN
      莱州
      LZ
      汕头
      ST
      厦门
      XM
      舟山
      ZS
      BH0.625 000.201 170.898 440.535 160.564 45
      JN−0.008 440.434 570.693 360.457 030.517 58
      LZ 0.000 930.012 840.638 670.018 550.021 48
      ST−0.005 06−0.011 29−0.002 920.355 470.654 30
      XM−0.001 54−0.004 820.044 380.003 940.826 17
      ZS−0.002 41−0.005 850.035 13−0.013 66−0.014 16

      表 3  基于ND2序列单倍型频率的少鳞群体FST (对角线下) 和相应P (对角线上)

      Table 3.  F-Statistics (below diagonal) and FST P value (above diagonal) from haplotype frequencies of S. japonica

      变异来源
      source of variation
      自由度
      df
      平方和
      SS
      方差分量
      variance components
      变异百分比
      percentage of variation
      种群间 among populations 5 9.720 0.000 78Va 0.04
      种群内 within populations 113 217.911 1.928 41Vb 99.96
      总计 total 118 227.630 1.929 19
      固定值 fixation index (FST) 0.000 41
       注:Va. 组群间方差、Vb. 组群内群体间方差
       Note: Va. variance between groups; Vb. variance between groups within the group

      表 4  少鳞群体线粒体ND2基因的AMOVA分析

      Table 4.  AMOVA analysis of S. japonica populations based on mtDNA ND2 gene

      从GenBank数据库下载近缘种多鳞 (S. sihama,GenBank登录号KR363150)和中国 (S. sinica,GenBank登录号KR363151) 作为外群,采用邻接法(neighbor-joining,NJ) 构建系统发育树 (图2),聚类树中不同群体的个体相互交织,没有形成明显的分支,种群之间存在广泛的基因交流。

      图  2  基于ND2序列构建的少鳞NJ树

      Figure 2.  NJ phylogenetic tree of haplotypes based on ND2 sequences of S. japonica

    • 在1 000次模拟抽样的情况下,用Arlequin软件对6个群体的少鳞进行中性检验(表5表6)和核苷酸不配对分析(图3)。当Tajima's D为负,且在统计学上达到显著水平,表明种群分化偏离中性突变理论模型,预示着群体在进化过程中不仅受到随机漂变因素影响,还有可能经历过大规模的群体扩张或受到过选择压力、瓶颈效应的影响[34]。Fu's Fs检验中,如果Fs大于零,表明种群趋于稳定,反之则表明种群趋于扩张[35]。结果表明,无论是从群体还是个体角度来看,Tajima's D和Fu's Fs均为负值且显著偏离中性(P<0.05),表明少鳞在历史进化过程中经历过群体扩张事件。此外,偏离方差(SSD)和Raggedness index (Rg)均较小,统计检验不显著,进一步表明群体扩增假说的成立[29]。核苷酸不配对分布图呈单峰且中性检验差异显著,表明少鳞群体近期可能受到了瓶颈效应的打击。根据核苷酸不配对分布得到的Tau,少鳞性成熟时间为1年,代时取1,计算得到少鳞群体发生扩张的时间大约在 (12~29) 万年前的第四纪更新世(Pleistocene)晚期。

      中性检验
      neutral test
      莱州
      LZ
      胶南
      JN
      舟山
      ZS
      厦门
      XM
      汕头
      ST
      北海
      BH
      Tajima's D−1.918 06−1.678 23−0.721 48−1.244 16−2.098 64−1.378 17
      P0.012 000.033 000.272 000.090 000.008 000.072 00
      Fu's Fs−7.368 12−6.792 29−3.610 54−12.466 49−4.380 38−6.687 65
      P0.000 000.001 000.050 000.000 000.030 000.003 00
      Tau2.23.73.54.35.33.4
      T (Mya)0.120.210.190.240.290.18

      表 5  少鳞群体中性检验及分化时间

      Table 5.  Neutral test and differentiation time of fit for S. japonica

      数量
      number
      Tajima's D Fu's Fs偏离方差
      SSD
      Rg
      DP FsP
      119 −2.223 54 0 −25.916 43 0 0.001 89 0.020 24

      表 6  少鳞种群中性检验

      Table 6.  Neutral test of population for S. japonica

      图  3  基于少鳞ND2基因单倍型的核苷酸不配对分布图

      Figure 3.  Mismatch distribution based on ND2 gene haplotypes of S. japonica

    • 单倍型多样性(h)和核苷酸多样性(π)是衡量物种遗传多样性丰富度的重要指标,也是物种生存适应和发展进化的前提。物种的遗传多样性越高,对环境的适应能力越强,对变化也有更强的调整能力[36]。本研究表明中国沿海少鳞群体遗传多样性较高(h=0.945 3,π=0.009 718),与薛泰强等[14]和王林燕[17]分别采用控制区序列和微卫星分子标记对不同地理群体少鳞遗传结构的分析结果相似,都呈现出相对丰富的遗传多样性水平。这可能与ND2基因本身具有较高的突变速率有关,另外种群的高单倍型多样性与种群的数量、所处环境的变化以及生活习性有密切联系[37],少鳞种群繁殖时间较长,单次后代产生数量大,能够弥补因捕食、捕捞等带来的种群缺失,使其存在较高的遗传多样性。

      莱州群体核苷酸多样性指数最低,且与厦门、舟山群体间的遗传差异显著。由于莱州位于采样点最北端,受渤海沿岸流和黄海冷水团的共同作用,所处海域水温相对较低,而浙闽沿岸受亚热带海洋性季风气候影响,水温普遍偏高,推测温度变化可能对少鳞遗传变异产生一定影响。鱼类作为变温性动物,温度变化对其分布和生长发育都具有重要影响。温度变化限制了鱼类成体的扩散并促进其适应局限的生存环境,使得彼此间交流的概率大幅降低,这种差异经长期积累并最终导致遗传分化。有文献资料表明,一些海洋鱼类的遗传多样性分布受到温度和纬度变化的影响,如王秀亮[38]对西北太平洋玉筋鱼(Ammodytes personatus) Cyt b基因研究发现,玉筋鱼单倍型频率在地理分布上与水温变化具有相关性,因而推测水域环境温度的改变以及洋流系统对玉筋鱼种群扩散及种群地理分布格局可能有着重要影响;Xu等[39]对中日沿海5个群体褐菖鲉(Sebastiscus marmoratus) ATP6和Cyt b基因所编码的氨基酸替换率进行比较时发现,温度可能是影响群体遗传结构和遗传分化的潜在环境因素。

      遗传多样性指数FST是评价物种群体间遗传分化尺度的标准。Wright[40]的研究发现,若FST>0.25,表明群体间存在高度分化,0.15<FST<0.25,表示群体间存在中度分化,0.05<FST<0.15表示群体间存在低度分化,而FST<0.05则表明群体间没有分化。本研究发现少鳞各群体间FST均小于0.05,结合平均遗传距离的大小来看,暗示群体间未发生明显的遗传分化,这与Gao等[16]的研究结果一致。推测少鳞这类小型鱼类群体数量大、分布范围广,加之早期营浮游生活时间长、本身具备较强迁移能力,基因交流频繁造成遗传的同质化,使得种群间遗传分化水平普遍偏低,阻碍了种群的分化。此外,基于邻接法构建的系统关系树显示来自不同群体的个体彼此相互交错分布,没有呈现显著的谱系结构,也反映了少鳞遗传分化水平较低这一结论。

      Hewitt[41]通过研究证实,更新世的剧烈气候波动对地球上动植物遗传多样性和分布格局产生了重要影响。基于线粒体控制区序列估算出中日韩不同地理种群少鳞发生群体扩张的时间大约在 (0.423~1.06)百万年前的晚更新世时期[16],本研究对ND2基因的核苷酸不配对分析得出少鳞群体分化时间大约在 (0.12~0.29)百万年前的第四纪冰川期-间冰期旋回[42]。尽管线粒体不同序列变异速率存在一定差异,但推算出扩张时间均处于第四纪更新世晚期,这一时期气候的急剧变化,全球性大幅度气温变冷,中高纬形成大面积冰盖,大气环流和洋流的变化直接影响动植物生长、演化和分布,导致大量生物种群的迁徙或者灭绝[43],使得少鳞等海洋鱼类的空间分布格局以及遗传结构发生了较大变化。

      鱼类种群结构的时空分布和动态变化与其所处水域的环境因子有着密切关系,过度捕捞、水体污染、工程建设以及外来物种入侵等都会使鱼类资源和遗传多样性遭到破坏。基于线粒体ND2基因分析的结果显示,我国沿海少鳞种群遗传多样性丰富程度较高,种质资源尚处于相对安全状态,具有一定的资源开发与利用潜力,但由于多遗传标记结合对于揭示种群进化历史和谱系地理格局的形成往往更全面、更准确。因此,在今后的渔业资源管理过程中,还需要综合采用多种分子生物学标记技术来监测少鳞的遗传多样性水平,以实现其资源多样性的合理保护和可持续利用。

参考文献 (43)

目录

    /

    返回文章
    返回