基于LC-MS技术的海、淡水养殖刀鲚卵巢的代谢组学比较分析

高淑芳, 张金鹏, 施永海, 袁新程, 刘其根

高淑芳, 张金鹏, 施永海, 袁新程, 刘其根. 基于LC-MS技术的海、淡水养殖刀鲚卵巢的代谢组学比较分析[J]. 南方水产科学, 2022, 18(3): 68-75. DOI: 10.12131/20210185
引用本文: 高淑芳, 张金鹏, 施永海, 袁新程, 刘其根. 基于LC-MS技术的海、淡水养殖刀鲚卵巢的代谢组学比较分析[J]. 南方水产科学, 2022, 18(3): 68-75. DOI: 10.12131/20210185
GAO Shufang, ZHANG Jinpeng, SHI Yonghai, YUAN Xincheng, LIU Qigen. Metabonomics analysis of ovaries of Coilia nasus in seawater and freshwater based on liquid chromatography-mass spectrometry[J]. South China Fisheries Science, 2022, 18(3): 68-75. DOI: 10.12131/20210185
Citation: GAO Shufang, ZHANG Jinpeng, SHI Yonghai, YUAN Xincheng, LIU Qigen. Metabonomics analysis of ovaries of Coilia nasus in seawater and freshwater based on liquid chromatography-mass spectrometry[J]. South China Fisheries Science, 2022, 18(3): 68-75. DOI: 10.12131/20210185

基于LC-MS技术的海、淡水养殖刀鲚卵巢的代谢组学比较分析

基金项目: 农业部鱼类营养与环境生态研究中心 (A1-3201-19-300204);上海市科委重点攻关项目(17391900300)
详细信息
    作者简介:

    高淑芳 (1994—),女,硕士研究生,研究方向为刀鲚分子生物学。E-mail: sfgao101@163.com

    通讯作者:

    刘其根 (1965—),男,教授,从事水域生态学研究。E-mail: qgliu@shou.edu.cn

  • 中图分类号: S 965.227

Metabonomics analysis of ovaries of Coilia nasus in seawater and freshwater based on liquid chromatography-mass spectrometry

  • 摘要: 为探究海、淡水养殖环境对刀鲚 (Coilia nasus) 卵巢发育的影响,采用非靶向代谢组学的方法,测定了海、淡水环境下刀鲚卵巢中代谢产物的差异情况,并与基因组百科全书 (KEGG) 数据库进行比对,找出相对应的代谢通路并分析其原因。结果表明,海水组和淡水组样品共鉴定出47种差异代谢物 (P<0.05、FC>1、VIP>1),与海水组相比,淡水组表达差异倍数最明显的为碳环血氧烷A2 (Carbocyclic thromboxane A2)、半乳糖神经酰胺(Galactosyl ceramide),差异倍数分别为10.40、2.78倍;与海水组相比,淡水组卵巢组织内皮质醇升高了1.61倍;对47种差异代谢物进行KEGG分析发现,变化显著的通路有氨酰-tRNA的生物合成和嘧啶代谢通路 (P<0.05),皮质醇、氨酰-tRNA的生物合成通路、嘧啶代谢通路和鞘磷脂代谢通路可能与刀鲚生殖洄游过程中卵巢发育有关。
    Abstract: In order to clarify the effects of seawater and freshwater on the ovary development of Coilia nasus, we analyzed their differences by using non-targeted metabolomics, and compared with database of KEGG directly to find out the corresponding metabolic pathways, then analyzed its causes. The results show that a total of 47 metabolites had significant difference between the two groups (P<0.05, FC>1, VIP>1). Compared with the seawater group, the most significant differences in the expression were carbocyclic thromboxane A2, Galactosylceramide, and their differences were 10.40 and 2.78 times, respectively. The cortisol in the ovarian tissue of the freshwater group increased by 1.61 times. According to the analysis on KEGG metabolic pathways of 47 different metabolites, the biosynthesis of aminoacyl-tRNA and pyridine metabolic pathways changed in the seawater and freshwater environments significantly (P<0.05). The biosynthesis pathway of cortisol, aminoyl-tRNA, pyrimidine metabolism pathway and sphingo-lipid metabolism pathway may be related to the ovarian development during the reproductive migration of C. nasus.
  • 鲇 (Silurus asotus) 属鲇形目、鲇科、鲇属,广泛分布于我国江河、湖泊,是洞庭湖流域主要鱼类之一[1],性凶猛,肉食性,多以鱼虾为食,具有适应性强、生命力旺、生长发育快、抗病能力强、肉质细腻鲜嫩、肌间刺少等特点[2],已被广泛推广养殖。在鲇形目鱼类繁殖生物学方面,国内已对南方鲇 (S. meriordinalis)[3-4]、怀头鲇 (S. soldatovi)[5]、兰州鲇 (S. lanzhouensis)[6]、洞庭湖鲶 (S. asotus)[7]、豹纹脂身鲶 (Pterygoplichthys pardalis)[8]、中华鲱鲇 (Clupisoma sinensis)[9]开展了一系列研究。

    鱼类的繁殖能力对其资源变化有一定影响[10]。沅江属长江流域洞庭湖支流水系,是湖南省第二大河流,全长1 033 km,起源于贵州,在湖南汉寿县汇入洞庭湖[11]。凌津滩水电站地处湖南省常德市桃源县,是沅江干流最下游的一个梯级,与上游五强溪水电站相距47.5 km。鱼类的生物学特性由遗传和环境共同作用决定[12],为适应环境,同一物种不同地理种群间生物学特性存在一定的差异[13-16]。目前尚未见对沅江干流水域鲇种群生物学相关的研究报道,选取凌津滩库区为本研究的取样及生态环境调查点,对沅江凌津滩库区鲇的生长与繁殖特征进行研究,一方面可为野生鲇的定种、品种改良以及人工繁殖等提供基础理论资料,另一方面通过探讨鲇对大坝隔断的生态响应,可为其种质资源保护与合理利用提供科学依据,并为凌津滩库区制定鱼类种群调控政策提供理论支撑。

    2018年5月至2020年4月共进行了24次采样,每月15日左右在湖南省常德市桃源县的凌津滩库区兴隆街河段采集鲇样本约30尾,共采集样本738尾。采集的鲇在鲜活状态下,采用带水增氧方式运输至环洞庭湖水产健康养殖及加工湖南省重点实验室进行分析。

    体长、体高等指标采用游标卡尺测量,精确到0.01 cm,体质量、去内脏质量 (空壳质量) 和性腺质量等采用电子天平称量,精确到0.01 g。怀卵量采用质量取样法,每年5—8月从IV期卵巢中称1.00 g卵,在Motic显微操作系统下用计数器对分散成熟的全部卵粒计数。参照文献[17]计算鲇绝对繁殖力和相对繁殖力,同时,随机选取每尾鱼的100粒卵子测定卵粒直径。

    体长与体质量的关系参考郜星晨等[17]的方法利用幂函数进行拟合,肥满度 (K) 的计算参照刘艳超[18]采用Fulton法,计算公式如下:

    $$ W=a L^b $$ (1)
    $$ K=100 \times W / L^3 $$ (2)

    式中:WLab分别代表体质量、体长、生长的条件因子、异速生长因子。

    性腺发育期根据黄小林等[19]通过目测比对区分。性别鉴定通过鱼体外形特征结合性腺解剖检查进行,并统计雌、雄个体数目计算性比,通过 $\chi^2 $ 检验性比差异显著性。雌性体长以10 mm为间距分组,统计各体长组性成熟个体比例,进行Logistic曲线拟合,并推算雌性鲇群体50% 性成熟体长,公式为:

    $$ P_i=\left[1+\mathrm{e}^{-{k}(L_ {i}-L _{50})}\right]^{-1} $$ (3)

    式中:Pi为各体长区段性成熟个体百分比;k为成熟度增长系数;Li为各体长区段中间值;L50为50% 性成熟体长。

    绝对繁殖力、相对繁殖力、性成熟系数 计算公式如下:

    $$ F=N \times W_{\mathrm{O}} $$ (4)
    $$ F_{\mathrm{W}}=F / W_{\mathrm{N}} $$ (5)
    $$ \mathrm{GSI}=100 \times\left(W_{\mathrm{O}} / W_{\mathrm{N}}\right) $$ (6)

    式中:F为绝对繁殖力;FW为相对繁殖力;GSI为性成熟系数 (%) ;N为每克卵巢卵数 (粒·g−1);WO 为性腺质量;WN为空壳质量。

    采用Excel 2016与SPSS 21.0软件进行数据整理、图表绘制和统计分析。

    外观上,可以根据鲇尾鳍分叉的深度来鉴别性别,雄性尾鳍分叉深而明显,超过尾鳍2/3以上,雌性尾鳍分叉浅,不到尾鳍的1/3 (图 1)。雌、雄胸鳍第一硬刺后缘不同,雌性较为光滑;雄性有粗壮的锯齿,用手摸有割手感。腹部感官不同,雌性膨大具有弹性,外生殖突短而粗大,生殖孔呈椭圆形且较突出、充血红肿;雄性腹部狭小,外生殖突细长,生殖孔狭小稍长、充血不明显。在繁殖盛期,雌性体型较雄性大,体表光滑多黏液,下腹膨胀柔软,卵巢轮廓隐约可见,挤压腹部有卵粒流出;雄性体表黏液少,躯体较狭长,腹部手感较硬。

    图  1  雌雄鲇尾鳍分叉对比
    Figure  1.  Differences of tail fins between female and male S. asotus

    沅江鲇的性比在不同时期有较大差异,但大多是雌性明显多于雄性。根据对738尾样本的统计,其中雌性404尾、雄性334尾,性比为1.21,通过$\chi^2 $检验,显著大于1 (P<0.05)。不同月份间性比有差异,7、8、10月最高;3、12月次之,高于4月;4月高于1、9月,1、9月高于2、5、11月,6月最低,且差异显著 (P<0.05) (表1)。

    表  1  鲇 1—12 月的性比
    Table  1.  Sex ratio of S. asotus from January to December
    月份
    Month
    雌 Female/尾 雄 Male/尾 雌∶雄 Female∶Male
    1 35 31 1.13c
    2 23 39 0.59b
    3 40 20 2.00e
    4 38 22 1.73d
    5 23 35 0.66b
    6 16 49 0.33a
    7 44 16 2.75f
    8 46 18 2.56f
    9 33 27 1.22c
    10 40 16 2.50f
    11 25 43 0.58b
    12 41 18 2.27e
    注:同列不同上标字母示差异显著(P<0.05)。表2表5同此。 Note: Different superscripts within the same column indicate significant differences (P<0.05). The same case in Table 2 and Table 5.
    下载: 导出CSV 
    | 显示表格

    按照采样月份对测量的738尾鲇的体长、体高、体质量进行统计 (表2),均在11月达到最高,分别为 (39.52±4.81) cm、(5.88±0.76) cm、 (381.60±158.40) g。从体长来看,5—11月除9月偏低外总体呈递增趋势,在11月达到最高,12月开始下降,次年1月最低,1—4月又呈现递增趋势。体高、体质量的变化趋势与体长基本一致。体长、体高、体质量周年变化与繁殖、越冬等生理活动相一致,9月偏低可能由取样数量偏少导致。雌性体长为19.20~44.00 cm,平均 (32.67±4.52) cm,优势体长为35.00~40.00 cm;体质量为125.00~1 042.00 g,平均 (381.78±105.34) g。雄性体长为14.80~37.20 cm,平均 (28.19±5.62) cm,优势体长为30~35 cm;体质量为55.00~514.00 g,平均 (276.43±88.72) g。

    表  2  鲇体长、体高、体质量的周年变化
    Table  2.  Annual changes in body length, body height and body mass of S. asotus
    月份
    Month
    体长
    Body length/cm
    体高
    Body height/cm
    体质量
    Body mass/g
    空壳质量
    Eviscerated mass/g
    范围
    Range
    均值
    Average
    范围
    Range
    均值
    Average
    范围
    Range
    均值
    Average
    范围
    Range
    均值
    Average
    1 23.50~38.50 30.02±4.45b 3.20~5.30 4.17±0.67a 92.00~499.00 229.21±96.08b 87.00~446.00 205.30±86.67b
    2 26.40~39.20 32.29±3.45b 3.30~6.15 4.69±0.52b 128.00~600.00 295.17±97.14c 121.00~529.50 272.00±85.28c
    3 14.80~43.20 34.38±3.74b 3.32~6.23 4.73±0.65b 55.00~691.00 306.97±116.53c 46.00~610.50 277.97±107.46c
    4 25.40~44.00 35.62±3.74b 3.46~7.25 5.03±0.73c 135.00~710.00 375.46±98.72e 119.50~628.50 337.23±88.92e
    5 14.80~39.30 26.71±5.62a 2.80~5.30 4.09±0.66a 55.00~583.00 249.63±124.11b 46.50~449.50 210.10±105.23b
    6 14.50~35.58 30.58±4.32b 3.40~5.60 4.46±0.49b 154.00~524.50 327.00±97.73d 142.00~473.50 290.84±84.29d
    7 27.10~39.40 33.14±3.01b 3.80~5.60 4.63±0.56b 168.00~449.50 288.70±88.70c 112.00~342.50 254.03±76.60c
    8 30.70~40.11 33.95±2.91b 3.60~5.75 4.77±0.56b 209.00~644.00 340.78±105.34e 195.50~564.50 300.72±90.08d
    9 21.52~36.70 26.77±4.09a 1.27~6.30 4.15±0.63a 91.50~583.00 177.59±101.51a 81.50~549.00 158.14±94.64a
    10 25.40~41.20 32.57±3.95b 3.47~8.60 4.93±1.21c 152.50~772.00 327.50±160.20d 174.00~661.50 296.07±139.70d
    11 21.50~49.73 39.52±4.81c 3.70~7.50 5.88±0.76d 184.00~1042.00 381.60±158.40e 124.00~917.50 341.43±137.50e
    12 26.80~47.32 35.07±5.34b 3.35~6.38 4.92±0.82c 149.50~925.00 370.90±182.30e 134.00~872.50 338.06±162.20e
    下载: 导出CSV 
    | 显示表格

    经协方差分析显示,雌雄间体长与体质量关系差异不显著 (F=0.39, P=0.57>0.05),不考虑性别因素,通过回归分析得到738尾鲇体质量和体长的幂函数关系为:W=0.056L2.57 (R=0.97,图2),b值2.57,用t检验分析显示低于3的理论匀速生长值,差异显著 (P<0.05),说明鲇体长按线性长度增加,体质量则按体长的幂次增加,群体生长为异速型。

    图  2  鲇体长-体质量关系
    Figure  2.  Relationship between body length and body mass of S. asotus

    鲇周年的发育分期 (表3) 与性成熟系数相一致。IV期即成熟期为5月,V期即生殖期为6—8月,VI—II期即疲弱期为9月,10—12月为II期即静息期,有部分是从VI期重新向II期过渡,处于I期的数量很少,卵巢中尚未排出的卵粒慢慢退化,并被亲鱼吸收利用。鲇在繁殖季节5—8月,各月性腺发育分期的组成存在一定差异 (图3),其中5月以IV期为主,占75.50%;6月以V期为主,占67.50%;7月以V和VI期为主,分别占52.90%和41.50%;8月以VI期为主,占72.50%。

    表  3  鲇性腺发育状况
    Table  3.  Gonadal development of S. asotus
    月份
    Month
    雄 Male 雌 Female
    发育分期
    Developmental stage
    颜色
    Color
    发育分期
    Developmental stage
    颜色
    Color
    5 IV, IV—V, V 白色,乳白色 IV, V 草绿色
    6 V, V—VI 白色,乳白色 IV, V 草绿色,黄绿色
    7 V, V—VI, VI 白色,淡红色 IV, V, VI 草绿色,黄绿色
    8 V—VI, VI 白色,淡红色 V—VI, VI 草绿色,黄绿色
    9 VI, VI—II 淡红色 VI, VI—II 浅绿色,黄绿色
    10 II 淡红色,浅肉红色 II 浅绿色,红黄色
    11 II, II—III 浅肉红色 II, II—III 红黄色,浅黄色
    12 II, II—III 肉红色 II, II—III 黄色,黄绿色
    1 III 肉红色 III, III—IV 黄绿色,浅绿色
    2 III—IV 深肉红色,白色 III, IV 黄绿色,浅绿色
    3 III, IV 白色,乳白色 III, IV 黄绿色,浅绿色
    4 III, IV 白色,乳白色 III, IV 黄绿色,浅绿色
    下载: 导出CSV 
    | 显示表格
    图  3  鲇卵巢发育分期的月份变化
    Figure  3.  Variation of ovaries developmental stages of S. asotus from Yuanjiang River in different months

    对404尾雌性鲇,体长以10 mm为间距分组,然后对各分组性成熟比例和体长进行Logistic回归拟合,获得性腺成熟度百分比与体长组的关系方程如下 (图4):

    图  4  鲇雌性性成熟比例与体长关系
    Figure  4.  Relationship between body length and sex mature frequency of female S. asotus in Yuanjiang River
    $$ P_i=\left[1+\mathrm{e}^{-0.092(L_{{i}}-258.85)}\right]^{-1}\left(R^2 = 0.98, n = 404\right) $$ (7)

    通过Logistic曲线回归法拟合估算得到,雌性鲇的群体50%性成熟体长为258.85 mm。

    鲇全年平均性成熟系数雌性高于雄性 (P<0.01),雌性平均7.51,雄性平均1.85。鲇雌性性成熟系数5月最高,达15.24,5—9月下降,9月达最低(2.04),9月至次年5月呈上升趋势,鱼类卵巢的成熟系数只在繁殖期达到最大,而周年中鲇性成熟系数仅出现一个峰值,说明鲇为一次产卵型鱼类;雄性性成熟系数也是5月最高(4.24),5—11月下降,11月达最低(0.74),11月至次年5月上升。说明鲇雌、雄性成熟系数的周年变化趋势基本同步,从而保障了鲇的繁殖率。

    鲇雌性性成熟系数平均值从1月一直增加到6月,之后逐渐下降,9月最低,10—12月较低;雄性性成熟系数变化规律与雌性相似 (表4)。鲇群体肥满度平均为0.88,雌性平均0.94,雄性平均0.81,雌性显著高于雄性 (P<0.01)。雌、雄肥满度周年变化呈相同趋势,12月至次年5月逐步上升,5—7月下降,7—10月上升,10—12月下降。雌、雄肥满度均在5月最高,分别为1.37、1.19;在11月最低,分别为0.65、0.56。肥满度周年变化与繁殖、越冬等生理活动相适应。

    表  4  鲇 1—12 月份的性成熟系数
    Table  4.  Gonadosomatic indexes of S. asotus from January to December
    月份
    Month
    体质量
    Body mass/g
    空壳质量
    Eviscerated mass/g
    性腺质量
    Gonad mass/g
    性成熟系数
    Gonadosomatic index
    肥满度
    Fullness

    Female

    Male

    Female

    Male

    Female

    Male

    Female

    Male

    Female

    Male
    1 264.97 193.43 238.56 171.50 16.05 1.82 6.73 1.06 0.79 0.79
    2 327.10 263.24 234.84 232.75 17.97 2.85 7.65 1.22 0.87 0.79
    3 345.19 268.75 306.13 249.80 26.34 4.67 8.60 1.87 0.89 0.68
    4 428.63 322.28 379.25 295.21 30.58 7.42 8.06 2.51 0.92 0.75
    5 264.45 234.81 215.14 206.07 32.79 8.74 15.24 4.24 1.37 1.19
    6 351.66 312.57 297.71 283.97 34.80 11.85 11.69 4.17 1.20 1.03
    7 298.60 286.63 257.36 250.69 24.51 4.08 9.52 1.63 0.83 0.72
    8 367.80 313.60 312.88 288.56 27.59 4.23 8.82 1.47 0.91 0.79
    9 194.47 160.90 171.68 144.52 3.50 1.63 2.04 1.13 0.97 0.84
    10 357.60 298.38 319.83 272.31 12.621 3.20 3.95 1.18 0.99 0.86
    11 421.08 321.36 367.71 315.17 13.66 2.34 3.71 0.74 0.90 0.77
    12 411.69 359.78 365.84 310.28 15.06 3.12 4.11 1.01 0.65 0.56
    下载: 导出CSV 
    | 显示表格

    共分析了5—8月67尾雌性性成熟个体的繁殖力,个体绝对繁殖力最大为125 324 粒,最小为4 259 粒,平均 (29 677±18 174) 粒;个体相对怀卵量最大为287 粒·g−1,最小为17 粒·g−1,平均 (107±65) 粒·g−1。最小的雌鱼个体体长为23.65 cm,体质量为109.52 g,绝对繁殖力 (F) 与体质量 (W) 和体长 (L) 均成直线相关,个体繁殖力随体质量、体长的增加而增加,关系式分别为:

    $$ F=437.36 W{\text{+}}945.09\left(R^2 = 0.93, P < 0.01\right) $$ (8)
    $$ F=1\;205.38 L{\text{−}}21\;285.27\left(R^2 = 0.87, P < 0.01\right) $$ (9)

    鲇平均卵径范围介于794.40~1 908.40 μm,平均 (1 444.50±97.45) μm,优势卵径组为1 200~1 400 μm (图5)。

    图  5  鲇 IV 期卵巢的卵径分布
    Figure  5.  Oocyte diameter-frequency distribution of S. asotus ovary at stage IV

    表3可知,鲇5—8月的性腺发育程度,5月以IV期为主,6月以V期为主,7—8月以VI期为主,结合表5得出,IV—VI期不同性腺发育期其卵径的大小不一样,表现为IV期高于V期高于VI期,独立t检验显示,IV—VI期的性腺发育程度与卵径呈负相关关系,且差异显著 (P<0.05)。

    表  5  鲇繁殖力和卵径
    Table  5.  Fecundity and egg diameter of S. asotus
    月份
    Month
    绝对繁殖力
    Absolute fecundity/粒
    相对繁殖力
    Relative fecundity/粒
    卵径范围
    Egg diameter range/μm
    平均卵径
    Average egg diameter/μm
    557 223266998.5~1 908.41 626.5±101.2a
    667 744228974.2~1 861.61 588.6±66.2b
    742 166164856.3~1 601.61 259.3±133.6c
    831 13199794.4~1 303.71 303.6±88.8c
    下载: 导出CSV 
    | 显示表格

    鲇总性别比例为雌性404尾、雄性334尾,雌∶雄=1.21∶1,通过 $\chi^2 $ 检验,大于 1∶1 的理论比值且差异显著 (P<0.05),与肖智[20]鲇繁殖习性研究结果一致,而与短尾高原鳅 (Triplophysa brevieauda)[18]、宽口裂腹鱼 (Schizothorax eurystomus)[21]等繁殖群体中雌性显著少于雄性不同,雄鱼明显多于雌鱼,可保障受精率,提高繁殖力。鲇雌、雄总比例大于1,但在2、5、6、11月低于1,呈现季节性差异,这种现象也出现在西藏哲古措异尾高原鳅 (T. stewarti)[22]、沙塘鳢 (Odontobutis obscura)[23]中。尤其是鲇6月的性比仅0.33,这可能与雄性个体小满足繁殖所需相关。鲇雌雄周年总比例大于 1∶1 的理论比值,可能由取样数较少导致,因此要进一步得出鲇的性别比例状况,还需加大取样量和采样的随机性;也可能是雌鱼生长快,具有竞争优势所致。雌性鲇生长速度明显快于雄性,因此,鲇雌性明显多于雄性,正好符合人工养殖需要,有利于提高养殖产量。但是雄性数量少,个体又较雌性小,难以保证受精率,对其繁殖极为不利。要解决这一矛盾需加强对鲇性别分化机制的探讨。何福林等[24]在沟鲶 (Ietalurus punetaus) 鱼苗生长初期使用甲基睾丸酮等激素处理,将沟鲶遗传上的雌性转化为生理上的雄性,从而提高沟鲇雄性的数量。对鲇能否利用甲基睾丸酮等激素诱导雌转雄,还有待进一步研究。若可行,便有望通过遗传雌性与转化雄性交配培育出全雌鲇,如同全雄黄颡鱼 (Pelteobagrus fulvidraco)[25]一样应用于养殖生产,提高养殖产量和经济效益。

    高小平等[26]认为对硬骨鱼类而言体质量 (W) 与体长 (L) 为幂函数关系:Wa Lbb的取值范围为2.5~4.0。本实验得到鲇体长与体质量关系式为:W=0.056L2.57R=0.97,b=2.57,在高小平等[26]报道的取值范围之内,但与低等甲壳动物中华绒螯蟹 (Eriocheir sinensis) 的直线线性相关不同[27]。杨丽亚等[28]研究显示鱼类是否处于等速生长可以用幂指数b的大小来判断,当b为3时,则表明鱼类的生长为等速型。鲇b=2.57,与3差异显著 (P<0.05),说明鲇的生长趋于异速生长,与网箱养殖福瑞鲤 (Cyprinus carpio)  (b=2.84)[29]、雅鲁藏布江下弧唇裂腹鱼 (Schizothorax curvilabiatus) (b=2.94)[30]、短尾高原鳅 (b=3.05)[31]、方氏云鳚 (Enedrias fangi) (b=3.10)[32]、鲢 (Hypophthalmichthys molitrix) (b=2.75)[33]等与3接近的研究结果不一致。影响b值变化的因素包括环境、发育阶段、年龄、季节、捕捞强度[34]、性别、海拔、种群[35]、食物、疾病[17]、个体大小[36]等,本实验并未考虑这些因素的影响,还有待进一步研究。

    本研究结果显示,鲇性成熟系数的周年变化表现出明显的周期性,雌性性成熟系数5月最高,5—9月下降,9月最低,9月至次年5月呈上升趋势;雄性性成熟系数也是5月最高,5—11月下降,11月最低,11月至翌年5月上升,雌、雄性成熟系数的周年变化规律基本一致,推测5月下旬至8月上旬为其繁殖期与实际的繁殖期一致,这一结果与郜星晨等[17]对长江宜昌段鳙 (Aristichthys nobilis) 的繁殖研究结果一致。鲇的卵巢大多处于II、III、IV期,在5月基本进入成熟期,部分开始进入排卵期即V期,GSI变化较大[37]。绝对怀卵量在4 259~125 324粒之间变化,差异很大,繁殖力受个体大小、环境因素、成熟期等影响[38-40]。鱼类的产卵类型通常可以通过卵径分布和成熟系数周年变化来判断[41]。本实验中鲇卵径的周年分布中只在5月达到高峰,这与长江宜昌段鳙[17]呈现双峰不同。鲇成熟卵径均值为1.63 mm,远低于长江宜昌段鳙卵径均值(3.75) mm[17],与同为鲇形目的其他鱼类比较,低于卵径均值为1.87 mm的豹纹脂身鲶,高于卵径均值为1.56 mm的光泽黄颡鱼 (P. nitidus)、卵径均值为1.50 mm的瓦氏黄颡鱼 (P. vachelli) 和卵径均值为1.05 mm的革胡子鲶 (Clariidae gariepinus)[8]。鲇的卵径较大,说明卵所含供孵化后仔鱼利用的营养物质较多,保证了孵化成活的机会,有利于提高受精卵的孵化率和仔鱼的成活率。在实验过程中发现,雌鲇产卵后卵巢内有明显的空卵泡,并残留大量处于III期后的未成熟期卵子,这可能是鲇与生俱来的一种正常自然现象,也可能是外来干扰造成的,其产生原因还有待进一步分析。鲇的卵巢发育是同步的,吴士平和李海洋[3]研究发现鲇可第二次产卵和再次进行人工催产。本实验中发现鲇在生殖期结束时,卵巢内残留较多的成熟卵母细胞,这种繁殖特性对鲇的养殖具有重要意义,可在鲇排卵后进行第二次催产,以提高其产卵量,这一现象与吴士平和李海洋[3]的研究结果一致,但是关于催产的时间间隔还有待进一步研究。

  • 图  1   海水组和淡水组的主成分分析 (a) 和正交偏最小二乘判别分析 (b) 得分图

    Figure  1.   PCA (a) and OPLS-DA (b) scores of seawater and freshwater

    图  2   置换检验图

    Figure  2.   Permutation test

    图  3   代谢产物Venn图

    Figure  3.   Venn map of metabolites

    图  4   差异代谢物筛选火山图

    Figure  4.   Volcanic map of different metabolites

    图  5   代谢物的通路富集分析

    Figure  5.   Metabolite pathway enrichment analysis

    表  1   差异代谢物信息表

    Table  1   Differential metabolites information sheet

    差异代谢物   
    Differential metabolite   
    分子式
    Molecular formula
    变量投影重要度
    VIP
    FC
    Fold change
    P变化趋势
    Variation trend
    磷脂酰肌醇 PI [20:4(5Z,8Z,11Z,14Z)/0:0] C29H49O12P 1.58 0.93 0.017 下降
    谷氨酰胺色氨酸 Glutaminyltryptophan C16H20N4O4 1.05 0.96 0.018 下降
    L-异亮氨酸 L-Isoleucine C6H13NO2 1.33 0.96 0.005 下降
    木麻黄6-α-D-葡萄糖苷 Casuarine 6-alpha-D-glucoside C14H25NO10 1.13 0.94 0.043 下降
    甘油一脂 MG(10:0/0:0/0:0) C13H26O4 3.59 0.71 0.025 下降
    皮质醇 Cortisol C21H30O5 3.28 1.61 0.008 上升
    碳环血氧烷A2 Carbocyclic thromboxane A2 C22H36O3 2.88 10.40 0.007 上升
    (±)9-十八碳二烯酸 (±)9-HPODE C18H32O4 2.90 0.58 0.020 下降
    N-棕榈酰蛋氨酸 N-palmitoyl methionine C21H41NO3S 1.86 0.75 0.022 下降
    13-羟基十八酸 13-hydroxyoctadecanoic acid C18H36O3 1.95 0.86 0.005 下降
    咖啡酰环戊醇 Caffeoylcycloartenol C39H56O4 2.07 0.81 0.036 下降
    11,13-二十碳二烯酸 15-OxoEDE C20H34O3 2.51 0.70 0.032 下降
    2,3-二氢苯并呋喃 2,3-dihydrobenzofuran C8H8O 1.19 1.08 0.033 上升
    4-甲酰基吲哚 4-formyl indole C9H7NO 1.34 0.83 0.040 下降
    红花素C Safflomin C C30H30O14 2.04 0.85 0.003 下降
    6-[(2-羧基乙酰基)氧]-3,4,5-三羟基氧烷-2-羧酸6-[(2-carboxyacetyl)oxy]-3,4,5-trihydroxyoxane-2-carboxylic acid C9H12O10 2.07 1.12 0.001 上升
    异戊二烯 Isoputreanine C7H16N2O2 2.33 0.84 0.001 下降
    脯氨酸 L-Proline C5H9NO2 1.80 0.91 0.003 下降
    羟脯氨酸 Hydroxyprolyl-hydroxyproline C10H16N2O5 3.13 0.72 0.007 下降
    5-O-阿魏酰黑精 5-O-Feruloylnigrumin C21H25NO10 1.0 0.95 0.037 下降
    壬二酸 Azelaic acid C9H16O4 1.12 0.96 0.012 下降
    尿苷 Uridine C9H12N2O6 1.35 0.95 0.001 下降
    视黄酯 Retinyl ester C20H30O2 1.24 0.96 0.004 下降
    顺-9,10-环氧硬脂酸 cis-9,10-Epoxystearic acid C18H34O3 1.66 0.92 0.025 下降
    异柠檬酸盐 Isocitrate C6H8O7 1.55 1.12 0.047 上升
    戊二酸 Glutaric acid C5H8O4 2.53 1.38 0.015 上升
    6-脱氧噬菌体胺 (6-脱氧花青) 6-Deoxyfagomine C6H13NO2 1.36 0.94 0.012 下降
    乙酰-L-酪氨酸 Acetyl-L-tyrosine C11H13NO4 1.38 0.87 0.010 下降
    γ-谷氨酰鸟氨酸 Gamma glutamyl ornithine C10H19N3O5 1.87 1.58 0.023 上升
    对茴香酸异戊酯 (异戊基异茴香酸酯) Isoamyl p-anisate C13H18O3 1.76 0.67 0.041 下降
    3-氧十二酸 3-Oxododecanoic acid C12H22O3 1.01 0.92 0.043 下降
    C-2神经酰胺 C-2 Ceramide C20H39NO3 2.21 0.75 0.028 下降
    磷脂酰乙醇胺 PE(15:0/16:1(9Z)) C36H70NO8P 1.42 1.10 0.030 上升
    半乳糖神经酰胺 Galactosyl ceramide (d18:1/14:0) C38H73NO8 1.97 2.78 0.045 上升
    N-棕榈酰甘氨酸 N-Palmitoyl glycine C18H35NO3 1.97 0.86 0.000 下降
    二十碳五烯酸 Eicosapentaenoic acid C20H30O2 1.01 0.98 0.004 下降
    溶血磷脂酰乙醇胺 LysoPE(0:0/20:2(11Z,14Z)) C25H48NO7P 1.33 1.12 0.037 上升
    9-羟基癸酸 9-Hydroxydecanoic acid C10H20O3 1.21 0.92 0.028 下降
    甘油一脂 MG(a-13:0/0:0/0:0)[rac] C16H32O4 1.92 0.79 0.040 下降
    反式-2-十二碳烯二酸 Traumatic acid C12H20O4 1.37 0.86 0.025 下降
    9-氧壬酸 9-Oxo-nonanoic acid C9H16O3 1.15 0.93 0.019 下降
    8-羟基-5,6-辛二烯酸 8-Hydroxy-5,6-octadienoic acid C8H12O3 1.04 0.94 0.016 下降
    (S)-3-磺酸盐 (S)-3-Sulfonatolactate C3H6O6S 2.57 3.05 0.046 上升
    γ-谷氨酰缬氨酸 gamma-Glutamylvaline C10H18N2O5 1.80 0.87 0.008 下降
    穗花牡荆苷 Agnuside C22H26O11 1.41 0.91 0.005 下降
    假尿苷 Pseudouridine C9H12N2O6 1.19 0.93 0.003 下降
    雌三醇7-(6-反式-对-香豆酰基葡萄糖苷)
    Eriodictyol 7-(6-trans-p-coumaroylglucoside)
    C30H28O13 1.44 0.92 0.020 下降
    注:FC表示某差异代谢物在淡水组相对于海水组的表达倍数变化,FC>1表示该代谢物上调,FC<1表示该代谢物下调;P<0.05表示差异显著,P<0.01表示差异极显著。 Note: FC value indicates the change of expression multiple of a differential metabolite in FOV group compared with SOV group. FC>1 indicates that the metabolite is upregulated; FC<1 indicates that the metabolite is down-regulated. P<0.05 indicates significant difference, and P<0.01 indicates extremely significant difference.
    下载: 导出CSV

    表  2   差异代谢产物的KEGG通路富集表

    Table  2   KEGG pathway enrichment table of differential metabolites

    富集的差异代谢产物名称 
    Enriched differential metabolite name 
    富集个数
    Numble
    富集通路
    ID Pathway
    通路描述  
    Pathway desciption  
    P
    半乳糖神经酰胺 Galactosyl ceramide (d18:1/14:0) 1 map00600 鞘脂代谢 0.059
    脯氨酸 L-Proline 1 map00330 精氨酸和脯氨酸代谢 0.176
    异柠檬酸盐 Isocitrate 1 map00020 柠檬酸循环 0.057
    戊二酸 Glutaric acid 1 map00310 赖氨酸降解 0.130
    皮质醇 Cortisol 1 map00140 类固醇激素生物合成 0.207
    尿苷 Uridine;假尿苷 Pseudo uridine 2 map00240 嘧啶代谢 0.011
    异柠檬酸盐 Isocitrate 1 map00630 乙醛酸和二羧酸代谢 0.137
    L-异亮氨酸 L-Isoleucine;脯氨酸 L-Proline 2 map00970 氨酰tRNA生物合成 0.008
    L-异亮氨酸 L-Isoleucine 1 map00290 缬氨酸、亮氨酸和异亮氨酸生物合成 0.064
    二十碳五烯酸 Eicosapentaenoic acid 1 map01040 不饱和脂肪酸生物合成 0.097
    L-异亮氨酸 L-Isoleucine 1 map00280 缬氨酸、亮氨酸和异亮氨酸降解 0.107
    (S)-3-磺酸盐 (S)-3-Sulfonatolactate 1 map00270 半膀氨酸和蛋氨酸代谢 0.141
    下载: 导出CSV
  • [1]

    SU M, DUAN Z, SHI H, et al. The effects of salinity on reproductive development and egg and larvae survival in the spotted scat Scatophagus argus under controlled conditions[J]. Aquac Res, 2019, 50(7): 1782-1794. doi: 10.1111/are.14056

    [2]

    MOHANTY B, GUPTA K, BABU G, et al. Exposure to salinity stress cause ovarian disruption in a stenohaline freshwater teleost, Heteropneustes fossilis (Bloch, 1794)[J]. Aquac Res, 2020, 51(5): 1964-1972. doi: 10.1111/are.14548

    [3]

    PHAM H Q, KJØRSVIK E, NGUYEN A T, et al. Reproductive cycle in female Waigieu seaperch (Psammoperca waigiensis) reared under different salinity levels and the effects of dopamine antagonist on steroid hormone levels[J]. J Exp Mar Biol Ecol, 2010, 383(2): 137-145. doi: 10.1016/j.jembe.2009.12.010

    [4]

    PHAM H Q, NGUYEN A T, KJØRSVIK E, et al. Seasonal reproductive cycle of Waigieu seaperch (Psammoperca waigiensis)[J]. Aquac Res, 2012, 43(6): 815-830. doi: 10.1111/j.1365-2109.2011.02894.x

    [5] 牛景彦, 刘占才. 影响鱼类性腺发育的生态因素研究[J]. 农业与技术, 2016, 36(16): 109.
    [6] 杜学芳. 盐度对凡纳滨对虾繁殖及家系生长、存活的影响[D]. 上海: 上海海洋大学, 2013: 9-17.
    [7] 吴旭干, 赵亚婷, 何杰, 等. 低盐度海水和淡水对中华绒螯蟹性腺发育及交配行为的影响[J]. 动物学杂志, 2013, 48(4): 555-561.
    [8] 吴建辉, 王家启, 戴小杰, 等. 基于概率模型的长江口鱼类空间共现模式分析[J]. 南方水产科学, 2019, 15(1): 1-9. doi: 10.12131/20180112
    [9]

    ZHU G L, WANG L J, TANG W Q, et al. Identification of olfactory receptor genes in the Japanese grenadier anchovy Coilia nasus[J]. Genes Genomics, 2017, 39(5): 521-532. doi: 10.1007/s13258-017-0517-8

    [10]

    DUAN J R, ZHOU Y F, XU D P, et al. Ovary transcriptome profiling of Coilia nasus during spawning migration stages by Illumina sequencing[J]. Mar Genomics, 2015, 21: 17-19.

    [11]

    LI W X, ZOU H, WU S G, et al. Richness and diversity of helminth communities in the JapaneseI grenadier anchovy, Coilia nasus, during its anadromous migration in the Yangtze River, China[J]. J Parasitol, 2012, 98(3): 449-452. doi: 10.1645/GE-2983.1

    [12]

    LI W X, SONG R, WU S G, et al. Seasonal occurrence of helminths in the anadromous fish Coilia nasus (Engraulidae): parasite indecators of fish migratory movements[J]. J Parasitol, 2011, 97(2): 192-196. doi: 10.1645/GE-2621.1

    [13]

    ZHANG H, WU G, XIE P, et al. Role of body size and temporal hydrology in the dietary shifts of shortjaw tapertail anchovy Coi-lia brachygnathus (Actinopterygii, Engraulidae) in a large floodplain lake[J]. Hydrobiologia, 2013, 703(1): 247-256. doi: 10.1007/s10750-012-1370-z

    [14] 代培, 严燕, 朱孝彦, 等. 长江刀鲚国家级水产种质资源保护区(安庆段)刀鲚资源现状[J]. 中国水产科学, 2020, 27(11): 1267-1276.
    [15]

    SHEN H, GU R, XU G, et al. In-depth transcriptome analysis of Coilia ectenes, an important fish resource in the Yangtze River: de novo assembly, gene annotation[J]. Mar Genomics, 2015, 23: 15-17. doi: 10.1016/j.margen.2015.03.002

    [16] 鲜博, 高建操, 徐钢春, 等. 盐度对刀鲚生长、抗氧化应激和渗透压调节能力的影响[J]. 海洋湖沼通报, 2020(2): 152-159.
    [17] 王武. 鱼类增养殖学[M]. 北京: 中国农业出版社, 2000: 193.
    [18]

    XU G, DU F, LI Y, et al. Integrated application of transcripto-mics and metabolomics yields insights into population-asynchronous ovary development in Coilia nasus[J]. Sci Rep, 2016, 6(1): 1-11. doi: 10.1038/s41598-016-0001-8

    [19]

    YIN D, LIN D, YING C, et al. Metabolic mechanisms of Coilia nasus in the natural food intake state during migration[J]. Genomics, 2020, 112(5): 3294-3305. doi: 10.1016/j.ygeno.2020.05.027

    [20]

    ZHAO H J, XU J K, YAN Z H, et al. Microplastics enhance the developmental toxicity of synthetic phenolic antioxidants by disturbing the thyroid function and metabolism in developing zebrafish[J]. Environ Int, 2020, 140: 105750. doi: 10.1016/j.envint.2020.105750

    [21]

    DING J, HUANG Y, LIU S, et al. Toxicological effects of nano-and micro-polystyrene plastics on red tilapia: are larger plastic particles more harmless?[J]. J Hazard Mater, 2020, 396: 122693. doi: 10.1016/j.jhazmat.2020.122693

    [22]

    ZHANG H, LIU Y, ZHOU L, et al. Metabonomic insights into the sperm activation mechanisms in ricefield eel (Monopterus albus)[J]. Genes, 2020, 11(11): 1259. doi: 10.3390/genes11111259

    [23]

    HUANG Y, ZHANG Y, ZHENG J, et al. Metabolic profiles of fish nodavirus infection in vitro: RGNNV induced and exploited cellular fatty acid synthesis for virus infection[J]. Cell Microbiol, 2020, 22(9): e13216.

    [24]

    OZEN G, NOREL X. Prostanoids in the pathophysiology of human coronary artery[J]. Prostag Oth Lipid M, 2017, 133: 20-28.

    [25]

    YAN H, ZHANG M Z, WONG G, et al. Mechanisms of U46619-induced contraction in mouse intrarenal artery[J]. Clin Exp Pharmacol Physiol, 2019, 46(7): 643-651. doi: 10.1111/1440-1681.13087

    [26]

    SMYTH E M. Thromboxane and the thromboxane receptor in cardiovascular disease[J]. Clin Lipidol, 2010, 5(2): 209-219. doi: 10.2217/clp.10.11

    [27]

    GU W, MADRID D M D, YANG G, et al. Unaltered influenza disease outcomes in swine prophylactically treated with α-galactosylceramide[J]. Dev Comp Immunol, 2021, 114: 103843. doi: 10.1016/j.dci.2020.103843

    [28]

    WANG T, LIU F, TIAN G, et al. Lineage/species-specific expansion of the Mx gene family in teleosts: differential expression and modulation of nine Mx genes in rainbow trout Oncorhynchus mykiss[J]. Fish Shellfish Immunol, 2019, 90: 413-430. doi: 10.1016/j.fsi.2019.04.303

    [29]

    HOEHN K, MARIEB E N. Human anatomy & physiology[M]. San Francisco: Benjamin Cummings, 2010: 1064-1094.

    [30]

    SOPINKA N M, CAPELLE P M, SEMENIUK C A D, et al. Glucocorticoids in fish eggs: variation, interactions with the environment, and the potential to shape offspring fitness[J]. Physiol Biochem Zool, 2017, 90(1): 15-33. doi: 10.1086/689994

    [31]

    MOUSA M A, IBRAHIM M G, KORA M F, et al. Experimental studies on the reproduction of the thin-lipped mullet, Liza ramada[J]. Egypt J Aquat Biol Fish, 2018, 22(3): 125-138. doi: 10.21608/ejabf.2018.9455

    [32]

    CUI W, MA A, WANG X, et al. Myo-inositol enhances the low-salinity tolerance of turbot (Scophthalmus maximus) by modula-ting cortisol synthesis[J]. Biochem Biophys Res Commun, 2020, 526(4): 913-919. doi: 10.1016/j.bbrc.2020.04.004

    [33] 胡静, 叶乐, 吴开畅, 等. 急性盐度胁迫对克氏双锯鱼幼鱼血浆皮质醇浓度和 Na+-K+-ATP 酶活性的影响[J]. 南方水产科学, 2016, 12(2): 116-120. doi: 10.3969/j.issn.2095-0780.2016.02.017
    [34]

    SARAVANAN M, RAMESH M, PETKAM R, et al. Influence of environmental salinity and cortisol pretreatment on gill Na+/K+-ATPase activity and survival and growth rates in Cyprinus carpio[J]. Aquac Rep, 2018, 11: 1-7. doi: 10.1016/j.aqrep.2018.04.002

    [35]

    KWON N H, FOX P L, KIM S. Aminoacyl-tRNA synthetases as therapeutic targets[J]. Nat Rev Drug Discov, 2019, 18(8): 629-650. doi: 10.1038/s41573-019-0026-3

    [36]

    GOMEZ M A R, IBBA M. Aminoacyl-tRNA synthetases[J]. RNA, 2020, 26(8): 910-936. doi: 10.1261/rna.071720.119

    [37]

    GARAVITO M F, NARVÁEZ-ORTIZ H Y, ZIMMERMANN B H. Pyrimidine metabolism: dynamic and versatile pathways in pathogens and cellular development[J]. J Genet Genomics, 2015, 42(5): 195-205. doi: 10.1016/j.jgg.2015.04.004

    [38]

    LAVIEU G, SCARLATTI F, SALA G, et al. Regulation of autophagy by sphingosine kinase 1 and its role in cell survival during nutrient starvation[J]. J Biol Chem, 2006, 281(13): 8518-8527. doi: 10.1074/jbc.M506182200

    [39]

    HANNUN Y A, OBEID L M. The ceramide-centric universe of lipid-mediated cell regulation: stress encounters of the lipid kind[J]. J Biol Chem, 2002, 277(29): 25847-25850. doi: 10.1074/jbc.R200008200

  • 期刊类型引用(1)

    1. 刘凤坤,赵吉臣,许敬轩,郑雅之,鲁慧杰,孙育平,阮灼豪,郭慧,杨俊贤,黄文. 耐盐好氧反硝化菌黏质沙雷氏菌HL4的分离鉴定及脱氮性能. 广东海洋大学学报. 2024(04): 38-46 . 百度学术

    其他类型引用(2)

图(5)  /  表(2)
计量
  • 文章访问数:  724
  • HTML全文浏览量:  319
  • PDF下载量:  100
  • 被引次数: 3
出版历程
  • 收稿日期:  2021-06-23
  • 修回日期:  2021-10-08
  • 录用日期:  2021-11-03
  • 网络出版日期:  2021-11-09
  • 刊出日期:  2022-06-04

目录

/

返回文章
返回