Population genetic structure of brushtooth lizardfish (Saurida undosquamis) based on mitochondrial cytochrome b gene sequences
-
摘要:
文章利用线粒体细胞色素b (cytochrome b, cytb)基因全序列作为分子标记,分析了中国近海和陆架的花斑蛇鲻(Saurida undosquamis)的遗传结构特征。从8个采样点266尾样本中共检测到142种单倍型,各采样点均表现出很高的单倍型多样性(0.925 1~0.992 9)和较低的核苷酸多样性(0.003 145~0.003 852)。单倍型的中间连接网络图呈现以4个优势共享单倍型为中心的星状发散结构,未发现与地理群体对应的谱系结构。分子方差分析表明花斑蛇鲻的遗传变异绝大部分(99.79%)来自种群内的个体之间,而群体之间几乎没有贡献遗传变异。成对遗传分化系数(FST)显示花斑蛇鲻群体间基因交流频繁,不存在明显的遗传差异,是一个随机交配群。中性检验表明种群历史动态显著偏离稳定种群模型,核苷酸错配分布表明花斑蛇鲻历史上曾经历过种群的快速扩张,扩张时间推算约在距今4万~10万年之前。研究结果表明,中国近海和陆架的花斑蛇鲻遗传分化不显著,在渔业上可以作为一个单元来管理。
Abstract:We examined the population genetic structure of brushtooth lizardfish (Saurida undosquamis) from the coastal waters and continental shelf of China based on the complete cytochrome b gene sequences. A total of 142 haplotypes have been derived from 266 individuals at eight sampling sites. Genetic diversity analysis shows that the S. undosquamis samples were characterized by rather high haplotype diversity (0.925 1−0.992 9) and relatively low nucleotide diversity (0.003 145−0.003 852). Median joining network shows that the haplotypes were connected in a star-like manner with four dominant haplotypes without clades of samples corresponding to sampling localities. Analyses of molecular variance suggest that nearly all the genetic variation (99.79%) was from individuals within population rather than from the populations. The pairwise FST reveals no genetic differentiation among different populations. A high rate of gene flow between populations implies a panmixia for S. undosquamis among the sampling sites. Neutrality tests show significant deviation of population historical dynamics from stable population model, and mismatch distribution analysis indicates a late Pleistocene expansion in S. undosquamis. The expansion time was estimated to be 40−100 thousands years ago. The results reveal that S. undosquamis from the coastal waters and continental shelf of China belong to the same population and can be managed as a unit in fisheries.
-
游泳能力对水生动物的生存能力有着巨大影响,而从力竭运动后状态恢复的能力反映了水生动物的持续极限运动能力。在野外环境中,力竭运动常出现在动物进行洄游、穿越激流、躲避天敌攻击和捕食等场景中;在人工养殖环境中,力竭运动常出现在动物抢食、逃避捕捞、同类相残、高密度暂养和运输等场景中[1-3]。对虾类具有较为出色的运动能力,其在低速游泳时使用游泳足摆动进行游泳,可适应较长距离的洄游活动;受到攻击或刺激时使用发达的腹部肌肉进行收缩,尾扇张开拨水,迅速向后弹跳逃离危险水域[4-5]。研究报道,鱼类和虾类在力竭运动时,无氧呼吸强度显著加强,无氧呼吸链的酶活性显著升高,有氧代谢活动被抑制[6-8];力竭运动结束后,肌肉中会积累乳酸,高浓度乳酸将会限制动物再次进行力竭运动的强度和持续运动的耐久力[9-12]。根据黑鲷 (Sparus macrocephlus) 幼鱼的血乳酸浓度恢复情况分析,黑鲷幼鱼在力竭运动后4 h即可再次进行力竭运动[13];斑马鱼 (Danio rerio) 在力竭运动后第8小时肌乳酸下降恢复至对照组水平,虽然肌糖原未恢复至对照组水平,但已能进行力竭运动[14],可见动物体内乳酸水平是限制动物重复进行力竭运动的关键因素。
对虾类在运动疲劳时,会出现肝、肌糖原浓度均显著下降、血糖和血乳酸浓度显著升高的现象[15-16],属于“乳酸释放型”动物。对虾类的内环境平衡比鱼类的更脆弱,失衡后恢复所需的时间亦更长[17-18]。凡纳滨对虾 (Litopenaeus vannamei) 和中国明对虾 (Fenneropenaeus chinensis) 虾苗在运动频率显著上升的背景下,出现生长率显著降低、免疫能力低下和死亡率升高的现象[19-20];凡纳滨对虾在弹跳运动中,体内蛋白质代谢供能、无氧和有氧代谢强度均显著提升[21],证明过量运动会对虾类的代谢健康造成强烈影响。斑节对虾 (Penaeus monodon) 隶属软甲纲、十足目、对虾科、对虾属,俗名草虾、黑虎虾等,具有适应性强、个体大、生长速度快等优点,是中国沿海重要的经济养殖和增殖放流品种[22]。在目前虾类增养殖业快速发展的背景下,虽然已经有关于虾类运动前后代谢变化的报道,但虾类力竭运动后的恢复代谢特征与调节机制尚不明晰。因此,本研究以斑节对虾为对象,在实验室水槽模拟斑节对虾幼虾力竭运动情况,探究其在力竭运动后肝胰腺与腹部肌肉中的代谢能力相关指标的恢复调节过程,以期为斑节对虾养殖生产和增殖放流提供科学参考。
1. 材料与方法
1.1 实验材料
实验使用的斑节对虾幼苗为中国水产科学研究院南海水产研究所深圳试验基地自行繁育的养殖虾苗。虾苗体长为0.5~0.9 cm时放入室外池塘中放养标粗3周,标粗后挑选体长为 (3.0±0.4) cm、体质量为 (0.2±0.07) g的约1 000尾规格相近的斑节对虾幼虾进行力竭运动胁迫。在力竭运动后的24 h内禁食,在之后的恢复阶段每日共投喂2.25 g (全部虾苗总质量的3%) 的斑节对虾配合0号饲料,分3次 (8:00、16:00、23:00) 投喂。力竭运动实验装置为自行设计的环形水道 (外径为1 004 mm,内置隔离装置直径为 604 mm,内置隔离装置高为395 mm,容量为200 L)。
1.2 实验方法
1.2.1 力竭状态幼虾组织样品收集
从育肥池中随机挑选30尾幼虾为一组,设3组平行,进行力竭运动胁迫。将幼虾放入环形水道中进行胁迫实验,期间水道中流速控制在12 cm·s−1 (在预实验时测得幼虾在该流速下可顶水游泳1.5 h),用小抄网从幼虾游泳方向前方捕捞攻击幼虾,使幼虾持续进行跳跃逃避运动,根据已有报道,当幼虾对刺激无法进行反应、游泳足摆动频率下降、无法再逆水流游泳的时候即进入力竭状态[21]。捕捞刺激持续约1 h,当大部分幼虾陷入力竭状态时,在其中所有进入力竭状态的幼虾中随机挑选12尾置于冰冷的取样盘中取样,收集其肝胰腺与腹部肌肉进行之后的指标测定,为力竭运动后第0小时的样品。取样后立即置入液氮中急冻,之后于−80 ℃保存。
1.2.2 恢复过程中幼虾组织收集
按恢复时间点 (第1、第4、第8、第12、第24、第48、第72小时) 设立7个处理组,每组设3组平行,每组约30尾。以1.2.1中方法使每组幼虾陷入力竭状态,将每组全部存活的力竭幼虾放入室内恢复池的网箱中,同时开始计算恢复时间,在设立的8个时间点取样,收集幼虾的肝胰腺与肌肉。实验水体信息见表1。
表 1 实验水体的物理化学参数Table 1. Physical and chemical parameters of experimental water body水体
Water body温度
Temperature/℃盐度
SalinitypH 溶解氧
Dissolved oxygen/(mg·L−1)育肥池 Fattening pond 25±1 23±2 8.0±0.3 4.2±0.5 环形水道 Circular waterway 25 23 7.8 4.5 恢复池 Recovering pond 25±1 23 7.8±0.3 5.6±0.5 1.2.3 对照组、空白组样品收集
力竭实验前从育肥池中随机挑选12尾幼虾,收集幼虾肝胰腺与肌肉组织作为空白组样品。在8个恢复时间点收集处理组样品前,在育肥池中随机选择12尾幼虾以收集其肝胰腺与肌肉组织,作为各个取样时间点的对照组样品。
1.2.4 样品生化指标测定
在所有样品采集完成后于实验室 (26 ℃) 测定肝胰腺和肌肉的储能物质含量:总蛋白质 (TP)、乳酸 (LA)、糖原 (GLY);呼吸代谢酶活指标:琥珀酸脱氢酶 (SDH) 活性、乳酸脱氢酶 (LDH) 活性;氨基酸代谢酶活指标:谷草转氨酶 (AST) 活性和谷丙转氨酶 (ALT) 活性。以上指标的测定均使用南京建成生物工程研究所的试剂盒。
1.3 数据处理
实验数据经Excel 2016软件处理后,使用SPSS 21.0软件进行数据分析,以力竭运动后的恢复时间为自变量、以取样时间序列排序的各项生理指标为因变量进行单因素方差分析 (One-Way ANOVA),若组间出现极显著 (P<0.01) 或显著差异 (P<0.05) 时,则用Duncan法进行组间差异显著性的多重比较,数据以“平均值±标准差 (
$\overline X \pm {\rm{SD}} $ )”表示。2. 结果
2.1 力竭运动后不同恢复时间对斑节对虾幼虾代谢能力相关指标的影响
根据表2,斑节对虾幼虾在力竭运动后,随恢复时间的延长,腹部肌肉中GLY、LA和LDH活性存在显著差异 (P<0.05),其中LDH活性差异极显著 (P<0.01)。斑节对虾幼虾在力竭运动后,随恢复时间的延长,肝胰腺中LA差异显著 (P<0.05);SDH、LDH、ALT和AST活性均有极显著差异 (P<0.01)。
表 2 斑节对虾幼虾在力竭运动后8个恢复时间点和空白组腹部肌肉、肝胰腺中代谢指标单因素分析Table 2. Analysis of metabolic capacity in abdominal muscles (M) and heaptopancreas (H) of juvenile P. monodon following strenuous activity by One-Way ANOVA因变量
Dependent variable自由度
df均方
MSF P 肌乳酸质量浓度 LA mass concentraion (M) 8 0.000 19 2.613 0.043 肝乳酸质量浓度 LA mass concentraion (H) 8 0.003 3.082 0.023 肌总蛋白质量分数 TP mass fraction (M) 8 3.152 1.294 0.307 肝总蛋白质量分数 TP mass fraction (H) 8 9.508 12.283 0.000 肌糖原质量分数 GLY mass fraction (M) 8 10.057 2.326 0.045 肝糖原质量分数 GLY mass fraction\ (H) 8 0.582 1.305 0.302 肌SDH活性 SDH activity (M) 8 0.007 0.491 0.847 肝SDH活性 SDH activity (H) 8 21.175 22.601 0.000 肌LDH活性 LDH activity (M) 8 0.232 6.318 0.001 肝LDH活性 LDH activity (H) 8 0.000 29 8.599 0.000 肌ALT活性 ALT activity (M) 8 0.937 2.151 0.085 肝ALT活性 ALT activity (H) 8 0.502 6.380 0.001 肌AST活性 AST activity (M) 8 0.560 0.943 0.507 肝AST活性 AST activity (H) 8 75.490 227.247 0.000 2.2 力竭运动后幼虾的代谢指标随恢复时间延长的变化
2.2.1 TP变化
斑节对虾幼虾肝胰腺TP水平在力竭运动前后无显著差异,随恢复时间延长,肝胰腺中TP水平在第4小时达到顶峰,显著高于对照组 (P<0.05),第4至第12小时间显著下降至低于对照组水平 (P<0.05),第48小时恢复至对照组水平 (图1-A)。腹部肌肉中的TP水平在力竭运动前后和之后的恢复中无显著变化 (图1-B)。
图 1 斑节对虾幼虾力竭运动后肝胰腺 (A) 与肌肉 (B) 中总蛋白质量分数变化不同字母代表各恢复时间点实验组间存在显著差异(P<0.05);“*”表示实验组和对照组在胁迫后同一时间点差异显著(P<0.05),后图同此Figure 1. Total protein mass fractions in hepatopancreas (A) and muscle (B) of juvenile P. monodon following strenuous activityDifferent letters represent significant differences between the experimental groups at different recovery time (P<0.05); *. Significant differences between experimental group and the control group at the same time (P<0.05). The same case in the following figures2.2.2 GLY变化
斑节对虾幼虾肝胰腺GLY水平在整个实验过程中并未出现显著差异(图2-A)。腹部肌肉GLY水平在力竭运动后出现显著下降 (P<0.05),第0至第4小时保持在低于对照组水平。腹部肌肉GLY水平在第1小时后开始恢复,保持上升,第8小时恢复至对照组水平,第72小时显著高于对照组水平 (P<0.05,图2-B)。
2.2.3 LA变化
力竭运动后斑节对虾幼虾肝胰腺中LA水平无显著变化,腹部肌肉LA水平显著上升 (P<0.05)。肝胰腺LA水平随恢复时间延长逐渐上升,第12小时显著高于对照组并达到顶峰 (P<0.05),第12至第24小时内显著回落,第72小时恢复至对照组水平 (图3-A)。腹部肌肉LA水平在第8小时恢复至对照组水平,但在第8至第12小时再次显著升高 (P<0.05),第24小时恢复至对照组水平 (图3-B)。
2.2.4 AST活性变化
斑节对虾幼虾肝胰腺中AST活性力竭运动前后无显著差异,第4至第12小时内显著上升 (P<0.05),第12至第72小时下降至对照组水平(图4-A)。腹部肌肉中的AST活性在力竭运动前后无显著差异,第0至第48小时内维持在对照组水平,但在第48至第72小时显著下降 (P<0.05,图4-B)。
2.2.5 ALT活性变化
斑节对虾幼虾肝胰腺中ALT活性力竭运动前后无显著差异,第4至第12小时内显著上升 (P<0.01),第12至第48小时下降至对照组水平(图5-A)。腹部肌肉中ALT活性在力竭运动后持续下降,第12小时时显著低于对照组水平 (P<0.05),在之后的恢复过程中保持较低水平,第72小时时仍显著低于对照组水平 (图5-B)。
2.2.6 LDH活性变化
斑节对虾幼虾肝胰腺中LDH活性在力竭运动前后无显著差异,第2至第12小时显著上升 (P<0.05),第12至第72小时降低至对照组水平 (P<0.05,图6-A)。腹部肌肉的LDH活性在力竭运动后出现显著升高 (P<0.05),第0至第1小时短暂下降后在第1至第12小时内继续显著升高,并在第12至第72小时内仍保持在较高水平 (P<0.05,图6-B)。
2.2.7 SDH活性变化
斑节对虾幼虾肝胰腺中SDH活性在力竭运动后升高,但无显著差异,第4小时显著低于第0小时水平 (P<0.05),第4至第12小时间显著升高,第24小时恢复至对照组水平(图7-A)。虾腹部肌肉力竭运动后的SDH活性显著降低,低于对照组水平,第12小时恢复至对照组水平,之后再次下降,第72小时仍显著低于对照组水平 (图7-B)。
3. 讨论
3.1 斑节对虾幼虾力竭运动后体内乳酸和糖原的变化
肌肉对乳酸的耐受性与糖原储备量被认为是限制动物进行力竭运动的关键因素,而乳酸代谢能力与糖原恢复能力则是限制动物重复力竭运动的关键因素[1]。具有出色游泳能力的虹鳟 (Oncorhynchus mykiss) 恢复时间短于瓦氏黄颡鱼(Pelteobagrus vachelli)等不喜游泳的鱼类[23],虹鳟和斑马鱼 (Danio rerio) 对乳酸的耐受能力也显著高于不喜游泳的瓦氏黄颡鱼和鲇 (Silurus asotus)[14, 24],可见乳酸清除和糖原恢复时间越短的水生动物游泳能力越强[24]。虾类力竭运动性能显著低于鱼类[19],在比较前人对不同体型对虾跳跃能力的研究中发现,体质量越大、体型越大的虾类跳跃距离越远[19],可认为体型大的成虾逃逸时所需的跳跃次数显著少于体型小的幼虾。
本研究中幼虾腹肌GLY、LA在力竭运动后第8小时即恢复至对照组水平,且肌GLY在恢复后第72小时仍显著高于对照组,力竭后的肌LA上升幅度、清除时间和GLY均高于体型较大的日本囊对虾 (Penaeus japonicus) 和凡纳滨对虾[5, 19],因此可认为幼虾耐受力竭的能力强于成虾,从力竭状态中恢复所需时间亦短于成虾。推测由于幼虾跳跃距离较短,需要更多次跳跃以逃离危险,因此腹部肌肉的糖原储备、乳酸耐受力适应性高于成虾。
研究表明,动物乳酸清除代谢有彻底氧化供能、合成葡萄糖和合成肌糖原3种途径,其中氧化供能的场所在肝胰腺,合成糖原和葡萄糖则在肌肉[25-26]。对虾属于乳酸释放型,力竭运动后通过血淋巴转运肌乳酸到其他组织进一步氧化分解[9, 21]。前人发现肝胰腺是水生动物清除体内乳酸的主要场所[27-28],本研究中幼虾肝LA水平在恢复阶段先升后降,说明幼虾肝胰腺承担清除乳酸的工作,保证机体在力竭运动后的内环境稳定;而肌LA在恢复后第8至第12小时再次上升,可能是肌肉组织将血乳酸重新摄取以合成糖原[29],但虾类体内是否存在此代谢途径还有待验证。
3.2 斑节对虾幼虾力竭运动后组织内转氨酶活性与蛋白质含量变化
谷丙转氨酶 (ALT) 和谷草转氨酶 (AST) 广泛存在于虾类的肝胰腺与心肌中,其活性可作为指示氨基酸代谢强度和肝功能是否健全的生理指标。两种酶联合介导的氨基酸转氨作用是细胞内糖异生作用的关键,在能源物质匮乏的情况下起到填补能量缺口的功能[30-31]。蛋白质作为生命活动的承担者,参与机体的绝大部分活动。
本研究中肌AST、ALT活性和TP水平在力竭运动前后并无显著变化,在恢复过程中AST、ALT活性持续降低,第72小时显著低于对照组水平,说明肌肉中氨基酸代谢被显著抑制,氨基酸并没有被作为能量物质而被分解供能。在恢复前4 h幼虾肝TP显著上升,转氨酶活性较低,可认为此阶段肝胰腺主要工作是合成蛋白。在恢复后第4小时转氨酶活性显著上升,推测是由于未补充外源营养,机体激活氨基酸代谢通路,通过脱氨作用产生丙酮酸和草酰乙酸供能或进入糖异生通路,进行代谢工作或补充能量缺口,与饥饿状态时的能量动员机制相似[31]。
虾类陷入缺氧状态时肝胰腺会合成并释放大量血蓝蛋白进入血淋巴中,以提高机体氧容量[32-34]。推测释放血蓝蛋白与糖异生作用是肝TP在第4至第8小时内显著下降的主要原因。第24小时开始摄入外源营养后,肝TP逐渐恢复到对照组水平。本研究中幼虾在恢复阶段的前24 h禁食,模拟幼虾在放流时遇到的极端情况,在幼虾进行剧烈运动后若无足够的能量摄入,其肝胰腺可在约24 h内通过氨基酸代谢供能以应对代谢能量不足的困境。
3.3 斑节对虾幼虾力竭运动后组织内呼吸代谢变化
呼吸作用是生物消耗能源物质,提供生命活动所需能量的唯一途径,各种代谢相关酶活性反映了呼吸代谢的强度[1-2]。LDH与SDH分别为无氧和有氧呼吸的关键酶,其活性标志了该组织进行无氧与有氧呼吸强度的高低,在氧充足环境中,LDH还能将乳酸逆转化为丙酮酸参与有氧呼吸,实现乳酸彻底分解[35-36]。鱼虾类在力竭运动或较长时间的持续运动后,肌肉中的无氧呼吸强度和能量动员能力显著升高,而生长率和储能物质则显著下降[37-38]。本研究中幼虾腹肌在陷入力竭时LDH活性显著升高,SDH活性显著降低,说明幼虾剧烈运动时肌肉中产生能量较快的无氧呼吸显著增强,而有氧呼吸作用被抑制。肝胰腺LDH与SDH活性在力竭运动前后无显著变化,但SDH活性在第1至第4小时内显著下降,LDH活性则显著上升,出现了与已有报道一致的对虾缺氧时的代谢特征[32]。
但本研究所有过程中幼虾均处于富氧环境,但幼虾在恢复前8 h内仍有陷入缺氧状态的迹象。Savant等[39]研究发现,持续的低pH胁迫会使泥蟹 (Ilyoplax gangetica) 的耗氧率下降;陶易凡等[40]研究发现,克氏原螯虾 (Procambarus clarkii) 在外源性低pH胁迫时,鳃细胞氧交换效率下降,导致机体陷入缺氧状态。推测由于对虾属于乳酸释放型动物[9-10],肌肉在运动后将乳酸释放至血淋巴,导致血淋巴pH下降,摄氧能力下降,使得幼虾在恢复的前4 h内氧交换能力下降并进入缺氧状态。在此阶段中,幼虾机体的能量主要由无氧呼吸提供,代谢平衡极其脆弱,再次出现胁迫或将导致幼虾猝死。
恢复后第8小时,幼虾肌肉中的LDH活性显著高于力竭运动后的水平,SDH活性在前12 h恢复至对照组水平时不再升高,之后下降。可能是由于肌肉组织将乳酸转化为丙酮酸重新利用,同时抑制有氧呼吸,使肌糖原过量恢复。肝胰腺在前24 h内SDH、LDH活性均显著升高,但糖原并无减少。说明幼虾肝胰腺通过转化乳酸、氨基酸参与三羧酸循环以完成清除力竭运动代谢产物的目的,与凡纳滨对虾经历无水运输后的恢复规律吻合[32]。
4. 小结
在不同恢复时间点上,斑节对虾幼虾肝胰腺中各因子指标均差异显著,而肌肉中仅呼吸相关指标差异显著。幼虾机体在力竭运动后恢复的前4 h内,由于内源性pH变化,导致机体摄氧能力降低,在肝胰腺中有氧呼吸强度显著升高的情况下,致使肌肉有氧呼吸被抑制。幼虾为了适应下次力竭运动作出调节:增强无氧、抑制有氧呼吸抑制肌糖原消耗,并通过糖异生作用进行储能,证明力竭运动训练可提高幼虾力竭刺激的生理适应性,但是否能提高其游泳能力还有待验证。
力竭运动会对虾类机体产生久远影响,其产生的刺激在幼虾经历72 h后仍使机体代谢平衡存在较显著差异。在斑节对虾的增养殖工作中,在经济条件允许的范围内可将幼虾尽量培养至较大体型,并对其进行一定强度的运动训练,以48~72 h为周期进行力竭训练,可有效提高幼虾对力竭刺激的生理适应性,有利于提高其养殖及野外的成活率;在进行操作前应加强幼虾的营养尤其是维生素B6和抗应激剂等补剂,使其体内有足够的能量储备和抗氧化能力以应对刺激。
-
图 2 花斑蛇鲻cytb基因序列单倍型的中间连接网络图
圆的大小代表单倍型频率,颜色代表所属群体,空心圆代表未检测到的单倍型,省略了单倍型之间的突变数
Figure 2. Median-joining network for cytb gene sequence haplotypes of S. undosquamis
The size of the circles is proportional to haplotype frequency and the colors represent the corresponding population. Small hollow circles represent missing haplotypes. The mutation steps between two connected haplotypes are omitted.
图 3 花斑蛇鲻cytb序列单倍型核苷酸错配分布曲线
柱状图表示观测值, 虚线为突然扩张模型下的期望值,实线为空间扩散模型下的期望值
Figure 3. Mismatch distribution of cytb haplotypes for S. undosquamis
The observed pairwise differences are shown in bars and the expected values under the sudden expansion model and spatial expansion model are in dash line and solid line, respectively.
表 1 花斑蛇鲻样本信息及cytb基因序列遗传多样性参数
Table 1 Specimen information of S. undosquamis and genetic diversity parameters based on cytb gene sequences
采样点
sampling
site经度/纬度
longitude/
latitude样本量 (N)
number of
samples单倍型数量 (H)
number of
haplotypes多态性位点数 (S)
number of polym-
orphic sites单倍型多样性 (h±SD)
haplotype
diversity核苷酸多样性 (π±SD)
nucleotide
diversity防城港 FCG 108°30'E/21°00'N 32 21 25 0.943 5±0.028 7 0.003 145±0.001 823 北部湾 BBW 107°15'E/19°15'N 34 22 26 0.953 7±0.021 9 0.003 205±0.001 850 西沙 XS 109°24'E/16°40'N 34 31 37 0.992 9±0.009 9 0.003 852±0.002 168 三亚 SY 109°46'E/17°58'N 35 27 41 0.968 1±0.020 3 0.003 724±0.002 103 海口 HK 111°18'E/20°18'N 34 21 32 0.925 1±0.035 8 0.003 154±0.001 824 珠江口 ZJK 114°05'E/21°41'N 34 27 38 0.985 7±0.010 8 0.003 692±0.002 089 汕头 ST 116°55'E/23°00'N 32 23 34 0.949 6±0.029 0 0.003 497±0.001 997 泉州 QZ 119°02'E/24°36'N 31 22 29 0.961 3±0.021 6 0.003 393±0.001 948 总计 total − 266 142 144 0.965 0±0.006 1 0.003 455±0.001 921 表 2 花斑蛇鲻8个地理群体cytb基因序列遗传变异的分子方差分析
Table 2 Analysis of molecular variance for eight populations of S. undosquamis based on cytb gene sequences
变异来源
source of variation自由度
degree of freedom变异百分比
percentage of variation分化系数
F statisticsP 群体间 among populations 7 0.21 0.002 1 0.572 6 群体内 within populations 258 99.79 所有样本 total samples 265 表 3 花斑蛇鲻两两地理群体间cytb基因序列的遗传分化系数 (对角线下方) 及显著性水平 (对角线上方)
Table 3 Pairwise FST (below diagonal) and P values (above diagonal) among geographic populations of S. undosquamis based on cytb gene sequences
群体 population FCG BBW XS SY HK ZJK ST QZ 防城港 FCG 0.773 8 0.372 0 0.342 6 0.288 2 0.540 9 0.779 1 0.807 2 北部湾 BBW −0.013 1 0.394 1 0.460 9 0.304 5 0.650 3 0.846 7 0.659 4 西沙 XS −0.000 9 −0.001 2 0.077 4 0.041 5 0.225 2 0.419 6 0.622 3 三亚 SY −0.000 2 −0.003 9 0.021 1 0.896 7 0.729 2 0.762 6 0.178 8 海口 HK 0.002 4 0.001 9 0.031 8 −0.012 9 0.747 6 0.588 1 0.118 5 珠江口 ZJK −0.006 7 −0.008 1 0.005 7 −0.008 9 −0.009 9 0.896 3 0.247 1 汕头 ST −0.012 5 −0.013 1 −0.002 4 −0.010 5 −0.007 9 −0.012 9 0.695 9 泉州 QZ −0.014 1 −0.009 9 −0.007 4 0.010 5 0.018 9 0.004 5 −0.010 3 表 4 花斑蛇鲻两两地理群体间随机交配假设检验的显著性水平
Table 4 P values of exact test of sample differentiation of S. undosquamis based on cytb gene haplotype frequencies
群体 population FCG BBW XS SY HK ZJK ST 北部湾 BBW 0.014 6 西沙 XS 0.067 7 0.109 4 三亚 SY 0.019 9 0.362 7 0.187 7 海口 HK 0.035 1 0.197 3 0.027 3 0.736 8 珠江口 ZJK 0.285 5 0.254 7 0.692 4 0.292 1 0.181 8 汕头 ST 0.390 6 0.236 1 0.269 5 0.704 2 0.873 0 0.674 6 泉州 QZ 0.527 0 0.310 6 0.128 8 0.303 3 0.486 4 0.392 6 0.945 3 表 5 花斑蛇鲻cytb基因序列核苷酸错配分布分析的参数估计值和中性检验统计值
Table 5 Mismatch distribution parameter estimates and neutrality tests statistics for S. undosquamis based on cytb gene sequences
群体
population错配分布 mismatch distribution 中性检验 neutrality test 突然扩张模型
sudden expansion model空间扩散模型
spatial expansion modelTajima' D Fu's FS 粗糙指数
HRIP 粗糙指数
HRIP D P FS P 防城港 FCG 0.782 9 0.936 2 0.590 0 0.950 1 −1.493 5 0.048 6 −14.339 4 0 北部湾 BBW 0.066 7 0.067 2 0.028 5 0.061 4 −1.492 5 0.048 5 −15.139 0 0 西沙 XS 0.075 9 0.313 3 0.214 5 0.173 7 −1.854 0 0.013 3 −25.721 9 0 三亚 SY 0.049 8 0.188 5 0.151 9 0.513 6 −2.068 2 0.004 2 −23.293 0 0 海口 HK 0.749 0 0.897 7 0.629 0 0.911 5 −1.926 2 0.010 6 −13.437 8 0 珠江口 ZJK 0.050 7 0.047 7 0.045 6 0.077 7 −1.973 3 0.008 2 −24.127 8 0 汕头 ST 0.890 7 0.834 6 0.909 0 0.934 4 −1.909 1 0.011 5 −17.081 6 0 泉州 QZ 0.136 3 0.217 8 0.156 5 0.352 5 −1.680 7 0.026 9 −15.934 9 0 合计 total 0.177 0 0.439 0 0.076 1 0.601 0 −2.555 7 0 −25.510 6 0 -
[1] FROESE R, PAULY D. FishBase[DB/OL]. [2019-04-25]. https://www.fishbase.de/summary/Saurida-undosquamis.html.
[2] 江艳娥, 许友伟, 范江涛, 等. 南海北部陆架水域多齿蛇鲻与花斑蛇鲻的年龄与生长[J]. 中国水产科学, 2019, 26(1): 82-90. [3] 陈再超, 刘继兴. 南海经济鱼类[M]. 广州: 广东科学与技术出版社, 1982: 184-188. [4] 舒黎明, 邱永松. 南海北部花斑蛇鲻生长死亡参数估计及开捕规格[J]. 湛江海洋大学学报, 2004, 24(3): 29-35. doi: 10.3969/j.issn.1673-9159.2004.03.007 [5] 舒黎明, 邱永松. 南海北部多齿蛇鲻生物学分析[J]. 中国水产科学, 2004, 11(2): 154-158. doi: 10.3321/j.issn:1005-8737.2004.02.012 [6] 黄梓荣. 休渔对南海北部多齿蛇鲻资源的影响[J]. 湛江海洋大学学报, 2002, 22(6): 26-31. doi: 10.3969/j.issn.1673-9159.2002.06.006 [7] 卢伙胜, 颜云榕, 侯刚. 等. 2009年度南海渔业资源调查报告[R]. 湛江: 广东海洋大学, 2010. [8] 孙典荣, 林昭进. 北部湾主要经济鱼类资源变动分析及保护对策探讨[J]. 热带海洋学报, 2004, 2(2): 62-68. doi: 10.3969/j.issn.1009-5470.2004.02.008 [9] 陈作志, 邱永松, 徐姗楠, 等. 北部湾花斑蛇鲻生物学特征的演化[J]. 中国水产科学, 2012, 19(2): 321-328. [10] COATES D J, BYRNE M, MORITZ C. Genetic diversity and conservation units: dealing with the species-population continuum in the age of genomics[J/OL]. Front Ecol Evol, 2018, 60 [2019-04-26]. https://www.frontiersin.org/articles/10.3389/fevo.2018.00165/full. doi: 10.3389/fevo.2018.00165.
[11] REISS H, HOARAU G, DICKEY-COLLAS M, et al. Genetic population structure of marine fish: mismatch between biological and fisheries management units[J]. Fish Fish, 2009, 10(4): 361-395. doi: 10.1111/j.1467-2979.2008.00324.x
[12] GOETHEL D R, BERGER A M. Accounting for spatial complexities in the calculation of biological reference points: effects of misdiagnosing population structure for stock status indicators[J]. Can J Fish Aquat Sci, 2017, 74(11): 1878-1894. doi: 10.1139/cjfas-2016-0290
[13] PALSBØLL P J, BERUBE M, ALLENDORF F W. Identification of management units using population genetic data[J]. Trends Ecol Evol, 2007, 22(1): 11-16. doi: 10.1016/j.tree.2006.09.003
[14] MALI K S, KUMAR M V, FAREJIYA M K, et al. Reproductive biology of Saurida tumbil (Bloch 1795) and Saurida undosquamis (Richardson 1848) inhabiting Northwest coast of India[J]. Int J Pure App Biosci, 2017, 5(6): 957-964. doi: 10.18782/2320-7051.6080
[15] CHHANDAPRAJNADARSINI E M, ROUL S K, SWAIN S, et al. Biometric analysis of brushtooth lizard fish Saurida undosquamis (Richardson, 1848) from Mumbai waters[J]. J Entomol Zool Stud, 2018, 6(2): 1165-1171.
[16] NAJMUDEEN T M, SEETHA P K, ZACHARIA P U. Stock dynamics of the brushtooth lizardfish Saurida undosquamis (Richardson, 1848) from a tropical multispecies fishery in the southeastern Arabian Sea[J/OL]. Aquat Living Resour, 2019, 32 [2019-04-26]. https://www.alr-journal.org/articles/alr/abs/2019/01/alr180110/alr180110.html. doi: 10.1051/alr/2019006.
[17] 张俊, 陈国宝, 陈作志, 等. 南沙南部陆架海域渔业资源声学评估[J]. 南方水产科学, 2015, 11(5): 1-10. doi: 10.3969/j.issn.2095-0780.2015.05.001 [18] 许友伟, 陈作志, 范江涛, 等. 南沙西南陆架海域底拖网渔获物组成及生物多样性[J]. 南方水产科学, 2015, 11(5): 76-81. doi: 10.3969/j.issn.2095-0780.2015.05.009 [19] HALL T, BIOSCIENCES I, CARLSBAD C. BioEdit: an important software for molecular biology[J]. GERF Bull Biosci, 2011, 2(1): 60-61.
[20] 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
[21] POSADA D. ModelTest Server: a web-based tool for the statistical selection of models of nucleotide substitution online[J]. Nucleic Acids Res, 2006, 34(sup 2): W700-W703.
[22] ROZAS J, FERRER-MATA A, SÁNCHEZ-DELBARRIO J C, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets[J]. Mol Biol Evol, 2017, 34(12): 3299-3302. doi: 10.1093/molbev/msx248
[23] POLZIN T, DANESCHMAND S V. On Steiner trees and minimum spanning trees in hypergraphs[J]. Oper Res Lett, 2003, 31(1): 12-20. doi: 10.1016/S0167-6377(02)00185-2
[24] BANDELT H J, FORSTER P, ROHL A. Median-joining networks for inferring intraspecific phylogenies[J]. Mol Biol Evol, 1999, 16(1): 37-48. doi: 10.1093/oxfordjournals.molbev.a026036
[25] 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/j.1755-0998.2010.02847.x
[26] TAJIMA F. Statistical-method for testing the neutral mutation hypothesis by DNA polymorphism[J]. Genetics, 1989, 123(3): 585-595.
[27] FU Y X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection[J]. Genetics, 1997, 147(2): 915-925.
[28] EXCOFFIER L. Patterns of DNA sequence diversity and genetic structure after a range expansion: lessons from the infinite-island model[J]. Mol Ecol, 2004, 13(4): 853-864. doi: 10.1046/j.1365-294X.2003.02004.x
[29] RAY N, CURRAT M, EXCOFFIER L. Intra-deme molecular diversity in spatially expanding populations[J]. Mol Biol Evol, 2003, 20(1): 76-86. doi: 10.1093/molbev/msg009
[30] HARPENDING H C. Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution[J]. Hum Biol, 1994, 66(4): 591-600.
[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] KOCHER T D, STEPIEN C A. Molecular systematics of fishes[M]. New York: Academic Press, 1997: 113-128.
[33] FRANKHAM R, BALLOU J D, BRISCOE D A. Introduction to conservation genetics[M]. Cambridge: Cambridge University Press, 2002: 78-104.
[34] 边力, 王鹏飞, 陈四清, 等. 基于线粒体Cytb基因序列的绿鳍马面鲀6个野生群体的遗传结构分析[J]. 中国水产科学, 2018, 25(4): 827-836. [35] 熊丹, 李敏, 李永振, 等. 南海短尾大眼鲷线粒体Cyt b基因序列及种群判别分析[J]. 中国水产科学, 2016, 23(1): 188-197. [36] 夏月恒, 章群, 高志远, 等. 中国近海鮸鱼遗传多样性的细胞色素b全序列分析[J]. 广东农业科学, 2013, 40(3): 101-105. doi: 10.3969/j.issn.1004-874X.2013.03.034 [37] 彭博, 章群, 赵爽, 等. 中国近海小黄鱼遗传变异的细胞色素b序列分析[J]. 广东农业科学, 2010, 37(2): 131-135. doi: 10.3969/j.issn.1004-874X.2010.02.045 [38] 沈朕, 关洪斌, 郑风荣, 等. 基于cytb 和D-loop的4个大泷六线鱼群体遗传多样性分析[J]. 海洋科学进展, 2017, 35(4): 524-534. doi: 10.3969/j.issn.1671-6647.2017.04.009 [39] PINSKY M L, PALUMBI S R. Meta-analysis reveals lower genetic diversity in overfished populations[J]. Mol Ecol, 2014, 23(1): 29-39. doi: 10.1111/mec.12509
[40] GRANT W S, BOWEN B W. Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation[J]. J Hered, 1998, 89(5): 415-426. doi: 10.1093/jhered/89.5.415
[41] WANG P. Response of western Pacific marginal seas to glacial cycles: paleoceanographic and sedimentological features[J]. Mar Geol, 1999, 156(1): 5-39.
[42] IMBRIE J, BOYLE E A, CLEMENS S C, et al. On the structure and origin of major glaciation cycles 1. Linear responses to Milankovitch forcing[J]. Paleoceanography, 1992, 7(6): 701-738. doi: 10.1029/92PA02253
[43] 黄小林, 李文俊, 林黑着, 等. 基于线粒体DNA D-loop序列的黄斑篮子鱼群体遗传多样性分析[J]. 热带海洋学报, 2018, 37(4): 45-51. [44] 李敏, 张鹏, 李玉芳, 等. 南海扁舵鲣种群遗传结构和遗传多样性评价[J]. 南方水产科学, 2015, 11(5): 82-89. doi: 10.3969/j.issn.2095-0780.2015.05.010 [45] 李敏, 李玉芳, 张鹏, 等. 基于线粒体控制区序列的南海圆舵鲣种群遗传结构分析[J]. 南方水产科学, 2016, 12(4): 88-95. doi: 10.3969/j.issn.2095-0780.2016.04.011 [46] HEWITT G M. The genetic legacy of the Quaternary ice ages[J]. Nature, 2000, 405(6789): 907-913. doi: 10.1038/35016000
[47] PALUMBI S R. Genetic divergence, reproductive isolation, and marine speciation[J]. Annu Rev Ecol Syst, 1994, 25: 547-572. doi: 10.1146/annurev.es.25.110194.002555
[48] 苏纪兰. 中国近海水文[M]. 北京: 海洋出版社, 2005: 1-367. [49] 孙冬芳, 董丽娜, 李永振, 等. 南海北部海域多齿蛇鲻的种群分析[J]. 水产学报, 2010, 34(9): 1387-1394. [50] LAIKRE L, PALM S, RYMAN N. Genetic population structure of fishes: Implications for coastal zone management[J]. Ambio, 2005, 34(2): 111-119. doi: 10.1579/0044-7447-34.2.111
-
期刊类型引用(1)
1. 高杰,郭华阳,刘明鉴,刘宝锁,朱克诚,张楠,郭梁,张殿昌. 无乳链球菌感染对卵形鲳鲹脾脏的损伤及细胞凋亡相关基因表达的影响. 南方水产科学. 2022(06): 44-51 . 本站查看
其他类型引用(2)