Sequence analysis of the internal transcribed spacer 1(ITS-1) of ribosomal DNA gene of Lutjanus fulviflamma
-
摘要:
以特异性引物扩增了金焰笛鲷(Lutjanus fulviflamma)的核糖体第一转录间隔区(ITS-1),扩增产物经克隆后测序,测得ITS-1长度为566 bp。其中A、T、G、C 4种碱基的含量分别为14.1%、16.1%、30.2%、39.6%,G+C(69.8%)含量明显高于A+T含量(30.2%)。将此引物在笛鲷属其他4种鱼类中扩增,发现该对引物有很好的通用性;比较发现在不同种中ITS-1存在着较大的差异,适合将其应用于分子系统学和种质资源方面的研究。
-
关键词:
- 金焰笛鲷 /
- 内转录间隔区(ITS-1) /
- 序列分析
Abstract:The PCR technique was used to amplify rDNA-ITS-1 of Lutjanus fulviflamma, then the purified PCR productions were cloned into T-vector and sequenced by M13+/-primers.As a result, 566 bp nucleotide sequences of rDNA-ITS1 were obtained. The average contents of A, T, G and C were 14.1%, 16.1%, 30.2% and 39.6%, respectively, the contents of GC (69.8%) were obviously higher than those of AT. Using these primers, ITS-1 region can be amplified in other four Lutjanus species, but the lengths of ITS-1 were largely different between different species.
-
Keywords:
- Lutjanus fulviflamma /
- ITS-1 /
- sequence analysis
-
初级生产力的流向(Fate)是当前海洋生态学研究的热点之一。初级生产力以颗粒有机物和溶解有机物2种方式存在于水体。细菌能迅速利用溶解有机物这部分初级生产力,但由于其体积小而不能被大中型浮游动物(体长>200 μm)有效摄食。而微型浮游动物(体长2~200 μm)能大量摄食细菌,自身又是meso级浮游动物(体长约200~2 000 μm)的饵料,从而使得这部分溶解态的初级生产力得以传递到食物网的高层[1]。因此,微型浮游动物在微食物网与传统食物链连接的能量通道中有着重要意义。研究海洋生态系统的营养动力学机制,微型浮游动物的摄食及能量转化是其研究重点之一[2-3]。在一些海区,微型浮游动物能大量摄食浮游植物,在避免赤潮的产生中也扮演重要角色[4-6]。我国学者利用稀释法或荧光标记法研究了胶州湾,东海,香港海域、厦门海域,厦门杏林虾池等水域微型浮游动物对浮游植物的摄食压力[7-14]。研究结果表明,在不同的海域,微型浮游动物对浮游植物的摄食压力存在较大差异。大亚湾是我国目前水域生物多样性保存良好的重要海湾,也是我国重要的亚热带物种种质资源库。“低营养,高生产力”是该湾不同于其它海湾的一个重要特点[15]。近20年来,大亚湾由贫营养状态发展到中营养状态,大亚湾营养限制因子由20世纪80年代的N限制过渡到目前的P限制[16]。在生态环境变化过程中,微型浮游动物在海湾生态系统中的重要性是否发生变化仍未知。有关大亚湾微型浮游动物的研究较少[17],其摄食生态学方面的研究还未开展。
本文利用稀释法,以叶绿素a的浓度来表示浮游植物的浓度,对大亚湾海域5个站位微型浮游动物对浮游植物的摄食压力进行了初步研究,为了解大亚湾生态系统的结构和能量流动提供基本资料和参数。
1. 材料与方法
1.1 调查站位
2005年3月13日至17日,在大亚湾湾口以内约500 km2的海域设置5个站位(图 1),采用稀释法研究微型浮游动物对浮游植物的摄食压力。
1.2 调查与分析方法
微型浮游动物摄食率应用LANDRY等[18]提出的稀释法测量。在各站位采集距表层0.5 m水样,按《海洋监测规范》操作,其中,pH、DO分别用pH计、YSL Model 55溶氧测定仪现场测定,盐度用比重法测定,化学需氧量(COD)用碱性高锰酸钾法测定,DIN用重氮-偶氮法、次溴酸钠氧化法、锌镉还原法测定,无机磷用磷钼兰分光光度法测定。
微型浮游动物摄食率实验用水样经200 μm筛绢过滤后,一部分水样用0.2 μm滤膜(孔径47 mm)Gelman滤器过滤以获得无颗粒水(particle-free water,PFW),另一部分水样与PFW按2∶0,2∶0.5,2∶1,2∶1.5比例混合,轻轻倒入2 L的磨口玻璃培养瓶中(培养瓶使用前经10%盐酸浸洗10 h,用自来水冲洗干净,每个稀释比例设2个平行样),放入培养箱,在甲板上利用自然海水流动循环培养24 h,使实验条件尽可能接近自然条件。培养期间多次晃动培养瓶,以使微型浮游动物和浮游植物混合均匀。
培养前,用孔径0.2 μm的滤膜过滤海水1 L,滤膜迅速冷冻保存。在实验室内,滤膜放入具塞刻度试管中用5 mL 90%丙酮,于20℃暗处萃取24 h。上清液用722光栅分光光度计测定吸光值,依据JEFFREY等[19]的改进公式计算叶绿素a的含量。培养24 h后用上述方法采集、测定培养后叶绿素a值。
1.3 微型浮游动物的摄食压力
稀释法最早由LANDRY等[18]提出,现已成为研究微型浮游动物摄食的常用方法之一。假设海水中浮游植物的生长率为μ,微型浮游动物的摄食率为g,浮游植物处于指数增长期,培养前的浓度为P0,培养后浓度为Pt ,那么Pt= P0e(μ-g) t,其中t为培养时间。方程可表示为:
$$ \frac{\ln \left(P_t / P_0\right)}{t}=\mu-g $$ (1) 计算每个培养瓶的浮游植物表观生长率(apparent growth rate,AGR) $\frac{\ln \left(P_t / P_0\right)}{t}$,然后计算实际稀释因子(actual dilution factor,ADF):
$$ \mathrm{ADF}=P_0\left(X_i\right) / P_0\left(X_0\right) $$ (2) 其中,P0(Xi)是初始培养处理中Xi组分的浮游植物现存量,Xi=1-无颗粒水占培养水的比例,P0(X0)是初始培养处理中未稀释的浮游植物现存量。实际稀释因子也称稀释度或稀释因子,也可表示为自然海水与混合海水的体积比。微型浮游动物的摄食率(g)和浮游植物的内禀生长率(μ)可用AGR和ADF的线性回归方程获得,其中截距为浮游植物的内禀生长率(μ),斜率为微型浮游动物的摄食率(g)
微型浮游动物的摄食影响可以用浮游植物净生长率(net growth rate of phytoplankton, NGR)、对浮游植物现存量的摄食压力(Ps)、对浮游植物潜在初级生产力的摄食压力(Pp)表示,计算公式如下:
$$ \text { NGR }=\mu-g $$ (3) $$ P_s=\frac{\left(C_0 e^\mu-C_0\right)-\left(C_0 e^{\mu-g}-C_0\right)}{C_0} \times 100 $$ (4) $$ P_s=\frac{\left(C_0 e^\mu-C_0\right)-\left(C_0 e^{\mu-g}-C_0\right)}{C_0 e \mu-C_0} \times 100 $$ (5) 另外,浮游植物加倍的时间(Td),每天的加倍数(n),公式如下:
$$ T_d=\ln 2 / \mu $$ (6) $$ n=\mu / \ln 2 $$ (7) 2. 结果
2.1 研究海区的水环境
各站位的水温变化较小,盐度总体上呈近岸向湾口递减的趋势。叶绿素a的高值区出现在湾北部的S1(6.55 μg·L-1),低值区则出现在湾东部的S7(0.70 μg·L-1),叶绿素与盐度的时空变化特点较相似。S1、S4站的COD值最高,这可能与沿岸施工及港口排污等等有关。各测站的环境参数见表 1。
表 1 各测站的表层环境参数Table 1. Environmental parameters at stations站位
stations水深/m
depth盐度*
salinity水温/℃
temperature叶绿素a/μg·L-1
Chl-a化学需氧量/mg·L-1
CODS1 5.8 32.21 15.4 6.55 6.73 S4 10.5 32.09 15.6 2.45 7.05 S7 10.0 31.82 15.2 0.70 6.20 S8 10.1 32.09 15.2 1.93 6.66 S11 18.0 31.69 15.3 1.70 6.80 注:* 统一校正到17.5℃时的盐度
Note:* Salinity is standardized to 17.5℃.2.2 浮游植物的生长率和微型浮游动物的摄食压力
稀释实验的结果见表 2,实验期间,S1和S11站浮游植物的生长率大于微型浮游动物的摄食率,而S4、S7和S8站浮游植物的生长率则小于微型浮游动物的摄食率。
表 2 大亚湾微型浮游动物的摄食结果分析Table 2. Analysis on the results of grazing of the microzooplankton in Daya Bay站位
stations采样水深/m
sampling depthR2 μ/d-1 G/d-1 NGR/d-1 %Ps %Pp Td/d n/d-1 S1 0.5 0.61 0.90 0.88 0.02 144 99 0.77 1.30 S4 0.5 0.90 0.59 0.96 -0.37 111 138 1.17 0.85 S7 0.5 0.76 1.14 1.74 -0.6 258 121 0.61 1.64 S8 0.5 0.53 0.16 0.67 -0.51 57 330 4.33 0.23 S11 0.5 0.45 0.51 0.33 0.18 47 70 1.36 0.74 注:R2. 由表观生长率[ln(Pt/P0)]与实际稀释因子(ADF)得到的回归直线中的相关因子;μ. 浮游植物内禀生长率;g. 微型浮游动物摄食率;NGR. 浮游植物净生长率;%Ps. 对浮游植物现存量的摄食压力;%Pp. 对浮游植物潜在初级生产力的摄食压力;Td. 浮游植物数量加倍所需的时间;n. 每天浮游植物加倍的数量
Note:R2. correlation coefficient;μ. instantaneous growth rate of phytoplankton;g. ingestion rate of microzooplankton;NGR. net growth rate of phytoplankton;%Ps. percentage of phytoplankton standing crop ingested by microzooplankton;%Pp. percentage of phytoplankton potential production ingested by microzooplankton;Td. time of doubling;n. times per dayS1站浮游植物现存量最高,其生长率也较高,且大于微型浮游动物的摄食率,浮游植物的实际生长率为正值(0.02),浮游植物的生物量很可能会增加。
S4站浮游植物的现存量、生长率、微型浮游动物的摄食率等各项指标值均处于5个站对应指标值的中间水平。
S7站浮游植物的生长率最高,但其现存量最低,微型浮游动物的摄食造成此现象的主要原因之一。在5个站位中,S7站微型浮游动物的摄食率是最高的,微型浮游动物对浮游植物现存量的摄食压力高达258%。S7站微型浮游动物的摄食对浮游植物潜在生产力的压力并不是最高,这预示着微型浮游动物对浮游植物的摄食作用有减小的趋势。浮游植物的生长是多种生物活动综合作用的结果,并不仅仅取决于微型浮游动物的摄食,该站浮游植物现存量仅需0.61 d(14 h左右)即可翻倍,但由于微型浮游动物的摄食,其现存生物量仍很低。可见,在S7站,微型浮游动物的摄食起着重要作用。
S8站浮游植物初始现存量较高,但浮游植物的生长率和微型浮游动物的摄食率都较低,微型浮游动物对浮游植物现存量的摄食压力最小(57%),但潜在初级生产力的摄食压力却最大。与S7站相反,该站浮游植物现存量翻倍所需的时间最长,为4.33 d。
S11站回归方程的相关系数R2较低,这可能与实验过程中一个培养瓶取样失败使实验数据缺失一个有关,该站浮游植物的现存量及生长率均处于中等水平,但微型浮游动物的摄食率最小,浮游植物的净生长率(NGR)为0.18,对浮游植物现存量和潜在生产力的摄食压力都是最小的,分别为47%,70%,但此站浮游植物现存量翻倍所需的时间却较长,这可能与该站的营养盐或其它因素有关,具体原因有待于今后进一步研究。
3. 讨论
浮游生态系中浮游植物与浮游动物的摄食关系研究是海洋生态系统中一个关键的过程研究。浮游动物的摄食受各营养阶层组成粒级结构的影响。在大细胞浮游植物占优势的海区,浮游植物被中型浮游动物摄食的经典食物链较重要,而在以小细胞植物为主的生态系统中,则以微型浮游动物对小细胞的浮游植物及细菌摄食的微食物环为主。近年来发现微型浮游生物在海洋生态系统中扮演重要角色,浮游动物摄食研究的重点就转向了微型浮游动物的摄食研究[18, 20-25]。评估微型浮游动物的摄食对浮游植物现存量的影响,对了解有害微藻水华发生的动力学机制有重要意义。
总体来讲,微型浮游动物的摄食速率、对浮游植物现存量的摄食压力由近岸区向湾中部呈增加趋势,在S7、S8站达到峰值,然后再向湾外呈减少趋势。若不考虑S7(由于浮游植物的初始浓度较低,微型浮游动物的摄食压力估计值可能偏高),其余4个站的研究结果与最优摄食理论(optimal foraging theory)相符合,摄食者对被摄食者的清除率(clearance rate,相当于本研究中的摄食率)随着被摄食者的浓度增加会达到一个最大值,但随后就降低(图 2)。在S1站,虽然浮游植物现存量最大,但摄食率却低于S4站。这反映微型浮游动物在控制浮游植物生长方面有重要作用,但若浮游植物的密度过高时,浮游动物是无法控制浮游植物生长的,具体阈值有待进一步研究。
由于缺乏本海区的有关微型浮游动物研究的历史资料,难以进行纵向比较。与其它海区比较(表 3),大亚湾微型浮游动物的摄食压力处于较高水平;与国内海区相比,比胶州湾、香港东部海域、厦门海域高,与渤海、东海海域大致相当。
表 3 不同海区微型浮游动物的对浮游植物潜在初级生产力摄食压力Table 3. Comparasion of microzooplankton grazing pressure on potential primary production海区
sea areas对浮游植物潜在初级生产力的摄食压力
grazing pressure on primary production (%Pp) /d-1参考文献
citation旧金山湾San Francisco Bay 44~722 [20] 墨西哥湾Gulf of Mexico 30 [21] 纽芬兰鲁及湾Logy Bay, Newfoundland 64~118 [22] 法国比斯开湾BiscayBay, France 73~136 [23] 加拿大北部琼斯海峡Jones Sound,NWT 40~114 [24] 加拿大北部巴芬湾Baffin Bay,NWT 37~88 [24] 华盛顿沿岸Washington coast 17~52 [22] 新西兰沿岸NZ coast 20~194 [25] 南极普里兹湾Prydz Bay,Antarctic 4~60 [11] 渤海Bohai sea, China 95~319 [10] 胶州湾内Jiaozhou Bay, China 53~93 [7] 胶州湾外Outside Jiaozhou Bay, China 74~84 [7] 东海East China Sea 74~203 [26] 香港东部海域East of Hong Kong water 79 [8] 香港西部海域West of Hong Kong water 127 [8] 厦门海域Xiamen water, China 52~102 [13] 厦门杏林虾池Xiamen Xinglin Shrimp Ponds, China 37~194 [12] 在大亚湾,微型浮游动物对初级生产有着重要的调控作用。近年来的研究表明,病毒在微食物环中也起着重要作用,病毒的存在使得微食物环的能流和物质循环更趋向复杂。今后,应进一步加强微食物环及大型浮游动物对微型浮游动物摄食的研究,因为微型浮游动物使碳通量减小,而大型浮游动物是碳通量的主要贡献者[27]。
-
-
[1] 刘丽, 刘楚吾, 许强. 运用RAPD技术对笛鲷属3种鱼的群体遗传学分析[J]. 湛江海洋大学学报, 2003, 23(3): 16-21. doi: 10.3969/j.issn.1673-9159.2003.03.004 [2] 王中铎, 刘楚吾, 郭昱嵩. 5种笛鲷mtDNA及Cytb基因片段的RFLP比较[J]. 水产学报, 2005, 29(3): 327-332. doi: 10.3321/j.issn:1000-0615.2005.03.008 [3] HILLIS M D M, DAVIS S K. Ribosomal DNA: Intraspeciefic polymorphism concerted evolution and phylogeny reconstruction[J]. Syst Zool, 1988, 37(1): 63-66. doi: 10.2307/2413191
[4] 郑雪松, 杨虹, 李道棠, 等. 基因间隔序列(ITS)在细菌分类鉴定和种群分析中的应用[J]. 应用与环境生态学报, 2003, 9(6): 678-684. http://www.cibj.com/#/digest?ArticleID=972 [5] 唐伯平, 周开亚, 宋大祥. 核rDNA ITS区在无脊椎动物分子系统学研究中的应用[J]. 动物学杂志, 2002, 37(4): 67-73. doi: 10.3969/j.issn.0250-3263.2002.04.017 [6] 中国科学院动物研究所, 中国科学院海洋研究所, 上海水产学院. 南海鱼类志[M]. 北京: 科学出版社, 1962: 469-470. [7] 萨姆布鲁克J, 弗里克E F, 曼尼阿蒂斯T. 分子克隆实验指南[M]. 北京: 科学出版社, 1995. [8] PHILLIPS R B, SAJDAK S L, DOMANICO M J. Relationships among charrs based on DNA sequences[J]. Nordic J Freshw Res; 1995, 71(2): 378-391. doi: 10.5433/1679-0359.2012v33Supl2p3287
[9] 喻达辉, 朱嘉濠. 珠母贝属6个种的ITS2分子标记研究[J]. 南方水产, 2005, 1(4): 6-12. doi: 10.3969/j.issn.2095-0780.2005.04.002 [10] 喻达辉, 朱嘉濠. 珠母贝属的系统发育: 核rDNA ITS序列证据[J]. 生物多样性, 2005b, 13(4): 315-323. doi: 10.1360/BIODIV.050050 [11] 喻达辉, 李有宁, 吴开畅. 中国、日本和澳大利亚珍珠贝的ITS2序列特征分析[J]. 南方水产, 2005, 1(2): 1-6. doi: 10.3969/j.issn.2095-0780.2005.02.001 [12] 喻达辉, 朱嘉濠, 贾晓平. 我国珠母贝属(Pinctada)主要种类亲缘关系的初步分析[J]. 海洋与湖沼, 2006, 37(3): 211-217. doi: 10.3321/j.issn:0029-814X.2006.03.004 [13] 何毛贤, 黄良民. 长耳珠母贝核rRNA基因ITS 2序列分析[J]. 热带海洋学报, 2004, 23(5): 81-84. doi: 10.3969/j.issn.1009-5470.2004.05.011 [14] 陈琳林, 孔晓瑜, 周立石, 等. 魁蚶核糖体DNA基因转录间隔区的序列特征[J]. 中国水产科学, 2005, 12(1): 104-107. https://www.fishscichina.com/zgsckx/article/abstract/3570?st=article_issue [15] 魏晓华. 栉孔扇贝和海湾扇贝的遗传多样性研究及扇贝科几种贝类的分子系统学研究[D]. 青岛: 中国海洋大学, 2004. [16] KONG Xiaoyu, ZHANG Liusuo, YU Ziniu, et al. Sequencing of ribosomal internal transcribed spacer regions and mitochondrial genef ragments in Crassostrea gigas[J]. J Fish Sci China, 2002, 9(4): 304-308. https://www.semanticscholar.org/paper/Sequencing-of-ribosomal-internal-transcribed-spacer-Kong-Zhang/97cabf934e3b697b3fe23d1139517cf7394be072
[17] YU Ziniu, KONG Xiaoyu, ZHUANG Zhimeng, et al. Sequence study and potential uses of ribosomal DNA internal transcribed spacers in scallop Chlamys farreri[J]. J Fish Sci China, 2001, 8(1): 6-9. https://europepmc.org/article/CBA/575294
[18] DOMANICO M J, PHILLIPS R B, OAKLEY T H. Phylogenetic analysis of Pacific salmon (genus Oncorhynchus) using nuclear and mitochondrial DNA sequences[J]. Can J Fish Aquat Sci, 1997, 54(8): 1865-1872. doi: 10.1139/f97-092
[19] LINN S S, PHILLIPS R B. Phylogenetic relationships among Coregonus species inferred from the DNA swquence of the first internal transcribed spacer (ITS1) of ribosomal DNA[J]. Can J Fish Aquat Sci, 1997, 54(8): 1494-1503. doi: 10.1139/f97-057