湖泊沉积岩芯采集与年代框架建立
获取覆盖过去28,000年的沉积序列，为时间序列分析提供年代基础。
从西伯利亚和阿拉斯加7个湖泊采集8个沉积岩芯，采用已有年龄模型，将记录划分为14个2000年时间切片。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
更新世猛犸草原丧失背景下的潜在植物灭绝
Potential plant extinctions with the loss of the Pleistocene mammoth steppe.
实验完整度
高多湖沉积岩芯sedaDNA（8个岩芯504样本）与GBIF合成数据集、全球数据库比对验证潜在灭绝，含功能性状与统计建模。
主要模型
西伯利亚和阿拉斯加7个湖泊的8个沉积岩芯 过去28,000年湖泊沉积物sedaDNA时间序列 SibAla_2023自定义DNA参考数据库 GBIF合成数据集
重点核对
8个湖泊沉积岩芯，504个样本，覆盖至少过去28,000年 sedaDNA扩增子序列变异（ASV）与SibAla_2023数据库90–100%相似性匹配 通过共现网络检测10个至少含5个ASV的社区以区分真实分类单元 2000年时间切片及1000次重采样平衡岩芯数、样本数和读数 广义线性混合模型评估植被周转、大型动物灭绝和7月温度变化的驱动作用
摘要
在更新世-全新世过渡期，欧亚大陆北部占主导地位的猛犸草原生态系统消失，同时大型动物灭绝。然而，由于缺乏可识别的化石记录，植物灭绝模式很少被检测到。在此，我们介绍一种方法，用于检测区域（局部灭绝）到潜在全球尺度（灭绝）的植物分类单元丧失及其原因，该方法基于过去28,000年西伯利亚和阿拉斯加湖泊沉积岩芯中的古代植物DNA宏条形码（sedaDNA）。总体而言，潜在的植物灭绝与更新世-全新世过渡期的温度变化、植被变化和大型动物灭绝同步。估计的潜在植物灭绝率为每百万物种年1.7–5.9次灭绝（E/MSY），高于背景灭绝率但低于现代估计值。主要的潜在植物灭绝事件检测到约在17,000和9,000年前，滞后于最大植被周转。我们的结果表明，草本分类单元和对β多样性贡献较小的分类单元更容易灭绝。虽然随着DNA参考文库和古代sedaDNA数据的扩展，估计的稳健性将增加，但现有数据支持植物比哺乳动物对环境变化更具韧性。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 更新世-全新世过渡期猛犸草原丧失过程中，植物分类单元是否发生区域性丧失或潜在全球灭绝，其速率和驱动因素是什么？
- 核心机制
- 潜在植物灭绝与温度变化、植被周转和大型动物灭绝相关，其中植被周转（滞后一个时间切片）解释力最强，其次为大型动物灭绝和7月温度变化；草本分类单元和对β多样性贡献较小的分类单元更易灭绝。
- 主要证据
- 基于西伯利亚和阿拉斯加8个湖泊沉积岩芯504个样本的sedaDNA宏条形码数据，鉴定359个植物分类单元，其中17–59个为潜在全球灭绝；估计潜在灭绝率为1.7–5.9 E/MSY，主要灭绝事件约在17,000和9,000年前。
- 研究意义
- 研究提供了检测古代植物潜在灭绝的方法框架，表明植物比哺乳动物对环境变化更具韧性，并强调当前灭绝事件的独特性；结果有助于理解生态系统丧失、植物组成变化和级联灭绝之间的相互作用。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
沉积物古代DNA提取与trnL P6loop扩增
从沉积物中获取植物DNA，用于分类单元鉴定。
在清洁实验室条件下对沉积物子样本进行DNA提取，使用trnL g和h引物扩增P6loop区域，产物纯化后测序。
ASV生成与自定义数据库SibAla_2023比对
将测序读数转化为分类单元信号，并允许检测数据库未覆盖的潜在灭绝分类单元。
使用OBITOOLS处理双端读数，生成ASV，与SibAla_2023数据库进行90–100%相似性比对，区分dbASV和非dbASV。
共现社区检测与分类单元确认
从ASV中区分真实分类单元信号与PCR/测序错误。
基于Spearman相关（>0.4）构建网络，使用Louvain算法检测社区，将同一社区内非dbASV合并到dbASV。
时间切片重采样与分类单元出现/消失分析
平衡不同时间切片的采样努力，识别现代时间切片中缺失的分类单元。
对岩芯数、样本数和读数进行1000次重采样，定义lost dbtaxa和lost non-dbtaxa为在现代时间切片（0–2000年）中缺失的分类单元。
合成数据集评估数据库覆盖影响
评估SibAla_2023数据库覆盖不足对潜在灭绝信号的影响。
基于GBIF数据模拟分类单元获得和丧失，比较模拟与观察的db/non-db比例，估计数据库覆盖可解释的lost non-dbtaxa数量。
全球数据库比对与再避难所评估
评估lost non-dbtaxa是否可能在全球其他地区存在，从而区分局部灭绝与全球灭绝。
使用MegaBLAST对NCBI nt数据库进行比对，并通过GBIF评估科级物种重叠，估计研究区外再避难所的可能性。
潜在灭绝率估计与时间动态分析
量化潜在灭绝率并识别主要灭绝事件时间。
基于17–59个潜在灭绝分类单元计算E/MSY，利用重采样估计每个时间切片的期望再现率，识别潜在灭绝事件。
驱动因素与脆弱性特征分析
识别潜在植物灭绝的主要驱动因素和易灭绝分类单元特征。
使用广义线性混合模型评估植被周转、大型动物灭绝和7月温度变化与灭绝率的关系；比较lost non-dbtaxa与其他分类单元的社区大小、SCBD、系统发育距离和功能型比例。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
生物信息学与分类单元鉴定
全球数据库比对与灭绝评估
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| DNA提取 | PowerMax®土壤DNA分离试剂盒 | Qiagen | 12988-10 |
| 蛋白酶K | Karl Roth | 7528.1 | |
| 二硫苏糖醇 | Omnilab | 1198463 | |
| GeneJET PCR纯化试剂盒 | Thermo Fisher Scientific | K0692 | |
| DNA定量 | ds-DNA BR检测试剂盒 | Invitrogen | Q32853 |
| Qubit 2.0荧光计 | Invitrogen | Q32866 | |
| PCR扩增 | Platinum® Taq高保真DNA聚合酶 | Invitrogen | 10675783 |
| 脱氧核苷三磷酸 | Life Technology | R1121 | |
| 牛血清白蛋白 | VWR | 11761797 | |
| PCR产物纯化 | MinElute PCR纯化试剂盒 | Qiagen | 28006 |
| 沉积物采样 | DNA Exitus PlusTM清洁剂 | VWR | -- |
| 实验室操作 | UVC/T-M-AR DNA/RNA清洁箱 | BIOSAN | -- |
| 沉积物采样 | 8 mL无菌管 | Sarstedt | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 湖泊沉积岩芯采集与年代框架 | 岩芯数量、湖泊位置、年龄模型来源、时间切片长度 阅读提示：Methods中Site description and timeframe及Table 1 |
| DNA提取与PCR扩增 | 样品量、提取试剂盒、PCR循环数、引物、重复数 阅读提示：Methods中Sampling and sediment DNA lab work |
| ASV生成与数据库比对 | 相似性阈值、最低读数、样本出现频率、数据库版本 阅读提示：Methods中Data analysis及Supplementary Note 1–2 |
| 共现社区检测 | 相关阈值、社区最小ASV数、网络构建参数 阅读提示：Methods中Data analysis及Supplementary Data 2 |
| 重采样与时间切片分析 | 重采样次数、时间切片划分、现代时间切片定义 阅读提示：Methods中Data analysis及Supplementary Note 1 |
| 潜在灭绝率估计 | 潜在灭绝分类单元数、E/MSY计算公式、置信区间方法 阅读提示：Methods中Data analysis及Supplementary Note 3 |
| 驱动因素分析 | GLM模型结构、植被周转滞后、大型动物灭绝时间框架来源 阅读提示：Methods中Data analysis及Supplementary Note 4 |