Evaluations and analyses of pollution factors of yogurt quality and safety in Lhasa
-
摘要:目的
全面分析评价影响西藏拉萨四县(区)(城关区、墨竹工卡县、林周县和当雄县)不同方式生产酸奶的质量安全关键污染因素，为提高拉萨酸奶质量安全水平、强化质量安全监管提供数据支撑。
方法2018年从拉萨四县(区)分别采集不同模式生产的酸奶样品106个，其中：林周县采样22个，城关区采样40个，墨竹工卡县采样23个，当雄县采样21个；按生产模式，企业采样8个，小作坊采样54个，家庭自产采样44个。分别使用电感耦合等离子体质谱法、HPLC-MS法、微生物培养与分离法对酸奶样品中5种重金属、黄曲霉毒素M1以及6种微生物进行检测分析与质量安全评价。
结果106个酸奶样品中霉菌污染严重，且超标率为46.23%；未检出志贺氏菌Shigella、沙门氏菌Salmonella和黄曲霉毒素；小作坊生产的酸奶中金黄色葡萄球菌Staphylococcus aureus检出率最高，检出率为25.93%；农牧民自身生产、民间小作坊和企业产业化生产酸奶中Pb是主要重金属污染物，污染物贡献率达42.5%；当雄县酸奶中重金属Cr含量较其他三县(区)高，超标率为19.0%；城关区酸奶中重金属Pb、As含量较高，超标率分别为15.0%和2.5%；Pb也是林周县、城关区、墨竹工卡县酸奶中的主要重金属污染物，其中城关区Pb的污染物贡献率最高，达54.02%。
结论拉萨地区酸奶重金属和霉菌污染状况突出，应当引起重视；西藏特殊的自然环境不利于饲料霉变，因此未检出黄曲霉毒素。
Abstract:ObjectiveTo evaluate and analyze the key pollution factors affecting the quality and safety of yogurt collected by different production methods in four different districts (Chengguan District, Maizhokunggar County, Linzhou County and Damxung County) of Lhasa in Tibet, and provide strong data support for improving the quality and safety of Lhasa yogurt and strengthening the quality and safety supervision.
MethodIn 2018, 106 yogurt samples with different production modes were collected from four counties (districts) of Lhasa, including 22 samples from Linzhou County, 40 samples from Chengguan District, 23 samples from Maizhokunggar County and 21 samples from Damxung County. Eight samples were collected from enterprises, 54 samples from individual workshops and 44 samples from families. Five kinds of heavy metals, aflatoxin M1 and six kinds of microorganisms in yogurt samples were detected and analyzed by ICP-MS, HPLC-MS, microbial culture and separation, respectively to evaluate the quality and safety.
ResultTotal 106 yogurt samples were seriously contaminated by mold, and the exceeding rate was 46.23%. Shigella, Salmonella and aflatoxin were not detected in yogurt. The detection rate of Staphylococcus aureus in yogurt produced by individual workshop was the highest at 25.93%. Pb was the main heavy metal pollutant in yogurt produced by farmers and herdsmen, individual workshops and enterprises, with the contamination contribution rate of 42.5%. The content of Cr in yogurt of Damxung County was higher than those in other three counties (districts), with the exceeding rate of 19.0%. The contents of Pb and As in yogurt of Chengguan District were higher, with the exceeding rates of 15.0% and 2.5%, respectively. Pb was also the main heavy metal pollutant in yogurt of Linzhou County, Chengguan District and Maizhokunggar County, among which the contamination contribution rate of Pb in Chengguan District was the highest reaching 54.02%.
ConclusionThe pollutions of heavy metals and molds in yogurt of Lhasa are prominent, which should be paid attention. The special natural environment in Tibet is not conducive to the mildew of feed, so no aflatoxin is detected.
-
Keywords:
- Lhasa /
- yogurt /
- toxic and harmful residue /
- heavy metal /
- aflatoxin M1 /
- microbial contamination
-
巴戟天Morinda officinalis别名巴戟、三蔓草、鸡肠风，是双子叶植物纲Dicotyledoneae、茜草目Rubiales、茜草科Rubiaceae多年生攀援木质藤本植物，是我国“四大南药”之一，以肉质根入药，其主要药效成分为寡糖、多糖蒽醌类、环烯醚萜类等活性物质[1-2]，具有补肝肾、强筋骨、祛风湿、抗衰老、抗肿瘤和治疗糖尿病、精神疾病以及细菌和病毒感染等作用[3-7]。巴戟天种植成本较高，且种植过程中存在开垦山林及乱采滥挖现象，这严重破坏了巴戟天的生长环境[8]。虽然巴戟天是第1批受到《广东省岭南中药材保护条例》保护的岭南药材，但其种植面积在逐年减少。土壤养分是影响巴戟天生长的重要环境因素之一，肉桂Cinnamomum cassia成林速度快，林内土壤动物、微生物种类更为丰富，可有效改善林地内土壤结构，从而提高土壤的肥力，营造出适宜巴戟天生长的种植环境[9]。
近年来，一些学者从巴戟天的药材化学成分及代谢组学[3, 10-12]、种植条件[13-14]、种植模式[15-16]等方面对巴戟天的有效成分含量进行比较研究，发现合理增加肉桂林种植地的光照量，可抑制藤苗生长，从而提高巴戟天肉质根的产量和有效成分积累，有效成分含量均随种植时间的延长呈上升趋势[16]。林下种植方式是控制种植地光照条件的重要途径之一，而且巴戟天林下种植具有较好的经济效益，解决了坡地连片种植带来的潜在山地环境破坏和水土流失问题，并且提高了巴戟天的产量和药效品质[8, 15]。这一林下生态种植方式为区域中药材的产业化发展提供了新思路。
广东省肇庆市德庆县是巴戟天的主产区，种植巴戟天有着悠久历史，是目前我国种植面积最大、产量最多和品质最好的巴戟天生产基地。本文研究了道地产区土壤养分和生长年限与巴戟天品质的关系，以及林下和非林下种植方式对巴戟天有效成分的影响，以期为巴戟天林下规范化种植提供基础数据和理论支持。
1. 材料与方法
1.1 研究区概况
本研究试验地位于广东省肇庆市德庆县，调查采样点分布于N 23°13′30″~N 23°17′37″，E 111°53′52″~E111°55′45″，海拔62~166 m。地处亚热带地区，热量丰富、雨热同期、夏长冬短，年平均温度20~25 ℃，年降雨量600~2 000 mm，土壤以红壤为主，其次是黄壤。德庆县当地农户种植的巴戟天有2个品种，分别是小叶品种‘黑幼籽’与大叶品种‘广宁特’，以小叶品种‘黑幼籽’最为常见。本次研究调查的品种均为‘黑幼籽’，种植密度普遍为株行距20 cm×30 cm，农家肥与大田复合肥为巴戟天施肥最常用肥料。大田复合肥的施肥量较低，通常为224~230 kg·hm−2，施肥频率为每年2~3次。
1.2 样品采集
于2018年10月，选取广东省德庆县高良镇以及邻近镇作为研究区域，根据该区域巴戟天种植面积和分布，选取18个具有代表性的非林下巴戟天种植点(表1)，调查分析土壤养分和种植年限对巴戟天寡糖含量的影响。另外，为了调查分析林下和非林下种植巴戟天的品质差异，在其中8个非林下巴戟天种植点的相邻林地上(相关信息见表2)，采集了林下种植点样品，巴戟天品种和种植年限以及土壤、坡度、坡向、海拔等均与对应的非林下种植点一致。在选定的地点，采用随机多点(5~8个点)采样法，在根系外围采集0~20 cm深的土壤，将土壤混匀并用四分法采集分析样品。土壤样品自然风干后，去除植物根系、落叶及石块等杂质，研磨，分别过1.00和0.25 mm筛，保存备用。同时在采样点随机整株挖取5株巴戟天。巴戟天肉质根用自来水清洗3次、去离子水润洗1次后，擦干。去除肉质根的木心部分，烘干，粉碎，过0.25 mm筛，待用。
表 1 道地产区非林下种植巴戟天采样点基本信息Table 1. Basic information of sample plot for no understory planting of Morinda officinalis in genuine producing areas采样点编号
Sample plot number经纬度
Longitude and latitude村或镇
Village or town生长年限
Growth year1 111°58′3″ E，23°18′25″ N 回村 Hui Village 1 2 112°2′12″ E，23°20′49″ N 莫村镇 Mocun Town 5 3 112°7′23″ E，23°20′51″ N 南田村 Nantian Village 4 4 112°2′20″ E，23°20′46″ N 莫村镇 Mocun Town 4 5 112°2′24″ E，23°20′47″ N 莫村镇 Mocun Town 2 6 112°2′22″ E，23°20′39″ N 莫村镇 Mocun Town 2 7 112°2′36″ E，23°20′44″ N 莫村镇 Mocun Town 5 8 112°5′49″ E，23°20′50″ N 荔枝村 Litchi Village 3 9 112°2′36″ E，23°20′53″ N 莫村镇 Mocun Town 5 10 112°2′39″ E，23°20′59″ N 莫村镇 Mocun Town 6 11 112°2′40″ E，23°21′0″ N 莫村镇 Mocun Town 2 12 112°2′54″ E，23°21′2″ N 莫村镇 Mocun Town 3 13 112°3′5″ E，23°21′23″N 莫村镇 Mocun Town 3 14 111°55′21″ E，23°15′55″ N 高良镇 Gaoliang Town 5 15 111°55′45″E，23°14′18″ N 榃浪村 Tanlang Village 5 16 111°55′11″ E，23°13′36″ N 中雄村 Zhongxiong Village 4 17 111°55′31″ E，23°13′42″ N 中雄村 Zhongxiong Village 4 18 111°54′54″ E，23°15′0″ N 年宅村 Nianzhai Village 5 表 2 林下种植巴戟天采样点基本信息Table 2. Basic information of sample plot for understory planting of Morinda officinalis采样点编号
Sample plot
number经纬度
Longitude
and latitude村或镇
Village
or town伴生植物
Associated plant透光度/%
Transparency生长年限
Growth
year1 111°55′45″E,
23°17′37″N降底村
Jiangdi Village杉木 Cunninghamia lanceolata、
马尾松 Pinus massoniana57.59 2 2 111°54′8″E,
23°17′34″N大揽村
Dalan Village火炬松 Pinus taeda、
山菅兰 Dianella ensifolia37.91 2 3 111°53′52″E,
23°15′29″N高良镇
Gaoliang Town肉桂 Cinnamomum cassia、
青皮竹 Bambusa textilis25.34 3 4 111°55′45″E,
23°14′18″N榃浪村
Tanlang Village肉桂 C. cassia、
马尾松 P. massoniana、
山苍子 Litsea cubeba35.89 5 5 111°55′31″E,
23°14′11″N榃浪村
Tanlang Village肉桂 C. cassia、火炬松 P. taeda、
五指毛桃 Ficus simplicissima36.15 5 6 111°55′11″E,
23°13′36″N中雄村
Zhongxiong Village杉木 C. lanceolata、
乌毛蕨 Blechnum orientale33.39 4 7 111°55′10″E,
23°13′30″N中雄村
Zhongxiong Village橄榄 Canarium album、
山麻黄 Psilopeganum sinense、
五指毛桃 F. simplicissima79.07 4 8 111°54′54″E,
23°15′0″N年宅村
Nianzhai Village阴香 Cinnamomum burmanni、
肉桂 C. cassia、荔枝 Litchi chinensis、
乌毛蕨 B. orientale61.38 5 1.3 样品分析
1.3.1 土壤养分含量分析
土壤样品的养分含量分析按照《土壤农化分析》[17]进行，具体为：称取10 g过1 mm筛的土壤，加入25 mL无二氧化碳去离子水浸提，利用pH计测定土壤pH；称取1 g过0.25 mm筛的土壤，采用凯氏定氮法测定土壤全氮含量；称取2 g过1 mm筛的土壤，采用碱解扩散法测定土壤碱解氮含量；称取1 g过0.25 mm筛的土壤，采用高氯酸−氢氟酸消煮−钼锑抗比色法测定土壤全磷含量；称取2 g过1 mm筛的土壤，采用氟化铵−盐酸浸提−钼锑抗比色法测定速效磷含量；称取0.25 g过0.25 mm筛的土壤，采用高氯酸−氢氟酸消煮−火焰光度法测定土壤全钾含量；称取5 g过1 mm筛的土壤，采用乙酸铵浸提–火焰光度法测定速效钾含量；称取0.5 g过0.25 mm筛的土壤，采用重铬酸钾氧化外加热法测定土壤有机质含量。
1.3.2 巴戟天肉质根品质成分分析
巴戟天肉质根全氮含量采用H2SO4−H2O2消煮−奈氏试剂比色法测定；全磷含量采用 H2SO4–H2O2 消煮−钼锑抗比色法测定；全钾采用H2SO4−H2O2消煮−原子吸收法测定。肉质根寡糖含量参照李倩[18]采用高效液相色谱–蒸发光散射法测定，通过峰面积计算样品中蔗糖、1−蔗果三糖、耐斯糖以及1F−果呋喃糖基耐斯糖的含量。
1.4 数据分析方法
试验数据采用Excel 2017和SPSS 25.0软件进行整理和统计分析，采用单因素方差分析(One-way ANOVA)，用Duncan’s多重极差检验法比较巴戟天品质指标的差异性；土壤养分与巴戟天品质指标的相关性则采用Pearson法进行分析；林下和非林下品质指标之间的差异性采用配对样本t检验法进行比较。
在采用主成分分析法进行综合评价时，选取反映巴戟天品质综合评估值的n个指标为X1，X2，…，Xn，m 个生长年限和种植环境的n项指标构成原始数据矩阵X=[Xij]n×m，其中，Xij为第i个模式的第j项指标数据(i=1, 2, …，n；j=1, 2, …，m)，按公式(1)对原始数据进行标准化，从而消除量纲和数量级的影响，然后确定主成分，从总方差分析表选取累积贡献率≥ 85%的前p个主成分，按公式(2)在建立p个主成分和标准化变量的关系后，按公式(3)确定权重，用第k个主成分的贡献率与选取的p个主成分的总贡献率的比值确定每个主成分的权重，最后按公式(4)构造综合评价函数，根据前p个主成分及其权重构造综合评价函数[19]。
$${X_{ij}}^* = ({X_{ij}} - \overline {{X_j}} )/{S_j},$$ (1) 式中，Xij*是 Xij的标准化数据，
$\overline {{X_j}} $ 和Sj是第j个指标的平均值和标准差。$${Y_k} = {u_{{k_1}}}{\boldsymbol{X}}_1^* + {u_{{k_2}}}{\boldsymbol{X}}_2^* + \cdots + {u_{{k_m}}}{\boldsymbol{X}}_m^*,$$ (2) 式中，Yk是第k个主成分(k=1, 2, ···，p)，uk1~ukm分别为第k1~km个主成分的因子荷载，X1*~Xm*均为标准化后的样本矩阵。
$${\omega _k} = {\lambda _k}/\sum\limits_{k = 1}^p {{\lambda _k}} ,$$ (3) 式中，ωk为第k个主成分的权重，λk为第k个主成分的贡献率。
$$F = \sum\limits_{k = 1}^p {{\omega _k}{Y_k}} ,$$ (4) 式中，F为不同生长年限和种植环境的巴戟天品质综合评价得分，分值越高表明该条件下的巴戟天品质越好。
2. 结果与分析
2.1 道地产区土壤养分含量
广东省德庆县是巴戟天道地产区。由表3可知，所收集到的土壤样品pH在4.08~4.90之间，差异较小，说明巴戟天适合种植在酸性土壤中。土壤有机质质量分数大部分在26~46 g·kg−1之间，11号样品有机质质量分数达到75.15 g·kg−1。土壤全氮、全磷、全钾的含量差异较大，这可能与当地的自然环境、施肥情况密切相关。
表 3 巴戟天道地产区土壤pH和养分含量Table 3. Soil pH and nutrient content in genuine producing areas of Morinda officinalis样品序号
Sample
number土壤类型
Soil typepH
w/(g·kg−1) w/(mg·kg−1) 全氮
Total N全磷
Total P全钾
Total K有机质
Organic
matter碱解氮
Available N速效磷
Available P速效钾
Available K1 红壤 Red soil 4.45 1.55 0.43 14.60 26.95 153.60 5.80 123.55 2 红壤 Red soil 4.25 1.91 0.62 11.18 27.72 157.52 46.26 154.32 3 红壤 Red soil 4.33 2.53 0.22 8.58 33.86 197.31 3.04 66.55 4 红壤 Red soil 4.50 2.35 0.25 9.13 43.09 188.70 4.67 68.00 5 红壤 Red soil 4.18 1.66 0.41 12.47 40.50 169.55 27.29 99.43 6 红壤 Red soill 4.10 2.00 0.20 12.89 30.25 199.93 3.62 71.67 7 红壤 Red soil 4.52 2.11 0.46 18.21 26.82 168.97 12.47 89.65 8 红壤 Red soil 4.51 2.06 0.26 10.05 37.04 202.02 11.27 131.75 9 红壤 Red soil 4.40 2.09 0.25 9.99 37.12 175.84 6.47 142.88 10 红壤 Red soil 4.25 1.40 0.21 17.32 36.23 342.07 3.45 74.70 11 红壤 Red soil 4.18 4.08 0.33 11.16 75.15 230.14 12.27 90.12 12 黄壤 Yellow soil 4.12 2.41 0.47 14.63 46.04 225.33 29.39 186.59 13 黄壤 Yellow soil 4.25 2.07 0.32 12.33 46.31 192.76 24.72 140.38 14 红壤 Red soil 4.08 1.06 0.27 6.78 37.05 162.21 14.09 90.78 15 红壤 Red soil 4.56 0.88 0.49 9.65 23.64 78.15 17.85 123.57 16 红壤 Red soil 4.66 2.70 0.66 7.17 33.68 118.69 38.14 146.22 17 红壤 Red soil 4.69 1.31 0.76 6.75 32.77 118.48 41.00 151.24 18 红壤 Red soil 4.68 1.47 0.47 4.99 28.03 136.20 32.77 136.34 平均值 Mean 4.37 1.98 0.39 10.99 36.79 178.75 18.59 115.99 最小值 Minimum 4.08 0.88 0.20 4.99 23.64 78.15 3.04 66.55 最大值 Maximum 4.69 4.08 0.76 18.21 75.15 342.07 46.26 186.59 标准误差 Standard error 0.01 0.04 0.01 0.20 0.65 3.13 0.79 1.97 变异系数/% Coefficient of variation 4.71 36.73 42.21 33.12 31.66 31.48 76.44 30.59 2.2 巴戟天品质成分含量
巴戟天药材中全氮质量分数为2.81~12.30 g·kg−1，全磷质量分数为0.12~0.71 g·kg−1，全钾质量分数为6.13~13.20 g·kg−1；巴戟天药材中4种寡糖成分含量差异较大，蔗糖质量分数为34.06~76.16 g·kg−1，1−蔗果三糖质量分数为9.99~59.30 g·kg−1，耐斯糖质量分数为27.93~69.41 g·kg−1，符合2015年版《中华人民共和国药典》巴戟天干燥品耐斯糖含量不得少于2.0%的要求[20]；1F−果呋喃糖基耐斯糖质量分数为46.60~97.12 g·kg−1(表4)。
表 4 道地产区巴戟天肉质根全氮、全磷、全钾和寡糖含量Table 4. Total N, total P, total K and oligosaccharide contents of succulent root of Morindus officinalis in genuine producing areasw/(g·kg−1) 样品序号
Sample
number全氮
Total N全磷
Total P全钾
Total K蔗糖
Sucrose1−蔗果三糖
1-Kestose耐斯糖
Nystose1F−果呋喃糖基耐斯糖
1F- Fructo-
furanosaccharides1 3.97 0.26 6.47 54.35 24.62 52.62 80.44 2 9.22 0.60 11.61 57.72 29.51 43.25 57.77 3 7.32 0.33 9.16 50.73 12.58 27.93 46.60 4 4.30 0.31 6.36 50.04 16.65 48.37 74.99 5 3.97 0.23 11.82 52.18 14.44 51.72 85.78 6 4.75 0.26 9.69 57.43 21.47 49.55 75.85 7 5.93 0.38 13.20 46.25 10.64 29.77 47.46 8 4.09 0.22 8.28 59.65 21.73 51.75 81.96 9 5.96 0.56 10.80 57.72 19.69 42.22 67.56 10 7.90 0.34 10.64 65.43 28.93 37.91 55.72 11 6.30 0.42 8.26 51.73 18.17 42.33 69.34 12 12.30 0.68 13.19 76.16 34.97 43.91 58.06 13 9.84 0.71 11.65 56.42 21.83 46.15 72.50 14 2.81 0.22 6.47 46.33 9.99 31.85 53.77 15 8.77 0.47 6.13 34.06 40.60 62.00 95.48 16 11.70 0.61 9.31 46.49 46.79 65.46 97.12 17 12.33 0.54 9.05 45.33 59.30 69.41 95.10 18 11.99 0.12 6.58 34.09 42.18 58.77 87.70 平均值 Mean 7.41 0.40 9.37 52.34 26.34 47.50 72.40 最小值 Minimum 2.81 0.12 6.13 34.06 9.99 27.93 46.60 最大值 Maximum 12.33 0.71 13.20 76.16 59.30 69.41 97.12 标准误差 Standard error 0.18 0.01 0.13 0.56 0.76 0.65 0.92 变异系数/% Coefficient of variation 43.51 44.07 25.36 19.27 51.83 24.64 22.75 2.3 巴戟天品质与土壤养分含量的相关性分析
从相关性分析结果(表5)可以看出，土壤养分含量与巴戟天营养元素含量密切相关，其中，土壤全磷、速效磷、速效钾含量与巴戟天全氮含量均表现出极显著正相关关系，土壤速效钾含量与巴戟天全磷含量之间以及土壤全钾含量与巴戟天全钾含量之间均具有极显著正相关关系。土壤养分含量与巴戟天寡糖含量密切相关，其中，土壤碱解氮含量与巴戟天蔗糖含量呈极显著正相关关系，但却与耐斯糖和1F−果呋喃糖基耐斯糖含量呈显著负相关关系。土壤全磷、速效磷、速效钾含量与巴戟天1−蔗果三糖含量均表现出极显著正相关关系，与耐斯糖含量表现出极显著或显著正相关关系。值得注意的是，土壤pH与巴戟天营养元素没有显著相关关系，但与巴戟天寡糖呈显著相关关系，其中，土壤pH与蔗糖含量表现出极显著的负相关关系，与1−蔗果三糖、耐斯糖以及1F−果呋喃糖基耐斯糖含量均有显著或极显著的正相关关系，说明随着土壤pH的增加蔗糖含量降低，而1−蔗果三糖、耐斯糖与1F−果呋喃糖基耐斯糖含量却呈增加趋势。另外，土壤有机质含量与巴戟天品质指标之间相关关系不显著，说明土壤有机质含量对巴戟天品质影响不大。
表 5 巴戟天品质与土壤养分含量的相关性分析1)Table 5. Correlation analysis between Morinda officinalis quality and soil nutrient contents指标
Index全氮
Total N全磷
Total P全钾
Total K蔗糖
Sucrose1−蔗果三糖
1-Kestose耐斯糖
Nystose1F−果呋喃糖基耐斯糖
1F- FructofuranosaccharidespH 0.373 −0.027 −0.349 −0.595** 0.578* 0.596** 0.570* 全氮 Total N 0.012 0.246 0.213 0.279 −0.246 −0.223 −0.193 碱解氮 Available N −0.192 −0.079 0.345 0.712** −0.384 −0.554* −0.570* 全磷 Total P 0.651** 0.424 0.105 −0.303 0.749** 0.609** 0.455 速效磷 Available P 0.700** 0.455 0.241 −0.156 0.650** 0.513* 0.366 全钾 Total K −0.228 0.111 0.593** 0.570* −0.363 −0.431 −0.458 速效钾 Available K 0.645** 0.610** 0.254 0.178 0.614** 0.471* 0.324 有机质 Organic matter −0.079 0.179 0.086 0.313 −0.258 −0.193 −0.131 1)“*”和“**”分别表示达0.05和0.01水平的显著相关(Pearson法)
1) “*” and “**” indicate significant correlations at 0.05 and 0.01 levels, respectively (Pearson method)2.4 不同生长年限巴戟天的品质指标含量
本次调查采集的是非林下种植环境的18个样地，巴戟天的生长年限包含1至6年生。由于1年生和6年生的仅有1个样点，故仅分析2~5年生巴戟天肉质根的营养元素和寡糖含量。
2.4.1 不同生长年限巴戟天药材中的营养元素含量
由图1可见，随着巴戟天生长年限的增长，其药用部分肉质根对全氮、全磷、全钾吸收积累量呈现先升高后降低的趋势，整体表现为全钾>全氮>全磷。2年生巴戟天的全氮含量显著低于其他生长年限的，其他处理间的全氮含量没有显著性差异。全磷含量、全钾含量在不同生长年限之间无显著差异，但均在3年生时达到最大值，分别为0.54、11.04 g·kg−1。
![]()
图 1 不同生长年限巴戟天药材营养元素含量各图中，柱子上方的不同小写字母表示生长年限间差异显著(P<0.05，Duncan’s法)Figure 1. Nutrient element contents in Morinda officinalis of different growth yearsIn each figure, different lowercase letters on the columns indicate significant differences among different growth years (P<0.05, Duncan’s method)2.4.2 不同生长年限巴戟天药材中的寡糖含量
由图2可以看出，巴戟天生长年限不同，寡糖类成分含量各异。随着巴戟天生长年限的增加，蔗糖、1−蔗果三糖含量变化均呈现出先增加后下降的趋势。其中，蔗糖含量为3年生巴戟天最高，质量分数达到64.08 g·kg−1，显著高于其他生长年限的，同时，2年生巴戟天蔗糖含量显著高于4年生和5年生的，而4年生与5年生巴戟天的蔗糖含量之间无显著差异。在1-蔗果三糖含量上，4年生巴戟天显著高于2年生巴戟天(高出76.96%)，4年生巴戟天达到32.68 g·kg−1，与3年生、5年生巴戟天没有显著差异，2年生、3年生与5年生巴戟天之间无显著差异。不同生长年限的巴戟天在耐斯糖和1F−果呋喃糖基耐斯糖含量上没有显著差异。
![]()
图 2 不同生长年限巴戟天寡糖含量各图中，柱子上方的不同小写字母表示生长年限间差异显著(P<0.05，Duncan’s法)Figure 2. Oligosaccharide contents in Morinda officinalisof different growth yearsIn each figure, different lowercase letters on the columns indicate significant differences among different growth years (P<0.05, Duncan’s method)2.5 林下与非林下种植巴戟天品质及其主成分分析
2.5.1 林下与非林下种植巴戟天的营养元素和寡糖含量
由于生产上巴戟天收获年龄为4年生，故本研究选择4年生的林下与非林下种植的巴戟天进行品质分析。林下与非林下种植的巴戟天全磷含量与全钾含量差异显著(图3)，林下种植的巴戟天全磷含量比非林下种植巴戟天高44.79%，而全钾含量比非林下巴戟天低13.78%，全氮含量无显著差异。
![]()
图 3 不同种植方式对4年生巴戟天营养元素含量的影响各图中，柱子上方的不同小写字母表示种植方式间差异显著(P<0.05，Duncan’s法)Figure 3. Effects of different planting methods on nutrient element contents of 4-year-old Morinda officinalisIn each figure, different lowercase letters on the columns indicate significant differences between different planting methods (P<0.05, Duncan’s method)由图4可以看出，林下与非林下种植的巴戟天除耐斯糖含量没有显著差异外，其他寡糖含量间均有显著差异。其中，林下种植的4年生巴戟天蔗糖和1−蔗果三糖质量分数分别达到53.14和60.22 g·kg−1，分别高出非林下种植巴戟天32.49%和18.68%。而1F−果呋喃糖基耐斯糖含量却低于非林下巴戟天3.97%。
![]()
图 4 不同种植方式对4年生巴戟天寡糖含量的影响各图中，柱子上方的不同小写字母表示种植方式间差异显著(P<0.05，Duncan’s法)Figure 4. Effects of different planting methods on oligosaccharide content of 4-year-old Morinda officinalisIn each figure, different lowercase letters on the columns indicate significant differences between different planting methods (P<0.05, Duncan’s method)2.5.2 林下与非林下种植巴戟天品质的主成分分析
主成分分析是对多个变量进行统计分析的一类定量分析方法，通过分析主成分载荷值和得分可以看出不同样品之间的关系，也可以解释样品的特点、分组、相似性及差别[21-22]。将巴戟天品质指标转化为7个主成分，分析结果如表6。7种指标可以提取出2种主成分，累计贡献率为85.776%。
表 6 巴戟天品质指标的主成分初始特征值及累计贡献率Table 6. Initial characteristic values and cumulative contribution rates of principal components of Morinda officinalis quality index变量1)
Variable初始特征值 Initial eigenvalue 提取载荷平方和 Squared sum of extracted load 特征值
Eigen
value贡献率/%
Contribution
rate累计贡献率/%
Cumulative
contribution rate特征值
Eigen
value贡献率/%
Contribution
rate累计贡献率/%
Cumulative
contribution rateX1 3.264 46.632 46.632 3.264 46.632 46.632 X2 2.740 39.144 85.776 2.740 39.144 85.776 X3 0.776 11.079 96.855 X4 0.166 2.370 99.225 X5 0.043 0.609 99.834 X6 0.011 0.159 99.993 X7 0.000 0.007 100.000 1)X1：全氮；X2：全磷；X3：全钾；X4：蔗糖；X5：1–蔗果三糖；X6：耐斯糖；X7：1F–果呋喃糖基耐斯糖
1)X1: Total N; X2: Total P; X3: Total K; X4: Sucrose; X5: 1-Kestose; X6: Nystose; X7: 1F- Fructofuranosaccharides为了更好地解释各品质指标与成分因子之间的关系，将所提取的主成分因子进行旋转处理，使一个变量在较少的几个因子上有较高的载荷，其载荷值的大小反映各有效成分指标在主成分中的重要程度(表7)。通过主成分分析可知，前2个主成分反映了7个指标的大部分综合信息，用2个主成分评价不同种植方式下巴戟天品质的优劣是可行的。因此，可用营养因子(F1)、寡糖因子(F2)代替原来7个指标(X1~X7)对巴戟天的品质进行分析，得到巴戟天品质前2个主成分的线性关系式分别为：
表 7 巴戟天品质指标前2个主成分的载荷值Table 7. Load values of the first two principal components of Morinda officinalis quality index评价指标
Evaluation index载荷值 Loading value 第一主成分
The first principal component第二主成分
The second principal component全氮 Total N 0.41 −0.15 全磷 Total P 0.52 0.17 全钾 Total K −0.42 −0.11 蔗糖 Sucrose 0.36 0.45 1−蔗果三糖 1-Kestose 0.10 0.59 耐斯糖 Nystose −0.28 0.51 1F−果呋喃糖基耐斯糖 1F- Fructofuranosaccharides −0.42 0.36 F1=0.41X1+0.52X2−0.42X3+0.36X4+0.10X5−0.28X6−0.42X7，
F2= −0.15 X1+0.17X2−0.11X3+0.45X4+0.59X5+0.51X6+0.36X7。
以每个主成分对应的特征值的方差提取贡献率建立综合评价模型：
F=0.15X1+0.36X2−0.28X3+0.40X4+0.32X5+0.08X6−0.06X7。
计算林下与非林下种植方式下巴戟天品质的综合评分(表8)，得到不同种植方式下巴戟天品质特征的优劣顺序：7号样地林下种植巴戟天>6号样地林下种植巴戟天>6号样地非林下种植巴戟天>7号样地非林下种植巴戟天。说明，林下种植的巴戟天品质优于同条件下非林下种植的巴戟天。
表 8 4年生巴戟天品质指标前2个主成分得分与排名Table 8. Scores and rankings of the first two principal components in the quality index of 4-year-old Morinda officinalis采样点编号
Sample
plot number种植方式
Planting
method第一主成分
The first principal component第二主成分
The second principal component综合得分
Synthesis
score排名
Rank得分 Score 排名 Rank 得分 Score 排名 Rank 6 非林 No understory −1.58 4 1.09 2 −0.36 3 林下 Understory −0.33 2 1.43 1 0.48 2 7 非林下 No understory −0.87 3 −2.58 4 −1.65 4 林下 Understory 2.78 1 0.07 3 1.54 1 3. 讨论与结论
3.1 讨论
土壤理化性质是反映土壤养分的重要指标，土壤养分是植物生长发育的主要营养来源，同时也是创造优质高产必备的物质条件之一，研究其对巴戟天品质的影响有非常积极的意义。本研究表明，巴戟天品质与土壤养分关系密切，这与刘瑾等[23-24]的研究结果一致。药理研究显示，巴戟天寡糖具有抗抑郁功效，其主要发挥功能的糖分为四~七糖[25]，如耐斯糖(简称：四糖)、1−果呋喃糖基耐斯糖(简称：五糖)，土壤碱解氮与巴戟天蔗糖含量呈极显著正相关关系，但与耐斯糖和1F−果呋喃糖基耐斯糖呈显著负相关关系，故减少土壤碱解氮的含量有利于巴戟天抗抑郁成分的积累。而土壤全磷、速效磷、速效钾对巴戟天1−蔗果三糖、耐斯糖均表现出显著的正相关关系，表明提高土壤磷和钾的有效性，有助于促进巴戟天寡糖累积。土壤的酸碱性对巴戟天品质的影响较大，巴戟天更适合种植于弱酸性土壤中[26]，与研究发现弱酸性的红壤或棕壤较适宜巴戟天的生长一致[15]。土壤pH与巴戟天根部营养元素含量无显著相关性，与蔗糖含量呈现出极强的负相关性，而与另外其他3种寡糖之间具有正相关关系，说明当土壤呈弱酸性时，会降低巴戟天根部的蔗糖含量，但会提高1−蔗果三糖、耐斯糖与1F−果呋喃糖基耐斯糖含量。土壤有机质是土壤肥力的重要组成部分，能够使植株正常生长发育，保证植物在生长过程中的养分需求[27]，但与巴戟天品质指标没有显著相关性。
不同生长年限的林下种植巴戟天的1−蔗果三糖和耐斯糖含量均大于同年生的非林下巴戟天，说明林下巴戟天从土壤中获取的养分较非林下巴戟天的多，从而使更多的单糖转化为1−蔗果三糖与耐斯糖。随着巴戟天生长年限的增加，巴戟天营养元素及寡糖含量呈现先增加后减少的变化趋势。这与玉竹Polygonatum odoratum、川党参Codonopsis pilosula[28]等药用植物有效成分积累的趋势一致，这主要是因为巴戟天生长前期(第1~4年)通过植物光合作用不断贮存营养物质，而进入衰老期(第5年)后营养物质的消耗量大于累积量，导致有效成分含量下降。此外，巴戟天药材中发挥抗痴呆、抗抑郁、保护生殖系统等功效的成分主要是低聚糖[29-31]，而低聚糖需要由单糖在一系列酶的催化下才能逐渐被合成，只有经过种植年限的积累，低聚糖含量才会升高，且巴戟天药材种植4年后其有效药用成分均能满足《中华人民共和国药典》[20]巴戟天干燥品耐斯糖含量不得少于2.0%的要求，故4年生巴戟天获得的药材质量最佳。
通过对不同种植方式下同年生巴戟天品质特征进行分析，林下种植的4年生巴戟天全磷、蔗糖、1−蔗果三糖含量显著高于非林下种植的，全钾、1F−果呋喃糖基耐斯糖含量却显著低于非林下种植的，可能是因为植物营养元素与有效成分并不是独立发挥作用，各成分之间具有一定相关性[32]，从而导致不同种植方式对巴戟天品质特征物质含量有不同的促进作用。树林的郁闭度较大，为林下巴戟天提供了很好的遮荫效果，起到降温保湿的作用[33]，这是林下种植的巴戟天有效成分含量优于非林下种植巴戟天的原因之一。
3.2 结论
土壤养分对巴戟天肉质根养分和寡糖含量有不同的影响，其中土壤pH与蔗糖含量表现出极显著的负相关性，与 1−蔗果三糖、耐斯糖以及1F−果呋喃糖基耐斯糖含量呈显著正相关性。土壤全磷、速效磷、速效钾含量与1−蔗果三糖和耐斯糖含量呈显著正相关关系，表明在巴戟天生长过程中可以通过增加土壤速效磷与速效钾含量以提高其品质。不同生长年限的巴戟天肉质根全氮、全磷、全钾含量没有显著差异，而3年生巴戟天蔗糖含量显著高于其他生长年限的巴戟天，4年生巴戟天的1−蔗果三糖含量显著高于2年生巴戟天。林下与非林下种植的巴戟天肉质根全氮、全磷、全钾、蔗糖、1−蔗果三糖、耐斯糖和1F−果呋喃糖基耐斯糖含量无一致规律，但主成分分析法综合评分结果显示，林下种植巴戟天品质优于非林下种植的巴戟天。
-
![]()
图 1 西藏4县(区)酸奶采样点分布示意图
Figure 1. Distribution diagram of yogurt sampling points in four counties (districts) of Tibet
表 1 乳制品质量分级标准
Table 1 Quality grading standards for dairy products
等级
Classification单因子污染指数
Single factor
pollution index综合污染指数
Composite
pollution index污染状况
Pollution status质量水平
Quality level1 ≤0.6 ≤0.6 有污染物残留，污染物含量接近或略高于背景
Residual pollutant was detected, and the content was close
to or slightly higher than background安全
Safety2 0.6~1.0 0.6~1.0 产品污染物含量较高，限量值范围内
The content of pollutant in product was high and
within the limit value range轻度污染
Mild contamination3 ≥1.0 ≥1.0 污染产品，污染物含量超过限量
The pollutant content in polluted product
exceeded the limit value重度污染
Severe contamination
下载: 导出CSV
表 2 不同区域生产酸奶中重金属检出情况
Table 2 Detection status of heavy metals in yogurt produced in different regions
区域
District检出率/%
Detection rate超标率/%
Exceeding standard rate污染指数及综合污染指数
Pollution index and
composite pollution index污染物贡献率/%
Pollutant contribution rateCd Cr Pb Hg As Cr Pb Hg As PCr PPb PAs PHg Pn KCr KPb KAs KHg 林周县
Linzhou County36.4 27.3 22.7 0 45.5 9.1 9.1 0 0.0 0.220 0.290 0.098 0.070 0.237 32.50 42.73 14.44 10.33 城关区
Chengguan
District30.0 27.5 22.5 2.5 47.5 7.5 15.0 0 2.5 0.213 0.541 0.177 0.070 0.421 21.32 54.02 17.67 6.99 墨竹工卡县
Maizhokunggar
County30.4 26.1 26.1 0 47.8 13.0 8.7 0 0 0.186 0.240 0.108 0.070 0.200 30.82 39.70 17.89 11.59 当雄县
Damxung County38.1 23.8 23.8 0 61.9 19.0 9.5 0 0 0.688 0.345 0.185 0.070 0.537 53.44 26.78 14.34 5.43 合计 Total 33.0 26.4 23.6 0.9 50.0 11.3 11.3 0 0.9 0.303 0.385 0.147 0.070 0.316 33.49 42.51 16.26 7.74
下载: 导出CSV
表 3 不同模式生产酸奶中重金属检出情况
Table 3 Detection status of heavy metals in yogurt produced by different modes
生产模式
Production modew/(mg·kg−1) 检出率/%
Detection rate超标率/%
Exceeding standard rate污染指数及综合污染指数
Pollution index and composite pollution index污染物贡献率/%
Pollutant contribution rateCd Cr Pb As Cd Cr Pb Hg As Cr Pb Hg As PCr PPb PAs PHg Pn KCr KPb KAs KHg 家庭自产Family production 0.008±0.007 0.248±0.246 0.052±0.053 0.022±0.018 31.8 25.0 25.0 0 45.5 11.4 9.1 0 0 0.207 0.270 0.101 0.070 0.223 31.99 41.63 15.61 10.77 企业Enterprise 0.004±0.005 0.017 0.031 0.016±0.013 25.0 12.5 12.5 0 37.5 0 0 0 0 0.008 0.088 0.066 0.070 0.074 3.49 37.78 28.64 30.09 小作坊Individual workshop 0.072±0.288 0.429±0.548 0.107±0.107 0.035±0.034 35.2 29.6 24.1 1.9 55.6 13.0 14.8 0 1.9 0.425 0.522 0.197 0.070 0.427 35.00 43.02 16.23 5.75 合计 Total 33.0 26.4 23.6 0.9 50.0 11.3 11.3 0 0.9 0.303 0.385 0.147 0.070 0.316 33.50 42.50 16.29 7.71
下载: 导出CSV
表 4 不同区域生产酸奶中重金属含量单因素方差分析
Table 4 One-way analysis of variance of heavy metal content in yogurt produced in different regions
重金属种类
Type of heavy metal源
OriginIII型平方和
Sum of squares of type IIIdf 均方
Mean squareF P Cd 组间 Interblock 0.148 3 0.049 1.104 0.362 组内 Intra-group 1.385 31 0.045 总数 Sum 1.534 34 Cr 组间 Interblock 1.668 3 0.556 3.539 0.030 组内 Intra-group 3.770 24 0.157 总数 Sum 5.438 27 Pb 组间 Interblock 0.023 3 0.008 0.989 0.417 组内 Intra-group 0.162 21 0.008 总数 Sum 0.185 24 As 组间 Interblock 0.002 3 0.001 0.945 0.426 组内 Intra-group 0.040 49 0.001 总数 Sum 0.042 52
下载: 导出CSV
表 5 不同模式生产酸奶中重金属含量单因素方差分析
Table 5 One-way analysis of variance of heavy metal content in yogurt produced by different production modes
重金属种类
Type of heavy metal源
OriginIII型平方和
Sum of squares of type IIIdf 均方
Mean squareF P Cd 组间 Interblock 0.037 2 0.018 0.394 0.678 组内 Intra-group 1.497 32 0.047 总数 Sum 1.534 34 As 组间 Interblock 0.003 2 0.001 1.697 0.194 组内 Intra-group 0.040 50 0.001 总数 Sum 0.042 52
下载: 导出CSV
表 6 酸奶中重金属含量相关性分析1)
Table 6 Correlation analysis of heavy metal content in yogurt
重金属种类
Type of heavy metalCd Cr Pb As Cd 0.243 −0.080 −0.054 Cr −0.466 0.582** Pb 0.030 As 1)“**”表示相关性达到0.01的显著水平(Pearson法)
1) “**” indicated significant correlation at 0.01 level (Pearson method)
下载: 导出CSV
表 7 不同区域生产酸奶中微生物状况分析
Table 7 Analyses of microorganisms in yogurt produced in different regions
区域
District样品数/个
No. of sample菌落总数
Aerobic plate count大肠菌群
Coliform霉菌
Mold金黄色葡萄菌
Staphylococcus aureus超标数/个
No. of sample exceeding standard超标率/%
Exceeding standard rate超标数/个
No. of sample exceeding standard超标率/%
Exceeding standard rate检出范围/
(cfu·g−1)
Range of detected dosage超标数/个
No. of sample exceeding standard超标率/%
Exceeding standard rate检出数/个
No. of sample detected检出率/%
Detection rate林周县
Linzhou
County22 0 0 1 4.55 0~183 11 50.00 4 18 城关区
Chengguan
District40 1 2.5 2 5.00 0~283 20 50.00 12 56 墨竹工卡县
Maizhokunggar
County23 0 0 2 8.70 0~260 15 65.22 5 23 当雄县
Damxung
County21 1 4.8 0 0 0~90 3 14.29 4 19 合计Total 106 2 1.9 5 4.72 0~283 49 46.23
下载: 导出CSV
表 8 不同模式生产酸奶中微生物状况分析
Table 8 Analysis of microorganisms in yogurt with different production modes
生产模式
Production mode样品数/个
No. of sample大肠菌群
Coliform霉菌
Mold金黄色葡萄菌
Staphylococcus aureus超标数/个
No. of sample exceeding standard超标率/%
Exceeding standard
rate超标数/个
No. of sample exceeding standard超标率/%
Exceeding standard rate检出数/个
No. of sample detected检出率/%
Detection rate家庭自产 Family production 44 3 6.82 25 56.82 9 20.45 企业 Enterprise 8 0 0 4 50.00 2 25.00 小作坊 Individual workshop 54 2 3.70 20 37.04 14 25.93
下载: 导出CSV
-
[1] 扎木苏. 西藏地区传统牦牛奶制品中乳酸菌的多样性研究[D]. 呼和浩特: 内蒙古农业大学, 2014. [2] 包秋华. 甘肃和四川省牦牛奶制品中乳酸菌的多样性研究[D]. 呼和浩特: 内蒙古农业大学, 2012. [3] 陶亮, 张乃明. 云南部分地区农产品重金属污染现状与分析[J]. 中国农学通报, 2017, 33(20): 83-89. doi: 10.11924/j.issn.1000-6850.casb16070064 [4] KEITH L, TELLIARD W. ES&T Special Report: Priority pollutants: I-a perspective view[J]. Environ Sci Technol, 1979, 13(4): 416-423. doi: 10.1021/es60152a601
[5] FAO. Principles and methods for the risk assessment of chemicals in food[M]. Geneva: Switzerland World Health Organization, 2009.
[6] 徐琳. 鲜牛奶中黄曲霉毒素M1检测的几种方法及控制措施[J]. 现代畜牧科技, 2018(12): 99. [7] 吴坤坦, 齐德生. 黄曲霉毒素B1在乳畜中的转化[J]. 饲料工业, 2018, 39(20): 59-64. [8] 翟晓瑞, 李燕. 乳品中黄曲霉毒素M1残留检测方法研究进展[J]. 食品安全质量检测学报, 2016, 7(6): 2309-2313. [9] 丁松乔, 肖志刚. 乳制品中黄曲霉毒素M1的常用检测方法比较分析[J]. 中国乳品工业, 2016, 44(5): 36-39. doi: 10.3969/j.issn.1001-2230.2016.05.009 [10] FRAZZOLI C, GHERARDI P, SAXENA N, et al. The hotspot for (global) one health in primary food production: Aflatoxin M1 in dairy products[J/OL]. Front Public Health, 2017, 4: 1-11. [2020−01−03]. https://doi.org/10.3389/fpubh.2016.00294.
[11] 许尨, 张冬洁, 李洪亮, 等. 酸奶发展的研究进展[J]. 农产品加工, 2019(12): 87-89. [12] 刘文文, 张靖晞, 李键, 等. 不同菌株对酸奶品质与感官特性的影响[J]. 中国酿造, 2019, 38(5): 95-100. doi: 10.11882/j.issn.0254-5071.2019.05.019 [13] 中华人民共和国国家卫生和计划生育委员会, 国家食品药品监督管理总局. 食品安全国家标准食品微生物学检验: GB 4789—2016[S]. 北京: 中国标准出版社, 2017. [14] 中华人民共和国国家卫生和计划生育委员会, 国家食品药品监督管理总局. 食品安全国家标准食品中黄曲霉毒素M族的测定: GB 5009.24—2016[S]. 北京: 中国标准出版社, 2017: 1-19. [15] 陈国友, 杜英秋, 李宛, 等. 应用ICP-MS、AFS、GF-AAS测定食品中As、Cd、Hg、Pb方法的对比研究[J]. 质谱学报, 2009, 30(4): 223-228. [16] 张唐伟, 吴雪莲, 郝治华, 等. ICP-MS法测定酸奶中重金属元素[J]. 农业与技术, 2019, 39(17): 35-37. [17] 祖文川, 汪雨, 李冰宁, 等. ICP-MS相关联用技术在食品元素形态分析中的应用及进展[J]. 质谱学报, 2013, 34(4): 247-256. doi: 10.7538/zpxb.2013.34.04.0247 [18] 史潜玉, 刘立, 柯润辉, 等. ICP-MS在食品质量安全领域应用研究进展[J]. 食品与发酵工业, 2012, 38(12): 118-123. [19] 中华人民共和国国家卫生和计划生育委员会, 国家食品药品监督管理总局. 食品安全国家标准 食品中污染物限量: GB 2762—2017[S]. 北京: 中国标准出版社, 2017. [20] 中华人民共和国农业部. 农、畜、水产品污染监测技术规范: NY/T 398—2000[S]. 北京: 中国标准出版社, 2000: 67-85. [21] 王海, 沈秋光, 邹明晖, 等. 各国乳品的生乳标准分析比对[J]. 乳业科学与技术, 2011, 34(6): 293-295. doi: 10.3969/j.issn.1671-5187.2011.06.012 [22] 李青冉, 高腾云, 凡复, 等. 牛奶质量评定与计价体系[J]. 中国牛业科学, 2007, 33(2): 37-39. doi: 10.3969/j.issn.1001-9111.2007.02.011 [23] 赵卓, 于冷, 陆骏飞. 我国原奶按质论价的原理、影响及政策含义[J]. 中国畜牧杂志, 2009, 45(4): 41-44. [24] 于炎湖. 饲料安全性问题: 4: 饲料中重金属元素污染的来源、危害及其预防[J]. 养殖与饲料, 2003(2): 3-5. doi: 10.3969/j.issn.1671-427X.2003.02.002 [25] LI S L, MIN L, WANG P P, et al. Aflatoxin M1 contamination in raw milk from major milk-producing areas of China during four seasons of 2016[J]. Food Control, 2017, 82: 121-125. doi: 10.1016/j.foodcont.2017.06.036
[26] 孙思, 王安波, 杨梅, 等. 生鲜乳中黄曲霉毒素M1检测技术研究进展[J]. 畜牧兽医科技信息, 2017(4): 12-13. doi: 10.3969/J.ISSN.1671-6027.2017.04.006 [27] 中华人民共和国国家卫生和计划生育委员会, 国家食品药品监督管理总局. 食品中真菌毒素限量: GB 2761—2017 [S]. 北京: 中国标准出版社, 2017: 2-3. [28] 张书义. 乳业生产安全概论: 二[J]. 中国乳业, 2011(2): 52-56. doi: 10.3969/j.issn.1671-4393.2011.02.014 -
期刊类型引用(4)
1. 陈海艺，徐悦，刘艺璇，张亚敏，郭红莲. 膜醭毕赤酵母蛋白复合膜对桃采后果腐病的抑制效果及机制. 食品与发酵工业. 2024(05): 274-279 . 
百度学术
2. 余生玲，申光辉，黄廷，杨露，张志清，侯晓艳，陈安均，黎杉珊. 藤椒精油对腐败解淀粉芽孢杆菌DY1a生物被膜的抑制作用. 微生物学通报. 2022(06): 2135-2151 . 百度学术
3. 韩松林，安才让吉. 大蒜素生物学功能及其在畜禽养殖中的应用研究进展. 饲料研究. 2022(20): 157-159 . 百度学术
4. 肖小年，盛丹梅，周洁，易孜成，朱碧泉，张俊，易醒. 大蒜油自微乳的抑菌作用. 中国食品学报. 2021(12): 73-80 . 百度学术
其他类型引用(4)
计量
- 文章访问数:
- HTML全文浏览量: 0
- PDF下载量:
- 被引次数: 8
