黄河高原鳅血红蛋白基因家族全基因组分析及低氧胁迫响应

谭瑾, 郭守全, 刘丹, 张存芳, 聂苗苗, 寇若彬, 姚占雯, 田菲, 祁得林

谭瑾, 郭守全, 刘丹, 张存芳, 聂苗苗, 寇若彬, 姚占雯, 田菲, 祁得林. 黄河高原鳅血红蛋白基因家族全基因组分析及低氧胁迫响应[J]. 水生生物学报, 2025, 49(5): 052506. DOI: 10.7541/2025.2024.0325
引用本文: 谭瑾, 郭守全, 刘丹, 张存芳, 聂苗苗, 寇若彬, 姚占雯, 田菲, 祁得林. 黄河高原鳅血红蛋白基因家族全基因组分析及低氧胁迫响应[J]. 水生生物学报, 2025, 49(5): 052506. DOI: 10.7541/2025.2024.0325
TAN Jin, GUO Shou-Quan, LIU Dan, ZHANG Cun-Fang, NIE Miao-Miao, KOU Ruo-Bin, YAO Zhan-Wen, TIAN Fei, QI De-Lin. GENOME-WIDE ANALYSIS OF THE TRIPLOPHYSA PAPPENHEIMI HEMOGLOBIN GENE FAMILY AND ITS RESPONSE TO HYPOXIA STRESS[J]. ACTA HYDROBIOLOGICA SINICA, 2025, 49(5): 052506. DOI: 10.7541/2025.2024.0325
Citation: TAN Jin, GUO Shou-Quan, LIU Dan, ZHANG Cun-Fang, NIE Miao-Miao, KOU Ruo-Bin, YAO Zhan-Wen, TIAN Fei, QI De-Lin. GENOME-WIDE ANALYSIS OF THE TRIPLOPHYSA PAPPENHEIMI HEMOGLOBIN GENE FAMILY AND ITS RESPONSE TO HYPOXIA STRESS[J]. ACTA HYDROBIOLOGICA SINICA, 2025, 49(5): 052506. DOI: 10.7541/2025.2024.0325
谭瑾, 郭守全, 刘丹, 张存芳, 聂苗苗, 寇若彬, 姚占雯, 田菲, 祁得林. 黄河高原鳅血红蛋白基因家族全基因组分析及低氧胁迫响应[J]. 水生生物学报, 2025, 49(5): 052506. CSTR: 32229.14.SSSWXB.2024.0325
引用本文: 谭瑾, 郭守全, 刘丹, 张存芳, 聂苗苗, 寇若彬, 姚占雯, 田菲, 祁得林. 黄河高原鳅血红蛋白基因家族全基因组分析及低氧胁迫响应[J]. 水生生物学报, 2025, 49(5): 052506. CSTR: 32229.14.SSSWXB.2024.0325
TAN Jin, GUO Shou-Quan, LIU Dan, ZHANG Cun-Fang, NIE Miao-Miao, KOU Ruo-Bin, YAO Zhan-Wen, TIAN Fei, QI De-Lin. GENOME-WIDE ANALYSIS OF THE TRIPLOPHYSA PAPPENHEIMI HEMOGLOBIN GENE FAMILY AND ITS RESPONSE TO HYPOXIA STRESS[J]. ACTA HYDROBIOLOGICA SINICA, 2025, 49(5): 052506. CSTR: 32229.14.SSSWXB.2024.0325
Citation: TAN Jin, GUO Shou-Quan, LIU Dan, ZHANG Cun-Fang, NIE Miao-Miao, KOU Ruo-Bin, YAO Zhan-Wen, TIAN Fei, QI De-Lin. GENOME-WIDE ANALYSIS OF THE TRIPLOPHYSA PAPPENHEIMI HEMOGLOBIN GENE FAMILY AND ITS RESPONSE TO HYPOXIA STRESS[J]. ACTA HYDROBIOLOGICA SINICA, 2025, 49(5): 052506. CSTR: 32229.14.SSSWXB.2024.0325

黄河高原鳅血红蛋白基因家族全基因组分析及低氧胁迫响应

基金项目: 

国家自然科学基金(32371578和31960127)资助

详细信息
    作者简介:

    谭瑾(1999—), 女, 硕士研究生; 研究方向为动物生态学。E-mail: 631999450@qq.com

    通信作者:

    祁得林(1973—), 男, 博士; 研究方向为动物分子遗传与生态学。E-mail: delinqi@126.com

  • 中图分类号: Q344+.1

GENOME-WIDE ANALYSIS OF THE TRIPLOPHYSA PAPPENHEIMI HEMOGLOBIN GENE FAMILY AND ITS RESPONSE TO HYPOXIA STRESS

Funds: 

Supported by the National Natural Science Foundation of China (32371578 and 31960127)

    Corresponding author:
  • 摘要:

    为探究血红蛋白基因家族在黄河高原鳅(Triplophysa pappenheimi)低氧适应中的作用, 研究基于黄河高原鳅基因组数据对血红蛋白基因家族进行了生物信息学分析, 并开展低氧胁迫和基因表达检测。结果显示, 黄河高原鳅血红蛋白基因家族由11个成员组成, 其中包括6个α珠蛋白基因(hbaa1hbaa2hbae1hbae3hbae4hbae5)和5个β珠蛋白基因(hbba1hbba2hbbe1.1hbbe1.2hbbe2)。Motif、Domain 及基因结构分析结果均表明该家族成员具有较高的保守性, 除hbbe2外, 其余基因结构相似。染色体定位分析结果显示, 血红蛋白基因家族成员分布在2条染色体上(Chr_07、Chr_17)。蛋白理化性质分析结果显示, 除hbbe1.2为疏水性不稳定蛋白, 其余基因产物均为疏水性稳定蛋白, 其中hbae3hbbe1.1hbbe2为酸性蛋白; 其余8个为碱性蛋白, 其中α-螺旋是主要的二级结构。亚细胞定位预测结果显示, 除hbae4基因产物位于细胞质中, hbae3位于细胞质和细胞外; 其余血红蛋白基因均位于线粒体中。基因表达研究表明, 在急性低氧胁迫12h (0.3±0.1 mg/L)后, hbaa1、hbaa2、hbba1hbba2基因在肝脏、鳃和血液中的表达量具有上调趋势, 但是差异不显著(P<0.05); 在慢性低氧24h至96h [(3.0±0.1) mg/L]时, hbaa1hbba1hbaa2hbba2基因在肝脏和血液中的表达量达到峰值(P<0.05), 随着胁迫时间延长, 上述基因的表达量均有所下调。在鳃组织中, hbaa1、hbaa2hbba1的最高表达量出现在慢性低氧后的196h, 而hbba2基因表达峰值出现在24h。研究展示了黄河高原鳅血红蛋白基因家族成员的低氧胁迫表达模式, 为黄河高原鳅的低氧适应机制研究积累科学数据。

    Abstract:

    Triplophysa pappenheimi is a representative species of Nemacheilinae, which can adapt to the aquatic environment of Tibetan Plateau. To investigate the role of hemoglobin gene family in T. pappenheimi adaptation to hypoxia, bioinformatics analysis of hemoglobin gene family was performed based on genome data, followed by an examination of expression changes of hemoglobin genes in response to hypoxia. The hemoglobin gene family to T. pappenheimi consists of 11 members, including six α-globin genes (hbaa1, hbaa2, hbae1, hbae3, hbae4, and hbae5) and five β-globin genes (hbba1, hbba2, hbbe1.1, hbbe1.2, and hbbe2). Motif, domain, and gene structure analysis showed high conservation among these genes; with the exception of hbbe2, the others exhibited similar structure features. Chromosomal localization analysis showed that the members of hemoglobin gene family were distributed on two chromosomes (Chr_07 and Chr_17). The results of physicochemical analysis indicated that all the gene products were stable hydrophobic proteins, except hbbe1.2, among which, hbae3, hbbe1.1, and hbbe2 were classified as acidic proteins, and the other eight were basic proteins, with α-helix is the main secondary structure. Subcellular localization prediction showed that hbae3 and hbae4 were located in the cytoplasm and the other hemoglobin genes were found in the mitochondria. The expression of hbaa1, hbaa2, hbba1 and hbba2 genes in liver, gill and blood showed an up-regulated trend after 12h of acute hypoxia stress, though the difference was not significant (P<0.05). The expression levels of hbaa1, hbba1, hbaa2, and hbba2 genes in liver and blood reached the peak at 24 h to 96h under chronic hypoxia (P<0.05). In gill tissue, the highest expression of hbaa1, hbaa2, and hbba1 appeared at 196h after chronic hypoxia, while hbba2 reached its peak at 24h. This study reveals the role of the hemoglobin gene family in hypoxic environment of T. pappenheimi, accumulating valuable scientific data for its conservation.

  • 随着野生渔业资源的日益枯竭, 以及我国内陆淡水池塘养殖的快速发展, 高密度放养成为当前池塘主要养殖方式, 然而由于养殖鱼类活动空间拥挤、范围小, 游泳运动减少, 造成商品鱼质量严重下降, 例如土腥味严重、脂肪堆积、肌肉松弛、鱼体免疫力下降及经济效益低等问题[1]。同时消费者对鱼类品质的要求不断提高, 流水养鱼这种不仅能够提高渔业产量又能改善肌肉营养品质的生态养殖模式正越来越受到人们的关注[2]。池塘内循环流水养殖(Inner-Pond Raceway Aquaculture, IPRA), 简称“跑道养殖”, 于2015年引入我国, 是对传统常规池塘养殖模式进行工程化改造, 形成高密度养殖水槽和生态净水区相结合的养殖系统, 是一种新型的生态高效环保养殖模式[3]。IPRA作为一种有效收集养殖对象残饵粪便和以养殖用水净化后循环利用为主要特征的高效养殖模式, 具有节水、环保、高产、智能化的特点, 符合绿色渔业发展的要求。截至2021年仅浙江省已建成2000余条的IPRA养殖水槽, 推广面积近10 km2 [4]。目前, 草鱼(Ctenopharyngodon idella)、花鲈(Lateolabrax maculatus)、大口黑鲈(Micropterus salmoides)、黄颡鱼(Pelteobagrus fulvidraco)、罗非鱼(Oreochromis mossambicus)及鲫(Carassius auratus)等多个品种已在国内开展了IPRA养殖试验[4], 相关研究主要总结了水质改善[5]、设备改进[6]、营养成分[7]和市场效益[8]等, 而有关水流诱导鱼类运动引起鱼类生长代谢的基础理论研究相对较少。因此, 探索水流因子对养殖对象生长性能、抗氧化能力、生理生化、能量代谢和组织结构等指标的影响, 为IPRA甚至流水养鱼提供重要的技术理论支撑。

    经笔者团队近几年的生产实践, 已证实大口黑鲈是适宜IPRA系统养殖的主要品种之一, IPRA的氮磷排放量及排污系数是传统池塘养殖(UPA)四倍以上, 且IPRA净利润率为(53.96±2.21)%, 显著高于UPA(48.04±2.05)%, 取得了较好的生态和经济效益[8]。目前该养殖模式仍处于探索起步阶段, 许多参数都还处于空白阶段, 而本研究前期养殖试验发现, IPRA养殖系统中当流速控制在(0.3±0.05) m/s时, 其养殖的大口黑鲈相对于UPA能增强能量代谢从而有效抑制大口黑鲈内脏脂肪的沉积。相关研究亦表明适当的游泳运动训练可以增强心血管功能, 促进鱼类生长, 改善其生理代谢机能, 提高鱼类生长率等[6], 但运动量过大又会导致代谢紊乱及免疫能力下降, 增加机体的易感性, 影响机体健康[9], 因此急需探索大口黑鲈最适水流速度。

    本项目拟通过开展IPRA模拟系统中不同流速大口黑鲈养殖的试验研究, 探明水流因子对大口黑鲈生长指标、体营养成分、形体指标、抗氧化指标、代谢指标及组织形态的影响, 探讨IPRA大口黑鲈养殖最适流速, 为IPRA大口黑鲈健康养殖提供技术支撑。

    在浙江省湖州市南浔区勤劳村现代渔业园区内选择1个水域面积约2.2 ha的IPRA模式池塘, 池塘安装有4条流水养殖槽, 生态净水区投放鲢(Hypophthalmichthys molitrix)和鳙(Aristichthys nobilis)等滤食性鱼类, 水面利用生态浮床种植铜钱草(Hydrocotyle vulgaris)和狐尾藻(Myriophyllum verticillatum L.)等耐低温水生植物, 并配有4台水车式增氧机(2.0 kW·h/台), 岸边靠近水槽的一侧建有1个集污池。在IPRA水槽进水端一侧放置4套自主研制的模拟水槽装置, 如图 1所示, 该套装置主要包括3条水槽、变频水泵、整流栅、PVC漂浮管及盖网。前期生产中IPRA养殖调查分析结果显示, 前端流水槽中推水速度基本处于0.3—0.4 m/s, 在该流速下大口黑鲈具有较高的生长速度, 因此参照此流速, 本试验设置3个水流速度处理组, 分别为低流速0.2 m/s (L)、中流速0.4 m/s (M)和高流速0.6 m/s (H), 另外设置空白对照组, 即静水组(C)。其中每条养殖水槽的水深为0.5 m, 长1.2 m, 宽0.3 m, 养殖体积为0.18 m3

    图  1  养殖模拟装置图
    A. 实际养殖装置; B. 养殖装置剖面图
    Figure  1.  Farming simulation device drawing
    A. actual breeding device; B. profile of breeding device

    模拟水槽投放鱼苗时间与IPRA系统中投放鱼苗保持一致, 即在12条模拟水槽内放养规格为(8.13±0.82) g的在IPRA系统中已驯食好的健康大口黑鲈苗种, 每条水槽初始投放苗种22尾, 未开启变频水泵, 以适应模拟水槽装置。另外可能由于本模拟装置为白色, 与IPRA系统养殖水槽为灰色有一定差异, 鱼类存在一定应激, 前2天不摄食, 第3天起开始摄食, 前4天内每条模拟水槽死亡3—5尾, 之后整个养殖实验周期无死亡现象。为了保持每条模拟水槽养殖密度同实际IPRA养殖密度(90尾/m3)保持一致, 在适应阶段的第7天每条水槽随机捞出1—3尾, 使得每条模拟水槽大口黑鲈养殖数量均为16尾/条, 养殖试验从7d适应阶段结束后即2020年5月30日开始, 10月30日结束, 养殖持续154d。每天投喂3次饲料(6:00、12:00和18:00), 达饱食水平, 并做好相关生产记录。

    在试验结束时, 4个水槽中各随机抽取10尾大口黑鲈, 现场测定其体重和体长, 并用消毒后的剪刀、镊子解剖取其内脏、肝脏和脂肪, 称量其重量, 计算大口黑鲈增重率(Weight Gain Rate, WGR, %)、特定生长率(Specific Growth Rate, SGR, %/d)、肝体指数(Hepatosomatic Index, HSI)、脏体指数(Viscerosomatic Index, VSI)和内脏脂肪指数(Visceral Adipose Index, VAI), 同时根据体长和体重计算肥满度(Condition Factor, CF, g/cm3), 并根据记录的饲料投喂量和每条水槽累计出鱼量, 计算饵料系数(Feed Coefficient, FC)和单位水体产鱼量(Fish Yield, FY, kg/m3)计算公式如下:

    $ {\rm{WGR}} = \frac{({W}_{t} - {W}_{0})}{{W}_{0}} \times 100 \qquad\quad$

    $ {\rm{SGR}} = \frac{\mathrm{L}\mathrm{n}\left({W}_{t}\right) - \mathrm{L}\mathrm{n}\left({W}_{0}\right)}{t} \times 100 $

    $ \text{}\text{HSI}=\frac{{W}_{\mathrm{肝}\mathrm{脏}}}{{W}_{t}}\times 100\qquad\qquad\;\;\; $

    $ \text{}\text{VSI}=\frac{{W}_{\mathrm{内}\mathrm{脏}}}{{W}_{t}}\times 100\qquad\qquad\;\;\; $

    $ \text{}\text{VAI}=\frac{{W}_{\mathrm{脂}\mathrm{肪}}}{{W}_{t}}\times 100 \qquad\qquad\;\;\;$

    $ {\rm{CF}}=\frac{{W}_{\mathrm{脂}\mathrm{肪}}}{{W}_{t}}\times 100\qquad\qquad\;\;\; $

    $ {\rm{FC}}=\frac{{W}_{\mathrm{饲}\mathrm{料}}}{{W}_{\mathrm{鱼}}} \qquad\qquad\qquad\;\;\;$

    $ {\rm{FY}}=\frac{{W}_{\mathrm{水}\mathrm{体}}}{{V}_{\mathrm{水}\mathrm{体}}}\qquad\qquad\qquad\;\;\; $

    式中, Wt表示养殖t时间时的体重; W0表示试验开始时体重; L表示体长; W肝脏表示肝脏的重量; W内脏表示内脏组织的重量; W脂肪表示内脏脂肪组织重量; W饲料表示单位产鱼量投喂饲料量; W表示单位产鱼量; W水体表示单位水体产鱼量; V水体表示单位水体。

    在测定生长指标之前, 随机采集各个水槽大口黑鲈各5尾, 尾柄抽血5 mL, 迅速注入10 mL的肝素钠抗凝管中, 于4℃保存, 12h内测定各血清指标, 然后每个水槽中的血清样品指标取平均值作为特定流速的1个重复。血清总蛋白(TP)、谷草转氨酶(AST)、碱性磷酸酶(ALP)、总胆汁酸(TBA)、总胆固醇(TC)和血糖(GLU)采用AURORA/Veras 10自动化工作站和SMT 100V生化自动分析工作站检测。甘油三酯(TG)、游离脂肪酸(EFA)、超氧化物歧化酶(SOD)、过氧化氢酶(CAT)、丙二醛(MDA)和总抗氧化力(T-AOC)均采用南京建成生物工程研究所的试剂盒, 具体测定方法按照试剂盒的说明进行。

    待抽完血后, 立即用75%酒精擦拭鱼体表后用蒸馏水漂洗, 消毒后的剪刀镊子解剖获取背部肌肉和肝脏, 置于冰面上的玻璃培养皿中, 迅速清除肝脏附着的脂肪, 并用PBS缓冲液反复冲洗肌肉和肝脏, 取约1 g肝脏和肌肉同一位置的组织分别放入2 mL的无菌离心管中, 立即放入液氮罐中速冻, 带回实验室后将肌肉和肝脏样品放入–80℃冰箱保存, 用于抗氧化酶、消化酶和脂肪代谢相关酶指标的测定。

    需要检测的抗氧化指标包括超氧化物歧化酶(SOD)、过氧化氢酶(CAT)、丙二醛(MDA)和总抗氧化力(T-AOC), 消化酶包括脂肪酶(LPS)、淀粉酶(AMS)和蛋白酶(PTS), 能量代谢相关酶包括谷氨酸脱氢酶(GDH)、激素敏感脂肪酶(HSL)和丙酮酸激酶(PK), 9种指标的测定均采用南京建成生物工程研究所的试剂盒, 具体测定方法按照试剂盒的说明进行。蛋白浓度用考马斯亮蓝染色法测定, 采用酶标仪测定吸光度值。每个水槽中的酶活指标取平均值作为特定流速的1个重复。

    在测定生长指标之后, 各组随机取5尾大口黑鲈, 首先解剖获取肝脏, 再剖离背部肌肉, 剔除附着杂物后及时称重, 记录数据。暂存于–80℃冰箱中用于测定常规营养成分组成及含量。水分采用直接干燥法测定, 参考GB 5009.3-2010; 粗蛋白采用凯式定氮法测定, 参考GB 5009.5-2010; 粗脂肪采用索氏抽提法测定, 参考GB/T 5009.6-2003; 灰分采用马弗炉灼烧法测定, 参考GB 5009.4-2010。每个水槽中的营养成分取平均值作为特定流速的一个重复。

    组织固定: 在测定生长指标时, 各组随机抽取3条大口黑鲈, 用75%酒精擦拭鱼体表后用蒸馏水漂洗, 消毒后的剪刀镊子获取背部统一位置直径0.5—0.8 cm的肌肉, 放入肌肉专用固定液中; 然后解剖获取肝脏, 取肝部同一位置直径0.5—0.8 cm的组织3份, 1份放入Bouin液中, 2份放入戊二醛电镜固定液中。

    组织切片制样: 肌肉固定液和Bouin液中的肌肉和肝组织固定24h后, 用乙醇逐级脱水, 石蜡包埋后切片, 苏木素-伊红(HE, hematoxylin-eosin)染色, 用Motic BA210光学显微镜观察并拍照。

    扫描电镜制样: 电镜固定液的肝脏经PBS缓冲液洗涤2次, 用4%(W/V)的蔗糖溶液洗涤1次, 用系列梯度酒精脱水, 加入100%酒精重悬后, 吸取少量悬液滴加在玻片上, 将玻片粘在导电胶上, 临界点干燥, 真空喷镀, 在SU8010(Hitachi公司)型扫描电镜下观察与摄片。

    利用数据处理软件Excel 2010进行数据处理, 采用SPSS 20.0软件对文中所得数据进行单因素方差分析(One-way ANOVA), 显著性水平设为P<0.05; 实验结果以平均值±标准差表示。

    表 1可知, 中流速组大口黑鲈的增重率和特定生长率显著高于其他3个处理组(P<0.05), 高流速组的特定生长率显著低于其他3个处理组; 中、高流速肥满度显著低于静水组和低流速组(P<0.05); 肝体指数和脏体指数处理组之间无显著性差异(P>0.05); 大口黑鲈内脏脂肪指数高低顺依次为静水组>低流速组>中流速组>高流速组, 且之间存在显著性差异(P<0.05)。

    表  1  不同处理组大口黑鲈生长性能和形体指标的差异
    Table  1.  Difference of growth indices of M. salmoides (n=10)
    项目Item静水组C低流速组L中流速组M高流速组H
    初重 Initial weight (g)8.82±0.828.85±0.798.80±0.818.82±0.84
    末重Final weight (g)210.1±10.6214.0±10.2225.9±9.8203.8±8.7
    饵料系数FC0.94±0.020.98±0.031.02±0.021.09±0.04
    单位水体产鱼量FY (kg/m3)18.67±1.2219.02±1.3120.08±1.3218.04±1.25
    增重率WGR (%)2382.2±107.8a2418.5±106.2a2567.4±112.8b2311.4±116.1a
    特定生长率SGR (%)1.89±0.035b1.78±0.032b2.14±0.035c1.59±0.029a
    肥满度CF (g/cm3)2.12±0.02b2.05±0.02b1.89±0.02a1.82±0.02a
    肝体指数HSI1.01±0.060.96±0.070.97±0.070.96±0.05
    脏体指数VSI6.89±0.576.79±0.566.95±0.656.84±0.63
    内脏脂肪指数VAI1.25±0.11a0.89±0.08b0.58±0.06c0.33±0.05d
    注: 同行数据肩标不同小写字母表示差异显著(P<0.05); 下同Note: In the same row, values with different small letter superscripts are significantly different (P<0.05). The same applies below
    下载: 导出CSV 
    | 显示表格

    表 2可知, 不同处理组肌肉和肝脏常规营养成分存在较大差异。对于肌肉而言, 中流速组粗蛋白显著高于其他3个处理组(P<0.05), 静水组和低流速组粗脂肪含量显著高于中、高流速组(P<0.05)。对于肝脏而言, 中、高流速组粗蛋白显著高于静水组和低流速组(P<0.05); 粗脂肪含量高低顺依次为静水组>低流速组>中流速组>高流速组, 且之间存在显著性差异(P<0.05); 高流速组粗灰分含量显著高于静水组。

    表  2  不同处理组大口黑鲈肌肉和肝脏常规营养成分比较
    Table  2.  Difference of routine nutritional components in muscle of M. salmoides (n=5; %)
    组织Tissue项目Item静水组C低流速组L中流速组M高流速组H
    肌肉 Muscle粗蛋白 Crude protein18.7±1.54a18.5±1.35a20.7±1.54b18.5±1.58a
    粗脂肪 Crude fat6.72±0.61c6.14±0.52c4.58±0.44b3.27±0.42a
    粗灰分 Ash1.39±0.08a1.33±0.10a1.35±0.11a1.87±0.18b
    水分 Moisture77.6±5.576.5±5.973.5±4.974.2±4.1
    肝脏 Liver粗蛋白 Crude protein23.24±1.84a24.58±2.18a26.14±2.71b27.79±2.32b
    粗脂肪 Crude fat11.28±1.21d9.54±1.05c7.37±0.89b5.19±0.47a
    粗灰分 Ash1.14±0.23a1.29±0.15ab1.25±0.16ab1.45±0.14b
    水分 Moisture62.54±5.8964.32±5.8964.52±5.7165.33±5.27
    下载: 导出CSV 
    | 显示表格

    运动对大口黑鲈血清生化和抗氧化指标产生较大影响(表 3), 血清中甘油三酯和总胆固醇变化趋势一致, 即中、高流速组显著低于静水组和低流速组(P<0.05); 高流速组游离脂肪酸显著低于其他3个处理组(P<0.05), 而总胆汁酸显著高于其他3个处理组(P<0.05); 高流速组总蛋白显著低于其他3个处理组(P<0.05), 而中流速组显著高于静水组和低流速组(P<0.05); 静水组血糖含量显著高于其他3个处理组(P<0.05)。中、高流速组谷草转氨酶显著高于静水组和低流速组(P<0.05), 且高流速组显著高于中流速组(P<0.05); 高流速组的碱性磷酸酶显著高于其他3组(P<0.05)。

    表  3  不同处理组大口黑鲈血清和肝脏生化指标比较
    Table  3.  Difference of biochemical and antioxidant indexes in serum of M. salmoides (n=5)
    项目
    Item
    生化指标
    Biochemical index
    静水组C低流速组L中流速组M高流速组H
    血清生化指标
    Serum biochemical index
    甘油三酯TG (mg/mL)2.35±0.18b2.21±0.17b2.05±0.15a1.93±0.20a
    游离脂肪酸EFA (μmol/L)1.25±0.11b1.12±0.10b1.08±0.09b0.87±0.08a
    总胆固醇TC (mmol/L)8.75±0.79b8.19±0.76b7.57±0.72a7.18±0.75a
    总胆汁酸TBA (μmol/L)14.21±1.21a14.82±1.32a14.65±1.21a17.84±1.41b
    总蛋白TP (g/L)35.6±2.41b34.6±2.50b38.7±2.50c27.98±2.22a
    血糖GLU (mmol/L)6.28±0.54b5.21±0.51a5.44±0.48a4.87±0.45a
    谷草转氨酶AST (U/L)192±11.9a207±12.8a218±11.80b245.73±12.7c
    碱性磷酸酶ALP (U/L)62.21±3.31a65.12±4.12a65.01±3.80a72.00±4.5.b
    血清抗氧化指标
    Serum antioxidant index
    超氧化物歧化酶SOD (U/mL)52.32±4.12a58.89±4.89b69.87±6.02b67.88±5.74b
    过氧化氢酶CAT (U/mL)4.21±0.52a5.04±0.41b5.15±0.48b5.77±0.67c
    丙二醛MDA (mmol/L)18.87±1.51c16.78±1.48b13.25±1.12a15.74±1.25b
    总抗氧化力T-AOC (U/L)178.21±12.91a188.97±13.8a223.14±18.40b189.95±14.87a
    肝脏抗氧化指标
    Liver antioxidant index
    超氧化物歧化酶SOD (U/mg prot)163.11±5.90168.92±6.61173.31±6.92159.89±6.51
    过氧化氢酶CAT (U/mg prot)14.88±1.21a15.81±1.44a18.73±1.52b17.65±1.51b
    丙二醛MDA (U/mg prot)1.72±0.33c1.43±0.41bc1.11±0.25a1.53±0.31b
    总抗氧化力T-AOC (U/mg prot)1.64±0.27a1.75±0.35ac2.12±0.33b1.89±0.32c
    下载: 导出CSV 
    | 显示表格

    在血清抗氧化指标中, 静水组SOD和CAT显著低于其他3个处理组(P<0.05), 且随着水流速度的增加而升高; 静水组MDA含量显著高于其他3个处理组(P<0.05), 而中流速组最低, 显著低于其他3个处理组(P<0.05); 中流速组T-AOC最高, 显著高于其他3个处理组(P<0.05)。

    在肝脏抗氧化指标中, 中、高流速组CAT活性显著高于静水组和低流速组(P>0.05); MDA含量顺序依次为中流速组<低流速组<高流速组<静水组, 且中流速组显著低于其他3个处理组(P<0.05); T-AOC顺序依次为中流速组>高流速组>低流速组>静水组, 且中流速组显著高于其他3个处理组(P<0.05)。

    表 4可知, 大口黑鲈不同组织消化酶活力变异较大, 肝脏和肠道是分泌脂肪酶的主要部位, 胃部是淀粉酶和蛋白酶的主要分泌部位, 其次是肠道, 最后是肝脏。不同处理组之间差异也较大, 肠道脂肪酶活性随着水流速度增强而增加, 中、高流速组显著高于静水组和低流速组(P<0.05); 流速组肠道蛋白酶显著高于对照组(P<0.05), 且中流速组显著高于低流速组和高流速组(P<0.05)。

    表  4  不同处理组下大口黑鲈消化酶活力比较
    Table  4.  Difference of digestive enzyme of M. salmoides (n=5; U/mg prot)
    组织Tissue项目Item静水组C低流速组L中流速组M高流速组H
    肠Intestine脂肪酶 LPS79.82±6.75a86.82±9.10a95.26±8.79b105.61±9.78c
    淀粉酶 AMS93.53±8.8196.81±8.2599.54±9.5497.26±8.91
    蛋白酶 PTS83.62±7.63a96.56±8.72b114.63±9.78c101.22±9.24b
    胃Stomach脂肪酶LPS31.25±3.45a35.57±3.97a42.61±4.15b67.87±6.24c
    淀粉酶AMS152.82±14.21a157.4±15.67a181.82±18.98c168.91±15.42b
    蛋白酶PTS112.35±11.24a116.23±10.14a128.47±11.44b118.60±10.81a
    肝Liver脂肪酶LPS90.61±8.47a92.34±8.89a99.89±8.98b96.83±8.97ab
    淀粉酶AMS68.74±6.21a85.15±8.12b96.52±9.12c108.79±9.59d
    蛋白酶PTS48.85±5.14a52.61±5.60ab58.74±5.47b56.25±6.57b
    下载: 导出CSV 
    | 显示表格

    胃组织中、高流速组脂肪酶和淀粉酶活性显著高于静水组和低流速组(P<0.05), 且高流速脂肪酶活性显著高于中流速组(P<0.05), 而中流速组淀粉酶活性显著高于高流速组(P<0.05); 中流速组蛋白酶活性最高, 显著高于其他3个处理组(P<0.05)。

    对于肝脏组织而言, 中流速组脂肪酶活性最高, 显著高于静水组和低流速组(P<0.05); 淀粉活性顺序依次为高流速组>中流速组>低流速组>静水组, 处理组之间存在显著性差异(P<0.05); 中、高流速组蛋白酶活性显著高于静水组(P<0.05)。

    表 5所示, 肌肉和肝脏中谷氨酸脱氢酶均是随着流速的增加而显著升高(P<0.05); 中流速组肌肉和肝脏中激素敏感脂肪酶最高, 且肌肉激素敏感脂肪酶显著高于其他3个处理组(P<0.05); 中流速组肌肉和肝脏中丙酮酸激酶均显著高于其他3个处理组(P<0.05)。

    表  5  不同处理组下大口黑鲈能量代谢酶活性比较
    Table  5.  Difference of energy metabolizing enzyme of M. salmoides (n=5; U/mg prot)
    组织Tissue项目Item静水组C低流速组L中流速组M高流速组H
    肌肉 Muscle谷氨酸脱氢酶GDH1.32±0.22a1.79±0.28b2.17±0.18c2.52±0.23d
    激素敏感脂肪酶HSL0.45±0.03a0.62±0.05b0.75±0.06c0.65±0.05b
    丙酮酸激酶PK26.54±1.21a35.87±2.32b54.21±3.87c43.5±3.21d
    肝Liver谷氨酸脱氢酶GDH2.25±0.31a2.99±0.33b3.22±0.42c3.45±0.31d
    激素敏感脂肪酶HSL1.25±0.21a1.41±0.19ab1.58±0.25bc1.49±0.26b
    丙酮酸激酶PK11.25±1.21a15.26±1.58b23.58±1.87c18.21±2.19b
    下载: 导出CSV 
    | 显示表格

    图 2所示, 肌纤维直径顺序依次为低流速组>中流速组>静水组>高流速组, 且低流速组显著高于其他3个处理组(P<0.05), 静水组和中流速组均显著高于高流速组(P<0.05)。肌纤维密度随着水流速度增加而显著增加(P<0.05; 表 6)。

    图  2  不同水流速度对大口黑鲈肌肉显微结构(HE)的影响
    D. 肌纤维直径 fiber diameter
    Figure  2.  Effects of different water velocity on the muscle microstructure (HE) of M. salmoides
    表  6  不同水流速度对大口黑鲈肌纤维直径和密度的影响
    Table  6.  The effect of different water velocity on the diameter and density of muscle fiber from M. salmoides (n=5)
    指标Index静水组C低流速组L中流速组M高流速组H
    肌纤维直径Fiber diameter (μm)112.36±16.54b132.14±14.56c113.58±12.39b98.84±8.75a
    肌纤维密度Fiber density (根/单位面积)34.61±2.63a49.51±4.22b66.89±5.98c75.98±6.47d
    下载: 导出CSV 
    | 显示表格

    图 3所示, 不同水流速度均明显影响大口黑鲈肝组织结构。对照组细胞空泡化程度明显高于处理组。随着水流速度增加, 空泡化现象越来越低, 至高流速组出现了部分肝细胞边缘模糊, 细胞排列不整齐、细胞核偏移现象, 甚至组织结构的网状支架被破坏, 一些细胞核聚集在一起。与静水组和高流速组相比, 低、中流速组具有细胞形态正常, 肝血窦、细胞核等结构清晰, 无细胞核空泡化和细胞核偏移聚集等现象。

    图  3  不同水流速度对大口黑鲈肝脏显微结构(HE)的影响
    A. 肝细胞hepatocyte; B. 细胞核nucleus; C. 肝血窦hepatic sinusoid; D. 细胞空泡化cell vacuolation; E. 细胞核偏移nuclear migration
    Figure  3.  Effects of different water velocity on the liver microstructure (HE) of M. salmoides

    水流, 是鱼类生活环境中的重要生态因子, 一定的水流条件会促进鱼类游泳运动, 增强活动代谢, 改善其生理功能, 影响鱼类的摄食、生长和鱼体营养组成[9]。Young等[10]研究发现与没有运动训练的相比, 条纹鲈(Morone saxatilis)在经过1.2—2.4 bl/s流速的逆流运动训练后, 其增重率和特定生长率都有所增加, 而超过2.4 bl/s, 则增重率和特定生长率会降低; 斜带石斑鱼(Epinephelus coioides)在1.0 bl/s流速下运动训练8周后, 增重率显著提高, 而在2.0 bl/s流速下增重率则显著下降[11]。这与本研究结果一致, 即中流速组大口黑鲈的增重率和特定生长率显著高于其他3个处理组, 说明在适宜的流速下鱼类的生长率达到最高, 低流速时鱼类的自发活动就会增多, 高流速时鱼类的游泳状态会不稳定, 导致鱼体疲劳而影响生长, 而持续性的适宜流水刺激一定程度上能够削弱鱼类种群中的攻击行为, 减少自主运动而产生的能量消耗, 有利于鱼体生长[12]

    相关研究表明运动对鱼体形态的影响主要与脂肪沉积有关, 适宜范围内的水流速度可以促进脂肪沉积, 导致肥满度、内脏脂肪指数上升[9], 但超过一定范围的运动, 鱼体会因脂肪分解供能而造成肥满度等指标出现下降[13, 14]。宋波澜等[15]在对多鳞四须鲃(Barbodes schwanenfeldi)的实验中发现, 设置3种流速(0、0.7和2.0 bl/s)条件下, 随流速增大, 鱼体粗脂肪含量显著下降; 而Li等[16]研究发现, 两种流速(1.0和1.5 bl/s)的运动锻炼有利于鱼体内脂肪的沉积, 其肥满度、脏体指数和肝体指数出现一定程度的上升。本研究中大口黑鲈内脏脂肪指数随着流速的增加而显著下降, 中、高流速组肥满度显著低于静水组和低流速组, 与上述研究存在一定差异。这可能是因为养殖品种与养殖环境之间的差异引起, 同时也说明在持续的游泳运动状态下, 大口黑鲈通过分解脂肪供能, 从而有效抑制内脏脂肪的沉积, 这一现象也与Yuan等[6]报道的IPRA养殖系统中大口黑鲈肥满度显著低于传统池塘养殖组的结果较一致。因此, 通过适宜水流刺激来促进鱼体运动, 可以提高鱼体代谢速率, 从而减少脂肪在鱼体内的过量累积, 提高对营养物质的吸收利用效率, 进一步提升商品鱼品质。

    水流运动会导致鱼体能量和物质代谢发生变化, 从而导致脂肪含量、营养成分发生变化, 同时对蛋白质、糖分等也产生影响[17]。目前, 国内外已有较多关于水流运动对鱼类营养成分的影响研究, 但结论不一。Yogata等[14]和Shin等[18]在研究流速对牙鲆(Paralichthysolivaceus)和虹鳟(Oncorhynchus mykiss)影响时发现水流运动未对两种鱼类的常规营养成分产生显著性差异; 而宋波澜等[15]和朱志明[17]发现适宜运动能促进多鳞四须鲃肌肉和肝脏蛋白质沉积但抑制了脂肪沉积。本研究结果表明中流速组肌肉粗蛋白显著高于其他处理组, 静水组和低流速组粗脂肪含量显著高于中、高流速组, 而高流速组肝脏粗蛋白显著高于静水组和低流速组, 粗脂肪含量随着流速的增加而显著降低, 说明不同流速对肌肉和肝脏产生了不同的影响, 中流速的大口黑鲈具有高蛋白低脂肪的特点, 较其他3个处理组肌肉品质有所改善。这主要是由于大口黑鲈随着运动量增加, 其代谢速率提高, 能量消耗加大, 首先大量消耗肝脏中的脂肪, 再分解蓄存在肌肉中的脂肪以供能, 因而造成肌肉和肝脏中脂肪随着运动量增加而降低; 同时肌肉蛋白质含量增加主要可能是由于持续性的运动锻炼促进了蛋白质的合成和转化效率[6]。大量研究也已证实适当运动会提高鱼体对糖类和脂肪的利用率[19], 从而减少蛋白质消耗, 有利于肌肉蛋白质的积累, 如逆流运动能明显提高吉富罗非鱼的蛋白质含量[20], 运动同样可提高黄尾鰤(Seriola lalandi)蛋白质含量[14], 但是当运动量过大时, 即当肌肉中的糖原和脂肪被消耗殆尽时, 则会进一步分解蛋白质以供能, 从而导致蛋白分解大于合成, 蛋白质含量有所降低[15]

    鱼类在持续有氧运动时, 鱼类代谢和心率加快, 血液循环加速, 血清生化指标能比较直观反映鱼类运动时的生理状况[16]。在本试验中, 血清中甘油三酯、总胆固醇、游离脂肪酸、血糖和总蛋白随着流速的增加而下降, 说明在本试验条件下, 运动对蛋白质、脂肪和糖类三大营养物质代谢有着明显的影响, 尤其是在高流速高强度运动状态下, 这些指标显著性降低, 机体通过分解血脂和血糖甚至血蛋白进行供能。这一结果与Li等[21]研究结果不一致, 该研究发现运动训练(0.36 m/s, 4周)会导致鲤(Cyprinus carpio)血清中甘油三酯和总胆固醇显著上升, 这可能是由不同的鱼类品种、饲料营养差异、不同的运动强度及时间等因素造成, 具体作用机理有待进一步探索。本试验结果显示血清中总胆汁酸、谷草转氨酶和碱性磷酸酶则随着流速的增加而升高, 且至高流速组达到显著性水平, 综合肝体比、肝脏脂肪含量、血清甘油三酯、胆固醇、脂肪酸及肝脏组织结构等结果, 表明适当的运动强度会降低肝脏脂肪含量, 而高强度的运动则会引起肝脏组织细胞产生损伤, 从而造成血清相关指标发生显著性变化。这一结果与Azuma等[22]研究结果较一致。

    运动调节免疫功能的主要是通过改变机体心血功能和代谢途径, 进而增加抵抗力, 因此运动训练也会对鱼类的免疫功能产生影响[18]。目前, 关于水流诱导的鱼类运动对鱼类免疫能力影响的研究还不多。Azuma等[22]发现大麻哈鱼(Oncorhynchus keta)在11周水流速度分别为2、13和22 cm/s运动训练后, 显著提高了血清SOD活性; 宋波澜[23]发现红鳍银鲫(Barbodes schwanenfeldi)在低流速(0.1 m/s)和高流速(0.3 m/s)的水流训练下其SOD、T-AOC随着水流速度的增加而升高; 而于丽娟[24]发现中华倒刺鲃(Spinibarbus sinensis)幼鱼免疫酶活性随着水流速度增加而升高, 当流速达到2 bl/s免疫酶活性达到最高, 其后随着流速的增加而又开始下降。本试验结果显示, 中流速组的血清和肝脏T-AOC均最高, 且肝脏中SOD、CAT活性最高, 而MDA含量最低, 说明适宜强度的水流运动可以增强抗应激能力, 改善鱼体机体状况, 提高抗氧化能力, 防止感染疾病, 而过量的运动会降低机体的抗氧化能力。

    水流作为重要的影响因子可以刺激鱼类运动, 增强能量代谢, 从而使消化组织器官产生相应的反应机制[25]。目前对消化能力和生长的关系国内外已进行了大量深入的研究报道[26], 而关于水流对鱼类消化酶影响的研究仍较少。肝脏是消化酶合成的重要场所, 肠和胃是鱼类的主要消化场所, 分泌有大量消化酶, 消化酶活性的提高可以促进鱼类对各种营养物质的吸收利用。赵璐琪等[25]发现罗非鱼在水流运动训练7d和30d后的消化酶活性都显著高于对照组, 且随着流速的增加其消化酶活性也增大, 在2 bl/s组达到最大, 其后则开始降低; 许亚琴[27]发现4个水流速度(0、2、5和8 cm/s)显著影响拉氏鱥(Rhynchocypris lagowskii)幼鱼淀粉酶和脂肪酶活力, 呈现先升后降的趋势, 在5 cm/s流速组活力达到最高。本研究中肠胃脂肪酶活性随着流速增加而显著升高, 至高流速组达到最大, 中流速组胃的淀粉酶和蛋白酶活性则最高; 而对于肝脏, 则是中流速组脂肪酶和蛋白酶活性最高, 高流速组淀粉酶活性最高, 这与以往的报道存在一定差异, 这主要可能是因为不同品种的鱼类及训练环境等因素造成的[16]

    运动代谢是指鱼类自发的或因由水流因子及其他外界环境因子的刺激而产生主动或被动运动所表现出来的能量消耗[17], 而其代谢快慢受相关限速酶控制, 其中谷氨酸脱氢酶(GDH)、激素敏感脂肪酶(HSL)和丙酮酸激酶(PK)分别是蛋白质、脂肪和糖原供能分解的关键酶, 在一定程度上可以反映能量代谢强度[28]。目前, 对GDH的研究较多集中在植物及哺乳动物上, 而对水产动物的研究相关报道较少[29]。在本试验中, 随着流速增加肌肉和肝脏中GDH活性显著增加, 说明肝脏中不断合成GDH, 并运送到肌肉组织中, 推测随着流速增加, 肌肉组织不断增强分解蛋白质能力以供能, 同时也反映出肝脏组织可能存在一定损伤[29]。同样运动训练后虹鳟肌肉中PK活性会显著增加, 对葡萄糖利用率提高了28倍[30]; Anttila等[30]发现运动训练能有效提高褐鳟(Salmo trutta)肝脏的HSL活性, 增强脂肪代谢能力; 魏小岚等[11]试验表明在1.0和1.5 bl/s运动训练下GDH等蛋白代谢酶活性显著低于对照组, 而在高流速组2.0 bl/s则蛋白代谢酶活性则显著增强, 蛋白质分解供能进一步增强。可见, 运动对鱼类糖类、脂肪和蛋白质代谢整体上具有显著影响, 但研究结果也存在一些差异。本试验结果显示, 运动能显著增加能量代谢相关酶活, 且高流速组肌肉和肝脏中蛋白脂酶和HSL活性最高, 而PK则是中流速组最高, 其原因可能是大口黑鲈超过一定运动强度后, 肌肉和肝脏中储存的糖原被快速分解殆尽, 主要依靠脂肪和蛋白质供能。综上, 持续且适度的水流运动虽然能增强鱼体蛋白质消耗, 但是合成大于分解而使肌肉蛋白质累积, 同时适度消耗鱼体储存的脂肪和糖原, 改善机体生理代谢。

    研究表明适宜的水流运动训练能有效改善肌肉的结构和形态, 增加鱼体肌纤维数目或直径[31], 促进其生长, 提高游泳能力[32]。如虹鳟在一定程度的游泳锻炼下肌纤维更加发达, 纤维数目和直径分别显著高于对照组[33], 欧鳟(Salvelinus confluentus)在3种水流速度(1.5、3.0和4.5 bl/s)运动训练4周后, 肌纤维直径比对照组分别提高40%、18%和20%, 数目分别提高14%、31%和15%[34], 这一现象也在溪红点鲑(Salvelinus fontinalus)的研究中发现[35]。而有报道显示金鱼在1.5 bl/s流速下, 肌纤维数目比对照组提高25%, 但在3.0和4.5 bl/s流速下则无显著性差异[36], 也有研究发现中华刺鲃经过1.0和2.0 bl/s的8周运动训练未对肌纤维直径和密度产生明显差异, 但在4.0 bl/s运动训练组显著降低了肌纤维直径, 显著增加了肌纤维密度[37]。综上可见有关运动对鱼肌纤维结构的影响并不一致, 主要受运动强度、运动时间及鱼类品种等因素的影响[35]。本试验结果显示低流速肌纤维直径最大, 其次为中流速、静水组和高流速组, 而纤维密度随着水流速度的升高而显著增加, 这与Harimana等[38]结果一致, 暗示适度的运动可以刺激肌纤维直径增大, 促进其增生, 提高生长率, 但是过高的水流刺激运动又会导致肌纤维直径缩小, 研究者认为由于这种较小的肌纤维更有利于游泳过程中氧气的扩散[36], 从而能够适应高强度运动下肌肉消耗氧气分解供能的需求。另外适度运动训练促进鱼类肌纤维直径和数目的增加是导致其生长率显著提高的主要原因[39, 40], 这与本试验中流速组最有最高的生长率结果一致。

    肝脏是鱼类进行各种代谢途径的主要器官, 极易受到外界环境的影响, 对外部应激极其敏感[41]。当鱼类肝脏受到外界环境胁迫时, 其遭到损伤后的典型变化包括水肿变性、空泡化和肝细胞坏死等现象[42], 根据肝脏形态学的变化可以评估其所受的损伤程度。目前关于运动对鱼类肝脏组织形态的影响研究报道还不多。在本试验中, 与静水组和高流速组相比, 低、中流速具有细胞形态正常, 肝血窦和细胞核等结构清晰, 无细胞核空泡化和细胞核偏移聚集等现象, 说明静水组和高流速组均会导致肝脏组织受损, 而长期适度运动则有利于预防脂肪肝, 改善肝脏生理机能[43], 这与本试验中流速组具有最高的消化酶和免疫酶活性相一致。

    本实验结果表明不同的运动强度会对大口黑鲈的生理生化状态产生不同的影响, 中流速组特定增重率和特定生长率最高, 而高流速组最慢; 中、高流速组具有条形好, 粗蛋白含量高的特点; 同静水组、低流速组相比, 中、高流速组可以增强机体消化能力和抗氧化能力, 且以中流速组改善效果最好; 同低、中流水组相比, 静水组和高流速组肝脏结构产生了轻微损伤。综上可知, 中流速组能够促进大口黑鲈生长, 提高机体抗氧化能力, 改善机体代谢, 增加粗蛋白含量, 抑制脂肪沉积, 有效改善肌肉营养品质。因此, 在本研究中, 从生理生化等指标及生产实际考虑, 0.4 m/s可以作为大口黑鲈IPRA养殖系统中最适养殖水流速度。

  • 图  1   斑马鱼和黄河高原鳅血红蛋白系统进化树

    Figure  1.   Phylogenetic tree of hemoglobin from Danio rerio and Triplophysa pappenheimi based on the protein sequences

    图  2   黄河高原鳅血红蛋白三级结构模型

    Figure  2.   Protein tertiary structure model of Hemoglobin in Triplophysa pappenheimi

    图  3   黄河高原鳅血红蛋白基因家族的Motif基序(a)、Domain结构域(b)及基因结构(c)

    Figure  3.   Motif (a), domain (b), and gene structure (c) of the hemoglobin gene family in Triplophysa pappenheimi

    图  4   黄河高原鳅血红蛋白基因在染色体上的分布

    Figure  4.   Distribution of hemoglobin gene on chromosomes in Triplophysa pappenheimi

    图  5   血红蛋白基因在不同组织的表达热图

    Figure  5.   Heat map of hemoglobin gene expression in different tissues

    图  6   血红蛋白各基因在各组织中的表达量

    a. 肝脏组织中的相对表达量; b. 鳃组织中的相对表达量; c. 血液中的相对表达量; 相同字母表示差异不显著(P>0.05); 不同字母表示差异显著(P<0.05)

    Figure  6.   Expression levels of hemoglobin gene in various tissues

    a. Relative expression in liver; b. Relative expression in gill; c. Relative expression in blood; The same letter indicates no significant difference (P>0.05), while different letters indicate significant differences (P<0.05)

    表  1   用于qRT-PCR 分析的引物序列

    Table  1   Primer sequences for qRT-PCR analysis

    基因
    Gene
    序列
    Sequence (5′—3′)
    退火温度
    Annealing
    temperature (℃)
    β-actin F GAACCCCAAGGCTAACAGAGAA 56.9
    β-actin R AGGCATACAGGGACAGCACA 59.1
    hbae1 F TGGTGTTGGTGAAGCTGTTGAGAA 58.8
    hbae1 R GCCAGAACGACGAGTATGTTGTGA 59.3
    hbae3 F GAGCCAGGAACTTATCCACAGACA 58.7
    hbae3 R CCTCAGACCAAGACCTACTTCTCC 58.8
    hbae4 F CATTTACGCTGCCACGGGAAGAA 60.4
    hbae4 R TGCTGAGGTGTGCCAGTGTAGT 61.2
    hbae5 F TCTCACTGGAGCGACCTGACAC 61.7
    hbae5 R GGAAGAGCATGGCAATCACAACAAG 59.0
    hbbe1.1 F TACAACGCCGCCGCCATCAT 63.2
    hbbe1.1 R GAGGAGACCACAACAGCGAGGAA 61.6
    hbbe1.2 F TACAACGCCGCCGCCATCAT 63.2
    hbbe1.2 R AGAGACCACAACAGCGAGGAACTT 60.5
    hbbe2 F ACTCAGCGGTACTTCGGCAGTT 61.3
    hbbe2 R GCAGCAATCACGATAGTCAGGCAAT 59.8
    hbaa1 F GACAAGGACAAGGCTGCCGTAAG 61.0
    hbaa1 R TGAGCGAAGTAGGTCTTGGTCTGAG 60.3
    hbaa2 F ACGATCTTGTTGGTGGATTGAGCAA 58.7
    hbaa2 R AGCAGATACAGCCGACAGGAACT 60.5
    hbba1 F TCCACCTTCGGCAACCTGTCA 61.8
    hbba1 R AGCGGACACGACCACAGACA 61.9
    hbba2 F GCCTTTACGCCAAACTAAGCAAACT 58.1
    hbba2 R GCCACGACGACAGACAGGAATT 60.1
    下载: 导出CSV

    表  2   黄河高原鳅血红蛋白基因家族成员理化性质

    Table  2   Physical and chemical properties of hemoglobin gene family members in Triplophysa pappenheimi

    基因
    Gene
    CDS长度
    CDS length
    (bp)
    氨基酸数目
    Number of
    amino acid
    分子质量
    Molecular
    weight (ku)
    理论等电点
    Theoretical
    isoelectric point
    脂溶性指数
    Aliphatic
    index
    亲水性平均值
    Grand average of
    hydropathicity
    不稳定性指数
    Instability
    index
    hbae1 444 143 15615.17 9.16 103.71 0.064 23.46
    hbae3 528 175 19411.69 6.37 112.57 0.245 24.81
    hbae4 426 141 15971.74 9.58 94.89 0.040 39.90
    hbae5 432 143 15796.51 9.04 101.68 0.121 24.35
    hbbe1.1 444 147 16417.94 6.96 96.87 0.083 16.63
    hbbe1.2 444 147 16426.95 7.01 96.87 0.085 96.87
    hbbe2 444 147 16786.39 6.90 94.22 –0.073 11.10
    hbaa1 432 143 15552.07 8.89 99.65 0.129 21.91
    hbaa2 432 143 15556.10 9.33 100.98 0.097 16.03
    hbba1 444 147 16113.64 7.70 88.98 0.086 20.53
    hbba2 444 147 16181.79 8.89 99.52 0.134 9.90
    下载: 导出CSV

    表  3   黄河高原鳅血红蛋白基因家族成员二级结构和亚细胞定位

    Table  3   Secondary structure and subcellular localization of hemoglobin gene family members in Triplophysa pappenheimi

    基因
    Gene
    α-螺旋
    Alpha-helix (%)
    β-折叠
    Beta-turn (%)
    无规卷曲
    Random coil (%)
    延伸链
    Extended strand (%)
    亚细胞定位
    Subcellular localization
    hbae1 66.43 4.20 19.58 9.79 Mitochondrion
    hbae3 58.29 6.29 22.29 13.14 Cytoplasm, Extracell
    hbae4 63.83 4.26 22.70 9.22 Cytoplasm
    hbae5 65.73 4.20 21.68 8.39 Mitochondrion
    hbbe1.1 65.31 6.80 17.69 10.20 Mitochondrion
    hbbe1.2 70.75 6.12 17.01 6.12 Mitochondrion
    hbbe2 63.27 3.40 24.49 8.84 Mitochondrion
    hbaa1 67.13 4.90 16.78 11.19 Mitochondrion
    hbaa2 61.54 4.20 25.87 8.39 Mitochondrion
    hbba1 73.47 6.80 14.97 4.76 Mitochondrion
    hbba2 70.07 6.12 19.05 4.76 Mitochondrion
    下载: 导出CSV

    表  4   基于转录组数据的血红蛋白基因在不同组织中的相对表达量

    Table  4   The relative expression of hemoglobin genes in different tissues based on RNA-seq data

    组织Tissue时间Time (h)基因相对表达量 Relative expression of gene
    hbae4hbbe1hbbe2hbaa1hbaa2hbba1hbba2
    肝脏Liver00.060.220.63447.34150.68447.7848.50
    120.070.220.52481.01204.95676.5470.12
    240.120.350.501333.31554.801336.64177.50
    960.160.310.331192.53448.941109.33113.44
    1680.130.300.62835.76354.51871.5871.14
    鳃Gill00.350.510.042253.08957.552217.45264.53
    120.582.1410.164603.072103.876968.25776.39
    241.261.344.466961.283612.676321.26996.59
    961.132.000.588994.123892.078518.22857.03
    1680.731.751.823693.781528.614072.82285.90
    血液Blood023.7742.620.40163575.3464712.88189796.5820481.05
    1215.1196.023.22187315.8082866.30331174.6639786.25
    2413.0116.849.21102190.1943706.21127182.4816116.07
    9621.5458.551.13222169.6987838.61226260.8524392.49
    16829.1286.1917.44201510.7666919.83217962.2419762.22
    下载: 导出CSV
  • [1] 陈宜瑜. 中国动物志·硬骨鱼纲·鲤形目(中卷) [M]. 北京: 科学出版社, 1998: 389-393.]

    Chen Y Y. Fauna Sinica, Osteichthyes, Cypriniformes Ⅱ [M]. Beijing: Science Press, 1998: 389-393. [

    [2] 朱松泉. 中国条鳅志 [M]. 南京: 江苏科学技术出版社, 1989: 101-102.]

    Zhu S Q. Chinese Nemacheilinae [M]. Nanjing: Jiangsu Science and Technology Press, 1989: 101-102. [

    [3] 曹亮, 张鹗, 臧春鑫, 等. 通过红色名录评估研究中国内陆鱼类受威胁现状及其成因 [J]. 生物多样性, 2016, 24(5): 598-609.] doi: 10.17520/biods.2015331

    Cao L, Zhang E, Zang C X, et al. Evaluating the status of China’s continental fish and analyzing their causes of endangerment through the red list assessment [J]. Biodiversity Science, 2016, 24(5): 598-609. [ doi: 10.17520/biods.2015331

    [4]

    Sidell B D. When bad things happen to good fish: the loss of hemoglobin and myoglobin expression in Antarctic icefishes [J]. Journal of Experimental Biology, 2006, 209(10): 1791-1802. doi: 10.1242/jeb.02091

    [5]

    Villota E D, Carmona M G, Rubio J J, et al. Equality of the in vivo and in vitro oxygen-binding capacity of haemoglobin in patients with severe respiratory disease [J]. British Journal of Anaesthesia, 1981, 53(12): 1325-1328. doi: 10.1093/bja/53.12.1325

    [6]

    Weber R E, Vinogradov S N. Nonvertebrate hemoglobins: functions and molecular adaptations [J]. Physiological Reviews, 2001, 81(2): 569-628. doi: 10.1152/physrev.2001.81.2.569

    [7] 齐小琼, 王艇. 植物血红蛋白基因家族的分子进化研究 [J]. 湖北农业科学, 2014, 53(18): 4447-4456.]

    Qi X Q, Wang T. Molecular evolution of hemoglobin families [J]. Hubei Agricultural Science, 2014, 53(18): 4447-4456. [

    [8] 冯蕾, 谭之京, 胡鸢雷, 等. 含血红蛋白基因(VHb)的铜绿假单胞菌重组菌的构建及鼠李糖脂表达条件的研究 [J]. 工业微生物, 2012, 42(3): 78-82.]

    Feng L, Tan Z J, Hu Y L, et al. Construction of recombinant strain containing vitreoscilla hemoglobin gene (VHb) of Pseudomonas aeruginosa and its rhamnolipid expressing conditions [J]. Industrial Microbiology, 2012, 42(3): 78-82. [

    [9] 孟琳, 陈良标. CRISPR/Cas9技术敲除hbae1.1基因对斑马鱼血红蛋白生成的影响 [J]. 生物学杂志, 2021, 38(5): 12-16.]

    Meng L, Chen L B. The effect of hbae1.1 gene on hemoglobin production in zebrafish knocked out by CRISPR/Cas9 [J]. Journal of Biology, 2021, 38(5): 12-16. [

    [10]

    Opazo J C, Butts G T, Nery M F, et al. Whole-genome duplication and the functional diversification of teleost fish hemoglobins [J]. Molecular Biology and Evolution, 2013, 30(1): 140-153. doi: 10.1093/molbev/mss212

    [11]

    Roesner A, Hankeln T, Burmester T. Hypoxia induces a complex response of globin expression in zebrafish (Danio rerio) [J]. Journal of Experimental Biology, 2006, 209(11): 2129-2137. doi: 10.1242/jeb.02243

    [12]

    Roesner A, Mitz S A, Hankeln T, et al. Globins and hypoxia adaptation in the goldfish, Carassius auratus [J]. The FEBS Journal, 2008, 275(14): 3633-3643. doi: 10.1111/j.1742-4658.2008.06508.x

    [13] 吕琳莉, 李朝霞, 崔崇雨. 高原河流溶解氧变化规律研究 [J]. 环境科学与技术, 2018, 41(7): 133-140.]

    Lü L L, Li Z X, Cui C Y. Study on the variation of dissolved oxygen in the plateau river [J]. Environmental Science& Technology, 2018, 41(7): 133-140. [

    [14] 刘壮, 魏峣, 陈强, 等. 长江黄河上游地表水溶解氧时空分布特征研究-以四川省为例 [J]. 环境保护科学, 2023, 49(3): 74-80.]

    Liu Z, Wei Y, Chen Q, et al. Spatial and temporal distribution characteristics of dissolved oxygen in surface water in the upper reaches of the Yangtze River and the Yellow River-Taking Sichuan Province as an example [J]. Environmental Protection Science, 2023, 49(3): 74-80. [

    [15]

    Shang E H H, Wu R S S. Aquatic hypoxia is a teratogen and affects fish embryonic development [J]. Environmental Science & Technology, 2004, 38(18): 4763-4767.

    [16]

    Wu R S S, Zhou B S, Randall D J, et al. Aquatic hypoxia is an endocrine disruptor and impairs fish reproduction [J]. Environmental Science & Technology, 2003, 37(6): 1137-1141.

    [17]

    Sollid J, Nilsson G E. Plasticity of respiratory structures-adaptive remodeling of fish gills induced by ambient oxygen and temperature [J]. Respiratory Physiology & Neurobiology, 2006, 154(1-2): 241-251.

    [18]

    Cossins A R, Crawford D L. Fish as models for environmental genomics [J]. Nature Reviews Genetics, 2005, 6(4): 324-333. doi: 10.1038/nrg1590

    [19]

    Zhu C D, Wang Z H, Yan B. Strategies for hypoxia adaptation in fish species: A review [J]. Journal of Comparative Physiology B, Biochemical, Systemic, and Environmental Physiology, 2013, 183(8): 1005-1013. doi: 10.1007/s00360-013-0762-3

    [20]

    Qi D L, Chao Y, Wu R R, et al. Transcriptome analysis provides insights into the adaptive responses to hypoxia of a schizothoracine fish (Gymnocypris eckloni) [J]. Frontiers in Physiology, 2018(9): 1326. doi: 10.3389/fphys.2018.01326

    [21]

    Tzaneva V, Bailey S, Perry S F. The interactive effects of hypoxemia, hyperoxia, and temperature on the gill morphology of goldfish (Carassius auratus) [J]. American Journal of Physiology Regulatory, Integrative and Comparative Physiology, 2011, 300(6): R1344-1351. doi: 10.1152/ajpregu.00530.2010

    [22]

    Vuori K A M, Soitamo A, Vuorinen P J, et al. Baltic salmon (Salmo salar) yolk-sac fry mortality is associated with disturbances in the function of hypoxia-inducible transcription factor (HIF-1α) and consecutive gene expression [J]. Aquatic Toxicology, 2004, 68(4): 301-313. doi: 10.1016/j.aquatox.2004.03.019

    [23]

    Øijordsbakken M. Effects of hypoxia on angiogenesis and angiogenic factors in crucian carp brain [D]. Oslo: University of Oslo, 2007: 35-47.

    [24]

    Xia M Z, Chao Y, Jia J L, et al. Changes of hemoglobin expression in response to hypoxia in a Tibetan schizothoracine fish, Schizopygopsis pylzovi [J]. Journal of Comparative Physiology B, 2016, 186(8): 1033-1043. doi: 10.1007/s00360-016-1013-1

    [25]

    Lei Y, Yang L, Jiang H, et al. Recent genome duplications facilitate the phenotypic diversity of Hb repertoire in the cyprinidae [J]. Science China Life Sciences, 2021, 64(7): 1149-1164.

    [26]

    Chen C, Chen H, Zhang Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data [J]. Molecular Plant, 2020, 13(8): 1194-1202. doi: 10.1016/j.molp.2020.06.009

    [27]

    Thompson J D, Gibson T J, Plewniak F, et al. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools [J]. Nucleic Acids Research, 1997, 25(24): 4876-4882. doi: 10.1093/nar/25.24.4876

    [28]

    Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11 [J]. Molecular Biology and Evolution, 2021, 38(7): 3022-3027. doi: 10.1093/molbev/msab120

    [29]

    Livak K J, Schmittgen T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method [J]. Methods, 2001, 25(4): 402-408. doi: 10.1006/meth.2001.1262

    [30]

    Bhattacharya M, Hota A, Kar A, et al. In silico structural and functional modelling of antifreeze protein (AFP) sequences of Ocean pout (Zoarces americanus, Bloch & Schneider 1801) [J]. Journal of Genetic Engineering and Biotechnology, 2018, 16(2): 721-730. doi: 10.1016/j.jgeb.2018.08.004

    [31]

    Ganis J J, Hsia N, Trompouki E, et al. Zebrafish globin switching occurs in two developmental stages and is controlled by the LCR [J]. Developmental Biology, 2012, 366(2): 185-194. doi: 10.1016/j.ydbio.2012.03.021

    [32] 王发艳. 花斑裸鲤血红蛋白基因组结构及其功能分化的时空表达调控 [D]. 西宁: 青海大学, 2022: 13-33.]

    Wang F Y. Genomic structure and spatiotemporal expression of functional differentiation of hemoglobin in Gymnocypris eckloni [D]. Xining: Qinghai University, 2022: 13-33. [

    [33] 王发艳, 刘丹, 高强等. 花斑裸鲤胚胎/仔鱼型血红蛋白基因家族成员鉴定与时序表达研究 [J]. 水生生物学报, 2022, 46(5): 718-724.]

    Wang F Y, Liu D, Gao Q, et al. The spatio temporal expression of embryo/larval hemoglobin gene in Gymnocypris eckloni [J]. Acta Hydrobiologica Sinica, 2022, 46(5): 718-724. [

    [34]

    Chen J S C, Wang T Y, Tzeng T D, et al. Evidence for positive selection in the TLR9 gene of teleosts [J]. Fish & Shellfish Immunology, 2008, 24(2): 234-242. doi: 10.1016/j.fsi.2007.11.005

    [35]

    Escriva H, Manzon L, Youson J, et al. Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution [J]. Molecular Biology and Evolution, 2002, 19(9): 1440-1450. doi: 10.1093/oxfordjournals.molbev.a004207

    [36]

    Hoffmann F G, Opazo J C, Storz J F. Whole-genome duplications spurred the functional diversification of the globin gene superfamily in vertebrates [J]. Molecular Biology and Evolution, 2012, 29(1): 303-312. doi: 10.1093/molbev/msr207

    [37]

    Shimeld S M, Holland P W H. Vertebrate innovations [J]. Proceedings of the National Academy of Sciences of the United Satats of America, 2000, 97(9): 4449-4452. doi: 10.1073/pnas.97.9.4449

    [38]

    Jaillon O, Aury J M, Brunet F, et al. Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype [J]. Nature, 2004, 431(7011): 946-957. doi: 10.1038/nature03025

    [39]

    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. doi: 10.1002/bies.20293

    [40]

    Quinn N L, Boroevich K A, Lubieniecki K P, et al. Genomic organization and evolution of the Atlantic salmon hemoglobin repertoire [J]. BMC Genomics, 2010(11): 539. doi: 10.1186/1471-2164-11-539

    [41] 武云飞, 康斌, 门强, 等. 西藏鱼类染色体多样性的研究 [J]. 动物学研究, 1999, 20(4): 258-264.]

    Wu Y F, Kang B, Men Q, et al. Chromosome diversity of Tibetan fishes [J]. Zoological Research, 1999, 20(4): 258-264. [

    [42]

    Lei Y, Yang L, Zhou Y, et al. Hb adaptation to hypoxia in high-altitude fishes: Fresh evidence from Schizothoracinae fishes in the Qinghai-Tibetan Plateau [J]. International Journal of Biological Macromolecules, 2021(185): 471-484.

    [43] 何云凌, 吴丽颖, 朱玲玲, 等. 线粒体自噬在低氧适应中的作用 [J]. 生物化学与生物物理进展, 2012, 39(3): 217-223.] doi: 10.3724/SP.J.1206.2011.00191

    He Y L, Wu L Y, Zhu L L, et al. The role of mitophagy in hypoxic adaptation [J]. Progress in Biochemistry and Biophysics, 2012, 39(3): 217-223. [ doi: 10.3724/SP.J.1206.2011.00191

    [44] 陈祺昌, 郑志琴, 刘丹, 等. 硬骨鱼类血红蛋白转换表达的新模式-以黄河裸裂尻鱼胚胎型血红蛋白研究为例 [J]. 水生生物学报, 2020, 44(6): 1199-1207.]

    Chen Q C, Zheng Z Q, Liu D, et al. A new pattern of hemoglobin switching in teleost fish-study of the embryonic hemoglobin in the Schizopygopsis pylzovi [J]. Acta Hydrobiologica Sinica, 2020, 44(6): 1199-1207. [

    [45]

    Wawrowski A, Gerlach F, Hankeln T, et al. Changes of globin expression in the Japanese medaka (Oryzias latipes) in response to acute and chronic hypoxia [J]. Journal of Comparative Physiology B, 2011(181): 199-208. doi: 10.1007/s00360-010-0518-2

  • 期刊类型引用(9)

    1. 侯兴蓉,陈攀,符圣波,毛海俊,林川,黄海. 两种养殖模式对银鼓鱼生长、消化酶及抗氧化性能的影响. 河北渔业. 2025(03): 1-3 . 百度学术
    2. 侯红霞,王禹奥,石兴伟,李川,唐伟斌. 大口黑鲈生长性能测定与分析. 邢台学院学报. 2025(02): 72-81 . 百度学术
    3. 张瑜,李辰钰,王涛,陈立婧,龚小玲. 上海市三大饮用水源地水库鳙肌肉特性比较. 上海海洋大学学报. 2024(03): 692-701 . 百度学术
    4. 丁兰,徐胜南,周聃,刘梅,李万祈,邹松保. 水产养殖固体废弃物有机肥对土壤肥效和蔬菜品质的影响. 中国农学通报. 2024(19): 102-108 . 百度学术
    5. 韩多彩,张静静,王亚冰,乔光德,王倩,彭士明. 大黄鱼不同抗流能力群体间能量代谢差异比较分析. 海洋渔业. 2024(05): 608-615 . 百度学术
    6. 张先炳,秦依雯,杨威,李耕,胡玉鹏,杨胜发,胡江,李文杰. 鱼类生物能量学模型研究及应用进展. 南方水产科学. 2024(06): 53-61 . 百度学术
    7. 顾玲玲,冯广朋,韩志强. 暗纹东方鲀运动训练后游泳能力和呼吸代谢的变化. 水生生物学报. 2023(12): 2011-2017 . 本站查看
    8. 汤保贵,陈秀灿,胡长圣,戴佳玥,王学锋,张静. 流速对黑鲷幼鱼游泳行为及运动生理的影响. 水生生物学报. 2023(12): 1993-2002 . 本站查看
    9. 蒲德成,郑吉澍,汪政希,翟旭亮,余米,李佩原,李脉,刘世仙. 不同地理大鲵子三代养殖群体肌肉营养成分比较. 饲料研究. 2023(23): 56-61 . 百度学术

    其他类型引用(7)

图(6)  /  表(4)
计量
  • 文章访问数:  89
  • HTML全文浏览量:  13
  • PDF下载量:  8
  • 被引次数: 16
出版历程
  • 收稿日期:  2024-08-28
  • 修回日期:  2024-10-31
  • 网络出版日期:  2024-12-11
  • 刊出日期:  2025-05-14

目录

/

返回文章
返回