薄荷是唇形科薄荷属多年生草本植物,药食两用植物。具有利咽、透疹等功效,用于咽喉肿痛及小儿麻疹等症[1-3]。现代药理研究表明,薄荷能够抗肿瘤、抗菌、消炎、有效缓解镇痛[4-7]。薄荷的药用价值很高,其营养成分也十分丰富,富含氨基酸、有机元素、维生素等多种营养成分,具有重要的营养价值和保健价值[8-10]。
目前,对于薄荷研究多以挥发油类成分为主[11-13],而作为一种药食两用的植物,薄荷作为常见的茶饮,其实际运用中不仅仅局限于挥发油, 其氨基酸和核苷类成分同样重要。而对薄荷中氨基酸、核苷类成分的研究较少,其营养价值尚未得到充分认识。氨基酸主要以游离及水解的形式存在于植物体内,氨基酸检测方法主要包括高效液相色谱法[14]、毛细管电泳法[15]、离子交换色谱法等[16]。本试验通过超快速液相色谱-质谱联用技术(ultra-fast liquid chromatography-mass spectrometry,UFLC-MS/MS),建立同时测定薄荷中21种氨基酸和10种核苷类成分含量方法,对3个不同产地的薄荷中氨基酸和核苷类成分差异进行分析比较,并应用PCA[17]和TOPSIS[18]方法对样品进行综合评价,为药食两用薄荷的质量评价和综合利用提供理论依据。
Trip Quad 4500串联四级杆质谱仪,美国 ABSciex公司;ExionLC AD 超快速液相色谱仪,日本Shimadzu公司;KQ-400KDE型高功率数控超声清洗器(100 kHz),昆山市超声仪器有限公司;离心机(HC-3018),安徽中科中佳科学仪器有限公司;SArtorius BT125D十万分之一电子天平,德国赛多利斯公司;ME36S 型百万分之一电子天平,德国赛多利斯公司;LGJ-10真空冷冻干燥机,北京松源华兴科技发展有限公司;Milli-Q超纯水系统,Millipore公司;0.22 μm 微孔滤膜,津腾,中国天津。
异亮氨酸(Ile)、组氨酸(His)、苯丙氨酸(Phe)、亮氨酸(Leu)、谷氨酸(Glu)、络氨酸(Tyr)、丙氨酸(Ala)、胱氨酸(Cys)、苏氨酸(Thr)、甲硫氨酸(Met)、脯氨酸(Pro)11种氨基酸,均购于中国药品生物制品检定所,批号为140624-200805。天冬氨酸(Asp,批号B21934)、天冬酰胺(Asn,批号B21935)、丝氨酸(Ser,批号140688)、甘氨酸(Gly,批号B21915)、缬氨酸(Val,批号140681)、精氨酸(Arg,批号B21920)、羟基脯氨酸(Hpro,批号B21928)、色氨酸(Tryptophan,批号140677)、盐酸赖氨酸(Lysine,批号140693)、γ-氨基丁酸(GABA,批号B21979)、尿嘧啶(Uracil,批号B20908)、腺嘌呤(Adenine,批号B21357)、胞苷(Cytidine,批号B20073)、肌苷(Inosine,批号B20582)、胸苷(Thymidine,批号B21915)、鸟苷(Guanine,批号B20905)、腺苷(Adenosine,批号B21356)、尿苷(Uridine,批号B20907)、次黄嘌呤(Hypoxanthine,批号B20211)、鸟嘌呤(Guanine,批号B20906),均购于上海源叶生物科技有限公司,质量分数均大于98%。水为超纯水;乙腈、甲酸为质谱纯,德国Merck公司;其余试剂为分析纯。
试验所研究的样品,采集于3个不同产地的薄荷植物新鲜叶片,经安徽中医药大学方成武教授鉴定为唇形科植物薄荷Mentha haplocalyx Briq.的干燥地上部分。样品信息如表1所示。
表1 薄荷主辅料样品来源及基本信息
Table 1 Source and basic information of main products of mint steamed vegetables
样品编号样品名称产地收集方式收集日期S1薄荷安徽省阜阳市实地采收2019.7S2薄荷安徽省阜阳市实地采收2019.7S3薄荷安徽省亳州市实地采收2019.7S4薄荷安徽省亳州市实地采收2019.7S5薄荷河南省商丘市实地采收2019.7S6薄荷河南省商丘市实地采收2019.7
1.3.1 色谱条件
Waters XBridge Amide专用柱(100 mm×2.1 mm,3.5 μm);流动相∶水(含0.2%甲酸)(A)-乙腈(含0.2%甲酸)(B),梯度洗脱(0~2.5 min,15% A;2.5~5 min,15%~50% A;5~7 min,50% A;7~8 min,50%~15% A;8~11 min,15% A);柱温 30 ℃,流速0.2 mL/min,进样量1 μL。其正离子多反映监测色谱图见图1。
图1 21种氨基酸和10种核苷的MRM图
Fig.1 UFLC-MS/MS multiple-reaction monitoring(MRM) chromatograms of 21 amino acids and 10 nucleosides
1.3.2 质谱条件
离子化模式ESI+下,采取多反应监测(MRM)检测;喷雾电压5.5 kV;离子化温度550 ℃;气帘气为241.3 kPa(35 psi),Gas1为379.2 kPa(55 psi),Gas2 为379.2 kPa(55 psi);同时对碰撞电压、去簇电压和检测离子对进行优化。31种目标成分的优化质谱检测条件参数见表2。
表2 31种目标成分的优化质谱检测条件参数
Table 2 Optimized MS/MS spectrometry detection condition parameters of 31 target components
化合物保留时间/min碰撞电压/eV去簇电压/VMRN离子对(m/z)天冬氨酸6.051059134.05/87.96天冬酰胺6.281346132.80/115.70谷氨酸6.781483147.08/83.92缬氨酸3.301054118.09/72.06甲硫氨酸3.411091150.06/104.03丝氨酸6.18867106.05/59.99甘氨酸5.6367376.04/30.00丙氨酸4.8107990.06/44.02异亮氨酸2.751566132.00/86.00亮氨酸2.621098132.10/86.05苯丙氨酸2.841456166.10/120.05胱氨酸5.611871240.80/151.90组氨酸7.131695156.08/110.03精氨酸6.681888175.12/70.02苏氨酸5.81154120.30/76.80酪氨酸3.881746182.16/136.08羟脯氨酸2.662552133.80/71.80脯氨酸4.381068116.07/70.02色氨酸2.7215202205.00/188.00γ—氨基丁酸3.181432103.70/86.90盐酸赖氨酸6.502452147.00/84.00尿嘧啶2.0321111113.04/70.00腺嘌呤3.352451136.06/136.00胞苷5.661061244.09/112.00肌苷3.721546269.00/137.05鸟苷5.231562284.30/152.00胸苷2.081361243.10/127.07腺苷3.602386268.10/136.07尿苷2.8413103244.90/113.00鸟嘌呤4.552551152.00/135.00次黄嘌呤2.782451137.05/137.05
1.3.3 溶液制备
1.3.3.1 对照品溶液制备
精密称取天冬氨酸、天冬酰胺、丝氨酸、丙氨酸、甘氨酸、缬氨酸、谷氨酸、甲硫氨酸、异亮氨酸、精氨酸、组氨酸、胱氨酸、苯丙氨酸、亮氨酸、苏氨酸、酪氨酸、羟脯氨酸、脯氨酸、色氨酸、γ-氨基丁酸、盐酸赖氨酸、尿嘧啶、腺嘌呤、胞苷、肌苷、鸟苷、胸苷、腺苷、尿苷、鸟嘌呤、次黄嘌呤对照品适量,分别置25 mL容量瓶中,加超纯水配制成0.16、0.18、0.12、0.12、0.14、0.13、0.11、0.14、0.18、0.12、0.13、0.13、0.13、0.12、0.11、0.15、0.13、0.11、0.17、0.17、0.15、0.17、0.13、0.11、0.15、0.12、0.15、0.16、0.11、0.13、0.13 mg/mL的对照品储备液。取各对照品储备液适量,加水配制成混合对照品溶液,逐级稀释,得到一系列不同质量浓度的混合对照品溶液,于4 ℃保存备用。
1.3.3.2 供试品溶液的配制
称取鲜品薄荷叶500 g,用水清洗干净,均匀铺散在托盘上,覆盖上已扎孔的薄膜,置于-20℃冰箱预冻,随后将其置于真空冷冻干燥机[19-21],直至水分含量达到8%以下,即干燥完成。取薄荷粉末(过40目筛)1.000 g,置于100 mL容量瓶中,加入20 mL超纯水,室温超声处理(100 kHz)提取30 min,吸取上清液于4 500 r/min转速下离心10 min,过0.22 μm 微孔滤膜,供分析。
1.3.4 线性关系考察、检测限和定量限
精密吸取“1.3.3”项下不同浓度系列对照品溶液及混合对照品储备液各1 μL,按照“1.3.1”和“1.3.2”项下色谱-质谱条件进样测定,以对照品的峰面积积分值(Y)为纵坐标,对照品的质量浓度(X)为横坐标,得出回归方程、相关系数(r)、线性范围。以信噪比(S/N)约为3 计算各成分的检测限(LOD),以信噪比(S/N)约为10,计算各成分的定量限(LOQ),结果见表3。
1.3.5 方法学考察
1.3.5.1 精密度试验
精密吸取一定浓度的混合对照品溶液,按照“1.3.1”和“1.3.2”项下色谱-质谱条件连续进样6次,得出21种氨基酸和10种核苷对照品峰面积的RSD在0.20%~3.36%。具体数据见表4。
1.3.5.2 重复性试验
精密称取S5—薄荷样品6份,每份1.000 g,按照“1.3.3.2”项下方法制备供试品溶液,按照“1.3.1”和“1.3.2”项下色谱-质谱条件进样测定,计算21种氨基酸和10种核苷的质量分数的RSD在0.21%~3.60%。具体数据见表4。
表3 21种氨基酸和10种核苷的线性关系考察
Table 3 Investigation of liner relationship of 20 amino acids and 10 nucleosides
化合物回归方程线性范围/(μg·mL-1)R2LOD/(ng·mL-1)LOQ/(ng·mL-1)天冬氨酸y=2×106x+875541.02~16.300.99942.759.16天冬酰胺y=386857x+102601.14~18.200.99984.0513.52谷氨酸y=4206x-319.640.71~11.300.99920.140.18缬氨酸y=2×107x-8532080.86~13.800.99961.695.65甲硫氨酸y=307086x+3472.20.87~13.900.99963.2010.65丝氨酸y=474974x-127520.73~11.600.99930.812.70甘氨酸y=6517.3x-472.370.90~14.40.99916.7922.64丙氨酸y=345588x+381110.78~12.500.99903.0710.25异亮氨酸y=2×107x+4×1061.12~17.900.99974.3114.38亮氨酸y=1×106x-991440.76~12.200.99906.9623.22苯丙氨酸y=5×107x-4×1060.78~12.500.99924.4714.90胱氨酸y=307.47x-20.2380.79~12.600.99970.160.55组氨酸y=230.35x-13.9440.79~12.600.99930.782.60精氨酸y=1×106x-1441770.77~12.300.99923.5011.66苏氨酸y=70795x-7347.70.71~11.40.99912.056.83酪氨酸y=7×106x-5917750.91~14.600.99922.518.36羟脯氨酸y=1048.6x-125.770.79~12.700.99900.220.73脯氨酸y=8×106x-4171700.71~11.300.99942.618.72色氨酸y=502153x-478101.03~16.500.99915.7319.11γ-氨基丁酸y=1×107x+7056901.06~16.900.99953.8212.73盐酸赖氨酸y=1×107x-2936970.93~14.800.99951.715.70尿嘧啶y=190353x-6272.11.06~17.000.99953.2410.80腺嘌呤y=5×107x-3×1060.81~13.000.99961.655.50胞苷y=8×106x-6714370.71~11.400.99904.5415.13肌苷y=1×107x-1541470.96~15.400.99944.8016.01鸟苷y=8×106x+5518270.78~12.400.99976.0220.07胸苷y=104245x-6960.20.92~14.700.99915.7719.23腺苷y=2×107x-7636390.99~15.900.99923.3111.04尿苷y=71315x-6299.80.71~11.300.99904.1313.78鸟嘌呤y=64087x-235.210.81~13.000.99943.9213.07次黄嘌呤y=1×107x+5384360.82~13.100.99952.387.95
1.3.5.3 稳定性试验
取S5—薄荷样品的供试品溶液,分别在 0、2、4、8、12、24 h时按照“1.3.1”和“1.3.2”项下色谱-质谱条件下进样分析,测定21种氨基酸和10种核苷成分峰面积的RSD值,RSD在0.65%~3.84%。具体数据见表4。
1.3.5.4 加样回收率试验
精密称取S5—薄荷样品6份,每份1.000 g,分别加入与1.000 g样品中各待测成分含量相当的对照品,按照“1.3.3.2”项下方法制备加样回收供试品溶液,按照“1.3.1”和“1.3.2”项下色谱-质谱条件下进样测定,计算加样回收率和相对标准偏差,结果见表4。31种成分的平均加样回收率在94.18%~98.76%,相对标准偏差在1.17%~3.56%。
表4 21种氨基酸和10种核苷方法学考察 单位:%
Table 4 Methodological study of 21 amino acids and 10 nucleosides
化合物重复性RSD稳定性RSD精密度RSD平均回收率平均回收率RSD天冬氨酸1.321.870.5795.022.18天冬酰胺0.211.240.2397.302.56谷氨酸2.873.541.3297.542.55缬氨酸1.781.110.9195.412.91甲硫氨酸1.463.150.4695.582.29丝氨酸0.843.562.7795.572.67甘氨酸2.192.582.3996.522.82丙氨酸2.482.392.1894.62.01异亮氨酸0.531.571.3194.381.92亮氨酸2.831.140.8394.392.18苯丙氨酸1.431.643.2997.391.69胱氨酸2.891.932.9297.752.44组氨酸2.642.332.1895.541.37精氨酸1.122.980.6195.212.61苏氨酸3.603.013.3696.33.03酪氨酸1.993.141.4196.562.65羟脯氨酸1.892.042.1195.032.89脯氨酸2.231.151.4797.022.59色氨酸2.141.650.8495.352.61γ-氨基丁酸2.012.670.4594.942.09盐酸赖氨酸0.780.650.4796.732.67尿嘧啶1.353.150.6295.582.88腺嘌呤1.071.800.3697.152.31胞苷2.012.651.0495.452.51肌苷2.071.970.7596.642.44鸟苷1.282.240.8395.902.61胸苷2.582.830.9194.892.11腺苷2.172.051.0994.182.47尿苷1.923.841.4195.262.06鸟嘌呤0.381.040.2098.761.17次黄嘌呤0.421.580.6493.241.36
按照“1.3.3.2”项下方法制备供试品溶液,按照“1.3.1”和“1.3.2”项下色谱-质谱条件进样得相应的峰面积,根据对应的标准曲线计算出供试样品中21种氨基酸和10种核苷类成分的含量,结果见表5。
表5 薄荷中21种氨基酸和10种核苷的含量测定结果 单位:μg/g
Table 5 Determination results of 21 amino acids and 10 nucleosides in peppermint
编号天冬氨酸天冬酰胺谷氨酸缬氨酸甲硫氨酸丝氨酸甘氨酸丙氨酸异亮氨酸亮氨酸苯丙氨酸胱氨酸组氨酸精氨酸苏氨酸S11090.651473.2426225.38230.905.59235.2628.97149.69253.84248.94121.81785.712006.81271.48841S2985.341164.3225774.07201.643.67145.5730.33123.46264.86265.69138.94804.631744.56220.76885.61S31032.063194.7724275.19278.516.06349.1223.42153.55380.48220.47109.67622.392131.19560.65751.35S4905.163312.1222639.26253.685.88364.3822.31166.29369.74212.4998.75583.861923.05542.69697.8S5723.623761.1019951.69293.179.91478.1522.7380.21546.97237.51111.55543.531107.25670.31532.45S6631.553645.6717395.57312.4311.43442.4720.0970.33563.43226.73110.29527.981004.21704.83510.68编号酪氨酸羟脯氨酸脯氨酸色氨酸γ-氨基丁酸盐酸赖氨酸尿嘧啶腺嘌呤胞苷肌苷鸟苷胸苷腺苷尿苷鸟嘌呤次黄嘌呤S152.82897.144.1555.82431.44218.89486.36121.54422.103.12106.697.780.85138.9921076.6524.89S248.34785.793.5264.87453.32223.43500.65138.76443.782.9596.456.540.78152.8022417.6931.04S3137.651339.818.38253.03337.44196.99236.9093.04230.562.08116.039.060.8872.3810109.533.26S4140.781275.469.64235.46318.90183.58253.6889.76254.591.86121.1210.870.9270.4214371.282.54S5101.941984.713.90215.83178.65217.49178.4658.3182.582.176.337.610.8458.124332.893.83S686.571795.282.06198.74153.84213.28185.7960.3280.641.545.465.970.7148.053754.412.97
由于薄荷样品中各成分在量纲和数量级上存在较大差异,为更好地分析结果,将原始数据进行标准化处理后,对21种氨基酸和10种核苷含量分别应用SPSS 23.0软件进行主成分(PCA)分析。以特征值大于1,累积贡献率大于70%为提取标准。
提取了PCA初始因子载荷矩阵,因子载荷矩阵主要说明主成分在各变量上的载荷,载荷越高,越相关。在氨基酸中,甘氨酸、谷氨酸、苏氨酸、胱氨酸、γ-氨基丁酸在PC1上有较高的载荷,说明PC1能反映这5种氨基酸的信息;同理,PC2能反映脯氨酸、丙氨酸的信息。在核苷中,尿嘧啶、腺嘌呤、胞苷、尿苷、鸟嘌呤在PC1上有较高的载荷,说明PC1能反映这5种核苷成分的信息;同理,PC2能反映胸苷的信息。综上,不同产地薄荷的氨基酸类成分中,甘氨酸、谷氨酸、苏氨酸、胱氨酸、γ-氨基丁酸、脯氨酸、丙氨酸为主要成分;同理,在核苷类成分中,尿嘧啶、腺嘌呤、胞苷、尿苷、鸟嘌呤、胸苷为主要成分。
分别以氨基酸和核苷类成分以及总成分含量为分析指标,采样TOPSIS分析方法对6份样品进行综合分析,得出排序结果。
2.3.1 归一化处理
将6份样品中 21种氨基酸和 10种核苷成分的含量数据矩阵分别进行归一化处理,并建立相应矩阵,计算α值,如公式(1)所示:
(1)
式中:αij表示第i个评价对象在第j个指标上归一化处理后的取值,Xij表示第i个评价对象在第j个指标上的取值。
2.3.2 最优与最劣方案
根据归一化处理后的样本情况,得到最优(Z+)和最劣(Z-)值向量。
2.3.3 相对贴近度(Ci)的计算
相对贴近度根据公式(2)计算:
(2)
式中:表示各个评价对象与最优方案的距离;表示各个评价对象与最劣方案的距离;D为每个评价指标与最优、最劣差值平方和的开方值。
2.3.4 TOPSIS分析结果
按照Ci大小将各评价对象排序,Ci值越大,表示综合品质越好。6份样品中 21种氨基酸和 10种核苷类成分 TOPSIS 分析结果见表6。
结果表明,以氨基酸为分析指标,S4样品最佳;以核苷为分析指标,S2样品最高;以总成分为分析指标,S1样品最高。
表6 6份样品中21种氨基酸和10种核苷类成分TOPSIS分析结果
Table 6 Results of TOPSIS analysis of 21 amino acids and 10 nucleosides in 11 samples
编号总成分排名氨基酸排名核苷排名D+iD-iCiD+iD-iCiD+iD-iCiS11.05841.36700.563611.02790.67040.394750.25261.19140.82512S21.18361.45090.550721.15570.62420.350760.25541.30980.83681S31.11781.17850.513240.52480.99340.654320.98700.63400.39114S41.07851.25520.537930.49641.02300.673310.95750.72740.43173S51.42600.97820.406950.75880.96500.559831.20730.16060.11745S61.64270.93880.363760.90650.93870.508741.36990.01620.01176
各产地样品在氨基酸和核苷酸样品的组成和比例均存在一定差异。从总氨基酸含量来看,亳州产样品总氨基酸含量高达48 337.11 μg/g,阜阳次之,商丘最低。氨基酸成分以谷氨酸平均含量最高(22 710.19 μg/g),天冬酰胺、组氨酸、羟脯氨酸、天冬氨酸等含量次之,脯氨酸平均含量最低(5.28 μg/g)。从核苷酸总量来看,阜阳产样品总核苷含量为2 309.02 μg/g,亳州次之,商丘最低。核苷类成分以鸟嘌呤平均含量最高(12 677.08 μg/g),尿嘧啶、胞苷等次之,腺苷平均含量最低(0.83 μg/g)。由此说明产区环境生态因子对薄荷中氨基酸和核苷类成分的合成积累影响较大。TOPSIS综合分析结果显示, 就总成分而言, S1(安徽阜阳) 样品薄荷相对较优。
本试验建立了UFLC-MS/MS同时测定3个不同产地薄荷中21种氨基酸和10种核苷含量的方法,比较了不同产地对薄荷中氨基酸和核苷类成分的影响,综合评价显示,阜阳产薄荷的综合质量最优,亳州其次,商丘综合质量最差。通过PCA和TOPSIS对药食两用薄荷质量进行评价,为优选薄荷最优产地提供理论依据,同时为进一步探索评价药食两用薄荷的内在质量与控制提供依据及新的质量综合评价参考方法。
[1] 华燕青.薄荷化学成分及其提取方法研究进展[J].陕西农业科学,2018,64(4):83-86.
[2] 崔培超.GC快速测定薄荷药材中薄荷酮和薄荷脑的含量[J].化工管理,2019(29):47-48.
[3] 王菲菲,刘杰,何风艳,等.薄荷属3种植物及变种的性状鉴别、化学成分分析及DNA条形码研究比较[J].药学学报,2019,54(11):2 083-2 088.
[4] 彭蕴茹,钱辉,石磊,等.薄荷非挥发性提取部位的药理活性研究[J].中药材,2008,31(1):104-106.
[5] ENCK P,JUNNE F,KLOSTERHALFEN S,et al.Therapy options in irritable bowel syndrome[J].Eur J Gastroenterol Hepatol,2010,22(12):1 402-1 411.
[6] DAVIES S J,HARDING L M,BARANOWSKI A P.A novel treatment of postherpetic neuralgia using peppermint oil[J].Clin J Pain,2002,18(3):200-202.
[7] 韩永胜,武松,袁晓,等.复方薄荷酊治疗瘙痒性皮肤病600例疗效观察[J].中医外治杂志,2003,12(3):46.
[8] 李祥,邢文峰.薄荷的化学成分及临床应用研究进展[J].中南药学,2011,9(5):362-365.
[9] 沈梅芳,李小萌,单琪媛.薄荷化学成分与药理作用研究新进展[J].中华中医药学刊,2012,30(7):1 484-1 487.
[10] 张荣发,杨宗发,江尚飞.薄荷油的药理毒理作用研究进展[J].中国药业,2012,21(19):1-3.
[11] 黄兴雨,杨黎燕,尤静.薄荷挥发油研究进展[J].化工科技,2019,27(3):70-74.
[12] 魏亮,方洪壮,吴比,等.东北野生薄荷挥发油成分及抑菌活性研究[J].安徽农业科学,2019,47(10):170-172;178.
[13] 束雅春,段煜,陈亚军,等.薄荷-荆芥药对与单味药挥发性成分分析及体外抑菌作用比较[J].中国实验方剂学杂志,2019,25(7):6-13.
[14] 唐晓雷,何慧楠,赵悦名,等.不同等级的梅花鹿鹿茸药材中氨基酸含量的柱前衍生-高效液相色谱法测定[J].时珍国医国药,2019,30(3):560-563.
[15] 林秀丽.毛细管电泳在某些中药和食品成分分析中的应用研究[D].福州:福州大学,2014.
[16] 柯丽群,王同珍,曹维强,等.微波水解-茚三酮柱后衍生离子交换色谱法测定鹌鹑蛋中的17种氨基酸[J].食品安全质量检测学报,2019,10(17):5 903-5 907.
[17] 吴振,李红,杨勇,等.基于无机元素的花椒产地溯源和品种聚类分析[J].食品科学,2019,40(16):213-219.
[18] 李月,胡鹏,张园娇,等.QTRAP-UFLC-MS/MS同时测定不同产地番荔枝子中氨基酸和核苷类成分的含量及PCA和TOPSIS分析[J].中草药,2018,49(11):2 671-2 677.
[19] 扶庆权,王海鸥,陈雨,等.不同干燥方式对白玉菇品质的影响[J].食品研究与开发,2019,49(17):148-152.
[20] 侯会,陈鑫,方东路,等.干燥方式对食用菌风味物质影响研究进展[J].食品安全质量检测学报,2019,10(15):4 877-4 883.
[21] 孙晓健,于鹏飞,李晨晨,等.HS-SPME结合GC-MS分析真空冷冻干燥香椿中挥发性成分[J].食品工业科技,2019,40(16):196-200.