CLONE AND EXPRESSION ANALYSIS OF TRANSCRIPTION COFACTOR LBH-B in GIBEL CARP
-
摘要: Lbh (Limb-bud and heart)基因是脊椎动物中高度保守的转录调控因子, 在早期胚胎发育及某些人类疾病的发病过程中发挥着重要作用。我们前期在银鲫(Carassius gibelio)垂体转录组中筛选到一个在垂体中大量表达的基因lbh-b。为了进一步研究lbh基因在银鲫的表达特征, 首先采用RACE方法克隆了银鲫lbh基因家族的成员lbh-b基因(Cglbh-b)。Cglbh-b的cDNA全长1526 bp, 开放阅读框549 bp, 共编码182个氨基酸。生物信息学分析表明CgLbh-b蛋白与其他脊椎动物的Lbh蛋白同源性在68%以上, 可能也是无序蛋白质家族的成员之一。成体组织RT-PCR分析表明Cglbh-b仅在银鲫的垂体、端脑、卵巢及眼睛中表达。不同胚胎发育时期的表达分析表明, 在受精卵至原肠胚中Cglbh-b转录产物是以母源形式存在的mRNA, 其合子转录起始于尾芽期。胚胎整体原位杂交结果显示从受精后2d到受精后3d, Cglbh-b大量表达于脑和眼睛。此外, 随着卵子成熟Cglbh-b在银鲫垂体中的表达上调。这些结果暗示, Cglbh-b可能在调控银鲫脑和眼睛的发育以及卵子成熟过程中发挥着重要作用。Abstract: LBH (limb-bud and heart) is a novel high-conserved transcription cofactor in vertebrates involved in embryonic development and pathogenesis of some human disease. We screened the lbh-b gene that was abundantly expressed in the pituitary of gibel carp, and cloned the lbh-b (Cglbh-b) cDNA sequence of lbh gene family from gibel carp cDNA library by RACE-PCR. The full-length of Cglbh-b cDNA was 1526 bp with a 549 bp long open reading frame (ORF) coding a 182 amino acid protein. Bioinformatics analysis showed that CgLbh-b protein shared high homology (>68%) with other vertebrate LBH and is one of intrinsic disordered proteins.Cglbh-b was abundantly expressed in pituitary, telencephalon, ovary and eye of adult gibel carp. The expression of Cglbh-b in female pituitaries was 4.66 times higher than that in male pituitaries, and Cglbh-b in pituitary was up-regulated during the process of oocyte maturation. During the early embryonic development, maternal mRNA of Cglbh-b was detected in embryos from 4-cells stage to gastrula stage, and its transcripts were synthesized at bud stage. The results of whole mount in situ hybridization showed that Cglbh-b was distributed on the brain and eyes of embryos from 2 days post fertilization (dpf) to 3dpf. These results suggest that Cglbh-b may play important roles in brain and retina development, and reproductive regulation in oocyte maturation of gibel carp.
-
Keywords:
- Gibel carp /
- lbh-b /
- Embryogenesis /
- Oocyte maturation
-
游泳运动是鱼类逃逸、捕食、繁殖等活动行为的重要实现方式, 为其生存生长及种群繁衍提供了基本保证。依据对氧气的需求, 鱼类游泳运动分为有氧运动和无氧运动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)。
不论是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)。
3. 讨论
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 基于NJ法构建的不同脊椎动物Lbh的系统发育树
下划线所示为CgLbh-b, Bootstraps 1000检验各分支置信度, 只显示置信度大于50的分支
Figure 1. Phylogenetic tree of Lbh in vertebrates based on NJ method
CgLbh-b is underlined. The percentages of replicate trees in the bootstrap test (1000 replicates) are calculated, and values large than 50 are shown next to the branches
图 2 Cglbh-b在银鲫成体组织(A)和胚胎发育过程中(B)的表达分析
β-actin作为阳性对照The β-actin was used as internal control; pituitary: 垂体, hypothalamus: 下丘脑, telencephalon: 端脑, mesencephalon: 中脑, cerebellum: 小脑, medulla oblongata: 延髓, spinal cord: 脊髓, liver: 肝脏, kidney: 肾脏, spleen: 脾脏, heart: 心脏, gut: 肠道, skin: 皮肤, gill: 鳃, ovary: 卵巢, testis: 精巢, retina: 视网膜, 256 cells: 256胞期, high: 高期, dome: 椭球期, 50% epiboly: 50%外胞期, bud: 尾芽期, 4s: 4体节期, 15s: 15体节期, 2 dpf: 受精后2d
Figure 2. Semi-quantitative RT-PCR detects the Cglbh-b expression in adult tissues (A) and in embryos during embryogenesis (B)
图 4 Cglbh-b在银鲫雌雄个体垂体(A)和卵子成熟过程中垂体(B)中的表达分析
♀代表雌性个体; ♂代表雄性个体。*表示与♀1或0h相比具有显著性差异(P<0.05)
Figure 4. The expression of Cglbh-b in pituitaries between female and male individuals (A) or during oocyte maturation (B) by real-time PCR
“♀” represents pituitary of female individual. “♂” represents pituitary of male individual. Asterisks (*) indicate significant differences (P<0.05) between other pituitaries of ♀1 or 0h
表 1 引物序列
Table 1 Primers sequences
引物
Primer引物序列(5′—3′, 下划线为T7启动子)
Sequence (5′-3′, T7 sequence are underlined)用途
Usage5′-RACE-R1 GATGACAAACTGACACTTTCTGTAGGCT 5′RACE-PCR 5′-RACE-R2 CACCCACTGAACTCCTCACGGTCCTTC 5′RACE-PCR 3′-RACE-F1 TTTCCCCTTGCCTTACCGTCTTTATCCT 3′RACE-PCR 3′-RACE-F2 GTGCAGAGCTCCACGCTGTCCATGCT 3′RACE-PCR 5′UPM AAGCAGTGGTATCAACGCAGAGTAC 5′RACE通用引物 3′UPM GTGGTATCAACGCAGAGTACTTTTTT 3′RACE通用引物 β-actin F AGCACGGTATTGTGACTAACTG 内参引物 β-actin R TCGAACATGATCTGTGTCATC 内参引物 lbh-b F1 AGCCTACAGAAAGTGTCAGTTTGTC RT-PCR lbh-b R1 AGGATAAAGACGGTAAGGCAAGG RT-PCR lbh-b F2 TTTGTCCCTTTTCTTTGTATTTCTG 荧光定量PCR lbh-b R2 TATTGCAAAGATGAAATCTGAATGA 荧光定量PCR Anti lbh-b F AGCATGGACAGCGTGGAGCTC 原位杂交 Anti lbh-b R TAATACGACTCACTATAGGGATGGTAGTGGTCCTGATGCAT 原位杂交 Sense lbh-b F TAATACGACTCACTATAGGGAGCATGGACAGCGTGGAGCTC 原位杂交 Sense lbh-bR ATGGTAGTGGTCCTGATGCAT 原位杂交 表 2 LBH基因同源DNA序列
Table 2 The homologous DNA sequences of LBH
物种
Species序列号
GenBank accession No.罗非鱼Tilapia Lbh XP_003442986 慈鲷科鱼Cichlids Lbh XP_005934338 河豚Fugu rubripes Lbh XP_011616972 青鳉Medaka Lbh-like XP_011488553 古比鱼Guppy Lbh XP_007576459 鸡Chick Lbh NP_001026209 非洲爪蟾African clawed toad Lbh NP_001081507 人Human LBH NP_112177 小鼠Mouse Lbh NP_084275 斑马鱼Zebrafish Lbh-a NP_956814 银鲫Gibel carp Lbh-b KT_696607 斑马鱼Zebrafish Lbh-like XP_001336471 古比鱼Guppy Lbh-like XP_007575999 青鳉Medaka Lbh-like isoform X3 XP_011481185 罗非鱼Tilapia Lbh-like XP_003441738 慈鲷科鱼Cichlids Lbh-like XP_005927379 斑马拟丽鱼Zebra mbuna Lbh-like XP_014265505 鸡Chick Lbh-like XP_421381 热带爪蟾Western clawed frog Lbh-like isoform X2 XP_002938037 古比鱼Guppy Lbh-like XP_007574801 比目鱼Southern platfish Lbh-like XP_014327879 斑马鱼Zebrafish Lbh-like XP_009305375 慈鲷科鱼Cichlids Lbh-like XP_005951713 斑马拟丽鱼Zebra mbuna Lbh-like XP_004554795 罗非鱼Tilapia Lbh-like XP_005453463 青鳉Medaka Lbh-like isoform X1 XP_004082740 河豚Fugu rubripes Lbh-like XP_011618576 古比鱼Guppy Lbh-likeisoform X1 XP_007540130 比目鱼Southern platfish Lbh-like XP_005811327 -
[1] Al-A H, Rieger M E, Seldeen K L, et al. Biophysical characterization reveals structural disorder in the developmental transcriptional regulator LBH [J]. Biochemicaland Biophysical Research Communications, 2010, 391(1): 1104—1109
[2] Ai J P, Wang Y Q, Tan K R, et al. A human homolog of mouse Lbh gene, hLBH, expresses in heart and activates SRE and AP-1 mediated MAPK signaling pathway [J]. Molecular Biology Reports, 2008, 35(2): 179—187
[3] Briegel K J, Joyner A L. Identification and characterization of Lbh, a novel conserved nuclear protein expressed during early limb and heart development [J]. Development Biology, 2001, 233(2): 291—304
[4] Briegel K J, Baldwin H S, Epstein J A, et al. Congenital heart disease reminiscent of partial trisomy 2p syndrome in mice transgenic for the transcription factor Lbh [J]. Development, 2005, 132(14): 3305—3316
[5] Ekwall A KH, Whitaker J W, Hammaker D, et al. The rheumatoid arthritis risk gene, LBH, regulates growth in fibroblast-like synoviocytes [J]. Arthritis Rheumatology, 2015, 67(5):1193—1202
[6] Rieger M E, Sims A H, Coats E R, et al. The embryonic transcription cofactor LBH is a direct target of the Wnt signaling pathway in epithelial development and in aggressive basal subtype breast cancers [J]. Molecularand Cellular Biology, 2010, 30(17): 4267—4279
[7] Powder K E, Cousin H, McLinden G P, et al. The transcriptional cofactor Lbh regulates angiogenesis andendochondral bone formation during fetal bone development [J]. Development Biology, 2009, 333(2): 348—358
[8] Powder K E, Cousin H, McLinden G P, et al. A nonsynonymous mutation in the transcriptional regulator lbh is associated with cichlid craniofacial adaptation and neural crest cell development [J]. Molecular and Cellular Biology, 2014, 31(12): 3113—3124
[9] Li W H, Zhou L, Li Z, et al. Zebrafish Lbh-like is required for Otx2-mediated photoreceptor differentiation [J]. International Journal Biological Sciences, 2015, 11(6): 688
[10] Zhou L, Gui J. Natural and artificial polyploids in aquaculture [J]. Aquaculture and Fisheries, 2017, 2(3): 103—111
[11] Zhou L, Gui J F. Karyotypic diversity in polyploid gibel carp, Carassius auratus gibelio Bloch [J]. Genetica, 2002, 115(2): 223—232
[12] Zhu H P, Ma D, Gui J F. Triploid origin of the gibel carp as revealed by 5S rDNA localization and chromosome painting [J]. Chromosome research, 2006, 14(7): 767—776
[13] Li X Y, Zhang X J, Li Z, et al. Evolutionary history of two divergent Dmrt1 genes reveals two rounds of polyploidy origins in gibel carp [J]. Molecular Phylogenetics and Evolution, 2014, 78: 96—104
[14] Gui J F, Zhou L. Genetic basis and breeding application of clonal diversity and dual reproduction modes in polyploid Carassius auratus gibelio [J]. Science China Life Sciences, 2010, 53(4): 409—415
[15] Liu X L, Jiang F F, Wang Z W, et al. Wider geographic distribution and higher diversity of hexaploids than tetraploids in Carassius species complex reveal recurrent polyploidy effects on adaptive evolution [J]. Scientific Reports, 2017, 7(1): 5395
[16] Zhou L, Wang Y, Gui J F. Genetic evidence for gonochoristic reproduction in gynogenetic silver crucian carp (Carassius auratus gibelio Bloch) as revealed by RAPD assays [J]. Journal of Molecular Evolution, 2000, 51(5): 498—506
[17] Yang L, Gui J F. Positive selection on multiple antique allelic lineages of transferrin in the polyploid Carassius auratus[J]. Molecular Biology and Evolution, 2004,21(7): 1264—1277
[18] Wang Z W, Zhu H P, Wang D, et al. A novel nucleo-cytoplasmic hybrid clone formed via androgenesis in polyploid gibel carp [J]. BMC Research Notes, 2011, 4(1): 82
[19] Zhang J, Sun M, Zhou L, et al. Meiosis completion and various sperm responses lead to unisexual and sexual reproduction modes in one clone of polyploid Carassius gibelio [J]. Scientific Reports, 2015, 5: 10898
[20] Li X Y, Zhang Q Y, Zhang J, et al. Extra microchromosomes play male determination role in polyploid gibel carp [J]. Genetics, 2016, 203(3): 1415—1424
[21] Yin J, Xia J H, Du X Z, et al. Developmental expression of CagMdkb during gibel carp embryogenesis [J].The International Journal Developmental Biology, 2007, 51(8): 761—769
[22] Xia W, Zhou L, Yao B, et al. Differential and spermatogenic cell-specific expression of DMRT1 during sex reversal in protogynous hermaphroditic groupers [J]. Molelularand Cellular Endocrinology, 2007, 263(1—2): 156—172
[23] Li X, Romero P, Rani M, et al. Predicting protein disorder for N-, C-and internal regions [J]. Genome Informatics, 1999, 10: 30—40
[24] Katoh K, Standley D M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability [J]. Molecular Biology and Evolution, 2013, 30(4): 772—780
[25] Hall T A. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT [J]. Nucleic Acids Symposium Series, 1999, 41: 95—98
[26] Tamura K, Stecher G, Peterson D, et al. MEGA6: molecular evolutionary genetics analysis version 6.0 [J]. Molecular Biology and Evolution, 2013, 30(12): 2725—2729
[27] Xiao Q, Xia J H, Zhang X J, et al. Type-IV antifreeze proteins are essential for epiboly and convergence in gastrulation of zebrafish embryos [J]. International Journal of Biological Sciences, 2014, 10(7): 715—732
[28] Kuraku S, Meyer A, Kuratani S. Timing of genome duplications relative to the origin of the vertebrates: did cyclostomes diverge before or after? [J].Molecular Biology and Evolution, 2009, 26(1): 47—59
[29] Meyer A, Van de Peer Y. From 2R to 3R: evidence for a fish-specific genome duplication (FSGD) [J]. Bioessays, 2005, 27(9): 937—945
[30] Patino R, Sullivan C V. Ovarian follicle growth, maturation, and ovulation in teleost fish [J]. Fish Physiology and Biochemistry, 2002, 26(1): 57—70
[31] Nagahama Y, Yoshikuni M, Yamashita M, et al. Regulation of oocyte growth and maturation in fish [J]. Current Topics in Developmental Biology, 1995, 30(30): 103—145
[32] Nagahama Y, Yamashita M. Regulation of oocyte maturation in fish [J]. Development Growth and Differentiation, 2008, 50(1): 195—219
[33] Cerda J, Mercade J, Lozano J J, et al. Genomic resources for a commercial flatfish, the Senegalese sole (Solea senegalensis): EST sequencing, oligo microarray design, and development of the Soleamold bioinformatic platform [J]. BMC Genomics, 2008, 9(1): 1
[34] Bobe J, Nguyen T, Fostier A. Ovarian function of the trout preovulatory ovary: new insights from recent gene expression studies [J]. Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology, 2009, 153(1): 63—68
[35] Dyson H J, Wright P E. Intrinsically unstructured proteins and their functions [J]. Nature Reviews Molecular Cell Biology, 2005, 6(3): 197—208
[36] Dunker A K, Oldfield C J, Meng J, et al. The unfoldomics decade: an update on intrinsically disordered proteins [J]. BMC Genomics, 2008, 9:1
[37] Uversky V N. Unusual biophysics of intrinsically disordered proteins [J]. Biochimica et Biophysica Acta (BBA)-Proteins and Proteomics, 2013, 1834(5): 932—951
[38] Chen T, Song J, Chan H S. Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding [J]. Current Opinion in Structural Biology, 2015, 30: 32—42
[39] Uversky V N. A decade and a half of protein intrinsic disorder: biology still waits for physics [J]. Protein Science, 2013, 22(6): 693—724
[40] Munro A D, Scott A P, Lam T. Reproductive seasonality in teleosts: environmental influences [M]. CRC Press. 1990, 264
[41] Migaud H, Davie A, Taylor J. Current knowledge on the photoneuroendocrine regulation of reproduction in temperate fish species [J]. Journal of Fish Biology, 2010, 76(1): 27—68
[42] Nagahama Y. Endocrine regulation of gametogenesis in fish [J]. International Journal of Developmental Biology, 1994, 38(2): 217—229
-
期刊类型引用(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)