Xantusiid "night" lizards: a puzzling phylogenetic problem revisited using likelihood-based Bayesian methods on mtDNA sequences

Historical Environment Is Reflected in Modern Population Genetics and Biogeography of an Island Endemic Lizard (Xantusia riversiana reticulata)

The restricted distribution and isolation of island endemics often produces unique genetic and phenotypic diversity of conservation interest to management agencies. However, these isolated species, especially those with sensitive life history traits, are at high risk for the adverse effects of genetic drift and habitat degradation by non-native wildlife. Here, we study the population genetic diversity, structure, and stability of a classic “island giant” (Xantusia riversiana, the Island Night Lizard) on San Clemente Island, California following the removal of feral goats. Using DNA microsatellites, we found that this population is reasonably genetically robust despite historical grazing, with similar effective population sizes and genetic diversity metrics across all sampling locations irrespective of habitat type and degree of degradation. However, we also found strong site-specific patterns of genetic variation and low genetic diversity compared to mainland congeners, warranting continued special management as an island endemic. We identify both high and low elevation areas that remain valuable repositories of genetic diversity and provide a case study for other low-dispersal coastal organisms in the face of future climate change.

Population genetics of the São Tomé caecilian (Gymnophiona: Dermophiidae: Schistometopum thomense) reveals strong geographic structuring

Islands provide exciting opportunities for exploring ecological and evolutionary mechanisms. The oceanic island of São Tomé in the Gulf of Guinea exhibits high diversity of fauna including the endemic caecilian amphibian, Schistometopum thomense. Variation in pigmentation, morphology and size of this taxon over its c. 45 km island range is extreme, motivating a number of taxonomic, ecological, and evolutionary hypotheses to explain the observed diversity. We conducted a population genetic study of S. thomense using partial sequences of two mitochondrial DNA genes (ND4 and 16S), together with morphological examination, to address competing hypotheses of taxonomic or clinal variation. Using Bayesian phylogenetic analysis and Spatial Analysis of Molecular Variance, we found evidence of four geographic clades, whose range and approximated age (c. 253 Kya – 27 Kya) are consistent with the spread and age of recent volcanic flows. These clades explained 90% of variation in ND4 (φ_CT = 0.892), and diverged by 4.3% minimum pairwise distance at the deepest node. Most notably, using Mismatch Distributions and Mantel Tests, we identified a zone of population admixture that dissected the island. In the northern clade, we found evidence of recent population expansion (Fu's F_s = −13.08 and Tajima's D = −1.80) and limited dispersal (Mantel correlation coefficient = 0.36, p = 0.01). Color assignment to clades was not absolute. Paired with multinomial regression of chromatic data, our analyses suggested that the genetic groups and a latitudinal gradient together describe variation in color of S. thomense. We propose that volcanism and limited dispersal ability are the likely proximal causes of the observed genetic structure. This is the first population genetic study of any caecilian and demonstrates that these animals have deep genetic divisions over very small areas in accordance with previous speculations of low dispersal abilities.

A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes

BACKGROUND

The extant squamates (>9400 known species of lizards and snakes) are one of the most diverse and conspicuous radiations of terrestrial vertebrates, but no studies have attempted to reconstruct a phylogeny for the group with large-scale taxon sampling. Such an estimate is invaluable for comparative evolutionary studies, and to address their classification. Here, we present the first large-scale phylogenetic estimate for Squamata.

RESULTS

The estimated phylogeny contains 4161 species, representing all currently recognized families and subfamilies. The analysis is based on up to 12896 base pairs of sequence data per species (average = 2497 bp) from 12 genes, including seven nuclear loci (BDNF, c-mos, NT3, PDC, R35, RAG-1, and RAG-2), and five mitochondrial genes (12S, 16S, cytochrome b, ND2, and ND4). The tree provides important confirmation for recent estimates of higher-level squamate phylogeny based on molecular data (but with more limited taxon sampling), estimates that are very different from previous morphology-based hypotheses. The tree also includes many relationships that differ from previous molecular estimates and many that differ from traditional taxonomy.

CONCLUSIONS

We present a new large-scale phylogeny of squamate reptiles that should be a valuable resource for future comparative studies. We also present a revised classification of squamates at the family and subfamily level to bring the taxonomy more in line with the new phylogenetic hypothesis. This classification includes new, resurrected, and modified subfamilies within gymnophthalmid and scincid lizards, and boid, colubrid, and lamprophiid snakes.

Rampant horizontal transfer of SPIN transposons in squamate reptiles

Transposable elements (TEs) are highly abundant in the genome and capable of mobility, two properties that make them particularly prone to transfer horizontally between organisms. Although the impact of horizontal transfer (HT) of TEs is well recognized in prokaryotes, the frequency of this phenomenon and its contribution to genome evolution in eukaryotes remain poorly appreciated. Here, we provide evidence that a DNA transposon called SPIN has colonized the genome of 17 species of reptiles representing nearly every major lineage of squamates, including 14 families of lizards, snakes, and amphisbaenians. Slot blot analyses indicate that SPIN has amplified to high copy numbers in most of these species, ranging from 2,000-28,000 copies per haploid genome. In contrast, we could not detect the presence of SPIN in any of the turtles (seven species from seven families) and crocodiles (four species) examined. Genetic distances between SPIN sequences from species belonging to different squamate families are consistently very low (average = 0.1), considering the deep evolutionary divergence of the families investigated (most are >100 My diverged). Furthermore, these distances fall below interfamilial distances calculated for two genes known to have evolved under strong functional constraint in vertebrates (RAG1, average = 0.24 and C-mos, average = 0.27). These data, combined with phylogenetic analyses, indicate that the widespread distribution of SPIN among squamates is the result of at least 13 independent events of HTs. Molecular dating and paleobiogeographical data suggest that these transfers took place during the last 50 My on at least three different continents (North America, South America and, Africa). Together, these results triple the number of known SPIN transfer events among tetrapods, provide evidence for a previously hypothesized transoceanic movement of SPIN transposons during the Cenozoic, and further underscore the role of HT in the evolution of vertebrate genomes.

Lizard auditory papillae: an evolutionary kaleidoscope

The evolutionary processes that modified the structure and function of lizard auditory papillae during the separation of the familial lineages during the Jurassic have resulted in a remarkable variety of family-typical papillae. These papillae vary structurally in their size, in the patterns of the distribution of hair-cell types, in the presence or absence of sub-papillae and in the configurations of the tectorial membranes. Functional differences, however, are much smaller than the structural variations might lead one to expect. To some extent, differences in innervation patterns and tectorial configurations compensate for 10-fold differences in papillar length. Nonetheless, although lizards with tiny papillae are able to maintain frequency-selective and relatively sensitive hearing, the best selectivity and most sensitive hearing is found in the largest and most complex papillae. Fundamental considerations of the tonotopic organisation of papillae leads to a likely scheme mapping the evolution of the hearing organs found in modern lizard families.

The relationship between the lizard eye and associated bony features: a cautionary note for interpreting fossil activity patterns

Activity pattern, the time of day when an animal is active, is associated with ecology. There are two major activity patterns: diurnal (awake during the day in a photopic environment) and nocturnal (awake at night in a scotopic environment). Lizards exhibit characteristic eye shapes associated with activity pattern, with scotopic-adapted lizard eyes optimized for visual sensitivity with large corneal diameters relative to their eye axial lengths, and photopic-adapted lizards optimized for visual acuity, with larger axial lengths of the eye relative to their corneal diameters. This study: (1) quantifies the relationship between the lizard eye and its associated bony anatomy (the orbit, sclerotic ring, and associated skull widths); (2) investigates how activity pattern is reflected in that bony anatomy; and (3) determines if it is possible to reliably interpret activity pattern for a lizard that does not have the soft tissue available for study, specifically, for a fossil. Knowledge of extinct lizards' activity patterns would be useful in making paleoecological interpretations. Here, 96 scotopic- and photopic-adapted lizard species are analyzed in a phylogenetic context. Although there is a close relationship between the lepidosaur eye and associated bony anatomy, based on these data activity pattern cannot be reliably interpreted for bony-only specimens, such as a fossil, possibly because of the limited ossification of the lepidosaur skull. Caution should be exercised when utilizing lizard bony anatomy to interpret light-level adaptation, either for a fossil lizard or as part of an extant phylogenetic bracket to interpret other extinct animals with sclerotic rings, such as dinosaurs.

Comparative analysis of the size and shape of the lizard eye

Lizards occupy both scotopic (light-limited) and photopic (light-rich) environments, thereby making this clade ideal for analyses of eye morphology adaptations. This study examines how in lizards the morphology of the eye varies according to activity in these different light environments. Measurements were collected on corneal diameters and axial lengths of the eye for 239 specimens of 116 lizard species (including Sphenodon) that include both species with scotopic and photopic visual adaptations. I show that the light level available to a lizard for vision has a significant effect on eye shape and size. Scotopic lizards have eye shapes that are optimized for visual sensitivity, with larger corneal diameters relative to axial lengths. However, photopic lizards do not exhibit absolutely larger axial lengths than do scotopic lizards, and the groups have the same absolute axial lengths of the eye. Results also indicate that the light level the lizard functions under is a more significant influence on eye shape, as defined by the relationship between corneal diameter and axial length of the eye, than is phylogeny.

Multi-locus DNA sequence data reveal a history of deep cryptic vicariance and habitat-driven convergence in the desert night lizardXantusia vigilisspecies complex (Squamata: Xantusiidae)

The lizard genus Xantusia of southwestern North America has received recent attention in relation to delimiting species. Using more than 500 lizards from 156 localities, we further test hypothesized species boundaries and clarify phylogeographical patterns, particularly in regions of potential secondary contact. We sequenced the entire mitochondrial cytochrome b gene for every lizard in the study, plus a second mitochondrial DNA (mtDNA) region and two nuclear introns for subsets of the total sample. Phylogenetic analyses of the mtDNA recover a well-resolved, novel hypothesis for species in the Xantusia vigilis complex. The nuclear DNA (nDNA) data provide independent support for the recognition of X. arizonae, X. bezyi and X. wigginsi. Differences between the respective mtDNA and nDNA topologies result from either the effects of lineage sorting or ancient introgression. Nuclear data confirm the inference that some populations of X. vigilis in northwestern Arizona converged on rock-crevice-dwelling morphology and are not X. arizonae with an introgressed X. vigilis mtDNA genome. The historical independence of ancient cryptic lineages of Xantusia in southern California is also corroborated, though limited introgression is detected. Our proposed biogeographical scenario indicates that diversification of this group was driven by vicariance beginning in the late Miocene. Additionally, Pleistocene climatical changes influenced Xantusia distribution, and the now inhospitable Colorado Desert previously supported night lizard presence. The current taxonomy of the group likely underestimates species diversity within the group, and our results collectively show that while convergence on the rock-crevice-dwelling morphology is one hallmark of Xantusia evolution, morphological stasis is paradoxically another.

Partitioned Bayesian analyses, partition choice, and the phylogenetic relationships of scincid lizards

Partitioned Bayesian analyses of approximately 2.2 kb of nucleotide sequence data (mtDNA) were used to elucidate phylogenetic relationships among 30 scincid lizard genera. Few partitioned Bayesian analyses exist in the literature, resulting in a lack of methods to determine the appropriate number of and identity of partitions. Thus, a criterion, based on the Bayes factor, for selecting among competing partitioning strategies is proposed and tested. Improvements in both mean -lnL and estimated posterior probabilities were observed when specific models and parameter estimates were assumed for partitions of the total data set. This result is expected given that the 95% credible intervals of model parameter estimates for numerous partitions do not overlap and it reveals that different data partitions may evolve quite differently. We further demonstrate that how one partitions the data (by gene, codon position, etc.) is shown to be a greater concern than simply the overall number of partitions. Using the criterion of the 2 ln Bayes factor > 10, the phylogenetic analysis employing the largest number of partitions was decisively better than all other strategies. Strategies that partitioned the ND1 gene by codon position performed better than other partition strategies, regardless of the overall number of partitions. Scincidae, Acontinae, Lygosominae, east Asian and North American "Eumeces" + Neoseps; North African Eumeces, Scincus, and Scincopus, and a large group primarily from sub-Saharan Africa, Madagascar, and neighboring islands are monophyletic. Feylinia, a limbless group of previously uncertain relationships, is nested within a "scincine" clade from sub-Saharan Africa. We reject the hypothesis that the nearly limbless dibamids are derived from within the Scincidae, but cannot reject the hypothesis that they represent the sister taxon to skinks. Amphiglossus, Chalcides, the acontines Acontias and Typhlosaurus, and Scincinae are paraphyletic. The globally widespread "Eumeces" is polyphyletic and we make necessary taxonomic changes.

Relict snakes of North America and their relationships within Caenophidia, using likelihood-based Bayesian methods on mitochondrial sequences

This paper focuses on the phylogenetic relationships of eight North American caenophidian snake species (Carphophis amoena, Contia tenuis, Diadophis punctatus, Farancia abacura, Farancia erytrogramma, Heterodon nasicus, Heterodon platyrhinos, and Heterodon simus) whose phylogenetic relationships remain controversial. Past studies have referred to these "relict" North American snakes either as colubrid, or as Neotropical dipsadids and/or xenodontids. Based on mitochondrial DNA ribosomal gene sequences and a likelihood-based Bayesian analysis, our study suggests that these North American snakes are not monophyletic and are nested within a group (Dipsadoidea) that contains the Dipsadidae, Xenodontidae, and Natricidae. In addition, we use the relationships proposed here to highlight putative examples of parallel evolution of hemipenial morphology among snake clades.

Testing Species Boundaries in an Ancient Species Complex with Deep Phylogeographic History: GenusXantusia(Squamata: Xantusiidae)

Identification of species in natural populations has recently received increased attention with a number of investigators proposing rigorous methods for species delimitation. Morphologically conservative species (or species complexes) with deep phylogenetic histories (and limited gene flow) are likely to pose particular problems when attempting to delimit species, yet this is crucial to comparative studies of the geography of speciation. We apply two methods of species delimitation to an ancient group of lizards (genus Xantusia) that occur throughout southwestern North America. Mitochondrial cytochrome b and nicotinamide adenine dinucleotide dehydrogenase subunit 4 gene sequences were generated from samples taken throughout the geographic range of Xantusia. Maximum likelihood, Bayesian, and nested cladogram analyses were used to estimate relationships among haplotypes and to infer evolutionary processes. We found multiple well-supported independent lineages within Xantusia, for which there is considerable discordance with the currently recognized taxonomy. High levels of sequence divergence (21.3%) suggest that the pattern in Xantusia may predate the vicariant events usually hypothesized for the fauna of the Baja California peninsula, and the existence of deeply divergent clades (18.8%-26.9%) elsewhere in the complex indicates the occurrence of ancient sundering events whose genetic signatures were not erased by the late Wisconsin vegetation changes. We present a revised taxonomic arrangement for this genus consistent with the distinct mtDNA lineages and discuss the phylogeographic history of this genus as a model system for studies of speciation in North American deserts.