基于形态分类学与土壤DNA宏条形码评估欧洲森林凋落物与土壤栖息跳虫群落的环境驱动因子

《Applied Soil Ecology》：Environmental drivers of litter- and soil-dwelling springtail communities across European forests, assessed by morphotaxonomy and soil DNA metabarcoding

　　

跳虫（Collembola）在欧洲森林凋落物与土壤中数量丰富，在有机质转化与分解中发挥关键作用。然而，多数跳虫群落数据来自地方与区域尺度研究，大尺度研究匮乏。为填补该空白，研究人员采用标准化方法，在涵盖17国25个样点的泛欧（Europe-wide）土壤-气候（pedo-climatic）梯度上调查森林跳虫群落，评估哪些环境变量影响跳虫密度及α与β多样性。研究人员分别采集森林凋落物与下层土壤，联合应用形态学鉴定与土壤DNA宏条形码（soil DNA metabarcoding）。平均而言，凋落物栖息跳虫约2500个体·m?2，土壤栖息约10500个体·m?2，北极与北方区最高，地中海区最低。形态分类学下凋落物物种丰富度均值5.2、土壤7.0，宏条形码显示土壤丰富度更高（10.2）。形态学表明，正向影响凋落物跳虫密度与α多样性的关键因子为凋落物干质量，年均温（MAT）负向影响；正向影响土壤跳虫密度的关键因子为土壤氮（N）与黏土含量，容重负向影响。宏条形码数据得出类似关键因子。结果表明，在需要简便可靠方法快速评估跳虫群落谱变化的大尺度研究中，宏条形码是一种颇具前景的途径。

研究背景方面，土壤生物是生态系统关键组分，跳虫（Collembola）作为中型土壤动物（mesofauna）在有机质形成、养分循环与碳（C）循环中作用显著，但既有跳虫数据多源于局部或区域研究，缺乏跨欧洲森林梯度的大尺度对照。已有泛欧研究或混合多种土地利用类型，或未同时联用形态分类学（morphotaxonomy）与土壤DNA宏条形码（soil DNA metabarcoding）两种方法，难以支撑欧盟土壤监测指令所需的快速评估。为此，研究人员依托EU COST Action KEYSOM，在38°N–69°N、?8°W–27°E跨度内选取17国25个森林样点，覆盖北极、北方、阿尔卑斯、大西洋、大陆与地中海6个生物群区（biome），采集凋落物与矿质土壤，同步开展形态学鉴定与基于CO1 minibarcode的宏条形码分析，检验方法一致性、生境差异、环境驱动及指示类群四类假设。该研究发表于《Applied Soil Ecology》。

关键技术方法上，研究人员建立跨欧洲纬度梯度森林样本队列（25站点×4样方，2017年5–7月采样），凋落物经Macfadyen提取器分离跳虫并称重干物质，土壤分三核混合提取动物、另三核混合过4 mm筛供DNA与理化分析；形态学由专家镜检定种，宏条形码用MiteMinBarF7/R4引物扩增CO1约200 bp区段，Illumina MiSeq双端测序，QIIME2+DADA2推断ASV（amplicon sequence variant）并以筛选后BOLD库归类；统计上采用秩转换线性混合效应模型、随机森林回归（Random Forest）、CAP（canonical analysis of principal coordinates）、PERMANOVA、Procrustes叠合、NMDS（non-metric multidimensional scaling）与LEfSe指示类群分析，站点作随机效应以避免伪重复。

研究结果部分，3.1节森林站点属性显示各生物群区凋落物干质量、C:N比、土壤有机碳（SOM）、pH、黏土含量差异明显且群区内变异大。3.2节形态学与宏条形码比对中，Procrustes分析（P<0.001）表明土壤跳虫群落结构两法高度吻合，z-score下α多样性无显著差异（混交林Poduromorpha除外），证明宏条形码可替代形态学评估群落框架。3.3节密度与α多样性显示凋落物均密度2531±385个体·m^?2、土壤10689±962个体·m^?2，形态学凋落物丰富度5.2、土壤7.0，宏条形码土壤达10.2；Folsomia quadrioculata为落叶与混交林优势种，Parisotoma notabilis等偏好常绿林凋落物，Onychiuridae在宏条形码中占26%–37%。3.4节驱动因子中，秩混合模型示凋落物密度正关纬度与凋落物C:N、负关MAT，α多样性正关凋落物干质量；土壤密度正关土N与黏土、负关容重；宏条形码α多样性正关纬度与土C:N、负关MAT与pH。随机森林示MAP与凋落物干质量主导凋落物密度，黏土、海拔、容重主导土壤密度。3.5节β多样性CAP与NMDS表明六生物群区群落显著分离，MAT、纬度、林型、土pH、黏土等envfit显著，形态与宏条形码对环境轴响应一致。3.6节LEfSe指示类群中，形态学识别常绿林凋落物指示种Anurophorus laricis、Parisotoma notabilis、Ceratophysella denticulata，土壤中断叶林Folsomia quadrioculata、混交林Sminthurinus signatus、常绿林Pseudosinella alba突出；宏条形码仅获Isotomiella、Mesaphorura（落叶）与Sminthurinus（常绿）三属级指示。

讨论与结论翻译部分，研究人员指出假设1未获清晰结论：Procrustes与z-score α多样性支持两法可比，但宏条形码土壤丰富度高出约50%，且指示类群与部分环境相关不一致，限制来自BOLD仅覆盖约17%已知种（约1600/9400），本研约40% ASV无法定属；但宏条形码免个体提取、可捕隐存种，宜作大尺度监测补充。假设2成立——土壤密度约为凋落物4倍、丰富度与Shannon指数显著更高。假设3成立——地理气候（纬度、MAT）主控密度，局地气候与林型介入α多样性，土壤N、黏土、容重、pH、C:N驱动土壤群落。假设4部分成立——形态学给出多指示种，宏条形码仅至属级。结论上，metabarcoding represents a promising approach for future large-scale studies when simple and reliable methods are needed to rapidly assess shifts in springtail community profiles（宏条形码在需简便可靠方法快速评估跳虫群落谱偏移的未来大尺度研究中是一种颇有前景的途径）。