《Biology》：Mitogenome-Informed Metabarcoding for Non-Invasive Identification of Sympatric Flamingos in Northern Chile
编辑推荐：
智利火烈鸟（Phoenicopterus chilensis）是一种广泛分布于南美洲的物种，栖息于从厄瓜多尔和秘鲁南至智利和阿根廷，东至巴西、乌拉圭和巴拉圭的多种湿地生态系统，分布范围从海平面到高海拔安第斯盐沼。尽管其分布广泛，但当地种群日益暴露于人为压力之下
智利火烈鸟（Phoenicopterus chilensis）是一种广泛分布于南美洲的物种，栖息于从厄瓜多尔和秘鲁南至智利和阿根廷，东至巴西、乌拉圭和巴拉圭的多种湿地生态系统，分布范围从海平面到高海拔安第斯盐沼。尽管其分布广泛，但当地种群日益暴露于人为压力之下，引发了区域尺度的保护关注。在本研究中，研究人员生成了线粒体基因组（mitogenome）资源，并开发了一种利用非侵入性样本进行物种鉴定的分子框架。研究人员对P. chilensis的基因组DNA进行了测序，并组装了其完整的线粒体基因组，从而能够对常用的元条形码（metabarcoding）标记进行计算机模拟评估。随后，研究人员从一个高安第斯泻湖的觅食地采集了八个粪便样本，该地P. chilensis与安第斯火烈鸟（Phoenicoparrus andinus）和詹姆斯火烈鸟（Phoenicoparrus jamesi）共存。作为阳性对照，研究人员从智利圣地亚哥国家动物园（Zoológico Nacional）采集了九个额外的粪便样本，该处仅饲养P. chilensis。提取DNA后，首先通过桑格测序（Sanger sequencing）进行初步筛选。为了稳健地分配分类学身份，研究人员实施了一个概率框架，结合最大似然法（maximum likelihood, ML）和贝叶斯推断（Bayesian inference, BI），基于观测到的遗传变异（genetic variants）估计物种成员的后验概率（posterior probability）。该方法揭示了两个无法归属于P. chilensis的遗传上独特的组，表明其他火烈鸟类群通过粪便DNA共存且可检测。研究结果证明了将线粒体基因组资源、元条形码和统计推断整合用于混合物种系统中非侵入性物种鉴定的价值。该框架为生物多样性监测提供了一个可扩展的工具，并对脆弱的安第斯高海拔生态系统中的保护策略和基于证据的政策制定具有直接影响。
**论文解读文章**
**研究背景与问题**
火烈鸟是南美洲湿地最具代表性的水鸟类群之一，作为特化消费者和环境变化的可见指示物种，它们发挥着重要的生态作用。在智利北部和安第斯山脉中部的阿尔蒂普拉诺高原，至少三种火烈鸟可能部分或局部共存：智利火烈鸟（Phoenicopterus chilensis）、安第斯火烈鸟（Phoenicoparrus andinus）和詹姆斯火烈鸟（Phoenicoparrus jamesi）。区域丰度与分布研究表明，这些种群中相当大一部分依赖于高安第斯和普纳湿地，这些生态系统对于觅食、季节性迁徙和繁殖至关重要。然而，这些湿地既高度独特且生产力丰富，又因依赖对气候变率和人为干扰敏感的水文平衡而脆弱。近年来，气候变化、地表水减少以及“锂三角”地区（包括智利北部）锂矿开采的扩张带来了日益严重的联合影响。火烈鸟已被提议作为高海拔湿地保护的旗舰物种，但其种群丰度可能因地表水减少和锂矿相关活动而下降，因而亟需开发稳健的监测工具。然而，在混合群体或从非侵入性样本中进行物种级别鉴定仍面临挑战：野外鉴定受限于空间共存、监测偏远湿地的后勤困难，而粪便等材料缺乏可明确分类的形态特征。尽管先前分子研究已澄清了火烈鸟科的进化关系，但可用的基因组参考资源仍相对有限，尤其是对于具有安第斯相关性的南美类群。完整的线粒体基因组（mitogenome）可提供高分辨率参考序列，用于系统发育推断、标记评估、分类分配以及解释通过元条形码（metabarcoding）产生的短扩增子，但高安第斯火烈鸟的线粒体基因组资源在文献和公共数据库中仍然不足。基于DNA的方法，如使用鸟类特异性元条形码引物组（MiBird）扩增线粒体12S rRNA基因短片段，已证明在直接观察不完整时可用于检测和鉴定鸟类。然而，针对南美火烈鸟的饮食研究在分类和地理覆盖上仍有限，分子方法评估饮食相对稀缺。因此，整合线粒体基因组参考资源、粪便DNA物种分配和饮食元条形码的分析框架，有望同时促进分类学、生态学和保护导向的监测。
**研究内容与结论**
本研究组装了P. chilensis的完整线粒体基因组，并以此为基础，结合分子方法分析了智利北部采集的未知火烈鸟粪便样本。首先，研究人员比较了传统桑格测序（Sanger sequencing）与基于线粒体引物的高分辨率元条形码方法。随后，对生成的数据进行评估，以识别材料中的遗传簇，包括一个可归属于P. chilensis的簇。总体而言，该研究为南美火烈鸟提供了新的基因组资源，并提出了一个可扩展的框架，用于对同域水鸟组合中的隐性遗传结构和营养生态进行非侵入性监测。论文发表在《Biology》。
**主要关键技术方法**
研究人员从P. chilensis肝脏组织提取基因组DNA，采用DNBSEQ-G400平台进行高通量测序，经参考引导和迭代映射策略组装完整线粒体基因组（17,450 bp）。随后在计算机中评估了MarVer3（16S rRNA）和MiBird-U（12S rRNA）两个元条形码引物组的结合位置和预期扩增子大小。粪便样本采集自两个地点：智利安托法加斯塔大区的Machuca泻湖（Laguna de Machuca）和智利圣地亚哥国家动物园（Zoológico Nacional），后者作为阳性对照。粪便DNA提取后，采用巢式PCR扩增12S和16S靶区域，并在Illumina iSeq 100系统上进行高通量测序。生物信息学处理使用DADA2流程进行扩增子序列变异（ASV）推断，通过最大似然法和贝叶斯推断的概率框架结合Dirichlet-多项分布模型进行物种归属分配。
**研究结果**
**3.1 P. chilensis线粒体基因组的组装及MarVer3标记的计算机评估**
高通量测序生成6,788,954对双端读长，组装出完整环状线粒体基因组，无缺口一致序列长度为17,450 bp，平均测序深度50.9×。计算机评估确认MarVer3和MiBird-U引物结合位点位于线粒体基因组内，支持其对粪便降解DNA的适用性。
**3.2 扩增子序列变异谱揭示三个主要火烈鸟遗传组**
粪便DNA的元条形码分析鉴定出火烈鸟相关的ASV，分为三个主要遗传单元：Phoenicopterus、Parina谱系A和Parina谱系B。国家动物园样本几乎全部由Phoenicopterus组主导，其中Phoenicopterus_01是主要ASV；Machuca泻湖样本中Parina谱系A占主导（ParinaA_01），而Parina谱系B由两个主要ASV（ParinaB_01和ParinaB_02）代表。
**3.3 遗传聚类支持Phoenicopterus与两个Parina谱系的分离**
基于ASV序列的近似最大似然树显示三个明确区分的遗传簇：Phoenicopterus、Parina A和Parina B。两个Phoenicopterus ASV聚集在一起，与Parina相关ASV明显分离；Parina谱系A和B也形成独立分支，表明粪便材料中存在至少两个遗传分化的Parina相关谱系。
**3.4 分层概率分配确认粪便样本的稳健物种级分类**
应用最大似然估计和贝叶斯Dirichlet模型，所有样本均以超过预定置信阈值（0.95）分配给单一主导遗传组，无混合、模糊或低信息样本。国家动物园样本全部归属Phoenicopterus；Machuca泻湖样本中LM-F01、F02、F04、F05、F06分配给Parina A，LM-F03、F07、F08分配给Parina B，表明该地点两个Parina谱系共存。
**3.5 次要次级信号被检测但未影响最终分配**
少数样本（如LM-F05和ZN-F12）含有微量次级组读长，但相对丰度低于1%，远低于混合样本阈值，被解释为技术或环境背景（如索引跳跃、低水平携带污染），而非混合粪便来源的证据。
**3.6 火烈鸟遗传组在采样地点的空间分布**
国家动物园样本仅检测到Phoenicopterus，验证了方法对预期P. chilensis信号的恢复能力。Machuca泻湖表现出最高遗传异质性，同时存在Parina A和Parina B，表明非侵入性粪便采样结合线粒体元条形码和概率分类能够检测到仅凭形态学难以区分的同域或局部共存的火烈鸟谱系。
**总结与讨论**
本研究开发了一个线粒体基因组指导的元条形码框架，用于智利北部火烈鸟粪便样本的非侵入性分子鉴定。通过生成P. chilensis的完整线粒体基因组，为评估元条形码标记和解释短粪便DNA序列建立了物种特异性参考。该方法成功恢复了火烈鸟相关的ASV，并利用ASV级分辨率和概率分类分配将Phoenicopterus与两个Parina相关遗传谱系区分开来。国家动物园样本一致归为Phoenicopterus，而Machuca泻湖样本揭示了两个操作性Parina谱系。这些结果表明，粪便DNA能够提供可靠的非侵入性证据，证明火烈鸟谱系在混合物种湿地系统中的存在。总体而言，该框架整合了线粒体基因组资源、高通量测序和统计分配，为脆弱的安第斯高海拔生态系统中的生物多样性监测、保护规划、栖息地利用研究和环境评估提供了有用工具。主要限制在于南美火烈鸟的参考数据库不完整，目前无法将Parina谱系A和B最终分配到P. andinus或P. jamesi。未来工作应纳入这些物种的完整线粒体基因组、验证过的标记序列和参考组织材料，扩大跨湿地、季节和年份的采样，并整合饮食元条形码标记，以将宿主身份与饮食和营养生态联系起来。这些改进将使所提出的框架从谱系级分配迈向对同域火烈鸟组合的更精确物种级监测。
