克氏原螯虾(
Procambarus clarkii)属于甲壳动物,是我国第四大淡水养殖虾类,截至2024年,我国克氏原螯虾年产量已达345×10
4 t,养殖面积达到2.03×10⁶ hm
2[1]。克氏原螯虾原产于北美洲中南部及墨西哥东北部,目前在我国主要分布于湖北、安徽、湖南、江苏、江西等省份,其为雌雄异体、体内受精的卵生动物,繁殖过程具有抱卵孵化、雌虾抱卵护幼的特点,通常一年产卵一次,繁殖活动多在洞穴内完成
[2]。在国内,克氏原螯虾长期采用自繁自养的养殖模式,该模式已引发系统性种质衰退
[3],其中低温不耐受特性尤为突出,克氏原螯虾在5 ℃的环境中可存活约1周,低于5 ℃时存活时间显著缩短
[4],这限制了其在北方地区的养殖推广,但目前克氏原螯虾低温适应相关研究还未见报道。水温是调控甲壳动物(含克氏原螯虾)地理分布、生长、蜕壳及繁殖的核心环境因子
[5]。不同甲壳动物对低温的生理阈值差异较大
[4,6],且其低温耐受能力具有生理可塑性,这种可塑性的遗传基础已在多种水产动物中得到广泛验证
[7]。通过持续环境胁迫驯化筛选,可快速积累适应性遗传变异,实现存活率等抗逆性状的代际增益
[8-10]。
肠道α-淀粉酶和脂肪酶可调控能量利用与重分配,在甲壳动物环境胁迫响应中发挥着关键作用
[11-12],甲壳动物的遗传背景会影响其在急性低温胁迫后的肠道淀粉酶活性。长期低温选择压力下,耐低温品系可能在肠道消化酶活力上产生遗传性定向选择,这种由环境耐受性引起的肠道消化酶活力变化已在日本沼虾(
Macrobrachium nipponense)
[13]、六斑刺鯚(
Diodon holocanthus)
[14]等多种水生物种中有所报道。肝胰腺作为甲壳动物的代谢枢纽,负责碳水化合物、脂肪和蛋白质的消化、吸收、储存与转运,低温环境下机体可通过调节代谢通路来平衡能量供需
[15-16],从而提高甲壳动物抵御低温胁迫、维持生理平衡的能力。
综上所述,甲壳动物低温耐受能力的生理可塑性依赖于消化功能调控与肝胰腺代谢重编程等机制。本课题组通过多代选择育种,获得了F
6代克氏原螯虾耐低温品系(cold resistant strain,RC),其可在6 ℃环境下正常生长摄食,且存活率、活动能力与抗氧化防御能力均显著高于普通养殖群体
[17-18]。但低温适应并非仅依赖抗氧化防御体系,急性低温胁迫下消化与代谢的能量重分配是决定克氏原螯虾低温耐受能力的关键。目前,低温适应能力不同的克氏原螯虾品系在急性低温胁迫下的消化酶活性调节及肝胰腺代谢重编程研究尚未见报道。为此,本研究以克氏原螯虾RC品系(F
6代)和普通养殖群体为试验对象,经短时低温胁迫后,比较两群体肠道α-淀粉酶、脂肪酶活性及肝胰腺代谢谱差异,以期揭示克氏原螯虾低温适应的消化与代谢调控机制,为其低温适应理论及耐低温种质培育提供科学依据。
1 材料与方法
1.1 材 料
采用本课题组在辽宁省沈阳市基地所选育的克氏原螯虾耐寒品系F6代(RC)与当地普通养殖克氏原螯虾种群(CT)为试验对象,活体样本均由沈阳玉兰农业有限公司捕捞后保鲜空运至无锡渔业学院,挑选活力较好、无明显外伤的个体作为试验虾。
1.2 样本处理
试验用克氏原螯虾在正式试验前均暂养2周,以消除运输及环境变化带来的应激影响,暂养期间水体温度维持在(26±1) ℃,溶解氧保持在(6.4±0.2) mg/L。试验开始前,对RC组与CT组虾个体进行标记。
为诱发可比的低温生理响应,在不引起大规模死亡的前提下,采用统一的短时低温胁迫作为试验处理条件。将标记后的RC组与CT组各50只虾个体同时置于6 ℃水体中,进行3 min的急性低温胁迫处理。低温胁迫处理期间,持续观察并记录2组克氏原螯虾的活动状态及是否出现异常反应或死亡。低温处理结束后,立即从RC组与CT组中各随机抽选3尾雄虾和3尾雌虾解剖,采集肝胰腺组织,于无RNA酶的冻存管中,置于液氮中低温过夜冻存后,转移至-80 ℃冰箱中保存。采集试验虾肠道组织,轻轻挤出内容物后将空肠剪开,先用载玻片刮去表层残留的食糜,再刮取黏膜,采集的黏膜样品立即放入液氮中速冻,-80 ℃冰箱冻存,以备后续测定。
1.3 肠道消化酶活性测定
取肠黏膜组织,采用淀粉酶和脂肪酶试剂盒(南京建成生化工程研究所)测定α-淀粉酶(α-amylase,AMS)和脂肪酶(lipase,LPS)活性,具体测定方法按照试剂盒的说明书进行。
1.4 肝胰腺组织样本处理
取约0.5 g肝胰腺组织样品置于研钵中,分次加入液氮,边加边研磨,使肝胰腺组织呈均匀粉末状,其间要确保样品始终处于冷冻状态不发生融化。取100 mg研磨后的肝胰腺组织样本置于EP管中,加入体积分数80%甲醇水溶液500 μL进行涡旋振荡,振荡后的样品冰浴静置5 min,4 ℃下15 000 r/min离心20 min,取上清液。上清液加质谱级水稀释至甲醇含量为53%,再次进行离心(条件同上),并收集上清液。每个样本取20 μL上清液混合,即得质控样本(quality control,QC)。
1.5 肝胰腺代谢组学分析
代谢组学检测试验在北京诺禾致源科技股份有限公司完成,采用液相色谱-串联质谱(liquid chromatograph-mass spectrometer,LC-MS/MS)技术对代谢物样本进行定性与定量分析,具体试验条件如下:液相色谱采用HPLC C18色谱柱(规格:100 mm×2.1 mm,粒径1.9 μm),分别设置阳离子与阴离子检测模式的流动相体系:阳离子模式流动相A为0.1%甲酸的水溶液,流动相B为纯甲醇;阴离子模式流动相A为pH 9.0的5 mmol/L醋酸铵水溶液,流动相B为纯甲醇。色谱分离的共性参数设定为:流速0.2 mL/min,柱温40 ℃,进样量5 μL。梯度洗脱程序如下:0~1.5 min,流动相B保持2%;1.5~3 min,流动相B由2%线性升至100%;3~10 min,维持100%流动相B进行洗脱;10~10.1 min,流动相B由100%快速降至2%;10.1~12 min,保持2%流动相B以平衡色谱柱,确保后续样本检测的稳定性。采用Q Exactive™ HF-X高分辨质谱仪采集质谱数据,结合上述色谱条件进行代谢物离子信号的采集与记录,为后续代谢组学分析提供基础数据支持。
1.6 肝胰腺代谢组学数据处理与统计分析
将原始数据导入Compound Discoverer 3.1 (CD 3.1)搜库软件中进行处理,对代谢物的保留时间与质荷比等参数进行简单筛选后,设置保留时间偏差0.2 min和质量偏差5×10
-6对样品进行峰对齐,使鉴定结果更加准确,随后进一步设置质量偏差5 ×10
-6、信号强度偏差30%、信噪比3、最小信号强度及加和离子等信息进行峰提取,并对峰面积进行定量,再对目标离子进行整合,通过分子离子峰和碎片离子进行代谢物分子式的预测,并在 mzCloud(
https://www.mzcloud.org/)、mzVault和Masslist数据库中进行比对,最后得到代谢物的定性和相对定量结果。为评估代谢组数据的稳定性与可靠性,分别在阳离子和阴离子模式下进行主成分分析(principal component analysis,PCA),以不同处理组样本在PCA图中的分离程度判断组间代谢差异,以质量控制(quality control,QC)样本在PCA图中的聚集程度验证仪器分析过程的稳定性及数据的重复性;使用KEGG数据库、HMDB数据库和Lipidmaps数据库对鉴定到的代谢物进行注释及分类;采用
t检验分析各代谢物在两组间的差异显著性。以变量投影重要度(variable importance in the projection,VIP)>1.0、差异倍数(fold change,FC)>1.2或<0.83且
P<0.05为标准对差异代谢物(differentially expressed metabolites,DEMs)进行筛选,对筛选到的差异代谢物进行聚类分析、偏最小二乘法判别分析(partial least squares discriminant analysis,PLS-DA)、KEGG 通路富集分析与代谢通路拓扑分析。
2 结果与分析
2.1 2组克氏原螯虾的耐低温性能及肠道酶活性差异
低温胁迫处理3 min内,2组克氏原螯虾均未出现即刻死亡现象。胁迫处理3 min后观察发现,RC组中仍有38%的个体可维持正常爬行等生理活动;而CT组所有个体均出现身体翻倒、活动暂停等典型的低温胁迫响应现象。由
图1可知,在急性低温胁迫处理条件下,RC组个体的肠道α-淀粉酶和脂肪酶活性均较CT组显著降低(
P<0.05),提示急性低温胁迫对RC组克氏原螯虾肠道消化酶活性的抑制作用更为明显。
2.2 肝胰腺代谢组数据质控
图2表明,试验数据整体质量良好。在阴离子模式(
图2-A)下,主成分1(PC1)和主成分2(PC2)的贡献率分别为44.48%和16.41%,累计贡献率达60.89%,能够较好地反映样本间的主要代谢差异;RC和CT样本在PCA图上明显分离,表明2组样本在代谢物组成上存在差异;QC 3个样本紧密聚集于坐标图中心区域,说明仪器分析过程稳定、数据重复性良好。在阳离子模式(
图2-B)下,PC1和PC2的贡献率分别为21.39%和18.81%,累计贡献率为40.20%;RC和CT样本表现出明显的分离趋势;QC 3个样本依旧集中分布,这进一步验证了检测系统的稳定性与试验结果的可靠性。综合来看,2种离子模式下的PCA结果均表明,QC样本聚集良好,RC和CT样本分离明显,表明试验仪器运行稳定、数据质量可靠,所得代谢组数据可用于后续差异代谢物筛选及生物通路分析。
2.3 肝胰腺代谢物鉴定及注释
对2组虾肝胰腺代谢物进行PCA,结果(
图3-A,B)显示,相较于CT组,RC组代谢物更为集中,且2组样本在阴阳离子模式下代谢物均可以较好地聚类区分。对2种群克氏原螯虾肝胰腺中的代谢物进行注释(
图3-C),共鉴定出1 216种代谢物,其中阳离子模式下848种,阴离子模式下368种;将代谢物分别进行KEGG、HMDB及Lipidmaps数据库注释,结果阳离子模式下分别注释到229,340和113个代谢物,阴离子模式下分别注释到113,156和105个代谢物。
2.4 肝胰腺代谢物分类
Lipidmaps 脂质分类结果(
表1)显示,克氏原螯虾肝胰腺脂质代谢物可归为7大类,其中脂肪酰类38种,占33.93%;甘油磷脂类29种,占25.89%;甾醇类19种,占16.96%;聚酮类化合物16种,占14.29%;异戊烯脂质7种,占6.25%;鞘脂2种,占1.79%;甘油脂1种,占0.89%。整体来看,克氏原螯虾肝胰腺脂质组成以脂肪酰类和甘油磷脂类为主。
KEGG通路注释分类结果(
表1)显示,代谢相关通路注释的代谢物最多,有253种,占比66.23%;其次为生物体系统,通路有79种代谢物,占比20.68%;环境信息处理相关通路注释到36种代谢物,占比9.42%;细胞过程与遗传信息处理通路注释到的代谢物数量较少,各有7种,占比均为1.83%。KEGG分类结果表明,注释代谢物主要富集于代谢通路,其次为生物体系统与环境信息处理相关通路,细胞过程和遗传信息处理通路的代谢物较少。
HMDB注释分类结果(
表1)显示,代谢物主要分为脂类和类脂质分子、有机酸及其衍生物、有机杂环化合物等3大类,其中脂类和类脂质分子类代谢物数量最多(92种),占比最大;其次为有机酸及其衍生物(78种)和有机杂环化合物(64种)类代谢物;苯类化合物,苯丙素类和聚酮类化合物,核苷、核苷酸及其类似物类代谢物分别为28,21和17种;有机含氧化合物、有机含氮化合物、生物碱及其衍生物类代谢物数量较少;烃类衍生物类代谢物仅有1种。
综合HMDB、KEGG和Lipidmaps等3个数据库的注释分类结果可知,本研究中的代谢物主要集中于脂质代谢及能量代谢等生物学过程,其中脂质和类脂质分子在数量与功能层面均占据主导地位。
2.5 肝胰腺DEMs筛选及PLS‐DA分析验证
以VIP>1.0,FC>1.2或<0.83且
P<0.05为标准,对DEMs进行筛选,阳离子模式下共获得16种上调DEMs,阴离子模式下共获得6种上调DEMs;阳离子模式下共获得71种下调DEMs,阴离子模式下共获得25种下调DEMs。由
图4可见,2种群克氏原螯虾的代谢物差异明显,且在阳离子模式下DEMs可以准确聚类至2个不同的种群。
进一步对DEMs进行PLS-DA分析,结果见图5。由图5可知,阳离子模式下模型拟合优度(R2)为0.91,显著高于阴离子模式的0.74;阳离子模式和阴离子模式的交叉验证预测率(Q2)分别为0.22和0.24,水平相近。综合聚类结果和模型参数表现可知,阳离子模式下RC和CT区分度较好且模型的拟合优度更好,故后续差异DEMs筛选及通路富集分析均基于阳离子模式结果进行。
2.6 肝胰腺差异代谢物功能富集分析
将阳离子模式下的87个DEMs导入Metabo Analyst 6.0数据库进行代谢通路富集分析,结果如
图6所示。以
P<0.05,即-lg
P >1.3、拓扑影响值>0.05作为重要代谢通路的筛选条件,共筛选出1条显著代谢通路,即甘油酸磷脂代谢通路,有2种DEMs出现在该通路中,分别是磷酸胆碱(phosphocholine,P-Cho)和甘油磷酸胆碱(glycerophosphocholine,GPC)。提取甘油酸磷脂代谢通路子网络图,
图6展示了克氏原螯虾在急性低温胁迫下DEMs所涉及的甘油酸磷脂代谢(glycerophospholipid metabolism)子通路,通路涉及的关键酶包括胆碱激酶(choline kinase,CHKA)、磷酸胆碱胞苷转移酶(phosphate cytidylyltransferase1A,PCYT1A)、神经元型磷脂酶(patatin-like phospholipase domain-containing protein6,PNPLA6)、磷脂酰转移酶样蛋白3(phospholipase A/acyltransferase3,PLAAT3)、甘油磷酸胆碱磷酸二酯酶1(glycerophospho-choline phosphodiesterase 1,GPCPD1)、胆碱磷酸乙醇胺磷酸转移酶1(choline/ethanolamine phosphotransferase 1,CEPT1)和磷酸酶1(phosphatase 1,PHOSPHO1),与CT组相比,RC组克氏原螯虾体内P-Cho和GPC均显著下调。
3 讨 论
几乎所有生物体都有代偿性调节机制,能够通过生理和代谢变化来适应外界环境的温度变化
[19]。为了能在低温环境下生存,甲壳动物会通过代谢抑制
[20]、调整膜脂组成
[14]等一系列的措施来应对低温环境胁迫,代谢抑制可减慢机体的代谢速率,降低能量消耗;调整膜脂组成可防止因低温造成的膜相变和渗透性破坏,增加糖原储备。然而,生理代谢调节等机制能够适应的温度范围有限。变温动物在应对低温胁迫时,普遍采取降低整体代谢率的策略,以减少能量消耗,将有限的能量重新分配至维持基本生理过程
[21]。在本研究中,RC组在急性低温下肠道α-淀粉酶与脂肪酶活性显著降低,这并非简单的酶活性因低温而被动降低,而是一种适应性的主动调节。消化是一个高能耗的生理过程,涉及多种消化酶的合成、分泌以及营养物质的吸收和同化,需要消耗大量ATP
[22]。在低温环境下,动物活动减少,摄食需求降低,此时维持高水平的消化功能会造成能量的浪费。RC品系通过主动抑制消化吸收功能,从而实现对整体代谢率的迅速下调,该结果与前人对罗非鱼的研究报道
[23]一致。
本研究代谢组学分析进一步揭示,RC组与CT组可基于代谢物谱清晰区分,表明连续多世代的选育使耐低温品系形成了稳定且可遗传的代谢特征。差异代谢物富集结果显示,甘油磷脂代谢通路发生显著富集,通路中磷酸胆碱和甘油磷酸胆碱含量均显著下调。该变化由CHKA、PCYT1A及磷脂酶家族PNPLA6、PLAAT3、GPCPD1等关键酶协同调控。低温胁迫通常会导致细胞膜脂质发生相变,从液晶态转为凝胶态,从而影响膜蛋白的功能和物质的跨膜运输
[24],甘油磷脂是构成细胞膜最基本的骨架结构,其组成和动态平衡对于维持细胞膜在低温下的流动性和完整性至关重要
[25],水生动物的脂肪酸代谢对温度变化很敏感
[26],膜脂肪酸去饱和常被认为是水生生物适应低温的重要机制,对维持膜的流动性、酶活力和细胞的正常功能具有重要作用
[27]。冷应激导致细胞中脂肪酸的组成发生改变,不饱和脂肪酸所占比例上升,这有利于维持细胞膜的流动性
[28-30],例如当吉富罗非鱼在15 ℃环境中时,其肝组织中的不饱和脂肪酸比例显著升高
[31]。甘油磷脂是细胞双层膜的主要组成成分,本研究代谢组的结果显示,低温胁迫下RC组甘油磷脂代谢通路被抑制,其中磷酸胆碱和甘油磷酸胆碱2种代谢物均较CT组下调。在代谢网络中,甘油磷酸二酯酶催化代谢反应生成的GPC与P-Cho,是后续胆碱代谢途径的直接前体物质,二者一旦进入该途径便被磷酸胆碱-磷酸胆碱胞苷转移酶活化为胞苷二磷酸(cytidine diphosphate‑choline,CDP)胆碱
[32],随后与二酰甘油结合生成磷脂酰胆碱(phosphatidylcholine,PC)
[32-33]。PC分子在sn-2位具有高度的酰基选择性,优先接受不饱和脂肪酸形成PC(16:0/PUFA)等不饱和磷脂
[34]。因而,当组织或细胞对多不饱和脂肪酸(polyunsaturated fatty acid,PUFA)的需求显著上升时,外源或内源 PUFA会被迅速转运至内质网并酰化到PC的sn-2位,驱动 CDP-胆碱途径的通量加速,以满足膜流动性、信号分子合成及脂肪酸储存的需求。该过程导致游离的P-Cho与GPC被快速消耗,因而在代谢物中检测到的自由P-Cho/GPC浓度下降,而不饱和磷脂及其下游PUFA产物则相对增多,这正是代谢流向从胆碱前体向不饱和脂肪酸-PC复合体转移的直接表现。RC品系通过这种高效的代谢重塑,在低温胁迫下迅速调整膜脂组成,从而维护细胞膜结构与功能的稳定。
中国水产科学研究院淡水渔业研究中心基本科研业务费项目(2025JBFR03)
江苏省种业振兴揭榜挂帅项目(JBGS [2021] 123)
盘锦市“带土移植”项目(DTYZ2024003)
中央公益性科研机构基础研究基金项目(2023TD39)