扩展功能
文章信息
- 程彦林, 刘京, 陈龙, 雷诗融, 邹家伟, 张霞, 李仕泽
- 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
蛙属 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 ,为贵州省两栖动物新记录。本文对其形态特征及生态习性进行描述,并初步探讨其形态特征的地理变异,为蛙属物种的系统学和谱系地理学研究提供了新信息。
物种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 |
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)。比较此次所采集的林蛙标本和其他蛙属物种的形态鉴定特征,发现本次采集的标本形态特征与大别山林蛙模式产地标本相符。
形态量度 | 雄性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 |
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)的描述吻合,系统发育分析也支持采自桐梓的林蛙标本与大别山林蛙关系较近。
物种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 等。
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 |