A phylogenetic and taxonomic assessment of the Cnemaspis alwisi group (Reptilia: Gekkonidae) in Sri Lanka with a description of two new species from isolated misty-mountains

a


Introduction
Sri Lanka, along with the Western Ghats of India, is ranked as one of the world's smallest biodiversity hotspots (Meegaskumbura et al. 2002).Of the 242 species of reptiles described, ~66% are endemic to the country, and thus Sri Lanka is also considered a reptile diversity hotspot (Roll et al. 2017;Karunarathna et al. 2020).Within this rich reptile assemblage, the diversity of geckos (Family Gekkonidae) is remarkable; 63 species (from eight genera) have been described so far, accounting for ~26% of the overall reptilian species-richness (de Silva et al. 2019;Amarasinghe and Karunarathna 2020;Amarasinghe et al. 2021b).Of these, 53 species (~84%) are endemic and 48 species (~76%) are threatened and vulnerable to extinction due to ongoing climate change effects, irresponsible development, loss of good quality habitats, and forest fragmentation (Amarasinghe and Karunarathna 2020;Dayananda et al. 2021).
Despite the recent resurrection of the genus Ancylodactylus Müller, 1907 to accommodate African species formerly assigned to Cnemaspis Strauch, 1887 (Malonza and Bauer 2022), more than 200 species of Cnemaspis are currently recognized, making it the second most speciose gekkonid genus in the world after Cyrtodactylus Gray, 1827 (Grismer et al. 2021;Uetz et al. 2022).The genus Cnemaspis comprises 40 species in Sri Lanka, all of which are endemic (Manamendra-Arachchi et al. 2007;Wickramasinghe et al. 2016;Batuwita et al. 2019;Karunarathna et al. 2021).Sri Lankan Cnemaspis represent about ~20% of global Cnemaspis richness despite occupying a small fraction of the global distribution area of the genus.Recent progress in studying the herpetofauna of the island has resulted in the discovery of many new species of Cnemaspis (de Silva et al. 2019;Karunarathna et al. 2019aKarunarathna et al. , 2019bKarunarathna et al. , 2019c) of which 50% have been described in the last decade.Most (~73%) of the Cnemaspis species in Sri Lanka occur within the wet bioclimatic zone (annual rainfall ≥ 2500 mm).However, the most recent species to be described from Sri Lanka (~27%) were recorded from small isolated forest areas or hillocks in the dry bioclimatic zone (annual rainfall between 1000-2000 mm) (Batuwita et al. 2019;Karunarathna and Ukuwela 2019;Vidanapathirana et al. 2014).
Sri Lankan Cnemaspis species represent two distinct evolutionary lineages, the kandiana and podihuna clades (Agarwal et al. 2017;Karunarathna et al. 2021;Pal et al. 2021).The species of the kandiana clade are morphologically characterized by small and irregularly shaped subcaudal scales, and the species of the podihuna clade are characterized by clearly enlarged, hexagonal or subhexagonal subcaudal scales (Karunarathna et al. 2019c).The high Cnemaspis species richness in Sri Lanka may be ac-counted for by possible multiple colonizations from the Indian mainland with isolated in-situ speciation and radiation (Agarwal et al. 2017(Agarwal et al. , 2020;;Pal et al. 2021).Most of the Sri Lankan Cnemaspis are point-endemics with distribution ranges limited to very small land areas.It is very likely that future studies on the taxonomy and biogeography of Cnemaspis in Sri Lanka will highlight the importance of these isolated habitats in generating and maintaining the diversity of these unique groups of geckos in the island (Deraniyagala 1944;Bauer et al. 2007;Wickramasinghe and Munindradasa 2007;Batuwita and Udugampola 2017).During recent fieldwork in lowland dry plains of northwestern and northeastern Sri Lanka, we discovered two unidentified populations of day geckos, which closely resemble species of the alwisi (comprising C. alwisi Wickramasinghe & Munindradasa, 2007;C. gunasekarai Amarasinghe et al., 2021a;C. hitihamii Karunarathna et al., 2019c;C. kohukumburai Karunarathna et al., 2019c;C. nilgala Karunarathna et al., 2019a;C. punctata Manamendra-Arachchi et al., 2007) and rammalensis (comprising C. gunawardanai Amara singhe et al., 2021a;C. rajakarunai Wickramasinghe et al., 2016;C. rammalensis Vidanapathirana et al., 2014) groups (Amarasinghe et al. 2021a) in the podihuna clade, yet had apparent morphological differences from all currently known species.Thus, on the basis of an integrative taxonomic approach, we describe them below as two new species.

Field sampling and specimens
We conducted field surveys in 173 different locations distributed across several bioclimatic regions (e.g.Semi-arid zone, dry zone, intermediate zone and wet zone) in Sri Lanka as a part of an on-going island-wide survey of lizards under permit number WL/3/2/42/18 (a-d), issued by the Department of Wildlife Conservation (DWC) and permit number R&E/RES/NFSRCM/2019-04, issued by the Forest Department of Sri Lanka to KDBU, ADS and SK.At each location, we surveyed and documented gecko species recorded with special attention to the focal genus Cnemaspis.On average, per location, we spent 12 man-hours per survey.Museum acronyms follow Uetz et al. (2019) and Sabaj (2020).The type materials reported in this paper is deposited in the National Museum of Sri Lanka (NMSL), Colombo.Specimens were caught by hand and photographed in life.They were subsequently euthanised using halothane and fixed in 10% formaldehyde for two days, washed in water and transferred to 70% ethanol for long-term storage.Tail tips were collected as tissue samples before fixation and were stored in 95% ethanol.For comparison, we examined 479 Cnemaspis specimens (catalogued and uncatalogued) representing all recognised Sri Lankan species, including all type specimens housed at the National Museum of Sri Lanka (NMSL), The Natural History Museum, London (BMNH), Museum of the Department of Wildlife Conservation, Giritale, Sri Lanka (DWC) and specimens collected by Anslem de Silva (bearing the field codes ADS), Aaron Bauer (bearing the field codes AMB) and Suranjan Karunarathna (bearing the field codes SSK and DMSSK), which have been deposited in the NMSL (see Appendix 1).Specimens that formerly belonged to the Wildlife Heritage Trust (WHT) collection, which bear WHT numbers, are currently deposited at the NMSL, and were catalogued under their original numbers.Additional information on the morphology and natural history of Sri Lankan Cnemaspis species was obtained from the relevant literature (Bauer et al. 2007;Manamendra-Arachchi et al. 2007;Wickramasinghe and Munindradasa 2007;Vidanapathirana et al. 2014;Wickramasinghe et al. 2016;Agarwal et al. 2017;Batuwita and Udugampala 2017;Batuwita et al. 2019;de Silva et al. 2019;Karunarathna et al. 2019aKarunarathna et al. , 2019bKarunarathna et al. , 2019c;;Karunarathna and Ukuwela 2019;Amarasinghe and Karunarathna 2020;Karunarathna et al. 2021;Amarasinghe et al. 2021a).Assignment of unidentified specimens to the new species was based on their morphometric and meristic characters, colour patterns and the level of geographic isolation, and the species limits were tested independently by molecular phylogenetic methods.These new species have not been included in previously published phylogenies (see Agarwal et al. 2017;Karunarathna et al. 2019cKarunarathna et al. , 2021)).

Morphometric characters
Forty morphometric measurements were taken from the two new populations of Cnemaspis from Ethagala mountain and Galgiriya mountain and representative individuals of the members of the alwisi and rajakarunai groups using a Mitutoyo digital Vernier calliper (to the nearest 0.1 mm).Detailed observations of scales and other structures were made through Leica Wild M3Z and Leica EZ4 dissecting microscopes.The following symmetrical meristic characters were taken on the left side of the body: eye diameter (ED, horizontal diameter of eye ball); orbital diameter (OD, the greatest diameter of orbit); eye to nostril length (EN, the distance between the anteriormost point of the orbit and the posterior border of nostril); snout length (ES, the distance between the anteriormost point of the orbit and the tip of snout); snout to nostril length (SN, the distance between the tip of snout and the anteriormost point of nostril); nostril width (NW, the maximum horizontal width of nostril); eye to ear distance (EE, the distance between the posterior border of eye and the anteriormost point of ear opening); snout to axilla distance (SA, the distance between axilla and the tip of snout); ear length (EL, the maximum length of ear opening); interorbital width (IO, the shortest distance between the left and right supraciliary scale rows); inter-ear distance (IE, the distance across the head between the two ear openings); head length (HL, the distance between the posterior edge of mandible and the tip of snout); head width (HW, the maximum width of head in-between the ears and the orbits); head depth (HD, the maximum height of head at the eye level); jaw length (JL, the distance between the tip of snout and the corner of mouth); internarial distance (IN, the smallest distance between the inner margins of nostrils); snout to ear distance (SED, the distance between the tip of snout and the anteriormost point of ear); upper-arm length (UAL, the distance between axilla and the angle of elbow); lower-arm length (LAL, the distance from elbow to wrist with palm flexed); palm length (PAL, the distance between wrist (carpus) and the tip of longest finger excluding the claw); lengths of digits I-V of manus (DLM, the distance between the juncture of the basal phalanx with the adjacent digit and the tip of the digit, excluding the claw); snout-vent length (SVL, the distance between the tip of snout and the anterior margin of vent); trunk length (TRL, the distance between axilla and groin); trunk width (TW, the maximum width of body); trunk depth (TD, the maximum depth of body); femur length (FEL, the distance between groin and knee); tibia length (TBL, the distance from knee to heel with ankle dorsiflexed); heel length (HEL, the distance between ankle (tarsus) and the tip of longest toe [excluding the claw] with both foot and tibia flexed); lengths of pedal digits I-V (DLP, the distance between the juncture of the basal phalanx with the adjacent digit and the digit tip, excluding the claw); tail length (TAL, the distance between the anterior margin of vent and the tail tip); tail base depth (TBD, the maximum height of the tail base); tail base width (TBW, the widest point of the tail base).

Meristic characters
Thirty discrete characters were observed and recorded using Leica Wild M3Z and Leica EZ4 dissecting microscopes on both the left (L) and the right (R) side of the body (reported in the L/R order): number of supralabials (SUP) and infralabials (INF) between the first labial scale and the corner of the mouth; number of interorbital scales (INOS) between the left and right supraciliary scale rows; number of postmentals (PM) bordered by chin scales, 1 st infralabial on the left and right and the mental; number of chin scales (CHS) touching medial edge of infralabials and mental between juncture of 1 st and 2 nd infralabials on the left and right; number of supranasal (SUN) scales between nostrils; presence of the postnasal (PON) scales posterior to the naris; presence of the internasal (INT) scale between supranasals; number of supraciliary scales (SUS) above the eye; number of scales between the eye and tympanum (BET) from posterior-most point of the orbit to anterior-most point of the tympanum; number of canthal scales (CAS), number of scales from posteri-or-most point of naris to anteriormost point of the orbit; total lamellae on manus I-V (TLM) counted from first proximal enlarged scansor, greater than twice the width of the largest palm scale, to distalmost lamella at tip of digits; number of dorsal paravertebral granules (PG) between pelvic and pectoral limb insertion points along a straight line immediately left of the vertebral column; number of mid-body scales (MBS) from the centre of mid-dorsal row diagonally towards the ventral scales; number of mid-ventral scales (MVS) from the first scale posterior to the mental to last scale anterior to the vent; number of belly scales (BLS) across the ventre between the lowest rows of granular dorsal scales; total lamellae on pes I-V (TLP), counted from first proximal enlarged scansor greater than twice the width of the largest heel scale, to distalmost lamella at tip of digits; number of femoral pores (FP) present on the femur; number of non-pored interfemoral scales (IFS), counted between pore-bearing femoral rows; numbers of non-pored distal femoral scales (DFS) counted from distal ends of femoral pore rows to knee.In addition, we also evaluated the texture [keeled (KD) or smooth (SM)] of the ventral scales, the texture [heterogeneous (HET) or homogeneous (HOM)] of the dorsal scales, the number of spinous scales on the flanks (FLSP) and characteristics, such as appearance of the caudal scales (except in specimens with regenerated tails).Colouration was determined from digital images of living specimens and also from direct observations in the field.

Distribution and natural history
During the surveys, behavioural and other aspects (e.g.habitat, microhabitat) of natural history of the focal species were observed through opportunistic and non-systematic means.Our surveys covered a wide variety of habitats, including open woodlands, primary and secondary forests, pine plantations, home gardens, and various agricultural lands (coconut, rubber, tea, cardamom, cinnamon, oil palm plantations, and paddy fields).Observations of geckos were made by the naked eye at a distance of 2 to 3 m without making any disturbance.
The ambient temperature and the substrate temperature were measured using a standard thermometer and a N19 Q1370 infrared thermometer (Dick Smith Electronics, Shanghai, China), respectively.The relative humidity and light intensity were measured with a QM 1594 multifunction environment meter (Digitek Instruments Co., Ltd., Hong Kong, China).To record elevation and georeference species locations, an eTrex 10 GPS (Garmin) was used (map datum WGS1984).Sex was determined by the presence of hemipenial bulges and femoral pores in males (M) or absence of the above in females (F).

DNA Based Species Delimitation
To determine genetic distinctiveness of the two new species from the already known Cnemaspis species, we examined the mitochondrial NADH dehydrogenase subunit 2 gene (ND2) including the adjacent tRNA regions.ND2 gene is commonly used as a genetic marker for geckos and other reptiles and a majority of DNA sequences available in GenBank for Cnemaspis species are of this gene.Genetic distinctiveness was determined by examining the haplotype clusters through phylogenetic analysis (Wiens and Penkrot 2002), estimation of the uncorrected pairwise genetic distances among taxa, and species delimitation analysis.

DNA isolation, PCR amplification and sequencing
Total genomic DNA was extracted from ethanol-preserved tail tissue samples using standard phenol-chloroform-proteinase K extraction procedures with consequent isopropanol precipitation (protocols followed Hillis et al. 1996, andSambrook andRussell 2001).The isolated total genomic DNA was visualized in 1% agarose gel electrophoresis in the presence of ethidium bromide.The concentrations of total DNA in the extracts were determined using a NanoDrop 2000 Nano-spectrophotometer (Thermo Scientific, USA), and consequently adjusted to ca. 100 ng DNA/μL.
PCR products were visualized through ethidium bromide stained 1% agarose gel electrophoresis and purified using 2 μl of a 1:4 dilution of ExoSapIt (Amersham, UK) per 5 μl of PCR product prior to cycle sequencing.Successful targeted PCR products were outsourced to Genetech Sri Lanka Pvt. Ltd.Colombo for purification and bidirectional sequencing.Consensus sequences from forward and reverse reads were assembled in GENEIOUS PRO 5.6 software (Drummond et al. 2009).The obtained sequences were deposited in GenBank under the accession numbers OL657119-OL657164 and OL741433 (Table S1).

Phylogenetic analyses
We downloaded most of the available ND2 sequences for Cnemaspis species of the South Asian radiation (Table S1) from GenBank.However, we did not include Cnemaspis sensu stricto from Southeast Asia as they are known to form a separate clade, unrelated from the South Asian Cnemaspis species except for C. modiglianii Das, 2005, C. tanintharyi Lee et al., 2019, and C. thayawthadangyi Lee et al., 2019 which are known to be nested within the South Asian radiation (Agarwal et al. 2017;Lee et al. 2019).Additionally, we included three Cnemaspis (C.alwisi [n = 8], C. molligodai [n = 3] and C. kandambyi [n = 1]) species generated in this study that have not been included in previous phylogenies.We used Calodactylodes illingworthorum Deraniyagala, 1953 as an outgroup to root the tree as it has been shown that Calodactylodes Strand, 1928 is the sister lineage to the South Asian Cnemaspis lineage (Agarwal et al. 2017).In total, we analysed 142 sequences of the fragment of mitochondrial ND2 gene and flanking tRNA genes.
Nucleotide sequences were initially aligned in GENEIOUS PRO 5.6 software (Drummond et al. 2009) with default parameters, and subsequently checked by eye and slightly adjusted.The final alignment contained 1328 bp.Partitioning schemes and optimal substitutional models for the alignment were estimated using the BIC criterion implemented in PartitionFinder (Lanfear et al. 2012).BIC indicated four partitions; each codon posi tion of the ND2 gene and the adjacent tRNA region.GTR+I+G substitution model was selected for each of the codon positions of the ND2 gene while GTR was selected for adjacent tRNA region.
We inferred the matrilineal genealogy using Bayesian Inference (BI) and Maximum Likelihood (ML) approaches.Partitioned Bayesian analysis was performed in MrBayes 3.2.6 (Ronquist and Huelsenbeck 2003) with unlinked model parameters using default priors for 80 million generations with two independent runs and four chains (one hot and three cold chains), sampling every 10000 generations.Convergence of the independent runs was assessed by examining split frequencies (<0.01) of clades across runs, effective sample sizes (ESS values) and likelihood plots in Tracer v1.4.1 (Rambaut et al. 2018).An all compatible consensus tree was built after first 25% of sampled trees were discarded as burn-in.We assessed confidence in tree topology by the frequency of nodal resolution (posterior probability; BI PP) (Huelsenbeck and Ronquist 2001).A partitioned ML analysis was implemented in RAxML 7.2.6.(Stamatakis et al. 2008) with 200 independent ML searches using the rapid hill climbing algorithm.Analysis involved same partitioning schemes and with GTR+G as the substitution model for each partition.Branch support was estimated using 1000 bootstrap (BS) pseudoreplicates.In both analyses, we regarded tree nodes with BI PP over 0.95 and BS values over 70 to be sufficiently resolved.Bayesian inference PP values between 0.95 and 0.90 were regarded as tendencies.Lower values were considered to indicate unresolved nodes (Huelsenbeck and Hillis 1993;Minh et al. 2013).Mean uncorrected genetic distances (p-distances) between sequences of mitochondrial ND2 gene (excluding the flanking tRNA) were calculated in MEGA X with an average site cut-off of 95% (Kumar et al. 2018).
Species delimitation analysis using Poisson Tree Process (PTP) (Zhang et al. 2013) was conducted using the rooted BI tree with 142 taxa as an input tree.The cal-culations were performed on the PTP web server (http:// species.h-its.org/ptp),with 500,000 MCMC generations, thinning set to 100 and burn-in set at 25%.The outgroup (C.illingworthorum) was removed from the tree, prior to the analysis to improve its performance.The probability of each node representing a species was calculated with both Maximum Likelihood (PTP) and Bayesian (bPTP) models.

Statistical analyses
All statistical analyses were conducted using R Core Team (2018).The 40 morphometric, 26 meristic, and six categorical color pattern characters used in statistical analyses are listed in (Table S2).To remove the effects of allometry in the morphometric characters (sec.Chan and Grismer 2022), 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(Thorpe , 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 potential intra-with interspecific variation (Reist 1986;McCoy et al. 2006).Precloacal and femoral pores were omitted from the multivariate analyses due to their absence in females.
Based on Levene's test for the normalized morphometric and meristic characters, equal variances were recovered for all species.Characters were then analyzed using an analysis of variance (ANOVA) and TukeyHSD post hoc test to search for statistically significant mean differences between all combinations of species pairs (Table S3).A principal coordinate analysis (PCoA) using a Gower (dis)similarity index was employed on a concatenated meristic and normalized morphometric data set.The dissimilarity matrix, not the original variables themselves, is used as the input to the analysis.Therefore, information concerning the original variables cannot be recovered.Because this unsupervised multivariate analysis is based on a (dis)similarity index, it is appropriate for data sets containing discrete characters (scale counts) because it does not require the data to fulfill the assumptions of linearity or unimodality (Marhold 2011;Paliy and Shankar 2016), thus allowing more flexible handling of mixed data sets.Average values were used for meristc data containing right and left counts.
Discriminant analysis of principal components (DAPC) from the ADEGENET package in R (Jombart 2021) was performed on the concatenated data set.The DAPC is a supervised analysis that places individuals from each predefined population into separate clusters (i.e., plots of points) bearing the smallest within-group variance that produce linear combinations of centroids having the greatest between-group variance (i.e., linear distance; Jombart et al. 2010).DAPC relies on scaled data from an 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 approximately 90% of the variation in the data set (Jombart and Collins 2015) as determined from a scree plot generated as part of the analysis.Retaining too many PCs forces false structure to appear in the data, while retaining too few runs the risk of missing true structure.
Morphospatial clustering and positioning among the species/populations were finally also analyzed using multiple factor analysis (MFA) on a concatenated data set comprised of 40 morphometric, 26 meristic, and six categorical color pattern characters (Table S3).The MFA was implemented using the mfa() command in the R package FactorMineR (Husson et al. 2017) and visualized using the Factoextra package (Kassambara and Mundt 2017).MFA is a global, unsupervised, multivariate analysis that incorporates qualitative and quantitative data (Pagès 2015), making it possible to analyze different data types simultaneously in a nearly total evidence environment.In an MFA, each individual is described by a different set of variables (i.e.characters), which are structured into different data groups in a global data frame-in this case, quantitative data (i.e.meristics and adjusted morphometrics) and categorical data (i.e.scale, tubercle, and caudal morphology).In the first phase of the analysis, separate multivariate analyses are carried out for each set of variables-principal component analyses (PCA) for the quantitative data sets and a multiple correspondence analysis (MCA) for categorical data.The data sets are then normalized separately by dividing all their elements by the square root of their first eigenvalues.For the second phase of the analysis, the normalized data sets are concatenated into a single matrix for a global PCA of the data.Standardizing the data in this manner prevents one data type from overleveraging another.In other words, the normalization of the data in the first phase prevents data types with the most number of characters or the greatest amount of variation from outweighing other data types in the second phase.This way, the contribution of each data type to the overall variation in the data set is scaled to define the morphospatial distance between individuals as well as calculating each data type's contribution to the overall variation in the analysis (Pagès 2015;Kassambara and Mundt 2017).
Non-parametric permutation multivariate analyses of variance (PERMANOVA) from the vegan package 2.5-3 in R (Oksanen et al. 2020) were used to determine if the centroid locations and group clustering of each species/ population in the multivariate analyses were statistically different from one another (Skalaski et al. 2018).Each analysis was based on the calculation of a Euclidean (dis) similarity matrix using 50,000 permutations of the input data, not on the output data of the multivariate analyses.A pairwise post hoc test calculates the differences between all combinations of population pairs, generating a Bonferroni-adjusted p value and a pseudo-F ratio (F statistic).A p < 0.05 is considered significant and larger F statistics indicate more pronounced group separation.A rejection of the null hypothesis (i.e., centroid positions and/or the spread of the data points (i.e., clusters) are no different from random) signifies a statistically significant difference between species/populations.

Assesssment of the conservation status
The conservation status of the newly described species of the alwisi group were assessed through the application of the IUCN Red List criteria, Extent of Occurrence (EOO) and Area of Occupancy (AOO).Extent of Occurrence and Area of Occupancy of each species was determined following the guidelines of the IUCN Red List Categories and Criteria (IUCN Standards and Petitions Committee 2022).

ND2 gene genealogy
Reconstructions using BI and ML methods resulted in highly similar topologies and branch lengths and hence only the Bayesian tree is shown (Fig. S1).Toplogy was identintical in both ML and BI analyses for the podihuna clade.However, there were slight differences in the topology in the kandiana clade but none of these differences were not strongly supported (Fig. S1).Our data confirm the results of both Agarwal et al. (2017) and Karunarathna et al. (2021) in the monophyly of podihuna (1.0/97; hereafter node support values are given for BI PP/BS respectively) and kandiana clades (1.0/99) (Fig. 1).The two novel species reported here were recovered with the Cnemaspis podihuna clade and hence only the C. podihuna clade is shown in Fig. 1 for clarity.Both new species were recovered in a strongly supported (1.0/100) monophyletic group comprising C. alwisi, C. gunasekarai, C. hitihamii, C. kohukumburai, C. nilgala, and C. punctata within the C. podihuna clade.The new species collected from Ethagala mountain, Ampara was monophyletic (1.0/100) and was recovered in a strongly supported clade (1.0/99) comprising C. nilgala and the undescribed species Cnemaspis sp. 9 (Fig. 1).The new species collected from Galgiriya mountain, Kurunegala was also monophyletic (1.0/100) and was recovered as the sister species (1.0/80) to C. alwisi, for which genetic data are presented here for the first time (Fig. 1).The newly included C. molligodai was monophyletic (1.0/100) and was recovered as sister to a clade (1.0/100) comprising C. podihuna and C. kandambyi.Cnemaspis kandambyi however, was recovered as nested inside C. podihuna (Fig. 1).

Sequence divergence
The average uncorrected pairwise genetic p-distance between new species from Ethagala mountain, Ampara and other taxa in the C. podihuna clade was 16.3% (range 11.7-20.6%)(Table 1) while it was 27.7% (range 20.8-32.1%) between the former and Indian Cnemaspis species (Table S6) exclusive of the members of the C. kandiana clade.The p-distance between new species from Ethagala Mountain, Ampara and other taxa in the C. kandiana clade was 28.3% (range 26.9-29.8%)(Fig. S1).The undescribed Cnemaspis sp. 9 has the lowest p-distance (11.7%) to the new species from Ethagala mountain, Ampara.The p-distance between the new species from Galgiriya mountain, Kurunegala and other taxa in the C. podihuna clade was 15.7% (range 8.9-19.5%)(Table 1) while it was 27.7% (range 20.1-32.1%) between the former and Indian Cnemaspis species, exclusive of the members of the C. kandiana clade.The p-distance between the new species from Galgiriya mountain, Kurunegala and other taxa in the C. kandiana clade was  29.5% (range 28.4-31.2%)(Fig. S1).Cnemaspis alwisi has the lowest p-distance (8.9%) to new species from Galgiriya mountain, Kurunegala.Species delimitation analyses was congruent among PTP and bPTP models, and confirmed all presently recognized species presented in the dataset, as well as the two newly described species (see the systematics section below).

Statistical ordination
The ANOVA and TukeyHSD post hoc tests recovered significant differences (p<0.05) between the new populations from Ethagala mountain and Galgiriya mountain forests and among them and all other species of Cnemaspis in the analyses across all meristic and normalized morphometric characters (Tables S2-S4).These analyses were corroborated by all three multivariate analyses that recovered notable morphospatial separation among several species.In the PCoA, the new population from Ethagala mountain was separated from all species except C. alwisi, C. gunasekarai, and C. punctata along PCoA axis 1 which accounted for 56.54% of the variation in the data set (Fig. 2A-D).It was separated from all species except C. gunasekarai, C. rajakarunai, C. gunawardanai, and C. rammalensis along PCoA axis 2 which accounts for an additional 10.96% of the data set.The new population from Galgiriya mountain was separated from all species except C. gunasekarai along axis 1, and the new population from Ethagala mountain, C. rajakarunai, C. gunawardanai, and C. rammalensis along axis 2. The PERMANOVA analysis of the PCoA data recovered statistically significant differences (p<0.05)among the centroid positions of the new taxa and all other species with a wide range of F.model values (Fig. 2B).These results were mirrored and corroborated in the DAPC, with the retention of the first two PCs accounting for 89.5% of the variation in the data set (Fig 2C).
The global MFA analysis recovered the greatest amount of separation among the species (Fig. 3A-C) with the first five dimensions accounting for 77.5% of the variation (Fig. 3C).Dimension 1 accounted for 41.6% of the variation and loaded most heavily for meristic data (Fig. 3B,C) followed by color pattern data for dimensions 2-4 and morphometric data for dimension 5 (Fig. 3C).The PERMANOVA analysis of the loadings of the first five dimension of the MFA data recovered statistically significant differences (p<0.05) in the position of nearly all species' centroids with a wide range of F.model values.However, the centroid position of the new population from Ethagala mountain was not significantly different from that of C. gunasekarai given their overlap along dimension 1.
Dorsal surface of the trunk with small, smooth, homogeneous granules, 126 paravertebral granules; 129 smooth midventral scales; 74 midbody scales across belly; 6/7  weakly developed tubercles on the flanks; ventrolateral scales small; granules on snout smooth and flat, larger than those on interorbital and occipital regions; canthus rostralis less pronounced, 12/13 smooth oval scales from eye to nostril; scales of the interorbital region oval and smooth; ear opening vertically oval, slanting from anterodorsal to posteroventral, 23/22 scales between anterior margin of the ear opening and the posterior margin of the eye.Supralabials 8/9 and infralabials 7/8, becoming smaller towards the gape.Rostral scale wider than long, partially divided (70%) by a median groove, in contact with first supralabial.Nostrils separated by 2/2 enlarged supranasals with 1 internasal and 1/1 postnasal; no enlarged scales behind the supranasals.Nostrils oval, dorsolaterally orientated, not in contact with first supralabials.
Color of living specimens.The dorsal color of the head, body and limbs generally from black-brown to ash with four or five paired black blotches; in between black blotches on the dorsum are small grey-white paravertebral blotches; tail is light-brown with 11-14 faded grey cross bands along the tail dorsum (Fig. 5); a small black spot is present in the occipital area; dorsal pattern on head reticulated; pupil circular and black with the surrounding scales yellowish white; limbs with black and grey band like dorsal patches; manus and pes with irregular, wide, black and white crossbands on the dorsum; supralabials and infralabials are dirty white and dusted with black; lateral body has a combination of 8-10 black and grey spots; chin, gular, and throat scales are bright yellow without dark spots in males; whole ventral side with dirty white without dark spots in females; pectoral region, arms, abdomen and cloacal scales dirty white; post-cloacal spur and subcaudal scales pale yellow, without spots.
Color of preserved specimens.Both males and females dorsally light brown, mixed with dirty white and dark paravertebral blotches; dark and grey cross bands on the dorsum of the tail; limbs with brown and dirty white band like dorsal patches; ventral surface dirty white without markings; subcaudals with dark brown margins; ventrolateral side of the trunk dusted with black.
Etymology.The specific epithet (jayaweerai) is a Latinized eponym in the masculine genitive singular, honouring Mr. Shanthasiri Jayaweera (a senior member and former president of the Young Zoologists' Association (YZA) of Sri Lanka; founder of the Junior YZA; senior instructor of the Fish study group of the YZA; renowned wildlife artist; a great educator; wildlife conservationist) for his friendship and valuable contribution to environmental conservation in Sri Lanka.
Distribution and natural history.The type locality, Ethagala forest (7.487114 to 7.514992N and 81.489325 to 81.525203E) is an isolated mountain with large pointed granite rock outcrops supporting a forest dominated by tall shade-bearing trees and bushes (e.g.Vatica obscura).
The dominant vegetation type is dry mixed semievergreen forests (Gunatileke and Gunatileke 1990), and the mean annual rainfall varies between 1000 and 1500 mm, received mainly during the northeast monsoon (November-February) (Somasekaram et al. 1988).The land area is approximately 450 ha and situated in the Ampara District, Eastern Province (dry bioclimatic zone) of Sri Lanka.The mean annual temperature of the area is 29.8-31.2°C,and its elevational range is 100-450 m a.s.l.Based on our preliminary investigations Cnemaspis jayaweerai sp.nov.appeared to be uncommon in Ethagala.Our survey of 20 ha revealed 16 (±0.3)geckos for 35 man-hours (field work done on 20 July 2019, 13-14 September 2019, and 18-19 March 2020).This species was restricted to granite caves (mainly those modified by natives 1,000 years ago) in closed canopy forest, and abandoned buildings associated with granite caves inside the forest (Fig. 6).Further, it prefers narrow (~6-10 mm), long (~120-390 mm) and deep (~80-260 mm) crevices as hiding sites.These microhabitats were poorly illuminated (light intensity: 0-621 Lux), relatively moist (relative hu-    Description of holotype.An adult male, 32.6 mm SVL, and 35.7 mm TAL.Body somewhat stout, relatively short (TRL/SVL ratio 39.0%).Head relatively small (HL/ SVL ratio 26.1% and HL/TRL ratio 66.9%), very narrow (HW/SVL ratio 14.8% and HW/HL ratio 56.7%), less depressed (HD/SVL ratio 9.5% and HD/HL ratio 36.6%) and distinct from neck.Snout relatively very long (ES/ HW ratio 74.0% and ES/HL ratio 41.9%), more than twice the eye diameter (ED/ES ratio 52.8%), more than half length of jaw (ES/JL ratio 62.2%), snout slightly concave in lateral view; eye relatively large (ED/HL ratio 22.1%), less than twice larger than the ear (EL/ED ratio 47.3%), pupil rounded; orbit length greater than eye to ear distance (OD/EE ratio 122.9%) and little longer than digit IV of the manus (OD/DLM IV ratio 103.2%); supraocular ridges not prominent; ear opening relatively large (EL/HL ratio 10.5%), deep, taller than wide, larger than nostrils; 2 rows of scales separate orbit from supralabials; interorbital distance is relatively broad (IO/ES ratio 99.7%), shorter than head length (IO/HL ratio 41.8%); eye to nostril distance greater than the eye to ear distance (EN/EE ratio 122.9%).
Dorsal surface of the trunk with small, smooth, homogeneous granules, 131 paravertebral granules; 122 smooth midventral scales; 88 midbody scales across belly; 6/6 weakly developed, large tubercles on the flanks; ventrolateral scales not enlarged; granules on snout smooth and flat, larger than those on interorbital and occipital regions; canthus rostralis not pronounced, 13/14 smooth oval scales from eye to nostril; scales of the interorbital region round and smooth; ear opening vertically oval, slanting from anterodorsal to posteroventral, 23/21 scales between anterior margin of the ear opening and the posterior margin of the eye.Supralabials 8/9 and infralabials 7/7, becoming smaller towards the gape.Rostral scale longer than wide, partially divided (90%) by a median groove, in contact with first supralabial.Nostrils separated by 2/2 enlarged supranasals with 3 internasal and 1/1 postnasal; no enlarged scales behind the supranasals.Nostrils oval, dorsolaterally orientated, not in contact with first supralabials.
Mental quadtriangular, longer than wide, posteriorly in contact with 4 enlarged postmentals (smaller than mental, and larger than chin scales); postmentals contact and bordered posteriorly by 6 smooth chin scales (larger than nostrils), contact with the 1 st infralabials; ventral scales larger than chin scales, and larger than nostrils.Smooth, rounded, juxtaposed granular scales on chin and gular region; pectoral and abdominal scales smooth, subimbricate towards precloacal region, abdominal scales larger than dorsals; 25 belly scales across venter; smooth, subimbricate scales around base of the tail; 8/7 femoral pores; 22 unpored interfemoral scales; 6/6 large posterior femoral scales.Original tail of holotype longer than snout-vent length (TAL/SVL ratio 109.6%),less spine-like tubercles along tail side, homogeneous scales on the dorsal aspect of the tail directed backwards, tail with 2 or 3 enlarged flattened obtuse scales forming whorls; hemipenial bulge greatly swollen (TBW 3.1 mm), a very small, round post-cloacal spur on each side; smooth subcaudals are arranged into a median series of clearly enlarged, subhexagonal scales.
Color of living specimens.The dorsal color of the head, body and limbs generally grey to brown, with 5-7 paired small black and white paravertebral blotches and 4-6 white crown-like markings along vertebra in both sexes; tail is cinnamon-brown on the dorsum with 12-14 faded cream-white and dark crossbands along its length (Fig. 9); a large blackish nuchal dash present on the head, reticulated pattern on the occipital area, snout yellowish; the pupil circular and black with the surrounding scales reddish brown in both sexes; supraciliaries grey; supralabials and infralabials yellowish and dusted with black; a straight, thin and dark postorbital stripe on either side running towards ear; chin and gular region pure white; throat, pectoral, abdomen, cloacal and subcaudal scales are bright-yellow without dark spots; limbs, pale-yellow ventrally, with white dorsal spots and dark irregular bands dorsally in both sexes; manus and pes with irregular black and white crossbands on dorsum.
Color of preserved specimens.Dorsum faded brown intermixed with pale crown-like markings; paravertebral blotches distinct light-brown and dirty white spots intermixed; reticulated pattern on the occipital clearly distinct in both sexes; limbs with dark and light dorsal spots; dorsal tail is pale-brown with dark cross bands; chin, gular, pectoral, cloacal and subcaudals dirty white with some scales on thigh, tail base and arms with irregular darkbrown margins in both sexes.
Etymology.The specific epithet (nanayakkarai) is a Latinized eponym in the masculine genitive singular, honouring Mr. Ananda Lal Nanayakkara (a senior member and a senior instructor of the reptile study group of the Young Zoologists' Association (YZA) of Sri Lanka; reputed lawyer; policy maker, researcher) for his friendship and valuable contribution to reptile conservation in Sri Lanka.
Distribution and natural history.The type locality, Galgiriya (8.0910 to 8.1514N and 80.6246 to 80.6896E), is an isolated mountain range with granite rock outcrop forests dominated by tall shade-bearing trees (e.g.Hopea brevipetiolaris), belonging to the dry mixed semievergreen forest type (Gunatileke and Gunatileke 1990).The land area is approximately 1000 ha and situated in the Kurunegala District, North Western Province (dry bioclimatic zone) of Sri Lanka.The mean annual rainfall varies between 1,000 and 1,500 mm, which is received mainly during the northeast monsoon (November-February), but a lesser amount during the southwest monsoon (May-September).The mean annual temperature of the area is 29.2-30.5°C,and its elevational range is 120-550 m a.s.l.Based on our preliminary investigations (22 December 2019, 10 February 2020, and 16 August 2020), Cnemaspis nanayakkarai sp.nov.appeared to be common in Galgiriya forest.Our survey of 12 ha revealed 32 (±0.2) geckos for 24 man-hours.This species was restricted to granite caves (mainly prehistoric granite caves modified 1,000 years ago by natives) in closed canopy forest, and old buildings associated with granite caves within the forest (Fig. 10).Furthermore, this species prefers narrow (~7-12 mm), long (~180-450 mm) and deep (~130-210 mm) crevices as hiding places.These microhabitats were poorly illuminated (light intensity: 246-692 Lux), relatively moist (relative humidity: 65-80% and canopy cover: 70-95%) and moderately cool (ambient temperature: 29.5-31.2°Cand substrate temperature: 28.2-29.4°C).The new species was sympatric with several other gecko species: Gehyra mutilata, Hemidactylus depressus, H. frenatus, H. parvimaculatus, and H. triedrus.Eggs, hatchlings and juveniles were not observed in the habitat, but gravid females were observed in year 2019 and 2020.
Conservation status.Application of the IUCN Red List criteria indicates that Cnemaspis nanayakkarai sp.nov. is Critically Endangered (CR) due to having an area of occupancy (AOO) <10 km 2 (six locations -single forest block, 0.27 km 2 in total assuming a 100 m radius around each georeferenced location) and an extent of occurrence (EOO) <100 km 2 (1.45 km 2 ) in the North Western Province [Applicable criteria B2-b (iii)] (Fig. 7) of Sri Lanka.

Discussion
Sri Lanka's terrestrial habitats are home to unique assemblages of floral and faunal communities with high endemism (Pethiyagoda and Sudasinghe 2021).It is also one of the smallest biodiversity hotspots in the world (IUCN-SL and MOE-SL 2012).With 244 known reptile species inhabiting the island, reptiles represent a major contributor to Sri Lankan vertebrate diversity.Nonetheless, they are also threatened, mainly due to loss of habitat, climate change, and pollution across all reptile distribution zones (MOMDE 2019; Dayananda et al. 2021).Amongst the diverse reptile fauna of the island, the diversity of geckos (Family Gekkonidae) is remarkable; reaching now 65 species, which accounts for ~27% of the overall reptilian richness.With the inclusion of these two new species described herein, the species richness of Cnemaspis rises to 42 (of which 20 species have been described from 1852-2007, and 22 species described from 2008-2023), making it the most species-rich reptile genus of Sri Lanka.The two new species described in this paper have not been included in any previous phylogenies of the genus (Agarwal et al. 2017;Karunarathna et al. 2021).
Taxonomy of the genus Cnemaspis recently started to receive special attention with several revisionary attempts providing a more of less detailed subgeneric system for the genus (Bauer et al. 2007;Manamendra-Arachchi et al. 2007;Wickramasinghe and Munindradasa 2007;Vidanapathirana et al. 2014;Wickramasinghe et al. 2016;Agarwal et al. 2017;Batuwita and Udugampala 2017;Batuwita et al. 2019;de Silva et al. 2019;Karunarathna et al. 2019aKarunarathna et al. , 2019bKarunarathna et al. , 2019cKarunarathna et al. , 2021;;Karunarathna and Ukuwela 2019;Amarasinghe and Karunarathna 2020;Amarasinghe et al. 2021a;Pal et al. 2021).Molecular phylogenies, based on mitochondrial genes (ND2, Cyt b) and nuclear genes have showed that Sri Lankan Cnemaspis belong to two unrelated lineages, namely the kandiana clade and the podihuna clade (Agarwal et al. 2017).The results of phylogenetic analyses comprising only mitochondrial ND2 gene in the current paper and previous studies (Karunarathna et al. 2019c) generally supported results of Agarwal et al. (2017) and Karunarathna et al. (2021).Other notable results of our phylogenetic analyses include phylogenetic placement of C. molligodai as sister species to C. podihuna (Fig. 1), and of C. gunasekarai as sister to C. alwisi and the new species from Galgiriya mountain, C. nanayakkarai sp.nov.. Also our results suggest the placement of the new species from Ethagala mountain, C. jayaweerai sp.nov.as a sister species to a lineage comprising C. nilgala and the undescribed species Cnemaspis sp. 9 (Fig. 1).
Three species (C.gunasekarai, C. kandambyi and C. molligodai) were described on the sole basis of morphological data, while "Cnemaspis alwisi", which appeared in previous works (Agarwal et al. 2017;Karunarathna et al. 2019cKarunarathna et al. , 2021) ) under the voucher number AMB7447 has recently been described as C. gunasekarai (Amarasinghe et al. 2021a).Our sequences of the specimens from the type locality of C. alwisi are notably distinct from AMB7447, so in the present paper we provide the molecular data for this species for the first time and, on the basis of their genetic distinctiveness, we thus confirm the species-level status of C. alwisi (p-distances > 8.98%), C. gunasekarai (p-distances > 9.52%) and C. molligodai (p-distances > 13.61%).Interestingly, our analysis recovered C. kandambyi to be nested within C. podihuna (Fig. 1) and this clade was strongly supported (1.0/100).Species delimitation analyses, however indicated that C. podihuna comprises four distinct 'species' and further validated the species level distinction of C. kandambyi.The genetic distinctiveness of these populations was shown to be greater than the proposed species-level divergence threshold of 3.7% for ND2 (Agarwal et al. 2017), and, as we show below, both taxa demonstrated significant morphological differences from all other congeners (Figs 2, 3).
Since the revisions of Agarwal et al. (2017) and Karunarathna et al. (2019c), the number of known Sri Lankan Cnemaspis has increased rapidly.The number of described species of the podihuna clade for which their phylogenetic placement was confirmed by molecular data arose from six to 14 (Agarwal et al. 2017;Karunarathna et al. 2019cKarunarathna et al. , 2021;;present paper).With several lineages with unclear status, the number of species in this clade is likely to increase from 20 in the coming years.Lineages and groupings for Sri Lankan Cnemaspis seem to be stable (Agarwal et al. 2017;Karunarathna et al. 2019c) and the current study further emphasises this fact.Currently being the most complete, our reconstruction of phylogenetic relationships between species of the podihuna clade is also quite strongly supported.These findings suggest that the mitochondrial ND2 gene to be extremely effective in species recognition and resolving interspecific relationships for this genus.
According to the morphological data and the molecular phylogeny, here we identify four groups (generally shown in phylogenetic tree, see Fig. 1) within the podihuna clade: (i) the alwisi group consisting of eight species (C. alwisi, C. gunasekarai, C. hitihamii, C. kohukumburai, C. nilgala, C. punctata, C. jayaweerai  The two new species were recovered within the alwisi group of the podihuna clade (Fig. 1) indicating speciation in the isolated mountains within the dry and intermediate bioclimatic zone of Sri Lanka.Most members of alwisi group have a cryptic morphology and coloration, both of which serve as camouflage.Our findings suggest that these geckos prefer narrow (~6-12 mm), long (~120-450 mm), and deep (~80-260 mm) crevices as microhabitats and refugia.Interestingly, the ventral and dorsal coloration of their tails is unique to each species and may serve as intraspecific communication.These crepuscular geckos are mostly rupicolous, found in granite caves, and are occasionally found in old buildings (cement walls or clay bricks walls) associated with granite caves.The members of the alwisi group are widely distributed in the dry and intermediate bioclimatic zones (~200-500 m, a.s.l) with the exception of C. kohukumburai, which is recorded from the wet bioclimatic zone which is quite unusual.Remarkably, both the rammalensis and scalpensis groups are restricted to the wet zone, and the podihuna group is found in both the dry and intermediate zones.It is possible and our unpublished data also suggest that the rupicolous species (alwisi group) in the podihuna clade may have higher speciation rates compared to the arboreal species (podihuna group) of this clade, probably due to the isolation in the clustered rocky habitats scattered in the intermediate and dry zones, as well as in the wet bioclimatic zone.
According to the estimations by Agarwal et al. (2021), the podihuna clade diverged from other South Asian Cnemaspis in the mid-Eocene (~40-50 mya).The clade then went through two independent dispersal events into Sri Lanka and two recent independent dispersals to SE Asia (Agarwal et al. 2020;Pal et al. 2021).The time of divergence in the alwisi group was estimated at ~12 mya in mid-Miocene (Agarwal et al. 2017(Agarwal et al. , 2020)), and the in-situ speciation was probably triggered by the loss of rainforest belt between India and Sri Lanka at that time (Morley 2018;Pethiyagoda and Sudasinghe 2021).Our research shows the distribution and divergence of the alwisi group has taken place along the 500 meter contour in Sri Lanka, and further emphasize the importance of isolated misty-mountains (Fig. 7).Recent studies (Agarwal et al. 2020;Pal et al. 2021) indicate that South Asian species of Cnemaspis originated in the rainforests of Western Ghats during the Eocene-Paleocene (~50-60 mya), and diversification is believed to have happened in the late Paleocene (~30-40 mya) and Miocene (23-25 mya) during a period of global warming.At that time, lowland rainforests of Malaysia, Philippines, India, Seychelles, Africa, Indonesia, and Sri Lanka were dominated by tree species of the family Dipterocarpaceae (Ashton and Gunatilleke 1987;Rust et al. 2010).
A few species of Dipterocarpaceae (e.g.Hopea sp. and Vatica sp.) that usually dominate the rainforest flora of Sri Lanka are found in several small isolated mountains (~400 m asl.) in the dry and intermediate bioclimat-ic zones of Sri Lanka.Soil erosion and weathering are believed to have contributed greatly to the isolation of these dry and intermediate zone mountains from central highlands and shaping these forest types with their unique biodiversity (Karunarathna et al. 2021;Pethiyagoda and Sudasinghe 2021) in this bioclimatic region.These isolated mountains are usually covered with mist after 1700 hrs to early morning (ca.0800 hrs) most of the year, and therefore, the forests in these mountains maintain wet microclimates which are well suited for moisture dependant herpetofauna as well as other faunal groups (Amarasinghe et al. 2021a; Karunarathna et al. 2019c).However, these moist microclimates are not found on isolated mountains that do not reach 300 m in elevation in the dry and intermediate bioclimatic zones, and most of them have dry conditions with low faunal diversity.
Almost all species of the alwisi group are found within relatively cool, moist habitats (substrate temperature: 25-29°C; relative humidity: 70-90%), which are well shaded (canopy cover: 60-90%; light intensity: 380-850 Lux) and have high-profile mature trees.The species in the alwisi group are point-endemic microhabitat specialists where distribution ranges are limited to <20 km 2 .Our studies suggest the possibility of speciation of alwisi group (currently comprising eight species) found in geographically isolated mountains with granite caves, rock outcrops and favourable environmental conditions; we therefore predict that the number of species will likely increase to more than 16 species in total.The high species richness of Cnemaspis in Sri Lanka may be accounted by the possibility of multiple, independent colonization events from the Indian mainland with subsequent, geographically-isolated in-situ speciation (Agarwal et al. 2020;Pal et al. 2021;Pethiyagoda and Sudasinghe 2021).It is also interesting to note that species like C. alwisi have genetically distinct populations on either sides of the mountains on which they occur (Fig. 1).Such findings need further investigation and may suggest the low dispersability of these small geckos, which prefer specific microhabitats.
Ethagala and Galgiriya are isolated misty-mountains and rock outcrops embedded within forest habitats having granite caves incorporated with historical Buddhist monasteries.These forest habitats with Buddhist monasteries serve as refugia for reptiles and other faunal groups, and it is imperative to conserve these habitats.Consequently, out of the 42 Cnemaspis species occurring in Sri Lanka, 35 (~84%) species are found inside Buddhist monasteries with granite caves (Manamendra-Arachchi et al. 2007;Wickramasinghe et al. 2016;Karunarathna et al. 2019c, Amarasinghe et al. 2021a).It is very likely that future studies on the biogeography of alwisi group in Sri Lanka will highlight the importance of these isolated habitats in generating and maintaining the diversity of these unique groups of geckos on the island.At the same time, it is important to note that the point-endemic and range-restricted species described herein, which are highly sensitive to changes in the habitat, would be severely affected by habitat degradation.Hence, past and present studies have emphasized the importance of conserving such iso-lated habitats throughout the country (Amarasinghe et al. 2021b).All these isolated habitats are susceptible to human-induced habitat degradation, including clear cutting and timber felling, forest fragmentation, granite mining, rubber cultivation, irresponsible vegetable farming, invasive species, human settlements, road and other infrastructure development, and waste disposal.It is evident that conservation of isolated misty mountains is crucial in conservation of the unique fauna and flora of the island.It is also imperative to inform the general public, educators, natural resource and land managers, and private entrepreneurs about environmental legislation through printed and electronic media and outreach activities (Karunarathna et al. 2017a(Karunarathna et al. , 2017b)).

Figure 1 .
Figure 1.Bayesian inference tree of the podihuna clade derived from the analysis of 1347 bp of the ND2 gene sequences representing 142 taxa (see Fig. S1 for the complete tree).Inset: Cnemaspis nanayakkarai sp.nov.from Galgiriya mountain.The black circles at nodes correspond to BI PP and ML BS support values greater than 0.95 and 70 respectively.For voucher specimen information and GenBank accession numbers see Table S1.The grey vertical bars indicate results of molecular species delimitation analyses (bPTP and PTP).
Figure 1.Bayesian inference tree of the podihuna clade derived from the analysis of 1347 bp of the ND2 gene sequences representing 142 taxa (see Fig. S1 for the complete tree).Inset: Cnemaspis nanayakkarai sp.nov.from Galgiriya mountain.The black circles at nodes correspond to BI PP and ML BS support values greater than 0.95 and 70 respectively.For voucher specimen information and GenBank accession numbers see Table S1.The grey vertical bars indicate results of molecular species delimitation analyses (bPTP and PTP).

Figure 2 .
Figure 2. A Principal coordinate analysis (PCoA) of the meristic and normalized morphometric characters showing the morphospatial relationships of Cnemaspis jayaweerai sp.nov.and C. nanayakkarai sp.nov. to the other species in an unsupervised analysis.B Distribution of the F.model values between all possible pairs of species analyzed.The higher the F.model value the more significantly different are the species of each species pair.C Discriminant analysis of principal components (DAPC) of the meristic and normalized morphometric characters showing the morphospatial relationships of Cnemaspis jayaweerai sp.nov.and C. nanayakkarai sp.nov. to the other species in a supervised analysis.D Color coding matches that of Eigenvalue scree plot of the PCoA showing the contribution of each axis to the variation in the data set.

Figure 3 .
Figure 3.A Global multiple factor analysis (MFA) of the color pattern, meristic, and normalized morphometric characters showing the complete separation in the morphospatial relationships of Cnemaspis jayaweerai sp.nov.and C. nanayakkarai sp.nov.and all other species in an unsupervised analysis.B Bar graph showing the percent contribution of each meristic and normalized morphometric character to the variation in the data set.C The percent contribution of each data type to the first five dimensions of the MFA.

Figure 4 .
Figure 4. Cnemaspis jayaweerai sp.nov.male holotype (NMSL.2019.23.01):A dorsal aspect of head, B lateral aspect of head, C ventral aspect of head, D homogeneous scales on dorsal surface of trunk, E lateral surface of trunk, F smooth ventral scales, G cloacal characters with only femoral pores, H subdigital lamellae of manus, I subdigital lamellae of pes, J dorsal side of tail, K lateral side of tail, L large subhexagonal subcaudals.Scale bar: 2 mm.Photos: Suranjan Karunarathna.

Figure 6 .
Figure 6.General habitat of Cnemaspis jayaweerai sp.nov. in Ethagala, Ampara District of Sri Lanka: A an isolated granite mountain with dense forest and canopy cover, B a granite cave habitat with lots of cracks and crevices, C historic granite caves with ancient inscriptions in a Buddhist monastery in Ethagala.Photos: Suranjan Karunarathna.

Figure 7 .
Figure 7. Currently known distribution of the members of the alwisi and rammalensis groups with new species -Cnemaspis jayaweerai sp.nov.and Cnemaspis nanayakkarai sp.nov.with the three major mountain habitats and seven elevational ranges in Sri Lanka (abbreviations: AZ Semi-arid zone, DZ dry zone, IZ intermediate zone, and WZ wet zone; Inset top to bottom: Knuckles mountain massif, Central highlands, and Rakwana mountain massif).Photos: Vimukthi Weeratunga and Ashan Geganage.Map: Buddhika Madurapperuma.

Figure 8 .
Figure 8. Cnemaspis nanayakkarai sp.nov.male holotype (NMSL.2021.06.01):A dorsal aspect of head, B lateral aspect of head, C ventral aspect of head, D homogeneous scales on dorsal surface of trunk, E lateral surface of trunk, F smooth ventral scales, G cloacal characters with only femoral pores, H subdigital lamellae of manus, I subdigital lamellae of pes, J dorsal side of tail, K lateral side of tail, L large subhexagonal subcaudals.Scale bar: 2 mm.Photos: Suranjan Karunarathna.

Figure 10 .
Figure 10.General habitat of Cnemaspis nanayakkarai sp.nov. in Galgiriya, Kurunegala District of Sri Lanka: A an isolated granite mountain ridge with a dense forest, B historic granite cave complex with lots of cracks and holes in a Buddhist monastry, C dense canopy cover with humid conditions inside the forest.Photos: Suranjan Karunarathna.
sp. nov.and C. nanayakkarai sp.nov.), (ii) rammalensis group consisting of three species (C.gunawardanai, C. rajakarunai, C. rammalensis), (iii) the podihuna group consisting of four species (C.kandambyi, C. manoae, C. molligodai and C. podihuna), and (iv) the scalpensis group consisting of five species (C.anslemi, C. gemunu, C. godagedarai, C phillipsi and C. scalpensis).Cnemaspis molligodai, C. alwisi and C. gunasekarai which were included in the phylogenetic analyses for the first time in the present study were recovered within the podihuna clade.This was expected because of their morphological similarity to other members of the podihuna clade characterised by the presence of enlarged hexagonal/subhexagonal subcaudal scales, and homogeneous dorsal scales.It is interesting to note the recovery of C. kandambyi within C. podihuna which rendered the latter paraphyletic.The confirmation from the species delimitation analyses of the presence of four genetically distinct species within C. podihuna including C. kandambyi may indicate the presence of cryptic species.However, additional samples are necessary to examine the validity of these putative taxa, given that C. podihuna is a widespread species that may comprise several genetically distinct lineages.

Table 1 .
Uncorrected pairwise genetic distances in the ND2 gene between Cnemaspis jayaweerai sp.nov.and Cnemaspis nanayakkarai sp.nov.and the members of the podihuna clade in which these two new species belong.

Table 2 .
Comparison of morphological and morphometric characters of C. jayaweerai sp.nov.(from Ethagala mountain) and C. nanayakkarai sp.nov.(from Galgiriya mountain) with the other congeners of the podihuna clade in Sri Lanka that can be used to diagnose two new species from the C. alwisi group and C. rammalensis group (Abbreviations: HOMG -homogenous, SMOO -smooth, SUBX -subhexagonal, mm -millimetres).