Limiting concentrations of nitrogen, phosphorus, iron on growth of Oscillatoria chlorine
-
摘要:
绿色颤藻(Oscillatoria chlorine)在虾池中大幅度增加时会释放微囊藻毒素, 给对虾的健康养殖带来危害。文章分析不同质量浓度的氮(N)、磷(P)、铁(Fe3+)对绿色颤藻叶绿素a和藻体质量浓度的影响, 探讨水中N、P、Fe3+的质量浓度对绿色颤藻生长的限制性条件。结果显示, 限制绿色颤藻叶绿素a质量浓度增加的条件是ρ(N)<0.06 mg·L-1, ρ(P)<9.93×10-2 mg·L-1和>1.99 mg·L-1, ρ(Fe3+)<2.58×10-2 mg·L-1和>51.67 mg·L-1; 限制藻体质量浓度增加的条件是ρ(N)<12.36 mg·L-1, ρ(P)<9.93×10-3 mg·L-1, ρ(Fe3+)<5.16×10-3 mg·L-1和>51.67 mg·L-1。绿色颤藻生长的阈值ρ(N)为0.06 mg·L-1、ρ(P)为9.93×10-3 mg·L-1、ρ(Fe3+)为5.16×10-3 mg·L-1。ρ(N)为7.42 mg·L-1、ρ(P)为4.96×10-2 mg·L-1和ρ(Fe3+) 为5.16×10-3 mg·L-1是限制绿色颤藻叶绿素和藻体质量浓度增加的最低条件组合; N是影响绿色颤藻叶绿素a和藻体质量浓度的主要因子, 其次是Fe3+和P。在对虾养殖过程中可以通过换水或循环水系统除去水体中的含N污染物, 降低养殖水体中溶解态N含量, 从而控制绿色颤藻的快速增长。
Abstract:The proliferation of Oscillatoria chlorine will release toxins which are harmful to prawn growth in the aquaculture ponds. We investigated the limiting concentrations of nitrogen (N), phosphorus (P) and iron (Fe3+) to the growth of O.chlorine as well as the effects of different N, P and Fe3+ concentrations on chlorophyll a content and dry weight of O.chlorine. The results show that N concentration < 0.06 mg·L-1, P concentration < 9.93×10-2 mg·L-1 and >1.99 mg·L-1, and Fe3+ concentration < 2.58×10-2 mg·L-1 and > 51.67 mg·L-1, limited the increase of chlorophyll a content; N concentration < 12.36 mg·L-1, P concentration < 9.93×10-3 mg·L-1 and the Fe3+ concentration < 5.16×10-3 mg·L-1 and > 51.67mg·L-1, limited the increase of dry weight. The threshold concentrations of N, P and Fe3+for the growth of O.chlorine were 0.06 mg·L-1, 9.93×10-3 mg·L-1 and 5.16×10-3 mg·L-1, respectively. The combined concentrations which limited the increase of chlorophyll a contents and dry weight of O.chlorine were 7.42 mg·L-1 of N, 4.96×10-2 mg·L-1 of P and 5.16×10-3 mg·L-1 of Fe3+, respectively. Nitrogen had the strongest effect on chlorophyll a content and dry weight, followed by Fe3+ and P. It is suggested that the rapid proliferation of the O.chlorine might be controlled through water exchange or circulation to reduce the dissolved nitrogen content in the aquaculture ponds.
-
Keywords:
- Oscillatoria chlorine /
- nitrogen /
- phosphorus /
- iron /
- chlorophyll a content /
- dry weight
-
近海海域在中国经济发展中具有重要的服务功能和战略地位。随着对近海海域的开发和利用,其环境逐渐面临诸多污染问题,其中石油类污染是主要污染问题之一。据国际油轮船东防污染联合会 (International Tanker Owners Pollution Federation, ITOPF) 的数据显示,2023年约有2 000 吨石油意外泄漏到环境中。国际海事组织 (International Maritime Organization, IMO) 和国际能源署 (International energy agency, IEA) 的报告显示,全球海上运输的油品中,原油和柴油的运输占比较高[1],因此原油和柴油污染对海洋生物的影响值得关注。此外,原油和柴油可以引起鱼类代谢紊乱、心脏速率减缓,甚至幼鱼发育畸形等问题[2];也会对多毛类[3]、虾类[4]、海洋藻类[5]、贝类[6]等生物的生长发育造成影响。因此,研究溢油污染对海洋生物生长发育的毒性效应,对保护海洋生态环境和渔业经济具有重要意义。
生物标志物是指环境因子对生物体产生的影响所导致的可测定变化,具有高敏感性的特点[7]。抗氧化酶系统和基因表达水平常用作生物标志物,指示污染及评估其生物毒性效应。蒋凤华等[8]使用超氧化物歧化酶 (SOD) 和过氧化氢酶 (CAT) 2种抗氧化酶生物标志物评估了不同浓度原油水溶性组分 (Water soluble fraction of crude oil, WSF) 对栉孔扇贝 (Chlamys farreri) 鳃和消化腺的毒性效应,结果显示SOD和CAT的活性随WSF暴露时间的增加呈降低-升高-降低的趋势;Xia等[9]通过测定菲律宾蛤仔 (Ruditapes philippinarum) 的抗氧化酶活性,发现石油类污染与该生物标志物之间存在剂量-效应和时间-效应关系。
综合生物标志物响应 (Integrated biological response, IBR) 指数通过将多个生物标志物与图形结合,从而更全面地判断生物的毒性响应,评估生态系统的健康状态[10]。Wang等[11]通过IBR指数评价重金属汞对鲢 (Hypophthalmichthys molitrix) 幼体的慢性毒性效应,发现汞会降低其抗氧化能力,导致机体功能受损;Jiang等[12]使用IBR指数评估缢蛏 (Sinonovacula constricta) 4种抗氧化酶对柴油和原油含水馏分 (Water accommodated fraction, WAF) 的氧化应激反应,结果显示不同浓度的石油污染物对4种抗氧化酶的诱导程度不同,且柴油的毒性大于原油。
凡纳滨对虾 (Litopenaeus vannamei) 又名南美白对虾,其生长快、抗逆性强、出肉率高且经济价值高,是全球养殖产量最高的经济虾类,2023年中国的养殖产量达220万吨[13]。凡纳滨对虾作为具有开放式血淋巴管系统的变温动物,其生理状态更易受到水体环境变化的影响[14]。本研究探讨了LH原油和0# 柴油对凡纳滨对虾肝胰腺抗氧化酶活性及相关功能基因表达的影响,以期为评估石油类污染对海洋生物的毒性效应提供基础数据支持。
1. 材料与方法
1.1 实验生物
凡纳滨对虾来自中国水产科学研究院南海水产研究所深圳试验基地,平均体质量为(6.50±0.21) g。正式实验前,将对虾置于室内300 L养殖桶中暂养7 d,养殖期间保持水温 (29.2±0.5) ℃、盐度30‰、pH 7.8~8.0,24 h连续充气,每天早晚各投喂1次,每日对虾配合饲料投喂量为虾体质量的5%,对虾日死亡率低于1‰ 时开始正式实验。
1.2 实验试剂及母液制备
1.2.1 试剂
南海LH原油来源于中国南海流花油田;0# 柴油采购自广州市加油站;富肯-3号溢油分散剂由广州富肯环保科技有限公司提供;正己烷、乙醇、乙醚等试剂购自国药控股广州有限公司,均为分析纯。
1.2.2 南海LH原油和0# 柴油母液及工作液制备
南海LH原油和0# 柴油乳化母液的配制方法参考文献 [15]:将分散剂分别按体积比1∶1与南海LH原油、0# 柴油混合后,置于超声波洗涤器中超声处理5 min,然后将混合液按体积比1∶10与过滤海水混合,连续超声处理3 h,再置于分液漏斗中静置3 h,取其下层水相液置于棕色瓶中,即为LH原油和0# 柴油母液。制备后的LH原油和0# 柴油母液参照GB 17378.4—2007《海洋监测规范 第4部分:海水分析》,采用紫外分光光度计法测定石油烃含量。在每次配制工作液前,先测定1次母液的石油烃浓度,随后使用分散剂梯度稀释母液用于制备实验用工作液,以保证每组工作液中分散剂浓度一致。工作液均为现用现配。
1.3 胁迫实验及样品采集
根据前期本实验室开展的急性毒性实验结果[15],LH原油和0# 柴油均设置3个实验浓度组,分别为低浓度组 (A组,25 μg·L−1)、中浓度组 (B组,50 μg·L−1) 和高浓度组 (C组,250 μg·L−1),和1个分散剂对照组 (D组,0 μg·L−1)。将240尾对虾分为4组,每组3个平行,每个平行20尾。连续胁迫15 d,每8 h换水1次,换水后使用工作液调节石油烃浓度至所设置浓度,分别于胁迫第1、第7和第15天在投喂及换水前取样,每组6尾,取肝胰腺迅速在液氮中冷冻,并储存于 −80 ℃冰箱。
1.4 样品处理及生理指标测定
将对虾肝胰腺取出,液氮研磨后分为2份,分别用于生理指标测定和基因表达分析。按照试剂盒说明书 (南京建成生物工程研究所) 测定SOD、CAT、过氧化物酶 (POD)、谷胱甘肽S转移酶 (GST)、谷胱甘肽过氧化物酶 (GPx) 的活性及还原型谷胱甘肽 (GSH) 的含量。
1.5 基因表达分析
根据已发表的文献,以凡纳滨对虾Lvactin为内参基因,选择细胞凋亡因子caspase2—caspase5、延伸因子 (EF)、SOD、CAT、GPx、过氧化物还原酶 (Prx)、酚氧化酶原 (proPO) 和溶菌酶 (LZM) 基因进行表达量的测定。引物序列见表1。
表 1 荧光定量PCR引物Table 1. Primers for fluorescent quantitative PCR引物名称
Primer引物序列 (5'—3')
Primer sequence (5'−3')退火温度
Annealing temperature/
℃参考文献
Reference引物名称
Primer引物序列 (5'—3')
Primer sequence (5'−3')退火温度
Annealing temperature/
℃参考文献
ReferenceLvcasp2-F ATGGCTCGTGGTTCATTCAG 57.8 [16] LvCAT-F ATCCATTCGACCTTACCA 52.7 [17] Lvcasp2-R CATCAGGGTTGAGACAATACAGG 60.2 [16] LvGPx-F GGCACCAGGAGAACACTAC 59.7 [17] Lvcasp3-F AGTTAGTACAAACAGATTGGAGCG 58.6 [16] LvCAT-R ACGCAATCTGCTCCACCT 57.3 [17] Lvcasp3-R TTGTGGACAGACAGTATGAGGC 60.1 [16] LvGPx-R CGACTTTGCCGAACATAAC 55.4 [17] Lvcasp4-F CATGCTTGACATACCCGATG 57.8 [16] LvPrx-F AGGAAGGTGGTCTTGGTACA 57.8 [17] Lvcasp4-R TGTCCGGCATTGTTGAGTAG 60.2 [16] LvPrx-R GCAGATCATTGATTGTTACCT 54.1 [17] Lvcasp5-F GAAGGAGTGAAGCTAAACGAGAC 60.3 [16] LvPro-F CGGTGACAAAGTTCCTCTTC 57.8 [18] Lvcasp5-R CAGTAGACCAGCAGATAAGGAAGT 61.9 [16] LvPro-R GCAGGTCGCCGTAGTAAG 59.6 [18] LvEF-1α-F GAAGTAGCCGCCCTGGTTG 61.9 [16] LvLZM-F GGACTACGGCATCTTCCAGA 59.8 [18] LvEF-1α-R CGGTTAGCCTTGGGGTTGAG 61.9 [16] LvLZM-R ATCGGACATCAGATCGGAAC 57.8 [18] LvSOD-F GAATTGGGTGAGGAACGA 55.0 [17] Lvactin-F CAAGATGTGTGACGACGAAG 57.8 [17] LvSOD-R ATGTAGCCCTGGTGATGC 57.3 [17] Lvactin-R AGCTCGTGGTGCATCGTC 59.6 [17] 1.6 荧光定量PCR
采用TaKaRa SYBR® Green I荧光定量试剂盒在罗氏LC480荧光定量PCR仪上进行荧光定量PCR分析,体系配制及PCR反应条件参照试剂盒说明书。
1.7 数据分析
荧光定量结果分析使用2−ΔΔCt法计算基因的相对表达量[19];IBR指数采用第二代综合生物标志物响应法 (Integrated biological response version 2, IBRv2),计算方法参照文献 [20-21];实验数据用Excel 2019软件进行统计整理,以“平均值±标准差 ($\overline { x}\pm s $)”表示;利用IBM SPSS 19.0软件进行单因素方差分析 (One-way ANOVA),差异显著性通过Tukey's HSD法检验,显著性水平α为 0.05;采用GraphPad Prism 9.0和Origin 2021软件制图。
2. 结果
2.1 生理指标分析
随LH原油胁迫时间的延长,中、高浓度组SOD、CAT、POD、GST和GPx活性总体呈升高趋势,并显著高于对照组 (p<0.05),表现出一定的时间-效应和剂量-效应 (图1)。低浓度组SOD、CAT、GST和GPx活性在第1和第7天均与对照组无显著性差异 (p>0.05)。GSH含量总体呈先上升后下降的趋势,各浓度组均在第7天达到最大值,在第15天时各浓度组GSH含量均显著低于对照组 (p<0.05)。
图 1 不同浓度LH原油乳化液对凡纳滨对虾肝胰腺抗氧化酶的影响注:* 表示组间存在显著性差异 (p<0.05);** 表示组间存在非常显著性差异 (p<0.01);*** 表示组间存在极显著性差异 (p<0.001)。Figure 1. Effects of LH crude oil emulsion on hepatopancreatic antioxidant enzyme activities of L. vannameiNote: * represents significant differences among the groups (p<0.05); ** represents very significant differences among the groups (p<0.01); *** represents extremely significant differences among the groups (p<0.001).随0# 柴油胁迫时间的延长,SOD、CAT、GST活性和GSH含量均显著低于对照组 (p<0.05),表现出一定的时间-效应关系 (图2)。中、高浓度组POD活性总体呈先升高后降低的趋势,至胁迫后期 (第15天) 各浓度组均显著低于对照组 (p<0.05)。GPx活性总体呈上升趋势,各浓度组均在第15天达最大值 (p<0.05)。
图 2 不同浓度0# 柴油乳化液对凡纳滨对虾肝胰腺抗氧化酶的影响注:* 表示组间存在显著性差异 (p<0.05);** 表示组间存在非常显著性差异 (p<0.01);*** 表示组间存在极显著性差异 (p<0.001)。Figure 2. Effects of No.0 diesel oil emulsion on hepatopancreatic antioxidant enzyme activities of L. vannameiNote: * represents significant differences among the groups (p<0.05); ** represents very significant differences among the groups (p<0.01); *** represents extremely significant differences among the groups (p<0.001).2.2 相关功能基因表达分析
随LH原油胁迫时间的延长,中、高浓度组SOD、CAT和GPx基因表达量逐渐上升并显著高于对照组 (p<0.05);中、高浓度组LZM、EF、proPO表达量呈先上升后下降的趋势;各浓度组caspase2、Prx表达量均显著低于对照组 (p<0.05,图3)。
图 3 不同浓度LH原油乳化液对凡纳滨对虾肝胰腺功能基因表达的影响注:* 表示组间存在显著性差异 (p<0.05);** 表示组间存在非常显著性差异 (p<0.01);*** 表示组间存在极显著性差异 (p<0.001)。Figure 3. Effects of LH crude oil emulsion on hepatopancreatic relative gene expression of L. vannameiNote: * represents significant differences among the groups (p<0.05); ** represents very significant difference among the groups (p<0.01); *** represents extremely significant differences among the groups (p<0.001).随0# 柴油胁迫时间的延长,各浓度组SOD表达量呈先升高后下降的变化趋势;高浓度组CAT、LZM、caspase2、EF表达量呈下降趋势;低、高浓度组Prx和各浓度组proPO表达量呈先下降后升高的趋势;至胁迫第15天,SOD、CAT、LZM、caspase2表达量被显著抑制 (p<0.05),而GPx、Prx、proPO和低、中浓度组EF表达量均显著高于对照组 (p<0.05,图4)。
图 4 不同浓度0# 柴油乳化液对凡纳滨对虾肝胰腺功能基因表达的影响注:* 表示组间存在显著性差异 (p<0.05);** 表示组间存在非常显著性差异 (p<0.01);*** 表示组间存在极显著性差异 (p<0.001)。Figure 4. Effects of NO.0 diesel oil emulsion on hepatopancreatic relative gene expression of L. vannameiNote: * represents significant differences among the groups (p<0.05); ** represents very significant differences among the groups (p<0.01); *** represents extremely significant differences among the groups (p<0.001).2.3 综合生物标志物响应 (IBR) 变化
不同浓度LH原油和0# 柴油胁迫下,凡纳滨对虾肝胰腺生物标志物IBRv2结果如图5所示。胁迫第1天时,高浓度组均具有最大覆盖面积;胁迫第7天时,中浓度组均具有最大覆盖面积;胁迫第15天时,原油中浓度组覆盖面积最大,柴油低浓度组覆盖面积最大。
图 5 LH原油和0# 柴油乳化液胁迫下凡纳滨对虾生物标志物IBRv2分析注:抗氧化酶指标包括SOD、CAT、POD、GST、GPx及GSH;基因表达指标包括SOD、CAT、GPx、LZM、caspase2、Prx、EF及proPO;a—c分别代表LH原油胁迫组第1、第7和第15天各生物标志物的IBRv2分析;d—f分别代表0# 柴油胁迫组第1、第7和第15天各生物标志物的IBRv2分析。Figure 5. IBRv2 analysis of biomarkers in L. vannamei after exposure to emulsion of LH crude oil and No.0 diesel oilNote: The antioxidant enzyme indicators include SOD, CAT, POD, GST, GPx and GSH, and the gene expression indicators include SOD, CAT, GPx, LZM, caspase2, Prx, EF and proPO. a−c and d−f represent the IBRv2 analysis of the biomarkers in the LH crude oil and NO.0 diesel oil exposure groups on 1st, 7th and 15th day, respectively.LH原油胁迫第1天时,低浓度组对Prx表达量有较大影响,中、高浓度组对EF、CAT表达量影响较大;第7天时,EF表达量在低浓度组中的响应水平较高,POD活性在中浓度组中的响应水平较高,CAT活性及其基因表达量在高浓度组中的响应水平较高;第15天时,GPx活性及其基因表达量在各浓度组中均有较高的响应水平。在0# 柴油胁迫的3个时间节点中,各浓度组GPx活性及其基因表达量均有较高的响应水平;胁迫第1天时,CAT表达量在低浓度组中的响应水平较高;第15天时,高浓度组caspase2基因表达量和CAT活性均有较高的响应水平。
对14种生物标志物的偏差系数绝对值进行求和,得到IBRv2值 (图6)。随LH原油和0# 柴油胁迫浓度的增加和时间的延长,IBRv2值呈剂量-效应和时间-效应关系。在胁迫的3个时间点中,其IBRv2值大小均为低浓度组<中浓度组<高浓度组;在相同浓度的LH原油和0# 柴油胁迫下,其IBRv2值随胁迫时间的延长呈升高趋势,均表现为第1天<第7天<第15天,且柴油组生物标志物基因表达量的IBRv2值变化范围均大于原油组。
2.4 生物标志物PCA分析
对LH原油胁迫下的14个生物标志物进行PCA分析 (图7-a),主成分1 (PC1) 和主成分2 (PC2) 分别解释了总变异的46.81%和15.32%。其中,PC1与SOD、CAT、GPx的活性和表达量具有较强的相关性,上述抗氧化酶活性及其基因表达量与中、高浓度组LH原油的长期胁迫 (15 d) 具有较强的相关性;PC2与EF、caspase2基因表达有关,说明低浓度组长期胁迫及高浓度组短期胁迫 (1 d) 与延伸因子、细胞凋亡因子有较强的相关性。
图 7 LH原油 (a) 和0# 柴油 (b) 乳化液胁迫下凡纳滨对虾肝胰腺中生物标志物的PCA分析注:抗氧化酶指标包括SOD、CAT、POD、GST、GPx和GSH;基因表达指标包括SOD、CAT、GPx、LZM、caspase2、Prx、EF及proPO。图中不同颜色的点分别代表各胁迫组,例如7-a中的红点“A15d”代表低浓度组 (A组,25 μg·L−1) 原油胁迫凡纳滨对虾15 d,其他颜色的点以此类推。Figure 7. Results of PCA analysis of biomarkers in hepatopancreas of L. vannamei exposed to emulsion of LH crude oil (a) and No.0 diesel oil (b)Note: The antioxidant enzyme indicators include SOD, CAT, POD, GST, GPx and GSH, the gene expression indicators included SOD, CAT, GPx, LZM, caspase2, Prx, EF and proPO. The different colored dots represent the stress groups. For example, the red dot "A15d" in Fig. 7-a means that the low concentration group (group A, 25 μg·L−1) of crude oil was stressed for 15 d, and the other colored dots are similar.对0# 柴油胁迫下的14个生物标志物进行PCA分析 (图7-b),PC1解释了40.56%,与SOD、CAT、GSH活性和caspase2、proPO表达量的相关性较强,表明低、中浓度组短期胁迫和高浓度组长期胁迫对抗氧化酶活性、酚氧化酶原系统和细胞凋亡基因的影响较大;PC2解释了19.02%,与POD活性和GPx、SOD表达量相关性较强,表明中浓度组长期胁迫 (7 d) 和高浓度组短期胁迫对抗氧化酶的影响较大。
3. 讨论
3.1 LH原油和0# 柴油对凡纳滨对虾肝胰腺抗氧化酶活性的影响
在环境胁迫下,生物机体会通过氧化代谢产生过量的活性氧自由基 (Reactive oxygen species, ROS),当ROS含量超过机体的清除能力将引发氧化应激反应,从而损伤细胞膜、使酶失活,并造成遗传物质等重要细胞成分氧化损伤,进一步诱导细胞凋亡[22]。研究表明,肝胰腺由于在消化、物质储备、脂质代谢、解毒等方面具有较高的代谢率,因此在环境胁迫中更容易受到氧化应激的影响[23]。
本研究中,SOD、CAT活性随着LH原油浓度的增加和胁迫时间的延长而逐渐上升,可能是由于机体在LH原油长期胁迫下,应激生成大量SOD、CAT用于清除过量的ROS,以减轻污染物胁迫所带来的伤害。相似的结果在李磊等[24]的研究中也有体现,在0.45和0.05 μg·L−1苯并[a]芘 (BaP) 胁迫下,脊尾白虾 (Exopalaemon carinicauda) 的SOD活性均逐渐上升,在第15天时2个实验组的SOD活性分别为第0天的3.04和2.38倍;王琪等[25]的研究也发现,随着菲胁迫浓度的增加,各浓度组泥鳅 (Misgurnus anguillicaudatus) 肝脏中的CAT活性显著上升。暴露于0# 柴油乳化液中第1天时,SOD、CAT活性与对照组无显著性差异,但随着胁迫时间延长和浓度增加,SOD、CAT活性逐渐降低,可能是因为在柴油胁迫初期,ROS积累远高于抗氧化酶的清除效率,造成机体损伤[26]。相似的结果也见于黑鲷 (Sparus macrocephalus) 肝脏[27],SOD和CAT活性随0# 柴油浓度的增加和胁迫时间的延长呈下降趋势,可能是由于ROS积累使SOD和CAT的酶蛋白结构、底物亲和力等受到损伤[28]。
研究表明,生物体在轻度胁迫时POD活性通常会升高,但长期胁迫会导致氧化损伤并引发细胞病变[29]。本研究中,随着LH原油胁迫时间的延长,各浓度组POD活性呈先升后降的变化趋势,说明原油短期胁迫可激活凡纳滨对虾的抗氧化防御系统,POD活性增强以减少ROS对机体的损伤,这与李传慧[30]的研究结果一致;但高浓度组的活性比中浓度组要低,可能是因为高浓度组胁迫超出了机体的生理调节范围,对机体的刺激作用增强[31]。与原油不同,0# 柴油长期胁迫时,各浓度组POD活性均被显著抑制,说明0# 柴油长期胁迫使凡纳滨对虾肝胰腺细胞严重损伤,导致其氧化应激能力下降。
从剂量-效应和时间-效应关系分析,LH原油胁迫使GST活性逐渐上升,以减少ROS的伤害,表明机体对胁迫产生了适应性反应。而高浓度组柴油的GST活性先升后降,可能是因为在相II反应酶系中,GST与各种亲电性化合物或某些相I代谢产物结合,并催化GSH反应生成水溶性的GSH硫结合物,起到脱毒作用。但高浓度组长期胁迫削弱了结合能力,破坏氧化还原平衡而抑制了GST活性。相似的结果在张先勇等[32]和任加云等[33]的研究中也有体现,BaP长期胁迫马氏珠母贝 (Pinctada martensii) 导致GST活性降低,多环芳烃 (PAHs) 长时间胁迫栉孔扇贝也会抑制其GST活性。
在本研究中,GPx活性与2种石油污染物之间呈显著的时间-效应和剂量-效应正相关关系。GSH含量在石油类短期胁迫中显著上升,可能是因为GPx催化GSH将过氧化氢 (H2O2) 分解为水,使得GPx活性增加,从而加速GSH对H2O2 的分解以清除过量的ROS[34]。随着柴油胁迫时间的延长,GSH含量显著下降,可能是由于GSH与石油烃中的有毒物质结合解毒,或是GSH在清除ROS时被GPx催化为了二硫化谷胱甘肽 (Glutathione disculfide, GSSG)[35-36]。
从2种石油对SOD、CAT、POD活性的时间影响来看,柴油各浓度组的最大诱导出现的时间 (第1天) 要早于原油 (第7和第15天),且在胁迫15 d后,柴油组抗氧化酶 (除GPx外) 的活性基本低于对照组,0# 柴油对酶活性的影响更显著,这与黄南建等[15]的研究结果一致。2种石油的毒性差异可能是由于其所含组分不同;通常,柴油中的饱和烷烃和芳香烃占比均高于原油[37],因此0# 柴油对凡纳滨对虾肝胰腺抗氧化酶活性的影响大于LH原油。
3.2 LH原油和0# 柴油对凡纳滨对虾肝胰腺功能基因表达的影响
本研究中,凡纳滨对虾在低浓度组石油污染物胁迫下,SOD和CAT基因表达的显著诱导时间 (1 d) 早于其酶活性,且在高浓度组LH原油胁迫下,抗氧化酶的基因表达量诱导倍数也大于酶活性的诱导倍数,说明相关功能基因表达量作为生物标志物可能具有更高的敏感性。在镉 (Cd)胁迫暗纹东方鲀 (Takifugu obscurus)[38]和锌离子 (Zn2+) 胁迫花斑裸鲤 (Gymnocypris eckloni)[39]的研究中也有体现,SOD、CAT、GPx表达量在肝脏中显著升高,表明相关功能基因表达量可能是氧化应激反应的重要生物标志物。
氧化应激会改变caspase家族等凋亡相关基因的表达,从而诱导细胞凋亡[40],其中的caspase2是关键基因。在本研究中,caspase2表达量在2种石油污染物胁迫下均呈现先升高后下降的趋势,高浓度组显著低于对照组,可能是胁迫初期机体上调基因表达应对初步损伤并诱发细胞凋亡[41-42],随着胁迫时间延长和毒素浓度增加,细胞的应激响应被抑制,caspase2表达量下降可能意味着细胞凋亡的终止或损伤修复机制的激活。本研究中,EF表达量在LH原油胁迫下先升高后下降,可能是胁迫初期机体诱导EF表达以提高蛋白质的合成速率,修复氧化损伤,但因长期胁迫无法维持正常代谢和生理状态,导致EF表达量下降。在0# 柴油胁迫下,低、中浓度组EF表达量均显著升高,说明EF可能参与了响应氧化应激的凋亡程序,通过肌动蛋白束缚等方式保护细胞免受caspase家族介导的细胞凋亡[43]。因此,推测caspase2表达量的下降是由于EF的抑制作用;高浓度组可能导致过量的细胞毒性,损伤细胞的代谢和蛋白合成路径。
在原油和柴油的胁迫下,各浓度组LZM表达量呈先升高后下降的趋势,这与柴油胁迫虹鳟 (Oncorhynchus mykiss)[44]的研究结果一致,可能是胁迫初期凡纳滨对虾增加LZM表达量以应答环境胁迫;随胁迫时间的延长,机体生理功能紊乱,免疫能力降低,LZM表达量下降。在原油胁迫下,各浓度组proPO表达量呈先升高后下降的趋势,低浓度组表达量在暴露第7天时达到峰值,高浓度组的表达量显著低于对照组。表明凡纳滨对虾在受原油胁迫时可以激活proPO系统,产生大量的酚氧化酶 (Phenoloxidase, PO) 应对环境变化,但随胁迫时间的延长,过量ROS积累导致细胞组织受损,proPO表达量下降。这与盐碱胁迫中华绒螯蟹 (Eriocheir sinensis) 的变化趋势相似[45]。而在柴油胁迫下,proPO表达量先降再升,可能是因为柴油直接抑制了proPO,导致细胞应激反应机制暂时失调,持续胁迫后机体产生了适应并激活了免疫系统。
本研究中,原油胁迫下的Prx表达量显著低于对照组,可能是因为原油中的PAHs损伤了相关的信号传递路径,破坏了细胞的抗氧化防御系统。而柴油胁迫下的Prx表达量呈上升趋势,可能是细胞通过增加Prx表达量应对氧化应激,保护自身免受损伤。相似的结果在Zhuang等[46]的研究中也有体现,将桡足类分别暴露在2种PAHs中,1, 2-NAPH胁迫下Prx表达量显著上调,芘胁迫下则显著下调,说明生物对不同的PAHs胁迫表现出敏感性差异,会导致不同的基因表达模式。
3.3 LH原油和0# 柴油暴露对凡纳滨对虾毒性效应的IBR及PCA分析
有研究表明,一般情况下,基因表达通常在氧化应激早期迅速变化,酶活反应则较为缓慢[47]。本研究中,SOD表达量在低浓度LH原油胁迫初期已显著上升,而SOD活性在高浓度胁迫下才显著上升。在低氧胁迫鲢的研究中也有相似的结果[48],SOD表达量在浮头组已显著上升,而SOD活性在窒息组才显著下降。说明基因表达变化可能比酶活性变化更早发生,而酶活性显著变化可能在较高浓度或长时间胁迫后才出现。
本研究的IBRv2分析结果显示,石油污染物短期胁迫凡纳滨对虾时,高浓度组对生物标志物的影响最大;随胁迫时间的延长,低、中浓度组柴油对凡纳滨对虾的影响更大,原油组则是中浓度组影响最大,高浓度组次之。在0硫柴油和平湖原油WAF分别胁迫缢蛏的研究中发现,2种油品的毒性差异源于组分不同,且柴油毒性大于原油[12]。因此,本研究中高浓度组柴油短期胁迫可能已对对虾造成了不可逆转的氧化损伤,而低、中浓度组长时间胁迫可能损伤细胞结构和功能。随2种油品胁迫浓度的增加和时间的延长,IBRv2值呈剂量-效应和时间-效应关系,类似结果在菲胁迫泥鳅的研究中也有体现[25],其综合毒性应激水平随着菲的暴露水平和时间的增加而上升。柴油组生物标志物基因表达量的IBRv2值变化范围均大于原油组,说明其表现出了更高的响应水平。Xu等[49]和Lu等[50]的研究也表明,生物标志物基因表达量在胁迫初期就显著变化,早于酶活性。
PCA分析结果显示,LH原油和0# 柴油的长期和短期胁迫均显著影响抗氧化酶活性及其基因表达,同时也影响细胞凋亡相关基因的表达。但2种油品所影响的生物标志物存在差异:低浓度组LH原油长期胁迫和高浓度组短期胁迫显著影响EF,而低、中浓度组0# 柴油短期胁迫和高浓度组长期胁迫显著影响proPO系统。可能是由于0# 柴油含有更多低分子量的化合物和芳香烃,毒性大于LH原油[15];EF通常与维持细胞蛋白质稳定性相关,其表达量增加是促氧化剂诱导细胞凋亡的早期事件[43],LH原油中的芳香烃等应激因素可能会引起蛋白质氧化等细胞损伤,导致EF表达量变化[51];而在PAHs类污染物屈 (Chrysene) 胁迫斑节对虾 (Penaeus monodon) 实验[52]中,proPO表达量显著下调,表明PAHs影响生物的生长和免疫系统。
研究表明,测定生物标志物CAT和GPx等的酶活性,可以显著区分不同养殖模式污染下的虾类[53]。在本研究中,PCA与IBRv2的分析结果相似,CAT和GPx的活性及表达量在2种油品胁迫下均受到显著影响,因此推测CAT和GPx的活性及表达量可以作为预测LH原油和0# 柴油污染暴露的早期预警生物标志物。
4. 结论
在LH原油和0# 柴油胁迫下,凡纳滨对虾肝胰腺抗氧化酶活性和抗氧化防御、细胞凋亡及免疫相关基因的表达均随油品浓度增加和胁迫时间延长而显著变化,其中0# 柴油的影响更为显著,表明凡纳滨对虾对0# 柴油的敏感性高于LH原油。IBRv2值的变化也显示出不同浓度和时间下的生物毒性差异,说明该指数可作为评估不同油品污染对凡纳滨对虾综合毒性效应的有效工具。此外,CAT和GPx的活性及表达量与油品暴露有较强的相关性,推测其可作为预测石油污染暴露的早期预警指标。
-
表 1 正交试验设计表
Table 1 Orthogonal design table
编号
No.ρ(氮)/mg·L-1
Nρ(磷)/10-2 mg·L-1
Pρ(铁)/10-2 mg·L-1
Fe3+1 2.47 (1) 0.99 (1) 0.52 (1) 2 2.47 (1) 4.96 (2) 1.68 (2) 3 2.47 (1) 9.93 (3) 2.58 (3) 4 4.94 (2) 0.99 (1) 1.68 (2) 5 4.94 (2) 4.96 (2) 2.58 (3) 6 4.94 (2) 9.93 (3) 0.52 (1) 7 7.42 (3) 0.99 (1) 2.58 (3) 8 7.42 (3) 4.96 (2) 0.52 (1) 9 7.42 (3) 9.93 (3) 1.68 (2) 表 2 不同氮质量浓度对绿色颤藻叶绿素a和藻体质量浓度的影响
Table 2 Effect of N concentrations on chlorophyll-a contents and dry weight of O.chlorine
编号
No.ρ(氮)/mg·L-1
Nρ(叶绿素a)/μg·L-1
Chl-a叶绿素a增加率/%
increase rate of Chl-aρ(藻体)/g·L-1
dry weight增重率/%
weight gain rate1 0 27.82±7.60a -0.48 0.16±0.01a 0.29 2 0.01 46.53±10.47ab -0.13 0.16±0.01a 0.29 3 0.06 65.75±11.24b 0.24 0.16±0.02a 0.24 4 0.25 71.23±8.89b 0.34 0.18±0.02a 0.41 5 1.24 69.33±9.06b 0.30 0.16±0.03a 0.29 6 2.47 69.16±11.16b 0.30 0.17±0.01a 0.35 7 6.18 99.69±14.21c 0.88 0.19±0.01a 0.47 8 12.36 227.75±15.87d 3.28 0.30±0.03b 1.35 9 24.72 294.44±22.82e 4.54 0.36±0.04c 1.83 10 61.80 297.73±19.98e 4.60 0.36±0.02c 1.83 11 247.21 292.18±25.38e 4.50 0.35±0.03c 1.76 12 494.42 288.73±7.14e 4.43 0.24±0.09c 0.91 表 3 不同磷质量浓度对绿色颤藻叶绿素a和藻体质量浓度的影响
Table 3 Effect of P concentrations on chlorophyll-a contents and dry weight of O.chlorine
编号
No.ρ(磷)/mg·L-1
Pρ(叶绿素a)/μg·L-1
Chl-a叶绿素a增加率/%
increase rate of Chl-aρ(藻体)/g·L-1
dry weight增重率/%
weight gain rate1 0 44.21±3.64a -0.17 0.24±0.02a 0.94 2 9.93×10-4 51.39±16.12ab -0.03 0.26±0.07ab 1.06 3 1.99×10-3 52.45±16.87ab -0.01 0.29±0.04abc 1.29 4 9.93×10-3 66.79±3.44ab 0.26 0.31±0.04bcd 1.47 5 1.99×10-2 74.95±15.58b 0.41 0.34±0.01cde 1.71 6 9.93×10-2 104.95±19.05c 0.97 0.34±0.01cde 1.71 7 0.20 115.77±2.05cd 1.18 0.36±0.04def 1.83 8 0.99 120.43±9.49cd 1.27 0.37±0.03defg 1.94 9 1.99 164.68±10.57e 2.10 0.38±0.02efg 2.00 10 4.96 139.39±14.96d 1.62 0.42±0.03fg 2.29 11 9.93 137.84±25.95d 1.59 0.42±0.02g 2.35 12 19.85 120.11±19.70cd 1.26 0.42±0.03fg 2.29 表 4 不同铁质量浓度对绿色颤藻叶绿素a和藻体质量浓度的影响
Table 4 Effect of Fe3+ concentrations on chlorophyll-a contents and dry weight of O.chlorine
编号
No.ρ(铁)/ mg·L-1
Fe3+ρ(叶绿素a)/μg·L-1
Chl-a叶绿素a增加率/%
increase rate of Chl-aρ(藻体)/g·L-1
dry weight增重率/%
weight gain rate1 0 72.50±3.25a 0.36 0.22±0.01a 0.71 2 5.16×10-5 74.25±1.24ab 0.40 0.22±0.02ab 0.76 3 2.58×10-4 78.93±1.95abc 0.48 0.24±0.03abc 0.88 4 5.16×10-3 82.56±8.36abc 0.55 0.27±0.01cde 1.18 5 2.58×10-2 90.43±1.17bcd 0.70 0.28±0.01de 1.23 6 5.16×10-2 103.57±4.37de 0.95 0.30±0.05def 1.35 7 0.26 117.57±7.54ef 1.21 0.33±0.02f 1.64 8 0.52 121.77±7.81f 1.29 0.34±0.01f 1.71 9 2.58 121.48±3.72f 1.29 0.33±0.03f 1.59 10 25.81 113.92±3.95ef 1.14 0.31±0.02ef 1.47 11 51.67 91.70±1.07cd 0.73 0.26±0.03bcd 1.06 表 5 绿色颤藻叶绿素a质量浓度的正交试验结果
Table 5 Chlorophyll-a contents result of orthogonal design test of O.chlorine
编号
No.ρ(氮)/mg·L-1
Nρ(磷)/10-2 mg·L-1
Pρ(铁)/10-2 mg·L-1
Fe3+ρ(叶绿素a)/μg·L-1
content of Chl-a1 2.47 (1) 0.99 (1) 0.52 (1) 25.81 2 2.47 (1) 4.96 (2) 1.68 (2) 18.87 3 2.47 (1) 9.93 (3) 2.58 (3) 22.33 4 4.94 (2) 0.99 (1) 1.68 (2) 50.54 5 4.94 (2) 4.96 (2) 2.58 (3) 48.77 6 4.94 (2) 9.93 (3) 0.52 (1) 47.15 7 7.42 (3) 0.99 (1) 2.58 (3) 45.21 8 7.42 (3) 4.96 (2) 0.52 (1) 51.04 9 7.42 (3) 9.93 (3) 1.68 (2) 60.50 表 6 绿色颤藻藻体质量浓度的正交试验结果
Table 6 Dry weight result of orthogonal design test of O.chlorine
编号
No.ρ(氮)/mg·L-1
Nρ(磷)/10-2 mg·L-1
Pρ(铁)/10-2 mg·L-1
Fe3+ρ(藻体)/ g·L-1
dry weight1 2.47 (1) 0.99 (1) 0.52 (1) 0.11 2 2.47 (1) 4.96 (2) 1.68 (2) 0.11 3 2.47 (1) 9.93 (3) 2.58 (3) 0.09 4 4.94 (2) 0.99 (1) 1.68 (2) 0.11 5 4.94 (2) 4.96 (2) 2.58 (3) 0.10 6 4.94 (2) 9.93 (3) 0.52 (1) 0.12 7 7.42 (3) 0.99 (1) 2.58 (3) 0.11 8 7.42 (3) 4.96 (2) 0.52 (1) 0.13 9 7.42 (3) 9.93 (3) 1.68 (2) 0.13 -
[1] 马祖友. 蓝藻的生长生理特征及其竞争优势研究[D]. 杨凌: 西北农林科技大学, 2005. [2] 戴瑾瑾, 陈德辉, 高云芳, 等. 蓝藻毒素的研究概况[J]. 武汉植物学研究, 2009, 27(1): 90-97. doi: 10.3969/j.issn.2095-0837.2009.01.014 [3] 王欣伊, 阚振荣, 王梅梅. 淡水藻类产毒研究进展[J]. 生物学杂志, 2005, 22(2): 5-9. doi: 10.3969/j.issn.2095-1736.2005.02.002 [4] 查广才, 麦雄伟, 周昌清, 等. 凡纳滨对虾低盐度养殖池浮游藻类群落研究[J]. 海洋水产研究, 2006, 27(1): 1-7. doi: 10.3969/j.issn.1000-7075.2006.01.001 [5] 彭聪聪, 李卓佳, 曹煜成, 等. 凡纳滨对虾半集约化养殖池塘浮游微藻优势种变动规律及其对养殖环境的影响[J]. 海洋环境科学, 2011, 30(2): 193-198. doi: 10.3969/j.issn.1007-6336.2011.02.010 [6] BRITTAIN S, MOHAMED Z A, WANG J, et al. Isolation and characterization of microcystins from a River Nile strain of Oscillatoria tenuis Agardh ex Gomont[J]. Toxicon, 2000, 38(12): 1759-1771. doi: 10.1016/S0041-0101(00)00105-7
[7] 金星. 藻类代谢产物中典型异味物质及其变化规律研究[D]. 南京: 南京理工大学, 2009. [8] 曹平, 黄翔鹄, 李长玲, 等. 颤藻对凡纳滨对虾生长和免疫相关酶活力的影响[J]. 渔业现代化, 2011, 38(5): 25-30. doi: 10.3969/j.issn.1007-9580.2011.05.006 [9] 梁伟峰, 陈素文, 李卓佳, 等. 虾池常见微藻种群温度、盐度和氮、磷含量生态[J]. 应用生态学报, 2009, 20(1): 223-227. [10] 王颜萍, 何鹰, 郑晓玲, 等. 铁限制条件下东海原甲藻分泌铁载体[J]. 生态学报, 2005, 25(10): 2788-2791. doi: 10.3321/j.issn:1000-0933.2005.10.047 [11] 蒲新明. 长江口区浮游植物营养限制因子的研究[D]. 青岛: 中国科学院海洋研究所, 2000. [12] SANO T, KAYA K. A 2-amino-2-butenoic acid (Dhb)-containing microcystin isolated from Oscillatoria agardhii[J]. Tetrahedron Lett, 1995, 36(47): 8603-8606. doi: 10.1016/0040-4039(95)01824-22095-1736.2005.02.002
[13] FUJⅡ K, SIVONEN K, NAGANAWA E, et al. Non-toxic peptides from toxic cyanobacteria, Oscillatoria agardhii[J]. Tetrahedron, 2000, 56(5): 725-733. doi: 10.1016/S0040-4020(99)01017-0
[14] JIA W, HUANG X H, LI C L. A preliminary study of the algicidal mechanism of bioactive metabolites of Brevibacillus laterosporus on Oscillatoria in prawn ponds[J]. Sci World J, 2014, doi: 10.1155/2014/869149.
[15] 贾雯, 黄翔鹄, 李长玲, 等. 侧孢短芽孢杆菌溶藻活性代谢产物对虾池颤藻溶藻效果研究[J]. 水产学报, 2013, 27(3): 465-472. https://www.cqvip.com/QK/90183X/201303/45244554.html [16] 贺春花, 黄翔鹄, 李长玲, 等. 温度、光照度、盐度和pH对颤藻生长的限制条件研究[J]. 渔业现代化, 2011, 38(6): 20-25. doi: 10.3969/j.issn.1007-9580.2011.06.005 [17] 寻涛. 高锰酸钾预氧化混凝去除水中颤藻的试验研究[J]. 安全与环境工程, 2009, 16(4): 45-48. doi: 10.3969/j.issn.1671-1556.2009.04.012 [18] MILLIE D F, GREENE D A, JOHNSEN P B. Effects of the carotenoid-inhihiting herbicide, flufidone on Oscillatoria agardhii Gomont(cyanobacteria)[J]. Aquat Toxicol, 1990, 16(1): 41-52. doi: 10.1016/0166-445X(90)90076-2
[19] 梁伟峰, 李卓佳, 陈素文, 等. 对虾养殖池塘微藻群落结构的调查与分析[J]. 南方水产, 2007, 3(5): 33-39. doi: 10.3969/j.issn.2095-0780.2007.05.006 [20] 刘孝竹, 李卓佳, 蓸煜成, 等. 低盐度养殖池塘浮游微藻的种类组成、数量及优势种群变动[J]. 南方水产, 2009, 5(1): 9-16. doi: 10.3969/j.issn.1673-2227.2009.01.002 [21] 曹煜成, 王少沛, 李卓佳, 等. 微绿球藻、隐藻和颤藻竞争氮磷营养的通径分析[J]. 渔业科学进展, 2012, 33(4): 99-106. doi: 10.3969/j.issn.1000-7075.2012.04.013 [22] WILBUR K M, ANDERSON N G. Electrometric and colorimetric determination of carbonic Anhydrase[J]. Biol Chem, 1948, 176(1): 147-154. doi: 10.1016/S0021-9258(18)51011-5
[23] 国家环保总局, 《水和废水检测分析方法》编委会. 水和废水监测分析方法[M]. 4版. 北京: 中国环境科学出版社, 2002: 670-671. [24] 石岩峻. 赤潮藻对营养盐的吸收及生长和相关特性研究[D]. 北京: 北京化工大学, 2004. [25] 李志伟, 崔力拓, 齐凤生, 等. 铜绿微囊藻生长特征及营养盐对其生长的影响[J]. 水利渔业, 2006, 26(1): 65-68. doi: 10.3969/j.issn.1003-1278.2006.01.031 [26] 江燕, 甘旭华, 唐欣昀, 等. 氮、磷营养因子对赤潮异弯藻生长的影响[J]. 应用生态学报, 2006, 17(3): 557-559. doi: 10.3321/j.issn:1001-9332.2006.03.040 [27] 黄翔鹄. 对虾高位池水环境养殖污染和浮游微藻生态调控机制研究[D]. 上海: 东华大学, 2013. [28] YUSOFF F M, ZUBAIDAH M S, MATIAS H B, et al. Phytoplankton succession in intensive marine shrimp culture ponds treated with a commercial bacterial product[J]. Aquac Res, 2002, 33(4): 269-278. doi: 10.1046/j.1355-557x.2002.00671.x
[29] 陆珠润, 蒋霞敏, 陆艳晶. 氮、磷、铁对原绿球藻生长的影响[J]. 宁波大学学报: 理工版, 2009, 22(3): 326-331. doi: 10.3969/j.issn.1001-5132.2009.03.006 [30] KRUSKOPF M, DU P S. Growth and filament length of the bloom forming Oscillatoria simplicissima (Oscillatoriales, Cyanophyta) in varying N and P concentrations[J]. Hydrobiologia, 2006, 556(1): 357-362. doi: 10.1007/s10750-005-1061-0
[31] 徐景亮. 东海原甲藻的氮、磷、铁、锰营养生态位研究[D]. 广州: 暨南大学, 2005. [32] MARTIN J H, FITZWATER S E. Iron deficiency limits phytoplankton growth in the northeast Pacific subartic[J]. Nature, 1988, 331(6154): 341-343. doi: 10.1038/331341a0
[33] 孙颖颖. 环境因子对球等鞭金藻生长的影响[D]. 大连: 大连理工大学, 2007. [34] 王扬才. 氮、磷、铁营养盐浓度对牟氏角毛藻生长的影响[J]. 海洋渔业, 2006, 28(2): 173-176. doi: 10.3969/j.issn.1004-2490.2006.02.016 [35] 吕秀平, 胡晗华, 张栩, 等. Fe3+对浮游颤藻生长和光合作用的影响[J]. 水生生物学报, 2005, 29(3): 318-322. doi: 10.3321/j.issn:1000-3207.2005.03.014 [36] 刘静. Fe3+对富营养化水体中三种常见淡水藻类生长的影响[D]. 扬州: 扬州大学, 2008. -
期刊类型引用(4)
1. 尹海晓, 王志康, 罗小勇, 王铭, 郭云, 李凯. 不同生长条件对藻源有机物化学特征的影响. 给水排水. 2019(02): 134-140 . 百度学术
2. 曾玲, 龙超, 文菁. 营养盐限制对利玛原甲藻生长和产毒的影响. 广东农业科学. 2018(02): 135-144 . 百度学术
3. 徐煜, 徐武杰, 文国樑, 胡晓娟, 苏浩昌, 杨铿, 李卓佳, 曹煜成. 颤藻浓度和水温对凡纳滨对虾响应颤藻粗提液毒性的影响. 南方水产科学. 2017(01): 26-32 . 本站查看
4. 孙志伟, 邱丽华, 段舜山, 蔡卓平, 李卓佳, 文国樑. 化感作用抑制有害藻类生长的研究进展. 生态科学. 2015(06): 188-192 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 3462
- HTML全文浏览量: 125
- PDF下载量: 1669
- 被引次数: 5