印度-西太平洋胡椒鲷亚科鱼类DNA条形码及分子系统进化研究
THE DNA BARCODING, CLASSIFICATION, AND PHYLOGENY OF INDO-WESTERN PACIFIC PLECTORHYNCHINAE
-
摘要: 为探讨COⅠ基因作为条形码在胡椒鲷亚科鱼类物种鉴定的可行性, 研究测定了胡椒鲷亚科8种鱼类51个个体线粒体COⅠ基因长度为651 bp的序列, 利用MEGA 5.0计算胡椒鲷亚科种内与种间的遗传距离, 基于最大似然法与贝叶斯法构建了系统进化树。结果显示: 胡椒鲷亚科鱼类种间平均遗传距离(0.142)显著大于种内平均遗传距离(0.003), 物种间遗传距离均大于Hebert推荐的物种鉴定最小种间遗传距离0.020(2%)。系统进化树上, 同一物种不同个体间均能形成独立的单系分支, 表明COⅠ基因可作为胡椒鲷亚科物种鉴定的有效条形码基因。研究同时揭示, 在形态分类上被认为是同种异名的两种胡椒鲷(条纹胡椒鲷与东方胡椒鲷)的COⅠ基因序列差异达到0.070, 有可能是两个独立的物种。此外, 在属级水平, 少棘胡椒鲷属与胡椒鲷属种类之间的平均遗传距离小于胡椒鲷属内部种间的遗传距离, 支持少棘胡椒鲷归属于胡椒鲷属的观点。Abstract: To explore whether COⅠsequences could be used as the DNA barcodes in the identification of Plectorhynchinae species, we analyzed the 651 bp COⅠsequences from 51 individuals of 18 Plectorhynchinae fish in this study. The genetic distances between and within species were calculated using MEGA 5.0, and the molecular phylogenetic trees were constructed using the maximum likelihood and Bayesian inference methods. The results showed that the average inter-species genetic distance (0.142) was 47 times higher than the average intra-species genetic distance (0.003). All the interspecies genetic distances were larger than the minimum value 0.020 (2%) suggested by Hebert for the species identification. In the molecular phylogenetic trees, each individual of the same species could form an independent monophyletic group, indicating that COⅠcould be potentially used as an effective barcode gene in the rapid and accurate identification of Plectorhynchinae species. According to the genetic distance values, we found that two Plectorhinchus species, P. orientalis and P. vittatus which were conventionally regarded as synonyms, could be actually two distinct species because of their large COⅠgenetic variations (0.070). The average genetic distance between Diagramma and Plectorhinchus was lower than the values of intra-Plectorhinchus species, which implied that the genetic variation between Diagramma and Plectorhinchus might be below the significant level. This supported the notion that Diagramma could be classified into the Plectorhinchus genus.
-
Keywords:
- DNA barcodes /
- COⅠgene /
- Plectorhynchinae /
- Phylogeny
-
游泳运动是鱼类逃逸、捕食、繁殖等活动行为的重要实现方式, 为其生存生长及种群繁衍提供了基本保证。依据对氧气的需求, 鱼类游泳运动分为有氧运动和无氧运动2种类型, 前者反映有氧代谢能力, 与日常活动(如巡游、自发运动)和洄游等有关, 常以临界游泳速度(Ucrit)[1, 2]或步法转换速度(Ugt)评价[2], 而后者则反映无氧代谢能力, 与捕食活动、逃逸敌害、穿越激流等相关, 常用匀加速游泳速度(Ucat)或爆发游泳速度评估[3]; 鱼类的有氧和无氧运动能力可能存在权衡[3—5]。自Ucat用于衡量鱼类无氧运动能力起, 不同研究采用加速度均是0.167 cm/s2 [即10 (cm/s)/min][3—8]。采用单一加速度考察不同鱼类Ucat的方式可能具有局限性, 因为不同鱼类的生物学特征(如发育阶段、体重大小、游泳能力及性别等)差异较大, 将影响游泳能力的测定结果。因此, 不同加速度与鱼类Ucat关系的相关研究十分迫切。
鱼类表型存在种内个体变异并且在不同环境条件下维持稳定, 具有重要进化生物学和生态学意义[9]。鱼类处于静止状态、无食物消化吸收时的最低能量消耗称作为静止代谢率(Resting metabolic rate, RMR), 代表了机体基本功能的维持水平和有氧代谢能力的下限, RMR与标准代谢率(Standard metabolic rate, SMR)有相近的内涵, 本研究统一使用RMR; 高强度有氧运动对鱼类施加巨大的生理负荷, 使其对氧的需求达到最大, 表现出最大代谢率(Maximum metabolic rate, MMR), 代表了机体有氧代谢能力的上限, 二者的差值为鱼类的有氧代谢空间(Aerobic scope, AS=MMR–RMR), 它表示为其他生理功能提供的最大代谢潜能[10]。鱼类MMR和AS的测量方法通常包括两种, 一种是Ucrit法[11—14], 别一种则是人为诱导运动力竭的过量耗氧(Excess post-exercise oxygen consumption, EPOC)法[11, 15, 16]。2种方法的本质都是通过高强度游泳运动对机体施加最大生理负荷, 诱导机体产生最强有氧代谢水平。然而, 由于实验方法的操作性及其相应的生态学意义不尽相同, 比较不同研究结果时可能面临一定的方法学局限性, 有关鱼类有氧代谢能力特征诱导方法的比较研究十分匮乏[11]。
鲫(Carassius auratus Linnaeus)是我国常见的鲤科鱼类, 广泛分布于江河湖泊、池塘沼泽等水域中。为考察鲤科鱼类运动能力的个体变异和表型特征的生态关联, 以及不同加速度对匀加速能力测定的影响, 本研究以鲫幼鱼为实验对象, 通过2种方法分别测定该种鱼的有氧代谢参数(RMR、MMR、AS及FAS)和有氧运动能力(Ucrit), 同时也测量在不同加速度(0.083、0.167、0.250、0.333 cm/s2)条件下鲫幼鱼的Ucat, 在Ucrit和Ucat测定过程中分别录像以分析获取相应的Ugt, 分析比较鲫幼鱼的有氧和无氧运动能力特征及表型关联, 比较匀加速游泳能力在不同加速度条件下的响应特征, 为鱼类匀加速游泳能力测定加速度的合理选择提供有益参考, 同时也为我国过渔设施(如鱼道)工程的设计开发提供重要的运动生理学资料。
1. 材料与方法
1.1 实验鱼及其驯化
鲫幼鱼购自于重庆市北碚区歇马人工养殖基地, 运回后将鱼放置于3个实验室循环控温水槽(约250 L)中驯化3周。驯化期间, 每天以商业饲料(中国通威公司)饱足投喂1次, 投喂0.5h后用虹吸法清除残饵和粪便以维持水体环境。驯化水槽的日换水量约为水体的10%, 实验用水为曝气后控温的自来水, 用充气泵不断向驯化水槽充气使溶氧水平大于7 mg/L, 驯化水温为(25.0±0.5)℃(与当时季节水温相近), 光周期为光照14 L﹕10 D。
1.2 实验设计
在驯化结束后, 挑选鱼体健康、体重相近的鲫幼鱼[n=30, 体重(13.08±0.41) g, 体长(7.52±0.08) cm]作为实验对象, 并将实验鱼移入已编号的30个养殖单元格中以便辨别实验鱼。实验之前, 对鱼禁食24h。在第1天先测量实验鱼的RMR, 随后采用游泳代谢仪测定鱼的Ucrit、Ugt和MMR并计算获取AS, 测定完毕后将鱼放回至养殖单元格中并投喂一次约0.5%体重的饵料[17], 以维持实验鱼的体重和生理恢复; 第3天测定EPOC法诱导鲫幼鱼的MMR并在运动力竭后连续测定鱼的恢复代谢率, 共计算10min的过量耗氧(EPOC), 在测定完毕后将鱼放回至养殖单元格中并投喂一次约0.5%体重的饵料; 从第5天开始随后的连续4d中, 于每天上午只测定一个加速度条件下实验鱼的Ucat和Ugt, 本研究设定的加速度包括0.083、0.167、0.250、0.333 cm/s2, 每天测定完毕后对单尾鱼投喂一次约0.5%体重的饵料。对于单尾鱼而言, 在Ucrit和EPOC测定时均有48h机体恢复时间, 在不同加速度条件下测量Ucat的恢复时间也达24h, 并在每次饵料投喂过程中所有实验鱼全部表现快速摄食, 但由于体重控制要求仅对每尾鱼投喂约0.5%体重的饵料。根据研究报道, 在运动力竭1h后鱼类白肌的三磷酸腺苷和磷酸肌酸恢复到运动前水平[17], 力竭运动后8h后机体的乳酸水平也恢复运动前水平[18]。本研究实验参数测定时间的间隔明显长于上述重要能量物质恢复和产物清除时间, 不影响鲫幼鱼游泳能力的测定。实验期间所有环境条件(如温度、溶氧、光周期等)与驯化期间保持一致。
1.3 RMR的测定
本研究采用流水式呼吸代谢仪测定实验鱼的RMR。将鱼从养殖单元格转移至呼吸代谢仪的呼吸室驯化24h, 以消除转移过程带来的影响。驯化结束后, 采用溶氧仪(HQ30d, 美国哈希公司)在第2 天的9:00、11:00、13:00、15:00、17:00、19:00和21:00测定实验鱼的耗氧率。呼吸代谢仪的环境条件(如水温、溶氧和光周期等)与驯化期间保持一致; 测定过程中实验室保持安静, 以减少声音对实验鱼生理状态造成影响。将7个时间点的耗氧率取平均值用来表征实验鱼的RMR [mg O2/(kg·h)], 计算公式如下:
${\rm{RMR}} = \Delta {{\rm{O}}_2} \times v/m$
式中, ΔO2是实验呼吸室(有鱼)出水口与空白呼吸室(无鱼)出水口的溶氧差值(mg O2/L), v是呼吸室流速(L/h), 由呼吸室出水口装满50 mL容量瓶所需时间计算获得, m为实验鱼的体重(kg)。由于实验鱼的体重十分相近, 无需进行代谢率的体重校正。
1.4 Ucrit和MMR的测定
(1) Ucrit测定法: 采用鱼类游泳代谢仪同时测定实验鱼的运动代谢率和Ucrit, 同时通过摄像系统拍摄实验鱼的游泳行为。先将实验鱼装入游泳代谢仪的游泳管中驯化1h, 驯化流速为6 cm/s; 随后每隔20min增加6 cm/s, 直至实验鱼运动力竭。力竭标准: 实验鱼因运动力竭不能抵抗水流而后退停靠于游泳管的尾端筛板, 并且在尾端筛板的停靠时间达20s以上[19]。Ucrit(cm/s)的计算公式如下:
${U_{{\rm{crit}}}} = V + \left( {t/\Delta T} \right) \times \Delta V$
式中, V是实验鱼能够完成20min游泳运动的最大流速, ΔV是流速增量(6 cm/s), ΔT是时间间隔(20min), t是在最高流水速度下实验鱼的游泳历时(持续时间短于20min)。
在测定Ucrit过程中, 同时测量实验鱼的运动代谢率。游泳代谢仪装有1个连接溶氧仪外部水循环系统(流量约为500 mL/min), 当溶氧仪探头连接该系统时游泳代谢仪形成封闭状态, 可监测游泳代谢仪内部水体溶氧水平。每隔20min, 打开游泳仪谢仪尾盖20s随后关闭, 该时间足以保证代谢仪内部与外部之间水体的充分置换, 溶氧水平可迅速恢复至正常状态。通过设定好的电脑程序控制目标流速, 用溶氧仪每隔2min测定一次游泳代谢仪内部水体的溶氧水平(mg/L), 连续测定10次。实验鱼的运动代谢率[MO2, mg O2/(kg·h)]的计算公式如下:
$M{{\rm{O}}_2} = \left( {{S_{\rm{t}}} - {S_{\rm{b}}}} \right) \times 60 \times V/m$
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
在Ucrit测定过程中, 在游泳管正上方安置一台高清摄像机, 录像单尾鱼Ucrit游泳行为。随流速的逐渐增加, 由于鱼类开始动用白肌(无氧代谢)提供游泳动力而打破体内乳酸平衡, 使得体内的乳酸积累一定程度时, 鱼类开始出现“爆发–滑行”的游泳行为变化并且摆尾频率变化明显, 该游泳行为被称为步法转换, 首次出现“爆发–滑行”游泳行为时的游泳速度称作为步法转换速度(Ugt), 评价最大有氧运动能力[2, 20]。虽然Ucrit测定过程中涉及到无氧代谢且饱受争议, 但现有研究仍采用Ucrit或Ugt评价鱼类的有氧运动能力[2]。
(2) EPOC测定法: 通过力竭运动诱导实验鱼的MMREPOC [mg O2/(kg·h)]并计算代谢空间(ASEPOC=MMREPOC–RMR)。EPOC的简要步骤[15]: 将鱼从养殖单元格中转移至力竭运动驱赶装置(Φ外=56 cm, Φ内=33 cm, 水深15 cm)的环行水道(流速~60 cm/s)中; 在环形水道中实验鱼逆流快速游泳, 一段时间后鱼体因高流速无法保持相对稳定游泳而后退, 此时人工连续驱赶实验鱼向前游泳直至其运动力竭; 处于力竭状态的鱼通常失去平衡且用手接触尾部时不再有明显应激反应, 整个过程一般在5min之内完成[21]; 然后将鱼迅速(10s之内完成)转移至呼吸代谢仪的呼吸室中进行30min的恢复代谢测定。溶氧数据采集时间分别为0、1、2、3、4、5、6、7、8、9、10、15、20、25、30min。在恢复代谢测定过程中, 呼吸室流速约为500 mL/min。经计算, 在此流速条件下呼吸室水体置换率达99%的历时小于1min[22]。因此, 将随后数分钟出现的最大耗氧率作为实验鱼的MMR, 其计算公式与RMR一样。鱼类的EPOC是指其在运动恢复期内超过静止状态耗氧水平的额外耗氧量, 在一定程度上反映了鱼类的最大无氧代谢能力, 计算恢复期内高于运动前的差值在时间上的总和即为该时期的EPOC[15]。为便于比较不同鲫个体的代谢恢复特征, 本研究仅选取前10min的EPOC用于数据分析。
1.5 Ucat的测定
本研究共采用了4个加速度(0.083、0.167、0.250、0.333 cm/s2), 每天只测定1个加速度下实验鱼的Ucat。将鱼放入游泳代谢仪的游泳管中驯化1h, 管内流速为6.0 cm/s。驯化结束后, 流速从6.0 cm/s开始按照某一加速度(每秒种递增)逐渐增加游泳代谢仪的流速, 直至实验鱼运动力竭, 该力竭状态被认为是鱼类的最大匀加速游泳能力, 此时对应的流速即为实验鱼的Ucat。游泳代谢仪的流速控制由电脑软件完成, 可即时观察游泳管中即时流速及其变化特征。在测定过程中, 在游泳管正上方安装一台高清摄像机, 录像单尾鱼整个的匀加速游泳过程。对获取的视频资料分析确定实验鱼首次“爆发–滑行”游泳行为的出现时间并结合摆尾频率的变化特征, 以确定在Ucat测定过程中实验鱼的Ugt。
1.6 数据处理与统计分析
采用Excel (2003)软件对所有实验数据作常规计算, 统计软件为SPSS (19.0)软件。实验数据均以“平均值±标准误”(Means±SE)表示, 显著水平为P<0.05。Ucrit测定中的能量代谢参数与游泳能力的关系及Ucat与无氧代谢(Ucat–Ugt)的关系均采用一般线性回归方程。两种方法诱导MMR和AS的比较采用T-检验, 能量代谢参数之间的相关性采用Pearson相关分析。采用单因素方差分析考察不同加速度对游泳能力参数的影响, 若组间存在显著差异则采用最小显著性差异法(Least significant difference, LSD)检验组间数据差异。采用组内相关系数(Intraclass correlation coefficient, ICC)和肯德尔和谐系数(Kendall concordance coefficient, KCC)考察匀加速游泳能力参数在不同加速度条件下的重复性。
2. 结果
2.1 Ucrit中能量代谢与游泳速度的关系
在Ucrit法中鲫幼鱼的MMR和AS均与Ucrit呈正相关, 但RMR与Ucrit不相关(图 1A)。此外, 鲫幼鱼的能量代谢参数(MMR、AS、RMR)与Ugt不相关(图 1B, 所有P>0.05)。
![]()
图 1 鲫幼鱼能量代谢与游泳速度的关系MMR=最大代谢率; AS=代谢空间; RMR=静止代谢率; 游泳速度与代谢率的回归方程如下: MMR=10.678×Ucrit+219.677, R2=0.445, P<0.001; AS=10.898×Ucrit+3.823, R2=0.443, P<0.001Figure 1. The relationships between metabolic rates and swimming velocity in juvenile crucian carpMMR=maximum metabolic rate; AS=aerobic scope; RMR=resting metabolic rate. Regressions equations between swimming velocities and metabolic rates as follows: MMR=10.678×Ucrit+219.677, R2=0.445, P<0.001; AS=10.898×Ucrit+3.823, R2=0.443, P<0.0012.2 两种方法中鲫幼鱼能量代谢参数的比较及相关性
本研究采用了2种常规诱导鱼类最大代谢率的方法, 分别为Ucrit法和EPOC法。结果发现, 两种方法获取鲫幼鱼的MMR和AS均无显著性差异(图 2A, MMR: T=0.200, P=0.843; AS: T=0.220, P=0.827); 同样, 2种方法获得的FAS (Ucrit法: 4.84±0.38; EPOC法: 4.98±0.36)也无明显差异(T=0.593, P=0.558)。虽然2种方法获得能量代谢参数的平均值无差异, 但2种方法的能量代谢参数(MMR、AS和FAS)仍保持较高的重复性(图 2B, MMR: r=0.434, P=0.017; 图 2C, AS: r=0.493, P=0.006; 图 2D, FAS: r=0.808, P<0.001)。
![]()
图 2 两种运动方式诱导鲫幼鱼有氧代谢能力的比较及相关性MMR=最大代谢率; AS=代谢空间; Ucrit=临界游泳速度; EPOC=运动后过量耗氧; FAS=相对代谢空间Figure 2. The comparison and relationship in aerobic capacities induced by two methods in juvenile crucian carpMMR=Maximum metabolic rate; AS =Aerobic scope; Ucrit=Critical swimming speed; EPOC=Excess post-exercise oxygen consumption; FAS=Factorial aerobic scope在Ucrit法中, 鲫幼鱼的RMR与MMR及AS均不相关, 但与FAS呈负相关; MMR与AS及FAS呈正相关, AS与FAS也呈正相关(表 1)。在EPOC法中, 鲫幼鱼的RMR与MMR不相关, 但与AS、FAS、EPOC均呈负相关, MMR与AS、FAS、EPOC均呈正相关(表 1)。AS、FAS、EPOC三者之间呈正相关(表 1)。
表 1 两种方法中鲫幼鱼能量代谢参数之间的相关性Table 1. Correlations among parameters of metabolic rate induced by two exercised methods in juvenile crucian carp测定方法Method 实验参数Parameter 相关性Correlation Ucrit测定法 最大代谢率MMR 代谢空间AS 相对代谢空间FAS 静止代谢率RMR r=0.043 P=0.820 r= –0.197 P=0.298 r= –0.630** P<0.001 最大代谢率MMR r=0.971** P<0.001 r=0.711** P<0.001 代谢空间AS r=0.850** P<0.001 EPOC测定法 最大代谢率MMR 代谢空间AS 相对代谢空间FAS 过量耗氧EPOC 静止代谢率RMR r= –0.326 P=0.079 r= –0.581** P=0.001 r= –0.800** P<0.001 r= –0.470** P=0.009 最大代谢率MMR r=0.959** P<0.001 r=0.775** P<0.001 r=0.544** P=0.002 代谢空间AS r=0.908** P<0.001 r=0.610** P<0.001 相对代谢空间FAS r=0. 519** P=0.003 注: RMR=静止代谢率; MMR=最大代谢率; AS=代谢空间; Ucrit=临界游泳速度; EPOC=运动后过量耗氧; FAS=相对代谢空间。*表示P<0.05, **表示P<0.01 Note: RMR=Resting metabolic rate; MMR=Maximum metabolic rate; AS =Aerobic scope; Ucrit=Critical swimming speed; EPOC=Excess post-exercise oxygen consumption; FAS=Factorial aerobic scope. * means P<0.05 and ** meansP<0.01 2.3 鲫幼鱼的有氧和无氧运动能力
鲫幼鱼的Ucrit低于不同加速度条件下的Ucat, 但加速度对Ucat测定无显著影响(图 3A, F=9.333, P<0.001)。然而, 随加速度的增加,Ugt呈现先增加后下降的变化规律(图 3B, F=87.386, P<0.001)。与Ugt类似, 运动能力比例(Ugt/Ucrit或Ugt/Ucat)也呈现先递增后降低的变化趋势(图 3C, F=58.797, P<0.001)。
![]()
图 3 鲫幼鱼游泳能力参数的比较Ucrit=临界游泳速度; Ucat=匀加速游泳速度; Ugt=步法转换速度Figure 3. The comparisons between parameters of swimming performance in juvenile crucian carp under two exercised methodsUcrit=Critical swimming speed; Ucat=Acceleration swimming speed; Ugt=Gait transition speed不论是Ucrit法还是Ucat法, 鲫幼鱼在两方法中的Ugt均与Ucrit或Ucat呈现正相关, 并且在不同加速度条件下Ugt与Ucat呈现的正相关方向保持不变(表 2)。此外, 鲫幼鱼的匀加速游泳能力参数[Ucat, Ugt, Ugt/Ucat, Ucat–Ugt, (Ucat–Ugt)/Ucat]在不同加速度条件下均保持较高的重复性(表 3)。
表 2 鲫幼鱼最大游泳能力与步法转换速度的回归方程Table 2. Regressions (±SE) of gait transition speed (x) as a function of swimming velocity (y, Ucrit or Ucat), y=a x+b测定方法Method 加速度
Accelerate rate (cm/s2)a b R2 P Ucrit法 / 0.198±0.080 16.380±5.250 0.177 0.021* Ucat法 0.083 0.367±0.062 11.244±5.120 0.552 <0.001** 0.167 0.514±0.141 22.784±11.637 0.335 0.001** 0.250 0.489±0.181 22.925±15.336 0.212 0.011* 0.333 0.773±0.120 –3.907±9.835 0.608 <0.001** 注: Ucrit=临界游泳速度, Ucat=匀加速游泳速度. *表示P<0.05, **表示P<0.01 Note: Ucrit=critical swimming speed, Ucat=Accelerated swimming speed. * denotes P<0.05, ** denotesP<0.01 表 3 鲫幼鱼匀加速游泳能力参数的重复性Table 3. Repeatability of parameters of acceleration swimming performance in juvenile crucian carp实验参数Parameters 组内相关系数Intraclass correlation coefficient (ICC) 肯德尔和谐系数Kendall concordance coefficient (KCC) Ucat ICC=0.790 P<0.001 KCC=0.028 P = 0.472 Ugt ICC=0.573 P=0.002 KCC=0.476 P <0.001 Ugt/Ucat ICC=0.627 P<0.001 KCC=0.577 P <0.001 Ucat–Ugt ICC=0.626 P<0.001 KCC=0.546 P <0.001 ( Ucat–
Ugt)/UcatICC=0.629 P<0.001 KCC=0.591 P <0.001 从加速度0.083 cm/s2逐渐增加至0.250 cm/s2, 鲫幼鱼Ucat的无氧代谢能力(Ucat–Ugt)与Ucat之间均呈现正相关(图 4A-C, 所有P<0.01), 而在加速度0.333 cm/s2下二者不相关(图 4D, r=0.342, P=0.069)。在匀加速能力测定中, 鲫幼鱼的有氧代谢能力(Ugt)与无氧代谢能力(Ucat–Ugt)呈负相关, 即二者存在权衡(图 5, r= –0.607, P<0.001)。
![]()
图 4 不同加速度条件下鲫幼鱼匀加速游泳能力与无氧代谢的关系Figure 4. Relationships between acceleration swimming performance and anaerobic metabolism in juvenile crucian carp at different accelerated rates![]()
图 5 鲫幼鱼的有氧代谢能力(Ugt)与无氧代谢能力(Ucat–Ugt)的相关性Figure 5. Relationship between aerobic metabolic capacity and anaerobic metabolic capacity in juvenile crucian carp3. 讨论
3.1 鲫幼鱼有氧运动能力与能量代谢的关系
随着游泳速度逐渐增加至Ucrit时, 鱼类呼吸循环系统的生理功能处于最大负荷状态, 表现出最大有氧代谢水平, 常用来评价鱼类的有氧运动能力[23]。本研究发现, 鲫幼鱼的运动代谢率随游泳速度的增加而呈指数性上升, 并且其MMR和AS与Ucrit均呈现正相关, 即有氧代谢能力越强的个体可以支撑更强水平的有氧运动, 反之亦然。这种现象不仅存在于鲤科鱼类的种间水平[13], 而且也发现于鲤科鱼类的种内水平[7, 24]。本研究还发现鲫幼鱼的RMR与Ucrit不相关, 这与鲤(Cyprinus carpio Linnaeus)[7]和青鱼(Mylopharyngodon piceusRichardson)[24]的正相关报道并不一致。在Ucrit测定过程中, 鲫幼鱼的Ucrit与Ugt呈正相关, 但Ugt与该种鱼的MMR、AS及RMR均不相关, 这可能是因为鲫幼鱼首次出现“爆发–滑行”游泳行为时(Ugt)的有氧代谢水平低于Ucrit时的MMR, 并且这两个代谢水平间的差距因个体而异, 具有明显个体变异; 另一方面由于游泳速度与有氧代谢之间的拟合方程指数不尽相同, 即相同游泳速度的增加, 指数越大的鲫个体其有氧代谢上升更为明显[11, 19]。
本研究采用了Ucrit法和EPOC法诱导鲫幼鱼的MMR和AS, 发现Ucrit法中MMR和AS的平均值与EPOC法的十分相近且无明显差异。然而, 这只是在组间水平的平均值差异检验, 还可以在组内个体水平上对两种方法进行比较。为此, 本研究将Ucrit法的能量代谢参数(MMR、AS和FAS)分别与EPOC法的进行重复性分析。结果发现, 3个能量代谢参数在两种方法之间保持较高的重复性, 即个体的生理表型保持较强的可重复测量水平并且在种群中的表型排序维持稳定, 表明在排除环境因素(如水温、体重和溶氧)后, 不论Ucrit法和EPOC法诱导鱼类的最大有氧代谢能力均可采用并进行比较。值得注意的是, 不论Ucrit法还是EPOC法均是排除实验鱼消化道中食物消化吸收的影响, 但先前研究发现鱼类最大代谢率的诱导与其代谢模式有关, 食物的消化吸收和游泳运动同时进行的方式可诱导某些鱼类(如鲤)的最大代谢率[23], 这种由2种生理活动叠加诱导的最大代谢率明显大于单一生理活动(运动或消化)的诱导值。因此, 在方法采用时还需要考虑实验鱼的代谢模式类型。
3.2 鲫幼鱼RMR与MMR或AS的关系
RMR是鱼类机体的最小维持代谢水平, 而MMR或AS代表了机体的最大有氧代谢能力, 二者的关系引起研究者的广泛关注[25]。如果自然选择倾向于动物既具有较低的能量维持消耗又具有较强的最大代谢能力的表型特征组合, 那么拥有这种最佳表型组合的动物在自然界中具有极高的适合度[25], 然而这种最佳表型关联的进化结果可能是不存在的[26, 27]。根据内温动物有氧代谢能力的进化理论, 自然选择在倾向于选择个体拥有较强的最大有氧代谢能力同时也会倾向于选择较高的维持代谢, 因为具备较强功能的器官(如发达的呼吸循环系统和运动系统)的机体相应具有更多收益(即具有更高的最大代谢率), 但这些器官也具有相应的代价(即相应较高的能量消耗)[28]。因此, 维持代谢消耗与最大代谢率在理论上将呈正相关。支持该假说的鱼类实验证据已见报道, 例如鲤[7]、青鱼[24]和鳟(Salmo trutta Linnaeus)[25]的RMR与MMR呈现正相关。
对于鱼类种内RMR与MMR呈正相关而言, 在排除体重差异和发育因素外迄今还未找到与之有关器官形态学的证据, 但在生化水平上发现证据: 如鳟肝脏上两种有氧代谢关键酶的活性与RMR及MMR均呈正相关[25]。然而, 也有研究发现南方鲇(Silurus meridionalis Chen)的RMR与MMR及AS不相关[14]; 欧鲈(Dicentrarchus labrax Linnaeus)[5]、鲫[15]和鳊(Parabramis pekinensisBasilewsky)[16]的RMR与MMR不相关, 与AS呈负相关; 大西洋鲑(Salmo salar Linnaeus)的RMR与FAS呈现负相关[29], 暗示该种鱼的RMR与MMR也不相关; 上述研究表明不同鱼类的研究结果差异甚大, 但可总结为两种类型: 一是正相关, 二是不相关。本研究发现, 两种运动法中鲫幼鱼RMR与MMR均不相关, EPOC法中RMR与AS呈负相关, 符合第二种相关类型。在自然界中, 假定自然选择会使鱼类同时具有较强的有氧代谢能力和较高的代谢维持, 并且较强的有氧代谢能力有助于个体的生存、生长和繁衍[30], 理论上拥有该表型个体会获得较高的适合度。因此, 在长进化时间尺度上, 自然选择结果是该鱼类种群的有氧代谢能力表型峰值更高且表型维度变异较窄。然而, 这种理论推测并未得到野外和实验室观测数据的支持, 因为野外和实验室表型特征的个体变异仍然较大[3—5]。另一方面, 在实验室条件下可以获取鱼类的最大有氧代谢能力, 但无法确定鱼类在野外条件下最大有氧代谢能力的使用频率[31], 该功能表型受到自然选择的强度及其与生活史的关联的研究也十分困难, 这也可能是导致维持代谢与最大代谢不相关的原因之一。
3.3 鲫幼鱼的匀加速能力及加速度影响
鱼类匀加速游泳能力的评价研究最早见于2000年报道, 该研究采用的加速度为0.167 cm/s2 [10 (cm/s)/min][3], 但该研究并未对加速度的采用原因进行解释, 即加速度的影响被忽视。本研究发现在4个加速度(0.083、0.167、0.250、0.333 cm/s2)下鲫幼鱼的匀加速能力Ucat不受加速度的影响。同样, 虹鳟(Oncorhynchus mykissWalbaum)[32]和鳊幼鱼的Ucat在较大加速度范围内也保持不变[33, 34], 与本研究结果类似。前期研究发现在较小加度速下的Ucat介于Ucrit与较大加速度下Ucat之间[33, 34], 即加速度越小, Ucat值越接近于Ucrit, 本研究鲫幼鱼Ugt也具有相似变化规律, 该现象可能是因为在较小加速度条件下实验鱼在大部分时间进行有氧运动, Ucrit很大程度上反映了鱼类有氧代谢能力并且可被认为是最小匀加速的特例, 鱼类的低匀加速游泳是一种非常接近有氧运动的游泳类型, 因而低加速度下Ucat值接近Ucrit值[33]。然而, 从另一方面考虑, 如果以过大加速度考察鱼类的匀加速游泳能力, 加速度越快, 实验鱼的运动力竭时间相应越短[33], 那么过高的加速度将影响鱼类匀加速游泳能力的准确评估, 因为过快频率的流速变化可能导致实验鱼提前停止游泳而造成假的运动力竭现象。值得注意的是, Farrell研究的加速度递增时间为每分钟[32], 本研究中及实验室前期研究[33, 34]的加速度递增时间为每秒钟。以每秒钟递增的加速度的设定方式对鱼类匀游泳运动的引导作用更为积极连贯, 可更好反映鱼类的匀加速游泳能力, 同时降低了因流速的突然陡增对鱼产生的生理胁迫影响。
研究认为, 鱼类Ucat包括有氧代谢组分(Ugt)和无氧代谢组分(Ucat–Ugt), 并且2个组分间可能存在权衡, 因为前者是由红肌有氧代谢供能而后者是白肌无氧代谢供能[4, 5]。本研究发现鲫幼鱼的Ugt与Ucat–Ugt呈现负相关, 并且Ucat与Ucat–Ugt呈正关, 表明该种鱼在匀加速游泳能力过程中有氧代谢和无氧代谢之间存在权衡。大西洋鳕(Gadus morhua Linnaeus)的研究数据支持上述权衡观点[3], 但欧鲈的实验结果并不支持[4, 5]。上述权衡观点在人类运动员研究中找到了肌肉形态结构的支持证据, 即短跑运动员的腿部快肌(无氧代谢)拥有比例高于慢肌(有氧代谢), 而马拉松运动员的腿部慢肌(有氧代谢)拥有比例高于快肌(无氧代谢)[35]。同样, 大部分鱼类的肌肉组成也包括常见的红肌和白肌两种类型, 二者的解剖学结构差异明显, 是上述功能权衡的形态学基础[36]。随着游泳速度的逐渐增加, 鱼类的游泳行为特征发生相应变化, 例如: 在较高游泳速度下稳定游泳的控制程度下降, 摆尾幅度和频率增加[20]及“爆发–滑行”行为的发生频率上升[10]; 在生化角度上, 无氧代谢大量参与的证据在于鱼类体内重要器官(如肌肉、肝脏和血液)的乳酸水平明显上升, 但是在静止状态下也保持一定的乳酸水平[18, 33, 34], 代谢底物和产物含量的个体差异可能有助于解释上述权衡现象。然而, 验证上述权衡的不同研究结果出现分歧, 仍有必要结合形态学与生化水平进一步探索权衡现象背后的内在机制。
综上所述, 本研究采用2种运动方法诱导鲫幼鱼有氧代谢能力的观测值不论在组间水平还是个体水平均无方法学的差异, 并且能量代谢参数之间存在表型关联; 加速度对鲫幼鱼的Ucat测定无影响, 建议选择较高的加速度, 加速度对Ugt测定有影响, 在不同加速度条件下鲫幼鱼的匀加速游泳能力保持较高的重复性, 并且有氧代谢与无氧代谢存在权衡。
-
[1] Cheng Q T, Zheng B S. Systematic Synopsis of Chinese Fishes [M]. Beijing: Science Press. 1987, 339341 [成庆泰,郑葆珊. 中国鱼类系统检索.北京:科学出版社. 1987, 339341]
[2] Shen S J. Fishes of Taiwan [M]. Taibei: Taiwan University Press. 1993, 360363 [沈世杰. 台湾鱼类志. 台北: 台湾大学出版. 1993, 360363]
[3] Huang Z G. Species and Distribution of Marine Organisms in China [M]. Beijing: China Ocean Press. 1994, 708709 [黄宗国. 中国海洋生物种类与分布. 北京: 海洋出版社. 1994, 708709]
[4] Meng Q W, Su J X, Miu X Z. Systematics of Fishes [M]. Beijing: Chinese Agricultural Press. 1995, 636640 [孟庆闻, 苏锦祥, 缪学祖. 鱼类分类学. 北京: 中国农业出版社. 1995, 636640]
[5] Mckay R J. Haemulidae. In: Fischer W, Bianchi G (Eds.), FAO Species Identification Sheets for Fishery Purposes. Western Indian Ocean (Fishing Area 51). Vol. 2. FAO, Rome. 1984
[6] Carpenter K E, Niem V H. FAO Species Identification Guide for Fishery Purposes.ISSN 10206868. The Living Marine Resources of the Western Central Pacific. Volume 5: Bony fishes part 3 (Menidae to Pomacentridae) [M]. FAO, Rome. 2001, 29612989
[7] Zhu S H, Zheng W J, Zou J X, et al. Molecular phylogenetic analysis of five Pomadasyidae fish based on mitochondrial cytochrome b sequences [J]. Journal of Tropical Oceanography, 2006, 25(4): 4245 [朱世华, 郑文娟, 邹记兴, 等. 5种石鲈科鱼类细胞色素b基因序列及分子系统分析. 热带海洋学报, 2006, 25(4): 4245]
[8] Ren G, Zhang Q, Qian K C, et al. Sequence analysis of twelve grunt fishes basedon 16S ribosomal RNA gene fragments [J]. Journal of Tropical Oceanography, 2007, 26(3): 4852 [任岗, 章群, 钱开诚, 等. 12种石鲈科鱼类线粒体16S rRNA基因的部分序列分析. 热带海洋学报, 2007, 26(3): 4852]
[9] Sanciangco M D, Rocha L A, Carpenter K E. A molecular phylogeny of the Grunts (Perciformes: Haemulidae) inferred using mitochondrial and nuclear genes [J]. Zootaxa, 2011, 2966: 3750
[10] Tavera J, Acero A, Balart E F, et al. Molecular phylogeny of grunts (Teleostei, Haemulidae), with an emphasis on the ecology, evolution, and speciation history of New World species [J]. BMC Evolutionary Biology, 2012, 12: 57
[11] Hebert P D N, Cywinska A, Ball S L et al. Biological identifications through DNA barcodes [J]. Proceedings of the Royal Society B-Biological Sciences, 2003a, 270: 313322
[12] Hebert, P D N, Ratnasingham S, deWaard J R. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species [J]. Proceedings of the Royal Society B-Biological Sciences, 2003b, 270: 596 599
[13] Emre K, Sevan A, Ali˙ S T DNA barcoding common non-native freshwater fish species in Turkey: Low genetic diversity but high population structuring [J]. Mitochondrial DNA, 2013, 24(3): 276287
[14] Keskin E, Atar H H. DNA barcoding commercially important fish species of Turkey [J]. Molecular Ecology Resources, 2013, 13(5): 788797
[15] Mao Y T, Gan X N, Wang X Z. DNA Barcodes and molecular phylogeny of Botiinae fishes based on the mitochondrial COⅠgene [J]. Acta Hydrobiologica Sinica, 2014, 34(4): 737744 [毛云涛, 甘小妮, 王绪祯. 基于线粒体COⅠ基因的沙鳅亚科鱼类DNA条形码及其分子系统发育研究. 水生生物学报, 2014, 34(4): 737744]
[16] Ward R D, Zemlak T S, Innes B H, et al. DNA Barcoding Australia's fish species [J]. Philosophical Transactions of the Royal Society of B-Biological Sciences, 2005, 360: 18471857
[17] Thompson J D, Higgins D G, Gibson T J. Clusial W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position specific gap penalties and weight matrix choice [J]. Nucleic Acids Research, 1994, 22(22): 46734680
[18] Tamura K, Peterson D, Peterson N, et al. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods [J]. Molecular Biology Evolution, 2011, 28: 27312739
[19] Swofford D P. Phylogenetic Analysis Using Parsimony (and Other Methods) Version 10 [M]. Sinauer, Sunderland, MA. 2003
[20] Huelsenbeck J P, Ronquist F. MrBayes: Bayesian inference of phylogeny [J]. Bioinformatics, 2001, 17: 754755
[21] Ronquist F, Huelsenbeck J P. MrBayes 3: Bayesian phylogenetic inference under mixed models [J]. Bioinformatics, 2003, 19: 15721574
[22] Peng J L, Wang X Z, Wang D. et al. Application of DNA barcoding based on the mitochondrial COⅠ gene sequences in classification of Culter (Pisces: Cyprinidae) [J]. Acta Hydrobiologica Sinica, 2009. 33(2): 271276 [彭居俐, 王绪祯, 王丁, 等. 基于线粒体COⅠ基因序列的DNA条形码在鲤科鲌属鱼类物种鉴定中的应用. 水生生物学报, 2009, 33(2): 271276]
[23] Randall J E, Lim K K P. A checklist of the fishes of the South China Sea [J]. Raffles Bulletin of Zoology Supplement, 2000, (8): 569667
[24] Gell F R, Whittington B. Diversity of fishes in seagrass beds in the Quirimba Archipelago, northern Mozambique [J]. Marine Freshwater Research, 2002, 53: 115121
[25] Gillibrand C J, Harris A R, Mara E. Inventory and Spatial Assemblage Study of Reef Fish in the Area of Andavadoaka, South-West Madagascar (Western Indian Ocean) [J]. Western Indian Ocean Journal of Marine Science, 2007, 6: 183197
[26] Unsworth R K F. 2010. Seagrass meadows of the Wakatobi National Park [M]. In: Clifton, J, R K F Unsworth, D J Smith (Eds.), Marine Conservation and Research in the Coral Triangle: The Wakatobi National Park. Nova Publishers. New York. 2010, 101126
[27] Satapoomin U, Randall J E. Plectorhinchus macrospilus, a new species of thicklip (Perciformes Haemulidae) from the Andaman Sea off southestern Thailand [J]. Phuket Marine Biological Center Research Bulletin, 2000, 63: 916
[28] Froese R, Pauly D. FishBase. World Wide Web Electronic Publication www.fishbase.org, version (9/2013)
[29] Randall J E and Lim K K P. A checklist of the fishes of the South China Sea [J]. The Raffles Bulletin of Zoology, 2000, 8(Supplement): 569667
[30] Ramesh R, Nammalwar P. Database on coastal information of Tamilnadu [M]. Institute for Ocean Management Anna University Chennal-600 025. 2008, 1132
[31] Shao K T, Ho H C, Lin P L, et al. A checklist of the fishes of southern Taiwan, Northern South China Sea [J]. The Raffles Bulletin of Zoology, 2008, 19(Supplement): 233171
-
期刊类型引用(7)
1. 刘铠源，段勇杰，程浩，吕俊贤，谢婷，李娟，关长涛，贾玉东. 斑石鲷(Oplegnathus punctatus)、云龙石斑鱼(Epinephelus moara♀×E.lanceolatus♂)、许氏平鮋(Sebastes schlegelii)游泳行为对比分析研究. 海洋与湖沼. 2024(03): 555-565 . 
百度学术
2. 段勇杰，谢婷，高云涛，王枫林，李梦迪，关长涛，贾玉东. 斑石鲷(Oplegnathus punctatus)游泳性能研究. 海洋与湖沼. 2023(01): 183-193 . 百度学术
3. 张静静，王亚冰，王倩，韩多彩，彭士明. 不同抗流能力大黄鱼(Larimichthys crocea)肌肉转录组学差异分析. 海洋与湖沼. 2023(05): 1507-1516 . 百度学术
4. 李新丹，肖新平，周卫国，罗江，周杨浩，荣义峰，祁洪芳，史建全，杜浩. 青海湖裸鲤的游泳能力及鱼梯通过效果. 渔业科学进展. 2022(02): 44-52 . 百度学术
5. 蔡露，Katopodis Christos，金瑶，黄应平，韩德举，胡望斌，陈小娟，陶江平，侯轶群. 中国鲤科鱼类游泳能力综合分析和应用. 湖泊科学. 2022(06): 1788-1801 . 百度学术
6. 孙春艳，曾令清. 颜色标记对两种鱼类有氧运动能力的影响. 重庆师范大学学报(自然科学版). 2020(03): 67-71 . 百度学术
7. 朱峰磊，张敬敬，龙静，闫玉莲，谢小军. 游泳加速模式对南方鲇和团头鲂运动后代谢恢复的影响. 淡水渔业. 2020(03): 65-71 . 百度学术
其他类型引用(4)
下载:
计量
- 文章访问数: 1571
- HTML全文浏览量: 1
- PDF下载量: 635
- 被引次数: 11