广西北部湾沿海市售贝类脂溶性毒素污染分析及食用风险评价

江姗姗, 张晨晓

江姗姗, 张晨晓. 广西北部湾沿海市售贝类脂溶性毒素污染分析及食用风险评价[J]. 南方水产科学, 2023, 19(4): 158-167. DOI: 10.12131/20230024
引用本文: 江姗姗, 张晨晓. 广西北部湾沿海市售贝类脂溶性毒素污染分析及食用风险评价[J]. 南方水产科学, 2023, 19(4): 158-167. DOI: 10.12131/20230024
JIANG Shanshan, ZHANG Chenxiao. Pollution analysis and dietary exposure risk assessment of lipophilic toxins in shellfish from Beibu Gulf seafood market in Guangxi Province[J]. South China Fisheries Science, 2023, 19(4): 158-167. DOI: 10.12131/20230024
Citation: JIANG Shanshan, ZHANG Chenxiao. Pollution analysis and dietary exposure risk assessment of lipophilic toxins in shellfish from Beibu Gulf seafood market in Guangxi Province[J]. South China Fisheries Science, 2023, 19(4): 158-167. DOI: 10.12131/20230024

广西北部湾沿海市售贝类脂溶性毒素污染分析及食用风险评价

基金项目: 广西重点研发计划项目 (桂科AB19110020)
详细信息
    作者简介:

    江姗姗 (1995—),女,硕士研究生,研究方向为海洋贝类毒素。E-mail: jss15056254097@163.com

    通讯作者:

    张晨晓 (1975—),女,教授,博士,研究方向为海产品加工与质量安全。E-mail: zhangcx@bbgu.edu.cn

  • 中图分类号: TS 254.2

Pollution analysis and dietary exposure risk assessment of lipophilic toxins in shellfish from Beibu Gulf seafood market in Guangxi Province

  • 摘要: 脂溶性贝类毒素分布广泛,可通过食物链传递,对消费者健康构成巨大威胁。广西北部湾海域是我国重要的海水养殖区,其贝类产量占比超过一半。为评估广西北部湾沿海居民因摄食贝类而引起的脂溶性贝类毒素摄入风险,并为建立毒素限量标准提供科学依据,采用液相色谱串联质谱法 (LC-MS/MS) 分析脂溶性贝类毒素含量,24 h 膳食回顾法调查居民连续7 d的海产品进食情况,应用海洋生物毒素点评估方法,计算膳食中每日毒素摄入量 (Estimated daily intake, ESI),对照欧洲食品安全局推荐的急性参考剂量 (Acute reference dose, ARfD) 值,评价其食用安全性。结果显示:检出的4种毒素含量均低于欧盟规定的安全限量值。检出率由高到低依次为环亚胺类毒素 (Gymnodimine, GYM)(83.02%)、大田软海绵酸 (Okadaic acid, OA)(51.16%)、鳍藻毒素 (Dinophysis toxin, DTX-2)(40.91%) 和虾夷扇贝毒素 (Homo-yessotoxin, Homo-YTX)(8.6%),检出的最高质量分数分别为105.4、31.39、38.19和159.66 μg·kg−1;毒素含量存在季节和物种差异,牡蛎秋季GYM含量最高,冬季OA含量最高。DTX-2更易在秋冬季样品中检出,Homo-YTX仅在夏秋季扇贝和贻贝样品中检出。膳食调查结果显示当地居民贝肉日平均消费量为45 g·d−1。以当地居民贝肉食用量和人体质量估算,OA组毒素ESI值为0.26 μg·(kg·d)−1,小于欧盟规定的ARfD值 [0.3 μg·(kg·d)−1]。但若以欧盟提供的贝肉食用量和人体质量数据估算,ESI值 [0.34 μg·(kg·d)−1] 则高于ARfD值。研究表明,虽然广西北部湾沿海市售贝类产品中脂溶性贝类毒素检出值低于欧盟规定的安全限量值,但当地居民存在OA组毒素膳食暴露风险。
    Abstract: Lipophilic shellfish toxins are widely distributed and will be transmitted to consumers through the food chain, posing a great threat to human health. The Beibu Gulf in Guangxi Province is an important mariculture area in China, accounting for more than half of the national output. To estimate the potential dietary exposure risk of coastal residents in that area, and to provide a scientific basis for the establishment of toxin limit standards, we analyzed the contents of lipophilic shellfish toxins by liquid chromatography tandem mass spectrometry. We have used a 24-hour dietary review method to investigate the consumption of marine products by coastal residents for seven consecutive days. Then we calculated the value of estimated daily intake (ESI) of toxins in the diet by using the marine biological toxin point assessment method, and evaluated the edible safety by comparing with the acute reference dose (ARfD) recommended by the European Food Safety Agency. The proportion of samples contaminated with toxins followed a descending order of gyrodimine (GYM) (83.02%), okadaic acid (OA) (51.16%), dinophysis toxin (DTX-2)(40.91%) and homo-yessotoxin (Homo-YTX) (8.6%), with the peak values of 105.4, 31.39, 38.19 and 159.6 µg·kg−1, respectively. The concentrations of toxins varied with different seasons and species. Among the six kinds of shellfish, the highest contents of OA and GYM were found in oyster samples, and the peak values were found in winter and autumn samples, respectively. The contents of DTX-2 in autumn and winter samples were significantly higher than those in spring and summer samples. Homo-YTX was only detected in scallop and mussel samples in summer and autumn. Dietary survey shows that the average daily shellfish consumption of coastal residents in Beibu Gulf was 45 g·d−1. The ESI value in OA toxin was estimated to be 0.26 μg·(kg·d)−1 based on the shellfish consumption and body masses of local residents, less than the ARfD value [0.3 μg·(kg·d)−1] set by the European Union (EU). However, based on the shellfish consumption and body masses of local residents provided by the EU, the ESI value [0.34 μg·(kg·d)−1] was higher than the ARfD value. The results suggest that although the detection rate of lipophilic shellfish toxin in shellfish samples was lower than the safe limit value set by the EU, there is a risk of dietary exposure to OA toxin for the Beibu Gulf residents.
  • 大型海藻主要分布在近岸海域的潮间带和潮下带区域,是海洋生态系统初级生产力的巨大贡献者,也是重要的海水养殖对象[1]。生活在潮间带的大型海藻经历着巨大的环境变化,如光照、温度、水分以及盐度变化等,这些环境因子共同影响着大型海藻的生长、生理和生态位等[1]。研究表明,温度胁迫通过影响植物光合机构的结构和功能,抑制光系统II (PSII) 供体侧水光解反应和PSII下游碳同化反应,进而影响光合电子传递,降低光合速率[2];在高光条件下,植物则会因吸收过量的光而产生光抑制现象,对光合系统造成损伤,导致光合能力下降、生长发育受阻[3];此外,高光胁迫也会影响叶绿体的形态结构和光合色素的组分[3]。由此可见,温度和光照对植物光合生理的影响尤为突出。

    大型海藻属于低等植物,光合生理是所有光合藻类积累有机物、增加生物量的关键生理反应,同时也是大型海藻重要的碳汇过程[4]。有研究表明,大型海藻经过长期的进化已衍生出一套抵御温度和光照胁迫的高效策略,通过降低叶绿素a等光合色素的含量以减少光能吸收是其应对高光胁迫的关键方式之一[3,5]:缘管浒苔 (Ulva linza) 为了避免高光胁迫,其叶绿素a和叶绿素b含量随光照强度的增加而显著降低[5];而通过热耗散机制应对过剩的产能则是大型海藻应对温度和光照胁迫的另一重要方式[3,6-7];海带 (Saccharina japonica) 孢子体在高温条件下,暗适应后的藻体PSII最大光化学效率 (Fv/Fm) 显著下降,但它可以通过非光化学淬灭 (NPQ) 的增加使其热耗散 (光保护) 能力增强,从而使所能利用的光合有效辐射范围更广[7]。然而,尽管目前关于温度和光照对大型海藻光合生理影响的研究已有报道,但不同种属大型海藻的适应范围不同,如:针叶蕨藻 (Caulerpa sertularioides) 在26 ℃时光合活性最好[8],而海带孢子体在环境温度达到14 ℃以上时,光合作用原初光反应过程会受抑制[7];皱紫菜 (Porphyra crispata) 的适合生长光照强度为4 000~7 500 lx,光照强度过低时其光合作用降低[9],而龙须菜 (Gracilaria lemaneiformis) 在光照强度3 000 lx以下生长较好,大于6 000 lx则生长相对较差[10]。因此,开展大型海藻对温度和光照的光合生理响应研究,有助于更科学地指导大型海藻在海水养殖、海洋生态修复等方面的应用。

    异枝江蓠 (Gracilariopsis bailinae) 属于大型经济红藻,其生长速度快,全年可生长,可用于鲍鱼养殖、琼胶提取[11]及养殖水处理等[12-14],在海水养殖和海洋生态修复上具有较好的应用前景,常见于我国的广东省和海南省[15]。目前,针对异枝江蓠的生长[11]、琼胶提取[14]与经济动物混养[16]等已有研究报道,但关于其光合生理方面的研究鲜有报道。为此,本研究利用叶绿素荧光技术,结合藻体光合色素含量和生长率的变化,探究了异枝江蓠对温度和光照强度的光合生理响应特征,以期丰富异枝江蓠在光合生理方面的研究,为其进一步的开发利用提供基础数据支撑。

    异枝江蓠于2021年7月采集自广东省湛江市雷州市乌石镇。采集后的鲜藻低温运回实验室,用灭菌海水和毛刷除去藻体表面的杂质后暂养于藻类培养室。暂养条件为温度25 ℃、光照强度5 000 lx、盐度30‰、pH 8.1、光周期12 L∶12 D。

    选取暂养后健康且生长状态相近的藻体,用吸水纸吸干表面水分后各称取1 g分别放入含有藻类营养液的不同烧杯 (500 mL) 中备用。根据前期预实验结果,设置了不同的温度梯度 (15、20、25、30和35 ℃) 和光照梯度 (1 000、3 000、6 000和9 000 lx),将选取好的藻体放入对应温度或光照强度的光照培养箱中。藻体培养周期为7 d,每个处理设置3个平行,每2 d更换一次新鲜培养液。7 d后称取藻体鲜质量,并取定量的样品用于藻体光合色素和叶绿素荧光参数的测定。

    用吸水纸吸干处理7 d后藻体表面水分并用分析天平 (FA1104,上海) 称量,根据Yong等[17]提出的公式计算藻体特定生长率 (SGR),公式如下:

    $$ {R}_{\mathrm{S}\mathrm{G}\mathrm{R}}=\frac{\mathrm{l}\mathrm{n}{M}_{7}-\mathrm{l}\mathrm{n}{M}_{0}}{7}\times 100\mathrm{\%} $$ (1)

    式中:RSGR表示特定生长率 (%·d−1);M0表示藻体初始鲜质量 (g);M7表示藻体 7 d后的鲜质量 (g)。

    叶绿素a  (Chl a)、类胡萝卜素 (Car) 使用甲醇法提取,称取处理7 d后的藻体0.1 g放入含有4颗研磨钢珠和2 mL甲醇的研磨管中 ,转入到组织细胞研磨仪 (JFstgrp-24,上海) 中,70 Hz研磨180 s,4 ℃黑暗提取24 h后转入高速离心机 (Heraeus Multifuge 3XR,德国) 5 000 r·min−1离心10 min,使用紫外分光光度计 (UV-2450,日本) 分别测量研磨管上清液在波长480、510、652、665、750 nm处的吸收值,根据以下公式计算Chl a和Car的质量分数[18]

    $$ {W}_{\mathrm{C}\mathrm{h}\mathrm{l}a}=\frac{[16.29\times ({A}_{665}-{A}_{750})-8.54\times ({A}_{652}-{A}_{750}\left)\right]\times V}{1\mathrm{ }000\times m} $$ (2)
    $$ {W}_{\mathrm{C}\mathrm{a}\mathrm{r}}=\frac{\left\{7.6\times \left[\left({A}_{480}-{A}_{750}\right)-1.49\times \left({A}_{510}-{A}_{750}\right)\right]\right\}\times V}{1\mathrm{ }000\times m} $$ (3)

    式中:WChla为叶绿素a质量分数 (mg·g−1);WCar为类胡萝卜素质量分数(mg·g−1);A 为吸光值;V 为色素提取液的总体积 (mL);m为藻体鲜质量 (g)

    藻红蛋白 (PE)、藻蓝蛋白 (PC) 采用0.1 mol·L−1磷酸盐缓冲液 (pH 6.8) 提取,提取方法与上述提取Chl a和Car的方法一致,按以下公式计算质量分数[19]

    $$ {W}_{\mathrm{P}\mathrm{E}}=\frac{\left[\left({A}_{564}-{A}_{592}\right)-\left({A}_{455}-{A}_{592}\right)\times 0.20\right]\times 0.12\times V}{F_{\mathrm{w}} }$$
    $$ {W}_{\mathrm{P}\mathrm{C}}=\frac{\left[\left({A}_{618}-{A}_{645}\right)-\left({A}_{592}-{A}_{645}\right)\times 0.51\right]\times 0.15\times V}{F_{\mathrm{w}}} $$

    式中:WPE为藻红蛋白的质量分数 (mg·g−1);WPC为藻蓝蛋白质量分数 (mg·g−1);A为吸光值;V为定容后的溶液体积 (mL);FW为藻体鲜质量 (g)

    处理7 d后的样品暗适应30 min后,使用便携式脉冲调制荧光仪 (FMS-2,英国) 测定其PSII最大光化学效率 (Fv/Fm)、实际光能转化效率 (ΦPSII)、电子传递速率 (ETR)、光化学淬灭系数 (qP) 和非光化学淬灭系数 (NPQ)。

    所有数据以“平均值±标准差 ($\overline x \pm {s} $)”表示,采用IBM SPSS Statistics 25软件进行Pearson相关性分析、单因素方差 (One-way ANOVA) 分析和Duncan's多重比较检验处理 (P<0.05作为差异的显著性水平)。利用双因素方差 (Two-way ANOVA) 分析温度和光照强度对特定生长率、光合色素含量及叶绿素荧光参数的交互作用。

    图1可知,不同温度条件下,光照强度对异枝江蓠的生长有不同程度的促进作用。除15 ℃外,异枝江蓠的特定生长率随光照强度的增加显著上升 (P<0.05),均在9 000 lx达到最大值 [(0.83±0.18)、(4.60±0.25)、(4.41±0.10)和(4.78±0.40)%·d−1],相对1 000 lx分别增加了76.60%、61.75%、173.91%和256.72%;在1 000和3 000 lx下,随温度的上升,特定生长率先上升后下降,特定生长率最大值均出现在25 ℃;而在6 000和9 000 lx下,特定生长率随温度的升高而上升,25~35 ℃的特定生长率[(1.33~4.78)%·d−1]显著高于15~20 ℃[(0.25~0.83)%·d−1, P<0.05]。相关性分析结果显示,在20~35 ℃和3 000~9 000 lx下,特定生长率分别与光照强度和温度呈正相关 (P<0.05),且随温度上升或光照强度增加相关性越强 (表1)。温度和光照强度对异枝江蓠特定生长率产生了极显著的影响 (P<0.01),且存在交互作用 (表2)。表明在一定范围内升高温度和增加光照强度均能促进异枝江蓠的生长。

    图  1  不同温度和光照强度下的异枝江蓠特定生长率
    注:不同小写字母表示在相同温度条件下不同光照强度处理间差异显著 (P<0.05),不同大写字母表示在相同光照强度条件下不同温度处理间差异显著 (P<0.05);后图同此。
    Figure  1.  Specific growth rate of G. bailinae at different temperatures and light intensities
    Note: Different lowercase letters indicate significant differences between different light intensity treatments under the same temperature condition (P<0.05), while different uppercase letters indicate significant differences between different temperature treatments under the same light intensity condition (P<0.05). The same case in the following figures.
    表  1  温度、光照强度与异枝江蓠生长、光合色素及叶绿素荧光参数的相关性分析
    Table  1.  Correlation analysis between temperature and light intensity and growth, photosynthetic pigments and chlorophyll fluorescence parameters of G. bailinae
    项目
    Item
    光照强度 Light intensity/lx温度 Temperature/℃
    1 0003 0006 0009 0001520253035
    特定生长率 SGR 0.501 0.668** 0.873** 0.889** 0.197 0.766** 0.898** 0.938** 0.958**
    叶绿素 a Chl a 0.242 −0.195 −0.704** −0.727** −0.130 −0.782** −0.910** −0.935** −0.950**
    类胡萝卜素 Car −0.042 0.041 −0.188 −0.215 −0.310 −0.736** −0.910** −0.101 −0.730**
    藻红蛋白 PE 0.978** 0.948** 0.886** 0.842** 0.011 −0.849** −0.910** −0.937** −0.950**
    藻蓝蛋白 PC 0.951** 0.969** 0.992** 0.923** −0.130 −0.882** −0.920** −0.947** −0.950**
    PSII最大光化学效率 Fv/Fm 0.804** 0.745** 0.752** 0.792** −0.710** −0.799** −0.420 −0.532 −0.230
    实际光能转化效率 ΦPSII 0.885** 0.776** 0.878** 0.895** −0.910** −0.853** −0.940** −0.913** −0.790**
    电子传递速率 ETR 0.859** 0.537* 0.784** 0.890** −0.910** −0.902** −0.900** −0.771** −0.600*
    光化学淬灭系数 qP 0.918** 0.866** 0.912** 0.929** −0.920** −0.835** −0.920** −0.755** −0.760**
    非光化学淬灭系数 NPQ 0.007 0.679** 0.952** 0.848** 0.912** 0.936** 0.937** 0.966** 0.902**
    注:** 表示在0.01水平 (双尾) 极显著相关,* 表示在0.05水平 (双尾) 显著相关。 Note: **. Extremely significant correlation at 0.01 level (Two-tailed); *. A significant correlation at 0.05 level (Two-tailed).
    下载: 导出CSV 
    | 显示表格
    表  2  温度和光照强度对异枝江特定生长率的双因素方差分析
    Table  2.  Two-way ANOVA analysis for effects of temperature and light intensity on specific growth rate of G. bailinae
    变异来源
    Source of variation
    自由度
    DF
    F显著性
    Sig.
    温度
    Temperature
    4 802.800 <0.000 1
    光照强度
    Light intensity
    3 233.200 <0.000 1
    温度×光照强度
    Temperature×Light intensity
    12 37.810 <0.000 1
    残差
    Residual
    40
    下载: 导出CSV 
    | 显示表格

    在温度为15 ℃时,异枝江蓠的Chl a、Car、PE和PC质量分数随光照强度的增加均无显著性变化 (P>0.05);而温度介于20~35 ℃时,4种色素质量分数整体上随光照强度的增加而显著下降 (P<0.05),且均在光照强度为1 000 lx时达到最大值、在9 000 lx时 (除30 ℃时的Car) 取得最小值 (图2)。在相同的光照强度下,随着温度的上升,Chl a和Car质量分数的变化不一 (图2-a和2-b),而PE和PC质量分数整体上呈增加趋势 (图2-c和2-d),且不同光照强度下两者均在35 ℃时达到最大值,分别为PE:(0.205±0.009)、(0.185±0.014)、(0.119±0.030) 和 (0.104±0.009) mg·g−1,PC:(0.349±0.015)、(0.336±0.020)、(0.243±0.008) 和 (0.168±0.015) mg·g−1。相关性分析结果显示 (表1),PE和PC质量分数在不同光照强度下与温度呈极显著正相关 (P<0.01),在20~35 ℃时与光照强度呈极显著负相关 (P<0.01);Chl a质量分数在6 000~9 000 lx和20~35 ℃时分别与温度和光照强度呈极显著负相关 (P<0.01);Car质量分数与温度无显著相关性 (P>0.05),在20、25和35 ℃时与光照强度呈极显著负相关 (P<0.01)。温度和光照强度对4种色素均有极显著影响 (P<0.01),除Car外 (P>0.05),温度和光照强度对Chl a、PE和PC含量的影响均存在交互作用 (P<0.01,表3)。

    图  2  不同温度和光照强度下的异枝江蓠色素质量分数
    Figure  2.  Pigment mass fraction of G. bailinae at different temperatures and light intensities
    表  3  温度和光照强度对异枝江蓠光合色素的双因素方差分析
    Table  3.  Two-way ANOVA analysis for effects of temperature and light intensity on photosynthetic pigments of G. bailinae
    光合色素 Photosynthetic pigment变异来源 Source of variation自由度 DFF显著性 Sig.
    叶绿素 a Chl a 温度 Temperature 4 36.650 <0.000 1
    光照强度 Light intensity 3 53.480 <0.000 1
    温度×光照强度 Temperature×Light intensity 12 4.415 0.000 2
    残差 Residual 40
    类胡萝卜素 Car 温度 Temperature 4 23.850 <0.000 1
    光照强度 Light intensity 3 14.830 <0.000 1
    温度×光照强度 Temperature×Light intensity 12 1.490 0.168 4
    残差 Residual 40
    藻红蛋白 PE 温度 Temperature 4 448.300 <0.000 1
    光照强度 Light intensity 3 219.400 <0.000 1
    温度×光照强度 Temperature×Light intensity 12 25.370 <0.000 1
    残差 Residual 40
    藻蓝蛋白 PC 温度 Temperature 4 356.000 <0.000 1
    光照强度 Light intensity 3 131.500 <0.000 1
    温度×光照强度 Temperature×Light intensity 12 14.400 <0.000 1
    残差 Residual 40
    下载: 导出CSV 
    | 显示表格

    PSII最大光化学效率在15和20 ℃时随光照强度的增加显著下降 (P<0.05),25~35 ℃时则无显著性变化 (P>0.05);在不同光照强度条件下,PSII最大光化学效率与温度呈极显著正相关 (P<0.01,表1),且25~35 ℃时PSII最大光化学效率 (0.738~0.803) 显著高于15~20 ℃ (0.513~0.610) (P<0.05,图3-a)。实际光能转化效率在不同温度下随光照强度的增加而下降,而电子传递速率在15~30 ℃时与实际光能转化效率的表现一致,在35 ℃时呈先下降后上升再下降的趋势;在相同光照强度下,实际光能转化效率和电子传递速率均随温度上升整体呈上升趋势,且均在35 ℃时达到最大值 (图3-b和3-c)。相关性分析结果显示,实际光能转化效率和电子传递速率与温度呈正相关,与光照强度呈负相关 (P<0.05,表1)。温度和光照强度对PSII的最大光化学效率、实际光能转化效率和电子传递速率均有极显著影响 (P<0.01),且对实际光能转化效率和电子传递速率有交互作用 (表4)。

    图  3  不同温度和光照强度下的异枝江蓠光化学效率参数
    Figure  3.  Photochemical efficiency parameters of G. bailinae at different temperatures and light intensities
    表  4  温度和光照强度对异枝江蓠 PSII 光化学效率参数的双因素方差分析
    Table  4.  Two-way ANOVA analysis for effects of temperature and light intensity on photochemical efficiency parameters of PSII of G. bailinae
    光化学效率参数
    Photochemical efficiency parameter
    变异来源
    Source of variation
    自由度
    DF
    F显著性
    Sig.
    PSII 最大光化学效率 Fv/Fm 温度 Temperature 4 290.800 <0.000 1
    光照强度 Light intensity 3 10.380 <0.000 1
    温度×光照强度 Temperature×Light intensity 12 1.864 0.069 9
    残差 Residual 40
    实际光能转化效率 ΦPSII 温度 Temperature 4 61.560 <0.000 1
    光照强度 Light intensity 3 106.000 <0.000 1
    温度×光照强度 Temperature×Light intensity 12 4.071 0.000 4
    残差 Residual 40
    电子传递速率 ETR 温度 Temperature 4 53.820 <0.000 1
    光照强度 Light intensity 3 70.460 <0.000 1
    温度×光照强度 Temperature×Light intensity 12 4.001 0.000 4
    残差 Residual 40
    下载: 导出CSV 
    | 显示表格

    图4为不同温度和光照强度下的异枝江蓠光化学淬灭系数和非光化学淬灭系数。结果显示:在不同温度下,光化学淬灭系数随光照强度的增加显著下降 (P<0.05),与光照强度呈极显著负相关 (P<0.01,表1);相反,在不同的光照强度条件下,光化学淬灭系数随温度的上升显著上升 (P<0.05),且与温度呈极显著正相关 (P<0.01,图4-a和表1)。非光化学淬灭系数在不同温度下随光照强度的增加显著上升 (P<0.05);在 6 000和9 000 lx时,非光化学淬灭系数随温度的上升显著上升 (P<0.05),在1 000和3 000 lx时则均表现为先上升后下降 (图4-b)。表明温度和光照强度对异枝江蓠荧光淬灭参数具有显著影响 (P<0.05),且存在交互作用 (表5)。

    图  4  不同温度和光照强度下的异枝江蓠荧光淬灭参数
    Figure  4.  Fluorescence quenching parameters of G. bailinae at different temperatures and light intensities
    表  5  温度和光照强度对异枝江蓠荧光淬灭参数的双因素方差分析
    Table  5.  Two-way ANOVA analysis for effects of temperature and light intensity on fluorescence quenching parameters of G. bailinae
    荧光淬灭参数
    Fluorescence quenching parameter
    变异来源
    Source of variation
    自由度
    DF
    F显著性
    Sig.
    光化学淬灭系数  qP 温度 Temperature 4 86.340 <0.000 1
    光照强度 Light intensity 3 73.660 <0.000 1
    温度×光照强度 Temperature×Light intensity 12 3.471 0.001 5
    残差 Residual 40
    非光化学淬灭系数 NPQ 温度 Temperature 4 204.600 <0.000 1
    光照强度 Light intensity 3 338.400 <0.000 1
    温度×光照强度 Temperature×Light intensity 12 33.670 <0.000 1
    残差 Residual 40
    下载: 导出CSV 
    | 显示表格

    不同大型海藻对温度和光照强度的适应性不同,生长率的变化可直观反映其生长状况的变化[8,20]。本研究中,温度 (15和20 ℃) 和光照强度 (1 000 lx) 较低时,异枝江蓠的生长受到抑制;而在20~35 ℃和3 000~9 000 lx下,异枝江蓠的特定生长率分别与光照强度和温度呈正相关 (P<0.05),并在35 ℃、9 000 lx时达到最大。而其他研究显示,龙须菜在水温19~27 ℃和光照强度3 000 lx以下时生长最佳;皱紫菜在17~20 ℃和4 000~7 500 lx时生长较好,高于23 ℃时其生长受到抑制,光照强度过低时其光合作用也会减弱[9]。由此可见,大型海藻对温度和光照强度的适应存在种属差异,过高或过低的温度和光照强度均会抑制藻体生长。大部分大型海藻生长的适宜温度为15~25 ℃[21],而本研究中,异枝江蓠生长最适温度为35 ℃,这与钟志海等[16]报道的结果 (33 ℃) 相近,表明异枝江蓠是一种喜高温的大型海藻。

    光合作用是绿色植物最基本的生命能量来源[22],叶绿素荧光包含了丰富的光合生理信息,是探究不同环境条件对藻类光合作用能力影响的最佳方法[8,23]。PSII最大光化学效率值的变化可以反映光抑制的程度,值越大其发生光抑制的程度越低[24]。相关研究表明,低温、高光条件可能使藻体PSII能量吸收和耗散平衡功能受损,导致PSII最大光化学效率降低,藻体光合系统吸收过剩光能后引起PSII反应中心关闭,即产生光抑制或光保护现象,以保护PSII免受光损伤[21,25]。本研究中,异枝江蓠PSII最大光化学效率在15~20 ℃时随光照强度的增加而显著下降 (P<0.05),且在15~20 ℃时的值显著低于25~35 ℃。可见,低温以及低温时给予的高光可能导致异枝江蓠PSII反应中心关闭,PSII活性降低,有效光能转化效率下降,具体表现为藻体生长受到一定程度的抑制。然而,在25~35 ℃时,高光并未导致异枝江蓠PSII最大光化学效率的显著下降,与缘管浒苔[5]、皱紫菜[6]等的研究结果不一致,这可能与温度和光照强度的共同作用效应相关,具体原因有待进一步研究。

    实际光能转化效率、电子传递速率和光化学淬灭系数与PSII电子传递活性和效率相关,其值越高说明电子传递活性越大[24]。研究表明,针叶蕨藻的生长率变化与其实际光能转化效率、电子传递速率及光化学淬灭系数的变化表现一致[8];然而,在本研究中,三种叶绿素荧光参数均随光照强度的增加而下降,但异枝江蓠的特定生长率却有相反的表现,说明光照强度增加所导致的异枝江蓠电子传递速率下降可能并未达到最终影响藻体生长的程度,异枝江蓠可能通过其他策略适应高光。目前认为,植物在响应高光时主要通过增加热耗散 (即非光化学淬灭系数) 来调节PSII的活性,避免光抑制,实现光保护[25-27]。大型海藻,如海带孢子体[7]、缘管浒苔[5]、皱紫菜[6]、裙带菜 (Undaria pinnatifida)[6] 等相关研究也表明非光化学淬灭系数的增加在藻体适应环境改变的过程中起重要作用。本研究中,异枝江蓠非光化学淬灭系数随光照强度的增加而上升,说明其可能通过增加热耗散的方式调节PSII的活性实现自我保护,整体表现为藻体特定生长率随光照强度的增加而上升。此外,异枝江蓠实际光能转化效率、电子传递速率和光化学淬灭系数在不同光照强度下均随温度的上升而升高,但低光照 (1 000和3 000 lx) 时,其特定生长率和非光化学淬灭系数却随温度的升高先升后降;这可能是因为在低光照条件下,过高的温度造成藻体吸收的多余能量无法耗散,最终导致藻体生长率下降。值得注意的是,目前有研究认为植物通过热耗散保护光系统的现象与光系统的循环电子流 (CEF) 相关,植物在干旱、强光、弱光、高温和低温等逆境胁迫下可以激发CEF来适应环境胁迫,但其响应机制仍不清楚[28],进一步深入研究可能有助于加深对大型海藻逆境响应机制的了解。

    光合色素是植物光合作用的关键物质[21]。研究表明,光合色素含量的下调也是植物实现光保护的一种重要手段[3]。本研究中,异枝江蓠4种色素含量在高光照时显著下降,说明其通过下调光合色素含量进行光保护;而色素含量减少生长率却上升,则可能与色素对能量利用效率的提高有关[29]。此外,PE和PC在高温时含量较高且与温度呈显著正相关 (P<0.05)。PE和PC既是红藻中重要的两种捕光色素蛋白[30],也是藻体细胞的贮存蛋白,有利于藻类在自然界中的生存竞争[31]。温度升高促进了异枝江蓠PE和PC的合成,说明异枝江蓠在高温条件下通过增加藻胆蛋白含量从而提高捕光效率及蛋白贮存,增强其生存能力,实现了在高温条件下的快速生长;这与刺枝鱼栖苔 (Acanthophora spicifera)[32]、皱紫菜[9]等的研究结果类似。

    综上所述,温度和光照强度对异枝江蓠生长及光合生理具有重大影响。异枝江蓠是一种喜高温的大型海藻,藻胆蛋白含量的增加有助于提高其在高温条件下的生存能力;高光照降低了异枝江蓠的电子传递速率,但其可通过增加热耗散和减少光合色素的合成来实现光保护。鉴于在实际的海水养殖和海洋生态修复中,随着纬度的不同以及季节的更替,温度和光合有效辐射均会有较大的波动,且海水深度的改变也会影响光合有效辐射;因此,异枝江蓠比较适合在我国热带、亚热带地区栽培或用作海洋生态修复的备选藻种。在应用过程中,应着重关注温度和光照强度的变化,优先选择水温25 ℃以上、光照强度9 000 lx左右的水域环境;温度和光照强度过低时,可能会使藻体光合生理受限而导致其养殖产量低、碳汇效果不佳、死亡等。

  • 图  1   贝类样品中脂溶性贝类毒素质量分数

    Figure  1.   Mass fractions of lipophilic shellfish toxins in different shellfish samples

    图  2   不同季节脂溶性贝类毒素质量分数

    注:同一图中标有不同字母表示有显著性差异 (P<0.05)。

    Figure  2.   Mass fractions of lipophilic shellfish toxins in shellfish in four seasons

    Note: Different letters within the same figure indicate significant differences (P<0.05).

    图  3   广西北部湾沿海居民海产品摄入种类分布

    Figure  3.   Distribution of seafood species intook by residents in Beibu Gulf in Gangxi Province

    图  4   广西北部湾沿海居民贝肉日均消费量

    Figure  4.   Mean daily shellfish meat consumption of residents in Beibu Gulf in Guangxi Province

    表  1   9 种脂溶性贝类毒素检测参数

    Table  1   Detection parameters of nine lipophilic shellfish toxins

    毒素
    Toxin
    电离模式
    ESI mode
    母离子
    Precursor ion/(m/z)
    子离子
    Product ion/(m/z)
    碰撞电压
    Collision energy/eV
    大田软海绵酸
    OA
    [M+H] 803.4 255, 563.2 −65, −55
    鳍藻毒素1
    DTX-1
    [M+H] 817.5 255.1, 113.1 −65, −100
    鳍藻毒素2
    DTX-2
    [M+H] 803.5 255.1, 563.2 −65, −55
    虾夷扇贝毒素
    YTX
    [M+H] 1 141.5 1 061.5, 855.5 −65, −50
    Homo-虾夷扇贝毒素
    Homo-YTX
    [M+H] 1 155.5 1 075.5, 869.5 −65, −50
    环亚胺毒素
    GYM
    [M+H]+ 508.3 174.3, 490.3 50, 35
    螺环内酯毒素
    SPX1
    [M+H]+ 692.5 444.4, 674.4 45, 45
    原多甲藻酸1
    AZA1
    [M+H]+ 842.5 806.5, 824.5 40, 50
    原多甲藻酸3
    AZA3
    [M+H]+ 828.5 810.4, 792.5 40, 50
    下载: 导出CSV

    表  2   贝中脂溶性贝类毒素检出率及含量

    Table  2   Detection rates and contents of lipophilic shellfish toxins in shellfish samples

    脂溶性贝类毒素
    Lipophilic shellfish
    toxin
    检出值范围
    Detection range/
    (μg·kg−1)
    平均检出值
    Mean/
    (μg·kg−1)
    检出率
    Detection
    rate/%
    环亚胺毒素
    GYM
    0.12~105.54 33.1 83.02
    螺环内酯毒素
    SPX1
    ND ND ND
    大田软海绵酸
    OA
    0.34~31.39 7.97 51.16
    鳍藻毒素1
    DTX-1
    ND ND ND
    鳍藻毒素2
    DTX-2
    0.29~38.19 7.65 40.91
    原多甲藻酸1
    AZA1
    ND ND ND
    原多甲藻酸3
    AZA3
    ND ND ND
    虾夷扇贝毒素
    YTX
    ND ND ND
    Homo-虾夷扇贝毒素
    Homo-YTX
    3.83~159.66 47.97 8.60
    注:ND表示未检出。 Note: ND. Not detected.
    下载: 导出CSV

    表  3   脂溶性贝类毒素的毒性等效因子与急性参考值

    Table  3   Toxic equivalent factors and ARfD of some shellfish toxins

    毒素组
    Toxin group
    毒素
    Toxin
    毒性等效因子
    Toxion equivalence factor
    急性参考剂量
    ARfD/[μg·(kg·d)−1]
    参考文献
    Reference
    大田软海绵酸
    OA
    大田软海绵酸
    TOA
    1 0.3 [6, 25]
    鳍藻毒素1
    TDTX-1
    1
    鳍藻毒素2
    TDTX-2
    0.6
    虾夷扇贝毒素
    YTXs
    Homo-虾夷扇贝毒素
    THomo-YTX
    1 25 [9, 25]
    45-OH-虾夷扇贝毒素
    T45-OH-YTX
    1
    45-OH-Homo虾夷扇贝毒素
    T45-OH-Homo-YTX
    0.5
    下载: 导出CSV

    表  4   广西北部湾沿海居民脂溶性贝类毒素急性暴露评估

    Table  4   Acute dietary exposure to some shellfish toxins of residents in Beibu Gulf in Guangxi Province

    毒素
    Toxin
    贝肉日P99消费量
    P99 daily shellfish
    consumption/g
    毒素检出的最大值
    Maximum toxin detection
    value/(μg·kg−1)
    每日摄入量
    ESI/[μg·(kg·d)−1]
    急性参考剂量
    ARfD/[μg·(kg·d)−1]
    第一组
    First group
    第二组
    Second group
    大田软海绵酸
    OA
    205 31.39 0.26 0.34 0.3
    38.19
    虾夷扇贝毒素
    YTX
    159.66 0.60 1.06 25
    下载: 导出CSV

    表  5   广西北部湾沿海居民脂溶性贝类毒素慢性暴露评估

    Table  5   Chronic dietary exposure to some shellfish toxins of residents in Beibu Gulf in Guangxi Province

    毒素
    Toxin
    贝肉日均消费
    Daily shellfish
    consumption/g
    毒素检出的均值
    Mean toxin detection
    value/(μg·kg−1)
    贝肉日P99消费量
    P99 daily shellfish
    consumption/g
    每日摄入量
    ESI/[μg·(kg·d)−1]
    第一组
    First group
    第二组
    Second group
    大田软海绵酸
    OA
    45 7.97 205 6.58×10−3 0.03
    鳍藻毒素2
    DTX-2
    7.65 3.79×10−3 0.017
    Homo-虾夷扇贝毒素
    Homo-YTX
    47.97 0.04 0.18
    下载: 导出CSV
  • [1]

    KONOKI K, ONODA T, WATANBE R, et al. In vitro acylation of okadaic acid in the presence of various bivalves' extracts[J]. Mar Drugs, 2013, 11(2): 300-315.

    [2] 王亚军, 余新威, 方力, 等. 营养环境对微小亚历山大藻C4生长和产毒的影响[J]. 水产学报, 2017, 41(10): 1588-1598.
    [3]

    VALDIGLESIAS V, PREGO-FARALDO M V, PÁSARO E, et al. Okadaic acid: more than a diarrheic toxin[J]. Mar Drugs, 2013, 11(11): 4328-4349. doi: 10.3390/md11114328

    [4]

    MARTELLI F, CIRLINI M, DELLAFIORA L, et al. Mitigation of marine toxins by interactions with bacteria: the case of okadaic acid and tetrodotoxin[J]. Food Control, 2021(1): 108428.

    [5]

    AUNE T, ESPENES A, AASEN J, et al. Study of possible combined toxic effects of azaspiracid-1 and okadaic acid in mice via the oral route[J]. Toxicon, 2012, 60(5): 895-906. doi: 10.1016/j.toxicon.2012.06.007

    [6]

    EFSA. Scientific Opinion of the Panel on Contaminants in the Food Chain on a request from the European Commission on marine biotoxins in shellfish-pectenotoxin group[J]. EFSA J, 2009, 1109: 1-47.

    [7]

    BLAANCO J, ARÉVALO F, CORREA J, et al. Effect of the industrial steaming on the toxicity, estimated by LC-MS/MS, of mussels exposed for a long time to diarrhetic shellfish poisoning (DSP) toxins[J]. Food Chem, 2015, 177: 240-247. doi: 10.1016/j.foodchem.2015.01.012

    [8]

    OTERO A, CHAPELA M J, ATANASSOVA M, et al. Cyclic imines: chemistry and mechanism of action: a review[J]. Chem Res Toxicol, 2011, 24(11): 1817-1829. doi: 10.1021/tx200182m

    [9]

    EFSA. Opinion of the Scientific Panel on Contaminants in the Food chain on a request from the European Commission on marine biotoxins in shellfish-yessotoxin group[J]. EFSA J, 2008, 907: 1-62.

    [10]

    EFSA panel on contaminants in the food chain (CONTAM). Scientific opinion on marine biotoxins in shellfish-cyclic imines (spirolides, gymnodimines, pinnatoxins and pteriatoxins)[J]. EFSA J, 2010, 8(6): 1628.

    [11] 徐轶肖, 韦光领, 王玉, 等. 钦州湾海水和香港牡蛎体内脂溶性贝类毒素污染特征[J]. 海洋与湖沼, 2021, 52(1): 144-152. doi: 10.11693/hyhz20200400126
    [12] 柳阳, 郭伟, 黎全江, 等. 广西北部湾近江牡蛎脂溶性藻毒素污染状况与安全风险评价[J]. 南方水产科学, 2020, 16(4): 108-113. doi: 10.12131/20190212
    [13]

    JI Y, YAN G W, WANG G X, et al. Prevalence and distribution of domoic acid and cyclici mines in bivalve mollusks from Beibu Gulf, China[J]. J Hazard Mater, 2022, 423: 127078. doi: 10.1016/j.jhazmat.2021.127078

    [14]

    WANG Z, DOUCETTE G J. Determination of lipophilic marine biotoxins by liquid chromatography-tandem mass spectrometry in five shellfish species from Washington State, USA[J]. J Chromatogr A, 2021, 1639: 461902. doi: 10.1016/j.chroma.2021.461902

    [15]

    ROURKE W A, JUSTASON A, MARTIN J L, et al. Shellfish toxin uptake and depuration in multiple Atlantic Canadian mollusca species: application to selection of sentinel species in monitoring programs[J]. Toxins, 2021, 13(2): 168. doi: 10.3390/toxins13020168

    [16]

    ZHENG R J, LIN S E, YANG Y, et al. Variability and profiles of lipophilic marine toxins in shellfish from southeastern China in 2017−2020[J]. Toxicon, 2021, 201: 37-45. doi: 10.1016/j.toxicon.2021.08.009

    [17]

    CHEN X, HUANG B Q, ZHAO Q H, et al. Shellfish contamination with lipophilic toxins and dietary exposure assessments from consumption of shellfish products in Shenzhen, China[J]. Ecotoxicol Environ Saf, 2021, 221: 112446. doi: 10.1016/j.ecoenv.2021.112446

    [18] 海洋讯. 自然资源部发布2021年度《中国海洋灾害公报》[J]. 自然资源通讯, 2022(10): 39-55.
    [19]

    MARTINS J C, DOMÍNGUEZ-PÉREZ D, AZEVEDO C B, et al. Molecular responses of mussel Mytilus galloprovincialis associated to accumulation and depuration of marine biotoxins okadaic acid and dinophysistoxin-1 revealed by shotgun proteomics[J]. Front Mar Sci, 2020, 7: 589822. doi: 10.3389/fmars.2020.589822

    [20] 刘红河, 秦逍云, 廖仕成, 等. 超高效液相色谱-四极杆串联线性离子阱质谱法测定贝类中的脂溶性贝类毒素[J]. 卫生研究, 2021, 50(6): 967-974.
    [21]

    QIU J B, CHEN H D, JI Y, et al. Evaluation of different strategies to minimize thematrix effects on LC-MS/MS analysis of multiple lipophilic shellfish toxins in both acidic and alkaline chromatographic conditions[J]. Toxicon, 2020, 188: 16-26. doi: 10.1016/j.toxicon.2020.10.002

    [22] 赵丽云, 马冠生, 朴建华, 等. 2010-2012中国居民营养与健康状况监测总体方案[J]. 中华预防医学杂志, 2016, 50(3): 204-207.
    [23]

    TWINER M J, El-ADKI R, KILCOYNE J, et al. Comparative effects of the marine algal toxins azaspiracid-1, -2, and -3 on Jurkat T lymphocyte cells[J]. Chem Res Toxicol, 2012, 25(3): 747-754. doi: 10.1021/tx200553p

    [24]

    ALVES R N, RAMBLA-ALEGRE M, BRAGA A C, et al. Bioaccessibility of lipophilic and hydrophilic marine biotoxins in seafood: an in vitro digestion approach[J]. Food Chem Toxicol, 2019, 129: 153-161. doi: 10.1016/j.fct.2019.04.041

    [25]

    EFSA. Scientific opinion of the panel on contaminants in the food chain on a request from the European commission on marine biotoxins in shellfish-summary on regulated marine biotoxins[J]. EFSA J, 2009, 1306: 1-23.

    [26] 白志毅, 温鹏超, 袁立, 等. 我国淡水贝类种质资源现状与保护利用[J]. 水产学报, 2022, 46(1): 149-157.
    [27] 丁剑楠, 张闪闪, 武旭跃, 等. 太湖贡湖湾水体中微囊藻毒素的时空分布及健康风险评估[J]. 生态环境学报, 2018, 27(11): 2095-2101.
    [28]

    CHEN J H, HAN T Z, LI X T, et al. Occurrence and distribution of marine natural organic pollutants: lipophilic marine algal toxins in the Yellow Sea and the Bohai Sea, China[J]. Sci Total Environ, 2018, 612: 931-939. doi: 10.1016/j.scitotenv.2017.08.304

    [29] 查道军, 李嫒芳, 丁任业, 等. 大亚湾大鹏澳海域贝类和浮游植物中脂溶性毒素及软骨藻酸研究[J]. 海洋环境科学, 2022, 41(5): 753-760. doi: 10.13634/j.cnki.mes.2022.05.016
    [30]

    WU H Y, YAO J H, GUO M M, et al. Distribution of marine lipophilic toxins in shellfish products collected from the Chinese market[J]. Mar Drugs, 2015, 13(7): 4281-4295. doi: 10.3390/md13074281

    [31] 韩蕾, 赵芮, 刘昭, 等. 大田软海绵酸和鳍藻毒素时间分辨荧光免疫层析试纸条的研制与应用[J]. 海洋环境科学, 2022, 41(5): 783-790.
    [32] 郑旭颖, 李兆新, 孙晓杰, 等. 渤海海域唐山贝类养殖区腹泻性和麻痹性贝类毒素的监测与风险评估[J]. 渔业科学进展, 2022, 43: 1-12.
    [33]

    BLANCO J, ARÉVALO F, CORREA J, et al. Lipophilic toxins in Galicia (NW Spain) between 2014 and 2017: incidence on the main molluscan species and analysis of the monitoring efficiency[J]. Toxins, 2019, 11(10): 612. doi: 10.3390/toxins11100612

    [34]

    LEITE I do P, SANDRINI-NETO L, SQUELLA F L, et al. Toxin accumulation, detoxification and oxidative stress in bivalve (Anomalocardia flexuosa) exposed to the dinoflagellate Prorocentrumlima[J]. Aquat Toxicol, 2021, 232: 105738.

    [35]

    LEI L, NING W C, YG C, et al. Seasonal variability of Protoceratium reticulatum and yessotoxins in Japanese scallop Patinopecten yessoensis in northern Yellow Sea of China[J]. Toxicon, 2017, 139: 31-34.

    [36]

    VILARIÑO N, LOUZAO M, ABAL P, et al. Human poisoning from marine toxins: unknowns for optimal consumer protection[J]. Toxins, 2018, 10(8): 324. doi: 10.3390/toxins10080324

    [37]

    EFSA Panel on Contaminants in the Food Chain (CONTAM). Statement on further elaboration of the consumption figure of 400 g shellfish meat on the basis of new consumption data[J]. EFSA J, 2010, 8(8): 1706. doi: 10.2903/j.efsa.2010.1706

    [38]

    MCCARRON P, KILCOYNE J, HESS P. Effects of cooking and heat treatment on concentration and tissue distribution of okadaic acid and dinophysistoxin-2 in mussels (Mytilus edulis)[J]. Toxicon, 2008, 51(6): 1081-1089. doi: 10.1016/j.toxicon.2008.01.009

    [39]

    DÍEZ-QUIJADA JIMÉNEZ L, GUZMÁNGUILLÉN R, CASCAJOSA LIRA A, et al. In vitro assessment of cyanotoxins bioaccessibility in raw and cooked mussels[J]. Food Chem Toxicol, 2020, 140: 111391. doi: 10.1016/j.fct.2020.111391

  • 期刊类型引用(6)

    1. 蒋利进,衣萌萌,可小丽,曹建萌,刘志刚,王章,谭泽加,卢迈新,王淼. 碱度对尼罗罗非鱼肝脏免疫、肠道菌群结构和代谢物的影响. 淡水渔业. 2025(01): 49-59 . 百度学术
    2. 江德泓,陈思宇,徐佳慧,张海洋,汪杨怀煜,李学军,石西. 淇河鲫幼鱼对碳酸盐碱度的耐受性研究. 河南水产. 2024(04): 23-26 . 百度学术
    3. 肖黄璇,刘天奇,鲁翠云,孙志鹏,那荣滨,常玉梅,郑先虎. 急性碳酸氢钠碱度胁迫对草鱼存活及相关酶活性的影响. 广东海洋大学学报. 2024(06): 1-9 . 百度学术
    4. 艾克达·热合曼,李雯,黄晶,梁利群,张立民,常玉梅. NaHCO_3碱度和NaCl盐度对两种规格高体雅罗鱼幼鱼的急性毒理研究. 水产学杂志. 2023(05): 32-41 . 百度学术
    5. 赵贺,赵年桦,李丽,聂竹兰,魏杰,沈建忠. 南疆地区四种土著裂腹鱼类幼鱼对氯化钠和碳酸氢钠的急性耐受能力. 淡水渔业. 2022(02): 34-41 . 百度学术
    6. 逯冠政,么宗利,来琦芳,高鹏程,周凯,朱浩拥,刘一萌,孙真. 高盐碱环境下大口黑鲈幼鱼生长性能、血液生理指标与质构特征研究. 渔业科学进展. 2022(04): 1-11 . 百度学术

    其他类型引用(5)

图(4)  /  表(5)
计量
  • 文章访问数:  437
  • HTML全文浏览量:  133
  • PDF下载量:  45
  • 被引次数: 11
出版历程
  • 收稿日期:  2023-02-16
  • 修回日期:  2023-02-23
  • 录用日期:  2023-03-01
  • 网络出版日期:  2023-03-05
  • 刊出日期:  2023-08-04

目录

/

返回文章
返回