中国花生黄曲霉毒素污染风险预警模型研究

发布时间：2020-08-01 17:03

【摘要】：花生(Arachis hypogaea L.)是我国主要油料作物和优势出口农产品,是食油兼用型经济作物,其营养价值高兼具食疗功能,在中国及世界油料生产中占有重要地位。然而花生易受真菌侵染产生黄曲霉毒素,严重影响人体健康,已成为影响我国花生消费安全和世界贸易的主要制约因素。花生黄曲霉毒素污染可以发生在种植、收获、储藏、运输、销售等各个环节,其中收获前田间污染过程既是花生黄曲霉毒素污染的主要原因,也是储藏期污染的源头。花生黄曲霉毒素污染防控是一个世界性难题,由于我国花生品种繁多,土壤类型和气候类型多样,因此了解花生黄曲霉毒素污染时空消长变化趋势,解析品种抗性、土壤类型、真菌产毒力及气候条件对其影响,开展花生黄曲霉毒素污染预警模型研究,构建农产品真菌毒素污染预警系统与黄曲霉毒素污染早期甄别技术,变事后处理为事前预防,对花生黄曲霉毒素污染科学控制具有重要意义。本文基于连续十年花生主产区的黄曲霉毒素污染数据和产毒真菌分布以及收集权威机构发布的品种信息、产地土壤数据与气候数据,探明了我国花生黄曲霉毒素污染消长变化及其与影响因素的关联,建立了适于中国花生黄曲霉毒素污染状况的规则分类—平衡取样—随机森林预警模型,实现了对黄曲霉毒素污染的预测,为从源头开展花生黄曲霉毒素风险治理提供关键技术支撑。主要研究方法与研究结果如下:1.明确了我国花生AFB_1的空间分布特征,绘制出中国花生AFB_1的污染分布地图。主要结果如下:(1)采用地统计分析对中国花生污染的AFB_1进行时空分析,获得了中国花生黄曲霉毒素的污染模式,表明中国花生AFB_1呈现显著的时间和空间差异;(2)利用GS+10.0软件获得全方位的半变异函数,表明花生AFB_1含量具有空间自相关,2009-2017年花生AFB_1半变异模型为球状模型;(3)利用普通克里金插值获得未取样地区的待测值,生成中国花生AFB_1的污染分布地图,表明我国花生AFB_1污染高风险区域呈分散分布,与这些地区花生品种、土壤类型、菌种产毒力、气候条件等因素相关。2.分析了2009-2017连续九年我国花生黄曲霉毒素污染消长变化,并阐述了花生品种抗性、土壤类型、菌株产毒力及气候条件对其产生的影响。主要结果如下:(1)黄曲霉毒素污染存在年度间差异,2017、2014、2015和2013年黄曲霉毒素污染水平高于其它年份,总量污染水平分别为26.67±116.53μg/kg、17.49±72.01μg/kg、15.33±87.19μg/kg和13.95±52.65μg/kg。2009年和2010年黄曲霉毒素污染水平最低,总量污染水平分别为2.40±32.48μg/kg和4.32±35.15μg/kg。(2)黄曲霉毒素污染存在地区间差异,长江流域主产区的花生黄曲霉毒素污染水平最高,不同省份总量污染范围为5.59-37.15μg/kg。其次是南方主产区(总量污染范围为4.23-19.23μg/kg)和北方主产区(总量污染范围为1.78-9.11μg/kg)。东北主产区的污染水平最低,黄曲霉毒素总量不超过1.00μg/kg。(3)黄曲霉毒素低抗品种、砂土土质、高产毒力菌株在超标样品中占比较大,高温干旱区域与黄曲霉毒素污染严重区域重合,验证了各影响因素对花生黄曲霉毒素污染消长变化的影响,为建立收获前花生黄曲霉毒素污染大尺度宏观预警模型奠定基础。3.建立了适合中国的规则分类—平衡取样—随机森林花生黄曲霉毒素预警模型,并对模型进行了验证。主要结果如下:(1)确定了中国花生黄曲霉毒素预警模型的主要气候变量,即纬度、8-20时降水量、平均气压和日平均气温。利用各气候参数与AFB_1之间的最大互信息系数值排序,初选得到8个变量:纬度、平均地表气温、日最低地表气温、8-20时降水量、平均气压、日最低气压、日平均气温和日最高气温;综合考虑初选变量间的多重共线性与全面性,本研究最终确定建模气候变量为:纬度、8-20时降水量、平均气压和日平均气温。(2)建立了模型分类规则,即高纬度区域(≥40°N)、低纬度区域(≤21°N)与低温区域(日平均气温≤16.8oC)适合制定规则。通过分析主要气候变量与总样本及AFB_1污染样本分布的关系,可知高纬度区域(≥40°N)、低纬度区域(≤21°N)与低温区域(日平均气温≤16.8oC)没有超限样本分布,即为了提高模型精度可将此区域样品直接认定为未超限样品。(3)建立了平衡取样—随机森林模型,对花生AFB_1是否超限预测精度达95%以上。以筛选得到的4个主要气候因子与花生品种抗性、土壤类型及真菌产毒力为模型输入参数,利用平衡取样—随机森林方法进行分类预测。结果表明当分类阈值为20μg/kg和5μg/kg时,未超限样本预测准确率分别为98.94%和96.97%,超限样本分类准确率均为100%。同其他常用分类器(决策树、支持向量机、K—近邻法、BP神经网络、径向基函数神经网络)相比较,随机森林能有效分离超限与未超限样品,具有较高的预测精度。(4)明确了模型输入参数对花生AFB_1分类预测的贡献率排序。进一步分析输入参数在分类模型中的贡献率排序,结果表明当以20μg/kg和5μg/kg为阈值时,7个输入参数在随机森林分类模型中的贡献率排序依次为:纬度日平均气温8-20时降水量平均气压品种抗性土壤类型真菌产毒力。(5)初步验证了模型预测效果,其预测准确率为86.89%。利用建立的规则分类—平衡取样—随机森林预警模型对2018年收获的花生AFB_1含量进行预测,结果表明综合规则分类样本,当分类阈值为20μg/kg和5μg/kg时,样本分类准确率均为86.89%。
【学位授予单位】：中国农业科学院
【学位级别】：博士
【学位授予年份】：2019
【分类号】：TS222.1
【图文】：

黄曲霉毒素,化学结构,分解温度,氯仿

1.3 黄曲霉毒素（AFB1、AFB2、AFG1、AFG2、AFM1和 AFM2）的化re 1.3 Chemical structures of aflatoxins (AFB1, AFB2, AFG1, AFG2, AFM1表 1.1 黄曲霉毒素的熔点和紫外吸收Table 1.1 Melting-points and ultraviolet absorption of aflatoxins分子式分子量熔点（oC）紫外线吸收（λmax(nm)ε (C17H12O6312.3 268-269（分解温度） 223 25（不溶于氯仿） 265 13362 21C17H14O6314.3 286-289（分解温度） 265 11（不溶于氯仿-正戊烷） 363 23C17H12O7328.3 244-246（分解温度） 243 11（不溶于氯仿甲烷） 257 9 9264 10362 16C17H14O7330.3 237-240（分解温度） 265 9 7（不溶于乙酸乙酯） 363 21C17H12O7328.3 299（分解温度） 226 23（不溶于甲醇） 265 11357 19 2001

市县,花生,主产区,样品

1 花生样品1.1 花生取样范围与取样量鉴于气候类型、土壤类型、花生品种和种植方式等的差异，将中国的花生主产区分为为东北主产区（包括黑龙江省、吉林省和辽宁省）、北方主产区（包括河北省、山）、长江流域主产区（包括安徽省、江西省、江苏省、湖南省、湖北省和四川省）（包括广东省、广西壮族自治区和福建省）（图 2.1）。花生抽样范围覆盖四大花生的花生样品覆盖中国花生生产总面积的 94%以及总产量的 96%。采用分层随机取主产县根据种植面积和产量抽取 10-30 份花生样品，每个样品 30 公斤。为抽取的花代表该地区花生黄曲霉毒素污染状况，采用四分法分取样，至少每点取样样品 3k中要及时准确地记录样品信息（包括花生样品的抽样地点、品种名称、土壤类型、理情况等）。将抽取的带壳花生样品在一周内晾晒至水分含量低于 11%之后运至实于通风干燥的环境，待测黄曲霉毒素含量。1.2 花生样品的制备经过去杂、去尘及脱壳后，用切片机切片（厚度为 0.5cm），再经样品粉碎机粉碎， 0.9mm 孔径的样品筛，然后将粉碎后的样品装于样品瓶内 4oC 储存，采用高效液gh performance liquid chromatography，HPLC）检测黄曲霉毒素含量。

变异函数,平均值

25图 2.3 2009-2017 年花生 AFB1含量平均值的变异函数Figure 2.3 Variograms of mean AFB1content in peanut during 2009-2017

【相似文献】