一种基于eDNA的黑土土壤动物空间变异监测方法

本发明涉及土壤生态监测，尤其涉及一种基于edna的黑土土壤动物空间变异监测方法。

背景技术：

1、黑土土壤动物是土壤食物网的核心组成部分，直接参与有机质分解、微生物群落调控、养分矿化及土壤团聚体形成等关键生态过程。不同地理位置及不同农业管理措施（如耕作方式、施肥制度）会导致土壤动物群落在组成和丰度上产生显著的空间异质性。因此，精准监测黑土土壤动物群落的空间变异，对于揭示农业管理措施与地理位置对土壤动物群落结构的影响、建立可重复的土壤动物监测技术体系具有重要意义。

2、现有监测手段主要依赖物理分离、形态学鉴定，存在通量低、周期长、对专家依赖度高且难以标准化等局限。环境dna（edna）技术提供了一种非破坏性的高通量替代方案，其常规流程包括样品采集、dna提取、目标区域扩增、高通量测序及下游的多样性分析。

3、环境dna技术能够通过检测土壤样品中残留的生物遗传物质，对土壤动物群落进行非破坏性和高通量监测。通常的edna分析流程包括样品采集、dna提取、扩增测序、序列质控、otu聚类、分类学注释以及α多样性、β多样性和差异分析。

4、然而，现有技术存在以下缺陷和不足：（1）常规edna分析通常将整体群落的α/β多样性作为核心结果，这往往会掩盖那些丰度不高、但在多个地点对特定处理具有一致响应趋势的关键指示物种（otu）；（2）现有方法在进行组间比较时，通常直接进行总体混合比较（如将所有地点的处理组与对照组直接合并），缺乏先同点配对、后跨点验证的分析逻辑，导致分析结果极易受到地点间固有环境噪声的干扰，降低了处理效应检测的敏感度和准确性；（3）在实际采样中，不同地点的处理设置往往不完全对称，部分地点可能缺失特定对照组，现有标准流程对此类非平衡数据缺乏有效处理机制，强行分析会导致样本配对失当；（4）现有技术无法输出兼具处理响应方向一致性和跨地点空间稳定性的指示性生物模块，难以直接为黑土碳循环研究提供空间直观的生物学证据；（5）当缺乏精确的经纬度坐标时，传统的距离衰减模型和空间自相关分析失效。

5、因此，如何提供一种无需依赖土壤理化指标和地理距离数据，仅基于edna otu丰度与分类信息，即可有效识别黑土土壤动物中具有稳健空间变异特征的指示类群的方法，是本领域亟待解决的技术问题。

技术实现思路

1、针对上述现有技术的不足，本发明的目的在于提供一种基于edna的黑土土壤动物空间变异监测方法，通过构建地点内配对比较—跨地点方向一致性筛选—空间响应模块聚合的分析链路，解决现有技术中低丰度指示类群易被掩盖、跨地点比较受环境噪声干扰、非平衡采样数据难以处理以及缺乏空间稳定指示性生物模块输出的技术问题。

2、为了实现上述目的，本发明是采用以下技术方案实现的：所述方法包括：

3、获取不同黑土样点、不同田间处理下的土壤样品的环境dna测序数据，构建otu丰度表、分类注释表和样品元数据，并对所述otu丰度表和所述样品元数据进行数据清洗；

4、依据所述分类注释表，采用分层筛选策略从所述otu丰度表中筛选目标土壤动物otu；

5、对筛选出的所述目标土壤动物otu的丰度进行基于hellinger转化的标准化处理，得到标准化丰度；

6、针对同一地点内同时存在处理组和对照组的样品，根据所述标准化丰度计算各地点内各所述目标土壤动物otu的处理响应值；

7、汇总各所述目标土壤动物otu在所有有效地点的所述处理响应值，从响应效应强度、响应方向一致性、有效地点覆盖率和留一法交叉验证稳定性四个维度计算空间稳定响应综合评分，并根据所述空间稳定响应综合评分筛选跨地点稳定指示otu；

8、基于所述跨地点稳定指示otu的分类信息和跨地点响应方向向量进行模块聚合，构建土壤动物空间变异模块；

9、输出土壤动物空间变异监测结果；

10、其中，所述处理响应值用于表征同一地点内处理组相对于对照组的相对丰度变化，所述有效地点为能够计算出有效处理响应值的地点。所述处理组为接受待评估田间管理措施的样品组，所述对照组为同一地点内接受常规田间管理措施的样品组。

11、在一个方案中，所述对所述otu丰度表和所述样品元数据进行数据清洗，包括：

12、强制要求所述样品元数据中的样品编号与所述otu丰度表的样品编号完全匹配；

13、对于同一地点下仅存在单一样品或处理标签不明晰的样品，将所述样品剔除出所述处理响应值的计算池，仅保留用于描述地点物种组成的地点效应分析。

14、在一个方案中，所述采用分层筛选策略从所述otu丰度表中筛选目标土壤动物otu，包括：

15、保留界级为后生动物metazoa的otu；

16、在所述后生动物metazoa内部锁定与黑土生态系统功能密切相关的目标类群，所述目标类群包括线蚓目enchytraeida、石蜈蚣目lithobiomorpha、小杆目rhabditida、啮虫目psocodea、中气门目mesostigmata及旋轮目philodinida中的至少之一；

17、筛除未注释到有效纲或目水平的otu，以及明确属于植物、真菌或原生生物的otu。

18、在一个方案中，所述对筛选出的所述目标土壤动物otu的丰度进行基于hellinger转化的标准化处理，包括：

19、基于每个样品的测序总丰度计算相对丰度，并在计算所述相对丰度时引入拉普拉斯平滑修正因子，所述拉普拉斯平滑修正因子为小于1的正数；

20、对所述相对丰度进行平方根转化，得到hellinger距离矩阵的基础值，作为所述标准化丰度。

21、在一个方案中，所述拉普拉斯平滑修正因子取0.65。

22、在一个方案中，所述根据所述标准化丰度计算各地点内各所述目标土壤动物otu的处理响应值，包括：

23、分别计算同一地点下处理组和对照组中各所述目标土壤动物otu的标准化丰度均值，以所述处理组的标准化丰度均值与所述对照组的标准化丰度均值的对数响应比作为该目标土壤动物otu在该地点的处理响应值；

24、所述处理响应值大于零表示该目标土壤动物otu在处理组富集，为正响应；所述处理响应值小于零表示该目标土壤动物otu在对照组富集，为负响应。

25、在一个方案中，所述响应效应强度为各所述有效地点的处理响应值绝对值的中位数；

26、其中，响应方向一致性为目标土壤动物otu的处理响应方向符号在其全部有效地点上求和后取绝对值、再除以该otu的有效地点数量得到的比值；所述响应方向一致性的取值范围为0至1，当所述响应方向一致性大于0.8时，定义对应的目标土壤动物otu具有强跨地点一致响应性；

27、所述有效地点覆盖率为能够计算出有效处理响应值的地点数量与全部地点数量的比值；

28、所述留一法交叉验证稳定性通过逐一剔除一个地点、基于剩余地点的处理响应值重新计算总体响应方向获得；每剔除一个地点记为一次留一检验，若该次留一检验中重新计算的总体响应方向与基于全部有效地点计算的方向一致，则记为一次稳定；所述留一法交叉验证稳定性为稳定次数与留一检验总次数的比值。

29、在一个方案中，所述空间稳定响应综合评分为归一化后的所述响应效应强度、所述响应方向一致性、所述有效地点覆盖率与所述留一法交叉验证稳定性的乘积；

30、所述根据所述空间稳定响应综合评分筛选跨地点稳定指示otu，包括：筛选所述空间稳定响应综合评分排名前10%的目标土壤动物otu作为所述跨地点稳定指示otu；对于仅在单一地点计算出有效处理响应值但处理响应值绝对值极高的目标土壤动物otu，标记为地点特异响应otu，且不纳入所述模块聚合。

31、在一个方案中，所述基于所述跨地点稳定指示otu的分类信息和跨地点响应方向向量进行模块聚合，包括：

32、将每个所述跨地点稳定指示otu在各地点中的处理响应方向转化为响应特征向量；

33、当两个otu在分类学上至少属于同一科或同一目，且二者的所述响应特征向量的杰卡德相似系数大于0.75时，将所述两个otu聚合为同一土壤动物空间变异模块；

34、聚合成功后计算模块丰度，所述模块丰度为模块内各otu相对丰度之和；

35、对聚合后的模块递归执行地点内处理响应值计算和跨地点稳定响应评估，并根据响应模式将输出结果划分为地点响应型、处理响应型、地点依赖型和稳定指示型四类。

36、本发明有益效果：

37、（1）仅需edna测序产出的otu表和基础元数据，无需额外测定昂贵的土壤理化指标或进行样点gps定位，实施成本低，适用范围广。

38、（2）通过严格的样品匹配和目标类群筛选，排除了非目标生物和异常标签样品的干扰，显著提高了下游分析的准确性。

39、（3）有效分离了处理效应和地点环境效应，相比传统直接混合比较的方法，能更纯净地提取出处理措施对生物群落的真实影响。

40、（4）引入了响应效应强度、响应方向一致性、有效地点覆盖率和留一法交叉验证稳定性四个评价指标，特别是独创的加权综合评分模型，能够筛选出在空间上高度稳健、生态指示意义明确的低丰度指示otu。

41、（5）基于分类地位与响应方向向量双重约束的聚合策略，将零散的otu转化为具有生态学意义的功能响应模块，使监测结果更宏观、更具解释力；在相同数据量下对处理效应的检测灵敏度更高，在非平衡采样下仍可完成地点内配对分析，监测结论的可重复性更强。

42、（6）明确提出了对缺失配对处理地点的处理方案（仅用于地点效应分析），避免了强行构造虚假对照导致的统计错误。

43、（7）通过留一法交叉验证的引入，本方法在评价指示otu时天然具有排除异常值干扰的能力，提高了监测结论的可靠性和鲁棒性。