Phylogeny and biogeography of Sumatra´s cloud forest lizards of the genus Dendragama and status of Acanthosaura schneideri

Abstract Lizards of the genus Dendragama are endemic to the highland cloud forests of Sumatra’s Barisan Mountain Range in western Indonesia, and recent studies have uncovered widespread diversity within the genus. Here, a suite of morphological characters and mitochondrial DNA are used to compare three geographically isolated populations of D. boulengeri from (1) Mount Kerinci in Jambi province, (2) Mount Marapi of west Sumatra, and (3) the Karo Highlands of north Sumatra. Additional phylogeographic analyses with two recently described sister species, D. australis and D. dioidema were conducted. Five genetically distinct clades of Dendragama, all distributed allopatrically of one another were identified and some are suspected to inhabit small distributions. Morphological and genetic data confirm the Karo Highlands population D. schneideri (previously Acanthosaura schneideri Ahl, 1926) should be revalidated from the synonymy of D. boulengeri. Dendragama schneideri is endemic to montane forests of the Karo Highlands surrounding Lake Toba in Sumatra Utara province. Pairwise genetic distances of 6–11% separate D. schneideri from congeners. Two distinct clades of D. boulengeri from Mount Kerinci and Mount Marapi were identified, which are 5.0% genetically distant from one another. Using morphological characters, we provide the first key for distinguishing between species of Dendragama. Based on biogeographic patterns and levels of genetic variation it is suspected that at least 18 other isolated cloud forest locations may hold new species or divergent populations of Dendragama but lack survey work. Collectively, these comparisons among populations of montane lizards further elucidate the complex biogeographic history of Sumatra’s montane forest species and the first phylogeny of the genus Dendragama.


Introduction
Uncovering tropical diversity is essential for conservation initiatives and understanding complex ecological and evolutionary processes (Ladle and Whittaker 2011). However, many regions and taxonomic groups across the globe remain largely unstudied, and for a multitude of reasons, Indonesia is particularly under-represented in this regard (O'Connell et al. 2018). Sumatra's Barisan Mountains are especially interesting because the entirety of the range is a volcanically active strip of the "Pacific Ring of Fire" which runs along the edge of Sumatra's west coast.
Volcanic activity and other historical biogeographic pressures throughout the Barisan Range have led to the development of a fascinating array of biodiversity (Demos et al. 2016), including a whole suite of montane forest species that live in high elevation cloud forest (temperate moist forest typically between 1,300 and 2,800 m of elevation in the tropics). However, Sumatra's cloud forest biodiversity has remained largely unexplored and the pressures that have driven the development of that diversity are not yet fully understood (O'Connell et al. 2018). Reptilian diversity is no exception in this regard despite the fact that Sumatra is known to harbor some of the most evolutionarily unique lineages of reptiles, including enigmatic species like the Modigliani's nose-horned lizard which was recently rediscovered (Putra et al. 2020).
Interestingly, despite the lack of herpetofaunal survey work throughout the region, Sumatra is already considered the most draconine (family Agamidae) diverse island in southeast Asia (Harvey et al. 2014). Given the complex geologic history, the presence of several isolated montane forest islands, and the wealth of diversity that has been discovered in nearby regions, it is clear that montane agamid lizards remain under studied and are an excellent group for better understanding the biogeographic history of the region (Harvey et al. 2014).
There has been considerable uncertainty regarding the taxonomic status of Dendragama boulengeri, which was considered to belong to a monotypic genus until only recently (Harvey et al. 2017). Early studies reported D. boulengeri from isolated montane forests of Jambi, west Sumatra (Sumatra Barat) and north Sumatra (Sumatra Utara) provinces (Doria 1874). However, Harvey et al. (2017) recently described two new species, D. australis from south Sumatra and D. dioidema from Aceh Province. They also provide a thorough redescription of D. boulengeri. In their paper, D. boulengeri is redescribed as a species distributed throughout much of Central Sumatra, including the type locality (Mount Singgalang in Sumatra Barat) and nearby Mount Marapi (paralectotype locality).
Previously, a population of Dendragama was described as Acanthosaura schneideri (Ahl 1926) from Sumatra Utara Province. The species was mentioned infrequently thereafter. In an unpublished dissertation, Moody (1980) conducted a thorough family-wide review of the Agamidae and split the genus Calotes into four genera: Bronchocela, Calotes, Dendragama, and Pseudocalotes. However, there is no mention of A. schneideri in his work. Only later did Manthey and Grossmann (1997) transfer A. schneideri to the synonymy of D. boulengeri. To date, there has been no phylogenetic analyses of Dendragama to confirm the status of the population from Sumatra Utara in relation to other populations.
We examined a series of Dendragama collected by MBH and ENS from their herpetofaunal inventory conducted between 2012 and 2014 throughout the Barisan Mountain Range of Sumatra. These samples include specimens from various isolated mountain peaks across much of central, northern, and southern Sumatra. Using an integrative analysis, we investigate species boundaries among these populations and biogeographic patterns revealed by the phylogenetic analyses.

Biological inventory
A thorough herpetofaunal survey was conducted across the highland forests of Sumatra's Barisan Mountain Range between 2012 and 2014. An international team of collaborators systematically targeted mountains based on geographic isolation from one another. We predominantly collected specimens at night. In total, we collected 195 specimens of Dendragama from localities throughout the highlands of Sumatra between 1231-2253 m elevation. We recorded GPS coordinates and ecological data on site. Animals were euthanized following appropriate IACUC protocols and DNA samples were taken for future identification prior to preservation in 10% formalin. Photographs were taken before and after euthanizing.

Counts and measurements
We scored 32 different morphological characters for each specimen from the three populations of D. boulengeri. To avoid systematic errors introduced by separate observers, K. Shaney collected all mensural and meristic characters. Sex was determined by examining the gonads. We examined 19 D. boulengeri collected from Mount Marapi, west Sumatra (Marapi population), nine specimens from Mount Kerinci and Mount Tujuh, Jambi (Kerinci population), and 15 specimens from various mountains across the Karo Highlands of Sumatra Utara, collected near the type locality of A. schneideri (Karo population; Fig. 1).
Because measurements and scores are often done differently depending upon the study, some of the morphological characters used here require further comment (Harvey et al. 2014). We measured eye-nostril distance from the anterior ocular angle (bony edge of the orbit) to the tip of the snout. We consider the last supralabial to be lower and more elongate than the supralabials in front of it, and the last infralabial is positioned directly below the center of the last supralabial. We counted gular scales beginning immediately behind the mental (or the first pair of infralabials when in medial contact) until the transverse crease where the gulars meet the pectoral region. From the crease, we counted ventral scales to the last scale on the anterior edge of the cloacal opening. We counted nasal-rostral scales as the number of scales between the nasal and rostral. Postrostral scales were counted as the number of scales contacting the rostral which were not supralabials. We counted canthals as the number of scales between the nasal scale and the first supraocular scale. The number of scales between the supralabials and the first canthal touching the nasal scale are described here as the "scales between first canthal and supralabial". Circumorbital scales include a canthal and the postciliary modified scale. We counted postmentals (or chin shields) if they contacted the infralabials. The number of midbody scales around the body included the ventrals and a scale of the dorsal crest. We counted lamellae on the fourth digit of the hands and feet, starting from the interdigital skin at the base of the digit and extending to the claw (i.e., including the elongate ungual scale). We counted nuchal crest scales from the first projecting scale to the last enlarged scale before the pectoral gap. Only projecting scales were counted and small flat vertebrals and paravertebrals that interrupt the crest were excluded.
To the nearest 0.1 mm with digital calipers, we measured snout-vent length (SVL) from the tip of the snout to the anterior lip of the vent and tail length by straightening the tail along the edge of a ruler and measuring from the posterior edge of the vent to the tip of the tail. We measured head length from skin covering the posterior edge of the mandible to the tip of the snout, trunk length as the distance from the axilla to the groin, and hand and foot length from the proximal margin of the palm or sole to the tip of the claw on the fourth digit. We measured the brachium as the length of the entire humerus and the antibrachium from skin covering the proximal end of the ulna (antebracheal fold) to the base of the palm. We measured the shank as the length of the tibia. The proximal and distal ends were determined with the elbow or leg flexed 90°. We measured internarial distance as the distance between the upper edge of each nostril, bony orbit, and tympanum width from the anterior to the posterior edge of each. Additional specimens examined from sister taxa and outgroups are provided in Suppl. material 1.

Statistical analyses
Using our mensural and meristic data we compared the Karo, Kerinci, and Marapi populations of Dendragama boulengeri. For meristic characters, we compared means between the three populations using Tukey's test after confirming assumptions of normality (using the Shapiro-Wilk test) and homoscedasticity (using Levene's test). Dendragama boulengeri populations from the Karo population represent the group that we hypothesize to be a distinct species, D. schneideri; however, we will continue to refer to this as the "Karo population" until the results section. Thus, for clarity between the text, figures, and tables, when we refer to D. boulengeri (Karo population) or D. schneideri we are referring to the same population.
When making comparisons among populations, we analyzed males only for head width and head length, because we found these traits to be sexually dimorphic in a preliminary study of our series from Marapi. To investigate sexual dimorphism and to compare mensural characters among populations, we used analysis of covariance treating SVL as a covariate. To avoid inflation of the type I error rate in our morphometric comparisons, we performed three additional calculations. First, we made Bonferroni corrections to the probability scores for the tests among populations (Glover and Mitchell 2016). Second, having identified several apparent morphometric differences in among populations, we then verified the difference by rerunning the analysis using a different measurement as covariate in each apparently different trait: eye-nostril distance as a covariate for comparisons among thigh lengths, length of brachium for comparisons among hand length and length of shank for comparisons among foot length. Third, as a final validation of these results, M. B. Harvey measured SVL and tail length of a separate sample of nine Dendragama housed in the MZB and compared them to his own measurements for specimens from the type locality and Marapi.

DNA extraction and amplification
We digested tissue in 100 μL of lysis buffer, then added 5 μL of proteinase K (20 mg/ ml) and incubated at 55 °C for 1-6 hours. After incubation, we added 1.8 μL of serapure beads (Rohland and Reich 2012) for every 1 μL of digested sample. DNA extraction was carried out following the same methods that are used in PCR cleaning protocols described in AMPure magnetic beads literature (Agencourt, Bioscience, Beverly, MA, USA).
Phylogenetic Analyses -We extracted genomic DNA from 17 specimens of Sumatran Dendragama, including: D. boulengeri from Mount Marapi (Marapi population), D. boulengeri from Mount Kerinci (Kerinci population), D. boulengeri representing the Karo Highlands population (Mount Sibuatan, Mount Pangururan and Vicinity of Tele), D. dioidema (Mount Kaba, Mount Patah and Mount Dempo) and D. australis (Berni Terlong and Takengan). It's important to note that individuals from the Karo Highlands were not collected precisely from the type locality of A. schneideri but were found in nearby geographic locations. We then combined new sequences from these specimens with sequences already published by Shaney et al. (2016) and Harvey et al. (2014Harvey et al. ( , 2017 on GenBank. The published sequences include the outgroup taxa Bronchocela cristatella (Kuhl, 1820), Lophocalotes ludekingi (Bleeker, 1860), and Pseudocalotes tympanistriga (Gray 1831) see Suppl. material 2.
ND4 provided sufficient data for us to generate a phylogenetic hypothesis of Dendragama. We sequenced a fragment of the NADH dehydrogenase subunit 4 (ND4) gene using the forward primer "ND4" (CACCTATGACTACCAAAAGCTCATGTA-GAAGC) and reverse primer "LEU" (CATTACTTTTACTTGGATTTGCACCA). The ND4 thermal cycle profile consisted of an initial denaturation at 94 °C for three minutes, followed by 30 cycles of denaturation at 94 °C for 30 seconds, a 50 °C annealing phase for 45 seconds and a 72 °C extension for one minute, followed by a 72 °C extension for seven minutes, then a holding phase at 4 °C.
All sequences were aligned using the Geneious aligner implemented within Geneious v. 6.1.8 (Kearse et al. 2012). ND4 sequences range in length from 616 to 934 bp. We identified the most likely model of evolution for each codon position using Bayesian information criteria implemented in PartitionFinder (Lanfear et al. 2012). We partitioned codon positions using GTR+ Γ. We conducted maximum likelihood analyses using raxmlGUI (Sylvestro et al. 2012). We utilized the thorough bootstrapping setting, sampling over 10 runs of 10,000 repetitions. We carried out Bayesian phylogenetic analysis using MrBayes v3.2.1 (Ronquist and Huelsenbeck 2003). We used four independent runs (nruns = 4) and four chains (three heated chains and one cold chain) for 10 million generations, sampling every 100 generations. We discarded the first 25% of samples as burn-in (Rambaut et al. 2014). We confirmed adequate mixing and assessed the appropriate amount of burn-in and convergence by inspecting the trace files in the program TRACER v1.6 (Rambaut et al. 2014). We conducted UPGMA analyses and calculated uncorrected pairwise distances using Mega 5.1 (Tamura et al. 2011).

Phylogenetics and biogeography
Both our Maximum Likelihood and Bayesian Analyses revealed the same relationships within Dendragama, the only difference being that the ML tree returned slightly lower support values (Fig. 2). Both analyses found Lophocalotes to be sister to Dendragama, followed by insular Pseudocalotes. We refer to Pseudocalotes from Java and Sumatra as "insular" because they are not closely related to Pseudocalotes from mainland Asia (Harvey et al. 2014).
Within Dendragama five clades have strong nodal support (posterior probabilities between 0.98 and 1.00): D. australis (Mount Dempo, Kaba and Patah region), D. boulengeri (Mount Kerinci), D. boulengeri (Mount Marapi), D. dioidema (Aceh Highlands) and D. boulengeri (Karo Highlands). For clarity, the Karo Highlands populations will be discussed in the redescription section of D. schneideri. Each of these clades is restricted to geographically isolated locations. Dispersal between populations is highly unlikely to occur given that all Dendragama rely on cloud forest habitat which is not continuous between any of these populations.
Dendragama schneideri from the Karo highlands represent a genetically distinct group from central Sumatra and includes all specimens from Mount Sibuatan, Mount Pangururan and the Vicinity of Tele (Fig. 2). Dendragama boulengeri populations from Mount Marapi (paralectotype locality of D. boulengeri) and D. boulengeri populations from Mount Kerinci are geographically isolated by elevation and breaks in montane cloud forest patches. Based on biogeographic patterns, there are 18 distinct localities that might hold new species of Dendragama.
2-4, 2 (5%), 3 (90%), 4 (5%) 2-3, 2 (17%), 3 (83%) 2-3, 2(37%), 9 (63%) Circumorbital Scales 13-16, 13 (37%), 14 (53%), 15 (5%), 16 (5%) Dendrama boulengeri specimens from Kerinci have relatively shorter tails, hands, and feet than specimens from the other two D. boulengeri populations. They also have shorter thighs than specimens from Karo and a smaller orbit than specimens from Marapi. Small specimens from Karo have a relatively smaller orbit than small specimens from Marapi; however, orbits are about the same size for larger specimens from the two populations. Our limited data suggests a different growth trajectory for the orbit at Karo vis-à-vis Marapi, but having violated the assumption of parallel regression lines, we do not report a probability for this comparison between Karo and Marapi. As detailed in  the methods, we confirmed each of these morphometric differences by treating other measurements as covariates. Moreover, a separate sample of nine Dendragama from Kerinci had relatively shorter tails (F 1,27 = 7.75, P = 0.010) than the sample of D. boulengeri from both Marapi and the type locality described by Harvey et al. (2017). We did not find differences among populations for eye-nostril distance, pectoral width, length of body, length of shank, length of brachium, length of antebrachium, or internarial distance (P > 0.12). Tables 2 and 3 provide statistics for the three populations compared.
In addition, Dendragama from the Marapi and Kerinci populations have a bright yellow buccal epithelium and tongue, whereas lizards from the Karo population have a pink to red buccal epithelium and tongue. Along their lower flanks, lizards from Karo have numerous distinctly enlarged tubercular scales. In contrast, specimens from the other populations lack these scales.

Species delineation
Our analysis revealed numerous differences between the Karo population on the one hand and the Marapi and Kerinci populations on the other. Numerous different means and high genetic divergence is evidence of an interruption in gene flow among these populations, but is of only limited diagnostic value. However, we also identified four fixed characters that distinguish the Karo population from the other two. Unlike these populations (characters in parentheses), the Karo specimens have pink to red buccal epithelia (yellow , Fig 4), numerous enlarged tubercles on the lower flank (scales of lower flanks homogenous or with few slightly enlarged scales), 59-68 scales around midbody (77-89), and 13-19 dorsals pointing upward and backward at midbody (20-25). Restricted to highland areas above 1,200 m, inhospitable lowlands isolate the Karo population from all other populations and species of Dendragama. Direct comparison of Ahl' (1926) type of Acanthosaura schneideri to the Karo specimens reveals that they are the same species. Accordingly, our analyses support the removal of A. schneideri Ahl, 1926 from the synonymy of Dendragama boulengeri Doria, 1888. Hereafter the Karo population should be recognized as a distinct species Dendragama schneideri.
Holotype. An adult male (ZMB 15664, Fig 5)  Diagnosis. A species reaching at least 201 mm in total length (SVL) and distinguished from congeners by the following characters: (1) midbody scales 58-67; (2) dorsal scales heterogeneous across flanks (Fig 6); (3) strongly keeled white/yellow scales randomly distributed along flanks (more numerous and distinct in females); (4) ventral scales 48-59; (5) banding pattern along flanks often muddled, but typically vertical when present; (6) mouth and tongue pink to red in life; (7) narrow, vertical black stripes across dorsal crest, limbs, digits and most of tail; (8) female color in life dark brown, yellow and black with amber coloration on underside, while males green and lacking amber coloration along ventral surface; (9) dorsal and nuchal crest clearly separated by 5-9 dorsal scales; (10) dorsal crest serrate, extending to base of the tail, comprised of 23-31 projecting, triangular scales; (11) a series of 3-4 enlarged tubercles present along the chin of males and females; (12) A series of 12-18 strongly keeled, white/yellow femoral spines present (combined count on both sides).
Description and variation. The description is based on the 19 referred specimens. Where appropriate we provide character state frequencies or means ± standard deviation in parentheses. When available and not subject to interobserver biases, we also provide data gathered by MBH for the holotype in brackets.
Nuchal crest clearly separated from dorsal crest and gap between crests spanning 5-9 scales; dorsal crest serrate, continuous down to tail; scales on dorsum, large and heterogeneous, with series of enlarged strongly keeled, yellow/white scales in row below dorsal crest; all other scales along dorsum and flank smooth to feebly keeled; scales along flank consistent with dorsum, with more enlarged strongly keeled scales in vertical rows along sides; midbody scales 58-67 (61.36 ± 2.8) [61], gulars smooth 32-44 (36.81 ± 3.51) [30]; ventral scales 48-59 (52.45 ± 3.14) [52], ventrals keeled from chest to lower abdomen before transitioning to smooth scales near precloacal area; precloacal scale width small 0.75-1.4 (1.02 ± 0.22); scales along limbs strongly keeled, with continuation of keeled scales down to fingers on both hands and feet; subdigital lamellae on finger IV 22-26 (24.27 ± 1.19) Figure 6. A flank of male Dendragama boulengeri (MZB 9825) and its thin, horizontal banding patterns, thick vertical bands along dorsal crest, small homogenous scales, and lack of enlarged, keeled scales B flank of male D. schneideri (UTA 62868) and its lack of horizontal banding along the flanks, thin vertical bands along dorsal crest, large heterogeneous scales, and enlarged, strongly keeled scales dispersed across the flanks (photographs by ENS).

Coloration in life.
There is distinct sexual dichromatism in this species and coloration changes in all Dendragama in response to rough handling. Females of Dendragama schneideri are typically shades of dark brown, green, and black with vertical black and yellow bands running along the extent of the dorsal crest. Bands extend almost to the end of the tail; tail bands 14-18 (12.2 ± 0.83), and enlarged green, yellow or white, strongly keeled scales are present intermittently along the flanks. Black and yellow/ green bands also extend along all limbs, hands, and feet. A black spot is present under the base of the nuchal crest as in other species of Dendragama. The throat has amber and brown coloration, which may or may not be broken up by small lateral brown lines. Brown and amber coloration extends along the lower flanks and all the way to the end of the tail. Yellow and black lines radiate around the eyes and across much of the face. Yellow, green or white enlarged tubercles are present below the eye and ear, and the mouth is pink to red.
Males may also be brown but are typically much lighter in coloration. They are often bright green and yellow with incomplete stripes of black scales, which zigzag vertically along the flanks. Black bands extend along the length of the dorsal crest and throughout the extent of the tail. Bands also cross the arms, legs, hands, and feet. A black prescapular blotch is present under the base of the nuchal crest, but may be less pronounced in some specimens. The venter is much lighter than in females, with a white or cream gular region, with some brown shading along the ventral side. Darker individuals may have some brown shading along the gular region as well. Green or yellow and black stripes radiate out from the eyes and, as with females, the mouth is pink to red.
Standard English name. Schneider's Tree Agamid Distribution and natural history. Dendragama schneideri occurs in high elevation, montane forest in north Sumatra's Bukit Barisan Mountain Range (Figs 1, 7). Our sampling encompasses Lake Toba (where the Toba blast occurred 71.6 kya) and the surrounding Karo Highlands. We hypothesize that the northern latitudinal limit of D. schneideri occurs at a break in the Barisan Range where elevation drops down into valley floors below 1000 m and does so continuously from the eastern to western edge of the range. That break in topography would certainly prohibit dispersal between montane sky islands today. Yet, maybe more importantly in terms of maintaining population differentiation, elevational ranges that drop below that point would have even prohibited dispersal during periods of glacial maxima during the Pleistocene epoch when montane forests retreated downslope 300-500 m (Hall 2009). The break occurs approximately at 3.196889°N 98.102583°E and we only found populations of D. dioidema north of that point. We hypothesize that the southern latitudinal limit of D. schneideri occurs approximately near the border of Sumatra Utara and Bengkulu. Similar to the northern limitation, there is a continuous elevational break that drops below 1000 m that extends across the width of the Barisan Range at the southern end of their range.
All the referred specimens were found sleeping in low vegetation 0.7-2.5 m above ground and between 1200-2800 m. All were found in montane cloud forest habitat and in higher elevations some individuals were found on moss covered vegetation in stunted forests where temperatures were slightly lower. Because the forests they inhabit receive high levels of precipitation relatively evenly distributed throughout cloud forests, D. schneideri are not dependent upon water bodies. Although some individuals were located along slopes near the edges of streams, those encounters did not seem to occur in any higher frequency than individuals being found in other areas away from bodies of water.
Many of the individuals collected were found along the edges of roadways and thus it seems they survive well along the edges of disturbed habitat, although they seem to be cloud forest obligates and require the presence of montane forest in some amount in order to persist in fragmented habitat.
Virtually nothing is known about the home range sizes and movement patterns of D. schneideri; however, given that they are relatively small arboreal lizards and live in cooler temperature habitats we would expect that they don't move long distances throughout the year and do not have large home ranges.
We documented distinct sexual dichromatism. Because the males tend to be more brightly colored, we suspect they use their brightly colored dewlaps to display during mating season. It is also likely they use their displays to defend territories, as do nearly all other brightly colored lizards (Sinervo and Lively 1996;Steffen and Guyer 2014). They may also use their colors when competing with other individuals. Although little is known about the life cycles of D. schneideri specifically, Harvey et al. (2017) suggest that other species of Dendragama reach maturity at around 60 mm, females lay 2-4 ovoid eggs and probably produce multiple clutches each year. This may also be true for D. schneideri given their other similarities in life histories. Little is known about the diet of D. schneideri; however, they are suspected to be insectivores.

Key to the species of Dendragama
We present a key to the species of Dendragama based on morphology and color pattern. Fig. 8 illustrates all four species. High supratemporal ridges enclosing a depressed parietal region, a row of white to yellow sublabial tubercles, and a visible tympanum immediately distinguish species of Dendragama from all other Sumatran agamids.  Shaney et al. (2020) breaks in species boundaries of Dendragama seem to occur consistently wherever elevation dips to 650-700 m (or lower) between mountains based on glacial periods during the Pleistocene. Thus, we hypothesize the following mountains will hold new D. dioidema sister lineages (Table 4).
Further south in latitude, the Karo highlands surrounding the Toba eruption site consist of topography that is continuously connected by higher elevation pieces of terrain and thus, there are likely few new sister species of D. schneideri throughout the Karo highlands between the latitudes 3.196889°N, 98.102583°E and 1.418139°N, 99.243389°E. However, there is distinct population subdivision among the populations that we sampled.
The break between populations of D. schneideri and D. boulengeri likely occurs around the point 1.418139N, 99.243389E which we hypothesize based on the dramatic drop-off in elevation that cuts through the entire width of the Barisan Range throughout the area. South of that point, the Barisan's topography becomes more variable in terms of having lower elevation valleys between mountains in a larger number of locations. Thus, we hypothesize the following nine mountains will have new sister species of D. boulengeri, including the mountains of the south Sapanuli Regency, Padang Lawas Regency, Rokan Hulu Regency, Mount Talakmau, Muaro Sungai Lolo of Pasaman Regency, Mount Sago, Mount Masurai, Bukit Daun, and Air Duku area of Rejang Lebong Regency. Not all of these ranges had clear names on the topographic maps, but for reference we provide GPS coordinates to the center point of the ranges in (Table 4).
The southern limitation of D. boulengeri populations probably occurs near the town of Curup where the mountains drop down very low in elevation, again separating D. boulengeri populations (and likely sister species) from D. australis populations. Among D. australis we hypothesize the following three mountains will hold new closely related species: Mount Kaba, Unknown strip of mountain in Empat Lawang Regency, Mountains in Sukabumi of Lahat Regency. The southern limitation of D. australis is not precisely known although we did not collect Dendragama further south in Latitude than Gunung Patah at -4.460662, 103.430325. We did collect Pseudocalotes cybelidermus and P. guttalineatus extensively throughout mountains further south, including Gunung Pesagi and montane forest above Ngarip and P. rhammanotus from Danau Ranau in Lampung Province. The absence of Dendragama from locations further south suggests populations of D. australis probably do not occur much further south than Gunung Patah. We hypothesize that each mountain probably has a new species, not just the possibility of one or two new species among all mountains mentioned. Based on Shaney et al. (2020) findings regarding the elevational shifts of Pleistocene cloud forests, it is likely that all of the 18 mountains mentioned remained isolated even during periods of glacial maxima. This latter point is crucial because Pleistocene connectivity is what would have facilitated gene flow historically and sympatric distributions contemporarily (which is not how they are distributed). Similar patterns have been seen in among peninsular Malaysia's herpetofauna as well (Loredo et al. 2013;Grismer et al. 2017). It is unlikely that Dendragama were able to disperse among any of the isolated mountains for millions of years, allowing for a whole array of species to diverge. Using these biogeographic patterns to identify likely locations for new species of Dendragama is certainly interesting on its own; however, perhaps more importantly, those same locations are probably where the distributions of many obligate cloud forest groups with limited dispersal (e.g. insects, rodents, amphibians) change in Sumatra. Thus, those mountains would be excellent candidates for other taxonomists conducting biological inventory throughout the region as well.
Additionally, a comparative biogeographic analysis among Dendragama and other agamid lizard groups would be a fascinating and informative study to better understand Sumatra's complex biological history. Grismer et al. (2016) provided a biogeographic study of draconid lizards of southeast Asia that yielded fascinating information; however, species of draconids from Sumatra remain underrepresented in biogeographic analyses (Inger and Voris 2001;Cannon et al. 2009). Furthermore, genetic splits between agamid lizard species seem quite old, and divergence dating of Sumatran agamid phylogenies may corroborate estimated geologic events from before the Pleistocene such as those put forth by (Voris 2000).

Phylogenetics
Phylogenetic analyses uncovered five distinct clades of Dendragama, among which are D. australis, D. dioidema, D. schneideri and two distinct clades of D. boulengeri, the first clade of D. boulengeri from Mount Kerinci and the second from near the type locality, Mount Marapi. These two clades are 5.0% pairwise genetically distant, which in many cases would constitute distinct species designations, if accompanied by readily identifiable morphological variation between populations. For example, Bradley and Baker (2009) show that sister species of mammals with readily distinguishable morphological characteristics typically had greater than 5.0% genetic variation in Cytochrome B sequences upon genetic evaluation later. Dendragama boulengeri populations examined in this study are on the edge of that genetic cutoff. Although means of some meristic and mensural characters differ between the Kerinci and Marapi populations of D. boulengeri, we did not find any fixed differences. Although geographically isolated from one another, we continue to recognize these two populations as D. boulengeri. However, we suspect that further studies may indeed identify these populations as distinct species.
These phylogenetic analyses further elucidate some of the patterns associated with cloud forest agamid lizard distributions in Sumatra. It seems that very few species have overlapping distributions and in contrast species, within Dendragama, species are distributed in a site-specific endemic pattern. Each species is distributed allopatrically with relatively small geographic ranges that match natural breaks in cloud forest habitat due to elevational changes. Furthermore, Lophocalotes is sister to the Dendragama group followed by insular Pseudocalotes, supporting Harvey et al. (2017).

Conservation
The conservation status of species of Dendragama has not yet been assessed by the IUCN (http://www.iucnredlist.org/), particularly because of the lack of population data. However, it is clear that Dendragama inhabits isolated cloud forest patches and Sumatra's rapid rates of deforestation are causing some of those patches to decline in size. Between 1900 and 2019 Sumatra's lowland forests were nearly wiped out for agriculture and timber, and because of forest fires (Poor et al. 2019). Deforestation is extending into highland areas that encompass remaining patches of D. schneideri habitat. Many forested areas have already been diminished beyond macroecological tipping points (Nowosad et al. 2019) and the ecological ramifications could cascade throughout the remaining cloud forest habitat.

Conclusions
The information from this paper contributes data on populations of D. boulengeri and D. schneideri, which we believe have small distributions, but are found in high abundance within those ranges. We estimated the latitudinal and longitudinal limits of D. boulengeri and D. schneideri (Fig. 1); however, we hypothesize that inventories of the rest of Sumatra's mountains would yield an array of new species and would show that the true distributional limits of all Dendragama species are actually small. These data may be used towards conducting IUCN Red List status assessments in the future; however additional information is needed to provide proper assessments. Regardless of their current conservation status, it is clear that rapidly increasing anthropogenic pressures throughout Sumatra are likely to have a significant impact on all species of Dendragama. Shaney et al. (2016) provide examples of how cryptic Indonesian lineages may be lost before being described and cryptic species may be overharvested due to poor taxonomic evaluation. Thus, continuation of biological inventory will be important in agamid lizard discovery and conservation of Sumatra's montane forest diversity in the near future. Given the rapid discovery of herpetofaunal diversity across Sumatra's highlands (Iskandar 2006;Kurniati et al. 2009) it is likely that an array of new agamid lizard species remains undiscovered throughout the region (Shaney et al. 2016).