A molecular phylogenetic hypothesis for the Asian agamid lizard genus Phrynocephalus reveals discrete biogeographic clades implicated by plate tectonics.

Phylogenetic relationships of the agamid lizard genus Phrynocephalus are described in the context of plate tectonics. A near comprehensive taxon sampling reports three data sets: (1) mitochondrial DNA from ND1 to COI (3' end of ND1, tRNAGln, tRNAIle, tRNAMet, ND2, tRNATrp, tRNAAla, tRNAAsn, tRNACys, tRNATyr, and the 5' end of COI) with 1761 aligned positional sites (1595 included, 839 informative), (2) nuclear RAG-1 DNA with 2760 aligned positional sites (342 informative), and (3) 25 informative allozyme loci with 213 alleles (107 informative when coded as presence/absence). It is hypothesized that Phrynocephalus phyletic patterns and speciation reflect fault lines of ancient plates now in Asia rejuvenated by the more recent Indian and Arabian plate collisions. Molecular estimates of lineage splits are highly congruent with geologic dates from the literature.  A southern origin for the genus in Southwest Asia is resolved in phylogenetic estimates and a northern origin is statistically rejected. On the basis of monophyly and molecular evidence several taxa previously recognized as subspecies are recognized as species: P. hongyuanensis, P. sogdianus, and P. strauchi as "Current Status"; Phrynocephalus bannikovi, Phrynocephalus longicaudatus, Phrynocephalus turcomanus, and Phrynocephalus vindumi are formally "New Status". Phylogenetic evaluation indicates a soft substrate habitat of sand for the shared ancestor of modern Phrynocephalus. Size diversity maximally overlaps in the Caspian Basin and northwestern Iranian Plateau. The greatest species numbers of six in sympatry and regional allopatry are found in the southern Caspian Basin and southern Helmand Basin, both from numerous phylogenetic lineages in close proximity attributed to tectonic induced events.

The Eurasian invasion: phylogenomic data reveal multiple Southeast Asian origins for Indian Dragon Lizards.

BACKGROUND: The Indian Tectonic Plate split from Gondwanaland approximately 120 MYA and set the Indian subcontinent on a ~ 100 million year collision course with Eurasia. Many phylogenetic studies have demonstrated the Indian subcontinent brought with it an array of endemic faunas that evolved in situ during its journey, suggesting this isolated subcontinent served as a source of biodiversity subsequent to its collision with Eurasia. However, recent molecular studies suggest that Eurasia may have served as the faunal source for some of India's biodiversity, colonizing the subcontinent through land bridges between India and Eurasia during the early to middle Eocene (~35-40 MYA). In this study we investigate whether the Draconinae subfamily of the lizard family Agamidae is of Eurasian or Indian origin, using a multi locus Sanger dataset and a novel dataset of 4536 ultraconserved nuclear element loci. RESULTS: Results from our phylogenetic and biogeographic analyses revealed support for two independent colonizations of India from Eurasian ancestors during the early to late Eocene prior to the subcontinent's hard collision with Eurasia. CONCLUSION: These results are consistent with other faunal groups and new geologic models that suggest ephemeral Eocene land bridges may have allowed for dispersal and exchange of floras and faunas between India and Eurasia during the Eocene.

Effect of body mass and melanism on heat balance in Liolaemus lizards of the goetschi clade.

The body temperature of ectotherms depends on the environmental temperatures and behavioral adjustments, but morphology may also have an effect. For example, in colder environments, animals tend to be larger and to show higher thermal inertia, as proposed by Bergmann's rule and the heat balance hypothesis (HBH). Additionally, dark coloration increases solar radiation absorption and should accelerate heat gain (thermal melanism hypothesis, TMH). We tested Bergmann's rule, the HBH and the TMH within the ITALIC! Liolaemus goetschilizard clade, which shows variability in body size and melanic coloration. We measured heating and cooling rates of live and euthanized animals, and tested how morphology and color affect these rates. Live organisms show less variable and faster heating rates compared with cooling rates, suggesting behavioral and/or physiological adjustments. Our results support Bergmann's rule and the HBH, as larger species show slower heating and cooling rates. However, we did not find a clear pattern to support the TMH. The influence of dorsal melanism on heating by radiation was masked by the body size effect in live animals, and results from euthanized individuals also showed no clear effects of melanism on heating rates. Comparison among three groups of live individuals with different degrees of melanism did not clarify the influence of melanism on heating rates. However, when euthanized animals from the same three groups were compared, we observed that darker euthanized animals actually heat faster than lighter ones, favoring the TMH. Although unresolved aspects remain, body size and coloration influenced heat exchange, suggesting complex thermoregulatory strategies in these lizards, probably regulated through physiology and behavior, which may allow these small lizards to inhabit harsh weather environments.

How Should Genes and Taxa be Sampled for Phylogenomic Analyses with Missing Data? An Empirical Study in Iguanian Lizards.

Targeted sequence capture is becoming a widespread tool for generating large phylogenomic data sets to address difficult phylogenetic problems. However, this methodology often generates data sets in which increasing the number of taxa and loci increases amounts of missing data. Thus, a fundamental (but still unresolved) question is whether sampling should be designed to maximize sampling of taxa or genes, or to minimize the inclusion of missing data cells. Here, we explore this question for an ancient, rapid radiation of lizards, the pleurodont iguanians. Pleurodonts include many well-known clades (e.g., anoles, basilisks, iguanas, and spiny lizards) but relationships among families have proven difficult to resolve strongly and consistently using traditional sequencing approaches. We generated up to 4921 ultraconserved elements with sampling strategies including 16, 29, and 44 taxa, from 1179 to approximately 2.4 million characters per matrix and approximately 30% to 60% total missing data. We then compared mean branch support for interfamilial relationships under these 15 different sampling strategies for both concatenated (maximum likelihood) and species tree (NJst) approaches (after showing that mean branch support appears to be related to accuracy). We found that both approaches had the highest support when including loci with up to 50% missing taxa (matrices with ~40-55% missing data overall). Thus, our results show that simply excluding all missing data may be highly problematic as the primary guiding principle for the inclusion or exclusion of taxa and genes. The optimal strategy was somewhat different for each approach, a pattern that has not been shown previously. For concatenated analyses, branch support was maximized when including many taxa (44) but fewer characters (1.1 million). For species-tree analyses, branch support was maximized with minimal taxon sampling (16) but many loci (4789 of 4921). We also show that the choice of these sampling strategies can be critically important for phylogenomic analyses, since some strategies lead to demonstrably incorrect inferences (using the same method) that have strong statistical support. Our preferred estimate provides strong support for most interfamilial relationships in this important but phylogenetically challenging group.

Thermal sensitivity of cold climate lizards and the importance of distributional ranges.

One of the fundamental goals in macroecology is to understand the relationship among species' geographic ranges, ecophysiology, and climate; however, the mechanisms underlying the distributional geographic patterns observed remain unknown for most organisms. In the case of ectotherms this is particularly important because the knowledge of these interactions may provide a robust framework for predicting the potential consequences of climate change in these organisms. Here we studied the relationship of thermal sensitivity and thermal tolerance in Patagonian lizards and their geographic ranges, proposing that species with wider distributions have broader plasticity and thermal tolerance. We predicted that lizard thermal physiology is related to the thermal characteristics of the environment. We also explored the presence of trade-offs of some thermal traits and evaluated the potential effects of a predicted scenario of climate change for these species. We examined sixteen species of Liolaemini lizards from Patagonia representing species with different geographic range sizes. We obtained thermal tolerance data and performance curves for each species in laboratory trials. We found evidence supporting the idea that higher physiological plasticity allows species to achieve broader distribution ranges compared to species with restricted distributions. We also found a trade-off between broad levels of plasticity and higher optimum temperatures of performance. Finally, results from contrasting performance curves against the highest environmental temperatures that lizards may face in a future scenario (year 2080) suggest that the activity of species occurring at high latitudes may be unaffected by predicted climatic changes.

Undersampling taxa will underestimate molecular divergence dates: an example from the South american lizard clade liolaemini.

Methods for estimating divergence times from molecular data have improved dramatically over the past decade, yet there are few studies examining alternative taxon sampling effects on node age estimates. Here, I investigate the effect of undersampling species diversity on node ages of the South American lizard clade Liolaemini using several alternative subsampling strategies for both time calibrations and taxa numbers. Penalized likelihood (PL) and Bayesian molecular dating analyses were conducted on a densely sampled (202 taxa) mtDNA-based phylogenetic hypothesis of Iguanidae, including 92 Liolaemini species. Using all calibrations and penalized likelihood, clades with very low taxon sampling had node age estimates younger than clades with more complete taxon sampling. The effect of Bayesian and PL methods differed when either one or two calibrations only were used with dense taxon sampling. Bayesian node ages were always older when fewer calibrations were used, whereas PL node ages were always younger. This work reinforces two important points: (1) whenever possible, authors should strongly consider adding as many taxa as possible, including numerous outgroups, prior to node age estimation to avoid considerable node age underestimation and (2) using more, critically assessed, and accurate fossil calibrations should yield improved divergence time estimates.

Chasing the Patagonian sun: comparative thermal biology of Liolaemus lizards.

The importance of the thermal environment for ectotherms and its relationship with thermal physiology and ecology is widely recognized. Several models have been proposed to explain the evolution of the thermal biology of ectotherms, but experimental studies have provided mixed support. Lizards from the Liolaemus goetschi group can be found along a wide latitudinal range across Argentina. The group is monophyletic and widely distributed, and therefore provides excellent opportunities to study the evolution of thermal biology. We studied thermal variables of 13 species of the L. goetschi group, in order to answer three questions. First, are aspects of the thermal biology of the L. goetschi group modelled by the environment or are they evolutionarily conservative? Second, have thermal characteristics of these animals co-evolved? And third, how do the patterns of co-evolution observed within the L. goetschi group compare to those in a taxonomically wider selection of species of Liolaemus? We collected data on 13 focal species and used species information of Liolaemus lizards available in the literature and additional data obtained by the authors. We tackled these questions using both conventional and phylogenetically based analyses. Our results show that lizards from the L. goetschi group and the genus Liolaemus in general vary in critical thermal minimum in relation to mean air temperature, and particularly the L. goetschi group shows that air temperature is associated with critical thermal range, as well as with body temperature. Although the effect of phylogeny cannot be ignored, our results indicate that these thermal biology aspects are modelled by cold environments of Patagonia, while other aspects (preferred body temperature and critical thermal maximum) are more conservative. We found evidence of co-evolutionary patterns between critical thermal minimum and preferred body temperature at both phylogenetic scales (the L. goetschi group and the extended sample of 68 Liolaemus species).

Evolution of extreme body size disparity in monitor lizards (Varanus).

Many features of species' biology, including life history, physiology, morphology, and ecology are tightly linked to body size. Investigation into the causes of size divergence is therefore critical to understanding the factors shaping phenotypic diversity within clades. In this study, we examined size evolution in monitor lizards (Varanus), a clade that includes the largest extant lizard species, the Komodo dragon (V. komodoensis), as well as diminutive species that are nearly four orders of magnitude smaller in adult body mass. We demonstrate that the remarkable body size disparity of this clade is a consequence of different selective demands imposed by three major habitat use patterns-arboreality, terrestriality, and rock-dwelling. We reconstructed phylogenetic relationships and ancestral habitat use and applied model selection to determine that the best-fitting evolutionary models for species' adult size are those that infer oppositely directed adaptive evolution associated with terrestriality and rock-dwelling, with terrestrial lineages evolving extremely large size and rock-dwellers becoming very small. We also show that habitat use affects the evolution of several ecologically important morphological traits independently of body size divergence. These results suggest that habitat use exerts a strong, multidimensional influence on the evolution of morphological size and shape disparity in monitor lizards.

Phylogeny of the freshwater copepod Mesocyclops (Crustacea: Cyclopidae) based on combined molecular and morphological data, with notes on biogeography.

We combined molecular and morphological characters in a copepod taxon for which obtaining a sufficiently high number of characters that evolve at different rates is a challenge. Few molecular markers are known to resolve evolutionary relationships in the copepods, and thus there is potential for morphology to contribute substantially to phylogenetic reconstruction. We used a morphology based tree of the entire Mesocyclops genus to guide our taxon sampling of 10 species for molecular and combined analyses. Morphology including polymorphic characters, 18S rDNA, and ITS2 sequences were analyzed using parsimony, ML, and Bayesian methods. Strong similarities among topologies were observed regardless of the character type or algorithm, with higher levels of support obtained in combined data analyses. In combined analyses Old World species formed a monophyletic group and New World species formed a paraphyletic group in this freshwater, predominantly (sub)tropical genus. Mesocyclops darwini was the single taxon whose relationships showed conflict among the previous reconstructions using only morphological characters and the tree inferred from the combined data set. Support for these alternative positions of M. darwini were compared using constraint tests, with the result supporting monophyly of Old World taxa.

Live birth among Iguanian lizards predates Pliocene--Pleistocene glaciations.

Among tetrapods, viviparity is estimated to have evolved independently within Squamata (lizards and snakes) more than 100 times, most frequently in species occupying cold climate environments. Because of this relationship with cold climates, it is sometimes assumed that many origins of squamate viviparity occurred over the past 2.5-4 Myr during the Pliocene-Pleistocene glaciations; however, this hypothesis is untested. Divergence-dating analysis on a 733-species tree of Iguanian lizards recovers 20 independent lineages that have evolved viviparity, of which 13 multispecies groups derived live birth prior to glacial advances (8-66 Myr ago). These results place the transitions from egg-laying to live birth among squamates in a well-supported historical context to facilitate examination of the underlying phenotypic and genetic changes associated with this complex shift in reproduction.

Phylogenetic relationships among iguanian lizards using alternative partitioning methods and TSHZ1: a new phylogenetic marker for reptiles.

We present phylogenetic hypotheses for the major iguanian lizard lineages and several squamate outgroups using a combined analysis of 4950 aligned base positions representing two intronless nuclear genes, TSHZ1 and RAG1. Bayesian analyses using reversible jump (RJ) mixture model selection are conducted and compared with a priori partitioned, mixed model maximum likelihood analyses. Bayesian credibility values and ML bootstraps are comparable with strong support at deep nodes and within acrodonts, but weak support for the twelve iguanid lineages. Accounting for pattern and rate heterogeneity is becoming commonplace and is essential for accurate phylogeny reconstruction.

Phylogenetic relationships and heterogeneous evolutionary processes among phrynosomatine sand lizards (Squamata, Iguanidae) revisited.

Phylogenetic analyses of DNA sequences were conducted to evaluate four alternative hypotheses of phrynosomatine sand lizard relationships. Sequences comprising 2871 aligned base pair positions representing the regions spanning ND1-COI and cyt b-tRNA(Thr) of the mitochondrial genome from all recognized sand lizard species were analyzed using unpartitioned parsimony and likelihood methods, likelihood methods with assumed partitions, Bayesian methods with assumed partitions, and Bayesian mixture models. The topology (Uma, (Callisaurus, (Cophosaurus, Holbrookia))) and thus monophyly of the "earless" taxa, Cophosaurus and Holbrookia, is supported by all analyses. Previously proposed topologies in which Uma and Callisaurus are sister taxa and those in which Holbrookia is the sister group to all other sand lizard taxa are rejected using both parsimony and likelihood-based significance tests with the combined, unparitioned data set. Bayesian hypothesis tests also reject those topologies using six assumed partitioning strategies, and the two partitioning strategies presumably associated with the most powerful tests also reject a third previously proposed topology, in which Callisaurus and Cophosaurus are sister taxa. For both maximum likelihood and Bayesian methods with assumed partitions, those partitions defined by codon position and tRNA stem and nonstems explained the data better than other strategies examined. Bayes factor estimates comparing results of assumed partitions versus mixture models suggest that mixture models perform better than assumed partitions when the latter were not based on functional characteristics of the data, such as codon position and tRNA stem and nonstems. However, assumed partitions performed better than mixture models when functional differences were incorporated. We reiterate the importance of accounting for heterogeneous evolutionary processes in the analysis of complex data sets and emphasize the importance of implementing mixed model likelihood methods.

A phylogenetic test for adaptive convergence in rock-dwelling lizards.

Phenotypic similarity of species occupying similar habitats has long been taken as strong evidence of adaptation, but this approach implicitly assumes that similarity is evolutionarily derived. However, even derived similarities may not represent convergent adaptation if the similarities did not evolve as a result of the same selection pressures; an alternative possibility is that the similar features evolved for different reasons, but subsequently allowed the species to occupy the same habitat, in which case the convergent evolution of the same feature by species occupying similar habitats would be the result of exaptation. Many lizard lineages have evolved to occupy vertical rock surfaces, a habitat that places strong functional and ecological demands on lizards. We examined four clades in which species that use vertical rock surfaces exhibit long hindlimbs and flattened bodies. Morphological change on the phylogenetic branches leading to the rock-dwelling species in the four clades differed from change on other branches of the phylogeny; evolutionary transitions to rock-dwelling generally were associated with increases in limb length and decreases in head depth. Examination of particular characters revealed several different patterns of evolutionary change. Rock-dwelling lizards exhibited similarities in head depth as a result of both adaptation and exaptation. Moreover, even though rock-dwelling species generally had longer limbs than their close relatives, clade-level differences in limb length led to an overall lack of difference between rock- and non-rock-dwelling lizards. These results indicate that evolutionary change in the same direction in independent lineages does not necessarily produce convergence, and that the existence of similar advantageous structures among species independently occupying the same environment may not indicate adaptation.

Assembly of the eastern North American herpetofauna: new evidence from lizards and frogs.

Darwin first recognized the importance of episodic intercontinental dispersal in the establishment of worldwide biotic diversity. Faunal exchange across the Bering Land Bridge is a major example of such dispersal. Here, we demonstrate with mitochondrial DNA evidence that three independent dispersal events from Asia to North America are the source for almost all lizard taxa found in continental eastern North America. Two other dispersal events across Beringia account for observed diversity among North American ranid frogs, one of the most species-rich groups of frogs in eastern North America. The contribution of faunal elements from Asia via dispersal across Beringia is a dominant theme in the historical assembly of the eastern North American herpetofauna.

Phylogenetic relationships of South American lizards of the genus Stenocercus (Squamata: Iguania): A new approach using a general mixture model for gene sequence data.

The South American iguanian lizard genus Stenocercus includes 54 species occurring mostly in the Andes and adjacent lowland areas from northern Venezuela and Colombia to central Argentina at elevations of 0-4000m. Small taxon or character sampling has characterized all phylogenetic analyses of Stenocercus, which has long been recognized as sister taxon to the Tropidurus Group. In this study, we use mtDNA sequence data to perform phylogenetic analyses that include 32 species of Stenocercus and 12 outgroup taxa. Monophyly of this genus is strongly supported by maximum parsimony and Bayesian analyses. Evolutionary relationships within Stenocercus are further analyzed with a Bayesian implementation of a general mixture model, which accommodates variability in the pattern of evolution across sites. These analyses indicate a basal split of Stenocercus into two clades, one of which receives very strong statistical support. In addition, we test previous hypotheses using non-parametric and parametric statistical methods, and provide a phylogenetic classification for Stenocercus.

Tempo and mode of evolutionary radiation in iguanian lizards.

Identification of general properties of evolutionary radiations has been hindered by the lack of a general statistical and phylogenetic approach applicable across diverse taxa. We present a comparative analytical framework for examining phylogenetic patterns of diversification and morphological disparity with data from four iguanian-lizard taxa that exhibit substantially different patterns of evolution. Taxa whose diversification occurred disproportionately early in their evolutionary history partition more of their morphological disparity among, rather than within, subclades. This inverse relationship between timing of diversification and morphological disparity within subclades may be a general feature that transcends the historically contingent properties of different evolutionary radiations.

Molecular phylogenetic evidence for ancient divergence of lizard taxa on either side of Wallace's Line.

Wallace's Line, separating the terrestrial faunas of South East Asia from the Australia-New Guinea region, is the most prominent and well-studied biogeographical division in the world. Phylogenetically distinct subgroups of major animal and plant groups have been documented on either side of Wallace's Line since it was first proposed in 1859. Despite its importance, the temporal history of fragmentation across this line is virtually unknown and the geological foundation has rarely been discussed. Using molecular phylogenetics and dating techniques, we show that the split between taxa in the South East Asian and the Australian-New Guinean geological regions occurred during the Late Jurassic to Early Cretaceous in two independent lizard clades. This estimate is compatible with the hypothesis of rifting Gondwanan continental fragments during the Mesozoic and strongly rejects the hypothetical origin of various members of the Australian-New Guinean herpetofauna as relatively recent invasions from South East Asia. Our finding suggests an ancient fragmentation of lizard taxa on either side of Wallace's Line and provides further evidence that the composition of modern global communities has been significantly affected by rifting and accretion of Gondwanan continental plates during the Middle to Late Mesozoic.

Rostral horn evolution among agamid lizards of the genus Ceratophora endemic to Sri Lanka.

The first phylogenetic hypothesis for the Sri Lankan agamid lizard genus Ceratophora is presented based on 1670 aligned base positions (472 parsimony informative) of mitochondrial DNA sequences, representing coding regions for eight tRNAs, ND2, and portions of ND1 and COI. Phylogenetic analysis reveals multiple origins and possibly losses of rostral horns in the evolutionary history of Ceratophora. Our data suggest a middle Miocene origin of Ceratophora with the most recent branching of recognized species occurring at the Pliocene/Pleistocene boundary. Haplotype divergence suggests that an outgroup species, Lyriocephalus scutatus, dates at least to the Pliocene. These phylogenetic results provide a framework for comparative studies of the behavioral ecological importance of horn evolution in this group.