《Molecular Ecology Resources》：Getting to the Core of the Matter—Assessing the Role of Replication in Metabarcoding-Based sedaDNA
编辑推荐：
重复设置是大多数实验和采样设计的核心，能够提高推断能力并捕获细尺度的数据异质性。然而，在一些生态学和进化研究背景中，其重要性仍缺乏充分评估。基于从沉积物档案中回收的DNA进行宏条形码研究即属此类情况——在这些研究中，生物信号通过沉积和埋藏过程整合生态信息，但通
重复设置是大多数实验和采样设计的核心，能够提高推断能力并捕获细尺度的数据异质性。然而，在一些生态学和进化研究背景中，其重要性仍缺乏充分评估。基于从沉积物档案中回收的DNA进行宏条形码研究即属此类情况——在这些研究中，生物信号通过沉积和埋藏过程整合生态信息，但通常每个位点仅依据单一沉积物岩芯进行推断。本研究利用嵌套采样设计，基于两个遗传标记（线粒体细胞色素c氧化酶I（COI）和核糖体18S rRNA基因）的沉积物DNA宏条形码数据，评估了不同类型重复设置的效果。该设计包含三个潮间带位点、每位点三个空间分离的沉积物岩芯（生物学重复）、每个岩芯两个沉积物层位，以及每个沉积物样本八个聚合酶链式反应（PCR）重复（技术重复）。方差分配结果显示，位点身份和沉积物年龄分组共同解释了β多样性中>70%的变异，表明位点间空间差异和地层差异是群落组成的主要驱动因素。PERMANOVA分析同样发现生物学重复的效应不显著。在来自同一沉积物样本的PCR重复之间，丰富度变化较大，而Shannon多样性更为一致。尽管存在这种变异性，技术重复间的群落组成差异仍小于与生物学重复或位点身份相关的差异，表明其对更宏观的生态格局影响有限。位点内重复岩芯之间的群落组成高度相似，与地层一致性相符。这些结果表明，位点内异质性有限，并提示在地层一致性条件下，增加生物学重复可能只能提供少量额外信息，而增强技术重复和地层分辨率则可以改善基于沉积物DNA宏条形码数据集的生态学推断。
**论文解读：基于宏条形码的沉积物古DNA研究中重复设置的作用评估**
**一、研究背景与科学问题**
实验和采样设计依赖于合理的重复设置，以确保生态学和进化推断的稳健性。充足的重复能够降低假阳性率、提高统计功效，并确保观察到的模式并非特定样本或条件的产物。在分子生态学中，重复设置是环境DNA（eDNA）和宏条形码数据集稳健性与可重复性的关键决定因素。重复主要分为两个层级：生物学重复（野外采样层面，如在同一采样点采集多个沉积物岩芯）和技术重复（实验室层面，如DNA提取、扩增和测序过程中引入的变异性）。尽管重复设置的重要性已被广泛认可，但在基于沉积物档案中回收的DNA（即沉积物古DNA，sedaDNA）开展的宏条形码研究中，不同重复类型的相对贡献仍知之甚少。沉积物古DNA代表了保存在沉积层中的、具有时间结构的eDNA亚集，能够记录过去生物群落的时空信息。然而，当前多数研究通常仅从每个位点采集单一沉积物岩芯，这忽略了位点内部可能存在的空间异质性，也使得生物学重复与技术重复的贡献难以区分。因此，系统评估不同层级重复在sedaDNA宏条形码研究中的作用，对于优化采样设计、提高生态推断的可靠性具有重要意义。
**二、研究设计与主要技术方法**
研究人员在西班牙北部圣托尼亚（Santo?a）河口系统的三个潮间带盐沼位点（Cicero, CCO；Colindres, CLR；Santo?a, STO）采集沉积物岩芯。于2023年11月，每个位点使用俄式泥炭采样器采集三个空间分离约50 cm的沉积物岩芯（生物学重复），每个岩芯选取浅层（4–5 cm）和深层（45–50 cm）两个沉积物层位（分别代表近期和较老沉积物）。研究人员采用嵌套设计：3个位点 × 3个岩芯 × 2个层位 × 8个PCR重复，每个遗传标记共144个样本。沉积物年龄通过
210Pb放射性测年法（以
210Po为代理，α能谱法测量）确定，并利用常量通量：常量沉积（CF:CS）模型估算沉积速率。使用手持式X射线荧光（XRF）扫描仪对每个生物学重复岩芯进行元素组成分析，以评估岩芯间的地层一致性。DNA提取采用DNeasy PowerMax Soil Kit。两个宏条形码标记分别为线粒体细胞色素c氧化酶I（COI）基因的313 bp片段（Leray-XT引物对）和核18S rRNA基因V7区（约110 bp，18S_allshort引物）。采用双标签（twin tagging）策略构建文库，在Illumina NovaSeq 6000平台测序。生物信息学流程包括Cutadapt去 multiplex、DADA2降噪、LULU错误ASV校正、BLAST+比对NCBI nt数据库进行物种注释。统计分析采用线性混合效应模型（lmer/glmer.nb）评估α多样性，PERMANOVA、距离冗余分析（db-RDA）和非度量多维尺度分析（NMDS）评估β多样性，并利用联合物种分布模型（jSDM，基于HMSC框架）进行方差分配。
**三、研究结果**
*3.1 定年与地球化学相关性：* 放射性测年显示所有位点均为近期沉积，但沉积速率差异显著：CLR为0.63 cm/年，CCO为0.53 cm/年，STO仅为约0.1 cm/年。据此推算，所分析的层位覆盖约最近500年，其中浅层对应近几十年（STO：1974年；CCO：2014年；CLR：2017年），深层则对应20世纪中叶（CCO：1939年；CLR：1952年）或更老的时期（STO约16世纪）。标准化XRF剖面显示，各位点内生物学重复岩芯之间存在中等到强的地层一致性，尤其是碎屑元素（Ti、Rb、Zr），为岩芯间的横向连续性提供了外部验证。
*3.2 测序产出与数据概览：* 测序获得COI 101,370个ASV、18S 57,415个ASV，质控后分别保留65,466和51,263个ASV。每个PCR重复的平均测序深度为COI 226,613 reads（SD=309,124）、18S 468,673 reads（SD=801,189）。基于覆盖度的分析（iNEXT）表明测序深度足以恢复每个PCR重复内可检测的多样性。稀释敏感性分析显示，尽管绝对丰富度数值降低（尤其COI），但总体模式与混合效应模型结论保持一致，表明生物学解释对测序深度变异不敏感。阴性对照读取数极低，表明污染控制有效。
*3.3 分类组成：* COI和18S均能回收后生动物、植物、原生生物和真菌，但相对丰度谱差异明显——COI中未分配ASV比例较高（53.08%），而18S中原生生物占最大比例（40.55%）。表层（4–5 cm）的后生动物和原生生物ASV比例较高，深层（45–50 cm）则相对富集植物。位点内生物学重复在界水平上表现出相似的图谱，表明位点内变异较低（STO除外）。在仅分析后生动物时，主要门类为环节动物、线虫、节肢动物和软体动物，但某些样本中这些门类的相对丰度在生物学重复间存在一定异质性。
*3.4 α多样性：* ASV丰富度的技术重复变异性大于Shannon多样性。Shannon多样性的组内相关系数（ICC）为0.43–0.69，而ASV丰富度的ICC仅为0.18–0.31，表明丰富度在PCR重复间可重复性较低。在考虑了技术变异性后，岩芯层面生物学重复对整体真核群落α多样性的影响有限，显著的岩芯间差异主要集中于STO位点。位点层面差异更为显著：STO的丰富度始终高于CCO和CLR，而Shannon多样性的位点差异在深层沉积物中更明显。在后生动物子集中，生物学重复的影响同样很小，位点效应则相对稳定。
*3.5 群落结构与β多样性：* NMDS排序显示，两个标记的群落组成均呈现显著的空间和地层结构，位点内生物学重复紧密聚集，分散程度极小，表明岩芯间群落组成高度相似。PCR重复间的Bray–Curtis相异度通常低于或相当于位点内生物学重复间的相异度，说明技术随机变异主要贡献于沉积物样本内部的变异性。jSDM方差分配显示，位点解释了63%–77%的变异，沉积物年龄组解释了14%–28%，而生物学重复相关的随机效应始终低于6%。db-RDA证实位点和沉积物年龄组对COI（位点F=3.33, p=0.001；年龄F=3.86, p=0.001）和18S（位点F=5.01, p=0.001；年龄F=5.18, p=0.001）均有显著效应。PERMANOVA分析中，位点（COI R2=0.272, p=0.001；18S R2=0.345, p=0.001）和沉积物年龄组（COI R2=0.157, p=0.001；18S R2=0.177, p=0.002）均显著，而生物学重复不显著（COI R2=0.032, p=0.712；18S R2=0.029, p=0.516）。后生动物子集分析中，生物学重复仍不显著，但位点和年龄组效应减弱或呈现标记依赖性。
**四、讨论与结论**
讨论部分指出，在所有位点和遗传标记中，岩芯层面的生物学重复对群落组成变异的贡献远小于位点身份和沉积物年龄组。这一模式在PERMANOVA、NMDS和jSDM中均得到一致验证。结果支持沉积物DNA档案具有信号整合特性，即生态信号通过沉积、再悬浮和埋后过程在空间上被混合，但仍保留足以解析时间生态格局的地层结构。
关于技术重复的讨论表明，其影响主要取决于多样性指标而非遗传标记：Shannon多样性对技术重复稳健，而ASV丰富度对技术重复敏感，这归因于稀有类群检测的随机性。因此，丰富度估计应谨慎解释，而丰度加权指标（如Shannon多样性）具有更强的鲁棒性。
关于生物学重复的讨论强调，XRF地球化学数据的横向一致性支持了低生物学方差并非由地层错位或采样伪迹所致的判断。位点内岩芯的群落组成高度相似，仅在STO位点（沉积速率较低）观察到局部异质性，可能反映了该位点沉积物输入减少且更异质的环境条件。
研究还讨论了生态驱动因素：位点间差异可能反映局部栖息地特征、周围植被、沉积物来源、水动力条件和与邻近环境的连通性等综合作用；沉积物年龄组则可能代表生态时间更替与年龄相关的DNA保存差异的叠加效应。
研究结论指出：在地层一致性条件下，生物学重复对基于宏条形码的sedaDNA分析贡献相对较小；位点身份和沉积物年龄组是变异的主要驱动因素。因此，在明确控制主要分析变异性来源的情况下，研究投入可更多导向空间覆盖、地层分辨率、测序深度和技术重复，而不是增加位点内生物学重复。然而，生物学重复不应被完全省略，因为一定水平的位点内重复仍是量化局部变异性并支持位点间稳健比较所必需的。该研究发表于《Molecular Ecology Resources》，为优化sedaDNA宏条形码研究的采样设计提供了实证依据，强调重复策略应结合沉积学背景和研究目标进行合理调整。
