枸杞是我国传统的药食同源食品[1],宁夏枸杞是唯一载入《2020年版中国药典》的品种,枸杞果实中富含微量元素、多糖、黄酮、甜菜碱、类胡萝卜素等多种功能营养成分[2-3],具有增强免疫力、抗肿瘤、抗疲劳、抗衰老、调节血脂、降血压[4-7]等多重功效,在我国保健行业和医药行业发挥着重要作用。枸杞作为宁夏地区特色优势产业,已成为我区生态治理、劳动力就业、脱贫富民的支柱产业,其产量、消费量和出口量均居世界首位。目前,枸杞中功能成分的研究主要集中在甜菜碱、多糖等有机成分,或者不同产区枸杞中功能营养成分差异性研究[8-9]。近年来,有很多有关中药材无机元素影响其质量的研究[10-12],这些研究表明,无机元素不仅直接或间接影响中药材的生长发育[13],还影响药材中次生代谢产物的积累[12],从而影响药材的药效作用,因矿物元素能够反映农产品生长的土壤环境、元素的种类和含量差异,具有一定地理标志特征,因此,矿物元素分析技术是被公认的植源性食品产地判别中最有效、最有前景的方法之一[14-17],已被广泛应用。但目前溯源指标稳定性及机理探讨的相关研究较为缺乏,不稳定的溯源指标不但增加了溯源的工作量,也可能影响稳定溯源模型的建立,因此,本研究以中宁枸杞为研究对象,检测分析不同成熟期枸杞中55种矿物元素含量,阐明不同成熟期枸杞矿物元素含量的变化特征,研究结果可为枸杞产地溯源提供数据支撑。
混合标准溶液(包括As、Ba、B、Bi、Be、Cd、Co、Cr、Cs、Cu、Fe、Li、Mg、Mn、Eu、Ni、Pb、Rb、Se、Sr、Tl、V、Zn、Ge、Mo、Nb、Ca、Ce、Dy、Gd、Ho、Nd、Pr、Sm、Tb、Lu、Tm、Y、Yb、Er、La、Sc、Sb、Sn、Ti、Zr、Pd、Ag、Ta、Hf、W、Ir、Pt、Ga、Hg元素,共计55种元素),美国 Perkin Elmer公司;硝酸为优级纯,德国 Merck 公司。
标准物质GBW10047(生物成分分析标准物质——胡萝卜),中国地质科学院地球物理地球化学勘查研究所研制。
ELAN DRC-e ICP-MS仪,美国 Perkin Elmer公司;Mars6 Xpress微波消解仪,美国 CEM 公司;AL104电子天平,梅特勒-托利多公司。
1.3.1 样品采集
枸杞样品采自宁夏中宁同3个基地(37°19′10″~37°28′19″ N,105°34′25″~105°55′8″ E)二茬果成熟时期的5个成熟阶段(7 d幼果期、14 d青果期、21 d黄变前期、28 d黄变后期和35 d红熟期)果实。鲜样晒干,置于50 ℃烘箱中烘干,然后进行充分研磨,过100目筛后备用。
1.3.2 样品前处理
枸杞中元素测定方法参考文献[18]。采用微波消解法消解样品,具体操作步骤:称取0.500 0 g(精确到0.000 1 g) 枸杞样品置于微波消解管中,加入硝酸10 mL,常温静置3 h后,放入微波消解仪中进行样品消解,选择温度控制,5 min爬升至120 ℃,保持10 min;5 min爬升至150 ℃,保持20 min;5 min爬升至185 ℃,保持30 min,消解完毕后冷却,轻轻拧开盖子,将微波消解管置于赶酸仪上120 ℃赶酸2 h,然后冷却至室温,用超纯水洗至50 mL刻度试管中,并用超纯水定容,摇匀;同时做试剂空白。
1.3.3 矿物元素含量测定
采用ICP-MS测定枸杞中元素含量。具体的工作条件:发生器功率为1 300 W;雾化器流量为0.95 L/min;等离子炬冷却气流量为17.0 L/min;辅助器流量为1.20 L/min;检测器模拟阶电压为-2 350 V;离子透镜电压为6.00 V。仪器测定枸杞样品及胡萝卜标准物质中55种元素,标准物质测定值在标准值范围内,仪器元素的检出限和定量限见表1。
表1 ICP-MS仪器测定矿物元素标准曲线、检出限和定量限
Table 1 ICP-MS instrument determination of mineral elements standard curve, detection limits and limits of quantification
元素工作方程相关系数检出限/(μg·kg-1)定量限/(μg·kg-1)元素工作方程相关系数检出限/(μg·kg-1)定量限/(μg·kg-1)AsY=4.209 3x0.999 8990.250.83PrY=18.756 2x0.999 9970.000 510.001 7BY=3.986 5x0.999 6391.555.2RbY=56.771 4x0.999 8020.000 0970.000 32BaY=32.081x0.999 9800.005 60.019SbY=4.434 08x0.999 9960.0350.12BeY=3.955 56x0.999 4820.391.3ScY=24.658x0.999 9080.0190.063BiY=30.170 9x0.999 9870.000 390.001 3SeY=0.987 727x0.999 9345.6919.0CdY=13.523 8x0.999 9230.003 30.011SmY=7.149 82x0.999 9950.008 60.029CeY=23.650 2x0.999 9940.000 430.001 4SnY=6.344 48x0.999 9890.0240.080CoY=12.295 7x0.999 9750.0110.037SrY=110.688x0.999 5760.000 210.000 70CrY=13.804 2x0.999 9540.0320.11TbY=71.122 7x0.999 9670.000 0280.000 093CsY=33.609 6x0.999 9860.000 110.000 37TiY=11.888 8x0.999 9530.0990.33CuY=11.755 4x0.999 8980.0140.047ThY=52.208 3x0.999 9851.274.23DyY=16.815 1x0.999 9770.001 00.003 3TlY=90.182 1x0.999 8420.000 0740.000 25ErY=27.982x0.999 9540.000 510.001 7TmY=136.841x0.999 8830.000 006 80.000 023EuY=35.385 6x0.999 970.000 150.000 50VY=17.996 9x0.999 9470.0850.283GdY=12.423 4x0.999 9810.001 70.005 7YY=97.176 9x0.999 7630.000 0270.000 090GeY=7.887 32x0.999 9130.0390.13YbY=46.533 1x0.999 8760.000 0760.000 25HgY=134.757x0.993 4260.000 0340.000 11FeY=2.397 89x0.998 462107.6359HoY=140.913x0.999 8730.000 006 30.000 021ZnY=6.887 32x0.999 6250.421.40LaY=16.892 8x0.999 9970.005 00.017TaY=134.574x0.999 8610.000 0410.000 14 LY=9.353 04x0.999 7990.0210.070IrY=38.915 5x0.999 9440.000 130.000 43 LuY=155.069x0.999 8500.000 004 10.000 014PtY=10.265 1x0.999 9680.005 70.019 MgY=11.382 8x0.999 8840.290.97WY=22.252 1x0.999 9500.000 790.002 6 MnY=14.900 2x0.999 9770.0170.057HfY=45.878 8x0.999 8980.000 140.000 47 MoY=21.107 2x0.999 7200.004 20.014PdY=7.593x0.999 9670.0220.073 NbY=70.750 8x0.999 7910.000 0560.000 19GaY=18.820 3x0.999 9290.0070.023 NdY=18.718 8x0.999 9650.001 30.004 3ZrY=32.594 9x0.999 8730.001 20.004 0 NiY=4.879 06x0.999 9080.301.0CaY=24.748x0.999 773103.7348.8PbY=17.102 2x0.999 9900.001 40.004 7
使用SPSS 25.0 Duncan进行显著性方差分析(P<0.05,差异显著);采用Origin 8.5进行绘图分析;采用SIMCA-P 进行主成分分析。
2.1.1 枸杞成熟过程中矿物元素含量分析
为了方便讨论,将55种元素按照其在枸杞生长过程中的作用分为微量元素、稀土元素、有害金属元素和其他元素,枸杞在5个生长阶段矿物元素含量值见表2。
枸杞果实中微量元素的含量与枸杞的药用价值密切相关,很多微量元素有害或无害是相对的,有益必需的微量元素,若改变其存在状态或浓度,就会产生毒害作用[19]。例如,V是植物生长非必需的微量元素,低浓度V会促进植物生长发育,但高浓度V对植物产生毒性作用[20]。对枸杞5个成熟阶段样品中的11种微量元素进行检测分析,并计算Fe/Mn值。由表2可知,各成熟期的干果中Cu、Co、B、Mg、Mn、Mo、Fe、Zn、Ca、Se、V元素含量范围分别为12.5~16.6、0.11~0.23、9.5~14.2、1 483~3 422、16.1~30.0、0.32~0.64、49.4~86.5、22.5~57.0、1 020~3 273、0.012~0.031、0.050~0.10 mg/kg。除Se元素外,其余10种微量元素在枸杞不同成熟阶段均具有显著性差异(P<0.05)。枸杞作为一种药食同源食品,是一种常用的中药材,中药的功效不是由单一成分的多少决定,而是由多种成分按照合理的比例共同作用的结果,中药中微量元素的生理作用也不能从单一元素的绝对含量来分析,而应该结合不同元素的比例来看[21],牛艳等[22]研究认为,适当的Fe/Mn值可能是宁夏中宁枸杞品质较好的重要原因之一,有研究表明,中药药性从温热到寒凉与Fe/Mn值呈正相关,热药的Fe/Mn值在0.7左右,温药为3.0左右[23]。根据本文中中宁枸杞果实中Fe/Mn值为2.35~3.70,红熟期的Fe/Mn值最高,可以确定中宁枸杞属于温药。
稀土元素是元素周期表中 15 种镧系元素、Sc和Y的总称,在植物体内相对比较稳定。稀土元素以不同形态广泛分布于土壤中,通常主要以难溶残渣态形式存在,但也存在可被植物吸收利用的水溶态、可交换态和有机态[24]。植源性农产品中稀土元素来源与地域环境具有较强的关联性,近年来,稀土元素已被广泛应用于农产品产地溯源研究中。本文对不同成熟期中宁干果枸杞中的16种稀土元素含量进行了分析(表2),可见各成熟阶段的枸杞中16种稀土元素含量均较低,其中La、Ce、Nd轻稀土元素及Sc、Y元素含量在枸杞中的分布较其他稀土元素较为丰富。La元素含量最高,为0.063~0.18 mg/kg,其次是Ce、Nd、Sc、Y元素,含量分别为0.035~0.095、0.014~0.035、0.015~0.033、0.013~0.020 mg/kg,之后是Pr、Eu、Yb、Dy、Er、Gd、Sm元素,含量为10-9数量级,Tb、Ho、Tm、Lu含量较低,为10-10数量级,枸杞中稀土元素显示出轻稀土元素富集、重稀土元素亏欠的特征,这与齐国亮等[25]的研究结果相同。利用ANOVA单因素方差分析比较了不同成熟阶段枸杞中稀土元素含量差异性,结果表明,除了Lu元素外,其余15种稀土元素在5个成熟阶段均存在显著性差异(P<0.05)。
对中宁枸杞中Pb、Cd、As、Hg、Cr、Ni、Ba、Tl、Sb、Sn 10种有害重金属元素含量分析结果表明,各成熟阶段的干果中Cr、Ni、Ba元素含量较高,为mg/kg级;Pb含量次之,含量为0.083~0.28 mg/kg,含量最低的元素为Hg、Tl含量较低,为μg/kg级,比Cd、As、Sb、Sn元素含量低10倍左右。不同成熟阶段枸杞中有害金属元素含量单因素方差分析结果表明,除了Hg、Tl元素,其他8种元素含量在不同成熟阶段枸杞中差异显著(P<0.05)。目前,中药材重金属超标现象依然严重,《中国药典》中规定了重金属残留限量:铅、镉、砷、汞、铜分别为5.0、1.0、2.0、0.2、20 mg/kg,中宁枸杞在整个成熟阶段铅、镉、砷、汞均未超过限量标准,Cu最高含量为16.6 mg/kg,也低于限量值,中宁枸杞质量属安全级。我国现行的有关枸杞重金属元素限量的相关标准,仅对As、Cd、Pb、Hg、Cu的限量值进行了规定,并未对其他元素含量做明确要求,但Ni、Cr、Tl等重金属的摄入对人体危害严重,如镍元素摄入过量会损害人的肝脏和心肺功能[26]。因枸杞为药食同源食品,日摄入量较大,需加强枸杞中重金属的监测。
中宁枸杞中Bi、Be、Sr、Li、Rb、Cs、Ti、Zr、Nb、Pd、Ag、Hf、Ta、W、Ir、Pt、Ga、Ge 18种元素含量分析结果显示,Sr、Rb元素含量较高,分别为7.2~20.6 mg/kg和11.1~15.2 mg/kg,Ti、Li含量分别为2.54~7.69 mg/kg和1.57~1.94 mg/kg,Ga、Cs、Zr、Ag元素含量较上述元素低100倍,Nb、Pt、Ir、W、Ta、Bi、Be、Pd、Ge、Hf含量较低,为μg/kg级。18种元素在5个时期含量的方差分析结果显示,5个阶段枸杞中Be、Nb、Pd、Ta、Hf、Pt、Ge元素含量不存在显著差异(P>0.05),其他11种元素均存在显著差异(P<0.05)。
根据以上分析结果,中宁地区5个成熟阶段枸杞中Cu、Co、B、Mg、Mn、Mo、Fe、Zn、Ca、V、La、Ce、Nd、Sc、Y、Pr、Eu、Yb、Dy、Er、Gd、Sm、Tb、Ho、Tm、Pb、Cd、As、Cr、Ni、Ba、Sb、Sn、Bi、Sr、Li、Rb、Cs、Ti、Zr、Ag、W、Ir、Ga 44种矿物元素在枸杞成熟阶段存在“时间差异”。
表2 不同成熟期枸杞中矿物元素含量 单位:mg/kg
Table 2 Mean values of mineral elements in different maturation stages of Lycium barbarum
元素分类成熟期幼果期青果期黄变前期黄变后期红熟期微量元素Cu15.6±0.4a16.2±0.4a13.4±0.000 6b12.7±0.02b13.3±0.6bCo0.21±0.008b0.22±0.009a0.13±0.003c0.12±0.004d0.11±0.002dB13.6±0.3b14.2±0.000 6a12.3±0.1c11.2±0.3d9.6±0.1eMg3 322±100a3 030±20b2 222±39c1 814±75d1 522±39eMn29.3±0.7a28.4±0.008a19.0±0.3b16.8±0.8c16.2±0.1cMo0.53±0.01b0.62±0.02a0.48±0.02c0.41±0.000d0.32±0.005eFe83.7±2.8a66.6±0.2b50.2±0.2d50.0±0.6d58.4±0.6cZn29.5±1.1b28.5±0.3b24.3±0.1c22.8±0.3c56.2±0.8aCa3 231±42a2 358±9.5b1 623±9.5c1 344±42d1 024±4.5eSe0.020±0.008a0.026±0.002a0.021±0.001a0.026±0.005a0.029±0.003aV0.10±0.004a0.080±0.005b0.053±0.002d0.056±0.002d0.073±0.004c稀土元素La0.16±0.02a0.068±0.005b0.069±0.02b0.069±0.01b0.074±0.01bCe0.089±0.000 6a0.068±0.004b0.039±0.004d0.039±0.001d0.058±0.004cNd0.034±0.001a0.027±0.004b0.017±0.002c0.016±0.002c0.024±0.000 0bSc0.026±0.001a0.026±0.004a0.029±0.004a0.019±0.004b0.015±0.000 0bY0.020±0.000 6a0.017±0.001b0.013±0.000 6c0.013±0.000 00c0.018±0.002bPr0.007 3±0.000 6a0.006 0±0.000 00b0.004 3±0.000 6c0.004 0±0.000 00c0.006 3±0.000 6bEu0.001 1±0.000 2a0.001 1±0.000 06a0.000 77±0.000 2a0.000 33±0.000 3b0.000 87±0.000 1aYb0.002 2±0.000 3a0.002 7±0.000 6a0.001 2±0.000 2b0.001 1±0.000 1b0.002 0±0.000 1aDy0.005 8±0.000 3a0.004 0±0.000 4b0.003 1±0.000 1c0.002 7±0.000 4c0.003 9±0.000 2bEr0.002 1±0.000 3a0.001 7±0.000 3a0.001 1±0.000 2b0.001 0±0.000 1b0.001 8±0.000 2aGd0.007 0±0.000 2a0.005 7±0.000 7b0.003 0±0.000 4d0.003 0±0.000 1d0.004 8±0.000 2cSm0.004 1±0.000 2a0.003 2±0.000 2b0.002 2±0.000 2c0.002 0±0.000 2c0.003 1±0.000 2bTb0.000 98±0.000 05a0.000 79±0.000 2b0.000 68±0.000 3c0.000 68±0.000 3c0.000 71±0.000 2bHo0.000 82±0.000 2a0.000 75±0.000 2a0.000 24±0.000 2b0.000 24±0.000 2b0.000 69±0.000 3abTm0.000 20±0.000 1a0.000 16±0.000 01ab0.000 09±0.000 08c0.000 10±0.000 01c0.000 14±0.000 01bLu0.000 17±0.000 2a0.000 17±0.000 2a0.000 12±0.000 1a0.000 10±0.000 06a0.000 13±0.000 1a稀土元素总量0.3610.2320.1840.1710.213 有害金属元素Pb0.28±0.000 7a0.16±0.007b0.084±0.01b0.12±0.04b0.18±0.07bCd0.072±0.005a0.074±0.002a0.058±0.000 1b0.052±0.001bc0.046±0.004cCr3.38±0.08a1.98±0.06b1.18±0.05c0.97±0.01d0.90±0.06dAs0.035±0.000a0.030±0.003ab0.020±0.001c0.020±0.003c0.026±0.005bHg0.002 8±0.000 5a0.002 2±0.000 9a0.002 2±0.000 2a0.001 9±0.000 01a0.002 3±0.000 1aNi1.74±0.03b1.94±0.11a1.72±0.01b1.54±0.04c1.38±0.04dBa2.36±0.09a1.84±0.06b0.99±0.06c0.85±0.04c0.90±0.01cSb0.012±0.001a0.009 7±0.001ab0.006 7±0.002b0.008 6±0.001ab0.008 6±0.001abTl0.006 6±0.000 1a0.006 6±0.000 4a0.006 0±0.000 01a0.005 6±0.000 3a0.006 0±0.000 2aSn0.018±0.000ab0.023±0.000a0.013±0.000 1bc0.008 7±0.000 1c0.012±0.000bc其他元素Bi0.002 4±0.000 3a0.002 0±0.000 7ab0.001 0±0.000 2b0.001 5±0.000 4ab0.001 4±0.000 5abBe0.002 1±0.000 1a0.001 0±0.000 1a0.001 2±0.001a0.000 6±0.000 01a0.002 1±0.001aSr20.0±0.6a15.6±0.1b14.1±0.8c12.1±0.003d7.4±0.2eLi1.68±0.03c1.57±0.00d1.93±0.007a1.84±0.01b1.60±0.007dRb14.6±0.6a14.5±0.4a12.1±0.3b11.5±0.1b11.3±0.2bCs0.023±0.001a0.021±0.001a0.014±0.001b0.013±0.001b0.016±0.001bTi7.65±0.4a6.52±2.4a3.20±0.8b3.00±0.4b2.64±0.1bZr0.022±0.008a0.023±0.003a0.012±0.004c0.011±0.001c0.015±0.006bNb0.009 5±0.002a0.012±0.008a0.005 5±0.001a0.004 5±0.001a0.006 0±0.001aPd0.001 9±0.000 4a0.002 0±0.000 5a0.002 2±0.002a0.002 3±0.000 7a0.002 0±0.000 4aAg0.012±0.000 7a0.009 8±0.000 9b0.006 8±0.000 1c0.007 8±0.000 5c0.010±0.000 7bTa0.004 1±0.001a0.003 0±0.000 8a0.005 0±0.000 4a0.004 4±0.000 4a0.005 3±0.001aHf0.001 9±0.000 4a0.002 0±0.000 5a0.002 2±0.002a0.002 3±0.000 7a0.002 0±0.000 4aW0.005 4±0.000 1a0.003 6±0.000 4b0.003 2±0.000 2bc0.002 8±0.000 1c0.002 7±0.000 3cIr0.006 2±0.001ab0.007 2±0.001a0.004 8±0.000 1ab0.005 4±0.002ab0.003 3±0.000 4bPt0.007 8±0.002a0.008 2±0.000 2a0.006 0±0.000 6a0.007 1±0.003a0.004 2±0.000 3aGa0.085±0.006a0.066±0.005b0.042±0.001c0.038±0.002c0.040±0.001cGe0.002 2±0.000 5a0.001 9±0.000 1a0.000 95±0.000 2a0.000 95±0.000 4a0.002 2±0.000 8a
注:同一行的不同小写字母及其组合(a~e)表示有显著性差异(P<0.05)
2.1.2 枸杞成熟过程中矿物元素变化规律
对部分元素含量进行标准化处理(Pr、Eu、Yb、Dy、Er、Gd、Sm、Pb元素含量乘以10,Co、Mo、V、Cd、As、Sb、Tb、Ho、Tm、Lu、Sn、Cs、Zr、Ga元素含量乘以100,Se、Hg、Tl、Bi、Be、Nb、Pd、Ag、Ta、Hf、W、Ir、Pt和Ge元素含量乘以1 000;Sr、Rb元素含量除以10,Ca、Mg元素含量除以100)后,做图,见图1。
微量元素中,Cu、Co、B、Mg、Mn、Mo、Ca元素在5个成熟期基本呈降低的趋势,在幼果期/青果期达到最大值,在红熟期达到最低;Fe、Zn、V元素含量呈先降低后升高的趋势,Fe、V在幼果期达到最高,随后逐渐降低,在黄变前期/后期达到最低值,在红熟期又有所升高,达到次高值,Zn在幼果期达到次高值,随后降低在黄变后期达到最低值,Zn在红熟期急剧升高,达到最高值;Se元素在5个成熟阶段含量变化趋势不明显。
稀土元素中,Sc元素呈先升高后降低的趋势,在黄变前期达到最高值;枸杞对其他15种稀土元素在整个成熟阶段具有相同的吸收累积规律,即先降低后升高的趋势,在黄变前期/黄变后期达到最低值,在幼果期/青果期达到最大值。稀土元素总量也呈先降低后升高的趋势,在幼果期含量达到最高,在黄变前/后期达到最低值。有研究表明,稀土元素具有生物活性,可提高作物的产量影响植物生理和生物化学反应的进程对植物的生命活动是有益的[27],而GB 2762—2005《食品中污染物限量》中规定了稻谷、玉米、小麦、绿茶等植源性食品中稀土元素总含量限量值为0.5~2.0 mg/kg,但在GB 2762—2017《食品中污染物限量》中又取消了食品中稀土元素的含量限值,中宁枸杞在整个成熟阶段稀土元素总含量在幼果期达到最高,为0.361 mg/kg,根据GB 2762—2005中的规定,中宁枸杞中稀土元素总量未超过限量值。
有害金属元素中,Pb、As、Hg、Sb、Tl、Sn元素在5个成熟期呈先降低后升高的趋势;Cd、Cr、Ni、Ba元素逐渐降低,并在后期趋于稳定,说明Cd、Cr、Ni、Ba、Sn元素的富集系数较低,或是从根茎部向枸杞果实的转运能力较弱。
其他元素中,枸杞对Bi、Be、Cs、Zr、Nb、Ag、Ta、Ge元素的吸收具有一致规律,即先降低后升高,枸杞中这些元素在黄变前期/黄变后期元素含量达到最低值;Sr、Rb、Ti、W、Ga元素具有一致的吸收规律,即枸杞在幼果期、青果期对这些元素具有较强的吸收作用,随着枸杞成熟,枸杞对这些元素的需求量减少,其含量逐渐降低;Li、Pd、Hf、Zr、Pt元素含量呈现先升高后降低的趋势,基本在黄变前期/黄变后期达到最大值。
A-微量元素;B-稀土元素;C-有害金属元素;D-其他元素
图1 不同成熟阶段枸杞中微量元素、稀土元素、有害金属元素、其他元素含量的动态变化
Fig.1 Dynamic changes of contents of trace elements, rare earth elements, harmful metal elements and other elements in Lycium barbarum at different maturation stages
2.2.1 主成分及主因子提取过程分析
主成分分析是通过降维以达到简化分析过程的目的,原始指标的因子载荷和主成分得分图分别表明其对主成分的贡献程度和不同样本之间的联系与差异[28]。
对中宁枸杞5个成熟阶段中55种元素含量进行PLS-DA分析,结果列于表3,根据特征值大于1的原则,可提取出2个主成分,2个主成分的特征值分别为9.76、2.38,方差贡献率分别为65.1%和15.8%,累计方差贡献率为80.9 %。因此可以将PLS-DA的前2个主因子作为评价枸杞样品质量的综合变量,达到了对评价枸杞矿物元素含量指标的降维目的。
表3 PLS-DA的累计贡献率
Table 3 Cumulative contribution rate of PLS-DA
主成分特征值方差贡献率/%累积方差贡献率/%19.7665.165.122.3815.880.930.9576.3887.340.5003.3390.650.4513.0193.6
2.2.2 基于PLS-DA的不同成熟阶段枸杞的分类
中宁枸杞5个成熟阶段中55种元素含量进行PLS-DA分析,构建不同成熟阶段枸杞的区分模型。该模型中R2X(cum)、R2Y(cum)和Q2(cum)分别为0.936、0.977和0.912,说明建立的PLS-DA模型中2个主成分能有效解释5个成熟阶段之间的差异,且该模型具有一定的预测能力。图2为PLS-DA模型第一、二主成分得分图。可以看出,各成熟阶段样品群体有明显的聚集趋势,幼果期(YG)枸杞样品聚在第一象限,红熟期(HS)枸杞样品聚在第二象限,黄变前期(HBQ)和黄变后期(HBH)枸杞样品聚在第三象限,青果期(QG)枸杞样品聚在第四象限。
根据VIP值(图3),可筛选出Sc(VIP:1.576 3)、Li(1.501 58)、Zn(1.450 72)、Ta(1.421 99)、Be(1.373 54)、La(1.325 38)、Se(1.305 46)、Sb(1.276 32)、Lu(1.251 62)、Pt(1.241 67)、Ni(1.199 67)、Ir(1.146 07)、Sr(1.135 21)、Hg(1.108 52)、Mo(1.100 46)、Bi(1.095 28)、Pb(1.046 19)、W(1.035 4)、Hf(1.034 67)、Yb(1.019 33)、Pd(1.008 18)、Sn(1.005 84)22种矿物元素作为中宁枸杞不同成熟阶段的有效区分指标。
图2 PLS-DA模型第一、二主成分得分图
Fig.2 Score plot of first versus second components in PLS-DA
图3 PLS-DA模型VIP值
Fig.3 The VIP plot of PLS-DA
追踪中宁枸杞中55种矿物元素在整个成熟过程中的动态变化,发现枸杞在不同成熟阶段元素含量差异较大,44种矿物元素含量在不同成熟阶段具有显著差异(P<0.05),且稀土元素显示出轻稀土元素富集、重稀土元素亏欠的特征。除Se元素外,其余54种矿物元素在整个成熟期有其各自的变化趋势,Cu、Co、B、Mg、Mn、Mo、Ca、Cd、Cr、Ni、Ba、Sr、Rb、Ti、W、Ga呈降低趋势,Fe、Zn、V、Pb、As、Hg、Sb、Tl、Sn、La、Ce、Nd、Y、Pr、Eu、Yb、Dy、Er、Gd、Sm、Tb、Ho、Tm、Lu、Bi、Be、Cs、Zr、Nb、Ag、Ta、Ge呈先降低后升高的趋势,Sc、Li、Pd、Hf、Zr、Pt呈先升高后降低的趋势。对中宁枸杞5个成熟阶段中55种元素含量进行PLS-DA分析,可提取出2个主成分,累计方差贡献率为80.9%,同时构建不同成熟阶段枸杞的区分模型,第一主成分、第二主成分得分图可以将5个成熟阶段的样品很好地区分开,同时可筛选出Sc、Li、Zn、Ta、Be、La、Se、Sb、Lu、Pt、Ni、Ir、Sr、Hg、Mo、Bi、Pb、W、Hf、Yb、Pd、Sn 22种矿物元素为中宁枸杞不同成熟阶段的有效区分指标。
[1] 路安民, 王美林.关于中药现代化中的物种鉴定问题——基于枸杞分类和生产问题的讨论[J].西北植物学报, 2003, 23(7):1 077-1 083.
LU A M, WANG M L.On the identification of the original plants in the modernization of Chinese herbal medicine——An example from the taxonomy and exploitation of ‘Gouqi’[J].Acta Botanica Boreali-Occidentalia Sinica,2003, 23(7):1 077-1 083.
[2] 如克亚·加帕尔, 孙玉敬, 钟烈州, 等.枸杞植物化学成分及其生物活性的研究进展[J].中国食品学报, 2013, 13(8):161-172.
RUKEYA·JIAPAER, SUN Y J, ZHONG L Z, et al.A Review of phytochemical composition and bio-active of Lycium barbarum fruit (Goji)[J].Journal of Chinese Institute of Food Science and Technology, 2013, 13(8):161-172.
[3] 葛影影, 何国戈, 郑经成, 等.枸杞的生物学功能及其在动物生产中的应用[J].安徽农业科学, 2021,49(3):25-27.
GE Y Y, HE G G, ZHENG J C, et al.The biological function of Lycium barbarum and its application in animal production[J].Journal of Anhui Agricultural Sciences, 2021,49(3):25-27.
[4] 杨永利, 明磊国, 林浩, 等.枸杞养肝明目功效研究进展[J].中国食物与营养, 2015, 21(7):75-78.
YANG Y L, MING L G, LIN H, et al.Research progress of liver and eyesight function of wolfberry[J].Food and Nutrition in China, 2015, 21(7):75-78.
[5] XU L, SHI J L, CHEN D F, et al.Optimization of fermentation conditions for production of Chinese wolfberry fruit wine with higher total flavonoid content[J].Food Science and Technology, 2011, 32:188-193.
[6] YE C H.Technological study on health food of wolfberry tofu[J].Food Research and Development, 2010, 31:113-115.
[7] 雷建武, 米佳, 罗青, 等.枸杞中类胡萝卜素及体外抗氧化活性研究[J].食品工业, 2015, 36(12):5-8.
LEI J W, MI J, LUO Q, et al.Composition and antioxidant activity of carotenoids from different varieties of Lycium[J].The Food Industry, 2015, 36(12):5-8.
[8] 杨春霞, 李彩虹, 杨静, 等.宁夏产区不同品种枸杞中总酚含量分析[J].食品研究与开发, 2020,41(8):183-187.
YANG C X, LI C H, YANG J, et al.Analysis of polyphenols content in different varieties of Lycium barbarum L.in Ningxia producing area[J].Food Research and Development, 2020, 41(8):183-187.
[9] 王亚军, 梁晓婕, 张波, 等.产地差异对宁夏枸杞果实形态及糖分含量的影响[J].干旱地区农业研究, 2018, 36(5):68-75;81.
WANG Y J, LIANG X J, ZHANG B, et al.Effects of different habitats on the fruit morphology and sugar content of Lycium barbarum L.[J].Agricultural Research in the Arid Areas, 2018, 36(5):68-75;81.
[10] 刘小云, 李拥军, 崔琴.甘肃红芪中5种微量元素含量及健康风险分析[J].疾病预防控制通报, 2018, 33(1):8-11.
LIU X Y, LI Y J, CUI Q.Analysis of content of five microelements in hedysari Radix from Gansu province and their associated health risks[J].Bulletin of Disease Control & Prevention(China), 2018, 33(1):8-11.
[11] 强正泽, 王燕, 王明伟, 等.红芪不同部位中13种无机元素分布规律研究[J].天然产物研究与开发, 2017,29(4):623-628;652.
QIANG Z Z, WANG Y, WANG M W, et al.Distribution of 13 inorganic elements in different parts of hedysari Radix[J].Natural Product Research and Development, 2017,29(4):623-628;652.
[12] KONIECZYNSKI P, ARCEUSZ A, WESOLOWSKI M.Essential elements and their relations to phenolic compounds in infusions of medicinal plants acquired from different European regions[J].Biological Trace Element Research, 2016,170(2):466-475.
[13] 杨波, 王振国.植物类中药寒热药性与无机元素相关性研究[J].南京中医药大学学报, 2011, 27(2):109-111.
YANG B, WANG Z G.Literature study on relationship between cold-heat nature of plant-based Chinese herbs and inorganic elements[J].Journal of Nanjing University of Traditional Chinese Medicine, 2011, 27(2):109-111.
[14] 钱丽丽, 邱彦超, 李殿威, 等.基于产地、品种和年份影响矿物元素含量的大米判别[J].食品科学, 2021,42(16):322-327.
QIAN L L, QIU Y C, LI D W, et al.Rice discrimination based on the influence of origin, variety and year on mineral element content[J].Food Science, 2021,42(16):322-327.
[15] SAYAGO A, GONZLEZ-DOMNGUZE R, BELTRN R, et al.Combination of complementary data mining methods for geographical characterization of extra virgin olive oils based on mineral composition[J].Food Chemistry, 2018, 261:42-50.
[16] BANDONIENE D, ZETTL D, MEISEL T, et al.Suitability of elemental fingerprinting for assessing the geographic origin of pumpkin (Cucurbita pepo var.styriaca) seed oil[J].Food Chemistry, 2013, 136(3-4):1 533-1 542.
[17] 连思雨, 谢瑜杰, 张紫娟, 等.多元素分析结合化学计量学方法快速判别宁夏和青海枸杞[J].食品与发酵工业, 2020, 46(13):250-254.
LIAN S Y, XIE Y J, ZHANG Z J, et al.Rapid discrimination of Lycium barbarum L.from Ningxia and Qinghai based on multi-element analysis combined with chemometrics[J].Food and Fermentation Industries, 2020, 46(13):250-254.
[18] 开建荣, 王彩艳, 王晓静, 等.微波消解-ICP-MS法测定枸杞中49种无机元素含量[J].宁夏农林科技, 2020, 61(2):55-59.
KAI J R, WANG C Y, WANG X J, et al.Determination of 49 inorganic elements in Lycium barbarum by microwave digestion-ICP-MS[J].Ningxia Journal of Agriculture and Forestry Science and Technology, 2020, 61(2):55-59.
[19] SNCHEZ T, MARTN S, SACO D.Some responses of two Nicotiana tabacum L.cultivars exposed to vanadium[J].Journal of Plant Nutrition, 2014, 37(5):777-784.
[20] HOU M, LI M Y, YANG X H, et al.Responses of nonprotein thiols to stress of vanadium and mercury in maize(Zea mays L.)seedlings[J].Bulletin of Environmental Contamination and Toxicology, 2019, 102(3):425-431.
[21] 傅文卫. 贵州道地药材吴茱萸的研究——微量元素与药材的道地性[J].贵阳中医学院学报, 1995(3):11-13.
FU W W.Study on Evodia rutaecarpa in Guizhou-trace elements and the genuine nature of the medicinal materials[J].Journal of GuiZhou University of Traditional Chinese Medicine, 1995(3):11-13.
[22] 牛艳, 王明国, 郑国琦, 等.宁夏不同地域枸杞子微量元素比较研究[J].干旱地区农业研究, 2005(2):100-103.
NIU Y, WANG M G, ZHENG G Q, et al.Comparison of the trace elements in the Licum bararum L.fruits of different regions of Ningxia[J].Agricultural Research in the Arid Areas, 2005(2):100-103.
[23] 彭励, 王林.不同产地枸杞子中微量元素和黄酮含量的比较[J].微量元素与健康研究, 2007, 24(5):14-16.
PENG L, WANG L.Comparison of trace elements and total flavone content in Chinese wolfberry in different regions[J].Studies of Trace Elements and Health, 2007, 24(5):14-16.
[24] 王立军, 王玉琦, 章申, 等.中国不同类型土壤中稀土元素的形态分布特征[J].中国稀土学报, 1997(1):64-70.
WANG L J, WANG Y Q, ZHANG S, et al.Speciation of rare earth elements in different types of soils in China[J].Journal of the Chinese Society of Rare Earths, 1997(1):64-70.
[25] 齐国亮, 苏雪玲, 王俊, 等.宁夏枸杞主产区土壤和果实中稀土元素含量及其相关性[J].南方农业学报, 2014,45(7):1 206-1 210.
QI G L, SU X L, WANG J, et al.Contents of rare earth elements in soils and fruits of Lycium bararum L.in different regions and their relationship[J].Journal of Southern Agriculture, 2014, 45(7):1 206-1 210.
[26] 朱美霖, 张亚红, 朱慧丽, 等.不同饮水途径下重金属暴露引起的健康风险比较[J].能源研究与管理, 2019(2):47-49;65.
ZHU M L, ZHANG Y H, ZHU H L, et al.Comparison of health effects based on heavy metals exposure in drinking water by different ways[J].Energy Research and Management, 2019(2):47-49;65.
[27] 吴晶, 冯秀娟.稀土元素对植物的生物有效性研究进展[J].中国资源综合利用, 2012, 30(4):37-40.
WU J, FENG X J.Progress in the plant ecological effects of rare earth elements[J].China Resources Comprehensive Utilization, 2012, 30(4):37-40.
[28] LIU Z C, WANG L P, LIU Y M.Analyzing differences in freshness of SA-1 hops by headspace solid-phase microextraction gas chromatography-mass spectrometry combined with chemometrics[J].Journal of the American Society of Brewing Chemists, 2017, 75(3):193-200.