四川动物  2021, Vol. 40 Issue (1): 59-65

扩展功能

文章信息

程彦林, 刘京, 陈龙, 雷诗融, 邹家伟, 张霞, 李仕泽
CHENG Yanlin, LIU Jing, CHEN Long, LEI Shirong, ZOU Jiawei, ZHANG Xia, LI Shize
贵州桐梓发现大别山林蛙
Identification of Rana dabieshanensis in Tongzi, Guizhou Province
四川动物, 2021, 40(1): 59-65
Sichuan Journal of Zoology, 2021, 40(1): 59-65
10.11984/j.issn.1000-7083.20200206

文章历史

收稿日期: 2020-06-02
接受日期: 2020-09-30
贵州桐梓发现大别山林蛙
程彦林 , 刘京 , 陈龙 , 雷诗融 , 邹家伟 , 张霞 , 李仕泽 *     
茅台学院资源环境系, 贵州仁怀 564500
摘要:2019年7月在贵州桐梓黄莲自然保护区(107°2'37.31″E,28°29'31.95″N,海拔1 666 m)采集到2号两栖动物标本,经形态特征比较,与大别山林蛙 Rana dabieshanensis 相似;基于线粒体 16S rRNA 构建的蛙属 Rana 部分物种贝叶斯系统发育树显示,此次采集的林蛙标本与大别山林蛙聚为一支,且具有较高的支持率(0.98);基于Kimura双参数模型估算,该标本与大别山林蛙模式产地标本间的遗传距离为0.6%~0.8%,远小于蛙属物种间的遗传距离。综合形态学和系统发育比较,确定此次采集到的标本为无尾目Anura蛙科Ranidae蛙属的大别山林蛙,系贵州省两栖动物分布新记录。
关键词大别山林蛙    新记录    贵州省    
Identification of Rana dabieshanensis in Tongzi, Guizhou Province
CHENG Yanlin , LIU Jing , CHEN Long , LEI Shirong , ZOU Jiawei , ZHANG Xia , LI Shize *     
Department of Resources and Environment, Moutai Institute, Renhuai, Guizhou Province 564500, China
Abstract: During the amphibian survey in the Huanglian Nature Reserve, Tongzi County, Guizhou Province in July 2019, 2 amphibian specimens were collected (107°2'37.31″E, 28°29'31.95″N, altitude 1 666 m). We compared the specimens with other species of Rana by using morphological and molecular methods. The specimens were morphologically similar with R. dabieshanensis , and the 16S rRNA gene based on Bayesian Inference tree showed that the 2 specimens were clustered with the R. dabieshanensis with a high support (0.98). The genetic distance between the 2 specimens and R. dabieshanensis from the type locality was 0.6%-0.8%, which was much lower than that between Rana species. Therefore, the specimens collected in this study were identified as R. dabieshanensis , and this is the new record of this species in Guizhou Province, China.
Keywords: Rana dabieshanensis    new record    Guizhou Province    

蛙属 Rana 广泛分布于欧亚大陆和美洲,目前该属已报道51种(Frost,2020),其中,中国记录分布25种,贵州省记录分布6种,分别为昭觉林蛙 R. chaochiaoensis 、徂徕林蛙 R. culaiensis 、峨眉林蛙 R. omeimontis 、威宁蛙 R. weiningensis 、镇海林蛙 R. zhenhaiensis 和寒露林蛙 R. hanluica (中国两栖类,2020)。昭觉林蛙分布在贵州西部(伍律等,1986),徂徕林蛙分布在贵州黄平等地(肖宁等,2019),峨眉林蛙分布在贵州东部和北部(费梁等,2012),威宁蛙分布在贵州西部威宁、水城等地(刘承钊等,1962),镇海林蛙分布在贵州梵净山和雷公山(张雷等,2010陈继军等,2014魏刚等,2017),寒露林蛙分布在贵州梵净山和雷公山(陈继军等,2017魏刚等,2017)。

2019年7月,在贵州桐梓黄莲自然保护区进行两栖动物调查时,于保护区内(107°2′37.31″E,28°29′31.95″N,海拔1 666 m)采集到蛙属成体标本2号(表 1),经形态学比较和系统发育分析,鉴定为大别山林蛙 R. dabieshanensis ,为贵州省两栖动物新记录。本文对其形态特征及生态习性进行描述,并初步探讨其形态特征的地理变异,为蛙属物种的系统学和谱系地理学研究提供了新信息。

表 1 样品信息表 Table 1 Information of samples used in this study
物种Species 标本号ID 采集地Locality 16S rRNA GenBank登录号
16S rRNA of GenBank accession number
本文标本 TZ20190714018
TZ20190714019
中国贵州桐梓 MT911420
MT911421
大别山林蛙 R. dabieshanensis AHU2016R001 中国安徽大别山 MF172963
黑龙江林蛙 R. amurensis MSUZP-SLK-RUS49 俄罗斯 KX269203
田野蛙 R. arvalis MSUZP-SLK-MKR21 俄罗斯 KX269197
中亚林蛙 R. asiatica KIZ-XJ0251 中国新疆47团 KX021945
昭觉林蛙 R. chaochiaoensis SCUM0405170CJ 中国四川昭觉 KX269192
中国林蛙 R. chensinensis KIZ-RD05SHX01 中国陕西户县 KX269186
徂徕林蛙 R. culaiensis KIZ-SD080501 中国山东徂徕山 KX021986
东北林蛙 R. dybowskii MSUZP-IVM-1d 俄罗斯 KX021949
寒露林蛙 R. hanluica SYNU07100490 中国湖南双牌 HQ228158
桓仁林蛙 R. huanrenensis 中国辽宁桓仁 GU225916
借母溪林蛙 R. jiemuxiensis KIZHUN0708013 中国湖南借母溪 KX269221
高原林蛙 R. kukunoris KIZCJ06102001 中国青海青海湖 KX021947
韩国林蛙 R. coreana KIZ-HUI040001 中国山东昆嵛山 KX269201
长肢林蛙 R. longicrus NMNS15022 中国台湾 KX269189
峨眉林蛙 R. omeimontis SCUM0405196CJ 中国四川洪雅 KX021946
镇海林蛙 R. zhenhaiensis KIZ0803271 中国浙江镇海 KX269218
梭德氏蛙 R. sauteri SCUM0405175CJ 中国台湾高雄 KX269204
明全蛙 R. zhengi SCUM0405190CJ 中国四川洪雅 KX269206
胫腺蛙 R. shuchinae CIB-HUI040009 中国四川昭觉 KX269210
威宁蛙 R. weiningensis SCUM0405171 中国贵州威宁 KX269217
黑斑侧褶蛙 Pelophylax nigromaculatus SCUM045199CJ 中国四川洪雅 KX269216
宝兴树蛙 Zhangixalus dugritei 中国四川宝兴 EU215541
1 材料与方法 1.1 材料

2号标本取肌肉组织保存于95%乙醇中,标本浸泡于10%福尔马林溶液中,保存于茅台学院生物学标本室。

1.2 形态鉴定与测量

对2号标本使用游标卡尺(200 mm,精度0.01 mm;无锡锡工量具有限公司)进行测量,测量依据《中国两栖动物检索及图鉴》(费梁等,2005), 并参照《中国两栖动物及其分布彩色图鉴》(费梁等,2012)进行形态鉴定。

1.3 分子鉴定

以高盐法提取肌肉组织的DNA。利用引物P7和P8(Simon et al., 1994)PCR扩增 16S rRNA 基因部分片段(528 bp)并测序。PCR扩增条件和步骤参照Chen等(2017)。扩增产物经电泳检测后,大小符合要求的产物交上海生工生物公司测序。所得序列上传到GenBank,登录号为MT911420和MT911421。

在GenBank中下载蛙属部分物种和蛙科Ranidae部分相近物种的 16S rRNA 基因,结合本次采集标本所测序列,进行分子系统发育分析。以宝兴树蛙 Zhangixalus dugritei 、黑斑侧褶蛙 Pelophylax nigromaculatus 为外群(Wang et al., 2017)。所有序列在Clustal X v1.8(Thompson et al., 1997)中进行比对,采用贝叶斯法(Bayesian-Inference,BI)重建系统发育关系。建树之前,利用jModeltest v2(David,2008)基于AIC标准选择最适核苷酸替代模型。贝叶斯法运算中马尔科夫链的蒙特卡罗方法(Markov chain,MonteCarlo process)设置为4条链(3条冷链1条热链)同时运行,共运行10 000 000代,每100代抽样一次,共得到100 001个样本及其相对的似然值和模型参数。在运行500 000代后,似然值趋于稳定,舍弃最初的500 000代(5 000个样本),用剩余的样本构建一致树,并获得系统树的支系结构和各支系的后验概率。用MEGA v6.0(Tamura et al., 2013)中的Kimura双参数模型计算表 1中所列蛙属物种间的遗传距离。

2 结果 2.1 形态描述

本次采集标本为成体雄性(TZ20160714018、TZ20160714019),体型较大,体长53.9 mm和54.8 mm,其他部位量度见表 2。吻端钝尖,头长略大于头宽,吻棱明显,瞳孔椭圆、横置,鼓膜明显,呈圆形,舌后端有明显缺刻。前臂粗壮,手指纤细,指端钝圆不膨大,指长序3>4>1≈2,内掌突椭圆形,小而明显,无外掌突;后肢细长,左右跟部重叠明显,胫跗关节贴体前伸达眼至鼻孔之间,趾长序4>3>5>2>1,趾端无沟,趾间具2/3蹼,趾关节下瘤明显:第四趾3个,第二、三、五趾各2个,第一趾1个。背部皮肤光滑,呈金色,背部有黑色斑点,背侧褶明显,从眼后角到腹股沟。生活时前肢、大腿、胫骨背侧面有灰褐色横纹,雄性婚垫灰褐色(图 1图 2)。比较此次所采集的林蛙标本和其他蛙属物种的形态鉴定特征,发现本次采集的标本形态特征与大别山林蛙模式产地标本相符。

表 2 贵州桐梓大别山林蛙成体标本形态量度 Table 2 Morphometric measurements of adult Rana dabieshanensis in Tongzi, Guizhou Province
形态量度 雄性Male ( n =2)
TZ20160714018 TZ20160714019
头体长/mm 53.9 54.8
头长/mm 18.0 18.0
头宽/mm 16.6 17.2
吻长/mm 8.2 8.5
鼻间距/mm 4.3 3.6
眼径/mm 5.4 5.1
眼间距/mm 2.7 3.6
上眼睑宽/mm 3.8 4.1
鼓膜径/mm 5.0 5.1
前臂及手长/mm 24.0 24.7
前臂宽/mm 5.4 5.2
后肢长/mm 100.9 103.9
股长/mm 28.5 29.4
胫长/mm 31.6 33.4
胫宽/mm 6.6 7.3
跗足长/mm 43.7 46.6
足长/mm 34.9 34.6

图 1 贵州桐梓大别山林蛙生活照 Fig. 1 Living individual of Rana dabieshanensis in Tongzi, Guizhou Province a.背面观, b.腹面观 a. dorsal view, b. ventral view

图 2 贵州桐梓大别山林蛙浸泡标本照 Fig. 2 Fixed specimen of Rana dabieshanensis in Tongzi, Guizhou Province a.侧面观, b.腹面观, c.手背面观, d.手正面观, e.脚腹面观 a. lateral view, b. ventral view, c. dorsal view of hand, d. facade view of hand, e. ventral view of foot
2.2 系统发育关系和遗传距离

PCR扩增获得线粒体 16S rRNA 基因片段长度为520 bp。T、C、A、G含量分别是25.9%、23.8%、29.7%、20.7%,A+T含量高于G+C。本文所用最适核苷酸替代模型为GTR+G+I。

基于贝叶斯算法构建系统发育树显示,采自桐梓的林蛙标本与大别山林蛙模式产地(安徽大别山)的标本聚为一支,且有较高的支持率(0.98)(图 3)。基于Kimura双参数模型估算本研究中所采用的蛙属部分物种的遗传距离(表 3)为1.2%~14.8%,平均遗传距离为8.1%。桐梓所采的林蛙标本与大别山林蛙模式产地标本间的遗传距离为0.6%~0.8%,远小于蛙属物种间的遗传距离。通过形态比较,在桐梓采集的林蛙标本与大别山林蛙原始文献(Wang et al., 2017)的描述吻合,系统发育分析也支持采自桐梓的林蛙标本与大别山林蛙关系较近。

图 3 基于 16S rRNA 基因序列构建的蛙属部分物种贝叶斯系统发育树 Fig. 3 Bayesian-Inference tree of Rana genus based on 16S rRNA gene sequences 枝上的数字代表重抽样支持率, TZ20190714018和TZ20190714019为本研究采集的样品, 其余物种序列信息见表 1 Numbers beside nodes indicate bootstrap values, sequences of TZ20190714018 and TZ20190714019 were from the Huanglian Nature Reserve, Tongzi, and the information of other species sequences were shown in table 1

表 3 基于16 S rRNA 基因序列估算的蛙属物种间Kamura双参数遗传距离 Table 3 Genetic distance of some Rana species using 16 S rRNA gene based Kamura 2-parameter model
物种Species 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
1 大别山林蛙 R. dabieshanensis
2 大别山林蛙 0.002
3 大别山林蛙 0.008 0.006
4 徂徕林蛙 R. culaiensis 0.012 0.010 0.017
5 镇海林蛙 R.zhenhaiensis 0.012 0.010 0.017 0.004
6 借母溪林蛙 R. jiemuxiensis 0.015 0.013 0.019 0.010 0.010
7 长肢林蛙 R. longicrus 0.017 0.015 0.021 0.008 0.008 0.015
8 寒露林蛙 R. hanluica 0.025 0.023 0.025 0.023 0.023 0.027 0.027
9 峨眉林蛙 R. omeimontis 0.038 0.036 0.038 0.036 0.036 0.038 0.041 0.041
10 昭觉林蛙 R. chaochiaoensis 0.038 0.036 0.036 0.036 0.036 0.043 0.040 0.038 0.023
11 桓仁林蛙 R. huanrenensis 0.047 0.045 0.047 0.043 0.043 0.038 0.045 0.049 0.056 0.058
12 黑龙江林蛙 R. amurensis 0.049 0.047 0.045 0.045 0.045 0.040 0.047 0.052 0.063 0.063 0.015
13 中亚林蛙 R. asiatica 0.052 0.049 0.052 0.047 0.047 0.047 0.052 0.063 0.063 0.060 0.043 0.050
14 田野蛙 R. arvalis 0.056 0.054 0.056 0.052 0.052 0.056 0.052 0.058 0.056 0.056 0.059 0.056 0.034
15 明全蛙 R. zhengi 0.065 0.063 0.065 0.060 0.060 0.056 0.060 0.067 0.063 0.074 0.063 0.065 0.060 0.045
16 胫腺蛙 R. shuchinae 0.074 0.072 0.079 0.065 0.065 0.061 0.065 0.072 0.072 0.074 0.072 0.074 0.063 0.061 0.058
17 梭德氏蛙 R. sauteri 0.074 0.072 0.079 0.067 0.067 0.063 0.067 0.074 0.065 0.074 0.058 0.065 0.043 0.054 0.067 0.072
18 东北林蛙 R. dybowskii 0.082 0.079 0.082 0.072 0.072 0.072 0.077 0.074 0.082 0.081 0.072 0.070 0.059 0.059 0.074 0.079 0.070
19 高原林蛙 R. kukunoris 0.086 0.084 0.091 0.081 0.081 0.082 0.081 0.089 0.074 0.086 0.067 0.077 0.063 0.045 0.060 0.088 0.065 0.050
20 中国林蛙 R. chensinensis 0.086 0.084 0.091 0.081 0.081 0.086 0.081 0.089 0.084 0.081 0.074 0.082 0.058 0.045 0.074 0.093 0.072 0.054 0.021
21 韩国林蛙 R. coreana 0.086 0.084 0.081 0.081 0.081 0.077 0.077 0.088 0.093 0.095 0.056 0.045 0.065 0.061 0.074 0.086 0.079 0.074 0.067 0.072
22 威宁蛙 R. weiningensis 0.148 0.146 0.146 0.148 0.148 0.149 0.143 0.153 0.135 0.148 0.142 0.140 0.143 0.140 0.130 0.135 0.134 0.148 0.143 0.164 0.132
注:1、2.大别山林蛙(桐梓黄莲), 3.大别山林蛙(安徽大别山)
Notes: 1, 2. R. dabieshanensis from the Huanglian Nature Reserve, Tongzi, 3. R. dabieshanensis from the Dabie Mountains, Anhui Province

因此,结合形态比较与分子系统学分析结果,确定采集于桐梓的林蛙标本为大别山林蛙,为贵州省两栖动物分布新记录。

2.3 生境描述

本次采集的大别山林蛙标本分布在桐梓黄莲自然保护区海拔1 650~1 700 m溪沟边树林茂密的草丛中,标本发现时在草丛边的道路上(图 4)。与其同域分布的两栖类有棘腹蛙 Quasipaa boulengeri 、峨眉树蛙 Zhangixalus omeimontis 、经甫树蛙 Zhangixalus chenfui 和斑腿泛树蛙 Polypedates megacephalus 等。

图 4 贵州桐梓黄莲自然保护区大别山林蛙生境 Fig. 4 Habitat of Rana dabieshanensis in the Huanglian Nature Reserve, Tongzi, Guizhou Province
3 讨论

蛙属物种在中国广泛分布,大多数物种间形态较为相似,因此在野外很难对其进行鉴定(费梁等,2012)。近年来,随着分子生物学的发展,采用形态比较结合分子系统学分析对物种进行鉴定的方法已广泛应用。

目前贵州省记录的蛙属物种6种,即昭觉林蛙、徂徕林蛙、峨眉林蛙、威宁蛙、镇海林蛙和寒露林蛙(中国两栖类,2020)。本文通过GenBank下载贵州已有报道的林蛙及其近缘种的 16S rRNA 序列,构建系统发育树,计算遗传距离,发现本文采集的标本和贵州省有记载分布的昭觉林蛙、徂徕林蛙、峨眉林蛙、威宁蛙、镇海林蛙和寒露林蛙的遗传距离分别为3.8%、1.2%、3.8%、14.8%、1.2%和2.5%,和大别山林蛙模式产地标本的遗传距离仅为0.6%和0.8%,仅较镇海林蛙与徂徕林蛙间的遗传距离大(0.4%),远小于其他蛙属物种的遗传距离。但有研究认为徂徕林蛙与镇海林蛙处于物种分化过程中的并系阶段,其物种分化还未完全完成(杨宝田等,2013)。因此,本文遗传距离比较不将徂徕林蛙和镇海林蛙纳入。

采自桐梓的林蛙标本与贵州已记载的6种蛙属物种在形态上也存在明显差异。头体长53.9~54.8 mm,可区别于威宁蛙(35 mm左右);无雄性线,区别于峨眉林蛙(有雄性线)和昭觉林蛙(有粉红色雄性线);颞褶明显,区别于徂徕林蛙(细弱)、镇海林蛙(弱)、寒露林蛙(弱)和峨眉林蛙(不显);趾端圆,趾端无沟,区别于威宁蛙的趾端吸盘小,具腹侧沟;后肢贴体前伸达眼至鼻孔之间,区别于徂徕林蛙、镇海林蛙、峨眉林蛙的后肢贴体前伸达鼻孔前方,威宁蛙的达鼻孔和吻端之间,寒露林蛙和昭觉林蛙的超过吻端。本次采集的大别山林蛙和模式产地标本在形态上也存在一些差异,桐梓的指长序为3>4>1≈2,大别山的为3>1>4>2。但本次样本数量有限,上述差异是否为适应不同的环境而导致的形态变异,还有待采集标本进行论证。

大别山林蛙为中国特有种,此前报道在我国的分布仅限于安徽岳西县鹞落坪自然保护区,海拔1 150 m左右的山区(Wang et al., 2017)。此次在贵州省桐梓县黄莲自然保护区内(海拔1 666 m)发现大别山林蛙,丰富了其分布信息,以及拓宽了对该物种分布海拔的认识。

参考文献
陈继军, 古定豪, 侯德华, 等. 2017. 贵州发现寒露林蛙[J]. 动物学杂志, 52(1): 155-157.
陈继军, 李仕泽, 吕敬才, 等. 2014. 贵州省雷公山国家级自然保护区两栖类新记录镇海林蛙[J]. 安徽农业科学, 42(28): 9797. DOI:10.3969/j.issn.0517-6611.2014.28.054
费梁, 叶昌媛, 黄永昭, 等. 2005. 中国两栖动物检索及图解[M]. 成都: 四川科学技术出版社.
费梁, 叶昌媛, 江建平. 2012. 中国两栖动物及其分布彩色图鉴[M]. 成都: 四川科学技术出版社.
刘承钊, 胡淑琴, 杨抚华. 1962. 贵州西部两栖类初步调查报告[J]. 动物学报, 14(3): 381-392.
魏刚, 张维勇, 郭鹏. 2017. 梵净山两栖爬行动物[M]. 贵阳: 贵州科技出版社.
伍律, 董谦, 须润华. 1986. 贵州两栖志[M]. 贵阳: 贵州人民出版社.
肖宁, 罗庆华, 罗涛, 等. 2019. 贵州分布长肢林蛙种组物种一新记录——徂徕林蛙[J]. 四川动物, 38(6): 616-622.
杨宝田, 周瑜, 李学龙. 2013. 基于线粒体 16S rRNA 基因徂徕林蛙( R. culaiensis )镇海林蛙( R. zhenhaiensis )的物种分化[J]. 沈阳师范大学学报(自然科学版), 31(2): 143-147. DOI:10.3969/j.issn.1673-5862.2013.02.003
张雷, 冉辉, 沈正雄, 等. 2010. 贵州省两栖动物新记录——镇海林蛙[J]. 安徽农业科学, 38(26): 14442. DOI:10.3969/j.issn.0517-6611.2010.26.100
中国两栖类. 2020. "中国两栖类"信息系统[EB/OL].[2020-05-20]. http://www.amphibiachina.org/.
Chen JM, Zhou WW, Nikolay A, et al. 2017. A novel multilocus phylogenetic estimation reveals unrecognized diversity in Asia toads, genus Megophrys sensu lato (Anura:Megophryidae)[J]. Molecular Phylogenetics and Evolution, 106: 28-43. DOI:10.1016/j.ympev.2016.09.004
David P. 2008. jModelTest:jModeltest phylogenetic model averaging[J]. Molecular Biology and Evolution, 25(7): 1253-1256. DOI:10.1093/molbev/msn083
Frost DR. 2020. Amphibian species of the world version 6.0, an online reference: names assigned to genus Megophrys [DB/OL].[2020-05-20]. http://research.amnh.org/vz/herpetology/amphibia/Amphibia/Anura/Ranidae/ Rana .
Simon C, Frari F, Beckenbaeh A, et al. 1994. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and compilation of conserved polymerase chain reaction primer[J]. Annals of the Entomological Society of America, 87(6): 651-701. DOI:10.1093/aesa/87.6.651
Tamura K, Stecher G, Peterson D, et al. 2013. MEGA6:molecular evolutionary genetics analysis version 6.0[J]. Molecular Biology and Evolution, 30(12): 2725-2729. DOI:10.1093/molbev/mst197
Thompson JD, Gibson TJ, Plewniak F, et al. 1997. The Clustal X windows interface:flexible strategies for multiple sequence alignment aided by quality analysis tools[J]. Nucleic Acids Research, 25(24): 4876-4882. DOI:10.1093/nar/25.24.4876
Wang CC, Qian LF, Zhang CL, et al. 2017. A new species of Rana from the Dabie Mountains in eastern China (Anura, Ranidae)[J]. ZooKeys, 724: 135-153. DOI:10.3897/zookeys.724.19383