Study of impacts of environmental factors on stock and recruitment relationship of the anchovy stock in the Yellow Sea
-
摘要:
鱼类年际资源量的波动可以归因于年间环境条件的变化和该种鱼类亲体数量的变化。文章根据1990~2001年间黄海中南部鳀鱼声学调查评估结果，以及黄海千里岩海区在此年间的表层水温和营养盐统计数据，以Ricker模型(R=αSe-βS)为基础对黄海鳀鱼(Engraulis japonicas)亲体-补充量关系进行了初步研究。结果表明，黄海千里岩水域表层水温、磷酸盐浓度等环境条件因素对补充量有重要影响。
Abstract:Changes in fish year-class strength have been attributed to year-to-year variability in environmental conditions, such as sea surface temperature (SST) and phosphate concentration, and spawning stock biomass(SSB).In this study, we examined the relationship between environmental factors, SSB and recruitment for the anchovy (Engraulis japonicas) stock in the Yellow Sea using a traditional stock-recruitment Ricker model. The recruitment fluctuation of the anchovy stock in the Yellow Sea was successfully reproduced by the Ricker stock-recruitment model with environmental factors. The findings suggest that the impacts of SST and phosphate concentration observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea on the anchovy stock in the Yellow Sea are very important.
-
Keywords:
- stock-recruitment model /
- environmental condition index /
- AIC /
- BIC
-
鱼类的生长与发育和哺乳动物相似，主要受生长激素/胰岛素样生长因子(growth hormone/insulin-like growth factor，GH/IGF)轴调节[1-4]。近年来，已有多种鱼类的GH和IGF基因被克隆，并且有关其重组蛋白的促生长活性也有一定的研究。研究表明，重组GH和IGF都能促进多数鱼类的生长，但是重组GH对虹鳟(Oncorhynchus mykiss)[5]并无显著促生长活性。最近研究表明，重组鲮(Cirrhinus molitorella)GH和IGF-Ⅰ均有促进罗非鱼(Oreochromis niloticus)生长活性，且重组鲮GH对罗非鱼稚鱼的促生长活性较重组鲮IGF-Ⅰ弱[6]，但是重组鲮GH和IGF-Ⅰ促生长作用的分子机制以及它们在鲮体内的相互作用方式并不清楚。
笔者在研究了重组鲮IGF-Ⅰ对鲮GH表达的影响后[7]，为更深一步了解重组鲮GH和IGF-Ⅰ促生长作用的分子机制以及它们在鲮体内的相互作用方式，该试验以我国南方主要养殖鱼类之一鲮为研究对象，在鲮GH和IGF-ⅠcDNA被克隆和体外重组表达的基础上[3, 8]，对鲮IGF-ⅠmRNA进行组织表达分析，并分析重组鲮GH处理后鲮IGF-Ⅰ表达变化情况。
1. 材料与方法
1.1 重组鲮GH
重组鲮GH的表达和纯化参照ZHANG等[8]的方法进行，纯化后的重组鲮GH充分透析，冷冻干燥，溶解在0.7%的生理盐水中待用(pH 7.0)。
1.2 实验动物的处理以及组织的提取
将鲮(500±10 g)从广州市番禺区鱼窝头镇养殖户的池塘中捕捞起后暂养于广东省鲮鱼原种场，待其状态稳定后，按鱼体重120 μg·g-1注射重组鲮GH或0.7%的生理盐水，12 h后取鲮鱼脑、肝和肌肉保存于液氮中，对照组取脑、肝、肌肉、肠、性腺、鳃、皮肤、脾脏、心脏和肾10种组织保存于液氮；并从每条鱼中抽取血液(大约3 mL)4℃冷凝30 min，4℃ 1 500 g离心15 min，取上清保存于-20℃，待检测GH和IGF-Ⅰ含量[9]。
1.3 总RNA的提取与RT-PCR
按Trizol Reagent(Invitrogen)操作说明提取各组织总RNA，用0.5×TBE配置1%琼脂糖凝胶，取2 μL总RNA与2 μL 2×DEPC处理后加样缓冲液混合，65℃加热15 min后，立即于冰上放置2 min，加样电泳，鉴定所提取总RNA的质量，取2 μL总RNA测定其光密度值OD260和OD280，计算总RNA浓度。
进行RT-PCR前必须测RNA浓度，逆转录体系对RNA量有一定要求，该研究按照逆转录试剂盒TakaRa RNA LA PCRTM Kit (AMV) Ver. 1.1说明书要求，所用总RNA量为500 ng。
按逆转录试剂盒TakaRa RNA LA PCRTM Kit (AMV) Ver. 1.1介绍方法进行逆转录，MgCl2 2μL，10×RNA PCR Buffer 1μL，RNase Free H2O 4.25 μL，dNTP Mixture(各10 mM) 1μL，Rnase Inhibitor 0.25 μL，AMV Reverse Transcriptase 0.5 μL, Oligo(dT)-Adaptor Primer 0.5 μL，Total RNA 0.5 μL(500 ng)。50℃45 min，99℃5 min，5℃5 min，1个循环。
1.4 PCR循环数确定及鲮IGF-ⅠmRNA组织表达
首先确定PCR扩增反应条件，即确定PCR扩增反应循环数[10]。用对照组肝cDNA作模板，引物为Actin-F，Actin-R(表 1)，94℃变性30 s，54℃退火30 s，72℃延伸45 s，每PCR循环数为15、20、25、30、35、40、45取4 μL，PCR产物在1.5%的琼脂糖凝胶电泳EB染色分析。循环数确定后利用看家基因Actin做内标，PCR检测鲮IGF-ⅠmRNA各组织表达情况以及重组鲮GH处理后IGF-ⅠmRNA在脑、肝和肌肉的表达，所有PCR反应体系包括1 μL(0.15 μg)cDNA，0.25 μL rTaq酶(TaKaRa公司)，正反向引物各1 μL，4 μL dNTP，5 μL 10×PCR缓冲液，双蒸水37.75 μL，总体积50 μL。PCR反应均94℃变性5 min，35个循环(94℃变性30 s，54℃退火30 s，72℃延伸45 s)，最后72℃延伸10 min。IGF-Ⅰ mRNA检测引物序列如表 1。
表 1 扩增IGF-Ⅰ、Actin的引物序列Table 1. Oligonucleotide primers used to amplify cDNA for mud carp IGF-Ⅰand Actin引物 primer 序列 sequence Actin-F 5′-GTGTTGGCG/ATACAGGTCCTTACG-3′ Actin-R 5′-CAGACTACCTC/GATGAAGATCCTGAC-3′ IGF-Ⅰ-F 5′-ATGGAAAACCAGCGCCTCTTC-3′ IGF-Ⅰ-R 5′-TGCATGTCCTTCTTGAAGCAAG-3′ 1.5 各组织mRNA表达检测结果分析
利用GeneTools软件分析，计算出步骤1.4中所有IGF-Ⅰ和Actin条带的相对质量，再将对应的IGF-Ⅰ与Actin质量相比，得出GH的相对表达量。
1.6 IGF-Ⅰ含量检测
利用GeneTools软件分析，计算出步骤1.4中所有IGF-Ⅰ和Actin条带的相对质量，再将对应的IGF-Ⅰ与Actin质量相比，得出IGF-Ⅰ的相对表达量。
2. 结果
2.1 总RNA提取与RT-PCR
从100 mg鲮组织中提取RNA，经1%琼脂糖凝胶电泳检测，可见清晰的28 s、18 s条带(图 1)，表明总RNA完整性很好，根据测定RNA的OD260和OD280值，计算其比值，比值均在1.9~2.0之间，也说明RNA纯度较高，可用于后续反应。
![]()
图 1 鲮组织总RNA 1%琼脂糖凝胶电泳1. 脑；2. 肝；3. 肾；4. 性腺；5. 肠；6. 鳃；7. 脾脏；8. 心脏；9. 皮肤；10. 肌肉Figure 1. 1% agarose gel electrophoresis of total RNA from tissues of mud carp1. brain; 2. liver; 3. kidney; 4. gonad; 5. intestine; 6. gill; 7. spleen; 8. heart; 9. skin; 10. muscle2.2 PCR循环数的确定
首先确定PCR反应条件，即确定PCR扩增反应循环数。通过鲮的分子内标Actin基因的相对表达量确定PCR扩增反应循环数，用于检测鲮IGF-ⅠmRNA组织表达。由图 2结果可知，循环数在25以下时，Actin扩增量较低，电泳条带不够清晰；循环数在30以上时结果较好，循环数为35最为理想(图 2)。
![]()
图 2 PCR循环数的确定M. 100 bp DNA分子量标准；15~45. 循环数分别为15、20、25、30、35、40、45Figure 2. Validation of PCR cyclesM. 100 bp DNA molecular weight marker; 15~45. Numbers on each lane represents the number of PCR cycles preformed.2.3 鲮IGF-ⅠmRNA组织表达与重组鲮GH对鲮IGF-ⅠmRNA组织表达的影响
通过半定量RT-PCR方法对经重组鲮GH处理前后鲮IGF-ⅠmRNA在各组织中表达检测结果显示，鲮IGF-ⅠmRNA在肝组织中表达最高，其次是肾和脑，另外在肠、鳃脾、性腺和心脏组织中也有表达，而在肌肉、皮肤中没有检测到鲮IGF-ⅠmRNA的表达；研究了经腹腔注射重组鲮GH(120 μg·g-1鱼体重) 12 h后，鲮IGF-Ⅰ表达变化情况，结果表明，经重组鲮GH处理后，鲮肝组织IGF-ⅠmRNA水平显著升高，脑组织IGF-ⅠmRNA水平也稍有升高，而在肌肉中仍未检测到IGF-ⅠmRNA的表达(图 3)。
![]()
图 3 半定量RT-PCR检测重组鲮GH处理前后鲮actin和IGF-ⅠmRNA的组织表达B1. 半定量RT-PCR检测对照组中鲮各组织中actin mRNA的表达；B2. 半定量RT-PCR检测对照组中鲮各组织中IGF-ⅠmRNA的表达；C1. 半定量RT-PCR检测处理组中鲮各组织中actin mRNA的表达；C2. 半定量RT-PCR检测处理组中鲮各组织中IGF-ⅠmRNA的表达；D. 处理前后鲮组织IGF-ⅠmRNA的相对表达；M. 100 bp DNA分子量标准；NC. 负对照；1~10. 肝、脑、肌肉、肾、肠、鳃、脾、性腺、心脏、皮肤Figure 3. Semi-quantitative RT-PCR analysis of expression of Actin and IGF-Ⅰ mRNA in various tissues in mud carpB1. Actin mRNA was detected in various tissues by semi-quantitative RT-PCR in control group; B2. IGF-ⅠmRNA was detected in various tissues by semi-quantitative RT-PCR in control group; C1. Actin mRNA was detected in various tissues by semi-quantitative RT-PCR in rmcGH in treated group; C2. IGF-Ⅰ mRNA was detected in various tissues by semi-quantitative RT-PCR in rmcGH in treated group; D. relative quantity of IGF-ⅠmRNA in various tissues in mud carp; M.100 bp DNA molecular weight marker; NC. a negative control (no template); 1~10. brain, liver, muscle, kidney, intestine, gill, spleen, gonad, heart, skin2.4 血清中IGF-Ⅰ的含量
血清中IGF-Ⅰ的含量，对照组为145.59±21.84 ng·mL-1，重组鲮GH处理组为247.71±2.83 ng·mL-1，经t检验，P=0.043 < 0.05，存在显著性差异；血清中GH的含量，对照组为0.037±0.0067 ng·mL-1，重组鲮IGF-Ⅰ处理组为0.04±0.0058 ng·mL-1，经t检验，P=0.508>0.05，故不存在显著性差异(表 2)。
表 2 对照组和重组鲮IGF-Ⅰ处理组鲮血清中GH和IGF-Ⅰ浓度Table 2. The GH and IGF-Ⅰconcentration in sera in control and rmc IGF-Ⅰtreated fish浓度±S.E.(n=3)/ng·mL-1 concentration±S.E. 对照组
control重组鲮IGF-Ⅰ处理组
rmcIGF-Ⅰtreated groupP值 胰岛素样生长因子-Ⅰ浓度
IGF-Ⅰconcentration145.59±21.84 247.71±2.83 P=0.043＜0.05 生长激素浓度
GH concentration0.037±0.0067 0.053±0.0067 P=0.508＞0.05 3. 讨论
采用半定量RT-PCR方法，检测到鲮肝、肾、脑、肠、鳃、脾、性腺、心脏组织都有IGF-ⅠmRNA表达，且在鲮肝组织中表达量最高，在鲮肾、脑、肠、鳃、脾、性腺、心脏组织中表达量较低；而在鲮皮肤和肌肉组织中未检测到有IGF-ⅠmRNA的表达。鲮IGF-ⅠmRNA组织表达的特异性，表明了鲮IGF-Ⅰ主要由肝组织合成并进入血液(内分泌功能)，其它组织也能合成少量的IGF-Ⅰ作用于自身细胞(自分泌或旁分泌功能)。这一结果与草鱼(Cirrhinus molitorella)、金鲷(Sparus aurata)以及其它一些脊椎动物中的结果一致[10-15]。
与哺乳类相似，生长激素是鱼类IGF-Ⅰ表达调控的主要影响因子之一。生长激素能诱导肝脏中IGF-Ⅰ和IGF-Ⅱ的合成[10, 16-17]。用放射免疫分析法分析经重组大麻哈鱼GH投喂或注射后，虹鳟的细胞内GH在较短时间内(12 h)升到最高值，并很快降低至正常水平，而IGF-Ⅰ逐渐升高并维持在较高水平至少96 h[18]。生长激素能提高鱼类IGF-ⅠmRNA表达水平，但生长激素对IGF-Ⅱ基因的转录调控可能在不同的物种具有不同的作用模式，例如，生长激素能使金鲷肝脏IGF-ⅠmRNA表达水平提高2.2倍，但IGF-ⅡmRNA表达水平不受生长激素的影响；而虹鳟IGF-Ⅰ和IGF-Ⅱ表达水平都受生长激素的影响。CAO[19]研究生长激素GH对鱼IGF-Ⅰ表达的调节作用，发现切除垂体后，肝组织IGF-ⅠmRNA丰度下降，而注射GH后IGF-ⅠmRNA水平又得到恢复。后来对GH家族对IGF-Ⅰ表达的调节作用进行的在体和离体试验均证明，GH能显著促进肝组织IGF-ⅠmRNA的表达，且具有剂量依赖性。该研究也证实了重组mcGH对鲮IGF-Ⅰ的表达和分泌都有促进作用。
在银大麻哈鱼(Oncorhynchus kisutch)和罗非鱼(O.mossambicus)中重组GH对肝组织IGF-ⅠmRNA表达有促进作用，但对肝外组织的IGF-ⅠmRNA表达水平没有影响[20-21]。该试验中，重组mcGH不仅促进了肝组织IGF-ⅠmRNA表达，而且促进肝外组织(脑组织)IGF-ⅠmRNA表达。BIGA等[10]利用重组牛生长激素(recombinant bovine somatotropin，rbST)处理虹鳟，发现不仅肝脏组织中IGF-ⅠmRNA水平升高，在鳃、性腺、脑、肠等一些肝外组织IGF-ⅠmRNA水平也有升高。TSE等[22]研究了用不同剂量的猪生长激素处理鲤(Cyprinus carpio)之后脑和肝脏中IGFs的表达情况，结果表明，生长激素能提高鲤脑和肝脏组织IGF-Ⅰ和IGF-ⅡmRNA的表达水平，且具有剂量依赖效应，脑中IGF-ⅡmRNA水平对生长激素反应更灵敏。GH能诱导肝细胞和其它细胞合成IGF-Ⅰ，并且能调节机体内IGF-Ⅰ浓度。PEREZ-SANCHEZ等[23]对虹鳟研究也表明重组GH刺激了肝和肝外组织IGF-ⅠmRNA表达。该试验中经重组mcGH处理，鲮脑组织中大于肝组织中IGF-ⅠmRNA增加量，这可能是器官特异性引起的不同的调节。
该试验中发现, 经外源GH处理，鲮血清中GH浓度变化不显著(表 2)。这可能由于重组mcGH不影响鲮内源GH的分泌。这一观点与LEEDOM等[24]利用重组牛生长激素刺激罗非鱼发现其内源GH变化并不明显的结果一致。
尽管重组mcGH刺激对鲮内源GH影响不大，但鲮内源IGF-Ⅰ显著升高(表 2)。这与对哺乳动物和其它一些脊椎动物中GH和IGF的研究结果一致[5, 25-26]。这也说明了在鲮中，GH对IGF-Ⅰ的分泌起到很重要的作用。
-
![]()
图 1 1987~2001年黄海鳀鱼亲体-补充量数据散点图
图中连线连接相邻年份的补充数据
Figure 1. Scatter plot of the stock and recruitment data for the anchovy (E.japonicas) stock in the Yellow Sea
The line joins recruitment values with the year indicated.
![]()
图 2 黄海千里岩1990~2001年4~7月份平均日表层水温数据(折线图)与黄海鳀鱼补充量数据(柱状图)变化趋势图
Figure 2. Variations of mean SST from April to July (line) and recruitment (bar) observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea in the period of 1990~2001
![]()
图 3 黄海千里岩1990~2001年平均磷酸盐浓度数据(折线图)与黄海鳀鱼补充量数据(柱状图)变化趋势图
Figure 3. Variations of mean phosphates concentration (line) and recruitment (bar) observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea in the period of 1990~2001
![]()
图 4 黄海千里岩环境条件指数αt和黄海鳀鱼补充量指数R变化趋势图
Figure 4. Variations of environmental conditions index and the recruitment observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea
![]()
图 5 黄海鳀鱼补充量观测值和计算值比较图
Figure 5. Comparison of observed and calculated recruitment of anchovy (E.japonicas) stock in the Yellow Sea
表 1 黄海千里岩水域环境因子与黄海鳀鱼补充量数据的相关系数
Table 1 The correlation of environmental factors and recruitment data observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea
环境因子
environmental factorsXi 相关系数
correlation coefficient年平均磷酸盐浓度
mean phosphates concentrationX1 0.4091 平均表层水温
average sea surface temperature1月 X2 0.2656 2月 X3 0.3707 3月 X4 -0.0761 4月 X5 -0.6213 5月 X6 -0.6614 6月 X7 -0.4987 7月 X8 -0.5290 8月 X9 -0.0451 9月 X10 -0.4223 10月 X11 -0.3174 11月 X12 -0.1593 12月 X13 -0.1838 全年 X14 -0.4679 4~7月 X15 -0.8096 
下载: 导出CSV
表 2 黄海鳀鱼SRR模型参数的估计值，AIC、BIC计算值
Table 2 Estimated parameters and AIC, BIC in stock recruitment models (SRR)
模型
model参数 estimated parameters AIC BIC α β α0 α1 α10 α15 Ricker模型
Ricker model1.5931 0.3058 - - - - 37.2062 38.176 加入αt的Ricker模型
Ricker model with αt- 0.3058 7.974 -0.0144 0.216 -0.713 9.3159 11.7404
下载: 导出CSV
-
[1] 张春霖, 成庆泰, 郑傈珊, 等. 黄渤海鱼类调查报告[M]. 北京: 科学出版社, 1955: 51-52. [2] 朱德山, IVERSEN S A. 黄、东海鳀鱼及其他经济鱼类资源声学评估的调查研究[J]. 海洋水产研究, 1990, 11: 18-31. [3] RICKER W E. Computation and interpretation of biological statistics of fish populations[J]. Bull Fish Res Board Can, 1975, 191(18): 1-382. https://www.semanticscholar.org/paper/Computation-and-interpretation-of-biological-of-Ricker/5a8d0d7094c356e3b851fd66bd929ed0e56aabfd
[4] HILBORN R, WALTERS C J. Quantitative fisheries stock assessment: Choice, dynamics and uncertainty[M]. New York: Chapman & Hall, 1992. doi: 10.1007/978-1-4615-3598-0
[5] BEVERTON R J H, HOLT S J. On the dynamics of exploited fish populations[J]. Fish Inverst, 1957, 19(2): 1-533. doi: 10.1007/bf00044132
[6] RICKER W E. Stock and recruitment[J]. J Fish Res Bd Can, 1954, 108(11): 599-623. doi: 10.1139/F54-039
[7] QUINN Ⅱ T J, DERISO R B. Quantitiative fish dynamics[M]. New York: Oxford University Press, 1999. doi: 10.1093/oso/9780195076318.001.0001
[8] ZHAO X, HAMRE J, LI F, et al. Recruitment, sustainable yield and possible ecological consquences of the sharp decline of the anchovy (Engraulis japonicus) stock in the Yellow Sea in the 1990s[J]. Fish Oceanogr, 2003, 12(4): 495-501. doi: 10.1046/j.1365-2419.2003.00262.x
[9] WANG Yanjun, LIU Qun, YE Zhenjiang. A bayesian analysis on the anchovy stock (Engraulis japonicus) in the Yellow Sea[J]. Fish Res, 2006, 82(1/3): 87-94. doi: 10.1016/j.fishres.2006.08.007
[10] 詹秉义. 渔业资源评估[M]. 北京: 中国农业出版社, 1995: 208-215. https://book.douban.com/subject/2048040/ [11] 邓景耀, 叶昌臣. 渔业资源学[M]. 重庆: 重庆出版社, 2000: 124-131. https://www.dushu.com/book/10041543/
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
粤公网安备 44010502001741号
