《BMC Genomics》：Genomics as a time capsule: insights from Oreobates chiquitanus type specimens
编辑推荐：
摘要背景：生物多样性面临日益严重的风险，两栖动物是最受威胁的脊椎动物类别，全球40.7%的物种处于危险之中。对于两栖动物而言，其数量下降的主要原因是栖息地丧失。这一全球性问题的规模在所谓的森林砍伐热点地区最为明显，例如玻利维亚的奇基塔诺旱林（Chiquitan
摘要背景：生物多样性面临日益严重的风险，两栖动物是最受威胁的脊椎动物类别，全球40.7%的物种处于危险之中。对于两栖动物而言，其数量下降的主要原因是栖息地丧失。这一全球性问题的规模在所谓的森林砍伐热点地区最为明显，例如玻利维亚的奇基塔诺旱林（Chiquitano Dry Forest），该地区拥有高度多样化的动植物群，许多物种仍有待正式描述。最近，青蛙Oreobates chiquitanus才在该森林内的一个地点被描述。此后，仅在另外两个地点发现了它。模式产地正在进行的伐木活动导致采样点被砍伐，甚至在该物种被正式描述之前就已发生。所有已知地点周边的持续伐木活动引发了关于该物种命运的疑问。在此，研究人员提供了该物种的全面遗传定义，并试图评估模式种群的遗传多样性和保护价值。结果：研究人员对所有可用的Oreobates chiquitanus模式标本进行了测序，提供了一个副模式产地标本（paratopotype）的参考基因组，以及所有模式标本的线粒体组（mitochondrial assemblies）。研究人员进一步表明，该物种所有已知种群都靠近近期伐木活动。种群历史动态的重建表明，该种群大约在30万年前（300kya）从一个低谷中恢复。亲缘关系和杂合度的估计暗示该种群在遗传上是活跃的。结论：本研究例证了自然历史收藏与基因组学计划（如TBG，即Translational Biodiversity Genomics）合作的重要性，以便更好地理解全球变化背景下的人为影响，识别并描述生物多样性，提高意识，并为潜在的保护措施提供信息。研究人员不仅为O. chiquitanus提供了全面且永久的分子定义，还评估了该种群的遗传多样性。因此，本研究为该物种未来遗传侵蚀（genetic erosion）的研究提供了重要基线。此外，它还说明了奇基塔诺旱林生物多样性丧失的风险，并强调了保护工作的必要性。
论文解读文章
**研究背景与问题**
当前，生物多样性丧失已成为全球性普遍现象。据估计，自上个世纪以来生物多样性下降了2%至10%，全球约有100万物种面临灭绝风险。驱动生物多样性变化的五大主要因素均与人类活动密切相关，其中以陆地和海洋利用变化影响最为显著。农业扩张作为最普遍的土地利用形式，影响了全球三分之一的地表，对陆地生态系统和物种产生巨大冲击。两栖动物是受威胁最严重的脊椎动物类群，全球40.7%的物种处于危险中，其数量下降的首要驱动因素正是栖息地丧失，而合法与非法的农业扩张导致的栖息地破坏是主要原因。这种区域性的生物多样性丧失在热带森林砍伐热点地区尤为突出，例如玻利维亚的奇基塔诺旱林（Chiquitano Dry Forest, CDF）。该森林是玻利维亚东部低地的特有种区，拥有高度多样化的动植物群，但研究尚不充分，许多物种仍等待正式描述。自1985年以来，由于农业扩张，近四分之一的区域面积减少，森林砍伐和野火的速率自2020年以来加速。正是在这样的背景下，青蛙Oreobates chiquitanus于近期被描述，其唯一已知模式产地即为该森林内的一处地点。然而，在物种描述正式发表之前，当地土地所有者发起的伐木活动已导致该采样点被彻底砍伐。此后，该物种仅在其他两个地点被发现，且所有已知地点周边均存在持续的、大规模的伐木活动。模式产地的毁灭使得从模式种群获取更多材料成为不可能，阻碍了未来的分类学或分子研究。模式标本是定义物种的物理参考，其保存对分类稳定性和物种鉴定至关重要。但模式产地本身会随外部影响而改变，其代表性无法长久保证。因此，本研究旨在利用基因组学手段，从可获得的O. chiquitanus模式标本中永久保存该物种的遗传定义，并评估模式种群在栖息地破坏前的遗传多样性，为未来研究遗传侵蚀和保护行动提供基线。
**研究人员开展的研究与结论**
研究人员对所有可用的Oreobates chiquitanus模式标本（共4份，存放于德国法兰克福森肯伯格博物馆）进行测序，利用一份副模式产地标本（SMF88498）通过长读长测序构建了参考基因组，并利用短读长数据对包括正模在内的所有模式标本进行重测序，同时组装了所有标本的线粒体基因组。研究进一步绘制了该物种已知三个分布点与伐木活动的关系图，显示所有种群均位于近期大规模伐木活动10公里范围内。通过重建种群历史动态（PSMC分析），发现模式种群在大约100万年前（Mya）经历了一次扩张，随后在30万年前（kya）左右降至最低有效种群大小，之后有所恢复。亲缘关系和杂合度分析表明，尽管所有标本采集自同一时间和地点，但线粒体单倍型差异暗示不同母系来源，且基因组水平杂合度较高（0.492%），与其他青蛙物种比较显示该种群遗传多样性较为健康。研究得出结论：研究人员为O. chiquitanus提供了全面且永久的分子定义，并评估了模式种群在破坏前的遗传状态，证明该种群在遗传上并非贫瘠或近交，但受持续伐木威胁，其生存前景堪忧。该研究发表于《BMC Genomics》。
**主要关键技术方法**
研究采用的长读长测序技术（PacBio Sequel II）对一份副模式产地标本（SMF88498）进行CLR测序，用于从头组装参考基因组；同时对所有四份模式标本（SMF88497–SMF88500）进行短读长测序（Illumina NovaSeq6000，150 bp双端文库，插入片段350 bp），平均测序深度约10倍（重测序标本）或30倍（参考基因组标本）。使用Flye v2.9进行基因组组装，并通过短读长数据进行抛光校正（DeepVariant、bcftools）。使用BUSCO、Quast和blobtoolkit评估组装质量。采用GetOrganelle从短读数据中组装每个标本的线粒体基因组，并用MitoZ进行注释和环化。种群分析方面，利用PSMC（Sequentially Markovian Coalescent）模型估算历史有效种群大小（假设世代时间3年，突变率0.3×10??每核苷酸每代）；使用ANGSD和ngsRelate计算亲缘关系和近交系数；通过随机降采样至10倍覆盖率后，采用GATK进行SNP calling并估计全基因组杂合度，并与三个额外蛙类物种比较。样本来源于森肯伯格博物馆的4份模式标本。
**研究结果**
- **Genome assembly（基因组组装）**：最终参考基因组长度约3.4 Gbp，contig N50约332 kbp。BUSCO分析显示84.3%的完整基因，6.4%破碎，9.3%缺失，符合两栖动物低分类覆盖背景下的较连贯组装。重复序列占72.92%，其中71.06%为散在重复，主要为未知重复家族。重复景观显示DNA重复在较远时期活跃，而长末端重复序列（LTRs）和长散在核重复序列（LINEs）在近期活跃。
- **Population analyses（种群分析）**：PSMC分析揭示种群大小约100万年前达峰值，之后下降至约30万年前最低点，随后恢复。所有标本的PSMC曲线相似，表明共享相同历史。Jacquard系数分析显示所有标本间无亲缘关系（J9=1，其余为0），但R0、R1和KING分析均将标本鉴定为三级亲属。全基因组杂合度在O. chiquitanus中为0.492%，高于H. boettgeri（平均0.157%）和L. fallax（0.041%），但低于S. bombifrons（1.119%）。线粒体基因组组装显示所有标本的线粒体单倍型不同（原始距离0.04%–0.09%），证实不同母系来源。每个线粒体基因组编码13个蛋白质编码基因、2个rRNA和23个tRNA（tRNA-Leu、tRNA-Ser和tRNA-His重复）。对于标本SMF88499，组装出两个不同的线粒体基因组，差异主要存在于控制区的插入缺失，可能源于异质性（heteroplasmy）而非NUMTs。伐木历史分析显示，所有三个已知O. chiquitanus种群距离近期大规模伐木活动均小于10公里。
**讨论与结论**
讨论部分指出：通过所有模式标本的测序，研究人员提供了该物种全面而永久的遗传特征，这对模式产地已遭破坏的情况尤其重要。虽然无法评估模式种群的当前状态，但当地的伐木活动已造成景观破碎化，可能也导致种群破碎化。PSMC分析无法反映近期种群历史，但亲缘关系和杂合度分析表明，在栖息地破坏前模式种群可能仍保持遗传多样性且非近交，尽管样本量小且存在背景亲缘关系的可能。线粒体差异和较高的杂合度支持这一判断。研究强调，在生物多样性丧失的时代，自然历史博物馆保存着现生和灭绝物种的珍贵收藏，基因组数据的使用极大地便利了保护生物学。结论部分翻译如下：研究结果提供了所有模式标本的永久而全面的基因组资源，即Oreobates chiquitanus的分子指纹。因此，本研究凸显了在生物多样性全球丧失的时代，自然历史博物馆的重要性，因为它们保存着现存和已灭绝物种的宝贵收藏。它强调了如何利用这类数据在栖息地退化和丧失之前捕捉O. chiquitanus模式种群的遗传结构和多样性信息。O. chiquitanus的案例强调了在保护生物学中使用基因组数据的好处，并例证了为减缓奇基塔诺旱林生物多样性丧失而采取保护行动的必要性。
