黄鳝生长差异基因发掘及其调控机制初步研究

霍欢欢, 刘瑜, 周秋白, 郭枫, 隗黎丽, 彭墨, 张燕萍, 陈文静

霍欢欢, 刘瑜, 周秋白, 郭枫, 隗黎丽, 彭墨, 张燕萍, 陈文静. 黄鳝生长差异基因发掘及其调控机制初步研究[J]. 南方水产科学, 2020, 16(1): 1-8. DOI: 10.12131/20190176
引用本文: 霍欢欢, 刘瑜, 周秋白, 郭枫, 隗黎丽, 彭墨, 张燕萍, 陈文静. 黄鳝生长差异基因发掘及其调控机制初步研究[J]. 南方水产科学, 2020, 16(1): 1-8. DOI: 10.12131/20190176
HUO Huanhuan, LIU Yu, ZHOU Qiubai, GUO Feng, WEI Lili, PENG Mo, ZHANG Yanping, CHEN Wenjing. Primary study on differentially expressed genes screening of Monopterus albus and their regulation mechanism[J]. South China Fisheries Science, 2020, 16(1): 1-8. DOI: 10.12131/20190176
Citation: HUO Huanhuan, LIU Yu, ZHOU Qiubai, GUO Feng, WEI Lili, PENG Mo, ZHANG Yanping, CHEN Wenjing. Primary study on differentially expressed genes screening of Monopterus albus and their regulation mechanism[J]. South China Fisheries Science, 2020, 16(1): 1-8. DOI: 10.12131/20190176

黄鳝生长差异基因发掘及其调控机制初步研究

基金项目: 国家特色淡水鱼产业体系建设专项资金 (CARS-46);国家自然科学基金项目 (31360641)
详细信息
    作者简介:

    霍欢欢 (1986—),男,博士,讲师,从事鱼类营养及免疫研究。E-mail: hhh16430@163.com

    通讯作者:

    周秋白 (1965—),男,博士,教授,从事鱼类营养研究。E-mail: zhouqiubai@163.com

  • 中图分类号: S 917.4

Primary study on differentially expressed genes screening of Monopterus albus and their regulation mechanism

  • 摘要: 为揭示黄鳝 (Monopterus albus) 生长相关基因的调控机制,对相同亲本具有显著性生长差异的黄鳝肝脏进行了转录组测序分析。结果显示,转录组测序共得到19 149个基因,其中差异基因598个,差异基因中有303个基因显著上调,295个显著下调。KEGG通路分析发现,598个差异基因分属262条通路中,其中有38条通路显著富集。GO功能注释发现,与生长相关的差异基因有7个,分别为col1α1、nkx6.1、nnosplexina4、igfbp1、pcgf1和h3.3,这些基因表达水平的变化可能对黄鳝神经、内分泌和消化系统的发育及生理活动产生了调节作用从而影响了黄鳝的生长。结合KEGG通路分析发现,col1α1所在的利什曼病通路和nnos所在的精氨酸和脯氨酸代谢通路富集显著,说明其对黄鳝生长具有重要调节作用。
    Abstract: In order to reveal the regulatory mechanism of growth related genes in eels (Monopterus albus), we carried out a transcriptome sequencing analysis on the liver of the eels that had significant growth difference with the same parent. A total of 19 149 genes were obtained by transcriptome sequencing, among which 598 were differentially expressed (303 and 295 were significantly up-regulated and down-regulated, respectively). KEGG pathway analysis shows that 598 differentially expressed genes belonged to 262 pathways, among which 38 pathways are significantly enriched. The GO functional annotation reveals that the seven differentially related genes were col1α1, nkx6.1, nnos, plexina4, igfbp1, pcgf1and h3.3, and variation in the expression levels of these genes had regulated the development and physiological activities of the nervous system, endocrine system and digestive system of the eels, which then affected their growth. Combined with KEGG pathway analysis, it is found that the leishmaniasis pathway, in which col1α1 is located, and the arginine and proline metabolic pathway, in which nnos is located, are significantly enriched, revealing their important influence on the growth of eels.
  • 鲇 (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   差异基因火山和统计图

    红色代表上调,蓝色代表下调,黑色代表无差异,每个点代表一个基因

    Figure  1.   Volcano-plots (a) and statistics (b) of differentially expressed genes

    Up-regulated genes are in red; down-regulated ones are in blue; black represents no difference; each dot represents a gene.

    图  2   差异表达基因GO功能分类图

    1. 单生物体过程;2. 细胞过程;3. 代谢过程;4. 生物调节;5. 生物过程调节;6. 定位;7. 刺激应答;8. 发展过程;9. 多细胞生物体过程;10. 生物体细胞组成或起源;11. 信号传导;12. 正调节生物过程;13. 负调节生物过程;14;生长;15. 运动;16. 免疫系统过程;17. 行为;18多生物体过程;19. 节律过程;20. 细胞;21. 细胞部分;22. 膜;23. 膜部分;24. 细胞器;25. 复杂大分子;26. 细胞器部分;27. 细胞外区域;28. 细胞外区域部分;29. 膜封闭腔;30. 细胞连接;31. 细胞外模型;32. 超分子纤维;33. 突触;34. 突触部分;35. 结合;36. 催化活性;37. 运输活性;38. 分子功能调节;39. 信号传感活性;40. 结构分子活性;41. 分子传感活性;42. 核苷酸结合转录因子活性

    Figure  2.   GO functional classification map of differentially expressed genes

    1. Single-organism process; 2. Cellular process; 3. Metabolic process; 4. Biological regulation; 5. Regulation of biological regulation; 6. Localization; 7. Response to stimulus; 8. Developmental process; 9. Multicellular organismal process; 10. Cell compent organization or biogenesis; 11. Signaling; 12. Positive regulation of biological process; 13. Negative regulation of biological process; 14. Gowth; 15. Locomotion; 16. Immune system process; 17. Behavior; 18. Multi-organism process; 19. Rhythmic process; 20. Cell; 21. Cell part; 22. Membrane; 23. Membrane part; 24. Organelle; 25. Macromolecular complex; 26. Organelle part; 27. Extracellular region; 28. Extracellular region part; 29. Membrane-enclosed lumen; 30. Cell junction; 31. Extracellular matrix; 32. Supramolecular fiber; 33. Synapse; 34. Synapse part; 35. Binding; 36. Catalytic activity; 37. Transporter activity; 38. Molecular function regulation; 39. Signal transducer activity; 40. Structural molecule activity; 41. Molecular transducer activity; 42. Nucleic acid binding transcription factor activity

    图  3   差异表达基因qPCR验证图

    Figure  3.   Validation of differentially expressed genes by qPCR

    表  1   样品信息

    Table  1   Information of samples

    样品Sample体质量Body mass/g全长Total length/cm
    实验组Treatment group AEG1 42.27 35.8
    AEG2 44.10 36.5
    AEG3 35.78 34.6
    对照组Control group ACG1 6.30 19.8
    ACG2 6.61 18.9
    ACG3 6.12 20.2
    下载: 导出CSV

    表  2   荧光定量PCR引物信息

    Table  2   Information of primers used for qPCR

    基因GeneqPCR引物序列 (5'–3')Primer sequence of qPCR退火温度Annealing temperature/℃片段大小Fragment size/bp
    col1α1 F:AGTTGTTTGCGGACCGAGAT 60.0 110
    R:GCAATCTGGCATTTCCTCACA 59.2
    nkx6.1 F:GGACAAAGATGGGAAACGAAA 56.7 96
    R:GCCAGGTATTTGGTCTGTTCA 58.2
    nnos F:CTATCAGTCTGGATGCCACAAC 58.8 115
    R:CAGAGCCCAACAGAAACATTAG 57.3
    plexina4 F:TGCTGAGAACCCTGAGTGGATA 60.6 159
    R:TAGCATTTGCGGTTGTCTTCAT 58.9
    pcgf1 F:CAGCCCTTACTCAACCTCAAA 57.9 167
    R:GCATCTGGCACAGCATCTACG 61.7
    igfbp1 F:CAGAGAGCCTTGGAAAAGATTG 57.3 171
    R:CTTGCCGTTCCAGGAGTGT 59.9
    h3.3 F:ATTTTGAGTTGCGGCGATTA 56.4 181
    R:GTAACGATGGGGCTTCTTCAC 59.0
    18S F:GTGGAGCGATTTGTCTGGTTA 57.8 162
    R:CGGACATCTAAGGGCATCAC 57.7
    下载: 导出CSV

    表  3   测序结果统计

    Table  3   Statistics of sequencing results

    样品Sample过滤后ReadsFiltered ReadsQ20/%Q30/%GC含量GC content/%匹配率Matching ratio/%
    对照组Control groupACG14120116898.0394.7546.7875.48
    ACG24104770498.0494.7746.6175.48
    ACG34092654098.0294.7848.1275.98
    实验组Treatment groupAEG14100432498.0994.9646.8774.32
    AEG24098404698.1295.0047.0974.44
    AEG34099936698.0794.8847.3375.21
    下载: 导出CSV

    表  4   生长相关差异表达基因

    Table  4   Growth-related differentially expressed genes

    基因名称Gene name基因IDGene ID差异倍数Fold change所属通路Belonged pathway通路IDPathway ID
    col1α1 109951101 −2.5 利什曼病 Ko05140
    nkx6.1 109952563 −3.3 青少年成熟性糖尿病 Ko04950
    nnos 109959109 4.0 精氨酸和脯氨酸代谢 Ko00330
    plexina4 109959815 3.4 细胞黏附分子 Ko04514
    igfbp1 109961253 2.3 p53信号通路 Ko04115
    pcgf1 109972045 3.3 干细胞潜能调节通路 Ko04550
    h3.3 109974907 2.2 癌症的转录失调 Ko05202
    下载: 导出CSV
  • [1] 周秋白, 朱长生, 吴华东, 等. 饲料中不同脂肪源对黄鳝生长和组织中脂肪酸含量的影响[J]. 水生生物学报, 2011, 35(2): 246-255.
    [2]

    VRTÍLEK M, ŽÁK J, POLAČIK M, et al. Rapid growth and large body size in annual fish populations are compromised by density-dependent regulation[J]. J Fish Biol, 2019, 95(2): 673-678. doi: 10.1111/jfb.14052

    [3]

    GAO Z, LUO W, LIU H, et al. Transcriptome analysis and SSR/SNP markers information of the blunt snout bream (Megalobrama amblycephala)[J]. PLoS One, 2012, 7(8): 1-10.

    [4]

    RIBAS L, ROBLEDO D, GÓMEZTATO A, et al. Comprehensive transcriptomic analysis of the process of gonadal sex differentiation in the turbot (Scophthalmus maximus)[J]. Mol Cell Endocrinol, 2016, 15(422): 132-149.

    [5]

    ROBLEDO D, HERMIDA M, RUBIOLO J A, et al. Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production[J]. Comp Biochem Phys D, 2016, 3(21): 41-55.

    [6]

    YIN F, GAO Q, TANG B, et al. Transcriptome and analysis on the complement and coagulation cascades pathway of large yellow croaker (Larimichthys crocea) to ciliate ectoparasite Cryptocaryon irritans infection[J]. Fish Shellfish Immun, 2016, 50: 127-141.

    [7]

    DU X, LI Y, LI D, et al. Transcriptome profiling of spleen provides insights into the antiviral mechanism in Schizothorax prenanti, after poly (I: C) challenge[J]. Fish Shellfish Immun, 2017, 62: 13-23.

    [8]

    MCCARTHY I D, CARTER C G, HOULIHAN D F. The effect of feeding hierarchy on individual variability in daily feeding of rainbow trout, Oncorhynchus mykiss (Walbaum)[J]. J Fish Biol, 1992, 41(2): 257-263. doi: 10.1111/j.1095-8649.1992.tb02655.x

    [9] 杨帆, 张世萍, 韩凯佳, 等. 投喂频率对黄鳝幼鱼摄食、生长及饵料利用效率的影响[J]. 淡水渔业, 2011, 41(3): 50-54. doi: 10.3969/j.issn.1000-6907.2011.03.009
    [10]

    XIAO S, HUANG Q, KUNG A W C. Genetics of osteoporosis in Chinese[J]. Int J Rheum Dis, 2010, 11(4): 359-365.

    [11]

    GUPTA S K, WESOLOWSKA-ANDERSEN A, RINGGAARD A K, et al. NKX6.1, induced pluripotent stem cell reporter lines for isolation and analysis of functionally relevant neuronal and pancreas populations[J]. Stem Cell Res, 2018, 29: 220-231. doi: 10.1016/j.scr.2018.04.010

    [12]

    YU S Y, ZHANG M, LUO J, et al. Curcumin ameliorates memory deficits via neuronal nitric oxide synthase in aged mice[J]. Prog Neuro-Psychoph, 2013, 45(4): 47-53.

    [13] 周立斌, 白茹, 马细兰, 等. 饲料维生素C对美国红鱼 (Sciaenops ocellatus) 组织NO含量、NOS活力和nNOS基因表达的影响[J]. 海洋与湖沼, 2014, 45(5): 1122-1126. doi: 10.11693/hyhz20140700193
    [14]

    GUTEKUNST C A, STEWART E N, FRANZ C K, et al. Plexina4 distribution in the adult rat spinal cord and dorsal root ganglia[J]. J Chem Neuroanat, 2012, 44(1): 1-13. doi: 10.1016/j.jchemneu.2012.03.002

    [15]

    SMOLKIN T, NIR-ZVI I, DUVSHANI N, et al. Complexes of plexin-A4 and plexin-D1 convey semaphorin-3C signals to induce cytoskeletal collapse in the absence of neuropilins[J]. J Cell Sci, 2018, 131(9): 5733-5783.

    [16]

    TAN C, LU N N, WANG C K, et al. Endothelium-derived semaphorin 3G regulates hippocampal synaptic structure and plasticity via neuropilin-2/PlexinA4[J]. Neuron, 2019, 101(5): 920-937. doi: 10.1016/j.neuron.2018.12.036

    [17]

    EMERSON S E, LIGHT S E, EBERT A M. Neuronal expression patterns of the PlexinA family during zebrafish development[J]. Gene Expr Patterns, 2018, 27: 56-66. doi: 10.1016/j.gep.2017.10.007

    [18]

    WANG J, LI Y C, DENG M, et al. Serum insulin-like growth factor 1 and its binding protein 3 as prognostic factors for the incidence, progression, and outcome of hepatocellular carcinoma: a systematic review and meta-analysis[J]. Oncotarget, 2017, 8(46): 81098-81108.

    [19]

    SHIMIZU M, DICKHOFF W W. Circulating insulin-like growth factor binding proteins in fish: their identities and physiological regulation[J]. Gen Comp Endocr, 2017, 252: 150-161. doi: 10.1016/j.ygcen.2017.08.002

    [20]

    FENG X, LIN J, XING S, et al. Higher IGFBP-1 to IGF-1 serum ratio predicts unfavourable survival in patients with nasopharyngeal carcinoma[J]. BMC Cancer, 2017, 17(1): 90-100. doi: 10.1186/s12885-017-3068-0

    [21]

    YUKIYA Y, AKIHIRO A, NOBUHIKO E. Clarifying the impact of polycomb complex component disruption in human cancers[J]. Mol Cancer Res, 2014, 12(4): 479-484. doi: 10.1158/1541-7786.MCR-13-0596

    [22]

    SHUKLA S, ZHAO Q, YING W, et al. Small molecule inhibitors of ring1B-Bmi1 E3 ligase target polycomb repressive complex 1 activity and regulate cell proliferation[J]. Cancer Res, 2015, 75(15): 3520-3520.

    [23] 关珊丽, 王健, 谢忆, 等. 多梳蛋白CBX家族在哺乳动物中结构与功能多样性的研究进展[J]. 中国科学:生命科学, 2017, 47(8): 847-862.
    [24]

    LIN C J, KOH F, WONG P, et al. Hira-mediated H3.3 incorporation is required for DNA replication and ribosomal RNA transcription in the mouse zygote[J]. Dev Cell, 2014, 30(3): 268-279. doi: 10.1016/j.devcel.2014.06.022

    [25]

    WEN D C, BANASZYNSKI L A, ROSENWAKS Z, et al. H3.3 replacement facilitates epigenetic reprogramming of donor nuclei in somatic cell nuclear transfer embryos[J]. Nucleus, 2014, 5(5): 369-375. doi: 10.4161/nucl.36231

    [26]

    WEN D C, BANASZYNSKI L A, LIU Y, et al. Histone variant H3.3 is an essential maternal factor for oocyte reprogramming[J]. P Natl Acad Sci USA, 2014, 111(20): 7325-7330. doi: 10.1073/pnas.1406389111

    [27]

    DELBARRE E, IVANAUSKIENE K, KUNTZIGER T, et al. DAXX-dependent supply of soluble (H3.3 –H4) dimers to PML bodies pending deposition into chromatin[J]. Genome Res, 2013, 23(3): 440-451. doi: 10.1101/gr.142703.112

    [28] 黄星卫, 程香荣, 王楠, 等. 组蛋白H3变体H3.3及其在细胞重编程中的作用[J]. 遗传, 2018, 40(3): 186-196.
  • 期刊类型引用(1)

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

    其他类型引用(2)

图(3)  /  表(4)
计量
  • 文章访问数:  5097
  • HTML全文浏览量:  2764
  • PDF下载量:  82
  • 被引次数: 3
出版历程
  • 收稿日期:  2019-09-08
  • 修回日期:  2019-11-02
  • 网络出版日期:  2019-11-29
  • 刊出日期:  2020-02-04

目录

/

返回文章
返回