A phylogeny-informed characterisation of global tetrapod traits addresses data gaps and biases.

Tetrapods (amphibians, reptiles, birds, and mammals) are model systems for global biodiversity science, but continuing data gaps, limited data standardisation, and ongoing flux in taxonomic nomenclature constrain integrative research on this group and potentially cause biased inference. We combined and harmonised taxonomic, spatial, phylogenetic, and attribute data with phylogeny-based multiple imputation to provide a comprehensive data resource (TetrapodTraits 1.0.0) that includes values, predictions, and sources for body size, activity time, micro- and macrohabitat, ecosystem, threat status, biogeography, insularity, environmental preferences, and human influence, for all 33,281 tetrapod species covered in recent fully sampled phylogenies. We assess gaps and biases across taxa and space, finding that shared data missing in attribute values increased with taxon-level completeness and richness across clades. Prediction of missing attribute values using multiple imputation revealed substantial changes in estimated macroecological patterns. These results highlight biases incurred by nonrandom missingness and strategies to best address them. While there is an obvious need for further data collection and updates, our phylogeny-informed database of tetrapod traits can support a more comprehensive representation of tetrapod species and their attributes in ecology, evolution, and conservation research.

Integrative species delimitation uncovers hidden diversity within the Pithecopus hypochondrialis species complex (Hylidae, Phyllomedusinae) and its phylogeography reveals Plio-Pleistocene connectivity among Neotropical savannas.

Despite their limited vagility and pronounced habitat heterogeneity in the tropics, many anuran species have unexpectedly extensive geographic ranges. One prominent example of this phenomenon is Pithecopus hypochondrialis, which is found in the Cerrado, Guianan savanna, and Llanos domains, as well as isolated tracts of savanna and open habitat within the Amazon Forest. The present study employs an integrative species delimitation approach to test the hypothesis that P. hypochondrialis is in fact a species complex. We also reconstruct the relationships among the lineages delimited here and other Pithecopus species. In this study, we employ Ecological Niche Modelling (ENM) and spatiotemporal phylogeographic reconstruction approaches to evaluate a multitude of scenarios of connectivity across the Neotropical savannas. We identified three divergent lineages, two of which have been described previously. The lineages were allocated to a lowland Pithecopus clade, although the relationships among these lineages are weakly supported. Both the ENM and the phylogeographic reconstruction highlight the occurrence of periods of connectivity among the Neotropical savannas over the course of the Pliocene and Pleistocene epochs. These processes extended from eastern Amazonia to the northern coast of Brazil. The findings of the present study highlight the presence of hidden diversity within P. hypochondrialis, and reinforce the need for a comprehensive taxonomic review. These findings also indicate intricate and highly dynamic patterns of connectivity across the Neotropical savannas that date back to the Pliocene.

Diversification of the terrestrial frog genus Anomaloglossus (Anura, Aromobatidae) in the Guiana Shield proceeded from highlands to lowlands, with successive loss and reacquisition of endotrophy.

Two main landscapes emerge from the Guiana Shield: the highlands to the west called the Pantepui region and the Amazonian lowlands to the east, both harbouring numerous endemic species. With 32 currently recognized species, the genus Anomaloglossus stands out among Neotropical frogs as one that diversified only within the Guiana Shield both in the highlands and the lowlands. We present a time-calibrated phylogeny obtained by using combined mitogenomic and nuclear DNA, which suggests that the genus originates from Pantepui where extant lineages started diversifying around 21 Ma, and subsequently (ca. 17 Ma) dispersed during the Miocene Climatic Optimum to the lowlands of the eastern Guiana Shield where the ability to produce endotrophic tadpoles evolved. Further diversification within the lowlands in the A. stepheni group notably led to an evolutionary reversal toward exotrophy in one species group during the late Miocene, followed by reacquisition of endotrophy during the Pleistocene. These successive shifts of reproductive mode seem to have accompanied climatic oscillations. Long dry periods might have triggered evolution of exotrophy, whereas wetter climates favoured endotrophic forms, enabling colonization of terrestrial habitats distant from water. Acquisition, loss, and reacquisition of endotrophy makes Anomaloglossus unique among frogs and may largely explain the current species diversity. The micro evolutionary processes involved in these rapid shifts of reproductive mode remain to be revealed.

Genetic structure and landscape effects on gene flow in the Neotropical lizard Norops brasiliensis (Squamata: Dactyloidae).

One key research goal of evolutionary biology is to understand the origin and maintenance of genetic variation. In the Cerrado, the South American savanna located primarily in the Central Brazilian Plateau, many hypotheses have been proposed to explain how landscape features (e.g., geographic distance, river barriers, topographic compartmentalization, and historical climatic fluctuations) have promoted genetic structure by mediating gene flow. Here, we asked whether these landscape features have influenced the genetic structure and differentiation in the lizard species Norops brasiliensis (Squamata: Dactyloidae). To achieve our goal, we used a genetic clustering analysis and estimate an effective migration surface to assess genetic structure in the focal species. Optimized isolation-by-resistance models and a simulation-based approach combined with machine learning (convolutional neural network; CNN) were then used to infer current and historical effects on population genetic structure through 12 unique landscape models. We recovered five geographically distributed populations that are separated by regions of lower-than-expected gene flow. The results of the CNN showed that geographic distance is the sole predictor of genetic variation in N. brasiliensis, and that slope, rivers, and historical climate had no discernible influence on gene flow. Our novel CNN approach was accurate (89.5%) in differentiating each landscape model. CNN and other machine learning approaches are still largely unexplored in landscape genetics studies, representing promising avenues for future research with increasingly accessible genomic datasets.

Diversity, biogeography, and reproductive evolution in the genus Pipa (Amphibia: Anura: Pipidae).

The genus Pipa is a species-poor clade of Neotropical frogs and one of the most bizarre-looking due to many highly derived anatomical traits related to their fully aquatic lifestyle. With their African relatives, they form the Pipidae family, which has attracted much attention, especially regarding its anatomy, reproductive biology, paleontology and biogeography. However, the actual diversity and phylogenetic relationships within Pipa remain poorly understood, and thus so do their historical biogeography and the evolution of striking features, such as the absence of teeth and endotrophy in some species. Using short mtDNA sequences across the distribution of the genus, we identified 15 main lineages (Operational Taxonomic Units - OTUs). This more than doubles the number of the currently seven valid nominal species. Several closely related OTUs do not share nuDNA alleles, confirming species divergence. Time-calibrated phylogenies obtained from mitogenomes and from 10 nuclear loci provide highly similar topologies but strikingly distinct node ages for Pipa. High dN/dS ratios and the variation of substitution rates across the trees suggest a strong effect of saturation on fast evolving positions of mtDNA, producing a substantially shorter stem branch of Pipa. Focusing on the nuDNA topology, we inferred an early Neogene Amazonian origin of the diversification of Pipa, with an initial split between the Guiana-Brazilian Shields and Western Amazonia, a pattern observed in many other co-distributed groups. All the western species are edentate, suggesting a single loss in the genus. Each of these groups diversified further out of Amazonia, toward the Atlantic Forest and toward trans-Andean forests, respectively. These events are concomitant with paleogeographic changes and match patterns observed in other co-distributed taxonomic groups. The two Amazonian lineages have probably independently acquired endotrophic larval development.

Diversification history of clown tree frogs in Neotropical rainforests (Anura, Hylidae, Dendropsophus leucophyllatus group).

General consensus emphasizes that no single biological process can explain the patterns of species' distributions and diversification in the Neotropics. Instead, the interplay of several processes across space and time must be taken into account. Here we investigated the phylogenetic relationships and biogeographic history of tree frogs in the Dendropsophus leucophyllatus species group (Amphibia: Hylidae), which is distributed across Amazonia and the Atlantic rainforests. Using Next Generation Sequencing (NGS) and double digest restriction-site associated DNA (ddRADseq), we inferred phylogenetic relationships, species limits, and temporal and geographic patterns of diversification relative to the history of these biomes. Our results indicate that the D. leucophyllatus species group includes at least 14 independent lineages, which are currently arranged into ten described species. Therefore, a significant portion of species in the group are still unnamed. Different processes were associated to the group diversification history. For instance, the Andes uplift likely caused allopatric speciation for Cis-Andean species, whereas it may also be responsible for changes in the Amazonian landscape triggering parapatric speciation by local adaptation to ecological factors. Meanwhile, Atlantic Forest ancestors unable to cross the dry diagonal biomes after rainforest's retraction, evolved in isolation into different species. Diversification in the group began in the early Miocene, when connections between Atlantic Forest and the Andes (Pacific Dominion) by way of a south corridor were possible. The historical scenario in Amazonia, characterized by several speciation events and habitat heterogeneity, helped promoting diversification, resulting in the highest species diversity for the group. This marked species diversification did not happen in Atlantic Forest, where speciation is very recent (late Pliocene and Pleistocene), despite its remarkable climatic heterogeneity.

The evolutionary history of Lygodactylus lizards in the South American open diagonal.

The Pleistocenic Arc Hypothesis (PAH) posits that South American Seasonally Dry Tropical Forests (SDTF) were interconnected during Pleistocene glacial periods, enabling the expansion of species ranges that were subsequently fragmented in interglacial periods, promoting speciation. The lizard genus Lygodactylus occurs in Africa, Madagascar, and South America. Compared to the high diversity of African Lygodactylus, only two species are known to occur in South America, L. klugei and L. wetzeli, distributed in SDTFs and the Chaco, respectively. We use a phylogenetic approach based on mitochondrial (ND2) and nuclear (RAG-1) markers covering the known range of South American Lygodactylus to investigate (i) if they are monophyletic relative to their African congeners, (ii) if their divergence is congruent with the fragmentation of the PAH, and (iii) if cryptic diversity exists within currently recognized species. Maximum likelihood and Bayesian phylogenetic analyses recovered a well-supported monophyletic South American Lygodactylus, presumably resulting from a single trans-Atlantic dispersal event 29 Mya. Species delimitation analyses supported the existence of five putative species, three of them undescribed. Divergence times among L. klugei and the three putative undescribed species, all endemic to the SDTFs, are not congruent with the fragmentation of the PAH. However, fragmentation of the once broader and continuous SDTFs likely influenced the divergence of L. wetzeli in the Chaco and Lygodactylus sp. 3 (in a SDTF enclave in the Cerrado).

Extinction risks forced by climatic change and intraspecific variation in the thermal physiology of a tropical lizard.

Temperature increases can impact biodiversity and predicting their effects is one of the main challenges facing global climate-change research. Ectotherms are sensitive to temperature change and, although predictions indicate that tropical species are highly vulnerable to global warming, they remain one of the least studied groups with respect to the extent of physiological variation and local extinction risks. We model the extinction risks for a tropical heliothermic teiid lizard (Kentropyx calcarata) integrating previously obtained information on intraspecific phylogeographic structure, eco-physiological traits and contemporary species distributions in the Amazon rainforest and its ecotone to the Cerrado savannah. We also investigated how thermal-biology traits vary throughout the species' geographic range and the consequences of such variation for lineage vulnerability. We show substantial variation in thermal tolerance of individuals among thermally distinct sites. Thermal critical limits were highly correlated with operative environmental temperatures. Our physiological/climatic model predicted relative extinction risks for local populations within clades of K. calcarata for 2050 ranging between 26.1% and 70.8%, while for 2070, extinction risks ranged from 52.8% to 92.8%. Our results support the hypothesis that tropical-lizard taxa are at high risk of local extinction caused by increasing temperatures. However, the thermo-physiological differences found across the species' distribution suggest that local adaptation may allow persistence of this tropical ectotherm in global warming scenarios. These results will serve as basis to further research to investigate the strength of local adaptation to climate change. Persistence of Kentropyx calcarata also depends on forest preservation, but the Amazon rainforest is currently under high deforestation rates. We argue that higher conservation priority is necessary so the Amazon rainforest can fulfill its capacity to absorb the impacts of temperature increase on tropical ectotherms during climate change.

Estimating synchronous demographic changes across populations using hABC and its application for a herpetological community from northeastern Brazil.

Many studies propose that Quaternary climatic cycles contracted and/or expanded the ranges of species and biomes. Strong expansion-contraction dynamics of biomes presume concerted demographic changes of associated fauna. The analysis of temporal concordance of demographic changes can be used to test the influence of Quaternary climate on diversification processes. Hierarchical approximate Bayesian computation (hABC) is a powerful and flexible approach that models genetic data from multiple species, and can be used to estimate the temporal concordance of demographic processes. Using available single-locus data, we can now perform large-scale analyses, both in terms of number of species and geographic scope. Here, we first compared the power of four alternative hABC models for a collection of single-locus data. We found that the model incorporating an a priori hypothesis about the timing of simultaneous demographic change had the best performance. Second, we applied the hABC models to a data set of seven squamate and four amphibian species occurring in the Seasonally Dry Tropical Forests (Caatinga) in northeastern Brazil, which, according to paleoclimatic evidence, experienced an increase in aridity during the Pleistocene. If this increase was important for the diversification of associated xeric-adapted species, simultaneous population expansions should be evident at the community level. We found a strong signal of synchronous population expansion in the Late Pleistocene, supporting the increase of the Caatinga during this time. This expansion likely enhanced the formation of communities adapted to high aridity and seasonality and caused regional extirpation of taxa adapted to wet forest.

Cryptic lineages and diversification of an endemic anole lizard (Squamata, Dactyloidae) of the Cerrado hotspot.

The Cerrado is a wide Neotropical savanna with tremendously high endemic diversity. Yet, it is not clear what the prevalent processes leading to such diversification are. We used the Cerrado-endemic lizard Norops meridionalis to investigate the main abiotic factors that promoted genetic divergence, the timings of these divergence events, and how these relate to cryptic diversity in the group. We sequenced mitochondrial and nuclear genes from 21 sites of N. meridionalis to generate species tree, divergence time estimations, and estimate species limits. We also performed population-level analysis and estimated distribution models to test the roles of niche conservatism and divergence in the group diversification. We found that N. meridionalis is composed by at least five cryptic species. Divergence time estimations suggest that the deepest branches split back into the early-mid Miocene, when most of the geophysical activity of the Cerrado took place. The deep divergences found in N. meridionalis suggest that beta anoles invaded South America much earlier than previously thought. Recent published evidence supports this view, indicating that the Panama gap closed as early as 15 mya, allowing for an early invasion of Norops into South America. The spatial pattern of diversification within N. meridionalis follows a northwest-southeast direction, which is consistent across several species of vertebrates endemic to the Cerrado. Also, we found evidence for non-stationary isolation by distance, which occurs when genetic differentiation depends on space. Our preliminary data in two out of five lineages suggest that niche conservatism is an important mechanism that promoted geographic fragmentation in the group.

Biogeographic history and cryptic diversity of saxicolous Tropiduridae lizards endemic to the semiarid Caatinga.

BACKGROUND: Phylogeographic research has advanced in South America, with increasing efforts on taxa from the dry diagonal biomes. However, the diversification of endemic fauna from the semiarid Caatinga biome in northeastern Brazil is still poorly known. Here we targeted saxicolous lizards of the Tropidurus semitaeniatus species group to better understand the evolutionary history of these endemic taxa and the Caatinga. We estimated a time-calibrated phylogeny for the species group based on two mitochondrial and two nuclear genes and jointly estimated the species limits and species tree within the group. We also devoted a denser phylogeographic sampling of the T. semitaeniatus complex to explore migration patterns, and the spatiotemporal diffusion history to verify a possible role of the São Francisco River as a promoter of differentiation in this saxicolous group of lizards. RESULTS: Phylogenetic analysis detected high cryptic genetic diversity, occurrence of unique microendemic lineages associated with older highlands, and a speciation history that took place during the Pliocene-Pleistocene transition. Species delimitation detected five evolutionary entities within the T. semitaeniatus species group, albeit with low support. Thus, additional data are needed for a more accurate definition of species limits and interspecific relationships within this group. Spatiotemporal analyses reconstructed the geographic origin of the T. semitaeniatus species complex to be located north of the present-day course of the São Francisco River, followed by dispersal that expanded its distribution towards the northwest and south. Gene flow estimates showed higher migration rates into the lineages located north of the São Francisco River. CONCLUSIONS: The phylogenetic and population structures are intrinsically associated with stable rock surfaces and landscape rearrangements, such as the establishment of drainage basins located to the northern and southern distribution ranges. The T. semitaeniatus complex preserved high genetic diversity during range expansion, possibly as a result of frequent long-distance dispersal events. Our results indicate that both the current course of the São Francisco River and its paleo-courses had an important role in promoting diversification of the Caatinga endemic T. semitaeniatus species group.

Toe pad morphology and adhesion in the miniaturized gecko, Chatogekko amazonicus (Gekkota: Sphaerodactylidae).

Chatogekko amazonicus is a miniaturized gecko from northern South America and is among the smallest of toe pad bearing lizards. The toe pads of C. amazonicus are miniscule, between 18% and 27% of the plantar surface area. We aimed to investigate the relationship between adhesive toe pad morphology, body size, and adhesive capabilities. Using scanning electron microscopy, we determine that the adhesive pads of C. amazonicus exhibit branched setae similar to those of other geckos, but that are generally much smaller. When compared with other gecko taxa, we show that C. amazonicus setae occupy a similar range of seta length: snout-vent length ratio and aspect ratio as other gekkonoid species (i.e. Gekkonidae, Phyllodactylidae, and Sphaerodactylidae). We demonstrate that C. amazonicus-even with its relatively small toe pads-is capable of climbing a smooth glass surface at a nearly vertical angle. We suggest that sphaerodactylids like C. amazonicus offer an excellent system for studying toe pad morphology and function in relation to miniaturization.

The role of environmental gradients and microclimates in structuring communities and functional groups of lizards in a rainforest-savanna transition area.

Environmental heterogeneity poses a significant influence on the functional characteristics of species and communities at local scales. Environmental transition zones, such as at the savanna-forest borders, can act as regions of ecological tension when subjected to sharp variations in the microclimate. For ectothermic organisms, such as lizards, environmental temperatures directly influence physiological capabilities, and some species use different thermoregulation strategies that produce varied responses to local climatic conditions, which in turn affect species occurrence and community dynamics. In the context of global warming, these various strategies confer different types of vulnerability as well as risks of extinction. To assess the vulnerability of a species and understand the relationships between environmental variations, thermal tolerance of a species and community structure, lizard communities in forest-savanna transition areas of two national parks in the southwestern Amazon were sampled and their thermal functional traits were characterized. Then, we investigated how community structure and functional thermal variation were shaped by two environmental predictors (i.e., microclimates estimated locally and vegetation structure estimated from remote sensing). It was found that the community structure was more strongly predicted by the canopy surface reflectance values obtained via remote sensing than by microclimate variables. Environmental temperatures were not the most important factor affecting the occurrence of species, and the variations in ecothermal traits demonstrated a pattern within the taxonomic hierarchy at the family level. This pattern may indicate a tendency for evolutionary history to indirectly influence these functional features. Considering the estimates of the thermal tolerance range and warming tolerance, thermoconformer lizards are likely to be more vulnerable and at greater risk of extinction due to global warming than thermoregulators. The latter, more associated with open environments, seem to take advantage of their lower vulnerability and occur in both habitat types across the transition, potentially out-competing and further increasing the risk of extinction and vulnerability of forest-adapted thermoconformer lizards in these transitional areas.

Quaternary range and demographic expansion of Liolaemus darwinii (Squamata: Liolaemidae) in the Monte Desert of Central Argentina using Bayesian phylogeography and ecological niche modelling.

Until recently, most phylogeographic approaches have been unable to distinguish between demographic and range expansion processes, making it difficult to test for the possibility of range expansion without population growth and vice versa. In this study, we applied a Bayesian phylogeographic approach to reconstruct both demographic and range expansion in the lizard Liolaemus darwinii of the Monte Desert in Central Argentina, during the Late Quaternary. Based on analysis of 14 anonymous nuclear loci and the cytochrome b mitochondrial DNA gene, we detected signals of demographic expansion starting at ~55 ka based on Bayesian Skyline and Skyride Plots. In contrast, Bayesian relaxed models of spatial diffusion suggested that range expansion occurred only between ~95 and 55 ka, and more recently, diffusion rates were very low during demographic expansion. The possibility of population growth without substantial range expansion could account for the shared patterns of demographic expansion during the Last Glacial Maxima (OIS 2 and 4) in fish, small mammals and other lizards of the Monte Desert. We found substantial variation in diffusion rates over time, and very high rates during the range expansion phase, consistent with a rapidly advancing expansion front towards the southeast shown by palaeo-distribution models. Furthermore, the estimated diffusion rates are congruent with observed dispersal rates of lizards in field conditions and therefore provide additional confidence to the temporal scale of inferred phylogeographic patterns. Our study highlights how the integration of phylogeography with palaeo-distribution models can shed light on both demographic and range expansion processes and their potential causes.

Seasonal flooding decreases fruit-feeding butterfly species dominance and increases spatial turnover in floodplain forests of central Amazonia.

The seasonal flood pulse in Amazonia can be considered a primary driver of community structure in floodplain environments. Although this natural periodic disturbance is part of the landscape dynamics, the seasonal inundation presents a considerable challenge to organisms that inhabit floodplain forests. The present study investigated the effect of seasonal flooding on fruit-feeding butterfly assemblages in different forest types and strata in central Amazonia. We sampled fruit-feeding butterflies in the canopy and the understory using baited traps in adjacent upland (unflooded forests-terra firme), white and blackwater floodplain forests (várzea and igapó, respectively) during the low- and high-water seasons. Butterfly abundance decreased in the high-water season, especially of dominant species in várzea, but the number of species was similar between seasons in the three forest types. Species composition differed between strata in all forest types. However, the flood pulse only affected butterfly assemblages in várzea forest. The ß-diversity components also differed only in várzea. Species replacement (turnover) dominated the spatial ß-diversity in igapó and terra firme in both seasons and várzea in the high-water season. Nonetheless, nestedness was relatively higher in várzea forests during the low-water season, mainly due to the effect of dominant species. These results emphasize the importance of seasonal flooding to structure butterfly assemblages in floodplain forests and reveal the idiosyncrasy of butterfly community responses to flooding in different forest types. Our results also suggest that any major and rapid changes to the hydrological regime could severely affect floodplain communities adapted to this natural seasonal hydrological cycle, threatening the existence of these unique environments.

Phylogeny and cryptic diversity in geckos (Phyllopezus; Phyllodactylidae; Gekkota) from South America's open biomes.

The gecko genus Phyllopezus occurs across South America's open biomes: Cerrado, Seasonally Dry Tropical Forests (SDTF, including Caatinga), and Chaco. We generated a multi-gene dataset and estimated phylogenetic relationships among described Phyllopezus taxa and related species. We included exemplars from both described Phyllopezus pollicaris subspecies, P. p. pollicaris and P. p.przewalskii. Phylogenies from the concatenated data as well as species trees constructed from individual gene trees were largely congruent. All phylogeny reconstruction methods showed Bogertia lutzae as the sister species of Phyllopezus maranjonensis, rendering Phyllopezus paraphyletic. We synonymized the monotypic genus Bogertia with Phyllopezus to maintain a taxonomy that is isomorphic with phylogenetic history. We recovered multiple, deeply divergent, cryptic lineages within P. pollicaris. These cryptic lineages possessed mtDNA distances equivalent to distances among other gekkotan sister taxa. Described P. pollicaris subspecies are not reciprocally monophyletic and current subspecific taxonomy does not accurately reflect evolutionary relationships among cryptic lineages. We highlight the conservation significance of these results in light of the ongoing habitat loss in South America's open biomes.

Two new species of geckos of the genus Phyllopezus Peters, 1878 (Squamata: Gekkota: Phyllodactylidae) from northeastern Brazil.

We describe two new species of Brazilian geckos of the genus Phyllopezus based on morphological and molecular data. The first species is currently known from a relictual Cerrado enclavecampos rupestres, in the mountains of the Serra do Espinhao in the state of Bahia. The second species is known from northeastern Atlantic Forest and transitional areas with Caatinga biome in the state of Alagoas. The two new species are sister taxa and together are the sister clade to the remaining species in the Phyllopezus pollicaris species complex. These new species can be morphologically distinguished from their congeners by meristic and morphometric characters, in addition to color pattern and genetic differentiation.

Phylogeographic model selection using convolutional neural networks.

The discipline of phylogeography has evolved rapidly in terms of the analytical toolkit used to analyse large genomic data sets. Despite substantial advances, analytical tools that could potentially address the challenges posed by increased model complexity have not been fully explored. For example, deep learning techniques are underutilized for phylogeographic model selection. In non-model organisms, the lack of information about their ecology and evolution can lead to uncertainty about which demographic models are appropriate. Here, we assess the utility of convolutional neural networks (CNNs) for assessing demographic models in South American lizards in the genus Norops. Three demographic scenarios (constant, expansion, and bottleneck) were considered for each of four inferred population-level lineages, and we found that the overall model accuracy was higher than 98% for all lineages. We then evaluated a set of 26 models that accounted for evolutionary relationships, gene flow, and changes in effective population size among the four lineages, identifying a single model with an estimated overall accuracy of 87% when using CNNs. The inferred demography of the lizard system suggests that gene flow between non-sister populations and changes in effective population sizes through time, probably in response to Pleistocene climatic oscillations, have shaped genetic diversity in this system. Approximate Bayesian computation (ABC) was applied to provide a comparison to the performance of CNNs. ABC was unable to identify a single model among the larger set of 26 models in the subsequent analysis. Our results demonstrate that CNNs can be easily and usefully incorporated into the phylogeographer's toolkit.