Specimens were collected in Suan Phueng District, Ratchaburi Province, Thailand by Parinya Pawangkhanant, Platon V. Yushchenko, Mali Naiduangchan, Chatmongkon Suwannapoom, and Nikolay A. Poyarkov during several field surveys in 2019 (Fig. 2). Geographic coordinates and elevations were obtained using a Garmin GPSMAP 60CSx and recorded in WGS 84 datum. Specimens were collected by hand, anaesthetized and euthanized with 20% benzocaine solution, femoral muscle or liver tissue samples removed, and then fixed in 10% buffered formalin before preserving in 70% ethanol. The tissue samples were stored in 95% ethanol. Specimens and tissues were subsequently deposited in the herpetological collections of the School of Agriculture and Natural Resources, University of Phayao (AUP, Phayao, Thailand) and of the Zoological Museum of Moscow University (ZMMU, Moscow, Russia).

Figure 2. 

A. Bayesian phylogeny of the Cyrtodactylus brevipalmatus group based on 1397 base pairs of ND2 with BBP and UFB nodal support values, respectively at the nodes. B. PCA based on concatenated meristic and adjusted morphometric data. C. DAPC based on concatenated meristic and adjusted morphometric data.

Specimen collection and animal use protocols were approved by the Institutional Ethical Committee of Animal Experimentation of the University of Phayao, Phayao, Thailand (certificate number UP-AE61-01-04-022, issued to Chatmongkon Suwannapoom) and were strictly complacent with the ethical conditions of the Thailand Animal Welfare Act. Field work, including collection of animals in the field and specimen exportation, was authorized by the Institute of Animals for Scientific Purpose Development (IAD), Bangkok, Thailand (permit numbers U1-01205-2558 and UP-AE59-01-04-0022, issued to Chatmongkon Suwannapoom).

Morphological data included both meristic and morphometric characters. Measurements were taken on the left side of the body when possible to the nearest 0.1 mm using dial calipers under dissecting microscope following Murdoch et al. (2019). Measurements taken were: snout-vent length (SVL), taken from the tip of the snout to the vent; tail length (TL), taken from the vent to the tip of the tail, original or partially regenerated; tail width (TW), taken at the base of the tail immediately posterior to the postcloacal swelling; humeral length (HumL), taken from the medial end of the humerus at its insertion point in the glenoid fossa to the distal margin of the elbow while flexed 90°; forearm length (ForL), taken on the ventral surface from the posterior margin of the elbow while flexed 90° to the inflection of the flexed wrist; femur length (FemL), taken from the medial end of the femur at its insertion point in the acetabulum to the distal margin of the knee while flexed 90°; tibia length (TibL), taken on the ventral surface from the posterior surface of the knee while flexed 90° to the base of the heel; axilla to groin length (AG), taken from the posterior margin of the forelimb at its insertion point on the body to the anterior margin of the hind limb at its insertion point on the body; head length (HL), the distance from the posterior margin of the retroarticular process of the lower jaw to the tip of the snout; head width (HW), measured at the angle of the jaws; head depth (HD), the maximum height of head measured from the occiput to base of the lower jaw; eye diameter (ED), the greatest horizontal diameter of the eye-ball; eye to ear distance (EE), measured from the anterior edge of the ear opening to the posterior edge of the bony orbit; eye to snout distance or snout length (ES), measured from anteriormost margin of the bony orbit to the tip of snout; eye to nostril distance (EN), measured from the anterior margin of the bony orbit to the posterior margin of the external nares; interorbital distance (IO), measured between the dorsal-most edges of the bony orbits; internarial distance (IN), measured between the nares across the rostrum; and ear length (EL), greatest oblique length across the auditory meatus.

Meristic characters evaluated were the number of supralabial scales (SL), counted from the largest scale immediately below the middle of the eyeball to the rostral scale; infralabial scales (IL), counted from the mental to the termination of enlarged scales just after the upturn of the mouth; the number of paravertebral tubercles (PVT), between limb insertions counted in a straight line immediately left of the vertebral column; the number of longitudinal rows of body tubercles (LRT), counted transversely across the center of the dorsum from one ventrolateral fold to the other; the number of longitudinal rows of ventral scales (VS,) counted transversely across the center of the abdomen from one ventrolateral fold to the other; the number of transverse rows of ventral scales (VSM), counted along the midline of the body from the postmentals to the cloacal opening; the number of expanded subdigital lamellae on the fourth toe proximal to the digital inflection (TL4E), counted from the base of the first phalanx where it contacts the body of the foot to the largest scale on the digital inflection—the large contiguous scales on the palmar and plantar surfaces were not counted; the number of small, generally unmodified subdigital lamellae distal to the digital inflection on the fourth toe (TL4U), counted from the digital inflection to the claw including the claw sheath (see Murdoch et al. 2019: Fig. 2); the total number of subdigital lamellae (TL4T) beneath the fourth toe (i.e. TL4E + TL4U = TL4T); the number of expanded subdigital lamellae on the fourth finger proximal to the digital inflection (FL4E) counted from the base of the first phalanx where it contacts the body of the foot to the largest scale on the digital inflection—the large contiguous scales on the palmar and plantar surfaces were not counted; the number of small, generally unmodified subdigital lamellae distal to the digital inflection on the fourth finger (FL4U) counted from the digital inflection to the claw including the claw sheath (see Murdoch et al. 2019: Fig. 2); the total number of subdigital lamellae (FL4T) beneath the fourth toe (i.e. FL4E + FL4U = FL4T); the total number of enlarged femoral scales (FS) from each thigh combined as a single metric; the number of enlarged precloacal scales (PCS); the number of precloacal pores in (PP) in males; the number of femoral pores in (FP) in males; the number of rows of post-precloacal scales (PPS) on the midline between the enlarged precloacal scales and the vent; and the number of dark body bands (BB) between the nuchal loop (the dark band running from eye to eye across the nape) and the hind limb insertions. Other morphological characters evaluated were the presence or absence of paravertebral tubercles, enlarged femoral scales, small tubercles on the forelimbs and flanks, paried or single enlarged subcaudal scales; large or small dorsolateral and ventrolateral caudal fringe scales (i.e. fringes), and the cross-section of the tail round or more square in shape.

Small sample sizes (n=1–3) from some of the populations/species precluded meaningful statistical analyses. However, the morphospatial clustering among the species for the meristic, morphometric, and concatenated datasets was visualized using principal component analysis (PCA) and discriminant analysis of principal components (DAPC) from the ADEGENET package in R (Jombart et al. 2010). We analyzed each data type (i.e. meristic and morphometric) separately in order to visualize the contribution of each to the variation in the group and then we concatenated the data types to visualize their overall contribution to the variation of the group. Previous studies have demonstrated that concatenated data sets generally outperform single data sets in capturing the overall morphological variation among closely related species (Grismer et al. 2018, 2020a, 2021).

Morphometric characters analyzed were SVL, AG, HumL, ForL, FemL, TibL, HL, HW, HD, ED, EE, EN, IO, EL, and IN. To remove potential effects of allometry in the morphometric dataset, size was normalized using the following equation: X_adj=log(X)-β[log(SVL)-log(SVL_mean)], where X_adj=adjusted value; X=measured value; β=unstandardized regression coefficient for each population; and SVL_mean=overall average SVL of all populations (Thorpe 1975, 1983; Turan 1999; Lleonart et al. 2000, accessible in the R package GroupStruct (available at https://github.com/chankinonn/GroupStruct). The morphometrics of each species were normalized separately and then concatenated so as not to conflate intra- with interspecific variation (Reist 1986). All data were scaled to their standard deviation to insure they were analyzed on the basis of correlation and not covariance.

Meristic characters analyzed were SL, IL, LRT, VS, VSM, TL4E, TL4U. TL4T, FL4E, FL4U, FL4T, PCS, PPS, and BB. Paravertebral tubercles (PVT), enlarged femoral scales (FS) were not included due to their absence in Cyrtodactylus elok. Precloacal and femoral pores were also omitted due to their absence in females. For corroboration of the PCA, a discriminant analysis of principal components (DAPC) was performed on the same data sets. DAPC is a supervised analysis that relies on scaled data calculated from its own internal PCA as a prior step to ensure that variables analyzed are not correlated and number fewer than the sample size. Dimension reduction of the DAPC prior to plotting, is accomplished by retaining the first set of PCs that account for 90–99% of the variation in the data set (Jombart and Collins 2015) as determined from a scree plot generated as part of the analysis. The ranges, means, medians, and 50% quartiles were visualized using boxplots for the meristic data and violin plots with embedded boxplots for the mensual data. All analyses were performed in R [v3.4.3].

For the molecular phylogenetic analyses, we extracted the total genomic DNA from ethanol-preserved femoral muscle tissue of six specimens of the new Thai populations using standard phenol-chloroform-proteinase K extraction procedures with consequent isopropanol precipitation, for a final concentration of about 1 mg/ml (protocols followed Hillis et al. (1996) and Sambrook and David (2001)). We visualized the isolated total genomic DNA in agarose electrophoresis in the presence of ethidium bromide. We measured the concentration of total DNA in 1 μl using NanoDrop 2000 (Thermo Scientific), and consequently adjusted to ca. 100 ng DNA/μL.

We amplified mtDNA fragments of ND2 and its flanking tRNAs using a double-stranded Polymerase Chain Reaction (PCR) to obtain a 1397 base pairs under the following conditions: 1.0 μl genomic DNA (10–30 μg), 1.0 μl light strand primer (concentration 10 μM), 1.0 μl heavy strand primer (concentration 10 μM), 15 μl Master Mix 2x (CWBIO, China), and 12 μl ultra-pure H2O. PCR reactions were executed on Bio-Rad T100™ gradient thermocycler under the following conditions: initial denaturation at 94 °C for 5 min, followed by a second denaturation at 94 °C for 60 s, annealing at 58 °C for 60 s, followed by a cycle extension at 72 °C for 60 s, for 35 cycles with the final extension step 72 °C for 10 min. For amplification of the full ND2 gene with parts of adjacent tRNAs we used the NADH2-Metf6 (5’-AAGCTTTCGGGCCCATACC-3’) and CO1H (5’-AGRGTGCCAATGTCTTTGTGRTT-3’) primers following Macey et al (1997).

PCR products were loaded onto 1.5% agarose gels in the presence of ethidium bromide and visualized in agarose electrophoresis. When distinct bands were produced, we purified the PCR products using 2 μl of a 1:4 dilution of ExoSapIt (Amersham) per 5 μl of PCR product prior to cycle sequencing. A 10 μl sequencing reaction included 2 μL of template, 2.5 μl of sequencing buffer, 0.8 μl of 10 pMol primer, 0.4 μl of BigDye Terminator version 3.1 Sequencing Standard (Applied Biosystems) and 4.2 μl of water. The cyclesequencing used 35 cycles of 10 sec at 96°C, 10 sec at 50°C and 4 min at 60°C. We purified the cycle sequencing products by ethanol precipitation. We carried out sequence data collection and visualization on an ABI 3730xl Automated Sequencer (Applied Biosystems).

Ingroup samples consisted of 14 individuals of the brevipalmatus group representing three nominal species including a sample from one of the new populations in Suan Phueng District plus five individuals of Cyrtodactylus cf. interdigitalis from Thailand and Laos. Four species of the linnwayensis group were used to root the tree following Grismer et al. (2021b). Sequence data for some specimens were acquired from GenBank and newly generated sequences were deposited in GenBank (Table 1).

Maximum Likelihood (ML) and Bayesian Inference (BI) were used to estimate phylogenetic trees. Best-fit models of evolution determined in IQ-TREE (Nguyen et al. 2015) using the Bayesian information criterion (BIC) implemented in ModelFinder (Kalyaanamoorthy et al. 2017) indicated that TNe+I+G4 was the best-fit model of evolution for the tRNAs and HKY+F+G4, TIM3+F+G4, and TPM3u+F+I+G4 were the best models of evolution for codon positions 1, 2, and 3, respectively. The ML analysis was performed using the IQ-TREE webserver (Trifinopoulos et al. 2016) with 1000 bootstrap pseudoreplicates using the ultrafast bootstrap (UFB) analysis (Minh et al. 2013; Hoang et al. 2018). The BI analysis was performed on CIPRES Science Gateway (Miller et al. 2010) using MrBayes v3.2.4 (Ronquist et al. 2012). Two independent runs were performed using Metropolis-coupled Markov Chain Monte Carlo (MCMCMC), each with four chains: three hot and one cold. The MCMCMC chains were run for 15,000,000 generations with the cold chain sampled every 1500 generations and the first 10% of each run being discarded as burn-in. The posterior distribution of trees from each run was summarized using the sumt function in MrBayes v3.2.4 (Ronquist et al. 2012). Stationarity was checked with Tracer v1.7 (Rambaut et al. 2018) to be sure effective sample sizes (ESS) for all parameters were well above 200. We considered Bayesian posterior probabilities (BPP) of 0.95 and above and ultrafast bootstrap support values (UFB) of 95 and above as an indication of strong nodal support (Huelsenbeck et al. 2001; Minh et al. 2013). Uncorrected pairwise sequence divergences (p-distance) were calculated in MEGA 7 (Kumar et al. 2016) using the complete deletion option to remove gaps and missing data from the alignment prior to analysis.

The ML and BI analyses recovered trees with identical topologies (Fig. 2A). Cyrtodactylus elok and C. brevipalmatus of the Thai-Malay Peninsula were recovered as poorly supported sister species (BPP 0.57/UFB 59) and the strongly supported (1.00/100) sister lineage to a strongly supported clade (1.00/100) containing the remaining taxa. Within that clade, the individual of the new population from Suan Phueng District was recovered as the sister lineage to a monophyletic group bearing the phylogenetic sequence of sp. 9, sp. 10, sp. 11, and terminating with the strongly supported (1.00/100) sister species C. ngati and C. cf. interdigitalis. Support values for these nodes in are shown in Figure 2A. The phylogenetic analyses strongly supported the recognition of the Suan Phueng population as a new species, which bears a 15.4–22.1% uncorrected pairwise sequence divergence from the other species in the brevipalmatus group (Table 2). The phylogenetic analyses also supported the separate species status of sp. 9 and sp. 10 as previously recovered by Chomdej et al. (2021).

The phylogeny indicates species recognition for ZMMU R-16492 from Phu Hin Rong Kla National Park, Phitsanulok Province, Thailand (Fig. 1; referred to here as sp. 11) in that it is recovered as the sister species to the strongly supported (1.00/99) lineage comprised of C. ngati and C. cf. interdigitalis (Fig. 2). Although the locality of ZMMU R-16492 at 1147 m in elevation is only 45 km to the northwest of the type locality of C. interdigitalis at 700 m in Nam Nao National Park, Petchabun Province (Ulber 1993), these upland sites are separated by a 20 km wide low-lying river basin of 136 m in elevation at its lowest point with no intervening hilly terrain (Fig. 1). Similarly, C. cf. interdigitalis FMNH 265806 from the Phu Luang Wildlife Sanctuary occurs at 850 m in elevation and is 65 km north of the type locality of C. interdigitalis but separated from it by an approximately15 km wide low-lying river basin at approximately 200 m at its lowest point. ZMMU R-16492 lies 50 km to the southwest of FMNH 265806 and it too is separated by a low-lying river basins. Given the close geographic proximity of these three topographically separated upland localities and the fact that ZMMU R-16492 and FMNH 265806 are not each other’s closest relatives, makes it premature to consider either as conspecific with C. interdigitalis in the absence of genetic data from the type locality or at least morphological data from the type series. Both options are currently being explored.

Although the monophyly of Cyrtodactylus cf. interdigitalis and C. ngati is strongly supported, the monophyly of each of them is not. Although it is difficult to conceive that C. ngati would not be monophyletic given that all specimens were found on the same karst formation approximately 200 km north of the nearest population of C. cf. interdigitalis (Fig. 1), nodal support for its monophyly is only 0.62/82. That of C. cf. interdigitalis ranging from northern Laos to central Thailand is even less—0.51/52. If these nodes were collapsed on the basis of low support values, it would be phylogenetically conceivable to consider C. ngati and C. cf. interdigitalis as conspecific. However, these two species are trenchantly different in morphology (see below) and again, the absence of any data from the type locality of C. interdigitalis, makes any taxonomic proposal premature.

The PCA of the concatenated meristic and adjusted morphometric data sets corroborated the phylogenetic analyses in that the Suan Phueng individuals are not embedded within the cluster of any other species along the combined ordination of principal component (PC) 1 and PC2 (Fig. 2B). PC1 accounted for 36.8% of the variation and loaded most heavily for AG, FemL, HW, HD, EL and FS. PC2 accounted for 15.3% of the variation and loaded most heavily for FL4T and BB (Table 3). These results were mirrored in the DAPC of the concatenated data (Fig. 2C). The PCA of the morphometric data also recovered separation of the Suan Phueng individuals along both axes but very close to C. cf. interdigitalis. In the DAPC, however, the Suan Phueng individuals were within the 66% confidence ellipsoid of C. cf. interdigitalis (Fig. 3). PC1 accounted for 50.2% of the variation and loaded most heavily for ES, HD, AG, FemL, and HL (Fig. 4) and PC2 accounted for an additional 16.2% of the variation, loading most heavily for IO, IN, and EN. The meristic data also recovered separation of the Suan Phueng populations along the combined axes in the PCA and DAPC but again, closest to C. cf. interdigitslis (Fig. 5). PC1 accounted for 35.0% of the variation and loaded most heavily for AG, FemL, HW, HD, EL and FS (Fig. 6) and PC2 accounted for 20.1% of the variation, loading most heavily for FL4E, TL4T, and BB.

Figure 3. 

A. PCA of the Cyrtodactylus brevipalmatus group based on the adjusted morphometric data. B. DAPC of same. C. Violin plots overlain with box plots showing the range, frequency, mean (white dot), and 50% quartile (black rectangle) for each character.

Figure 4. 

A. PCA summary statistics of the adjusted morphometric data of the Cyrtodactylus brevipalmatus group B. Factor loadings of PC1. C. Factor loadings of PC2.

Figure 5. 

A. PCA of the Cyrtodactylus brevipalmatus group based on the meristic data. B. DAPC of same. C. Box plots showing the range, frequency, mean (blue dot), and 50% quartile (black rectangle) for each character. White dots are values on y-axis.

Figure 6. 

A. PCA summary statistics of the meristic data of the Cyrtodactylus brevipalmatus group B. Factor loadings of PC1. C. Factor loadings of PC2.

Based on the phylogenetic relationships, PCA, DAPC, and several discrete morphological differences between the Suan Phueng individuals and all other species of the brevipalmatus group (see comparisons below and Table 4), we hypothesize that they comprise a discretely diagnosable lineage that is not reticulating with any other lineage and as such, should be accorded species status.

Cyrtodactylus rukhadeva sp. nov.

http://zoobank.org/B979D28D-D18F-4B8A-A071-062904D491CC

Figures 7, 8 Suggested common name: tree spirit bent-toed gecko

Cyrtodactylus brevipalmatus Ulber 1993:198 (partim).

Holotype

Adult male ZMMU R-16851 (field tag NAP-09743, tissue sample ID HLM0372) from Thailand, Ratchaburi Province, Suan Phueng District, Khao Laem Mountain (13.53846N, 99.20071E, elevation 994 m a.s.l.), collected by Platon V. Yushchenko and Kawin Jiaranaisakul on 19 June 2019.

Paratype

Adult female ZMMU R-16852 (field tag NAP-09744) from Thailand, Ratchaburi Province, Suan Phueng District, Hoop Phai Tong (13.56210N, 99.20670E, elevation 593 m a.s.l.), collected by Platon Yushchenko on 11 July 2019.

Diagnosis

Cyrtodactylus rukhadeva sp. nov. can be separated from all other species of the brevipalmatus group by having 9–11 supralabials, 10 or 11 infralabials, 27–30 paravertebral tubercles, 19 or 20 rows of longitudinally arranged tubercles, 34–43 transverse rows of ventrals, 152–154 longitudinal rows of ventrals, nine expanded subdigital lamellae on the fourth toe, 11 unexpanded subdigital lamellae on the fourth toe, 18–20 total subdigital lamellae on the fourth toe, eight or nine expanded subdigital lamellae on the fourth finger, nine or 10 unexpanded subdigital lamellae on the fourth finger, 17–19 total subdigital lamellae on the fourth finger, 16–17 enlarged femorals, 20 femoral pores in the male; 17 precloacal pores in the male; 13–17 enlarged precloacals; 16 post-precloacals; enlarged femorals and enlarged precloacals not continuous; proximal femorals not less than one-half the size of the distal femorals; small tubercles on forelimbs and flanks; small dorsolateral caudal tubercles and ventrolateral caudal fringe; paired enlarged subcaudals; and maximum SVL 79.4 mm (Tables 1, 4).

Description of holotype

(Figs 7 and 8) Adult male SVL 74.9 mm; head moderate in length (HL/SVL 0.27), width (HW/HL 0.72), depth (HD/HL 0.46), distinct from neck, triangular in dorsal profile; lores concave slightly anteriorly, weakly inflated posteriorly; prefrontal region slightly concave; canthus rostralis rounded; snout elongate (ES/HL 0.41), rounded in dorsal profile; eye large (ED/HL 0.23); ear opening round, small; eye to ear distance greater than diameter of eye; rostral rectangular, partially divided dorsally by inverted Y-shaped furrow, bordered posteriorly by large left and right supranasals and one small azygous internasal, bordered laterally by first supralabials; external nares bordered anteriorly by rostral, dorsally by large supranasal, posteriorly by two smaller postnasals, bordered ventrally by first supralabial; 11(R,L) rectangular supralabials extending to below midpoint of eye, second supralabial same size as first; 10(R,L) infralabials tapering smoothly to just below and slightly past posterior margin of eye; scales of rostrum and lores flat to domed, larger than granular scales on top of head and occiput; scales of occiput and interorbital region intermixed with distinct, small tubercles; superciliaries subrectngular, largest anteriorly; mental triangular, bordered laterally by first infralabials and posteriorly by large left and right trapezoidal postmentals contacting medially for 60% of their length posterior to mental; one row of slightly enlarged, elongate sublabials extending posteriorly to fifth infralabial; gular and throat scales small, granular, grading posteriorly into slightly larger, flatter, smooth, imbricate, pectoral and ventral scales.

Figure 7. 

A. Dorsal view of the head of the holotype of Cyrtodactylus rukhadeva sp. nov. (ZMMU R-16851) from Khao Laem Mountain, Suan Phueng District, Ratchaburi Province, Thailand. B. Gular region of the holotype. C. Right lateral view of the head of the holotype. D. Pores in the femoral and precloacal region of the holotype. E. Dorsal views of the holotype (left) and paratype (ZMMU R-16852) from Hoop Phai Tong, Suan Phueng District, Ratchaburi Province, Thailand. F. Subcaudal region of the holotype showing the single median row of transversely enlarged subcaudal scales. G. Dorsal view of body showing the tubercles intermixed with granular scales. Photos by Roman A. Nazarov.

Figure 8. 

In life coloration and pattern of (A) the paratype of Cyrtodactylus rukhadeva sp. nov. (ZMMU R-16852) from Hoop Phai Trong, Suan Phueng District, Ratchaburi Province, Thailand and (B) the holotype (ZMMU R-16581) from Khao Laem Mountain, Suan Phueng District, Ratchaburi Province, Thailand. Photos by Platon V. Yushchenko. C. In life coloration and pattern of a hatchling of Cyrtodactylus rukhadeva sp. nov. (not collected); photo by Mali Naiduangchan.

Body relatively short (AG/SVL 0.46) with defined ventrolateral folds; dorsal scales small, granular interspersed with larger, conical, semi-regularly arranged, moderately keeled tubercles; tubercles extend from interorbital region onto base of tail; smaller tubercles extend anteriorly onto nape and occiput, diminishing in size and distinction in interorbital region; approximately 19 longitudinal rows of tubercles at midbody; approximately 27 paravertebral tubercles; 34 longitudinal rows of flat, imbricate, ventral scales much larger than dorsal scales; 153 transverse rows of ventral scales; 17 large, pore-bearing, precloacal scales; no deep precloacal groove or depression; and nine rows of post-precloacal scales on midline.

Forelimbs moderate in stature, relatively short (ForL/SVL 0.11); granular scales of forearm slightly larger than those on body, interspersed with small tubercles; palmar scales rounded, slightly raised; digits well-developed, relatively short, inflected at basal interphalangeal joints; digits narrower distal to inflections; subdigital lamellae wide, transversely expanded proximal to joint inflections, narrower transverse lamellae distal to joint inflections; claws well-developed, claw base sheathed by a dorsal and ventral scale; nine expanded and 11 unexpanded lamellae beneath first finger; hind limbs more robust than forelimbs, moderate in length (TibL/SVL 0.13), covered dorsally by granular scales interspersed with moderately sized, conical tubercles and anteriorly by flat, slightly larger scales; ventral scales of thigh flat, subimbricate, larger than dorsals; subtibial scales flat, imbricate; one row of 9,8(R,L) enlarged femoral scales not continuous with enlarged precloacal scales, terminating distally at knee; proximal femoral scales smaller than distal femorals, the latter forming an abrupt union with smaller, rounded, ventral scales of posteroventral margin of thigh; 8,9(R,L) femoral pores; plantar scales flat; digits relatively long, well-developed, inflected at basal interphalangeal joints; 8(R,L) wide, transversely expanded subdigital lamellae on fourth toe proximal to joint inflection that extend onto sole; nine expanded and 10 unexpanded lamellae beneath first toe; and claws well-developed, sheathed by a dorsal and ventral scale at base.

Tail long (TL/SVL 1.29) 6.8 mm in width at base, tapering to a point; dorsal scales flat, square bearing tubercles forming paravertebral rows and slightly larger tubercles forming a dorsolateral longitudinal row; enlarged, posteriorly directed semi-spinose tubercles forming small but distinct ventrolateral caudal fringe; median row of paired, transversely expanded subcaudal scales, larger than dorsal caudal scales; base of tail bearing hemipenal swellings; three conical postcloacal tubercles at base of hemipenal swellings; and postcloacal scales flat, imbricate.

Coloration in life (Fig. 8)

Ground color of the head body, limbs, and tail brown; faint, diffuse mottling on the top of the head; wide, dark-brown post-orbital stripe; lores dark- brown; nuchal band faint, bearing two posterior projections; three wide faint body bands edged in dark-brown between limb insertions; band interspaces bearing irregularly shaped, dark markings; faint dark speckling on limbs and digits; five wide slightly darker caudal bands separated by five lighter caudal bands; posterior tip of regenerated tail unbanded; ventral surfaces beige, generally immaculate; gular region lacking pigment; poaterior subcaudal region faintly banded; iris orangish gold in color.

Variation (Fig. 8)

The female paratype differs significantly from the male holotype in coloration and pattern. The holotype has a much less contrasting dorsal pattern with only faint indications of dorsal body and caudal bands whereas the paratype is boldly marked, bearing a one dark jagged-edged nuchal band and three dark jagged-edged body bands, all with darkened borders. The caudal pattern of the paratype consists of six dark and seven light-colored contrasting bands on a full original tail whereas the holotype has four dark and light-colored weakly contrasting bands and a regenerated tail tip. The top of the head of the holotype is essentially unicolor whereas that of the paratype is darkly mottled. With only one specimen of each sex we cannot say if these differences in color pattern are sexually dimorphic. An uncataloged juvenile had a generally unicolor tan dorsal ground color overlain with a dark-brown nuchal loop bearing two posteriorly oriented projections, four complete to incomplete dark-brown body bands, and nine dark-brown caudal bands separated by eight tan to white caudal bands (Fig. 8)

Distribution

Cyrtodactylus rukhadeva sp. nov. is known only from the type locality from Khao Laem Mountain, Suan Phueng District, Ratchaburi Province, Thailand and Hoop Phai, Suan Phueng District, Ratchaburi Province, approximately 7.7 km to the west of the type locality (Fig. 1).

Etymology

The specific epithet “rukhadeva” is given as a noun in apposition and refers to the spirits or gods residing in trees in Thai mythology, known as Rukha Deva (literally “Tree Nymphs”). According to Thai folklore, these sylvan spirits live on tree branches and on large older trees wearing traditional Thai attire, usually in reddish or brownish colours, and are believed to protect the forest. The new arboreal species of Cyrtodactylus resides in one of the remaining fragments of the north Tenasserim montane forests. We want to underscore the need for the immediate assessment of herpetofaunal diversity surveys and implementation of adequate conservation measures for these relic forests.

Comparisons

Given the low sample size (n=2) of Cyrtodactylus rukhadeva sp. nov., meaningful statistical comparisons of meristic were not possible. However, some character ranges—at this point—are diagnostic in that they are widely separated from one another. We are well aware that a larger sample size may preclude some of these characters as being diagnostic but it is equally probable that they will demonstrate that they are truly diagnostic. Cyrtodactylus rukhadeva sp. nov. can be separated from C. ngati in having fewer paravertebral tubercles (PVT; 27–30 vs 38–40), from C. elok in having more longitudinal rows of dorsal tubercles (LRT; 19 or 20 vs 4–7), from C. brevipalmatus and C. elok in having fewer transverse rows of ventral scales (VSM; 152–154 vs 172–180 and 190–234, respectively), from C. brevipalmatus in having more enlarged femoral scales (FS; 16 or 17 vs 10–16), from all species in having more femoral pores in males (FP; 20 vs 0–17 collectively), from C. cf. interdigitalis, C. brevipalmatus, and C. elok in having males having more precloacal pores (PP; 17 vs 13, 7, and 8, respectively), and from all species except C. elok in having more post-precloacal scales (PPS; 16 vs 2–10, cllectively. All ranges are presented in Table 5. Summary statistics and comparisons among the adjusted morphometric characters are presented in Table 6.

Table 5.

Download as

CSV

XLSX

Summary statistics of meristic characters. sd = standard deviation and n= sample size. * = males only. Shaded cells denote diagnostic character differences from C. rukhadeva sp. nov. Abbreviations are in the Materials and methods.

rukhedeva sp. nov.	SL	IL	PVT	LRT	VS	VSM	TL4E	TL4U	TL4T	FL4E	FL4U	FL4T	FS	FP*	PP*	PCS	PPS	BB
mean	10	10.5	28.5	19.5	38.5	153	9	11	19	8.5	9.5	18	16.5	20	17	15	16	3
±1 sd	1.41	0.71	1.41	0.71	6.36	1.41	0	00.00	1.41	0.71	0.71	1.41	0.71	0.00	0.00	2.83	0	0
minimum	9	10	27	19	34	152	9	11	18	8	9	17	16	20	17	13	16	3
maximum	11	11	30	20	43	154	9	11	20	9	10	19	17	20	17	17	16	3
n	2	2	2	2	2	2	2	2	2	2	2	2	2	1	1	2	2	2
ngati
mean	10	9	39.5	18.8	35.3	167.0	9.0	10.3	19.25	7.25	8.5	15.8	18.3	14	0	12.3	2	6
±1 sd	0	0	1.32	2.22	2.50	8.41	0.82	0.58	1.73	0.50	0	1.5	1.5	0.00	0	0.5	0	0
minimum	10	9	38	17	32	158	8	9	18	7	8	15	17	14	0	12	2	6
maximum	10	9	40	22	38	178	10	11	20	8	9	18	20	14	0	13	2	6
n	4	4	4	4	4	4	4	4	4	4	4	4	4	1	1	4	4	4
cf. interdigitalis
mean	12	10.1	28.1	19.1	35.1	162.0	9.3	10.9	20.1	8.1	9.6	17.7	16.3	14.0	13	13.1	7.1	3.1
±1 sd	2.05	1.28	1.77	2.14	1.51	6.87	0.89	0.76	1.16	0.64	0.74	0.89	1.06	1.00	0.00	0.35	1.49	0.35
minimum	9	8	26	17	33	146	8	10	18	7	8	16	15	13	13	12	5	3
maximum	14	12	32	24	37	166	10	12	21	9	10	19	18	15	13	13	10	4
n	7	7	7	7	7	7	7	7	7	7	7	7	7	2	2	7	7	7
brevipalmatus
mean	11	9	38.0	16	38	176	8	11.7	19.7	8.0	10.0	18.0	12.3	7	7	7.0	10.3	4.3
±1 sd	1.00	1.00	1.00	1.00	0	6	1.00	1.15	0.58	0	1.00	1.00	3.21	0.00	0.00	0	0.58	1.53
minimum	10	8	37	15	38	170	7	11	19	8	9	17	10	7	7	7	10	3
maximum	12	10	39	17	38	182	9	13	20	8	11	19	16	7	7	7	11	6
n	3	3	3	3	3	3	3	3	3	3	3	3	3	1	1	3	3	3
elok
mean	10.3	9.8	0	5.3	43.3	210.3	9.3	10.3	19.5	8.8	10.5	19.3	0	0	8	7.8	13.8	4.0
±1 sd	2.22	1.50	0	1.50	4.92	22.54	0.50	0.96	1.29	0.50	2.38	2.22	0	0	0.00	0.50	3.40	1.15
minimum	8	8	0	4	36	190	9	9	18	8	8	17	0	0	8	7	11	3
maximum	13	11	0	7	47	234	10	11	21	9	13	22	0	0	8	8	18	5
n	4	4	4	4	4	4	4	4	4	4	4	4	4	4	1	4	4	4
sp. 11
mean	8	9	30	18	34	160	9	10	19	10	9	19	17	17	13	13	8	3
±1 sd	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0
minimum	8	8	30	18	34	160	9	10	19	10	9	19	17	17	13	13	8	3
maximum	8	9	30	18	34	160	9	10	19	10	9	19	17	17	13	13	8	3
n	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1

Table 6.

Download as

CSV

XLSX

Summary statistics of adjusted morphometric characters. sd = standard deviation and n= sample size. Abbreviations are in the Materials and methods.

rukhadeva sp. nov.	SVL	AG	HumL	ForeaL	FemurL	TibL	HL	HW	HD	ED	EE	ES	EN	IO	EL	IN
mean	1.865	1.526	1.023	0.916	1.086	0.987	1.294	1.146	0.947	0.648	0.792	0.903	0.778	0.505	0.040	0.332
±1 sd	0.0134	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
minimum	1.8555	1.5263	1.0233	0.9163	1.0863	0.9866	1.2944	1.1459	0.9468	0.6483	0.7924	0.9030	0.7778	0.5051	0.0400	0.3324
maximum	1.8745	1.5263	1.0233	0.9163	1.0863	0.9866	1.2944	1.1459	0.9468	0.6483	0.7924	0.9030	0.7778	0.5051	0.0400	0.3324
n	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2
ngati
mean	1.832	1.471	0.912	0.990	1.061	1.050	1.310	1.084	0.849	0.575	0.758	0.865	0.806	0.743	-0.116	0.426
±1 sd	0.0090	0.0017	0.0043	0.0074	0.0000	0.0059	0.0005	0.0027	0.0096	0.0350	0.0157	0.0153	0.0042	0.0091	0.0192	0.0078
minimum	1.8228	1.4702	0.9070	0.9846	1.0607	1.0435	1.3090	1.0806	0.8422	0.5513	0.7401	0.8546	0.8010	0.7325	-0.1292	0.4166
maximum	1.8407	1.4732	0.9148	0.9980	1.0607	1.0543	1.3100	1.0855	0.8596	0.6148	0.7687	0.8824	0.8086	0.7489	-0.0943	0.4308
n	3	3	3	3	3	3	3	3	3	3	3	3	3	3	3	3
cf. interdigitalis
mean	1.883	1.568	0.954	0.998	1.129	1.080	1.318	1.121	0.934	0.706	0.787	0.926	0.790	0.683	0.051	0.390
±1 sd	0.0361	0.0249	0.0617	0.0225	0.0137	0.0181	0.0123	0.0167	0.0259	0.0574	0.0329	0.0094	0.0138	0.1168	0.0636	0.0332
minimum	1.8331	1.5100	0.8481	0.9638	1.1106	1.0549	1.2949	1.0981	0.8869	0.6455	0.7191	0.9134	0.7688	0.4949	-0.0339	0.3367
maximum	1.9400	1.5854	1.0400	1.0268	1.1482	1.1039	1.3351	1.1458	0.9648	0.8125	0.8257	0.9384	0.8098	0.7966	0.1584	0.4354
n	8	8	8	8	8	8	8	8	8	8	8	8	8	8	8	8
brevipalmatus
mean	1.831	1.519	0.953	0.979	1.080	1.058	1.283	1.117	0.888	0.668	0.740	0.865	0.727	0.692	0.063	0.304
±1 sd	0.0222	0.0235	0.0208	0.0092	0.0055	0.0017	0.0011	0.0027	0.0125	0.0384	0.0046	0.0007	0.0200	0.0338	0.0650	0.0579
minimum	1.8069	1.5000	0.9368	0.9718	1.0743	1.0559	1.2825	1.1141	0.8783	0.6372	0.7360	0.8645	0.7106	0.6653	-0.0098	0.2393
maximum	1.8500	1.5452	0.9768	0.9895	1.0848	1.0593	1.2846	1.1194	0.9023	0.7110	0.7449	0.8659	0.7490	0.7303	0.1151	0.3504
n	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4
elok
mean	1.905	1.589	0.803	1.053	1.136	1.127	1.337	1.195	0.996	0.709	0.834	0.947	0.794	0.703	0.187	0.434
±1 sd	0.0161	0.0121	0.3410	0.0259	0.0223	0.0237	0.0010	0.0116	0.0082	0.0171	0.0199	0.0078	0.0101	0.0699	0.0619	0.0327
minimum	1.8932	1.5730	0.3367	1.0175	1.1118	1.0955	1.3363	1.1811	0.9833	0.6836	0.8063	0.9361	0.7790	0.6088	0.1532	0.3965
maximum	1.9284	1.6007	1.0886	1.0790	1.1629	1.1529	1.3387	1.2094	1.0008	0.7201	0.8527	0.9544	0.8001	0.7583	0.2797	0.4762
n	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4
sp. 11
mean	1.940	1.571	1.031	0.979	1.145	1.071	1.311	1.139	0.921	0.645	0.793	0.915	0.788	0.508	0.015	0.337
±1 sd	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
minimum	1.940	1.571	1.031	0.979	1.145	1.071	1.311	1.139	0.921	0.645	0.793	0.915	0.788	0.508	0.015	0.337
maximum	1.940	1.571	1.031	0.979	1.145	1.071	1.311	1.139	0.921	0.645	0.793	0.915	0.788	0.508	0.015	0.337
n	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1

Cyrtodactylus rukhadeva sp. nov. can be readily separated from all other species on the basis of discrete morphological differences (Table 4). Cyrtodactylus rukhadeva sp. nov. has paravertebral tubercles and enlarged femoral scales which are absent in C. elok. Cyrtodactylus rukhadeva sp. nov. has small tubercles on the forelimbs which are absent in C. cf. interdigitalis, C. ngati, and C. elok. Cyrtodactylus rukhadeva sp. nov. has small tubercles on the flanks, femoral pores in males, enlarged dorsolateral tubercles, and a large ventrolateral caudal fringe that are absent in C. elok. Cyrtodactylus rukhadeva sp. nov. has a square tail in cross-section versus a round tail in C. cf. interdigitalis, C. ngati, and C. brevipalmatus. Cyrtodactylus rukhadeva sp. nov. has single enlarged subcaudal scales versus paired enlarged subcaudal scales in C. cf. interdigitalis and C. brevipalmatus. Cyrtodactylus elok and C. ngati lack distinctly enlarged subcaudal scales. Cyrtodactylus rukhadeva sp. nov. has a prehensile tail apparently lacking in C. ngati. Figure 9 illustrates the general morphological and color pattern similarity of the other members of the brevipalmatus group to C. rukhadeva sp. nov.

Figure 9. 

Illustration showing the general similarity body shape and color pattern among five species of the Cyrtodactylus brevipalmatus group. A. Adult female C. elok (LSUDPC 2589) from Negeri Sembian, Peninsular Malaysia. Photo by L. Lee Grismer. B Adult male sp. 11 (ZMMU R-16492) from Phu Hin Rong Kla National Park, Petchabun Province, Phitsanulok, Thailand. Photo by Nikolay A. Poyarkov. C. Adult female C. interdigitalis from Tak Mal National Park, Phetchabun Province, Thailand. Photo from Creative Commons Attribution Share Alike. Photo by Evan S. H. Quah. D. Adult female C. cf. interdigitalis (NCSM 79472) from Ban Pha Liep, Houay Liep Stream, Xaignabouli Province, Laos. Photo by Bryan L. Stuart. E. Adult female C. ngati (VNUF R.2020.12) from Pa Thom Cave, Pa Xa Lao Village, Pa Thom Commune, Dien Bien District, Dien Bien Province, Vietnam. Photo by Dzung T. Le. F. Adult female C. brevipalmatus from the type locality at Khao Luang National Park, Nakon Si Thammarat Provice, Thailand. Photo from Creative Commons Attribution Share Alike. G. Adult male C. cf. ruhkhadeva sp. nov. form Kaeng_Krachan National Park, Petchaburi Province, Thailand. Photo from Creative Commons Attribution Share Alike.

Natural History

Cyrtodactylus rukhadeva sp. nov. inhabits montane and submontane evergreen polydominant tropical forests, and was most often recorded in bamboo forests mixed with dry dipterocarp forests at elevations ranging from 600 to 1100 m a.s.l. and dominated by the large trees Anisoptera costata Korth, 1841 and Dipterocarpus gracilis Blume, 1825 (Dipterocarpaceae) (Fig. 10A). Lizards were usually observed near large tree holes or hiding among the roots of strangler fig trees (Ficus sp.; Moraceae) (Fig. 10B). All were active from approximately 19:30–23:50 hrs and approximately 2 m above the ground. One lizard was observed ca. 6 m above the ground on the branch of Syzygium sp. (Myrtaceae) above a rocky streambed. Given its morphology, we believe C. rukhadeva sp. nov. is an arboreal specialist that generally resides in the upper canopy, but may be forced down to lower portions of the canopy or to the understory layer by rain storms or strong winds. Gravid females bearing two eggs were recorded in March; two eggs of the new species were found attached in a tree hole; total length of the hatchlings was approximately 7 cm. Cyrtodactylus rukhadeva sp. nov. is a highly aggressive species that opens its mouth and waves its tail from side to side in a defensive posture when threatened. At the type locality, C. rukhadeva sp. nov. was recorded in sympatry with Cnemaspis selenolagus Grismer, Yushchenko, Pawangkhanant, Nazarov, Naiduangchan, Suwannapoom & Poyarkov, 2020, Cyrtodactylus cf. oldhami (Theobald, 1876), and Gekko (Ptychozoon) kaengkrachanense (Sumontha, Pauwels, Kunya, Limlikhitaksorn, Ruksue, Taokratok, Ansermet & Chanhome, 2012).

Figure 10. 

Habitat of Cyrtodactylus rukhadeva sp. nov. at the type locality: Khao Laem Mountain, Suan Phueng District, Ratchaburi Province, Thailand. (A) Submontane forest on the slopes of Khao Laem Mountain. (B) A male of the new species hiding among the roots of a strangler fig (Ficus sp.). Photos by Parinya Pawangkhanant.