中国河北晚侏罗世蝾螈亚目（两栖纲，有尾目）基干属种干沟青龙螈（Qinglongtriton gangouensis）骨骼高精度μCT三维数据集

贾佳，高克勤

1.北京大学，地球与空间科学学院，北京 100871

2.山东科技大学，山东省沉积成矿作用与沉积矿产重点实验室，山东青岛266590

3.中国科学院南京地质古生物研究所，中国科学院资源地层学与古地理学重点实验室，南京 210008

数据库（集）基本信息简介

数据库（集）名称 蝾螈亚目（两栖类，有尾目）基干属种干沟青龙螈骨骼高精度μCT扫描源数据集数据作者 贾佳、高克勤数据通信作者 贾佳（jia_jia@pku.edu.cn）数据时间范围青龙螈化石层位属于晚侏罗世髫髻山组，U-Pb同位素测年显示距今约（160.889±0.069）Ma～（160.254±0.045）Ma[1]；化石标本采集于2009年夏，并于2015年春扫描。

?

引 言

在传统古生物学研究中，当把一件标本从野外运输到实验室内，研究者为了研究其形态特征，通常会在显微镜下首先对其进行机械修理或者酸处理，以去除其围岩从而展露出化石的真实面貌[2]。然而这些修理手段很难完整地揭示出化石本身的三维形态特征，反而常常会对标本造成损坏，带来不可挽回的损失。如若需要观察化石的内部结构，在标本充足的前提下，研究者通常选择切片技术（磨片法、揭片法）来获取标本的部分截面，从而进行观察[3]。但是这些方法并不能实现完全的连续切片，致使无法完整地观察化石的内部结构。

高精度计算机断层扫描（high-resolution X-ray tomography，high-resolution micro-CT或μCT）是一种利用高能 X光射线在短时间内获得物体内部三维结构的技术。与上述传统手段相比，该技术以其高精度、高效率和无损探测等优势迅速获得了古生物研究人员的青睐，并广泛应用在许多门类的关键过渡属种中，如肉鳍鱼类Eusthenopteron[4-5]，低等四足动物Ichthyostega[6]，原始锄足蟾类Prospea[7]，原始蛇类Dinilysia[8]，原始翼龙类Liaodactylus[9]，霸王龙类恐龙Alioramus[10]，早期鸟类孔子鸟[11]以及原始哺乳类Vilevolodon[12]等。基于μCT扫描所得到的虚拟模型极大地方便了研究者对标本的测量和观察，同时也扩展了对生物本身的结构功能学和生物力学的研究，如对低等四足动物Ichthyostega行走能力的研究[6]。除了辅助科学研究以外，μCT扫描所生成的虚拟模型还可以通过3D打印机制成实体模型，用于博物馆展览、公众科普等。

本文所提供的原始蝾螈类干沟青龙螈（Qinglongtritongangouensis）的骨骼μCT三维源数据集便是μCT的众多应用实例之一。蝾螈类（有尾目）是一类有尾两栖动物，包括隐鳃鲵亚目（2科，11属，70种）、蝾螈亚目（7科，55属，633种）和鳗螈亚目（1科，2属，4种），它们与蛙类和蚓螈类共同构成滑体两栖动物亚纲[13]。蝾螈类是一支古老的演化支系，分子生物学数据显示该支系起源于约262.5～357.0 Ma之间（如[14-15]）。自上世纪90年代末以来，中国河北北部、辽宁西部和内蒙古自治区东南部的中生代地层出产了数以万计的蝾螈类标本，为了解蝾螈类在中生代时期的演化提供了珍贵的化石材料。这些化石标本不仅数目众多，而且保存完整，同时分类多样性也非常丰富。目前已识别出隐鳃鲵亚目的8属10种和蝾螈亚目的2属2种[16-18]。本文对青龙螈多件标本的高精度断层扫描和3D打印是在我国化石蝾螈类研究中的首次尝试[17]。与以往研究相比，该项研究提供了更为详尽的特征描述，揭示了这一原始蝾螈物种的多个独特的骨骼结构，为更好地理解蝾螈亚目的早期特征演化提供了信息[17]。

共享高质量的μCT扫描源数据将为获取某些重要标本的虚拟模型提供便利，促进科研合作。我们在研究过程中发现多个可用共享平台，如Dryad、DigiMorph等已经提供了少数标本的源数据共享和无偿下载功能，但是共享力度仍然不足。因此，本文提供了我国晚侏罗世髫髻山组产出的干沟青龙螈的正型标本（PKUP V0226）和2件归入标本（PKUP V0228、PKUP V0254）的μCT扫描源数据，正型标本的上半身骨骼可用于3D打印的stl文件以及3件标本的复原图和相关视频。

1 数据采集和处理方法

1.1 获取数据

在进行扫描以前，首先需要对标本进行初步的镜下机械修理，使得标本能够更好地固定在扫描仪内部的载物台上，进而获得较高的分辨率。2009年的野外工作共发现干沟青龙螈标本46件[17]。所有标本的修理工作在北京大学地球与空间科学学院完成。对正型标本PKUP V0226和归入标本PKUP V0228、PKUP V0254的扫描工作在中国地质大学（北京）科学研究院工业微焦CT实验室的扫描仪（Nikon XT H 320 LC）上操作完成。后续对标本进行的图像分割（segmentation）和三维重建（3D reconstruction）在该实验室工作站上利用VG Studio Max 2.2（Volume Graphics，Heidelberg，Germany）软件完成。尽管该软件应用广泛，但并非仅此一类软件可供处理图像堆栈（image stack），读者可参照[19]关于其他付费或免费软件的相关介绍。

1.1.1 标本预处理和μCT扫描

干沟青龙螈的正型标本PKUP V0226及其他2件归入标本PKUP V0228、PKUP V0254在野外采集过程中断裂为大小不一的几块。在解剖显微镜下，使用碳钢针小心去除每一块标本覆盖在骨骼上的围岩。待骨架暴露完全后，使用锯条或小型齿轮机沿骨架的长轴将多余的围岩去除，以减少标本整体的宽度，从而提高扫描成像效果。将清漆与丙酮溶液的混合液用于化石标本加固，以免骨骼碎裂脱落；之后，使用苏打粉与 2-氰基丙烯酸乙酯将同一标本碎裂的部分粘连在一起。

将粘连好的标本固定在扫描仪的载物台上，并保持整体骨骼的长轴方向与载物台垂直。通过电脑调整载物台与X光射线源的距离，使标本在X光接收器上的投影尽量大而且完整；旋转载物台 360°以确保化石标本在每个角度上的投影都出现在X光接收器上。由于不同标本厚度不一，调节电流与电压以强化X光射线的穿透力，从而提高化石与围岩之间的对比度。

1.1.2 数据生成

在扫描青龙螈标本的过程中，载物台每旋转0.11457670273711°，X光接收器采集一张包含化石标本灰度值（gray value）的二维X光图像（tiff格式）。扫描过程中，每件标本共生成3 142张tiff格式二维图像和一个xtekct格式的参数文本文件。此后利用Nikon扫描仪自带的3D Pro软件读取3 142张tiff文件，通过矫正在扫描过程中可能出现的标本晃动或X光射线硬化等造成的标本边缘模糊现象，进而生成垂直于标本长轴方向（Z轴）的1 998张16位 tiff格式的图像堆栈，即本文所提供的扫描源数据（PKUPV0226rawdata.zip，PKUPV0228rawdata.zip，PKUPV0254rawdata.zip），流程见图1。

图1 干沟青龙螈标本修理、扫描及数据处理流程示意图

1.2 处理数据

扫描数据处理过程中，将1 998张16位tiff格式图像堆栈导入VG Studio Max 2.2软件，通过弱化背景和强化化石骨骼的灰度值，从而提高二者在二维截屏面上的对比度。完成源数据的导入之后，VG Studio Max软件分为4个窗口展示该标本，其中3个窗口为相互垂直的3个截平面（X-Y轴，X-Z轴和Y-Z轴），另一窗口显示标本的三维画面。在某一截平面上，利用软件的画笔工具将特定的骨骼轮廓一帧一帧地画出，生成感兴趣区（region of interest，简称ROI；图2）。将ROI提取即可生成某一骨骼的三维虚拟模型。倘若需要对某一结构进行3D打印，只需将该三维模型进行表面测定，导出stl格式的文件即可，例如本文提供的青龙螈正型标本上半身骨骼模型的stl文件（PKUPV0226stl.zip）。此外，本文将青龙螈的3件标本分别进行了二维截面的图像分割，揭示出了化石标本的完整结构，并将骨骼的背、腹面视图保存为高分辨率的tiff格式图片和avi格式的视频（pics.zip和videos.zip）。

图2 干沟青龙螈正型标本头骨渲染图及感兴趣区（ROI）的建立

2 数据样本描述

干沟青龙螈的正型标本（PKUP V0226）和2件归入标本（PKUP V0228、PKUP V0254）均为成年个体（图3～4）。其中PKUP V0226的头长（前颌骨至外枕骨后缘距离）为45.11 mm，宽（头骨与两侧下颌关节之间最大距离）为40.35 mm，吻臀距（前颌骨至坐骨后缘距离）为155.71 mm；PKUP V0228的头长40.79 mm，宽35.78 mm，吻臀距148.47 mm；PKUP V0254的头长46.18 mm，宽44.23 mm，该标本躯干后部保存不完整无法获得吻臀距。其中，正型标本为腹面向上保存，其余两件为背面向上保存（图 3）。通过对标本的观察发现干沟青龙螈具有双头肋骨（bicapitate ribs），而且下颌中的隅骨（angular）与前关节骨（prearticular）愈合，属于有尾目中蝾螈亚目的典型特征。分支分析研究发现该属种与此前发现的建平北燕螈（Beiyanerpetonjianpingensis）[20]构成姐妹群，共同位于蝾螈亚目的最基干位置[17]。

青龙螈的系统分类信息如下：

两栖纲 Amphibia Linnaeus, 1758

滑体两栖亚纲 Lissamphibia Haeckel, 1866

有尾超目 Caudata Scopoli, 1777

有尾目 Urodela Duméril, 1806

蝾螈亚目 Salamandroidea Dunn, 1922

科未定 Family Incertae Sedis

干沟青龙螈QinglongtritongangouensisJia and Gao, 2016

图3 干沟青龙螈3件标本照片

对青龙螈 3件标本的高精度断层扫描识别出了传统镜下机械修理所不能揭露的结构特征。举例来说（图2），紫色和橙色部分分别为头骨背面的额骨（frontal）和顶骨（parietal），墨绿色部分为头骨腭面的副蝶骨（parasphenoid），通过对 μCT扫描的多个横截面观察，并未发现眶蝶骨（orbitosphenoid）。在大部分现在存活的蝾螈类中，眶蝶骨作为脑颅的一部分，起到保护软体脑颅和支撑头骨顶部的作用。该骨仅在较为进步的洞螈科（Proteidae）和无肺螈科（Plethodontidae）Eurycea属的个别种类中缺失。眶蝶骨在青龙螈所有已扫描标本中的缺失现象说明这一特征在蝾螈亚目演化初期就已有出现，而在洞螈科以及Eurycea中的缺失很有可能是由于同塑演化（homoplasy）导致的。

图4 干沟青龙螈归入标本PKUP V0228的μCT重建图

3 数据质量控制和评估

如前所述，青龙螈骨骼高精度μCT扫描源数据集是使用中国地质大学（北京）科学研究院工业微焦CT实验室的扫描仪（Nikon XT H 320 LC）完成。该扫描仪的X光接收器的分辨率为2000×2000像素。对不同标本的扫描均未使用任何滤片，其电压、电流、分辨率汇总如表1，其余参数见数据集中的xtekct参数文本文件。

表1 干沟青龙螈（Qinglongtriton gangouensis）标本高精度μCT扫描参数表

干沟青龙螈正型标本的上半身骨骼模型stl格式文件，是研究者基于对化石材料的图像分割和光滑处理后生成的。对每一块骨骼的原始ROI使用光滑力度（smooth strength）为1，深度（depth）为8位的光滑处理，之后提取成体积文件并做表面测定再导出成stl格式文件。

3件标本的视频文件均由VG Studio Max 2.2的动画（Animation）选项中的关键帧（key frame）完成，视频导出时无任何压缩处理。

4 数据价值

与传统的连续切片技术以及耗时耗力的化石修理相比，高精度X光扫描可以无损而且高效揭示出化石标本内部结构，极大方便了科研人员对标本的观察。对青龙螈骨骼的高精度μCT扫描是我国中生代蝾螈类研究中的首次尝试。包含扫描源数据和成果图片及视频的数据集的共享将为研究干沟青龙螈的骨骼特征提供完整的信息，为日后μCT扫描在蝾螈类研究中的应用提供参照。干沟青龙螈正型标本上半身骨骼模型stl文件的共享将为获取3D打印青龙螈的实体模型提供便利。