南海短尾大眼鲷的种群遗传结构分析

熊丹, 李敏, 陈作志, 李永振, 李玉芳, 黄梓荣

熊丹, 李敏, 陈作志, 李永振, 李玉芳, 黄梓荣. 南海短尾大眼鲷的种群遗传结构分析[J]. 南方水产科学, 2015, 11(2): 27-34. DOI: 10.3969/j.issn.2095-0780.2015.02.004
引用本文: 熊丹, 李敏, 陈作志, 李永振, 李玉芳, 黄梓荣. 南海短尾大眼鲷的种群遗传结构分析[J]. 南方水产科学, 2015, 11(2): 27-34. DOI: 10.3969/j.issn.2095-0780.2015.02.004
XIONG Dan, LI Min, CHEN Zuozhi, LI Yongzhen, LI Yufang, HUANG Zirong. Genetic structure of Priacanthus macracanthus population from the South China Sea[J]. South China Fisheries Science, 2015, 11(2): 27-34. DOI: 10.3969/j.issn.2095-0780.2015.02.004
Citation: XIONG Dan, LI Min, CHEN Zuozhi, LI Yongzhen, LI Yufang, HUANG Zirong. Genetic structure of Priacanthus macracanthus population from the South China Sea[J]. South China Fisheries Science, 2015, 11(2): 27-34. DOI: 10.3969/j.issn.2095-0780.2015.02.004

南海短尾大眼鲷的种群遗传结构分析

基金项目: 

农业部财政重大专项 NFZX2013

国家科技支撑计划项目 2013BAD13B06

中央级公益性科研院所基本科研业务费专项资金(中国水产科学研究院南海水产研究所资助项目) 2014TS23

农业部财政专项“南海北部近海渔业资源调查” 20141005

详细信息
    作者简介:

    熊丹(1989-),女,硕士研究生,从事鱼类种群分析研究。E-mail:xiongdan_2013@163.com

    通讯作者:

    黄梓荣(1963-),男,副研究员,从事渔业资源评估和渔业生态研究。E-mail:hzr0114@163.com

  • 中图分类号: Q347;S931.5

Genetic structure of Priacanthus macracanthus population from the South China Sea

  • 摘要:

    在南海北部陆架和南海南沙西南陆架6个地理位点采集224尾短尾大眼鲷(Priacanthus macracanthus)样品,采用线粒体控制区(mtDNA D-loop区)序列分析了其种群遗传结构,探讨了6个采样位点之间短尾大眼鲷种群归属关系。所得224条同源序列D-loop片段(729 bp)中检测到101个突变位点和94个核苷酸多态位点,定义了172种单倍型;遗传多样性表现出高单倍型多样性(0.980 5~0.997 1)和低核苷酸多样性(0.048 2~0.060 9)的特点;系统发育分析、分子方差分析和成对的Fst分析显示南海各群体之间有较高的遗传同质性,遗传分化不显著,也没有出现明显的地理分支或聚簇;中性检测和不对称分布分析发现南海的短尾大眼鲷群体发生过种群扩张。结果表明,南海海域的6个短尾大眼鲷群体属于同一种群,该结果为今后中国与南海周边国家渔业资源共享和争取捕捞配额提供了依据。

    Abstract:

    We investigated the genetic structure of red bigeye (Priacanthus macracanthus) population in six sampling sites off the South China Sea including five coastal waters of Beibu Gulf, Sanya, Maoming, Zhuhai, Shantou and the Southwest Nansha. Sequences of the mitochondrial control region (D-loop) were amplified by polymerase chain reaction. Fragments of 729 bp of the mitochondrial D-loop were sequenced from 224 individuals, and 101 mutated sites, 94 nucleotide polymorphism sites as well as 172 haplotypes were identified. All the six groups were characterized with high haplotype diversity (0.980 5~ 0.997 1) and low nucleotide diversity (0.048 2~0.060 9). Analyses of phylogeny, molecular variances and pairwise Fst showed no significant genealogical clades or clusters of samples corresponding to sampling localities. Both neutrality tests and mismatch distribution analyses revealed a recent population expansion in P.macracanthus. Thus, the six P.macracanthus groups in the South China Sea belong to one population. The results provide references for the fishery resources sharing and quotas winning between China and the South China Sea surrounding countries.

  • 罗氏沼虾(Macrobrachium rosenbergii)又名马来西亚大虾、淡水长臂大虾,隶属节肢动物门、甲壳纲、十足目、长臂虾科、沼虾属。原产于印度和东南亚,是重要的淡水经济虾类[1]。20世纪70年代,中国水产科学研究院珠江水产研究所从日本、马来西亚引进罗氏沼虾,先后采用“浓缩海水”与“人工海水”配方,人工繁殖虾苗获得成功,并向全国多个省市推广养殖。现在中国罗氏沼虾年产量为13.5×104 t,产值达30×109[2]。近年在罗氏沼虾养殖中出现了生长速度减慢、性成熟时间提早、抗病力下降等问题,严重影响罗氏沼虾的养殖效益,也危及产品的质量安全。因此了解现有罗氏沼虾养殖种群的遗传多样性可为培育罗氏沼虾优良品种奠定基础。

    近十年来针对罗氏沼虾种质遗传变异已有不少的研究报道。张海琪等[3]采用同工酶和RAPD分析缅甸罗氏沼虾野生群体和浙江养殖群体的生化遗传变异;蒋钦杨等[4]用线粒体16S rRNA基因部分片段和RFLP比较缅甸引进种子代、广西养殖群体及江苏养殖群体的遗传结构;杨学明等[5]用线粒体DNA的COI基因比较了缅甸原种F1代、江苏养殖群体和广西选育群体的遗传多样性;姚茜等[6]用线粒体控制区(D-loop)基因序列比较浙江湖州养殖群体和“南太湖1号”的遗传多样性;陈雪峰等[7]用核糖体转录间隔区2基因(ITS2)序列进行孟加拉和缅甸野生群体、浙江和广西养殖群体、以及“南太湖2号”选育群体的遗传多样性的比较。

    应用微卫星标记进行罗氏沼虾遗传多样性的研究报道尚不多见,仅见CHAREONTAWEE等[8]用6个微卫星标记对泰国5个养殖群体和2个野生群体的遗传多样性进行了分析;朱其建等[9]用8个微卫星标记对16个来自于上海养殖群体的选育群体进行遗传多样性分析。广东与江苏是罗氏沼虾的两大主养区,两省的罗氏沼虾养殖产量占全国总养殖产量的近80%,同时也均为罗氏沼虾苗种主要培育基地[10]。由于罗氏沼虾引进中国的时间较长,虽然也有多次小规模的重新引种,但多数虾苗场尚缺乏科学的保种与选育技术,因此笔者研究选取广东、江苏以及泰国的罗氏沼虾养殖群体,采用微卫星标记技术进行种群遗传结构和遗传多样性的分析,旨在为罗氏沼虾的选育种工作提供基础数据。

    试验样品分别采自6个养殖场的普通养殖群体,均属非选育群体。其中在中国境内共有4个采样点,广东养殖群体1(GD1,32尾),广东养殖群体2(GD2,30尾),广东养殖群体3(GD3,32尾),江苏养殖群体(JS,34尾)。在泰国境内有2个采样点,泰国养殖群体1(TG1,29尾),泰国养殖群体2(TG2,25尾)。6个群体共计180尾。每个群体均为随机取样。分别剪取虾尾部肌肉组织,保存于95%的乙醇中。

    根据GenBank上登录的罗氏沼虾微卫星序列[11-13],选取、设计合成30对引物并进行预扩增,从中挑选12对特异性强、重复性好和多态性较高的微卫星引物(表 1),由生工生物工程(上海)股份有限公司合成。各对引物分别带有FAM、TAMRA和HEX的特异性荧光标记用于测算扩增目的条带的大小。

    表  1  用于罗氏沼虾微卫星扩增的引物序列及其退火温度
    Table  1.  Primer sequence and annealing temperature for microsatellite marker amplification
    序号No. 位点loci 重复序列repeat sequence 引物序列primer sequence 片段大小production size 复性温度/℃ annealing
    1 EMR-31B (CA)7 F: FAM-GCTGTGCTCCAAAATCTCTCTC
    R: CTCACCCATACTTGACAACGAC
    210 58
    2 EMR-55 (TTA)13 F: FAM-GAAGTCATCCGACAAACTTCAC
    R: AGTAATCATGTGGCCTAGCCTAG
    200 58
    3 EMR81 (CT)15 F: TAMRA-GGAACCAGTGAAAAAGCAATGA
    R: GGGGTGCATTCAAAAATAGGT
    240 58
    4 EMR85 (AG)12 F: FAM-GACGGACAGACATTCATTAGCC
    R: ATTCACCCCACACTTTGACATT
    220 58
    5 Mbr1 (GA)24 F: HEX-CCCACCATCAATTCTCACTTACC
    R: TCCTTTTCACATCGTTTCCAGTC
    272~320 60
    6 Mbr-2 (GT)22 F: HEX-TTCCCGACCAATTTCTCTTTCTC
    R: GGCAAAAATGATCTTGGATTCAC
    298~336 60
    7 Mbr-3 (AG)14 F: TAMRA-CAACTCTATGTTTCGGCATTTGG
    R: GGGGAATTTTACCGATGTTTCTG
    232~284 62
    8 Mbr-4 (GT)29 F: HEX-CCACCTACCGTACATTCCCAAAC
    R: CGGGGCGACTTTTAGTATCGAC
    288~326 62
    9 Mbr-5 (AG)25 F: HEX-CAAGGCTCGTGTCTCTTGTTTC
    R: GCTTGTACTTGTTCAGCTTTTGC
    286~328 62
    10 Mbr-9 (TG)5(AG)17 F: TAMRA-TTGTTTGCTTGTTTAGTGTCAAGG
    R: CTCCAAAACCGAAAAATCCTCAC
    240~284 60
    11 Mbr11 (AG)31 F: TAMRA-GTATTGAGAACAAAGGCGCACAG
    R: ATCTCTTTCCAAAACAGGGCACA
    263~291 60
    12 MRMB25 (CT)17 F: FAM-CCGGTCCAAAGGAATACAGA
    R: GTGGGTGGTGTCCCTATGTT
    194~216 46
    下载: 导出CSV 
    | 显示表格

    按照天根生化科技(北京)有限公司的TIANamp Genomic DNA Kit试剂盒说明书提取罗氏沼虾肌肉组织的总DNA,并用0.8%的琼脂糖凝胶电泳和分光光度计检测DNA质量和浓度,-20 ℃保存备用。

    引物分为3类,分别用3种不同的荧光进行标记(FAM、TAMRA和HEX)(表 1)。PCR反应总体积为20 μL,含有10×buffer 2.0 μL、氯化镁(MgCl2)(25 mmol·L-1)0.8 μL、dNTP(2.5 mmol·L-1)0.4 μL、上下游引物(20 μmol·L-1)各0.4 μL、基因组DNA 40 ng、Taq酶[天根生化科技(北京)有限公司出品]1 U。PCR扩增程序为94 ℃预变性5 min后进入35个循环,94 ℃,30 s;55 ℃退火,30 s;72 ℃,30 s,循环结束后72 ℃再延伸10 min。

    STR基因分型委托生工生物工程(上海)股份有限公司完成。使用DYY-8型稳压稳流电泳仪(上海琪特分析仪器有限公司出品)、凝胶成像系统(Gene Genius公司出品)和3730XL测序分析仪(美国ABI公司出品)进行STR序列分析,根据每个扩增条带分子量的差异性,判断每个个体各基因座的基因型。

    通过STR分型确定每个标记在每个个体中的基因型,根据条带(等位基因)的大小,按从大到小按字母排序,用POPGENE 32进行统计分析,计算各群体每一个微卫星位点的等位基因数(Na)、有效等位基因数(Ne)、Shannon指数(I)、观测杂合度(Ho)、期望杂合度(He)、群体间遗传距离(Da)、遗传相似度(S)。FSTAT软件计算遗传分化指数(Fst)。根据Nei氏遗传距离通过MEGA 5软件利用非加权配对算术平均法(unweighted pair group method with arithmetic mean,UPGMA)构建系统树[14]。用CERVUS 3.0软件计算多态信息含量(polymorphism information content,PIC)[15]。群体间Fst、群体分子方差分析(AMOVA)采用ARLEQUIN 3.1软件进行[16]。利用STRUCTURE 2.3软件分析群体的遗传结构,再用最大似然性法分析最佳K值,即理论群体数[17]

    根据预扩增结果从30个微卫星位点中选取特异性强、重复性好和多态性较高的12对引物,对6个罗氏沼虾群体共180个个体的基因组DNA进行扩增,结果显示,各对引物在6个群体中均能扩增出目的条带,且表现出不同程度的多态性(表 2)。其中Na为13~28个,Ne为6.61~18.33,I为2.124~3.045,Ho为0.488~0.994,He为0.852~0.948;PIC为0.853~0.941,所检测的12个位点均属于高度多态位点(PIC>0.5)。

    表  2  罗氏沼虾12个微卫星位点的等位基因数、杂合度及多态信息含量
    Table  2.  Number of alleles, heterozygosity and polymorphic information content of 12 microsatellite loci of M.rosenbergii
    位点locus 等位基因数Na 有效等位基因数Ne Shannon指数I 观测杂合度Ho 期望杂合度He 多态信息含量PIC
    EMR-31B 22 12.71 2.766 0.673 0.924 0.915
    EMR-55 14 10.79 2.506 0.906 0.910 0.908
    EMR81 23 13.72 2.853 0.980 0.930 0.919
    EMR85 21 14.70 2.840 0.853 0.935 0.933
    Mbr1 19 12.96 2.725 0.671 0.926 0.897
    Mbr-2 13 6.61 2.124 0.691 0.852 0.853
    Mbr-3 16 8.68 2.449 0.590 0.888 0.887
    Mbr-4 24 18.33 3.024 0.723 0.948 0.941
    Mbr-5 19 11.40 2.666 0.891 0.915 0.914
    Mbr-9 28 16.57 3.045 0.713 0.943 0.939
    Mbr11 17 8.41 2.439 0.488 0.884 0.886
    MRMB25 14 7.88 2.321 0.994 0.876 0.876
    均值mean 19.17 11.90 2.646 0.765 0.911 0.906
    下载: 导出CSV 
    | 显示表格

    6个罗氏沼虾群体的遗传多样性见表 3。GD3群体的平均Na、平均NeI均为最高(Na=13.833,Ne=8.947,I=2.348),而TG2的平均Na、平均NeI均为最低(Na=10.917,Ne=6.558,I=2.046)。Ho以GD1的最高(Ho=0.821)、TG2的最低(Ho=0.683)。He和多态信息含量则以GD3的最高(He=0.896,PIC=0.880)、TG2的最低(He=0.848,PIC=0.806)。综合上述各参数,可认为所检测的6个群体均具有较高的遗传多样性,并以广东群体3(GD3)的遗传多样性最高、泰国群体TG2的遗传多样性最低。

    表  3  6个罗氏沼虾养殖群体的遗传多样性
    Table  3.  Genetic diversity of six populations of M.rosenbergii
    参数index 广东1 GD1 江苏JS 广东2 GD2 广东3 GD3 泰国1 TG1 泰国2 TG2
    等位基因数Na 13.000 12.750 12.000 13.833 12.500 10.917
    有效等位基因数Ne 7.409 7.204 6.791 8.947 8.274 6.558
    Shannon指数I 2.207 2.159 2.135 2.348 2.239 2.046
    观测杂合度Ho 0.821 0.751 0.720 0.761 0.808 0.683
    期望杂合度He 0.871 0.858 0.861 0.896 0.881 0.848
    多态信息含量PIC 0.850 0.846 0.845 0.880 0.857 0.806
    下载: 导出CSV 
    | 显示表格

    6个罗氏沼虾养殖群体间的Fst为0.006 1~0.093 1(表 4),属于中低水平遗传分化(Fst<0.15)。其中泰国群体1和泰国群体2之间遗传分化水平最低(Fst=0.006 1),属于低程度的分化(Fst<0.05)。而泰国群体2和广东群体1之间遗传分化水平最高(Fst=0.093 1),两者间的分化为中等程度的分化(0.05<Fst<0.15)(表 4)。

    表  4  6个罗氏沼虾养殖群体的遗传分化指数
    Table  4.  Fixation index (Fst) of six populations of M.rosenbergii
    广东1 GD1 江苏JS 广东2 GD2 广东3 GD3 泰国1 TG1 泰国2 TG2
    广东1 GD1
    江苏JS 0.023 0
    广东2 GD2 0.024 5 0.006 3
    广东3 GD3 0.031 0 0.046 9 0.041 0
    泰国1 TG1 0.071 3 0.073 4 0.074 9 0.058 4
    泰国2 TG2 0.093 1 0.077 8 0.082 5 0.072 9 0.006 1
    下载: 导出CSV 
    | 显示表格

    AMOVA分析结果表明,群体中仅有5.27%的遗传变异来自于群体间,94.73%的遗传变异来自于群体内(P<0.01),表明遗传变异主要存在于个体之间,个体间的遗传变异远大于群体间的遗传变异(表 5)。

    表  5  6个罗氏沼虾养殖群体的AMOVA分析
    Table  5.  AMOVA analysis among six populations of M.rosenbergii
    变异来源source of variation 自由度df 平方和sum of squares 方差组分variance component 百分率/% percentage
    群体间among populations 5 108.518 0.296 76 5.273 51
    群体内within populations 312 1663.064 5.330 52 94.726 49**
    总变异total 317 1771.582 5.627 28
    注:* *. 1 023次模拟检验后显示为极显著(P<0.01)
    Note:* *. very significance (P<0.01) after 1 023 permutation tests
    下载: 导出CSV 
    | 显示表格

    6个群体间Nei氏Da为0.160 8~0.695 8,S为0.311 7~0.851 5(表 6),其中江苏群体和广东群体2的Da最近(0.160 8)、S最高(0.851 5),而广东群体1和泰国群体2间的Da最远(0.695 8)、S最低(0.311 7)。根据遗传距离构建的UPGMA聚类树显示6个罗氏沼虾群体聚为2个大支,3个广东群体与江苏群体聚为一支,2个泰国群体聚为另一支(图 1)。

    表  6  6个罗氏沼虾养殖群体的遗传距离(对角线下)和遗传相似度(对角线上)
    Table  6.  Nei′s genetic distance (below diagonal)and genetic similarity (above diagonal)of six populations of M.rosenbergii
    广东1 GD1 江苏JS 广东2 GD2 广东3 GD3 泰国1 TG1 泰国2 TG2
    广东1 GD1 - 0.760 4 0.737 4 0.664 0 0.393 9 0.311 7
    江苏JS 0.273 9 - 0.851 5 0.572 2 0.407 7 0.432 5
    广东2 GD2 0.304 6 0.160 8 - 0.597 4 0.382 2 0.390 0
    广东3 GD3 0.409 4 0.458 2 0.415 2 - 0.423 5 0.379 7
    泰国1 TG1 0.611 7 0.597 1 0.661 8 0.559 3 - 0.814 2
    泰国2 TG2 0.695 8 0.568 1 0.611 7 0.628 4 0.205 5 -
    下载: 导出CSV 
    | 显示表格
    图  1  基于Nei′s遗传距离构建的6个罗氏沼虾群体的UPGMA聚类树
    Figure  1.  UPGMA dendrograrn of six populations of M.rosenbergii based on Nei′s genetic distance

    利用STRUCTURE软件执行了2-7的假设K,设10次重复,根据K对应的参数、采用基于群体结构的贝叶斯分析法发现,K=4时出现平台期,推断笔者研究的所有参试个体最佳分组为4个理论群(图 2)。广东群体1、广东群体2和江苏群体群体聚为一组,广东群体3和2个泰国群体则各自成为一组。

    图  2  参试罗氏沼虾群体在K=4时的遗传结构图
    Figure  2.  Cluster analysis based on 14 microsatellite loci for M.rosenbergii specimens from STRUCTURE (K=4)

    微卫星标记是分布于真核生物基因组中的简单重复序列(simple sequence repeats,SSRs),也称短串联重复序列(simple tandem repeats,STRs)。因其在基因组中分布广泛、具较高的多态性、呈共显性遗传,且由于分析方法相对简便快捷也被广泛应用于水产生物种群遗传多样性分析、种质资源保护、遗传图谱建立和QTL定位等研究中[18-20]。微卫星标记分析通常是通过PCR扩增、电泳检测和片段大小分离分析各等位基因。传统的片段分离采用聚丙烯酰胺凝胶电泳结合放射显影或银染的方法,其分辨率和效率均较低。笔者研究利用近年发展起来的STR分型新方法,通过荧光标记的PCR引物扩增微卫星位点区域序列,利用ABI遗传分析仪对荧光标记的DNA片段进行毛细管电泳检测,结合分子量内标计算DNA片段长度,使STR分型更高效与准确[21]

    多态信息含量(PIC)起衡量基因位点多态性的作用,通常PIC能反映出某个群体的遗传变异程度、位点多样性等[22]。笔者研究所用的12个微卫星位点的PIC介于0.853~0.941,根据BOTSTEIN等提出的划分标准,PIC>0.5时该等位基因座位为高度多态性座位[22],所采用的12个位点均属于高度多态(PIC>0.5),可用于罗氏沼虾群体的遗传多样性及遗传结构的评估。

    杂合度(heterozygosity)是评价群体遗传变异的重要指标。Ho容易受样本大小的影响,而He更能够反映群体的遗传多样性,杂合度越高物种遗传多样性越丰富,对环境适应能力则越强[24-25]。CHAREONTAWEE等[8]用6个微卫星标记分析了泰国的5个养殖群体和2个野生群体,认为这些养殖群体仍具有较高的遗传多样性水平,其中2个养殖群体的遗传多样性(He均为0.73)甚至高于2个野生群体(He=0.71和0.68)。朱其建等[9]用8个微卫星标记分析采自上海地区的16个罗氏沼虾选育群体,结果显示各群体的He相差较大(0.551 9~0.733 2)。在应用微卫星标记进行生物遗传多样性分析时其结果可能受选用的标记在基因组中所处的位置、位点的多态性以及分型手段等的影响。相比于STR基因分型技术,聚丙烯酰胺凝胶电泳的分辨率相对较低,可能会造成等位基因缺失、杂合度偏低[26]。笔者研究采用STR基因分型,结果显示6个罗氏沼虾养殖群体的平均He介于0.848~0.896,说明这6个群体均具有较高的遗传多样性与选育潜力。

    遗传距离和遗传相似度是衡量群体间遗传关系的指标[27]。笔者研究结果发现,江苏群体和广东群体2的Da最近(0.160 8)、S最高(0.851 5),说明其亲缘关系最近;而广东群体1和泰国群体2之间的遗传距离最远(0.695 8)、S最低(0.311 7),说明其亲缘关系最远,且广东群体1和泰国群体2间的遗传分化水平也最高(Fst=0.093 1)。同时系统进化树更直观地表现出这种关系,中国的4个群体聚为一支,泰国的2个群体聚为另一支,说明泰国罗氏沼虾养殖群体和中国养殖群体的亲缘关系相对较远。

    STRUCTURE软件是基于个体的遗传组成进行群体模拟分析,不受各群体样本数的影响,是群体遗传结构分析的理想工具[26, 28]。笔者研究的所有参试个体被划分为4个地理群,支持UPGMA树的聚类结果。同时遗传结构图的直观显示结果与遗传分化分析结果相吻合。广东群体1、广东群体2和江苏群体群体聚为一组,说明广东群体1、广东群体2和江苏群体之间的遗传结构相似,3个群体之间存在着一定程度的基因交换,这可能与江苏和广东两个罗氏沼虾养殖大省存在着一定程度的亲虾或苗种互相供给与交流有关。

    文章所分析的罗氏沼虾养殖群体均具有较高的遗传多样性;中国养殖群体间遗传分化水平较低,但中国与泰国群体间的遗传分化程度达到中等水平,且遗传距离大、亲缘关系远,在下一步的罗氏沼虾选育种中可考虑引入泰国群体以提高选育群体的遗传多样性。

  • 图  1   短尾大眼鲷标本采集地点图

    Figure  1.   Sample sites for Priacanthus macracanthus

    图  2   短尾大眼鲷线粒体D-loop序列的NJ发育树

    Figure  2.   Neighbour-joining tree for mitochondrial D-loop sequences of P.macracanthus

    图  3   短尾大眼鲷线粒体D-loop序列单倍型的MJ网络图

    Figure  3.   Median-joining network for mitochondrial D-loop haplotypes of P.macracanthus

    表  1   不同群体短尾大眼鲷的碱基频率、中性检验和不对称分布参数

    Table  1   Statistics of Tajima′s D, Fu′s Fs and mismatch distribution parameter of P.macracanthus

    群体group 编号ID T/% C/% A/% G/% Tajima′s D Fu′s Fs mismatch distribution
    D P Fs P τ θ0 θ1
    南沙Nansha NS 34.2 19.1 31.8 15.0 -1.798 0 0.014 2 -25.346 8 0.000 0 3.675 7 1.453 7 99 999.000 0
    北部湾Beibu Gulf BG 34.1 19.1 31.8 15.0 -1.696 4 0.022 7 -25.345 2 0.000 0 5.718 7 0.000 0 99 999.000 0
    三亚Sanya SY 34.1 19.1 31.8 15.0 -1.801 9 0.015 0 -25.642 2 0.000 0 4.500 0 0.000 0 99 999.000 0
    茂名Maoming MM 34.1 19.1 31.8 15.0 -1.641 3 0.031 0 -25.779 2 0.000 0 4.343 7 0.001 7 99 999.000 0
    珠海Zhuhai ZH 34.1 19.1 31.8 15.0 -1.730 1 0.020 9 -25.597 2 0.000 0 4.679 6 0.014 0 143.359 3
    汕头Shantou ST 34.1 19.1 31.7 15.0 -1.704 2 0.024 5 -25.104 1 0.000 0 4.187 5 1.221 4 260.625 0
    平均average 34.1 19.1 31.8 15.0 -2.062 6 0.000 3 -25.101 2 0.000 9 4.931 6 0.008 7 99 999.000 0
    下载: 导出CSV

    表  2   短尾大眼鲷样本信息及遗传多样性参数

    Table  2   Sampling information and descriptive statistics of genetic diversity of P.macracanthus

    编号ID 群体group 经度longitude 纬度latitude 采样数量quantity of sample 单倍型数(私有单倍型数) No.of haplotypes (private) 多态位点(S) segregating site 单倍型多样性(h) haplotype diversity 核苷酸多样性(π) nucleotide diverity 平均核苷酸差异(k) mean number of nucleotide differences
    NS 南沙Nansha 109.650E 6.037N 40 37(28) 49 0.994 9±0.007 5 0.060 6±0.032 9 5.698 7±2.788 8
    BG 北部湾Beibu Gulf 107.483E 20.427N 46 42(37) 49 0.997 1±0.005 7 0.060 9±0.032 9 5.725 6±2.792 8
    SY 三亚Sanya 110.135E 18.075N 37 30(26) 39 0.989 5±0.009 2 0.048 9±0.027 3 4.603 6±2.311 9
    MM 茂名Maoming 111.548E 21.132N 34 31(23) 32 0.992 9±0.009 9 0.044 7±0.025 3 4.210 3±2.142 9
    ZH 珠海Zhuhai 113.774E 22.353N 37 29(23) 37 0.980 5±0.013 5 0.048 2±0.026 9 4.531 5±2.280 2
    ST 汕头Shantou 116.932E 23.062N 30 28(21) 42 0.993 1±0.011 8 0.060 7±0.033 3 5.710 3±2.813 9
    平均average 224 172(158) 94 0.990 3±0.002 3 0.094 2±0.048 1 9.140 7±4.218 3
    下载: 导出CSV

    表  3   短尾大眼鲷的遗传分化系数Fst(对角线下)和P值(对角线上)

    Table  3   Pairwise Fst (below diagonal) and P value of population differentiation (above diagonal) of P.macracanthus

    南沙NS 北部湾BG 三亚SY 茂名MM 珠海ZH 汕头ST
    南沙NS 0.040 0 0.135 7 0.036 1 0.324 2 0.958 9
    北部湾BG 0.013 2 0.607 4 0.669 9 0.652 3 0.044 9
    三亚SY 0.008 1 -0.002 9 0.793 9 0.976 5 0.338 8
    茂名MM 0.017 1 -0.003 4 -0.006 9 0.832 0 0.015 6
    珠海ZH 0.002 4 -0.003 7 -0.011 4 -0.008 0 0.238 2
    汕头ST -0.012 6 0.014 6 0.001 8 0.024 0 0.004 9
    下载: 导出CSV

    表  4   不同群体短尾大眼鲷遗传差异的AMOVA多样性分析

    Table  4   Analysis of molecular variance (AMOVA) among different groups of P.macracanthus

    变异来源source of variation 自由度degrees of freedom 平方和sum of squares 变异组成variance components 变异百分比percentage of variation F P
    群体间among populations 5 14.087 0.007 1 0.28
    群体内within populations 218 556.654 2.553 5 99.72
    总计total 223 570.741 2.560 6 100 Fst=0.002 8 0.207 9
    下载: 导出CSV
  • [1]

    SOKOLOVSKAYA T G, SOKOLOVSKII A S, SOBOLEVSKII E I. A list of fishes of Peter the Great Bay (the Sea of Japan)[J]. J Ichthyol, 1998, 38(1): 1-11. https://www.semanticscholar.org/paper/A-LIST-OF-FISHES-OF-PETER-THE-GREAT-BAY-(THE-SEA-OF-Sokolovskaya-Sokolovskii/c65b189aad74d9abcdf1cd8052fe67600840f0b1

    [2] 卢占晖, 朱增军, 林显鹏. 东海区短尾大眼鲷食性分析[J]. 海洋渔业, 2010, 32(2): 199-203. doi: 10.3969/j.issn.1004-2490.2010.02.014
    [3] 张仁斋, 陆穗芬. 短尾大眼鲷卵及仔、稚鱼的发育[J]. 水产学报, 1982(3): 243-251. http://scxuebao.ijournals.cn/scxuebao/article/abstract/19820306
    [4] 何筱洁. 大眼鲷匹里虫病及其病原体的研究——包括匹里虫属一新种的描述[J]. 水产学报, 1982(2): 97-106. http://scxuebao.ijournals.cn/scxuebao/article/abstract/19820202
    [5] 孙典荣, 邱永松. 北部湾长尾大眼鲷生长和死亡参数估计[J]. 海洋湖沼通报, 2004(3): 27-34. doi: 10.3969/j.issn.1003-6482.2004.03.005
    [6]

    LIU K M, CHENG Y L. Virtual population analysis of the big eye Priacanthus macracanthus in the waters off northeastern Taiwan[J]. Fish Res, 1999, 41(3): 243-254. doi: 10.1016/S0165-7836(99)00025-9

    [7]

    LIU K M, HUNG K Y, CHEN C T. Reproductive biology of the big eye Priacanthus macracanthus in the north-eastern waters off Taiwan[J]. Fish Sci, 2001, 67(6): 835-838. doi: 10.1046/j.1444-2906.2001.00355.x

    [8]

    BENJAKUL S, OUNGBHO K, VISESSANGUAN W, et al. Characteristics of gelatin from the skins of bigeye snapper, Priacanthus tayenus and Priacanthus macracanthus[J]. Food Chem, 2009, 116(2): 307-325. doi: 10.1016/j.foodchem.2009.02.063

    [9] 孙冬芳, 李永振, 董丽娜. AFLP技术在鱼类种群分析中的应用[J]. 南方水产, 2009, 5(6): 71-75. doi: 10.3969/j.issn.1673-2227.2009.06.013
    [10] 宁平, 俞存根, 虞聪达, 等. 浙江南部外海渔场鱼类区系特征研究[J]. 浙江海洋学院学报, 2008, 27(3): 267-270. https://www.zhangqiaokeyan.com/academic-journal-cn_journal-zhejiang-ocean-university-natural-science_thesis/0201248287722.html
    [11]

    KOICHIRO T, GLEN S, DANIEL P, et. al. MEGA6: Molecular evolutionary genetics analysis version 6.0[J]. Mol Biol Evol, 2013, 30(12): 2725-2729. doi: 10.1093/molbev/mst197

    [12]

    LIBRADO P, ROZAS J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data[J]. Bioinformatics, 2009, 25(11): 1451-1452. doi: 10.1093/bioinformatics/btp187

    [13]

    EXCOFFIER L, LAVAL G, SCHNEIDER S. Arlequin ver. 3.0: an integrated software package for population genetics data analysis[J]. Evol Bioinformatics, 2005, 1: 47-50. https://www.researchgate.net/publication/24238444_ARLEQUIN_Ver_30_An_integrated_software_package_for_population_genetic_data_analysis

    [14]

    NEI M. Molecular evolutionary genetics[M]. New York: Columbia University Press, 1987: 153-187. doi: 10.7312/nei-92038

    [15]

    TAJIMA F. Evolutionary relationship of DNA sequences in finite populations[J]. Genetics, 1983, 105(2): 437-467. doi: 10.1093/genetics/105.2.437

    [16]

    WRIGHT S. The genetical structure of populations[J]. Annals of Eugenics, 1951, 15(1): 323-354. doi: 10.1111/j.1469-1809.1949.tb02451.x

    [17]

    WRIGHT S. Evolution and the genetics of populations[M]//Variability within and among natural populations Vol. 4, Chicago: University of Chicago Press, 1978: 179-182. https://www.semanticscholar.org/paper/Evolution-and-the-Genetics-of-Populations%2C-Volume-Curnow-Wright/7ffa4dd7ecfb20e985c073866956be97764d86a1

    [18]

    FREELAND J R. Molecular ecology[M]. Chichester: John Wiley & Sons Ltd, 2005: 112-116. https://onlinelibrary.wiley.com/toc/1365294x/2005/14/1

    [19]

    BANDELT H J, FORSTER P, ROHL A. Median joining networks for inferring intraspecific phylogenics[J]. Mol Biol Evol, 1999, 16(1) : 37-48. doi: 10.1093/oxfordjournals.molbev.a026036

    [20]

    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. doi: 10.1093/oxfordjournals.molbev.a040727

    [21]

    LI W H. Molecular evolution[M]. Sunderland, Massachusetts, USA: Sinauer Associates, Inc. Publishers, 1997: 59-63.

    [22]

    HARPENDING H C. Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution[J]. Hum Biol, 1994, 66(4): 591-600. https://pubmed.ncbi.nlm.nih.gov/8088750/

    [23]

    LEE W J, CONROY J, HOWELL W H, et al. Structure and evolution of teleost mitochondrial control region[J]. Mol Evol, 1995, 41(4): 54-66. doi: 10.1007/BF00174041

    [24]

    RAMOS-ONSINS S E, ROZAS J. Statistical properties of new neutrality tests against population growth[J]. Mol Biol Evol, 2002, 19(12): 2092-2100. doi: 10.1093/oxfordjournals.molbev.a004034

    [25]

    TAJIMA F. Statistical-method for testing the neutral mutation hypothesis by DNA polymorphism[J]. Genetics, 1989, 123(3): 585-595. doi: 10.1093/genetics/123.3.585

    [26]

    FU Y X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection[J]. Genetics, 1997, 147(2): 915-925. doi: 10.1093/genetics/147.2.915

    [27] 牛素芳, 苏永全, 王军, 等. 福建近海蓝圆鲹群体遗传结构分析[J]. 厦门大学学报: 自然科学版, 2012, 51(4): 759-766. https://www.cnki.com.cn/Article/CJFDTotal-XDZK201204024.htm
    [28] 牛素芳, 苏永全, 王军, 等. 福建近海竹筴鱼线粒体DNA控制区和细胞色素b遗传多态性[J]. 中国水产科学, 2011, 19(1): 66-74. https://www.doc88.com/p-9734747774803.html
    [29]

    LIU J X, GAO T X, WU S F, et al. Pleistocene isolation in the Northwestern Pacific marginal seas and limited dispersal in a marine fish, Chelon haematocheilus (Temminck & Schlegel, 1845)[J]. Mol Ecol, 2007, 16(2): 275-288. doi: 10.1111/j.1365-294X.2006.03140.x

    [30] 苏天凤, 江世贵, 马之明, 等. 二长棘犁齿鲷线粒体DNA控制区结构和进化[J]. 南方水产科学, 2012, 8(6): 1-8. doi: 10.3969/j.issn.2095-0780.2012.06.001
    [31] 罗旭光, 齐景伟, 安晓萍, 等. 内蒙古达里湖及岗更湖东北雅罗鱼D-loop基因序列及多样性研究[J]. 南方水产科学, 2014, 10(2): 92-96. doi: 10.3969/j.issn.2095-0780.2014.02.014
    [32]

    TUDELA S, GARCÍA-MARÍN J L, PLA C. Genetic structure of the European anchovy, Engraulis encrasicolus, in the north-west Mediterranean [J]. J Exp Mar Biol Ecol, 1999, 234(1): 95-109. doi: 10.1016/S0022-0981(98)00142-7

    [33] 张东亚, 汪登强, 刘绍平, 等. 怒江濒危鱼类缺须盆唇鱼基于线粒体Cyt b序列的群体遗传结构分析[J]. 中国水产科学, 2009, 17(4): 477-486. doi: 10.3321/j.issn:1005-8737.2009.04.002
    [34] 刘进贤, 高天翔, 吴世芳, 等. 梭鱼的分子系统地理学研究——晚更新世西北太平洋边缘海隔离分化及其有限的扩散能力[J]. 中国海洋大学学报: 自然科学版, 2007, 37(6): 931-938. doi: 10.16441/j.cnki.hdxb.2007.06.014
    [35]

    AVISE J C. Phylogeography: the history and formation of species[M]. Cambridge, Massachusets: Harvard Unversity Press, 2000. https://www.nature.com/articles/6887654

    [36]

    HEWITT G M. The genetic legacy of the Quaternary ice age[J]. Narture, 2000, 405(6789): 907-913. https://www.nature.com/articles/35016000

    [37]

    DYNESIUS M, JANSSON R. Evolutionary consequences of changes in species′ geographical distribution driven by Milankovitch climate oscillations[J]. Proc Natl Acad Sci USA, 2000, 97(16): 9115-9120. doi: 10.1073/pnas.97.16.9115

    [38]

    LAMBECK K, HONZA E. Global tectonics, and formation of marginal basins: role of the westerm Pacific[J]. Episodes, 1991, 14: 224-230. doi: 10.18814/epiiugs/1991/v14i3/005

    [39] 杨海军, 刘秦玉. 南海海洋环流研究综述[J]. 地球科学进展, 1998(4): 47-51. doi: 10.3321/j.issn:1001-8166.1998.04.007
    [40] 苏纪兰. 南海环流动力机制研究综述[J]. 海洋学报: 中文版, 2005, 27(6): 3-10. http://www.hyxb.org.cn/aos/ch/reader/view_abstract.aspx?file_no=20050601&st=alljournals
    [41] 管秉贤, 袁耀初. 中国近海及其附近海域若干涡旋研究综述Ⅰ. 南海和台湾以东海域[J]. 海洋学报: 中文版, 2006, 28(3): 1-16. https://qikan.cqvip.com/Qikan/Article/Detail?id=21970714
    [42] 孙冬芳, 董丽娜, 李永振, 等. 南海北部海域多齿蛇鲻的种群分析[J]. 水产学报, 2010, 34(9): 1387-1394. doi: 10.3724/SP.J.1231.2010.06975
    [43] 孙冬芳, 卢伟华, 艾红, 等. 中沙和南沙群岛密点胡椒鲷的种群分析[J]. 广东农业科学, 2010(5): 4-7. doi: 10.3969/j.issn.1004-874X.2010.05.002
    [44] 董丽娜. 南海北部陆架区金线鱼(Nemipterus virgatus)种群分子分析[D]. 上海: 上海海洋大学, 2011. https://cdmd.cnki.com.cn/Article/CDMD-10264-1011304986.htm
    [45] 李玉芳, 李永振, 陈国宝, 等. 基于线粒体Cyt b基因的黄鳍马面鲀种群分析[J]. 水产学报, 2014, 38(1): 23-32. doi: 10.3724/SP.J.1231.2014.48864
  • 期刊类型引用(30)

    1. 唐峰华,吴祖立,杨嘉樑,黄洪亮,蒋科技,ESTEFANIA Gonzalez Maiez Violeta,ARREOLA Juan Pedro Vela. 加利福尼亚湾茎柔鱼(Dosidicus gigas)渔场声学调查方法. 中国农业科技导报. 2022(04): 185-192 . 百度学术
    2. 王炜祺,童剑锋,薛铭华. 商业渔船渔业声学数据采集及应用研究进展. 电声技术. 2022(12): 50-53 . 百度学术
    3. 李哲,朱文斌,陈峰,李德伟,刘连为,张洪亮. 近年我国渔业资源声学评估研究进展. 浙江海洋大学学报(自然科学版). 2021(01): 80-85+92 . 百度学术
    4. 王新元,李建龙,章翔,陈孟玲,吴沛霖,马文刚,彭思颖,张玺,谢建军,高菲,许强,王爱民. 三亚蜈支洲岛热带海洋牧场渔业资源现状及季节变动. 海洋与湖沼. 2021(06): 1557-1566 . 百度学术
    5. 黄卉,魏涯,杨贤庆,岑剑伟,赵永强,陈胜军,胡晓,王悦齐,郝淑贤. 肌肉改性处理对鸢乌贼品质的影响. 食品与发酵工业. 2020(22): 42-47 . 百度学术
    6. 范江涛,张俊,冯雪,陈作志. 南沙海域鸢乌贼渔场与海洋环境因子的关系. 上海海洋大学学报. 2019(03): 419-426 . 百度学术
    7. 招春旭,邱星宇,何雄波,王国易,顾银娜,钟亚娜,颜云榕,康斌. 南海春季月相、水深、作业时间与鸢乌贼CPUE的关系. 水产学报. 2019(11): 2372-2382 . 百度学术
    8. 王欢欢,张俊,陈作志,汤勇,孙铭帅,黄佳兴,邱永松. 鸢乌贼目标强度绳系控制法测量. 水产学报. 2019(12): 2533-2544 . 百度学术
    9. 江淼,马胜伟,吴洽儿. 南海鸢乌贼资源探捕与开发. 中国渔业经济. 2018(02): 65-70 . 百度学术
    10. 张俊,邱永松,陈作志,张鹏,张魁,范江涛,陈国宝,蔡研聪,孙铭帅. 南海外海大洋性渔业资源调查评估进展. 南方水产科学. 2018(06): 118-127 . 本站查看
    11. 王东旭,陈国宝,汤勇,李斌,王志超. 大亚湾南部海域渔业资源水声学评估. 安徽农业科学. 2017(06): 95-98+159 . 百度学术
    12. 董淑华,黄卉,李来好,杨贤庆,郝淑贤,荣辉,李春生. pH调节法优化鸢乌贼分离蛋白制备的工艺研究. 食品工业科技. 2017(09): 181-185 . 百度学术
    13. 范江涛,张俊,冯雪,陈作志. 基于地统计学的南沙海域鸢乌贼渔场分析. 生态学杂志. 2017(02): 442-446 . 百度学术
    14. 张俊,江艳娥,陈作志,龚玉艳,陈国宝. 南海中南部中层鱼资源声学积分值及时空分布初探. 中国水产科学. 2017(01): 120-135 . 百度学术
    15. 张鹏,张俊,李渊,张然,林龙山,晏磊,邱永松,孙典荣,陈森. 秋季南海中南部海域的一次灯光罩网探捕调查. 南方水产科学. 2016(02): 67-74 . 本站查看
    16. 张俊,张鹏,陈作志,陈国宝,张魁,许友伟,孙铭帅. 南海外海鲹科鱼类资源量及其分布. 南方水产科学. 2016(04): 38-48 . 本站查看
    17. 粟丽,陈作志,张鹏. 南海中南部海域春秋季鸢乌贼繁殖生物学特征研究. 南方水产科学. 2016(04): 96-102 . 本站查看
    18. 范江涛,陈作志,张俊,冯雪. 基于海洋环境因子和不同权重系数的南海中沙西沙海域鸢乌贼渔场分析. 南方水产科学. 2016(04): 57-63 . 本站查看
    19. 江艳娥,方展强,林昭进,张鹏,陈作志. 南海鸢乌贼耳石微量元素研究. 南方水产科学. 2016(04): 71-79 . 本站查看
    20. 杨贤庆,吴静,胡晓,李来好,吴燕燕,林婉玲,黄卉,裘超颖,马海霞. 响应面法优化酶解鸢乌贼制备抗氧化肽的工艺研究. 食品工业科技. 2016(11): 242-248 . 百度学术
    21. 杨丽芝,杨贤庆,黄卉,李来好,邓建朝,赵永强,杨少玲. 响应面法优化南海鸢乌贼墨汁多糖提取工艺. 食品工业科技. 2016(18): 292-296 . 百度学术
    22. 张俊,陈国宝,陈作志,于杰,范江涛,邱永松. 南沙南部陆架海域渔业资源声学评估. 南方水产科学. 2015(05): 1-10 . 本站查看
    23. 范江涛,张俊,冯雪,陈作志. 基于栖息地模型的南沙海域鸢乌贼渔情预报研究. 南方水产科学. 2015(05): 20-26 . 本站查看
    24. 许友伟,陈作志,范江涛,张俊. 南沙西南陆架海域底拖网渔获物组成及生物多样性. 南方水产科学. 2015(05): 76-81 . 本站查看
    25. 江艳娥,张鹏,林昭进,邱永松,方展强,陈作志. 南海外海鸢乌贼耳石形态特征分析. 南方水产科学. 2015(05): 27-37 . 本站查看
    26. 杨贤庆,杨丽芝,黄卉,李来好,邓建朝,赵永强,杨少玲. 南海鸢乌贼墨汁营养成分分析与评价. 南方水产科学. 2015(05): 138-142 . 本站查看
    27. 黄卉,沈玉,吴燕燕,杨贤庆,胡晓,王锦旭. 南海中南部海域鸢乌贼甲醛本底质量分数调查分析. 南方水产科学. 2015(05): 132-137 . 本站查看
    28. 李杰,晏磊,陈森,杨炳忠,谭永光,张鹏. 灯光罩网网口沉降与闭合性能研究. 南方水产科学. 2015(05): 117-124 . 本站查看
    29. 张俊,陈国宝,江艳娥,龚玉艳,陈作志,邱永松,梁沛文. 南海中南部长钻光鱼渔业生物学初步研究. 热带海洋学报. 2015(03): 62-67 . 百度学术
    30. 张俊,陈丕茂,房立晨,陈国宝,黎小国. 南海柘林湾—南澳岛海洋牧场渔业资源本底声学评估. 水产学报. 2015(08): 1187-1198 . 百度学术

    其他类型引用(11)

图(3)  /  表(4)
计量
  • 文章访问数:  3697
  • HTML全文浏览量:  209
  • PDF下载量:  1787
  • 被引次数: 41
出版历程
  • 收稿日期:  2014-12-09
  • 修回日期:  2015-01-08
  • 刊出日期:  2015-04-04

目录

/

返回文章
返回