尼罗罗非鱼对碳酸盐碱度耐受性研究

王艳玲，赵 岩，陈 铭，施郁松，徐 俊，李帅帅，唐首杰，赵金良

（上海海洋大学/水产种质资源发掘与利用教育部重点实验室/水产动物遗传育种中心上海市协同创新中心，上海 201306）

我国盐碱水面积大，广泛分布于东北、西北、华北等内陆地区，水质类型多样，以碳酸盐型居多[1]。过高的碳酸盐碱度对鱼类有强烈的毒害作用，如环境中氢氧根离子 (OH−)、碳酸氢根离子(HCO3−)、碳酸根离子 (CO32−) 等直接作用于鱼类的鳃上皮细胞，造成器质性损伤[2]；鱼类鳃内二氧化碳分压 (pCO2) 随碳酸盐碱度的升高维持在较低水平，导致低碳酸血症和呼吸性碱中毒[3-4]；鱼体内的氨气 (NH3) 不能与氢离子 (H+) 结合排出体外，发生氨中毒等[5]。鱼类可以通过改变鳃组织的形态、调整含氮产物的代谢和排泄等方式适应高碱度环境[6-8]，基因的差异表达也是鱼类应答碱胁迫的途径之一[9-11]。不同鱼类对碳酸盐碱度的耐受性不同，目前国内外已有对包括青海湖裸鲤 (Gymnocypris przewalskii) 在内的20多种鱼类碱耐受性的报道[7-8,12-23]，但多数是在极端高碱环境下经过进化具有生存优势的鱼类，如马加迪湖罗非鱼[7](Oreochromis alcalicus grahami)、金字塔湖克拉克大麻哈鱼[8](Oncorhynchus clarki henshawi)、凡湖塔氏卡拉白鱼[12](Chalcalburnus tarichi) 等，且研究过程中往往胁迫方式比较单一、统计处理也比较简单。

尼罗罗非鱼 (O. niloticus) 生长速度快，是世界上最广泛养殖的鱼类之一，2019年我国养殖产量为1.64×104t[24]。尼罗罗非鱼也是一种广盐性鱼类，具有在盐碱水域中养殖推广的潜力。有关尼罗罗非鱼对氯化钠 (NaCl) 盐度耐受性的研究较多[25-28]，但对碳酸盐碱度耐受性的报道较少。雷衍之等[29]发现当水体总碱度为 48.2 me·L−1时，相较于鲢 (Hypophthalmichthys molitrix)、鳙 (H. nobilis)、草鱼(Ctenopharyngodon idella) 等淡水鱼，尼罗罗非鱼存活时间更长，但该碱度计量单位目前使用较少；赵岩等[30]发现，尼罗罗非鱼在 4和 6 g·L−1的碳酸盐碱水体中胁迫 48 h 无死亡个体，10 g·L−1碱环境下的开始死亡时间明显早于8 g·L−1，但并未计算准确的半致死浓度 (LC50)。此外，在对尼罗罗非鱼进行耐盐碱新品种的选育过程中，为了优化选育和驯化方案，也需要对其盐碱耐受性进行全面深入的评估。

本研究以体质量差异3倍以上的大、小2种规格的尼罗罗非鱼幼鱼为对象，用碳酸氢钠 (NaHCO3)配置不同浓度的碱性水体，通过高浓度碱胁迫致死、急性碱胁迫和慢性碱驯化3种不同的胁迫方法，系统研究了尼罗罗非鱼对碳酸盐碱度的耐受性能，为其盐碱选育和盐碱水养殖提供参考资料。

1 材料与方法

1.1 实验材料

实验自2020年10月1日起于上海海洋大学罗非鱼种质资源试验站进行。实验用鱼为新吉富罗非鱼品系后代。挑选同年繁殖且大、小2种规格 [ 大：(21.56±0.27) g，小：(5.92±0.03) g] 的鱼进行实验。正式实验前暂养24 h，暂养期间不投喂，选取体质健康、活力较强、大小均匀的个体。实验用水提前配置，根据水族箱 (46 cm×38 cm×31.5 cm) 体积计算所需NaHCO3(分析纯) 的质量，将相应NaHCO3加至完全曝气的自来水中，均匀搅拌，连续充氧，待其稳定后使用。

1.2 实验设计

1.2.1 高浓度碱胁迫致死实验 将大、小规格实验鱼由淡水分别直接移入碳酸盐碱度为30和22 g·L−1水体中。发现死鱼立即捞出，记录每尾鱼的死亡时间。每组每种鱼设3个重复，每个重复30尾。

1.2.2 急性碱胁迫实验 根据预实验，大规格实验组设置 8、9、10 g·L−13个碱度梯度，小规格实验组设置 6、7、8 g·L−13个碱度梯度。将暂养后的实验鱼直接放入不同碱度组的水箱内，每组40尾，设3个平行，充分曝气的自来水组为对照，进行96 h急性碱胁迫实验。实验开始后，每4 h观察记录死亡情况。

1.2.3 慢性碱驯化实验 水体起始碱度为 0 g·L−1，分别以+2、+4、+6 g·(L·d)−1的速度增加碳酸盐碱浓度来进行实验。每组40尾鱼，3个重复。实验开始后，每12 h观察记录死亡情况。实验期间用气泵不间断充氧，水温19～24 ℃，不喂食不换水，用虹吸管吸出粪便以保持水体清洁。以实验鱼静卧水箱底部、呼吸停止、用玻璃棒触碰无反应为死亡标准。

1.3 数据分析

利用 MATLAB 8.4、Systat Sigma Plot 10.0和 DPS软件对所得实验数据进行统计分析。

1.3.1 高浓度碱胁迫致死实验 计算平均存活时间 (Mean surviving time, MST, min)，以“平均值±标准差 ()”表示。

1.3.2 急性碱胁迫实验 绘制随时间-碱度变化的累计死亡率观察值图；将累计死亡率转换成概率单位后，用线性回归模型拟合碱度与累计死亡率的回归方程，求出96 h LC50；用互补重对数模型(Complementary log-log model, CLL) 拟合不同碱度处理下大、小规格鱼各时间区间的累计死亡率曲面，模型拟合在DPS数据处理系统下进行[31]。

其中实验鱼在时间 tj(j=12, 24,···, J) 被作用因子 (NaHCO3) 的剂量 di(i=6, 7,···, I) 致死的概率可表达为：

式中：β为描述剂量效应的待估参数；τj为至时间tj的时间效应参数；pij为剂量di在0到j个时间单位内产生的累计死亡概率 (Cumulative mortality)；剂量 di使受试个体在时间区间 [tj−1,tj] 内可能遭受的死亡率称为条件死亡概率 (Conditional mortality)，可表示为：

式 (2) 中 β 的含义与式 (1) 相同。γj为描述时间区间 [tj−1, tj] 内时间效应的待估参数，与 τj在时间上的含义有所区别。γj和β通过二项分布最大似然函数逼近获得。对模型的拟合值与实测值之间的差异，分别采用Pearson卡方检验及Hosmer-Lemoshow统计量来检验模型的拟合度高低。

1.3.3 慢性碱驯化实验 绘制随时间-碱度变化的累计死亡率观察值图；用Logistic模型拟合累计死亡率和胁迫时间的关系曲线，分别求出致死碱度范围和LC50。

2 结果

2.1 高浓度碱胁迫致死实验

将实验鱼分别移入碳酸盐碱水体中，3 h内均全部死亡，但大规格鱼在高胁迫浓度下，MST为 (136.92±3.72) min，仍长于小规格鱼 [(114.08±4.36) min]。

2.2 急性碱胁迫实验

实验期间，对照组呼吸与游动保持正常，水质较清澈，全程未出现死亡个体，而碱胁迫组实验鱼陆续死亡。大规格鱼在碳酸盐碱度为8 g·L−1时，无强烈应激表现，前24 h未出现死亡，随时间延长累计死亡率上升较平缓，96 h内不足20%；碳酸盐碱度为 9 g·L−1时，前 72 h 累计死亡率为 22.5%，72 h后迅速上升，胁迫96 h达到52.5%；碳酸盐碱度为 10 g·L−1时，第 12小时即出现死亡，前 60 h累计死亡率小幅增长，60～72 h骤增，72 h死亡数达一半以上，96 h后仅剩少量存活 (图1-a)；小规格鱼在碳酸盐碱度 6 g·L−1时，胁迫 48 h 未出现死亡，96 h累计死亡率不足50%；碳酸盐碱度为7 g·L−1时，前 24 h 累计死亡率为 0，24～72 h 呈匀速变化 (6%, 14%, 22%, 30%)，72 h 后迅速增至64%；碳酸盐碱度为8 g·L−1时，胁迫前期即开始死亡，出现急游冲撞、窜出水面、眼睛突出、眼白出血、呼吸运动加快等现象，胁迫60 h死亡已超半数，96 h后累计死亡率超过85% (图1-b)，且死亡个体大多具有眼部凸起发黑、身体僵硬、腹腔积水等特征；同时水质浑浊较快，含有大量泡沫。

图1 尼罗罗非鱼随时间-碱度变化的累计死亡率观察值Figure 1 Observed value of cumulative mortality of Nile tilapia with time-alkalinity change

基于线性回归模型拟合碱度与累计死亡率的回归方程，大、小规格急性胁迫96 h的LC50分别为8.95和 6.25 g·L−1(表 1)。

表1 基于线性回归模型尼罗罗非鱼累计死亡率的相关性分析Table 1 Correlation analysis of cumulative mortality of Nile tilapia by linear regression model

进行各碱度各时间段累计死亡率的CLL模型拟合 (图2)。在实验设置的时间与碱度范围内，模拟的死亡概率是一个曲面，曲面的高度和变化均与观察值曲面相近。CLL模型下，大、小规格实验鱼急性胁迫 96 h 的 LC50分别为 9.01和 6.35 g·L−1(表 2)。条件死亡率效应参数 (γ) 和累计死亡率效应参数 (τ) 估计值均随胁迫时间的延长而增大；模型经Pearson卡方检验及Hosmer & Lemeshow统计量检验均显著相关；模型系数β、γj经t检验相关极显著 (P<0.01，表 2)。

表2 尼罗罗非鱼互补重对数模型的参数估计及其显著性检验Table 2 Estimated and tested parameters of complementary log-log model of Nile tilapia

图2 互补重对数模型拟合尼罗罗非鱼随时间-碱度变化的累计死亡率Figure 2 Cumulative mortality of Nile tilapia estimated by CLL model with time-alkalinity change

2.3 慢性碱驯化实验

大规格鱼在3种碱变化速度下，累计死亡率波峰不显著 (图3-a)；小规格鱼在3种碱变化速度下，其累计死亡率出现3个明显的峰值，按时间先后分别为+6、+4和+2 g·(L·d)−1组 (图 3-b)。实验中以 2 g·(L·d)−1的幅度增加碱度，大、小规格实验鱼胁迫 96 h均未出现死亡，碳酸盐碱度达 12 g·L−1时累计死亡率均为 100%；+4 g·(L·d)−1组在碱度达12 g·L−1时，大规格鱼累计死亡半数以上，小规格鱼全部死亡；+6 g·(L·d)−1组在碱度达 12 g·L−1时，大规格鱼累计死亡24%，小规格鱼累计死亡率从0升至100%。

图3 尼罗罗非鱼慢性碱驯化下的累计死亡率观察值Figure 3 Observed value of cumulative mortality of Nile tilapia with gradual alkalinity

利用Logistic模型将尼罗罗非鱼的累计死亡率和胁迫时间按照不同碳酸盐碱度变化速度进行统计拟合 (图4)，得出慢性碱驯化过程中实验鱼的致死碱度范围及半致死碱度 (表3)。

图4 Logistic模型拟合慢性碱驯化下尼罗罗非鱼的累计死亡率Figure 4 Cumulative mortality of Nile tilapia estimated by Logistic model with gradual alkalinity

表3 慢性碱驯化下尼罗罗非鱼累计死亡率相关分析Table 3 Correlation Analysis of cumulative mortality of Nile tilapia with gradual alkalinity

3 讨论

LC50是衡量盐碱水体对鱼毒害影响的指标[9]。常见鱼类的 96 h LC50[13-21]见表 4。本实验急性碱胁迫96 h大规格尼罗罗非鱼的LC50为106.55×10−3mol·L−1，小规格的为 74.40×10−3mol·L−1(不同文献中碱度单位有差异，此处统一按 1 g·L−1=84×10−6mol·L−1计算)。可见尼罗罗非鱼的耐碱能力虽弱于青海湖裸鲤、卡拉白鱼等耐高碱鱼类，但仍强于大多淡水养殖鱼类。另外，当慢性驯化实验中碱度为 12 g·L−1时，+2 g·(L·d)−1组累计死亡率均达 100%；+4 g·(L·d)−1组大规格鱼累计死亡半数以上，小规格鱼全部死亡；+6 g·(L·d)−1组小规格鱼累计死亡率从0升至100%。推测碳酸盐碱度12 g·L−1可能是尼罗罗非鱼碱耐受能力的拐点。

表4 常见鱼类急性胁迫96 h的半致死浓度Table 4 96 h LC50 values of alkalinity of common fish species

大规格鱼的高浓度致死MST、急性碱胁迫96 h LC50[(136.92±3.72) min, 8.95 g·L−1] 均高于小规格[(114.08±4.36) min, 6.25 g·L−1]；慢性碱驯化下，大规格鱼 LC50均高于 10 g·L−1，小规格鱼的低于 10 g·L−1；大规格的生存碱度上限比小规格高1～4 g·L−1，死亡碱度下限与小规格差异不显著。上述结果说明体质量会影响鱼类对环境的耐受性。在幼鱼阶段，体质量增加3倍以上的罗非鱼对碳酸盐碱胁迫的耐受性显著增强。值得注意的是，慢性碱驯化实验中大规格鱼LC50随碱变化速度的增大而上升，小规格鱼则随碱变化速度的增大而下降，究其原因有2个：1) 慢性碱驯化实验中，小规格鱼累计死亡率出现3个明显的峰值，而大规格鱼死亡变化平缓，在碳酸盐碱度达12 g·L−1时，小规格鱼较大规格死亡情况更集中；2) 实际操作中碱度以脉冲式增加，碱度的拟合值较实测值准确度可能较低，拟合曲线并未穿过全部数据点。

本研究在慢性碱驯化实验中发现，同规格鱼在不同碱变化速度下的存活时间不同。+2 g·(L·d)−1组实验鱼的存活时间较长 (大：144 h，小：140 h)，+6 g·(L·d)−1组实验鱼较短 (大：68 h，小：52 h)。+4 g·(L·d)−1组实验鱼存活时间 (大：92 h，小：68 h)和+6 g·(L·d)−1组差异不显著。尽管在较高的碱变化速度下，生存碱度上限、死亡碱度下限发生变化，但实际该数值变化并不会提高耐受性。因此认为在碱驯化过程中，碱度日增加2 g·L−1对尼罗罗非鱼的存活更加友好，养殖中可通过 2 g·(L·d)−1的碱增加使其耐受更高浓度的碳酸盐碱水体。另外，可以推测碱度日增加1 g·L−1时罗非鱼存活时间更长，但驯化时间和成本也会相应增加。

碱度对鱼类的毒性具有剂量效应和时间效应[17]。CLL模型最先应用在不同昆虫的化学及生物杀虫剂毒力测定的数据分析上，优点在于可以揭示时间与剂量间的互作效应，使实验数据更充分更完整[31]。本研究引入CLL模型用于急性碱胁迫下尼罗罗非鱼的“时间-碱度-死亡率”分析。结果发现，和线性回归相比，CLL模型可直接反映不同时间区间的死亡变化情况，拟合程度更高，且胁迫过程中的时间效应在数值上表现为逐渐增强 (大：−23.677 6<−21.714 4，−22.710 7<−20.800 2；小：−12.630 2<−10.923 7，−11.353 7<−9.524 9)，不仅解释了胁迫后期实验鱼死亡数量骤增现象，同时也说明在有效的胁迫时间、碱度范围内，CLL模型适于此类数据的分析。但对于慢性碱驯化而言，时间因子和碱度因子并非2个独立变量，存在线性关系，两者共同作用导致累计死亡率的特定变化，不满足CLL模型基于双独立变量的要求。此外，笔者也试图利用双变量Logistic模型[32]解释慢性碱驯化下碱度、时间对尼罗罗非鱼的致死影响，但由于双变量Logistic模型多应用于温度等连续变化的作用因子，实验点较多，和本实验浓度变化明显不同，因此并不适合。