Molecular Phylogenetic Status of Rhacophorus laoshan and Zhangixalus yinggelingensis (Anura: Rhacophoridae) from China

Limei Yuan1, Kong Yang1* and Dechun Jiang2*

1Institute of Qinghai-Tibetan Plateau, Southwest Minzu University, Chengdu, Sichuan, P.R. China.

2Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu, Sichuan, P.R. China.

ABSTRACT

Recently, Rhacophorus sensu lato was divided into three genera. In addition, Rhacophorus laoshan and Zhangixalus yinggelingensis were placed in the genus Rhacophorus sensu stricto and Zhangixalus, respectively, only based on morphological characters. The research on the genus Rhacophorus sensu lato faces major challenges due to their complex interspecies relationships, and the phylogenetic status of many species are still unclear, which hampered the taxonomy and protection of these species. In this study, we investigated the molecular phylogenetic status of R. laoshan and Z. yinggelingensis using mitochondrial DNA (12S rRNA, tRNAVal, 16S rRNA and Cyt b) fragments and nuclear DNA (RAG-1, RHOD, TYR) fragments. Our results revealed that Rhacophorus laoshan is closely related to R. verrucopus, R. orlovi and R. calcaneus, and Zhangixalus yinggelingensis should indeed be placed under Zhangixalus. This research clarified the phylogenetic positions of two species within the genus Rhacophorus sensu lato, and it also has an impact on the protection and biogeographic analysis of these species.

Article Information

Received 11 June 2020

Revised 27 September 2020

Accepted 10 November 2020

Available online 30 December 2021

(early access)

Published 15 July 2022

Authors’ Contribution

DCJ and KY conceived the study. LMY performed the analysis and prepared the initial manuscript draft. DCJ approved the final version of the manuscript.

Key words

Rhacophorus, Zhangixalus, Tree frog

DOI: https://dx.doi.org/10.17582/journal.pjz/20200611050629

* Corresponding author: lx-yk@163.com, jiangdc@cib.ac.cn

0030-9923/2022/0005-2417 $ 9.00/0

Copyright 2022 by the authors. Licensee Zoological Society of Pakistan.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

INTRODUCTION

Old world tree frogs of the family Rhacophoridae include 428 species that are widely distributed across Subsaharan Africa and southern Asia (from Sri Lanka, Nepal, and India, to Japan, the Philippines and Sulwesi) (Frost, 2020). The genus Rhacophorus sensu lato is one of the most diverse genera of the family Rhacophoridae, which contains 94 species and is widely distributed across China, Japan, and India, and from the Philippines to Sulawesi (Frost, 2020).

Many new species of the genus Rhacophorus sensu lato have been discovered recently (Dehling, 2015; Hamidy and Kurniati, 2015; Kropachev et al., 2019; Matsui et al., 2013; Nguyen et al., 2017; Pan et al., 2017b; Streicher et al., 2014; Yu et al., 2019), and these studies have made valuable taxonomic contributions of the genus Rhacophorus sensu lato. In addition, in the last ten years, there have been many reports that have inferred the phylogenetic relationships within the genus Rhacophorus sensu lato (Chen et al., 2019; Li et al., 2008, 2009, 2012a,b, 2013; O’Connell et al., 2018a,b; Pan et al., 2017a; Yu et al., 2008). However, the molecular data of many species of Rhacophorus sensu lato are still poorly known and not publicly available, such as these cryptic species R. laoshan and Z. yinggelingensis from southern China. A brown tree frog species, R. laoshan, which was first described in Mo et al. (2008) based on seven adult male specimens collected from Cenwanglaoshan Nature Reserve, Guangxi, China. This frog can be distinguished from all other Asian Rhacophorus Kuhl and van Hasselt, 1822 by the combination of: skin brown and smooth; Y-shaped cartilage visible dorsally on tips of fingers and toes; outer fingers one-third webbed; distinct dermal ridges present on forearms, above vent, and calcars present on heels; anterior and posterior surface of thighs tangerine in color without distinct dark or light spots; tympanum distinct and large, about 6.6% of SVL; dorsum brown with wide dark cross-shaped mark (Mo et al., 2008). To date, R. laoshan is still known from type series only, the validity of this species is based on morphological evidence only (Mo et al., 2008). Similarly, another species, Z. yinggelingensis, was first described in Chou et al. (2007), and it is known from the type locality in the Yinggeling Mountains of Hainan, China (Frost, 2020). It can be distinguished from other species by the following combination of characteristics: snout rounded, without protruding process; green dorsum with few fine white spots, dorsolateral folds absent; skin on dorsum smooth, flanks and venter granulate; metatarsal skinfolds faint; fingers and toes webbed; front of thigh yellow and red-tinged, rear of thigh and inner side of tibia red; white-tipped tuberculous dermal ridge on cloaca present, not well developped; iris silvery, partly light golden in upper half (Chou et al., 2007). In addition, no further specimens or new distribution record has been reported for these two poorly known tree frogs, to the best of our knowledge (Fei et al., 2009, 2012; Frost, 2020).

Recently, Jiang et al. (2019) divided Rhacophorus sensu lato into three genera: genus Rhacophorus sensu stricto Kuhl and Van Hassalt, 1822, genus Leptomantis Peters, 1867, and a new genus erected, namely Zhangixalus Jiang et al., 2019. Furthermore, R. laoshan and Z. yinggelingensis were placed in the genus Rhacophorus sensu stricto and Zhangixalus, respectively, only based on morphological characters (Jiang et al., 2019). Due to the lack of molecular data for the two species, their phylogenetic status remains unclear. Previous studies suggested that it is essential to clarify species’ phylogenetic status using molecular data. We thus conducted molecular phylogenetic analyses of the two Rhacophorus sensu lato species to clarify their taxonomic status.

MATERIALS AND METHODS

Ethical statement

This study was performed in accordance with the approval of Experimental Animal Ethics Committee of Chengdu Institute of Biology, Chinese Academy of Sciences.

Taxon sampling and data collection

Two specimens were collected in this study, including one specimen of Z. yinggelingensis (HN2018002), which was collected in Hainan, China on 2018 (19°3’3.84N, 109°32’20.60″E) at 1,569 m altitude (Fig. 1), and one specimen of R. laoshan (1705014), which was collected from Cenwanglaoshan Nature Reserve (Guangxi, China) (Fig. 1). After euthanization, liver tissues were preserved in 95% ethanol for DNA extraction. Except for the two individuals, sequences of 39 individuals were collected from GenBank in this study, including sequences for eight species of Rhacophorus sensu stricto, one species of Leptomantis, and 18 species of Zhangixalus as ingroups and two species of Polypedates as outgroups according to previous studies (Jiang et al., 2019; Li et al., 2008, 2009). All novel sequences obtained in this study were deposited in GenBank (Accession No. are shown in Supplementary Table S1).

DNA extraction, amplification, and sequencing

Genomic DNA of R. laoshan was extracted from ethanol-preserved liver tissue samples using an Ezup Column Animal Genomic DNA Purification Kit (Sangon Biotech, China), according to the protocols of the manufacturer. A total of 1890 bp of the 12S rRNA, tRNAVal, and 16S rRNA was targeted and amplified using the primers that are listed in Table I. Polymerase chain reaction (PCR) were performed at 25 μL volume, including 12.5 µl Taq PCR Master Mix (2 X blue dye), 9.5 µl ddH2O, 1µl FS01, 1µl Rend (C=10µl/L), and 1 µl DNA template. PCR cycles consisted of initial denaturation at 95°C for 3 min followed by 35 cycles: denaturation at 94°C for 1 min, annealing at 55°C for 1 min, and elongation at 72°C for 1 min. A final elongation step of 10 min at 72°C was performed.

 

Table I. Primers used in this study.

Primer	Primer sequence	Resource
FS01	5’-AAC GCT AAG ATG AAC CCT AAA AAG TTC T-3’	Wilkinson et al., 2002
Rend	5'-GAC CTG GAT TAC TCC GGT CTG A-3'	Wilkinson et al., 2002
F0483	5'-GAA GAG GCA AGT CGT AAC ATG G -3'	Wilkinson et al., 2002
F1209	5’-CTC GTA CCT TTT GCA TCA TG-3’	Yang, 2018
R1209	5’-CAT GAT GCA AAA GGT ACG AG-3’	Yang, 2018

 

In order to conduct further research, we used high-throughput sequencing technology to extract genes of Z. yinggelingensis. A total of 0.88 μg DNA extracted from Z. yinggelingensis muscle tissue was sheared into short fragments (<500bp) using the Covaris system to construct pair-end DNA libraries. DNA libraries were constructed and sequenced using Illumina-Hi-Seq 2500. The sample was sequenced to a target 10×depth (30Gbp) of its genome. Reads containing more low-quality bases (phred score <30) than 20% were filtered out. Gene fragments of other species in this study were aligned using MEGA 7.0 (Kumar et al., 2016) and adjusted manually. Then, we used the MEANGS (https://github.com/YanCCscu/MEANGS) to get the purpose genes. The process are as follows: The alignments of genes that aligned using MEGA 7.0 were used in construction profile HMM(s) with hmmer-build in HMMER V3.16 (Eddy, 2008). Hereafter, nhmmer were used to search reads that fitting the hmm profiler for each gene from high-throughout sequencing data. Finally, we assembled these matched reads with the de novo assembler SSAKE v2.2.1 (Warren et al., 2007), and the assembled scaffolds were identified by blastting to the purpose DNA fragments, including four mitochondrial fragments, 12S rRNA, tRNAVal, 16S rRNA, Cytochrome b (Cyt b), and three nuclear genes, Recombination activating gene 1 (RAG-1), Exon 1 of tyrosinase (RHOD) and exon 1 of tyrosinase (TYR). In order to evaluate the accuracy of the assembled genes of Z. yinggelingensis, we mapped all reads back to the assembled genes, and statistic the mapping quality and average depth of these genes.

Phylogenetic analyses

Phylogenetic trees were reconstructed by only mitochondrial fragments (12S rRNA, tRNAVal, 16S rRNA) and combined mitochondrial and nuclear fragments (12S rRNA, tRNAVal, 16S rRNA, Cyt b, RAG-1, RHOD, TYR), respectively. Sequences were aligned using MEGA 7.0 (Kumar et al., 2016) with default settings and rechecked in MAFFT version 7 (Katoh and Standley, 2013), then adjusted manually. Different gene fragments of each individual were concatenated by using Phylosuite (Zhang et al., 2020), and positions of indels were treated as missing data for all datasets. The most optimum substitution model of evolution (GTR+I+G for mitochondrial fragments, HKY+I for RAG-1 gene and RHOD gene, and SYM+G for TYR gene) was calculated by Partition Finder v2.1.1 (Lanfear et al., 2017) under Akaike Information Criterion (AIC) (lnL= -18812.57, AICc= 37858.23). We then reconstructed phylogenetic trees using the Bayesian inference (BI) and maximum likelihood (ML) methods, respectively.

The Bayesian phylogenetic tree was constructed using MrBayes 3.2.0 (Ronquist et al., 2012). Runs of Markov chains for 5,000,000 generations were summarized and sampled every 100 generations. The first 25% of samples as burn-in, and convergence was investigated using the parameter of average standard deviation of split frequencies ≤ 0.01 (Ronquist et al., 2012). Nodes were considered strongly supported when BPP ≥ 0.95. Software RAxML v8.2.10 (Stamatakis, 2014) was used to implement Maximum Likelihood (ML) analysis under the best-fit model of evolution (GTRGAMMA) based on the AIC criterion. And bootstrap proportions (BSP) were evaluated with 1000 nonparametric bootstrap pseudoreplications, and it were considered significantly supported when the node’s BSP ≥ 70. In addition, uncorrected pairwise distance (p-distances) of the four mitochondrial fragments (12S rRNA, tRNAVal, 16S rRNA and Cyt b) and three nuclear genes (RAG-1, RHOD and TYR) dataset among species in Rhacophorus sensu lato were calculated using MEGA 7 (Kumar et al., 2016), with the “P-distance” model and “Transitions + Transversions” substitutions to include.

RESULTS AND DISCUSSION

The mapping average depth of 12S rRNA, tRNAVal, and 16S rRNA, Cyt-b, RAG-1, RHOD, and TYR for the Z. yinggelingensis is 3331, 2315, 13.95, 5.847, 9.681, respectively. Simultaneously, the mitochondrial fragments consisted of 1890 bp, and the combined mitochondrial and nuclear fragments consisted of 3746 bp, including 1890 bp for 12S rRNA, tRNAVal, and 16S rRNA, 546 bp for Cyt-b, 507 bp for RAG-1, 279 bp for RHOD, 524 bp for TYR. The present molecular phylogeny is concordant with previous molecular studies regarding the monophyly of Rhacophorus sensu lato (Chan et al., 2018, 2019; Li et al., 2008, 2009, 2012a, 2013; O’Connell et al., 2018b; Pan et al., 2017a). Both mitochondrial fragments (Supplementary Fig. S1) and combined mitochondrial and nuclear fragments (Fig. 2 and Supplementary Fig. S2) resulted in virtually identical topology. Although interspecific relationships were not fully resolved, the phylogenetic trees from both the BI and ML methods showed that all species contained in Rhacophorus sensu lato were clustered to three clades as monophyletic groups and each clade consisted of the same species. According to phylogenetic trees from combined mitochondrial and nuclear fragments, clade A was congruent with genus Zhangixalus and was strongly supported by a bootstrap proportion of 83 and a posterior probability of 0.96. Clade B was congruent with genus Leptomantis and it was sister group to Zhangixalus that supported by a bootstrap proportion of 80 and a posterior probability of 0.96, which was consistent with Jiang et al. (2019), but was different from Chen et al. (2019). Chen et al. (2019) indicated that Leptomantis was sister group to Rhacophorus sensu stricto. They used anchored hybrid enrichment (AHE) targeted sequencing and their sequences were long larger than this study. Clade C was congruent with genus Rhacophorus sensu stricto and supported by a bootstrap proportion of 100 and a posterior probability of 1.00.

The phylogeny (Fig. 2 and Supplementary Fig. S2) depicted that R. laoshan belongs to Rhacophorus sensu stricto and Z. yinggelingensis was nested within Zhangixalus. Jiang et al. (2019) propose the morphological diagnosis and difference of each genus of Rhacophorus sensu lato. They shed light on Zhangixalus differs from Rhacophorus sensu stricto by the absence of dermal folds along limbs and tarsal projections (vs. present), the absence of supracloacal fold (vs. present or not). And we found that R. laoshan and Z. yinggeligensis corresponds with the morphological diagnosis of Rhacophorus sensu stricto and Zhangixalus, respectively. On the one hand, the morphological characteristics of R. laoshan similar to the characteristics diagnosis of Rhacophorus sensu lato in the presence of dermal folds along limbs (Mo et al., 2008; Jiang et al., 2019). On the other hand, the morphological characteristics of Z. yinggelingensis similar to the following combination of characteristics diagnosis of Zhangixalus: Dermal folds along limbs absent; tarsal projections absent; body size relatively moderate or large, SVL 30–120 mm (mostly above 50 mm); dorsal coloration mostly green (Chou et al., 2007; Jiang et al., 2019). Consequently, based on our molecular results and published morphological data, we agree with Jiang et al. (2019), place R. laoshan and Z. yinggelingensis into the genera Rhacophorus and Zhangixalus, respectively.

Furthermore, in the clade Rhacophorus sensu stricto, R. laoshan was the sister-taxon to the assembly of three species (BPP = 1.00; BSP = 99), including R. verrucopus (6254 RAO), R. orlovi (AMNH A161405) and R. calcaneus (FMNH 256465; Fig. 2). However, we cannot confirm the sister species of R. laoshan, and it may be related to R. calcaneus, which was recovered as polyphyletic and formed two distinct clades. Specifically, between the two R. calcaneus clades, one taxon (FMNH 256465, from Laos) was sister to R. orlovi (AMNH A161405; BPP = 1.00; BSP = 100) possessed a genetic divergence of 0.2% (Supplementary Table S2). Additionally, the second R. calcaneus (AMNH A163749, from Vietnam) and R. robertingeri (VNMN:4123, from Vietnam) formed a well-supported clade with short branch length (BPP = 1.00; BSP = 100), and this clade possessed a genetic divergence of 0.3% (Supplementary Table S2). However, on the one hand, the localities of the species R. calcaneus were certainty only from three places (Frost, 2020), and the localities of the two specimens of R. calcaneus (FMNH 256465, from Laos; AMNH A163749, from Quang Nam, Vietnam) didn’t fall within the known distribution sites of R. calcaneus. On the other hand, some species are poorly diagnosed yet and characters highly convergent in previous studies (David et al., 2007, 2015a; Guo et al., 2014). Consequently, we suspect that the two specimens (FMNH 256465, AMNH A163749) of R. calcaneus were misidentifications. Due to we do not have the two specimens, it is need to confirm whether the two specimens were misidentified in future study. Moreover, it may be that some sequences of the species R. calcaneus from the NCBI have problems, because taxonomic identification problems of existing Genbank sequences are very common (Guo et al., 2014; Ren et al., 2018).

 

For Z. yinggelingensis, the p-distances varied from 3.0% to 7.0% between it and other species of Zhangixalus that collected in this study [minimum distance to Z. duboisi; maximum distance to Z. achantharhena (Supplementary Table S2)]. Additionally, Chou et al. (2007) suggested that Z. yinggelingensis seems to be closely related to Z. arvalis, Z. aurantiventris, Z. chenfui, Z. dorsoviridis, Z. hungfuensis, Z. taipeianus, and Z. yaoshanensis, Z. moltrechti, Z. nigropunctatus on the basis of morphological evidences. However, there are many morphologically similar and conservative species in Zhangixalus (Li et al., 2012), and we can’t distinguish them well. Thus, we suggest that it not only use morphological and molecular datasets, but also their distributions when we identified the species, which can reduce the probability of misidentifications.

Although phylogenetic status of Z. yinggelingensis is not resolved well, we recognize the validity of it using the molecular data for the first time. However, why phylogenetic status of Z. yinggelingensis is not resolved well? The possible reasons are as follows: (1) The sequence of Z. yinggelingensis in this study is short; (2) There are only some species of Zhangixalus in this study; (3) the genus Zhangixalus evolution history is complex. As genomic data have become more common in phylogenetic studies, many studies have been shown that it can use genomic data resolving both shallow and deep-scale evolutionary relationships (Brandrud et al., 2019; Burbrink and Gehara, 2018; Chen et al., 2019; Godefroid et al., 2019; Kuntner et al., 2019; Spriggs et al., 2019). Thus, it seems possible that we can use genomic data to resolve the phylogenetic status of Z. yinggelinensis in future study. In addition, the research on amphibians in Southern China was scarce, and there have been some cryptic species were discovered in the past three years (Liu et al., 2017; Pan et al., 2017b; Yu et al., 2019). Finally, we suggest that we can strengthen the investigation of amphibian diversity in southern China, which lay a foundation for the study of biogeography in southern China.

CONCLUSIONS

In this study, we have clarified the phylogenetic status of R. laoshan and Z. yinggelingensis using mitochondrial fragments (12S rRNA, tRNAVal, 16S rRNA and Cyt b) and nuclear DNA (RAG-1, RHOD, TYR) fragments. Our results suggest that R. laoshan and Z. yinggelingensis belongs to Rhacophorus sensu stricto and Zhangixalus, respectively. Finally, this research play an important role in the further studies, such as protection, biogeographical patterns and the evolution of life history strategies of these species.

ACKNOWLEDGMENTS

This work was supported by the National Natural Science Foundation of China (31900322; 32070410); Sichuan Science and Technology Bureau (2020YFH0005); and Innovative Research Project for Postgraduates of Southwest Minzu university (CX2019SZ89).

Supplementary material

There is supplementary material associated with this article. Access the material online at: https://dx.doi.org/10.17582/journal.pjz/20200611050629

Statement of conflict of interest

The authors have declared no conflict of interest.

REFERENCES

Brandrud, M.K., Paun, O., Lorenz, R., Baar, J. and Hedrén, M., 2019. Restriction-site associated DNA sequencing supports a sister group relationship of Nigritella and Gymnadenia (Orchidaceae). Mol. Phylogenet. Evol., 136: 21-28. https://doi.org/10.1016/j.ympev.2019.03.018

Burbrink, F.T. and Gehara, M., 2018. The biogeography of deep time phylogenetic reticulation. Syst. Biol., 67: 743-755. >https://doi.org/10.1093/sysbio/syy019

Chan, K.O., Grismer, L.L. and Brown, R.M., 2018. Comprehensive multi-locus phylogeny of old world tree frogs (Anura: Rhacophoridae) reveals taxonomic uncertainties and potential cases of over and underestimation of species diversity. Mol. Phylogenet. Evol., 127: 1010-1019. >https://doi.org/10.1016/j.ympev.2018.07.005

Chen, J.M., Prendini, E., Wu, Y.H., Zhang, B.L., Suwannapoom, C., Chen, H.M., Jin, J.Q., Moriarty Lemmon, E., Lemmon, A.R., Stuart, B.L., Raxworthy, C.J., Murphy, R.W., Yuan, Z.Y. and Che, J., 2019. An integrative phylogenomic approach illuminates the evolutionary history of old world tree frogs (Anura: Rhacophoridae). Mol. Phylogenet. Evol., 145: 106724. >https://doi.org/10.1016/j.ympev.2019.106724

Chou, W., Lau, M. and Chan, B.P., 2007. A new treefrog of the genus Rhacophorus (Anura: Rhacophoridae) from Hainan Island, China. Raffles Bull. Zool., 55: 157-165.

David, P., Agarwal, I., Athreya, R., Mathew, R., Vogel, G. and Mistry, V.K., 2015a. Revalidation of Natrix clerki Wall, 1925, an overlooked species in the genus Amphiesma Duméril, Bibron and Duméril, 1854 (Squamata: Natricidae). Zootaxa, 3919: 375–395. >https://doi.org/10.11646/zootaxa.3919.2.9

David, P., Bain, R.H., Nguyen, T.Q., Orlov, N.L., Vogel, G., Vu, T.N. and Ziegler, T., 2007. A new species of the natricine snake genus Amphiesma from the Indochinese Region (Squamata: Colubridae: Natricinae). Zootaxa, 1462: 41–60. >https://doi.org/10.11646/zootaxa.1462.1.2

Dehling, J.M., 2015. A new species of Rhacophorus (Anura: Rhacophoridae) from Gunung Kinabalu, Borneo. Salamandra, 51: 1-11.

Eddy S.R., 2008. A probabilistic model of local sequence alignment that simplifies statistical significance estimation. PLoS Comput. Biol., 4: e1000069. >https://doi.org/10.1371/journal.pcbi.1000069

Fei, L., Hu, S.Q., Ye, C.Y. and Huang. Y.Z., 2009. Fauna sinica. Amphibia. Volume 2, Chinese academy of Science. Science Press, Anura, Beijing.

Fei, L., Ye, C.Y. and Jiang, J.P., 2012. Colored atlas of Chinese amphibians and their distributions. Sichuan Publishing House of Science and Technology, Sichuan, China.

Frost, D.R., 2020. Amphibian species of the world: An online reference. Version 6.0. American Museum of Natural History, New York, USA. Partners Available at: https://amphibiansoftheworld.amnh.org/ (accessed 10 May 2020).

Godefroid, M., Meseguer, A.S., Sauné, L., Genson, G., Streito, J.C., Rossi, J.P., Riverón, A.Z., Mayer, F., Cruaud, A. and Rasplus, J.Y., 2019. Restriction-site associated DNA markers provide new insights into the evolutionary history of the bark beetle genus Dendroctonus. Mol. Phylogenet. Evol., 139: 106528. >https://doi.org/10.1016/j.ympev.2019.106528

Guo, P., Zhu, F., Liu, Q., Zhang, L., Li, J.X., Huang, Y.Y. and Pyron, R.A., 2014. A taxonomic revision of the Asian keelback snakes, genus Amphiesma (Serpentes: Colubridae: Natricinae), with description of a new species. Zootaxa, 3873: 425–440. >https://doi.org/10.11646/zootaxa.3873.4.5

Hamidy, A. and Kurniati, H., 2015. A new species of tree frog genus Rhacophorus from Sumatra, Indonesia Amphibia, Anura). Zootaxa, 3947: 49-66. >https://doi.org/10.11646/zootaxa.3947.1.3

Jiang, D.C., Jiang, K., Ren J. L., Wu, J. and Li, J.T, 2019. Resurrection of the genus Leptomantis, with description of a new genus to the family Rhacophoridae (Amphibia: Anura). Asian Herpetol. Res., 10: 1–12.

Katoh, K. and Standley, D.M., 2013. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol., 30: 772-780. >https://doi.org/10.1093/molbev/mst010

Kropachev, I.I., Orlov, N.L., Ninh, H.T. and Nguyen, T.T., 2019. A new species of Rhacophorus genus (Amphibia: Anura: Rhacophoridae: Rhacophorinae) from Van Ban District, Lao Cai Province, Northern Vietnam. Rus. J. Herpetol., 26: 325-334. >https://doi.org/10.30906/1026-2296-2019-26-6-325-334

Kuhl, H. and Van Hasselt, J.C., 1822. Uittreksels uit breieven van de Heeren Kuhl en van Hasselt, aan de Heeren C. J. Temminck, Th. van Swinderen en W. de Haan. Algemeene Konst-en Letter-Bode, 7: 99–104.

Kumar, S., Stecher, G. and Tamura, K., 2016. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol., 33: 1870-1874. >https://doi.org/10.1093/molbev/msw054

Kuntner, M., Hamilton, C.A., Cheng, R.C., Gregorič, M., Lupše, N., Lokovšek, T., Lemmon, E.M., Lemmon, A.R., Agnarsson, I. and Coddington, J.A., 2019. Golden orbweavers ignore biological rules: phylogenomic and comparative analyses unravel a complex evolution of sexual size dimorphism. Syst. Biol., 68: 555-572. >https://doi.org/10.1093/sysbio/syy082

Lanfear, R., Frandsen, P.B., Wright, A.M., Senfeld, T. and Calcott, B., 2017. Partition finder 2: New methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. Evol., 34: 772-773. >https://doi.org/10.1093/molbev/msw260

Li, J.T., Li, Y., Murphy, R.W., Rao, D.Q. and Zhang, Y.P., 2012a. Phylogenetic resolution and systematics of the Asian tree frogs, Rhacophorus (Rhacophoridae, Amphibia). Zool. Scr., 41: 557-570. >https://doi.org/10.1111/j.1463-6409.2012.00557.x

Li, J.T., Liu, J., Chen, Y.Y., Wu, J.W., Murphy, R.W., Zhao, E.M., Wang, Y.Z. and Zhang, Y.P., 2012b. Molecular phylogeny of treefrogs in the Rhacophorus dugritei species complex (Anura: Rhacophoridae), with descriptions of two new species. Zool. J. Linn. Soc., 165: 143-162.

Li, J.T., Che, J., Bain, R.H., Zhao, E.M. and Zhang, Y.P., 2008. Molecular phylogeny of Rhacophoridae (Anura): A framework of taxonomic reassignment of species within the genera Aquixalus, Chiromantis, Rhacophorus and Philautus. Mol. Phylogenet. Evol., 48: 302-312. >https://doi.org/10.1016/j.ympev.2008.03.023

Li, J.T., Che, J., Murphy, R.W., Zhao, H., Zhao, E.M., Rao, D.Q. and Zhang, Y.P., 2009. New insights to the molecular phylogenetics and generic assessment in the Rhacophoridae (Amphibia: Anura) based on five nuclear and three mitochondrial genes, with comments on the evolution of reproduction. Mol. Phylogenet. Evol., 53: 509-522. >https://doi.org/10.1016/j.ympev.2009.06.023

Li, J.T., Li, Y., Klaus, S., Rao, D.Q., Hillis, D.M. and Zhang, Y.P., 2013. Diversification of Rhacophorid frogs provides evidence for accelerated faunal exchange between India and Eurasia during the Oligocene. Proc. Natl. Acad. Sci. U.S.A., 110: 3441-3446. >https://doi.org/10.1073/pnas.1300881110

Liu, B.Q., Wang, Y.F., Jiang, K., Chen, H.M., Zhou, J.J., Xu J.N. and Wu C.H., 2017. A new treefrog species of the genus Rhacophorus found in Zhejiang, China (Anura: Rhacophoridae). Chinese J. Zool., 52: 361-372.

Matsui, M., Shimada, T. and Sudin, A., 2013. A new gliding frog of the genus Rhacophorus from Borneo. Curr. Herpetol., 32: 112-124. >https://doi.org/10.5358/hsj.32.112

Mo, Y., Jiang, J.P., Xie, F. and Ohler, A., 2008. A new species of Rhacophorus (Anura: Ranidae) from China. Asiatic Herpetol. Res., 11: 85-92.

Nguyen, T.T., Pham, C.T., Nguyen, T.Q., Ninh, H.T. and Ziegler, T., 2017. A new species of Rhacophorus (Amphibia: Anura: Rhacophoridae) from Vietnam. Asian Herpetol. Res., 8: 221-234.

O’Connell, K.A., Hamidy, A., Kurniawan, N., Smith, E.N. and Fujita, M.K., 2018a. Synchronous diversification of parachuting frogs (Genus Rhacophorus) on Sumatra and Java. Mol. Phylogenet. Evol., 123: 101-112. >https://doi.org/10.1016/j.ympev.2018.02.003

O’Connell, K.A., Smart, U., Smith, E.N., Hamidy, A., Kurniawan, N. and Fujita, M.K., 2018b. Within-island diversification underlies parachuting frog (Rhacophorus) species accumulation on the Sunda Shelf. J. Biogeogr., 45: 929-940. >https://doi.org/10.1111/jbi.13162

Pan, T., Zhang, Y., Wang, H., Wu, J., Kang, X., Qian, L., Chen, J., Rao, D., Jiang, J. and Zhang, B., 2017a. The reanalysis of biogeography of the Asian tree frog, Rhacophorus (Anura: Rhacophoridae): geographic shifts and climatic change influenced the dispersal process and diversification. Peer J., 5: e3995. >https://doi.org/10.7717/peerj.3995

Pan, T., Zhang, Y., Wang, H., Wu, J., Kang, X., Qian, L., Li, K., Zhang, Y., Chen, J., Rao, D., Jiang, J. and Zhang, B., 2017b. A new species of the genus Rhacophorus (Anura: Rhacophoridae) from Dabie Mountains in East China. Asian Herpetol. Res., 8: 1-13.

Peters, W.C.H., 1867. Herpetologische Notizen. Monatsberichte der Königlichen Preussische Akademie des Wissenschaften zu Berlin, 1867: 13–37.

Ren, J.L., Wang, K., Nguyen, T.T., Van Hoang, C., Zhong, G.H, Jiang, K., Guo, P. and Li, J.T., 2018. Taxonomic re-evaluation of the monotypic genus Pararhabdophis Bourret, 1934 (Squamata: Colubridae: Natricinae) with discovery of its type species, P. chapaensis, from China. Zootaxa, 4486: 31-56. >https://doi.org/10.11646/zootaxa.4486.1.2

Ronquist, F., Teslenko, M., Van, d.M.P., Ayres, D.L., Darling, A., Höhna, S., Larget, B., Liu, L., Suchard, M.A. and Huelsenbeck, J.P., 2012. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol., 61: 539-542. >https://doi.org/10.1093/sysbio/sys029

Spriggs, E.L., Eaton, D.A., Sweeney, P.W., Schlutius, C., Edwards, E.J. and Donoghue, M.J., 2019. Restriction-site-associated DNA sequencing reveals a cryptic Viburnum species on the North American coastal plain. Syst. Biol., 68: 187-203. >https://doi.org/10.1093/sysbio/syy084

Stamatakis, A., 2014. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics, 30: 1312-1313. >https://doi.org/10.1093/bioinformatics/btu033

Streicher, J.W., Hamidy, A., Harvey, M.B., Anders, B., Shaney, K.J., Kurniawan, N. and Smith, E.N., 2014. Mitochondrial DNA reveals a new species of parachuting frog (Rhacophoridae: Rhacophorus) from Sumatra. Zootaxa, 3878: 351-356. >https://doi.org/10.11646/zootaxa.3878.4.2

Warren, R.L., Sutton, G.G., Jones, S.J. and Holt, R.A., 2007. Assembling millions of short DNA sequences using SSAKE. Bioinformatics, 23: 500-501. >https://doi.org/10.1093/bioinformatics/btl629

Wilkinson J.A., Drewes R.C. and Tatum O.L., 2002. A molecular phylogenetic analysis of the family Rhacophoridae with an emphasis on the Asian and African genera. Mol. Phylogenet. Evol., 24: 265-273. >https://doi.org/10.1016/S1055-7903(02)00212-9

Yang T., 2018. Phylogenetic, biogeography and evolution of Rhacophorus (Rhacophoridae: Amphibia). Master’s thesis of the Chinese Academy of Sciences, Beijing, China.

Yu, G.H., Hui, H., Hou, M., Wu, Z., Rao, D. and Yang, J., 2019. A new species of Zhangixalus (Anura: Rhacophoridae), previously confused with Zhangixalus smaragdinus (Blyth, 1852). Zootaxa, 4711: 275-292. >https://doi.org/10.11646/zootaxa.4711.2.3

Yu, G.H., Rao, D.Q., Yang, J.X. and Zhang, M.W., 2008. Phylogenetic relationships among rhacophorinae (Rhacophoridae, Anura, Amphibia), with an emphasis on the Chinese species. Zool. J. Linn. Soc., 153: 733-749. >https://doi.org/10.1111/j.1096-3642.2008.00404.x

Zhang, D., Gao, F., Jakovlić, I., Zou, H., Zhang, J., Li, W.X. and Wang, G.T., 2020. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour., 20: 348-355. >https://doi.org/10.1111/1755-0998.13096