雷山髭蟾蝌蚪形态学分析

摘要 为了解雷山髭蟾蝌蚪的形态学特征，2020年12月在雷公山国家级自然保护区内格头村、仙女塘、康利水厂附近小溪采集雷山髭蟾蝌蚪，并测量其体重和形态指标；采用多变量形态学分析方法对74只雷山髭蟾蝌蚪的体重、全长、眼前头长、体高、尾长等18项指标进行方差分析、主成分分析和聚类分析。结果表明：不同全长组蝌蚪的体重、吻肛距、全长等18项指标差异均极显著。在主成分分析中，有1个主成分的特征值大于1，其累计方差贡献率为76.351%；体重、全长、眼间距、体宽、尾长等指标对其有较大的影响。进行聚类分析的74只雷山髭蟾蝌蚪可聚为3个大类，全长0～50 mm组和全长101～150 mm组的雷山髭蟾蝌蚪关系较远；全长51～100 mm组的蝌蚪处于生长发育的过渡阶段，与全长0～50 mm组的雷山髭蟾蝌蚪有较近的关系。雷山髭蟾蝌蚪外部形态分析可为雷山髭蟾的生物学特征研究积累基础数据。

关键词 雷山髭蟾蝌蚪；形态指标；主成分分析；聚类分析

中图分类号 Q 954文献标识码 A

文章编号 0517-6611（2023）18-0096-07

doi：10.3969/j.issn.0517-6611.2023.18.022

Morphological Analysis of Vibrissaphora leishanensis Tadpoles

MENG Li-xia，WANG Zhao，YANG Jian et al

（Kaili University，Kaili，Guizhou 556011）

Abstract In order to understand the morphological characteristics of Vibrissaphora leishanensis tadpoles，V. leishanensis tadpoles were collected in the streams near Getou Village，Xiannvtang and Kangli Water Plant in Leigongshan National Nature Reserve in December 2020，and their body weight andmorphological indices were measured. Using the method of multivariable morphological analysis，18 indices such as body weight，total length，anterior head length，body height and tail length of 74 tadpoles of V. leishanensis were analyzed based on variance analysis，principal component analysis and cluster analysis. The results showed that there were extremely significant differences in 18 indices of tadpoles among different full-length groups，such as body weight，snout anal distance，total length. In principal component analysis，the eigenvalue of one principal component was greater than 1，and its cumulative variance contribution rate was 76.351%. The body weight，total length，eye spacing，body width，tail length and other morphological indices had greater impacts on it.The cluster analysis showed that 74 tadpoles of V. leishanensis were clustered into three categories： the relationship of V. leishanensis tadpoles betweentotal length 0-50 mm group andtotal length 101-150 mm group was farther. The tadpoles in total length 51-100 mm group were in the transitional stage of growth and development，which was closely related to the tadpoles in total length 0-50 mm group. The external morphological analysis of V.leishanensis tadpoles can accumulate the basic data for the study on the biological characteristics of V. leishanensis tadpoles.

Key words Vibrissaphora leishanensis tadpole;Morphological indices;Principal component analysis;Cluster analysis

髭蟾属（Vibrissaphora）隶属两栖纲（Amphibia）无尾目（Anura）角蟾科（Megophryidae），原隶属锄足蟾科（Pelobatidae），为中国特有属。刘承钊等［1］通过形态学特征将髭蟾属划入角蟾科，雄性在繁殖期上唇缘生有角质刺，俗称“角怪”“胡子蛙”。目前已知髭蟾属有峨眉髭蟾（V.boringii Liu）、崇安髭蟾指名亚种（V.liui liui Pope）、崇安髭蟾瑶山亚种（V.liui yaoshanensis Liu and Hu）、雷山髭蟾（V.leishanensis Liu and Hu）、哀牢髭蟾（V.ailaonica Yang，Chen and Ma 1983）、原髭蟾（V.promustache Rao，Wilkinson et Zhang）6个种（或亚种）［2］。雷山髭蟾仅分布于贵州雷公山，分布区域非常狭窄，主要栖息在海拔800 ～2 100 m、植被类型保存完好、溪流水源充足的常绿阔叶林、落叶阔叶林、常绿落叶阔叶混交林、针叶混交林以及针阔混交林下的山间小溪附近［3］。成体一般产卵在林间小溪扁平石块下；卵40 d左右孵出蝌蚪；蝌蚪完成变态过程要越过2个冬天，从卵发育成成体需要3年时间［4］。成蟾多爬行，不善跳跃，行动缓慢，容易遭天敌捕食；其蝌蚪期长，易遭天敌吞食，成活率低。近年来，雷山髭蟾的栖息环境遭到破坏，当地居民大量捕食雷山髭蟾成体、蝌蚪甚至卵，使得雷山髭蟾的数量呈现急剧减少的趋势，面临灭绝的危险。国家林业和草原局、农业农村部公布的2021年2月1日起实施的《国家重点保护野生动物名录》中，雷山髭蟾被列为国家二级重点保护野生动物［5］。

目前已发表有关雷山髭蟾研究论文20多篇，主要关于雷山髭蟾的生物学基本特征［6］、雷山髭蟾成体及幼体的生态环境［3，7-8］、卵群发育［9-10］、种群动态［11］、种群历史［12］、保护与发展［4，13］、繁殖期雷山髭蟾性二态分析及求偶鸣声［14］、微卫星DNA标记［15］、MHC Ⅰ类基因的多态性与进化分析［16］、线粒体全基因组［17］、营养成分分析及评价［18-19］等。目前尚未见到采用多变量形态学分析方法对雷山髭蟾蝌蚪形态特征进行分析的报道。雷山髭蟾是我国特有的珍稀保护动物，关于其生活史方面的研究很少。为了解雷山髭蟾蝌蚪的形态学特征，笔者采用多变量形态学分析方法对雷山髭蟾蝌蚪的全长、体高、尾长、尾高等18项指标进行测量，并对其进行方差分析、主成分分析和聚类分析，分析雷山髭蟾蝌蚪不同全长段间形态特征的差异，了解雷山髭蟾蝌蚪在野外自然状态下的生长发育情况及其适应水生生活的形态特征，旨在为雷山髭蟾蝌蚪的生物学特征积累基础数据。

1 研究方法

1.1 样本的采集

2020年12月在雷公山国家级自然保护区内格头村、仙女塘和康利水厂附近的溪流采集雷山髭蟾蝌蚪（图1），采集样本数74只。为减少该研究对雷山髭蟾蝌蚪的影响，对雷山髭蟾蝌蚪体重和形态指标进行测量后及时将其放回溪流，所测量蝌蚪均在自然条件下生长发育。

1.2 体重和形态指标测量

使用电子天平称量雷山髭蟾蝌蚪的体重，精确到0.01 g。使用电子数显卡尺对雷山髭蟾蝌蚪的形态指标进行测量，精确到0.01 mm，共测量全长、吻肛距、尾长、尾高、体宽、体高、上尾鳍高、下尾鳍高、吻出水孔距、眼前头长、吻鼻距、鼻间距、眼径、眼间距、口宽、吻肛距、尾基宽17项形态指标。依据组距式分组法，结合雷山髭蟾蝌蚪全长数据及样本数，以50 mm为组距，将全部样本分为3个组：0～50 mm组的样本有3只（1～3号样本），51～100 mm 组的样本有62只（4～65号样本），101～150 mm组的样本有9只（66～74号样本）。

1.3 数据处理

为消除雷山髭蟾蝌蚪个体间异速增长的影响，先使用Excel软件将雷山髭蟾蝌蚪各外部形态指标的原始数据转换为以10为底数的对数值，再使用SPSS 20.0软件进行方差分析、主成分分析和聚类分析。

2 结果与分析

2.1 雷山髭蟾蝌蚪的外部形态

雷山髭蟾蝌蚪体形肥大，皮肤光滑，呈红棕色，体尾有小黑点；头部呈淡红色且有少许褐色块斑，尾部具有淡黄色“V”形斑（图2a、c、e），尾鳍半透明，尾肌强壮，尾末端钝尖；吻端平圆，鼻孔近于吻两侧、靠近眼。眼在体前面的背侧；两眼相距较远。口为下位，具有发达的下颌肌及坚硬的锯齿状下颌鞘，唇齿发达。腹部呈灰白色，嘴呈开“W”形（图2b、d、f），盘曲的肠道等清晰可见，肛孔开口于腹部尾基中央。全长0～50 mm组雷山髭蟾蝌蚪的全长41.75～49.63 mm，平均全长45.66 mm；全长51～100 mm组雷山髭蟾蝌蚪全长53.55～99.19 mm，平均全长75.33 mm；全长101～150 mm组雷山髭蟾蝌蚪全长100.44～118.30 mm，平均全长106.73 mm。全长0～50 mm组雷山髭蟾蝌蚪尾部褐色斑点较为集中，且主要集中在基部，且后肢芽均未发育（图2a、b）。全长51～100 mm组和全长101～150 mm组雷山髭蟾蝌蚪褐色斑点分布比较分散，下尾鳍和上尾鳍均有褐色斑点分布（图2c～f）。全长51～100 mm组少部分蝌蚪后肢芽已生长，大部分蝌蚪未生长后肢芽。全长101～150 mm组的全部蝌蚪均具后肢芽，且生长程度不同。雷山髭蟾蝌蚪的四肢先后长出，后肢芽慢慢发育成后肢，后肢长出后再长出前肢。此次采集到的蝌蚪样本不管有无后肢芽，均无前肢。

2.2 雷山髭蟾蝌蚪体重和形态指标分析

2.2.1 雷山髭蟾蝌蚪体重和形态指标测量。雷山髭蟾蝌蚪的体重和形态指标见表1。全长0～50 mm组雷山髭蟾蝌蚪均无后肢，全长51～100 mm组雷山髭蟾蝌蚪部分有后肢，全长101～150 mm组雷山髭蟾蝌蚪均有后肢。

2.2.2 方差分析。

先将3组雷山髭蟾蝌蚪体重和形态指标的原始数据转换成以10为底数的对数值，再使用SPSS 20.0软件进行单因素方差分析，结果见表2。由表2可知，雷山髭蟾蝌蚪体重、吻肛距、全长、尾长、尾高、体宽、体高、上尾鳍高、下尾鳍高、吻出水孔距、眼前头长、吻鼻距、鼻间距、眼径、眼间距、口宽、后肢芽长、尾基宽的P值均小于0.01，表明3组雷山髭蟾蝌蚪的18项指标差异均极显著。

2.2.3 主成分分析。先将雷山髭蟾蝌蚪体重和形态指标的原始数据转换成以10为底数的对数值，再进行主成分分析，KMO值接近1（KMO=0.930），Bartlett球度检验的P值小于0.01（P=0.000），表明所测数据适合进行因子分析。从18项指标中提取特征值大于1的主成分有1个（图3），因为这个主成分的累计方差贡献率为76.351%，所以这个主成分保留了原有指标的大部分信息。通过主成分得分系数矩阵（表3）可知，对主成分影响较大的指标为体重（0.965）、全长（0.950）、眼间距（0.928）、体宽（0.907）、尾长（0.918）。

主成分的得分值=体重×0.070+吻肛距×0.068+全长×0.071+尾长×0.070+尾高×0.069+体宽×0.070+体高×0.056+上尾鳍高×0.064+下尾鳍高×0.068+吻出水孔距×0.063+眼前头长×0.053+吻鼻距×0.046+鼻间距×0.061+眼径×0.050+眼间距×0.070+口宽×0.061+后肢芽长×0.061+尾基宽×0.069。

2.2.4 聚类分析。

将雷山髭蟾蝌蚪体重和形态指标的原始数据转换成以10为底数的对数值，再进行聚类分析。采用单个样本进行聚类分析，使用欧氏距离法对体重、全长、眼间距、体宽、尾长等18个指标进行聚类分析，聚类树状图见图4。74只雷山髭蟾蝌蚪共进行73步聚类。以标线5为基准，可分成五大类：5号样本单独为一大类（规定为第二类），1、2、3、4、6、7号样本为一大类（第三类），51、54、55、56为一大类（第四类），57、58、59、60、61、62、63、64、65、66、67、68、69、70、71、72、73、74号样本为一大类（第五类），其他样本为第一类，第一类样本最多（45个）；第一类和第二类聚为一支，再与第三类聚在一起；第四类和第五类聚在一起；前三类和后两类的距离最远。全长0～50 mm组（1～3号样本）与全长101～150 mm组（66～74号样本）的雷山髭蟾蝌蚪关系最远；全长51～100 mm组（4～65号样本）的雷山髭蟾蝌蚪处于生长发育的过渡阶段，与全长0～50 mm组雷山髭蟾蝌蚪有较近的关系，所以较多蝌蚪与该组蝌蚪聚为一类，而较少蝌蚪与全长101～150 mm组雷山髭蟾蝌蚪聚为一类。