• 中文核心期刊
  • 中国科技核心期刊
  • CSCD来源期刊
  • DOAJ 收录
  • Scopus 收录

液相色谱-高分辨质谱法在中国东北地区农作物有机硒形态分析中的应用

朱帅, 沈亚婷, 贾静, 劳昌玲

朱帅, 沈亚婷, 贾静, 劳昌玲. 液相色谱-高分辨质谱法在中国东北地区农作物有机硒形态分析中的应用[J]. 岩矿测试, 2021, 40(2): 262-272. DOI: 10.15898/j.cnki.11-2131/td.202005130070
引用本文: 朱帅, 沈亚婷, 贾静, 劳昌玲. 液相色谱-高分辨质谱法在中国东北地区农作物有机硒形态分析中的应用[J]. 岩矿测试, 2021, 40(2): 262-272. DOI: 10.15898/j.cnki.11-2131/td.202005130070
ZHU Shuai, SHEN Ya-ting, JIA Jing, LAO Chang-ling. Determination of Organic Selenium Compounds in Crops by Liquid Chromatography-Quadrupole/Electrostatic Field Orbitrap High-resolution Mass Spectrometry[J]. Rock and Mineral Analysis, 2021, 40(2): 262-272. DOI: 10.15898/j.cnki.11-2131/td.202005130070
Citation: ZHU Shuai, SHEN Ya-ting, JIA Jing, LAO Chang-ling. Determination of Organic Selenium Compounds in Crops by Liquid Chromatography-Quadrupole/Electrostatic Field Orbitrap High-resolution Mass Spectrometry[J]. Rock and Mineral Analysis, 2021, 40(2): 262-272. DOI: 10.15898/j.cnki.11-2131/td.202005130070

液相色谱-高分辨质谱法在中国东北地区农作物有机硒形态分析中的应用

基金项目: 

中国地质调查局中国地质科学院基本科研业务费项目 CSJ201905

国家自然科学基金项目(21607033);中国地质调查局中国地质科学院基本科研业务费项目(CSJ201905)

国家自然科学基金项目 21607033

详细信息
    作者简介:

    朱帅, 硕士, 助理研究员, 主要从事有机污染物与环境行为研究。E-mail: zhu15131215153@126.com

  • 中图分类号: O657.63

Determination of Organic Selenium Compounds in Crops by Liquid Chromatography-Quadrupole/Electrostatic Field Orbitrap High-resolution Mass Spectrometry

  • 摘要: 随着公众对富硒产品的关注,对富硒农作物中硒含量的研究越来越受到人们的重视。针对现有方法测定农作物有机硒的基体效应和多原子离子质谱干扰问题,本研究利用液相色谱与抗干扰能力强的四极杆/静电场轨道阱高分辨质谱相结合,建立了农产品中有机硒的形态分析方法。经优化以超声波辅助蛋白酶XIV提取农作物中有机硒,用Waters Symmetry RP18液相色谱柱分离,以纯水和乙腈作为流动相进行梯度洗脱,采用正、负离子切换同时测定了硒代蛋氨酸、硒代胱氨酸、甲基硒代半胱氨酸。三种有机硒在负离子模式下采用Full MS/dd-MS2扫描模式对目标物进行定性筛查和确证,三种有机硒在5~500μg/L的浓度范围内线性关系良好,相关系数均大于0.997,检出限为1.0~4.7μg/kg,回收率为76.4%~116.0%,相对标准偏差(RSD)≤6.5%。与现有方法比较,本方法线性范围宽,适用于多种谷物和蔬菜类样品的分析。将建立的方法应用于测定中国粮食主产区黑龙江地区多种农作物的可食部分,结果发现谷类中以硒代蛋氨酸为主,蔬菜类以甲基硒代半胱氨酸为主,并且该地区农作物的有机硒含量处于较低水平,亟待强化补硒。
    要点

    (1) 利用液相色谱-四级杆/静电场轨道阱高分辨质谱法建立了农产品中硒形态高效、高灵敏的方法。

    (2) 该方法与现有方法有机硒的检出限和回收率相当,但本方法线性范围更宽,应用范围更广,可用于多种蔬菜谷物等农作物的分析,解决了现有方法只可应用于某种或某类样品分析的难题。

    (3) 黑龙江地区农作物可食部分中硒形态的组成,谷类以硒代蛋氨酸为主,蔬菜类以甲基硒代半胱氨酸为主,并且该调查区的农作物硒含量处于缺硒水平。

    HIGHLIGHTS

    (1) An efficient and high-sensitivity method for determination of selenium forms in agricultural products by liquid chromatography-quadrupole/electrostatic field orbitrap high-resolution mass spectrometry was established.

    (2) This method had the same detection limit and recovery rate compared to the existing determination method for organic selenium. However, the method established in this study had a wider linear range and application range. This method can be used for the analysis of a variety of vegetables, grains and other crops, and solves the problem that the existing methods can only measure a specific type of sample.

    (3) In crops of Heilongjiang area, the selenium species in cereals were mainly selenomethionine, and in vegetables were mainly methylselenocysteine. The selenium content of crops in the survey area was at a relatively low level.

  • 硒(Se)是人类必需的营养元素,也是植物的有益元素,它以多种天然形式存在于人类、动植物甚至蓝细菌中[1]。Rayman[2]于2012年指出硒是有效的免疫刺激剂,对机体免疫非常重要,硒还具有降低氧化应激反应、延缓衰老,以及促进细胞分化等作用[3]。当人们硒摄入量不足时会导致多种疾病的发生[4],例如当糖尿病人每天Se的摄入量不足20μg时会增加患癌症、阿兹海默氏症或帕金森氏症等其他疾病的风险[5]。在食物链中,植物硒的代谢优于动物硒,植物是补充硒的重要来源[6]。中国存在富硒土地[7],可用于种植富硒作物,食用富硒农作物是预防硒缺乏症最有效的途径[8]。因此,对农产品中硒研究具有现实意义。硒的代谢途径、生物利用度、生理功能和毒性不仅取决于摄入量,更取决于元素的化学形式[9]。通常,硒在农作物中主要以有机硒的形式存在(占比40%~90%)。例如有机硒的形态有硒代蛋氨酸(SeMet)、硒代胱氨酸(SeCys2)、甲基硒代半胱氨酸(MeSeCys)、硒代乙硫氨酸(SeEt)、硒脲(SeUr)、亚硒基全蛋氨酸等[10]。与无机硒相比,有机硒具有更高的生物利用度,SeMet的生物利用度高达90%以上,而人体摄入过量的无机硒会引起中毒[11-12]。人们普遍在土壤或者叶面施用含硒肥料用于生产富硒农产品,主要通过转运蛋白或磷酸盐转运蛋白的转运作用,外源无机硒被农作物根吸收,在芽和树叶中转化为有机硒,然后在叶绿体作用下进一步代谢合成硒代半胱氨酸(SeCys)、硒代蛋氨酸(SeMet)、硒代胱氨酸(SeCys2)等有机硒[13]。因此,在评估富硒农作物时,主要考虑这三种有机硒的含量与组成。

    分离分析农产品中硒形态主要包括提取和分析。提取过程中需特别注意不能使农产品中硒形态发生变化,近些年酶提取法在硒的形态分析研究广泛应用[14],超声微波结合酶提取法可以有效地缩短提取时间[15]。硒在农产品中含量较低,因此对分析方法的灵敏度和准确度要求比较高。硒化物的分离技术主要有气相色谱法(GC)[16]、电泳法[17]和高效液相色谱法(HPLC)等。气相色谱具有分离速度快、检出限低等优点,但分析时间长,且不适用分离低聚肽的硒化合物。电泳法具有检出限低、分辨率高、保留样品的完整性等优点,但是接口技术的滞后阻碍了CE-ICP-MS在分析农作物硒形态上的应用。高效液相色谱技术广泛应用于分离植物体中各种硒形态[18-19]。硒形态的检测器主要有电感耦合等离子体质谱(ICP-MS)[19-20]、原子吸收光谱(AAS)[21-22]、原子荧光光谱(AFS)[23-25]等。ICP-MS具有灵敏度高、检出限低、选择范围宽等优点,与阴离子交换色谱联合使用在农作物硒化物形态分析方面有广泛的应用,但氩的多原子离子会产生质谱干扰(如40Ar和36Ar+76Se,40Ar和37Ar+77Se,40Ar和38Ar+78Se等),基质干扰严重,现有的分析方法应用范围窄,不能同时适用于多种复杂基质样品的测定分析。同时Se属于难电离物质,电离不完全,影响测定的准确度[26],此外,ICP-MS仅适用于拥有标准品的目标分析物的鉴定。因此,质谱干扰问题、电离不完全、不能同时测定多种复杂基质样品仍是分析工作面临的挑战。

    静电场轨道离子阱(Orbitrap)具有较强的抗干扰能力、高精度和高通量等优点,可在很大程度上解决质谱干扰的问题[27]。四极杆/静电场轨道阱高分辨质谱(Q-Exactive)是将四极杆与Orbitrap相结合,降低了复杂基质背景的干扰,测定准确,适用于复杂基质样品的测定分析。该质谱拥有高质量精度和高分辨率,从而能够识别未知物的结构信息[28],在有机硒分析方面初见成效。Bierla等[29]利用电喷雾Orbitrap MS3对葵花籽油中复杂的含Se化合物进行高质量精确度(<10-6)分子特异性检测。本文将超声辅助酶提法与液相色谱-四极杆/静电场轨道阱高分辨质谱相结合,初次尝试同时测定中国粮食主产区东北黑龙江地区多种农产品中三种有机硒形态。对样品提取和仪器分析各参数进行优化,建立一种前处理过程快速简单、检出限低、准确度高和适合多种复杂基质农产品中硒形态分析方法。

    农作物样品(包括水稻13份、玉米13份、大豆9份、马铃薯2份、白菜2份、尖椒2份、白萝卜2份、茄3份)于2019年9月采自黑龙江地区(图 1)。该地区主要为黑土,pH值为7.21~8.24,有机质含量18.9~26.2mg/g,氮含量95.8~137.4mg/kg,磷含量233.7~321.0mg/kg,钾含量183.1~211.7mg/kg,含水率3.2%~7.5%。样品在野外先用自来水冲洗去除表面灰尘和残留土壤,再用矿泉水润洗后用纱布包裹晾样。待样品运回实验室后离开取出,在室温下自然风干。风干的样品放在托盘中于55℃烘24h,将烘干的样品用万能粉碎机粉碎后储存于干燥器中,用于测定总硒和硒形态。

    图  1  采样点位图
    Figure  1.  Locations of sampling sites

    液相色谱-四极杆/静电场轨道阱高分辨质谱系统(LC-Q Exactive Orbitrap MS,美国ThermoFisher公司);电子天平(德国Sartorius公司);离心机(日本Hitachi公司);B5510-E超声波清洗仪(美国Branson公司)。

    硒代胱氨酸和硒代半胱氨酸(SeCys2和MeSeCys,纯度98%,北京百灵威公司);硒代蛋氨酸(SeMet,纯度99%,美国Alfa公司);脂肪酶和淀粉酶(分析纯,北京索莱宝科技有限公司);蛋白酶K和蛋白酶XIV[天根生化科技(北京)有限公司];甲醇和乙腈(色谱纯,美国J.T.Baker公司);甲酸(色谱纯,德国Fluka公司);Milli-Q水:由Millipore纯水仪制备。

    总硒含量测定:用四位天平精确称取约0.5g样品置于消解罐中,同时做试剂空白。总硒含量参照国家标准《食品安全国家标准食品中硒的测定》(GB 5009.93-2017),用ICP-MS进行测定,以72Ge作为内标,同位素78Se和80Se进行测定。

    有机硒形态含量测定:用四位天平精确称取0.5g样品置于离心管中(15mL),做试剂空白,加入5mL超纯水和20mg蛋白酶XIV,涡旋30s混匀,在37℃水浴中超声1.5h,超声后移入离心机离心分离,离心机转速6000r/min,时间10min,静置后取上清液过0.22μm滤膜,上机测定。Q-Exactive气质的主要仪器参数见表 1

    表  1  梯度洗脱程序
    Table  1.  Gradient elution procedure
    时间(min) 流速(mL/min) 流动相A (%) 流动相B (%)
    初始 0.300 97.0 3.0
    8.0 0.300 85.0 15.0
    9.0 0.300 50.0 50.0
    10.0 0.300 10.0 90.0
    12.0 0.300 10.0 90.0
    12.5 0.300 97.0 3.0
    16.5 0.300 97.0 3.0
    下载: 导出CSV 
    | 显示表格

    高效液相色谱条件:Waters Symmetry RP18液相色谱柱(3.9mm×150mm, 5μm);进样量3μL;柱温20℃;流动相:A水,B乙腈,流速0.3mL/min,梯度洗脱程序如表 1所示。质谱条件:加热电喷雾离子源(HESI);喷雾电压3.0kV;毛细管温度350℃;加热器温度300℃;扫描模式为Full MS/dd-MS2;质谱扫描范围m/z 50~500;一级全扫描(Full MS)分辨率为70000,C-trap最大容量(AGC target)为3×106,C-trap最大注入时间为200ms;数据依赖二级子离子扫描(dd-MS2)分辨率为17500,C-trap最大容量(AGC target)为1×105,C-trap最大注入时间为50ms,二级质谱的碰撞能量(CE)为30eV;鞘气压力35arb;辅助气压力10arb。

    样品提取、净化、分析过程中严格按照1.3节分析条件对样品进行质量控制。在样品分析时,每7个样品设置一个流程空白作为实验室背景干扰。实验过程中,有机硒在流程空白样品未检出,三种有机硒的方法检出限(S/N=3)为1.0~4.7μg/kg,回收率为76.4%~116.0%,相对标准偏差(RSD)≤6.5%,通过标准品测定证实实验过程不会造成有机硒的形态发生改变。

    液相色谱柱的选择对分离至关重要,为获得更好的灵敏度和分离效果,本研究比较了Hamilton PRP-X100阴离子交换柱(4.1mm×250mm, 10μm)、Waters Symmetry RP18液相色谱柱(3.9mm×150mm, 5μm)、Agilent ZORBAX Eclipse Plus C18(RRHD, 2.1mm×150mm, 18μm)、Waters BEH C18液相色谱柱(2.1mm×50mm, 1.7μm)来分离有机硒。如图 2图 3a所示,三种有机硒在Waters Symmetry RP18液相色谱柱上的峰形最好,响应值最高,具有更好的保留特性,该色谱柱特有的极性基团嵌入技术使有机硒在该色谱柱上的峰形比其余的三个色谱柱的峰形更优。此外, 该液相色谱柱可用于分离多种亲水化合物,对实际样品的检测适用性更好、更耐用,因此本实验选择Waters Symmetry RP18液相色谱柱。

    图  2  有机硒在不同液相色谱柱上的色谱图
    a-Hamilton PRP-X100阴离子交换柱(4.1mm×250mm, 10μm);b-Waters Symmetry RP18液相色谱柱(3.9mm×150mm, 5μm);c-Agilent ZORBAX Eclipse Plus C18液相色谱柱(RRHD,2.1mm×150mm, 18μm);d-Waters BEH C18液相色谱柱(2.1mm×50mm,1.7μm)。
    Figure  2.  Chromatograms of organic selenium on different liquid chromatography columns
    图  3  有机硒在不同仪器参数下的响应值
    a-有机硒在不同液相色谱柱上的响应值;b-不同流动相对硒形态的影响;c-不同温度对硒形态的影响。
    Figure  3.  Values of organic selenium under different instrument parameters

    依据优先选择简单和无毒试剂作为流动相的原则,本实验选择A水、B乙腈和A水、B甲醇作为流动相进行优化分析。如图 3b所示,有机相采用乙腈的信号强度优于甲醇相,当流动相水中加入0.1%的甲酸或1%的乙酸铵后,有机硒的信号强度骤减,故采用纯水和乙腈作为流动相。如图 3c所示,柱温在20℃时,三种有机硒的响应值最高,优化后的色谱条件见1.3节所述,各有机硒组分在10min内获得较好的分离效果,目标物峰形对称性好、尖锐度和灵敏度高。

    本方法对三种目标物在负离子模式下采用Full MS/dd-MS2扫描模式对目标物进行定性筛查和确证。对于复杂基质中痕量组分的定性,分辨率起到了非常关键的作用,本研究发现一级全扫描分辨率R=70000时,所有待测物可以与基质中的干扰物实现基线分离,响应值也明显提高,基质干扰大大降低。因此,考虑到定性和定量的准确性,一级质谱全扫描的分辨率设置为70000,二级质谱扫描设置为17500。进一步对二级质谱的碰撞能量(CE)进行优化。结果显示,CE值为30eV时能够获得最佳的二级质谱定量离子响应。三种有机硒的保留时间、精确质量数、相对质量偏差等信息见表 2。由表中数据可知,三种有机硒精确质量数的相对质量偏差均低于3×10-6,证明目标化合物的质量准确度均很好,满足定性需求。

    表  2  三种有机硒的保留时间、精确质量数和相对质量偏差
    Table  2.  Retention times, accurate masses and relative deviations of three kinds of organic selenium
    硒形态 保留时间(min) 精确质量数理论值(m/z) 精确质量数测定值(m/z) 相对质量偏差(10-6)
    SeMet 5.64 195.98822 195.98862 2.04
    SeCyS2 3.41 334.90547 334.90593 1.37
    MeSeCys 3.41 167.95692 167.95669 1.37
    下载: 导出CSV 
    | 显示表格

    选取硒含量较高的水稻籽粒样品进行提取优化实验,采用37℃水浴超声1.5h后,保持37℃水浴过夜的提取方法分别考察水、10%甲酸、脂肪酶、α-淀粉酶、蛋白酶K、蛋白酶E和蛋白酶XIV对水稻籽粒样品中硒形态提取效果的影响。如图 4a所示,水稻籽粒样品中主要的硒形态为硒代蛋氨酸,采用蛋白酶E和蛋白酶XIV的提取效果最佳,三种形态硒的响应值总和最高;水、10%甲酸、脂肪酶和α-淀粉酶无法将SeCyS2从水稻籽粒中提取出来;蛋白酶K的提取效率仅为蛋白酶E或蛋白酶XIV的50%。蛋白酶E或蛋白酶XIV的提取效率相当,但是蛋白酶XIV对SeCys2的提取效率比蛋白酶E的提取效率高(图 4b)。本实验选择蛋白酶XIV作为提取剂考察其对硒形态提取效果的影响。

    图  4  有机硒在不同前处理条件下的响应值
    a、b-不同提取剂对有机硒提取效率的影响;
    c-蛋白酶XIV用量对有机硒提取效率的影响;
    d-提取时间对有机硒提取效率的影响。
    Figure  4.  Values of organic selenium under different instrument parameters

    水浴超声处理的温度为37℃,提取时间设为3h,分别考察10、20、30、40、50和60mg蛋白酶XIV对水稻籽粒样品中硒形态提取效果的影响。如图 4c所示,当蛋白酶XIV用量为20mg时,提取效率最佳;当蛋白酶XIV用量大于20mg时,提取效率不再增加。随着蛋白酶XIV用量的增加,提取出的SeMet含量先升高后降低,当蛋白酶XIV用量为20mg时,提取出的SeMet含量最高。SeCys2含量随蛋白酶XIV用量的增加而逐渐增加。MeSeCys含量随蛋白酶XIV用量增加先升高后下降。最终选择蛋白酶XIV用量为20mg。

    在水浴超声处理的温度为37℃,蛋白酶XIV用量为20mg的条件下,分别考察提取时间为1.0、1.5、2.0、2.5和3.0h对水稻籽粒样品中硒形态提取效果的影响。如图 4d所示,提取效率随提取时间增长而增加,当提取时间≥1.5h时,提取效率基本保持一致;提取出的SeMet和SeCys2含量随提取时间的增加而逐渐减少,这是由于这两种有机硒随时间的增长逐渐酶解导致其含量降低,MeSeCys含量随时间增长变化不大,选择提取时间为1.5h最优。

    分别配制5.0、10.0、50.0、100.0、200.0、500.0 μg/L的SeCys2、MeSeCys和SeMet混合标准溶液系列,以各浓度色谱峰面积(y)对应质量浓度(x)绘制标准曲线。在最佳实验条件下,在5.0~500.0μg/L浓度范围内,相关系数(r)均小于0.997,线性关系良好。采用逐级稀释法,三倍的信噪比(S/N=3)确定该化合物的检出限,SeCys2、MeSeCys和SeMet方法检出限为1.0~4.7μg/kg(表 3)。

    表  3  三种硒的线性方程、相关系数和检出限
    Table  3.  Regression equations, correlation coefficients and detection limit of three selenium species
    硒形态 线性方程 相关系数(r) 检出限(μg/kg)
    水稻籽粒 黄豆 玉米 白菜 马铃薯 白萝卜
    SeMet y=302208+ 45923.1x 0.9997 1.0 1.7 1.2 1.3 1.8 1.4
    SeCys2 y=-2987460+ 163256x 0.9982 2.5 3.4 2.3 2.6 3.5 2.7
    MeSeCys y=1231460+ 63033.2x 0.9970 3.3 4.7 3.5 3.1 3.7 3.3
    下载: 导出CSV 
    | 显示表格

    分别选取一空白水稻籽粒样品,添加适量的有机硒标准溶液至浓度为50、1000和5000μg/kg三个水平的三种硒形态混合标准溶液,每个浓度水平做7个平行样,进行回收率实验。表 4中的结果表明,SeMet的加标回收率在89.2%~106.0%,SeCys2的加标回收率在81.2%~107.1%,MeSeCys的加标回收率在87.2%~116.0%。在所选取的农作物样品中三种硒形态的相对标准偏差(RSD)均≤6.5%(表 4),方法的重复性良好。

    表  4  水稻籽粒中三种有机硒的加标回收率和精密度
    Table  4.  Standard recovery rate and precision tests of three kinds of organic selenium in rice grain
    硒形态 添加量(μg/kg) 测量值(μg/kg) 加标回收率(%) RSD (%)
    SeMet 50. 0 42.7~51.3 89.2~102.6 5.1
    1000.0 937.2~1060.1 93.7~106.0 3.7
    5000.0 4930.0~5110.6 98.6~102.2 1.3
    SeCys2 50.0 40.6~49.0 81.2~98.0 6.5
    1000.0 998.3~1070.8 99.8~107.1 3.1
    5000.0 4820.3~4960.5 96.4~99.2 0.9
    MeSeCys 50.0 43.6~51.6 87.2~103.2 6.2
    1000.0 1020.7~1160.2 102.1~116.0 4.7
    5000.0 5070.6~5190.4 101.4~103.8 1.3
    下载: 导出CSV 
    | 显示表格

    选取有机硒未检出的黄豆、玉米、白菜、马铃薯和白萝卜参考水稻籽粒样品进行回收率实验,黄豆、玉米、白菜、马铃薯和白萝卜中三种有机硒的加标回收率分别为76.4%~93.8%、86.9%~108.2%、81.0%~93.8%、79.9%~88.7%和80.9%~92.3%, 表明该方法适用于不同类型的样品。

    表 5将本研究建立的测定有机硒形态的方法与文献中报道的HPLC-ICP-MS、HPLC-HG-AFS等方法进行比较。可以看出,本方法测定的化合物检出限和回收率可与文献报道的方法相媲美。传统的测定硒形态的方法只适用于某种或某类样品的分析,本方法可降低基质干扰,线性范围更宽,应用范围更广,可用于多种蔬菜谷物等农作物的分析。因此,LC-Q Exactive Orbitrap MS可作为HPLC-ICP-MS的替代仪器来分析农作物的有机硒形态。

    表  5  本方法与文献中相关分析方法的比较
    Table  5.  Comparison of performance with other methods in the references
    硒形态 样品 样品提取方法 分析方法 线性范围(μg/L) 回收率(%) 检出限 来源文献
    SeMet, SeCys2, MeSeCys, Se(Ⅳ), Se(Ⅵ) 谷类食品 超声酶提取 HPLC-ICP-MS 0.5~200.0 60.5~120.8 2.5~10.0μg/kg [30]
    SeMet, SeCys2, MeSeCys, Se(Ⅳ), Se(Ⅵ) 水稻籽粒 超声酶提取 HPLC-HG-AFS 10.0~200.0 85.4~105.5 2.27~3.89μg/L [31]
    SeMet 富硒面包 微波酶提取 HPLC-IDA-ICP-MS 20.0~100.0 68.0~84.0 0.7μg/L [32]
    SeMet, SeCys2, MeSeCys, Se(Ⅳ), Se(Ⅵ) 白菜 超声提取 HPLC-ICP-MS 0.0~200.0 52.1~111.6 0.85~4.68mg/mL [33]
    SeMet, SeCys2, MeSeCys 蔬菜和谷物 超声酶提取 LC-Q Exactive Orbitrap MS 5.0~500.0 76.4~116.0 1.0~4.7μg/kg 本文
    下载: 导出CSV 
    | 显示表格

    采用优化好的前处理及检测方法对采集的农作物样品中硒形态进行分析。检出率为30.43%,表 6列出了部分检出农作物可食用部分有机硒含量。SeCys2、MeSeCys和SeMet在黄豆、玉米样品中均有不同程度的检出。水稻籽粒样品中除MeSeCys外,其余两种有机硒化合物均有检出。SeCys2在白菜、马铃薯和白萝卜这三种蔬菜中均未检出。在黄豆、玉米和水稻籽粒谷物类样品中SeMet占总有机硒含量的53.1%~100.0%,占总硒含量的64.1%~80.0%,说明谷类食品中的硒以硒代蛋氨酸为主,此结果与文献中报道的大米和小麦中硒以硒代蛋氨酸为主一致[14, 34]。与谷物类样品不同的是,蔬菜类样品中的硒以甲基硒代半胱氨酸为主,占总有机硒含量的58.8%~100.0%,占总硒含量的35.7%~37.4%。蔬菜中有机硒的含量与谷物产品相比较低,这可能是因为蔬菜中的硒主要通过土壤、根施硒肥或喷施叶面硒肥的方式获得的,在生长周期不足的情况下,没有得到很好的生物转化,大部分硒还是以无机硒的形式存在。部分检出有机硒的谷物和蔬菜样品对硒的转化率从强到弱依次为:玉米>黄豆>水稻籽粒>白菜>马铃薯>白萝卜。所检出谷物和蔬菜样品玉米对有机硒转化能力最强,硒转化率达72.73%。由此可见,谷物尤其是玉米可作为人体补充硒的重要来源。

    表  6  部分检出农作物可食用部分有机硒含量
    Table  6.  Organic selenium content in edible parts of some plants
    样品 有机硒含量(μg/kg) 无机硒转化率(%) 有机硒检出率(%)
    SeMet SeCys2 MeSeCys 总有机硒 总硒
    水稻籽粒-1 92.0 12.7 ND 104.7 162.5 64.4 30.8
    水稻籽粒-2 84.0 10.3 ND 94.4 135.2 69.8
    水稻籽粒-3 11.7 ND ND 11.7 19.2 60.9
    水稻籽粒-4 109.3 27.8 17.3 154.4 195.6 78.9
    黄豆-1 43.4 27.3 11.1 81.8 112.7 72.6 44.4
    黄豆-2 32.5 ND ND 32.5 45.3 71.7
    黄豆-3 87.2 34.8 24.2 146.2 200.4 73.0
    黄豆-4 45.6 12.5 ND 58.1 85.7 67.8
    玉米-1 84.7 21.3 8.6 114.6 143.2 80.0 23.1
    玉米-2 47.0 9.8 ND 57.3 89.4 64.1
    玉米-3 58.0 34.7 ND 93.1 125.7 74.1
    白菜 ND ND 8.0 8.0 21.4 37.4 50.0
    马铃薯 5.6 ND 10.3 15.9 43.7 36.4 50.0
    白萝卜 6.3 ND 9.0 15.3 42.8 35.7 50.0
    注:ND表示未检出。
    下载: 导出CSV 
    | 显示表格

    由于采集样品时蔬菜未完全成熟,本研究仅将谷物作物的硒含量与国内外研究作对比。本研究测定黑龙江水稻样品有机硒含量为0.0~154.4μg/kg,整体上低于文献报道湖北恩施富硒地区富硒水稻的有机硒含量80.0~1540.0μg/kg[35],与陕西省安康市大米有机硒含量24.0~373.0μg/kg处于同一水平[30]。黑龙江玉米样品有机硒含量为0.0~143.2μg/kg,远远低于陕西省安康市玉米有机硒含量3668.0~10738.0μg/kg[30]。黑龙江大豆样品硒含量为0.0~200.4μg/kg,与西班牙东南部豆类硒含量17.9~268.7μg/kg相当[36],低于江苏省苏州大豆硒含量1.5~123.6mg/kg[37]。通过对比发现调查区的谷物硒含量处于缺硒水平,植物体内的硒含量与形态种类除了受植物本身对硒的迁移转化机制和土壤硒含量的影响外,还受成土母质、土壤类型、土壤pH值和有机质含量等因素的影响。采集农作物样品的研究区域以黑土为主,还有少量的黑钙土,该区域土壤呈碱性,pH值为7.21~8.24,土壤有机质含量(18.9~26.2mg/g)比较高,研究区域农作物硒含量较低,但该地区所采集的土壤硒总量(0.29~1.16mg/kg)整体处于足硒甚至富硒水平。因此,有必要开展东北地区农作物的富硒机制研究,为农作物科学种植、合理补硒、改善及预防缺硒性地方病提供科学的理论依据。

    本研究建立了高效液相色谱-四极杆/静电场轨道阱高分辨质谱法测定多种农作物中有机硒形态的方法。该方法的线性范围广,可广泛应用于富硒和非富硒多种农作物有机硒的测定。此外,该方法的适用范围广,可用于多种蔬菜谷物等农作物的分析,解决了现有方法只可测定某种或某类样品的难题。

    本研究对采集于东北黑龙江地区的水稻籽粒、玉米、大豆、马铃薯、白菜、尖椒和白萝卜等46份样品中的硒形态进行分析研究,发现谷类以硒代蛋氨酸为主,蔬菜类以甲基硒代半胱氨酸为主。研究区域的农作物中硒含量处于较低水平,有必要开展农作物的富硒机制研究,为提升该地区人民健康水平提供科学支撑。

  • 图  1   采样点位图

    Figure  1.   Locations of sampling sites

    图  2   有机硒在不同液相色谱柱上的色谱图

    a-Hamilton PRP-X100阴离子交换柱(4.1mm×250mm, 10μm);b-Waters Symmetry RP18液相色谱柱(3.9mm×150mm, 5μm);c-Agilent ZORBAX Eclipse Plus C18液相色谱柱(RRHD,2.1mm×150mm, 18μm);d-Waters BEH C18液相色谱柱(2.1mm×50mm,1.7μm)。

    Figure  2.   Chromatograms of organic selenium on different liquid chromatography columns

    图  3   有机硒在不同仪器参数下的响应值

    a-有机硒在不同液相色谱柱上的响应值;b-不同流动相对硒形态的影响;c-不同温度对硒形态的影响。

    Figure  3.   Values of organic selenium under different instrument parameters

    图  4   有机硒在不同前处理条件下的响应值

    a、b-不同提取剂对有机硒提取效率的影响;
    c-蛋白酶XIV用量对有机硒提取效率的影响;
    d-提取时间对有机硒提取效率的影响。

    Figure  4.   Values of organic selenium under different instrument parameters

    表  1   梯度洗脱程序

    Table  1   Gradient elution procedure

    时间(min) 流速(mL/min) 流动相A (%) 流动相B (%)
    初始 0.300 97.0 3.0
    8.0 0.300 85.0 15.0
    9.0 0.300 50.0 50.0
    10.0 0.300 10.0 90.0
    12.0 0.300 10.0 90.0
    12.5 0.300 97.0 3.0
    16.5 0.300 97.0 3.0
    下载: 导出CSV

    表  2   三种有机硒的保留时间、精确质量数和相对质量偏差

    Table  2   Retention times, accurate masses and relative deviations of three kinds of organic selenium

    硒形态 保留时间(min) 精确质量数理论值(m/z) 精确质量数测定值(m/z) 相对质量偏差(10-6)
    SeMet 5.64 195.98822 195.98862 2.04
    SeCyS2 3.41 334.90547 334.90593 1.37
    MeSeCys 3.41 167.95692 167.95669 1.37
    下载: 导出CSV

    表  3   三种硒的线性方程、相关系数和检出限

    Table  3   Regression equations, correlation coefficients and detection limit of three selenium species

    硒形态 线性方程 相关系数(r) 检出限(μg/kg)
    水稻籽粒 黄豆 玉米 白菜 马铃薯 白萝卜
    SeMet y=302208+ 45923.1x 0.9997 1.0 1.7 1.2 1.3 1.8 1.4
    SeCys2 y=-2987460+ 163256x 0.9982 2.5 3.4 2.3 2.6 3.5 2.7
    MeSeCys y=1231460+ 63033.2x 0.9970 3.3 4.7 3.5 3.1 3.7 3.3
    下载: 导出CSV

    表  4   水稻籽粒中三种有机硒的加标回收率和精密度

    Table  4   Standard recovery rate and precision tests of three kinds of organic selenium in rice grain

    硒形态 添加量(μg/kg) 测量值(μg/kg) 加标回收率(%) RSD (%)
    SeMet 50. 0 42.7~51.3 89.2~102.6 5.1
    1000.0 937.2~1060.1 93.7~106.0 3.7
    5000.0 4930.0~5110.6 98.6~102.2 1.3
    SeCys2 50.0 40.6~49.0 81.2~98.0 6.5
    1000.0 998.3~1070.8 99.8~107.1 3.1
    5000.0 4820.3~4960.5 96.4~99.2 0.9
    MeSeCys 50.0 43.6~51.6 87.2~103.2 6.2
    1000.0 1020.7~1160.2 102.1~116.0 4.7
    5000.0 5070.6~5190.4 101.4~103.8 1.3
    下载: 导出CSV

    表  5   本方法与文献中相关分析方法的比较

    Table  5   Comparison of performance with other methods in the references

    硒形态 样品 样品提取方法 分析方法 线性范围(μg/L) 回收率(%) 检出限 来源文献
    SeMet, SeCys2, MeSeCys, Se(Ⅳ), Se(Ⅵ) 谷类食品 超声酶提取 HPLC-ICP-MS 0.5~200.0 60.5~120.8 2.5~10.0μg/kg [30]
    SeMet, SeCys2, MeSeCys, Se(Ⅳ), Se(Ⅵ) 水稻籽粒 超声酶提取 HPLC-HG-AFS 10.0~200.0 85.4~105.5 2.27~3.89μg/L [31]
    SeMet 富硒面包 微波酶提取 HPLC-IDA-ICP-MS 20.0~100.0 68.0~84.0 0.7μg/L [32]
    SeMet, SeCys2, MeSeCys, Se(Ⅳ), Se(Ⅵ) 白菜 超声提取 HPLC-ICP-MS 0.0~200.0 52.1~111.6 0.85~4.68mg/mL [33]
    SeMet, SeCys2, MeSeCys 蔬菜和谷物 超声酶提取 LC-Q Exactive Orbitrap MS 5.0~500.0 76.4~116.0 1.0~4.7μg/kg 本文
    下载: 导出CSV

    表  6   部分检出农作物可食用部分有机硒含量

    Table  6   Organic selenium content in edible parts of some plants

    样品 有机硒含量(μg/kg) 无机硒转化率(%) 有机硒检出率(%)
    SeMet SeCys2 MeSeCys 总有机硒 总硒
    水稻籽粒-1 92.0 12.7 ND 104.7 162.5 64.4 30.8
    水稻籽粒-2 84.0 10.3 ND 94.4 135.2 69.8
    水稻籽粒-3 11.7 ND ND 11.7 19.2 60.9
    水稻籽粒-4 109.3 27.8 17.3 154.4 195.6 78.9
    黄豆-1 43.4 27.3 11.1 81.8 112.7 72.6 44.4
    黄豆-2 32.5 ND ND 32.5 45.3 71.7
    黄豆-3 87.2 34.8 24.2 146.2 200.4 73.0
    黄豆-4 45.6 12.5 ND 58.1 85.7 67.8
    玉米-1 84.7 21.3 8.6 114.6 143.2 80.0 23.1
    玉米-2 47.0 9.8 ND 57.3 89.4 64.1
    玉米-3 58.0 34.7 ND 93.1 125.7 74.1
    白菜 ND ND 8.0 8.0 21.4 37.4 50.0
    马铃薯 5.6 ND 10.3 15.9 43.7 36.4 50.0
    白萝卜 6.3 ND 9.0 15.3 42.8 35.7 50.0
    注:ND表示未检出。
    下载: 导出CSV
  • Germ M, Stibilj V. Selenium and plants[J]. Acta Agriculturae Slovenica, 2007, 89(1): 65-71.

    Rayman M P. Selenium and human health[J]. The Lancet, 2012, 379: 1256-1268. doi: 10.1016/S0140-6736(11)61452-9

    Wichman J, Winther K H, Bonnema S J, et al. Selenium supplementation significantly reduces thyroid auto-antibody levels in patients with chronic autoimmune thyroiditis: A systematic review and meta-analysis[J]. Thyroid, 2016, 26(12): 1681-1692. doi: 10.1089/thy.2016.0256

    Wei J, Zeng C, Gong Q Y, et al. The association between dietary selenium intake and diabetes: Across-sectional study among middle-aged and older adults[J]. Nutrition Journal, 2015, 14(1): 18. doi: 10.1186/s12937-015-0007-2

    Labunskyy V M, Hatfield D L, Gladyshev V N. Selenoproteins: molecular pathways and physiological roles[J]. Physiological Reviews, 2014, 94(3): 739-777. doi: 10.1152/physrev.00039.2013

    Achibat H, Alomari N A, Messina F, et al. Organoselenium compounds as phytochemicals from the natural kingdom[J]. Natural Product Communications, 2015, 10(11): 1885-1892. http://smartsearch.nstl.gov.cn/paper_detail.html?id=19e6514ec141eb022dfdd3c135175b29

    Toan D Q, Cui Zewei C, Huang J H, et al. Selenium distribution in the Chinese environment and its relationship with human health: A review[J]. Environment International, 2018, 112: 294-309. doi: 10.1016/j.envint.2017.12.035

    Turło J, Gutkowska B, Herold F, et al. Biological availability and preliminary selenium speciation in selenium-enriched mycelium of Lentinula edodes (Berk. )[J]. Food Biotechnology, 2011, 25(1): 16-29. doi: 10.1080/08905436.2011.547113

    Liu K, Zhao Y, Chen F, et al. Purification and identification of Se-containing antioxidative peptides from enzymatic hydrolysates of Se-enriched brown rice protein[J]. Food Chemistry, 2015, 187: 424-430. doi: 10.1016/j.foodchem.2015.04.086

    Liu L, Yun Z, He B, et al. Efficient interface for online coupling of capillary eectrophoresis with inductively coupled plasma-mass spectrometry and its application in simultaneous speciation analysis of arsenic and selenium[J]. Analytical Chemistry, 2014, 86(16): 8167-8175. doi: 10.1021/ac501347d

    Foster L, Sumar S. Selenium in health and disease: A review[J]. Critical Reviews in Food Science & Nutrition, 1997, 37(3): 211-228.

    Hadrup N, Ravn-Haren G. Acute human toxicity and mortality after selenium ingestion: A review[J]. Journal of Trace Elements in Medicine and Biology, 2020, 58: 126435. doi: 10.1016/j.jtemb.2019.126435

    Zhou F, Dinh Q T, Yang W, et al. Assessment of spe-ciation and in vitro bioaccessibility of selenium in Se-enriched Pleurotus ostreatus and potential health risks[J]. Ecotoxicology and Environmental Safety, 2019, 185: 109675. doi: 10.1016/j.ecoenv.2019.109675

    Guzmán Mar J L, Reyes L H, Mizanur Rahman G M, et al. Simultaneous extraction of arsenic and selenium species from rice products by microwave-assisted enzymatic extraction and analysis by ion chromatography-inductively coupled plasma-mass spectrometry[J]. Journal of Agricultural and Food Chemistry, 2009, 57(8): 3005-3013. doi: 10.1021/jf803598k

    Vale G, Pereira S, Mota A, et al. Enzymatic probe sonication as a tool for solid-liquid extraction for total selenium determination by electrothermal-atomic absorption spectrometry[J]. Talanta, 2007, 74(2): 198-205. doi: 10.1016/j.talanta.2007.05.051

    Zhang X, Yang L, Mester Z. Determination of amino acids in selenium-enriched yeast by gas chromatography-mass spectrometry after microwave assisted hydrolysis[J]. Analytica Chimica Acta, 2012, 744: 54-59. doi: 10.1016/j.aca.2012.07.011

    Kannamkumarath S S, Wrobel K, Wuilloud R G. Study-ing the distribution pattern of selenium in nut proteins with information obtained from SEC-UV-ICP-MS and CE-ICP-MS[J]. Talanta, 2005, 66(1): 153-159. doi: 10.1016/j.talanta.2004.10.010

    Khanam A, Platel K. Bioavailability and bioactivity of selenium from wheat (Triticum aestivum), maize (Zea mays), and pearl millet (Pennisetum glaucum), in selenium-deficient rats[J]. Journal of Agricultural and Food Chemistry, 2019, 67(22): 6366-6376. doi: 10.1021/acs.jafc.9b02614

    do Nascimento da Silva E, Aureli F, D'Amato M, et al. Selenium bioaccessibility and speciation in selenium-enriched lettuce: Investigation of the selenocompounds liberated after in vitro simulated human digestion using two-dimensional HPLC-ICP-MS[J]. Journal of Agricultural and Food Chemistry, 2017, 65(14): 3031-3038. doi: 10.1021/acs.jafc.7b01188

    Hildebrand J, Greiner A, Drexler H, et al. Determination of eleven small selenium species in human urine by chromatographic-coupled ICP-MS methods[J]. Journal of Trace Elements in Medicine and Biology, 2020, 61: 126519. doi: 10.1016/j.jtemb.2020.126519

    Heinrichs H, Keltsch H. Determination of arsenic, bismuth, cadmium, selenium, and thallium by atomic absorption spectroscopy with a volatilization technique[J]. Analytical Chemistry, 1982, 54(7): 1211-1214. doi: 10.1021/ac00244a049

    Mohammed E, Mohammed T, Mohammed A. Optimization of instrument conditions for the analysis for mercury, arsenic, antimony and selenium by atomic absorption spectroscopy[J]. Methods, 2018, 5: 824-833. http://www.ncbi.nlm.nih.gov/pubmed/30112290

    Ma Y, Zeng K, Duan T. Fast and clean determination of total selenium in biological materials by an improved oxygen flask combustion method combined with hydride generation atomic fluorescence spectrometry[J]. Microchemical Journal, 2019, 148: 743-747. doi: 10.1016/j.microc.2019.04.064

    Tsopelas F N, Ochsenkühn-Petropoulou M T, Mergias I G, et al. Comparison of ultra-violet and inductively coupled plasma-atomic emission spectrometry for the on-line quantification of selenium species after their separation by reversed-phase liquid chromatography[J]. Analytica Chimica Acta, 2005, 539(1): 327-333. http://www.sciencedirect.com/science/article/pii/S0003267005003697

    Huang C, Li Q, Mo J, et al. Ultratrace determination of tin, germanium, and selenium by hydride generation coupled with a novel solution-cathode glow discharge-atomic emission spectrometry method[J]. Analytical Chemistry, 2016, 88(23): 11559-11567. doi: 10.1021/acs.analchem.6b02807

    Klencsár B, Li S, Balcaen L, et al. High-performance liquid chromatography coupled to inductively coupled plasma-mass spectrometry (HPLC-ICP-MS) for quantitative metabolite profiling of non-metal drugs[J]. TrAC Trends in Analytical Chemistry, 2018, 104: 118-134. doi: 10.1016/j.trac.2017.09.020

    D'Amato R, Regni L, Falcinelli B, et al. Current knowledge on selenium biofortification to improve the nutraceutical profile of food: A comprehensive review[J]. Journal of Agricultural and Food Chemistry, 2020, 68(14): 4075-4097. doi: 10.1021/acs.jafc.0c00172

    韩梅, 侯雪, 邱世婷, 等. 超高效液相色谱-四极杆/静电场轨道阱高分辨质谱测定蔬菜中的百草枯[J]. 分析测试学报, 2019, 38(9): 1126-1131. doi: 10.3969/j.issn.1004-4957.2019.09.016

    Han M, Hou X, Qiu S T, et al. Determination of paraquat in vegetables by ultra performance liquid chromatography-quadrupole/electrostatic field orbitrap high resolution mass spectrometry[J]. Journal of Instrumental Analysis, 2019, 38(9): 1126-1131. doi: 10.3969/j.issn.1004-4957.2019.09.016

    Bierla K, Flis-Borsuk A, Suchocki P, et al. Speciation of selenium in selenium-enriched sunflower oil by high-performance liquid chromatography-inductively coupled plasma mass spectrometry/electrospray-orbitrap tandem mass spectrometry[J]. Journal of Agricultural and Food Chemistry, 2016, 64(24): 4975-4981. doi: 10.1021/acs.jafc.6b01297

    陈绍占, 唐德剑, 李晓玉, 等. 谷类食品中硒形态超声酶提取-高效液相色谱-电感耦合等离子体质谱法测定[J]. 中国公共卫生, 2020, 36(1): 130-136. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGW202001033.htm

    Chen S Z, Tang D J, Li X Y, et al. Determination of selenium species in cereal food with ultrasonic enzyme extraction and high performance liquid chromatography-inductivel coupled plasma mass spectrometry[J]. China Journal Public Health, 2020, 36(1): 130-136. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGW202001033.htm

    胡文彬, 贾彦博, 魏琴芳, 等. 应用液相色谱-原子荧光联用仪测定富硒大米中的5种硒形态[J]. 分析仪器, 2019(1): 120-124. doi: 10.3969/j.issn.1001-232x.2019.01.024

    Hu W B, Jia Y B, Wei Q F, et al. Determination of selenium form in selenium-rich rice by liquid chromatography-hydride generation atomic fluorescence spectrometry[J]. Analytical Instrumentation, 2019(1): 120-124. doi: 10.3969/j.issn.1001-232x.2019.01.024

    Sánchez-Martínez M, Pérez-Corona T, Cámara C, et al. Preparation and characterization of a laboratory scale selenomethionine-enriched bread. Selenium bioaccessibility[J]. Journal of Agricultural and Food Chemistry, 2015, 63(1): 120-127. doi: 10.1021/jf505069d

    Thosaikham W, Jitmanee K, Sittipout R, et al. Evaluation of selenium species in selenium-enriched pakchoi (Brassica chinensis Jusl var parachinensis (Bailey) Tsen & Lee) using mixed ion-pair reversed phase HPLC-ICP-MS[J]. Food Chemistry, 2014, 145: 736-742. doi: 10.1016/j.foodchem.2013.08.116

    Cubadda F, Aureli F, Ciardullo S, et al. Changes in sele-nium speciation associated with increasing tissue concentrations of selenium in wheat grain[J]. Journal of Agricultural and Food Chemistry, 2010, 58(4): 2295-2301. doi: 10.1021/jf903004a

    邵鹏威, 路国慧, 郑宇, 等. 高效液相色谱-电感耦合等离子体质谱测定大米粉中的硒形态[J]. 环境化学, 2020, 39(5): 1434-1441. https://www.cnki.com.cn/Article/CJFDTOTAL-HJHX202005031.htm

    Shao P W, Lu G H, Zheng Y, et al. Determination of selenium species in rice flour using high performance liquid chromatography-inductively coupled plasma mass spectrometry[J]. Environmental Chemistry, 2020, 39(5): 1434-1441. https://www.cnki.com.cn/Article/CJFDTOTAL-HJHX202005031.htm

    Díaz-Alarcón J P, Navarro-Alarcón M, Serrana H G, et al. Determination of selenium in cereals, legumes and dry fruits from southeastern Spain for calculation of daily dietary intake[J]. Science of the Total Environment, 1996, 184(3): 183-189. doi: 10.1016/0048-9697(96)05079-6

    Lu X, He Z, Lin Z, et al. Effects of Chinese cooking methods on the content and speciation of selenium in selenium bio-fortified cereals and soybeans[J]. Nutrients, 2018, 10(3): 1-10. http://www.ncbi.nlm.nih.gov/pubmed/29518925

  • 期刊类型引用(5)

    1. 任昀,廖美林,陈明,黄波,陈杨,杨斌. HPLC-ICP-MS法测定富硒大米中5种硒形态组分. 食品工业. 2025(01): 303-306 . 百度学术
    2. 刘晓杰,杨霄,郭玉斐,张驰,张思远. 富硒农产品中有机硒形态分析研究进展. 现代食品. 2024(13): 75-77+87 . 百度学术
    3. 李乾玉,姚晓慧,刘丽萍,陈绍占,刘洋,何洪巨. 高效液相色谱-电感耦合等离子体质谱法分析研究西兰花中硒形态. 岩矿测试. 2023(03): 523-535 . 本站查看
    4. 张川,李双龙,李春霖,刘薇祎,刘林. 不同形态外源硒对铁皮石斛硒含量的影响研究. 中国林副特产. 2023(05): 67-69 . 百度学术
    5. 宋延斌,王喜宽,夏炎,侯进凯,赵瑞. 河南洛阳市土壤和农作物中硒分布及富集特征. 岩矿测试. 2022(04): 652-662 . 本站查看

    其他类型引用(1)

图(4)  /  表(6)
计量
  • 文章访问数:  4080
  • HTML全文浏览量:  1026
  • PDF下载量:  54
  • 被引次数: 6
出版历程
  • 收稿日期:  2020-05-12
  • 修回日期:  2020-09-28
  • 录用日期:  2021-01-24
  • 发布日期:  2021-03-27

目录

/

返回文章
返回