《Molecular Ecology Resources》：Substrate-Dependent Variation in Environmental DNA Persistence and Degradation From a Small Mammal
编辑推荐：
环境DNA（eDNA）分析应用于陆地环境已有十多年历史，然而关于eDNA沉积、迁移和降解动态的认识仍然有限。为解决这一问题，研究人员开展了一项为期一年的受控实验，评估了一种树栖哺乳动物的eDNA在四种底物（水体、沉积物、遮蔽巢箱表面和暴露巢箱表面）上的持久性。
环境DNA（eDNA）分析应用于陆地环境已有十多年历史，然而关于eDNA沉积、迁移和降解动态的认识仍然有限。为解决这一问题，研究人员开展了一项为期一年的受控实验，评估了一种树栖哺乳动物的eDNA在四种底物（水体、沉积物、遮蔽巢箱表面和暴露巢箱表面）上的持久性。研究人员在移除圈养个体后反复采样其沉积的eDNA，并使用物种特异性定量PCR（qPCR）和宏条形码（metabarcoding）进行分析。研究人员采用定量衰减、检测概率和条件衰减框架对降解过程进行建模，并应用随机断裂模型估算理论片段持久性。在水体和巢箱表面上，衰减遵循一级指数动力学，但在沉积物中则不然，后者表现出高度随机性。研究人员观察到不同底物间降解速率存在巨大差异，指数衰减模型显示半衰期从水体的7.3天到遮蔽内部表面的276天不等。水体中的快速衰减将检测限制在约77天内，而持久性模型预测遮蔽表面可保留可检测DNA超过4年（>1700天），表明在受保护环境中可能存在长期遗留DNA积累。虽然物种特异性qPCR比宏条形码具有更高的灵敏度，但两种方法的衰减模式平行。研究结果表明，在设计eDNA调查时，底物选择是一个关键决策，因为需要在最大化检测概率与保持时间分辨率之间进行权衡。
研究背景与科学问题
环境DNA（environmental DNA，eDNA）分析作为一种高灵敏度、非侵入性的物种检测手段，已广泛应用于生物多样性监测。然而，eDNA结果的解读受到对其沉积、迁移和降解动态认知不足的严重制约。eDNA的持久性是指生物体脱落的DNA片段在沉积后可被检测到的时间长度，反映了沉积、物理保护和环境降解三者相互作用的净结果。降解过程涉及微生物活动和酶解等生物机制，以及水解、紫外线（UV）暴露、温度、pH值和湿度等非生物过程与条件。不同生态系统条件对降解的影响各异，使得建模变得困难。例如，eDNA在被保护性土壤颗粒或生物膜结合时可能持久更久，而在永久冻土沉积物中甚至可保存数百万年；相比之下，水生环境中由于稀释作用和相对较高的微生物活性，DNA降解更为迅速。尽管已有少数研究比较了不同底物的降解速率，但尚无研究直接在陆地环境中对不同底物的降解速率进行比较。
若要将eDNA有效用于生态管理，必须理解不同底物间的降解速率差异。缺乏可靠的eDNA持久性模型会导致物种检测数据解读的不确定性，进而影响管理决策并危及物种保护。例如，快速降解可能导致假阴性，低估物种占有率或未能检测到脆弱种群；反之，过长的持久性则会引发歧义，使检测信号可能代表历史存在的残留信号，导致假阳性。这两种情况都会造成生态和经济后果——未能准确界定检测窗口可能导致无法发现濒危物种的存在，或耗费数百万美元去根除实际上已不存在的入侵物种。因此，将比较性持久性估计纳入eDNA监测的最佳实践指南至关重要。
此外，分子检测方法的灵敏度在eDNA数据解读中也扮演着重要角色。目前主要依赖两种方法：针对单个物种的物种特异性定量PCR（quantitative PCR，qPCR）和同时针对多个物种的宏条形码（metabarcoding）。宏条形码通常灵敏度较低，其原因包括扩增偏倚和引物效率低下。尽管已有研究比较了qPCR与宏条形码的灵敏度差异，但尚无研究直接对比两种方法随时间变化的eDNA信号持久性，这构成了时间数据解读的关键知识空白。
本研究以红尾帚尾袋鼩（Phascogale calura）为研究对象，考察其eDNA在四种陆地环境底物（水体、沉积物、遮蔽巢箱表面和暴露巢箱表面）中长达一年的降解动态，并利用物种特异性qPCR和宏条形码评估其时间可检测性。
主要技术方法
研究人员于西澳大利亚Dryandra国家公园捕获10只成年红尾帚尾袋鼩，将其单独圈养于塑料围栏内7至12天，使其与四种实验底物（沉积物、木制巢箱内外表面、雨水）充分接触以沉积eDNA。随后将围栏转移至红尾帚尾袋鼩分布区外的科廷大学本特利校区户外鸟舍，模拟自然环境条件进行为期380天的降解实验。在第0、1、3、9、18、27、50、81、140、240和380天采集样品，共获得484份样品。沉积物样品通过随机采集并均质化五个表层沉积物子样品获得；巢箱表面使用滚筒拭子取样；水样经0.45μm醋酸纤维注射器过滤器过滤后加入DNA Shield溶液保存。DNA提取分别采用改良磷酸盐缓冲液法（沉积物）、DNeasy Blood and Tissue试剂盒（水样）和磁珠提取法（表面拭子）。物种特异性检测使用针对细胞色素B（cytochrome B，CYTB）区域的RTP_CYTB_02探针法qPCR检测133 bp片段；宏条形码检测使用脊椎动物通用12S-V5引物扩增98 bp的12S rRNA基因片段，并在Illumina NextSeq 2000平台测序。数据分析采用线性混合效应模型（linear mixed-effects model，LMEM）、二项式广义线性混合效应模型（generalized linear mixed-effects model，GLMM）和贝叶斯广义线性模型（Bayesian GLM）进行多层面降解建模。
研究结果
3.1 eDNA衰减
指数衰减模型验证显示，暴露表面（R2=0.99）、水体（R2=0.85）和内部表面（R2=0.81）显著符合一级动力学，而沉积物因高度随机性（R2=0.05）不符合该模型。对于133 bp片段，水体半衰期最短（7.3天），暴露表面居中（43.7天），内部表面最长（276.0天）。线性混合效应模型表明，时间、底物及其交互作用均显著影响eDNA浓度。内部表面衰减率不显著；沉积物与内部表面无显著差异；而水体（β=-1.16，p<0.001）和暴露表面（β=-0.33，p=0.044）的降解显著更快。预测达到检测限（limit of detection，LOD）的时间分别为：水体77天、暴露表面254天、沉积物629天（外推值），内部表面超过1700天（1773天，外推值）。
3.2 eDNA持久性
内部巢箱表面在整个研究期间维持100%检测率，故排除于模型之外。广义线性混合效应模型显示，时间、底物及其交互作用显著影响检测概率。水体的持久性下降速度显著快于沉积物（β=-4.71，p<0.001）；暴露表面的总体检测概率显著低于沉积物（β=-1.91，p<0.001），但下降速率无显著差异。
3.3 条件衰减（阳性检测样本）
对阳性检测样本的线性混合效应模型显示，时间和底物均显著影响eDNA浓度，但两者无交互作用。各底物与内部表面相比，衰减率无统计学显著差异（β<0.78，p≥0.066），表明一旦检测到DNA，其在各底物中的降解速率相似。
3.4 宏条形码与qPCR的比较
在qPCR检测阳性的164份样品中，12S-V5宏条形码成功检出119份，真阳性率为72.6%，阳性预测率为97.5%。McNemar检验证实qPCR灵敏度显著高于宏条形码（χ2=37.53，p<0.001）。qPCR阳性但宏条形码阴性的样品平均浓度为40±7.7拷贝/反应，而两种方法均阳性的样品平均浓度为395±47.1拷贝/反应。两种方法均阳性的样品中，qPCR浓度与宏条形码读长计数呈显著正相关（ρ=0.5，p<0.001）。时间与方法无交互作用，表明两种方法的降解模式平行；三因素交互模型亦未发现底物间的方法差异。
讨论与结论
本研究首次提供了多种陆地底物间eDNA持久性的长期定量比较，证明底物类型是决定红尾帚尾袋鼩检测窗口的主要驱动因素。水体中半衰期仅7.3天，而遮蔽内部表面半衰期达276天，近40倍的差异导致检测窗口从77天延伸至超过4年。这一发现挑战了eDNA检测的二元解读方式，表明有效采样窗口显著受底物稳定性的影响。水体提供较短的检测窗口和高时间分辨率，而内部巢箱表面的多年代持久性虽延长了检测机会，但也增加了检测到不再存在于采样点的遗留DNA的风险。
沉积物eDNA的高随机性归因于土壤基质中eDNA沉积的空间异质性——动物行为导致的有机物质“热点”（如粪便、毛发、皮肤细胞）呈离散分布，点状采样难以一致捕获。相比之下，大面积连续采样的表面滚轮方法更易捕获剩余热点，维持较高检测率。暴露表面与内部表面的六倍差异进一步凸显了紫外线、风和降雨等环境暴露因素对eDNA稳定性的关键影响。
物种特异性qPCR灵敏度显著高于宏条形码（多检出约28%的样品），但两种方法展现出平行的衰减模式和显著的定量相关性，表明宏条形码在目标DNA浓度高于随机检测阈值时同样能准确反映相对时间衰减趋势。
研究提出，底物选择是一项战略性决策，决定了监测数据的时间分辨率和生态相关性。对于树栖物种管理，遮蔽巢箱可作为有效的“eDNA陷阱”长期积累分子证据，极大提高检测灵敏度。然而，为克服遗留信号的时序模糊性，建议采用多方法筛查验证策略：对eDNA阳性位点优先进行传统活捕或红外相机实地验证，同时探索使用多标记物测定不同片段长度比率来估算信号年龄的“eDNA时钟”方法。在实现更好的时间分辨率之前，研究建议生态管理者在采用eDNA技术进行存在/不存在调查时必须充分考虑底物动态，并将eDNA采样整合到长期巢箱监测计划中，以补充传统视觉和相机调查，提升对低密度种群或干扰后恢复状况的推断能力。
