本文所述方案旨在指导用户对从甲虫中获得的细胞色素c氧化酶亚基I(COI)基因序列进行基于计算机的处理,从而能够从DNA数据中生成称为分子操作分类单元(MOTUs)的物种聚类假说。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
本文所述方案旨在指导用户对从甲虫中获得的细胞色素c氧化酶亚基I(COI)基因序列进行基于计算机的处理,从而能够从DNA数据中生成称为分子操作分类单元(MOTUs)的物种聚类假说。
生物多样性衰退正以难以想象的速度发生,然而生物多样性监测却受到诸多因素的阻碍,包括基线多样性本底数据不完整以及熟练分类学家数量不断减少。弥合物种编目中的这一障碍 "隐存分类群",或在分类学中常被忽视的高多样性类群,在面临巨大环境压力的热带地区尤其重要。尽管基于形态学的初级分类学仍是物种鉴定的金标准,基于DNA的方法在加速物种界定与发现方面已展现出巨大潜力。本文介绍了一种用于生成分子操作分类单元(MOTUs)的整合分析流程,该流程整合了在 两个极为不起眼的甲虫属,即 Byrrhinus Motschulsky,1858,属于扁甲科 Anacaena Thomson, 1859(隶属于水龟甲科 Hydrophilidae),采自菲律宾。由于本实验方案聚焦于DNA条形码的处理,且这些遗传数据既可通过测序获得,也可从在线数据库下载,因此本方案的起始点为DNA序列。在获得COI条形码后,该分析流程结合使用三种基于阈值的分子物种划分方法,即TaxonDNA、K2P和ASAP,这些方法以序列比对结果作为输入数据。此外,该流程还进一步采用三种基于溯祖理论的方法,即PTP-ML、PTP-BI和mPTP,这些方法需要以系统发育树文件作为输入数据。部分方法需借助网络服务器完成,而其他方法则使用本地软件。通过所选的算法组合,可基于共识或多数原则识别出分子操作分类单元(MOTUs)。值得注意的是,该分析流程成功区分了来自菲律宾的多种显著不同的甲虫物种,包括那些采自相同或相邻地区的物种。
鉴于全球不同生物类群的物种正以前所未有的速度衰退1,2,3,4,全面清查生物多样性的紧迫性日益凸显,这一工作包括发现和描述此前未知的物种,已成为与时间赛跑的任务2,5。一方面,生物多样性保护只有在具备充分的分类学知识基础上才具有实际意义6,7。因此,在物种水平上鉴定生物多样性十分必要,因为亲缘关系相近但不同的物种可能具有独特的基础生态位8。
在甲虫分类学中,与大多数(即使不是全部)动物类群类似,分类学中的初级分类或传统分类是物种鉴定的金标准9该方法依赖形态学特征对生物进行高阶分类学归类,并实现物种水平的鉴定10,11在许多情况下,家庭12 或甚至属水平的鉴定13 仅通过外部形态特征(如体型和特化适应结构)即可完成。然而,甲虫的种级鉴定通常需要 通过 雄性外生殖器或阳茎结构细节的比较14,15,16.
尽管这种基于形态学的分类学方法被广泛接受,但该物种发现途径受到多种“分类学障碍”的限制6,17,18,尤其是在无脊椎动物中更为突出19。在甲虫系统学中,利用分类学进行物种鉴定和(或)发现面临诸多挑战,包括熟练的分类学家数量有限20,21、早期描述简略且信息量不足22,以及获取文献和模式标本方面的后勤难题23。
当所研究的分类单元具有隐存性或极难被察觉时,这一问题尤为突出。隐存种是指同一属甚至不同属中的一类物种,由于种间表型差异极为细微,难以通过比较形态学方法进行明确区分24。一些隐存类群被称为“黑暗类群”25,26,27,它们不仅物种多样性极高,还普遍存在研究严重不足的问题。遗憾的是,这种情况在河岸带和水生甲虫中并不罕见8,28。若未能解决隐存性“黑暗类群”中的物种界定问题,将对生物多样性保护构成威胁,因为具有不同生态位和生态需求的生物可能被统一施加相同的保护措施,从而导致保护策略失效8,29。
鉴于这些挑战,在过去二十年中,“整合分类学”23,即在传统形态分类(α分类学)基础上结合其他证据来确立新物种的方法,已逐渐得到广泛应用30,31,32。在甲虫33,34,35,36,37,38,39及其他昆虫目40,41,42,43,44,45的整合分类学研究中,常用的一种辅助证据是DNA序列。其中,线粒体细胞色素c氧化酶亚基I(COI)基因因其具有适度保守性,正被广泛用于对多种动物进行DNA条形码鉴定46,47,48。尽管其658 bp长的5′端(COI-5′),也被称为“条形码片段”或“Folmer片段”,已被提议作为“诊断性”序列49,但由于存在一系列针对3′端的高效引物,723 bp长的COI-3′也正被广泛使用50,51,52。
尽管利用COI进行DNA条形码分析最初是为物种水平的鉴定而设计的,但其应用受限于以下因素 先验的 序列数据库53,例如 GenBank54 和 BOLD55,分子物种界定方法仅利用提供的DNA序列来确定物种界限,无需预先建立的参考数据库56. 根据输入文件的不同,分子物种划分方法可分为两类。第一类以序列比对作为输入,基于阈值的方法利用设定的截断值来确定物种间和物种内的遗传分歧度48,57,58. 作为基于距离的方法,这些方法依赖于昆虫观察到的种内最大距离为3%,该值被称为“条形码间隙”51,59一些基于阈值的方法包括在 TaxonDNA 中实现的 SpeciesIdentifier60,遗传距离 通过 Kimura双参数(K2P)61,以及通过自动划分组装物种(ASAP)62.
其次,以系统发育树作为输入文件,基于溯祖理论的方法通过确定谱系分支速率的变化来识别物种边界并划分独立的谱系63,64,66。这类方法通常依赖于系统发育物种概念67,包括泊松树过程(Poisson tree processes, PTP)68及其变体——多速率PTP69。无论采用何种方法,分子物种划分方法所形成的聚类统称为分子操作分类单元(molecular operational taxonomic units, MOTUs)。
本文介绍了一种用于生成分子操作分类单元(MOTUs)的整合分析流程,该流程包含三种基于阈值的方法和三种基于合并模型(coalescent-based)的方法(图1)。该流程在两个数据集上进行了测试:(1)河岸带Byrrhinus Motschulsky, 1858 鞘翅目昆虫的 COI-3′ 序列,该数据集包含两个先前通过整合分类学方法确立的已发表物种35,以及来自未描述物种的额外序列;(2)Anacaena Thomson, 1859 鞘翅目昆虫的 COI-5′ 序列,该数据集包含两个先前通过α分类学方法确立的已发表物种70,以及来自未描述物种的额外序列。

图1:生成MOTU的整合流程图。请点击此处查看该图的放大版本。
访问受限。请登录或开始试用以查看此内容。
该组合分析流程的起始材料为DNA序列,这些序列可从在线数据库下载(例如BOLD、GenBank),或通过Sanger测序获得71。表1和表2列出了本文所展示代表性结果所对应数据集的登录号。
1. 基于阈值方法的准备步骤:序列比对

图 2:MEGA X 中的序列比对。(A)手动校正后的最终 DNA 序列比对结果,以及(B)Byrrhinus 数据集的翻译后蛋白质序列比对结果,显示无插入、缺失和空位。请点击此处查看该图的放大版本。
2. 在物种鉴定器 1.8 的分类单元 DNA 模块中进行界定
注意:物种标识符可从 https://taxondna.sourceforge.net/ 下载

图3:Species Identifier 1.8 软件中 TaxonDNA 界面显示用户设定的3%聚类阈值用于物种划分。请点击此处查看该图的高清版本。
在 MEGA X 中使用 Kimura 2-参数(K2P)进行界定

图4:MEGA X 中显示 Kimura-2 参数法结果的窗口。 请点击此处查看该图的放大版本。
4. 通过自动划分组装物种进行界定(ASAP)

图5:ASAP中显示可能物种划分结果的窗口。 第一个结果(蓝色框)被接受,因其ASAP得分最低(红色框)且P值排名最高(黄色框)。请点击此处查看该图的放大版本。
5. 基于溯祖方法的准备步骤:系统发育树构建

图6:在 MEGA X 中进行系统树推断,作为基于溯祖理论分析的前期步骤。(A)用于确定最优核苷酸替换模型的设置。(B)对于 Byrrhinus 数据集,由于 BIC 和 AICc 值最低(红框),最优模型为 GTR + I。(C)在运行最大似然法构建系统树时,参数设置为 GTR(黄框)+ I(蓝框)。(D)所得原始 ML 树,将用于后续分析。请点击此处查看该图的高清版本。
6. 基于泊松树过程(Poisson Tree Processes, PTP)的物种划分
注意:以下步骤将产生 PTP-ML 和 PTP-BI 的结果。

图7:基于泊松树过程的分子物种界定。(A)显示输入参数的窗口。(B)PTP网络服务器中的结果窗口,其中包含物种界定结果的链接(红色框)。请点击此处查看该图的放大版本。
7. 基于多速率泊松树过程(mPTP)的物种界定

图 8:通过多速率泊松树过程进行分子物种界定。(A)显示 mPTP-read 关于已上传树信息的窗口。(B)mPTP 网络服务器中的结果窗口,包含物种界定结果的链接(红色框)。请点击此处查看此图的放大版本。
8. 生成 MOTU
访问受限。请登录或开始试用以查看此内容。
本文采用了一种分子物种界定复合分析流程,包括TaxonDNA、K2P、ASAP、两种变体的PTP以及mPTP方法,应用于水生及河岸带甲虫。六种方法各自所得结果称为“分子聚类”,而其综合结果则称为“分子操作分类单元”(MOTUs)75。
该分析流程在两个数据集上进行了测试。第一个数据集(表1,图9)包含24条来自菲律宾采集的河岸带Byrrhinus Motschulsky, 1858甲虫的COI-3'序列,以外群为两个未定种的Limnichus Latreille, 1829物种76。Byrrhinus序列包括两条来自先前已发表物种的序列35,这些物种基于形态学的α分类学并结合K2P划分方法建立,此外还包括来自五个不同省份的未知物种的序列。
访问受限。请登录或开始试用以查看此内容。
尽管DNA序列(尤其是COI49)可与在线公共数据库54,55,79,80中的参考序列22,53,78进行比对,以实现生物鉴定46,但该方法受限于已有的测序数据沉积情况。相比之下,分子物种界定方法可在无需预先建立序列数据库的前提下,确定物种边界并实现分子聚类56。鉴于分子物种界定方法可分为两类,在分析流程中存在两个关键步骤,即输入文件的准备,具体包括序列比对的构建和基因树的生成。
如本研究所示,当将组合...
访问受限。请登录或开始试用以查看此内容。
我们没有需要披露的内容。
我们感谢菲律宾渔业和水产资源局(BFAR)以及相关地方政府部门授予我们事先知情同意,以开展生物多样性采样工作。感谢 Hendrik Freitag 博士在分类学研究方面提供的指导。同时感谢审稿人和编辑们提出的富有建设性的宝贵意见。最后,我们谨向马尼拉雅典耀大学大学研究委员会重大项目资助(URC 2022-09)和开放获取出版资助、雅典耀科学与工程研究所 SC Johnson 学生科研与环境领导力基金,以及伦敦林奈学会与系统分类学会联合设立的 LinnéSys:系统分类研究基金,对本研究提供的经费支持表示最诚挚的感谢。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| ASAP | Spart 探索者 | https://bioinfo.mnhn.fr/abi/public/asap/asapweb.html | |
| 计算机 | 任意处理单元 | ||
| MEGA X | MEGA 软件 | https://www.megasoftware.net/ | |
| mPTP 网络服务器 | https://mptp.h-its.org/#/tree | ||
| PTP 网络服务器 | 海德堡理论研究学会(HITS) | https://species.h-its.org/ptp/ | |
| 物种鉴定器 | https://taxondna.sourceforge.net/ |
访问受限。请登录或开始试用以查看此内容。