Methodological congruence in phylogenomic analyses with morphological support for teiid lizards (Sauria: Teiidae)

Central African dwarf crocodiles found in syntopy are comparably divergent to South American dwarf caimans

Recent molecular taxonomic advancements have expanded our understanding of crocodylian diversity, revealing the existence of previously overlooked species, including the Congo dwarf crocodile (Osteolaemus osborni) in the central Congo Basin rainforests. This study explores the genomic divergence between O. osborni and its better-known relative, the true dwarf crocodile (Osteolaemus tetraspis), shedding light on their evolutionary history. Field research conducted in the northwestern Republic of the Congo uncovered a locality where both species coexist in sympatry/syntopy. Genomic analysis of sympatric individuals reveals a level of divergence comparable to that between ecologically similar South American dwarf caimans (Paleosuchus palpebrosus and Paleosuchus trigonatus), suggesting parallel speciation in the Afrotropics and Neotropics during the Middle to Late Miocene, 10-12 Ma. Comparison of the sympatric and allopatric dwarf crocodiles indicates no gene flow between the analysed sympatric individuals of O. osborni and O. tetraspis. However, a larger sample will be required to answer the question of whether or to what extent these species hybridize. This study emphasizes the need for further research on the biology and conservation status of the Congo dwarf crocodile, highlighting its significance in the unique biodiversity of the Congolian rainforests and thus its potential as a flagship species.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Nuclear phylogenomics, but not mitogenomics, resolves the most successful Late Miocene radiation of African mammals (Rodentia: Muridae: Arvicanthini)

The tribe Arvicanthini (Muridae: Murinae) is a highly diversified group of rodents (ca. 100 species) and with 18 African genera (plus one Asiatic) represents probably the most successful adaptive radiation of extant mammals in Africa. They colonized a broad spectrum of habitats (from rainforests to semi-deserts) in whole sub-Saharan Africa and their members often belong to most abundant parts of mammal communities. Despite intensive efforts, the phylogenetic relationships among major lineages (i.e. genera) remained obscured, which was likely caused by the intensive radiation of the group, dated to the Late Miocene. Here we used genomic scale data (377 nuclear loci; 581,030 bp) and produced the first fully resolved species tree containing all currently delimited genera of the tribe. Mitogenomes were also extracted, and while the results were largely congruent, there was less resolution at basal nodes of the mitochondrial phylogeny. Results of a fossil-based divergence dating analysis suggest that the African radiation started early after the colonization of Africa by a single arvicanthine ancestor from Asia during the Messinian stage (ca. 7 Ma), and was likely linked with a fragmentation of the pan-African Miocene forest. Some lineages remained in the rain forest, while many others successfully colonized broad spectrum of new open habitats (e.g. savannas, wetlands or montane moorlands) that appeared at the beginning of Pliocene. One lineage even evolved partially arboricolous life style in savanna woodlands, which allowed them to re-colonize equatorial forests. We also discuss delimitation of genera in Arvicanthini and propose corresponding taxonomic changes.

Unusual lizard fossil from the Miocene of Nebraska and a minimum age for cnemidophorine teiids

Teiid lizards are well represented in the fossil record and are common components of modern ecosystems in North and South America. Many fossils were referred to the cnemidophorine teiid group (whiptails, racerunners and relatives), particularly from North America. However, systematic interpretations of morphological features in cnemidophorines were hampered by the historically problematic taxonomy of the clade, and the biogeography and chronology of cnemidophorine evolution in North America is poorly understood from the fossil record. Few fossil cnemidophorines were identified with an apomorphy-based diagnosis, and there are almost no fossil cnemidophorines that could be used to anchor node calibrations. Here, I describe a cnemidophorine from the Miocene Ogallala Group of Nebraska and diagnose the fossil using apomorphies. In that process, I clarify the systematic utility of several morphological features of cnemidophorine lizards. I refer the fossil to the least inclusive clade containing Aspidoscelis, Holcosus and Pholidoscelis. The most conservative minimum age of the locality of the fossil is a fission-track date of 6.3 Ma, but mammal biochronology provides a more refined age of 9.4 Ma, which can be used as a minimum age for the crown cnemidophorine clade in divergence time analyses. The fossil indicates that a cnemidophorine lineage that does not live in Nebraska today inhabited the area during the Miocene. I refrain from naming a new taxon pending discovery of additional fossil material of the lineage to which the fossil belonged.

Aspidoscelis costatus costatus (Squamata, Teiidae): high elevation clutch production for a population of whiptail lizards

Clutch size and number of clutches per reproductive cycle are important life history traits that can be influenced by anatomical, physiological, evolutionary, and ecological factors. This report on the clutch size and number of clutches of an endemic Mexican whiptail lizard, Aspidoscelis costatus costatus (Cope, 1878), is based on a study of population at an unusually high elevation for a member of this genus. The study site is located in Ixtapan de la Sal, southeastern Estado de México, Central Mexico, at 2090 m a.s.l. Lizards were sampled in June 2006, and from May to July 2007, where females of Aspidoscelis costatus costatus were collected by hand along a drift fence. Female reproductive condition was evaluated based on abdominal palpation for presence of developing eggs; clutch size was determined by actual counts of either vitellogenic follicles or oviductal eggs. The smallest reproductive female was 77 mm snout-vent length; females produced a minimum of two clutches during the breeding season, the mean clutch size of 6.5 eggs (n = 33) was one of the largest reported for the genus. However, both length and width of its eggs, and the relative clutch mass have not been diminished by development of a large clutch. Additionally, comparisons of clutch size were undertaken within the polytypic A. costatus complex, within the genus Aspidoscelis, and between certain genera of whiptail lizards. This apparently represents the first study of whiptail lizards (genus Aspidoscelis), assessing the aforementioned reproductive characteristics, in a population above 2000 m.

Optimizing Phylogenomics with Rapidly Evolving Long Exons: Comparison with Anchored Hybrid Enrichment and Ultraconserved Elements

Marker selection has emerged as an important component of phylogenomic study design due to rising concerns of the effects of gene tree estimation error, model misspecification, and data-type differences. Researchers must balance various trade-offs associated with locus length and evolutionary rate among other factors. The most commonly used reduced representation data sets for phylogenomics are ultraconserved elements (UCEs) and Anchored Hybrid Enrichment (AHE). Here, we introduce Rapidly Evolving Long Exon Capture (RELEC), a new set of loci that targets single exons that are both rapidly evolving (evolutionary rate faster than RAG1) and relatively long in length (>1,500 bp), while at the same time avoiding paralogy issues across amniotes. We compare the RELEC data set to UCEs and AHE in squamate reptiles by aligning and analyzing orthologous sequences from 17 squamate genomes, composed of 10 snakes and 7 lizards. The RELEC data set (179 loci) outperforms AHE and UCEs by maximizing per-locus genetic variation while maintaining presence and orthology across a range of evolutionary scales. RELEC markers show higher phylogenetic informativeness than UCE and AHE loci, and RELEC gene trees show greater similarity to the species tree than AHE or UCE gene trees. Furthermore, with fewer loci, RELEC remains computationally tractable for full Bayesian coalescent species tree analyses. We contrast RELEC to and discuss important aspects of comparable methods, and demonstrate how RELEC may be the most effective set of loci for resolving difficult nodes and rapid radiations. We provide several resources for capturing or extracting RELEC loci from other amniote groups.

Steroid hormones, energetic state, and immunocompetence vary across reproductive contexts in a parthenogenetic lizard

Reproduction is energetically expensive and investing in this life history trait is likely accompanied by significant changes in physiological activity. Investment strategy necessary for achieving reproductive success in reptiles can vary with reproductive form and pattern, potentiating different consequences for competing fitness-related traits such as those key to survival. The goal of this study was to assess if and how energetic state (i.e., energy metabolites) and self-maintenance (i.e., immunocompetence) are hormonally modulated across reproductive contexts in an oviparous, parthenogenetic lizard, the Colorado Checkered Whiptail Aspidoscelis neotesselata. Here blood plasma samples were collected from lizards within the US Army Fort Carson Military Installation near Colorado Springs, CO, USA, during seasons of reproductive activity (i.e., June) and inactivity (i.e., August). Measures of reproductive (i.e., estradiol) and energy-mobilizing (i.e., corticosterone) hormones, energy metabolites (i.e., glucose, triglycerides, and free glycerol), and innate immunity (i.e., bactericidal ability) were compared by season and reproductive stage. Levels of energy metabolites and bactericidal ability were compared to levels of E₂ and CORT. Bactericidal ability was also compared to levels of energy metabolites. Corticosterone and glucose levels were lower during the reproductive season while triglyceride levels and bactericidal ability were higher, but both estradiol and free glycerol levels did not differ between seasons. Throughout vitellogenesis, corticosterone and glucose levels as well as bactericidal ability did not differ, but estradiol levels were higher during early and mid-stage and both triglyceride and free glycerol levels were lower during gravidity. Corticosterone levels were negatively associated with circulating triglycerides and bactericidal ability, but were not related to glucose nor free glycerol levels. Estradiol levels were positively associated with free glycerol levels and bactericidal ability, but were not related to glucose nor triglyceride levels. Finally, bactericidal ability was negatively associated with glucose, but positively associated with triglycerides. Differences in energetic state and immunocompetence are thus reflected by shifts in hormone secretion across reproductive investment. These findings provide partial support for the hypothesis that energetic state is differentially regulated by steroid hormones to afford reproduction, potentially at the cost of future survival.

The combined role of dispersal and niche evolution in the diversification of Neotropical lizards

Ecological requirements and environmental conditions can influence diversification across temporal and spatial scales. Understanding the role of ecological niche evolution under phylogenetic contexts provides insights on speciation mechanisms and possible responses to future climatic change. Large-scale phyloclimatic studies on the megadiverse Neotropics, where biomes with contrasting vegetation types occur in narrow contact, are rare. We integrate ecological and biogeographic data with phylogenetic comparative methods, to investigate the relative roles of biogeographic events and niche divergence and conservatism on the diversification of the lizard genus Kentropyx Spix, 1825 (Squamata: Teiidae), distributed in South American rainforests and savannas. Using five molecular markers, we estimated a dated species tree, which recovered three clades coincident with previously proposed species groups diverging during the mid-Miocene. Biogeography reconstruction indicates a role of successive dispersal events from an ancestral range in the Brazilian Shield and western Amazonia. Ancestral reconstruction of climatic tolerances and niche overlap metrics indicates a trend of conservatism during the diversification of groups from the Amazon Basin and Guiana Shield, and a strong signal of niche divergence in the Brazilian Shield savannas. Our results suggest that climatic-driven divergence at dynamic forest-savanna borders might have resulted in adaptation to new environmental niches, promoting habitat shifts and shaping speciation patterns of Neotropical lizards. Dispersal and ecological divergence could have a more important role in Neotropical diversification than previously thought.

The Karyotype of Salvator merianae (Squamata, Teiidae): Analyses by Classical and Molecular Cytogenetic Techniques

In this study, we analyzed the karyotype of Salvator merianae (Teiidae) from the Brazilian semiarid region using different cytogenetic markers. Chromosomes were examined by classical (Giemsa and AgNOR staining) and molecular (FISH with ribosomal, telomeric, and microsatellite probes) cytogenetic approaches. S. merianae showed a diploid chromosome number of 2n = 38 (10 biarmed macrochromosomes + 28 microchromosomes). No sex-linked chromosome heteromorphisms were observed. Clusters of 18S/28S rDNA were localized in the terminal region of the long arm of pair 2. In addition to the typical telomeric signals, (TTAGGG)n repeats were detected in the pericentromeric region of some macrochromosome pairs, which might indicate the occurrence of chromosomal rearrangements via chromosome fusions. Hybridization signals of the microsatellite probes (GA)n, (GAA)n, and (GAG)n were uniformly distributed across all chromosomes, while (CA)n, (CAA)n, and (CAC)n produced brighter signals in the telomeric and pericentromeric regions of specific chromosome pairs. The comparison with previous studies demonstrates that, despite the wide distribution of S. merianae, the macrostructure organization of the karyotype remained unchanged, showing stability in diploid number and chromosome morphology.

Interrogating Genomic-Scale Data for Squamata (Lizards, Snakes, and Amphisbaenians) Shows no Support for Key Traditional Morphological Relationships

Genomics is narrowing uncertainty in the phylogenetic structure for many amniote groups. For one of the most diverse and species-rich groups, the squamate reptiles (lizards, snakes, and amphisbaenians), an inverse correlation between the number of taxa and loci sampled still persists across all publications using DNA sequence data and reaching a consensus on the relationships among them has been highly problematic. In this study, we use high-throughput sequence data from 289 samples covering 75 families of squamates to address phylogenetic affinities, estimate divergence times, and characterize residual topological uncertainty in the presence of genome-scale data. Importantly, we address genomic support for the traditional taxonomic groupings Scleroglossa and Macrostomata using novel machine-learning techniques. We interrogate genes using various metrics inherent to these loci, including parsimony-informative sites (PIS), phylogenetic informativeness, length, gaps, number of substitutions, and site concordance to understand why certain loci fail to find previously well-supported molecular clades and how they fail to support species-tree estimates. We show that both incomplete lineage sorting and poor gene-tree estimation (due to a few undesirable gene properties, such as an insufficient number of PIS), may account for most gene and species-tree discordance. We find overwhelming signal for Toxicofera, and also show that none of the loci included in this study supports Scleroglossa or Macrostomata. We comment on the origins and diversification of Squamata throughout the Mesozoic and underscore remaining uncertainties that persist in both deeper parts of the tree (e.g., relationships between Dibamia, Gekkota, and remaining squamates; among the three toxicoferan clades Iguania, Serpentes, and Anguiformes) and within specific clades (e.g., affinities among gekkotan, pleurodont iguanians, and colubroid families).

Environmental temperatures shape thermal physiology as well as diversification and genome-wide substitution rates in lizards

Climatic conditions changing over time and space shape the evolution of organisms at multiple levels, including temperate lizards in the family Lacertidae. Here we reconstruct a dated phylogenetic tree of 262 lacertid species based on a supermatrix relying on novel phylogenomic datasets and fossil calibrations. Diversification of lacertids was accompanied by an increasing disparity among occupied bioclimatic niches, especially in the last 10 Ma, during a period of progressive global cooling. Temperate species also underwent a genome-wide slowdown in molecular substitution rates compared to tropical and desert-adapted lacertids. Evaporative water loss and preferred temperature are correlated with bioclimatic parameters, indicating physiological adaptations to climate. Tropical, but also some populations of cool-adapted species experience maximum temperatures close to their preferred temperatures. We hypothesize these species-specific physiological preferences may constitute a handicap to prevail under rapid global warming, and contribute to explaining local lizard extinctions in cool and humid climates.

A conservation checklist of the amphibians and reptiles of Sonora, Mexico, with updated species lists

Sonora has a rich natural diversity, including reptiles and amphibians. Sonora's location on the United States-Mexico border creates some unique conservation challenges for its wildlife. We compiled a list of the amphibian and reptile species currently known for Sonora, summarized the conservation status of these species, and compared our list of species with known species lists for adjacent states. The herpetofauna of Sonora comprises 200 species of amphibians and reptiles (38 amphibians and 162 reptiles). Overall, Sonora shares the most species with Chihuahua, Sinaloa, and Arizona. Approximately 11% of the amphibian and reptile species are IUCN listed, but 35.5% are placed in a protected category by SEMARNAT, and 32.6% are categorized as high risk by the Environmental Vulnerability Score.

Tracking down the lizards from Gravenhorst's collection at the University of Wrocław: type specimens of Callopistes maculatus Gravenhorst, 1838 and three Liolaemus species rediscovered

Johann Ludwig Christian Gravenhorst's herpetological collection at the Museum of Natural History, University of Wrocław included numerous important specimens of amphibians and reptiles. The majority, if not the entirety, of this collection has long been thought to be lost. However, we were able to rediscover some type specimens of lizards. The rediscovered specimens include the holotypes of Liolaemus conspersus and L. hieroglyphicus, one syntype of Callopistes maculatus (here designated as the lectotype) and two syntypes of L. lineatus (one of which is herein designated as the lectotype). Reexamination of these specimens indicates that previous synonymies proposed for L. conspersus and two syntypes of L. hieroglyphicus are problematic; furthermore, more complex taxonomic work is needed to resolve this issue. Two rediscovered syntypes of L. lineatus differ in several scalation traits and are possibly not conspecific. The type specimens of several other species of lizards from Gravenhorst's collection (Liolaemus marmoratus, L. unicolor and two other syntypes of L. lineatus, Leiocephalus schreibersii and Chalcides viridanus) were not found and are probably lost.

Naturally clonal vertebrates are an untapped resource in ecology and evolution research

Science requires replication. The development of many cloned or isogenic model organisms is a testament to this. But researchers are reluctant to use these traditional animal model systems for certain questions in evolution or ecology research, because of concerns over relevance or inbreeding. It has largely been overlooked that there are a substantial number of vertebrate species that reproduce clonally in nature. Here we highlight how use of these naturally evolved, phenotypically complex animals can push the boundaries of traditional experimental design and contribute to answering fundamental questions in the fields of ecology and evolution.

A longitudinal study of attitudes toward evolution among undergraduates who are members of the Church of Jesus Christ of Latter-day Saints

Polling data reveal a decades-long residual rejection of evolution in the United States, based on perceived religious conflict. Similarly, a strong creationist movement has been documented internationally, including in the Muslim world. Members of the Church of Jesus Christ of Latter-day Saints (LDS, Mormon), a generally conservative denomination, have historically harbored strong anti-evolution sentiments. We report here a significant shift toward acceptance, compared to attitudes 30 years earlier, by students at Brigham Young University, which is owned and operated by the LDS church. This change appears to have multiple explanations. Students currently entering the university have been exposed to a much-improved introduction to evolution during high school. More importantly, there has been a significant decrease in negative messaging from Church authorities and in its religious education system. There is also evidence that current students have been positively influenced toward evolution by their parents, a large percentage of whom were BYU students, who earlier were given a strong science education deemed compatible with the maintenance of religious belief. A pre-post comparison demonstrates that a majority of current students become knowledgeable and accepting following a course experience focused on evolutionary principles delivered in a faith-friendly atmosphere. Elements of that classroom pedagogy, intended to promote reconciliation, are presented. Our experience may serve as a case-study for prompting changes in acceptance of evolution in other conservative religious groups.

The mitochondrial genomes of five frog species of the Neotropical genus Ischnocnema (Anura: Brachycephaloidea: Brachycephalidae)

We report the mitogenomes for five species of the Ischnocnema guentheri series, being the first described for this genus of brachycephalid frogs. We assembled mitogenomes from anchored hybrid enrichment data and recovered the 13 protein-coding genes, 22 tRNA genes, and two rRNA genes for all species. The general structure agrees with most previously sequenced neobatrachians, with two exceptions: the origin of replication of L-strand (OL) was found between tRNA-A and tRNA-N, and the position of tRNA-L and tRNA-T, which are dispersed in the control region. We provide a phylogenetic tree with outgroups, which is consistent with previous phylogenetic hypotheses.

Anchored phylogenomics illuminates the skipper butterfly tree of life

BACKGROUND

Butterflies (Papilionoidea) are perhaps the most charismatic insect lineage, yet phylogenetic relationships among them remain incompletely studied and controversial. This is especially true for skippers (Hesperiidae), one of the most species-rich and poorly studied butterfly families.

METHODS

To infer a robust phylogenomic hypothesis for Hesperiidae, we sequenced nearly 400 loci using Anchored Hybrid Enrichment and sampled all tribes and more than 120 genera of skippers. Molecular datasets were analyzed using maximum-likelihood, parsimony and coalescent multi-species phylogenetic methods.

RESULTS

All analyses converged on a novel, robust phylogenetic hypothesis for skippers. Different optimality criteria and methodologies recovered almost identical phylogenetic trees with strong nodal support at nearly all nodes and all taxonomic levels. Our results support Coeliadinae as the sister group to the remaining skippers, the monotypic Euschemoninae as the sister group to all other subfamilies but Coeliadinae, and the monophyly of Eudaminae plus Pyrginae. Within Pyrginae, Celaenorrhinini and Tagiadini are sister groups, the Neotropical firetips, Pyrrhopygini, are sister to all other tribes but Celaenorrhinini and Tagiadini. Achlyodini is recovered as the sister group to Carcharodini, and Erynnini as sister group to Pyrgini. Within the grass skippers (Hesperiinae), there is strong support for the monophyly of Aeromachini plus remaining Hesperiinae. The giant skippers (Agathymus and Megathymus) once classified as a subfamily, are recovered as monophyletic with strong support, but are deeply nested within Hesperiinae.

CONCLUSIONS

Anchored Hybrid Enrichment sequencing resulted in a large amount of data that built the foundation for a new, robust evolutionary tree of skippers. The newly inferred phylogenetic tree resolves long-standing systematic issues and changes our understanding of the skipper tree of life. These resultsenhance understanding of the evolution of one of the most species-rich butterfly families.

Amphibians and reptiles of the state of Durango, Mexico, with comparisons with adjoining states

A summary of the species of amphibians and reptiles of Durango, as well as their geographic distributions, habitat, and conservation status have been compiled. The herpetofauna of Durango consists of 36 species of amphibians and 120 species of reptiles. Durango shares the most species with Chihuahua (74.0%), and shares fewer species with Sinaloa (48.0%), Nayarit (48.7%), and Coahuila (48.0%). Arid-semiarid and Sierras habitat types have the most species, with valleys and Quebradas habitat types having fewer species. In Durango, there are several taxa of particular conservation concern including eleutherodactylid frogs, eublepharid, iguanid, phrynosomatid, and xantusid lizards, boid, colubrid, and natricid snakes, and emydid and testudinid turtles.

Climatic suitability, isolation by distance and river resistance explain genetic variation in a Brazilian whiptail lizard

Spatial patterns of genetic variation can help understand how environmental factors either permit or restrict gene flow and create opportunities for regional adaptations. Organisms from harsh environments such as the Brazilian semiarid Caatinga biome may reveal how severe climate conditions may affect patterns of genetic variation. Herein we combine information from mitochondrial DNA with physical and environmental features to study the association between different aspects of the Caatinga landscape and spatial genetic variation in the whiptail lizard Ameivula ocellifera. We investigated which of the climatic, environmental, geographical and/or historical components best predict: (1) the spatial distribution of genetic diversity, and (2) the genetic differentiation among populations. We found that genetic variation in A. ocellifera has been influenced mainly by temperature variability, which modulates connectivity among populations. Past climate conditions were important for shaping current genetic diversity, suggesting a time lag in genetic responses. Population structure in A. ocellifera was best explained by both isolation by distance and isolation by resistance (main rivers). Our findings indicate that both physical and climatic features are important for explaining the observed patterns of genetic variation across the xeric Caatinga biome.

The biomechanical role of the chondrocranium and sutures in a lizard cranium

The role of soft tissues in skull biomechanics remains poorly understood. Not least, the chondrocranium, the portion of the braincase which persists as cartilage with varying degrees of mineralization. It also remains commonplace to overlook the biomechanical role of sutures despite evidence that they alter strain distribution. Here, we examine the role of both the sutures and the chondrocranium in the South American tegu lizard Salvator merianae We use multi-body dynamics analysis (MDA) to provide realistic loading conditions for anterior and posterior unilateral biting and a detailed finite element model to examine strain magnitude and distribution. We find that strains within the chondrocranium are greatest during anterior biting and are primarily tensile; also that strain within the cranium is not greatly reduced by the presence of the chondrocranium unless it is given the same material properties as bone. This result contradicts previous suggestions that the anterior portion (the nasal septum) acts as a supporting structure. Inclusion of sutures to the cranium model not only increases overall strain magnitudes but also leads to a more complex distribution of tension and compression rather than that of a beam under sagittal bending.

Genomic timetree and historical biogeography of Caribbean island ameiva lizards (Pholidoscelis: Teiidae)

The phylogenetic relationships and biogeographic history of Caribbean island ameivas (Pholidoscelis) are not well-known because of incomplete sampling, conflicting datasets, and poor support for many clades. Here, we use phylogenomic and mitochondrial DNA datasets to reconstruct a well-supported phylogeny and assess historical colonization patterns in the group. We obtained sequence data from 316 nuclear loci and one mitochondrial marker for 16 of 19 extant species of the Caribbean endemic genus Pholidoscelis. Phylogenetic analyses were carried out using both concatenation and species tree approaches. To estimate divergence times, we used fossil teiids to calibrate a timetree which was used to elucidate the historical biogeography of these lizards. All phylogenetic analyses recovered four well-supported species groups (clades) recognized previously and supported novel relationships of those groups, including a (P. auberi + P. lineolatus) clade (western + central Caribbean), and a (P. exsul + P. plei) clade (eastern Caribbean). Divergence between Pholidoscelis and its sister clade was estimated to have occurred ~25 Ma, with subsequent diversification on Caribbean islands occurring over the last 11 Myr. Of the six models compared in the biogeographic analyses, the scenario which considered the distance among islands and allowed dispersal in all directions best fit the data. These reconstructions suggest that the ancestor of this group colonized either Hispaniola or Puerto Rico from Middle America. We provide a well-supported phylogeny of Pholidoscelis with novel relationships not reported in previous studies that were based on significantly smaller datasets. We propose that Pholidoscelis colonized the eastern Greater Antilles from Middle America based on our biogeographic analysis, phylogeny, and divergence time estimates. The closing of the Central American Seaway and subsequent formation of the modern Atlantic meridional overturning circulation may have promoted dispersal in this group.

The mitochondrial genomes of Atlas Geckos (Quedenfeldtia): mitogenome assembly from transcriptomes and anchored hybrid enrichment datasets

The nearly complete mitogenomes of the two species of North African Atlas geckos, Quedenfeldtia moerens and Q. trachyblepharus were assembled from anchored hybrid enrichment data and RNAseq data. Congruent assemblies were obtained for four samples included in both datasets. We recovered the 13 protein-coding genes, 22 tRNA genes, and two rRNA genes for both species, including partial control region. The order of genes agrees with that of other geckos.