Evolutionary origins and diversification of dragon lizards in Australia's tropical savannas

Australia's monsoonal tropics are dominated by the largest and least modified savanna woodlands in the world, and they are globally significant for their high biodiversity and regional endemism. Despite this, there have been very few molecular studies of the evolutionary origins and diversification of vertebrates in this region. The semi-arboreal dragon lizards of Lophognathus and Amphibolurus are widely distributed in the savanna and dry sclerophyll woodlands of Australasia, including the monsoon tropics. We sequenced a ~1400 bp region of mitochondrial DNA and a ~1400 bp nuclear gene (RAG1) to investigate the phylogenetic relationships and phylogeographic structuring of all seven species of Lophognathus and Amphibolurus. Our analyses show that there is a higher level of species and generic diversity in the monsoon tropics than previously thought, and a full morphological review and taxonomic revision of these genera is required. Relaxed molecular clock analyses indicate that species across both genera originated in the late Miocene and early Pliocene, with significant phylogeographic structure within species. We did not find any evidence that the monsoon tropics species were a monophyletic group that had diversified within the region; instead Amphibolurus and Lophognathus represent at least three independent evolutionary colonizations of the monsoon tropics. It is probable that the origins and phylogeographic patterns of the northern Lophognathus species have evolved under the climatic influence of the Australian monsoon, rather than being either an ancient Gondwanan lineage that pre-dates the monsoon or the result of a more recent dispersal event across Wallace's Line.

Habitat use affects morphological diversification in dragon lizards

Habitat use may lead to variation in diversity among evolutionary lineages because habitats differ in the variety of ways they allow for species to make a living. Here, we show that structural habitats contribute to differential diversification of limb and body form in dragon lizards (Agamidae). Based on phylogenetic analysis and ancestral state reconstructions for 90 species, we find that multiple lineages have independently adopted each of four habitat use types: rock-dwelling, terrestriality, semi-arboreality and arboreality. Given these reconstructions, we fit models of evolution to species' morphological trait values and find that rock-dwelling and arboreality limit diversification relative to terrestriality and semi-arboreality. Models preferred by Akaike information criterion infer slower rates of size and shape evolution in lineages inferred to occupy rocks and trees, and model-averaged rate estimates are slowest for these habitat types. These results suggest that ground-dwelling facilitates ecomorphological differentiation and that use of trees or rocks impedes diversification.

The importance of phylogenetic scale in tests of Bergmann's and Rapoport's rules: lessons from a clade of South American lizards

We tested for the occurrence of Bergmann's rule, the pattern of increasing body size with latitude, and Rapoport's rule, the positive relationship between geographical range size and latitude, in 34 lineages of Liolaemus lizards that occupy arid regions of the Andean foothills. We tested the climatic-variability hypothesis (CVH) by examining the relationship between thermal tolerance breadth and distribution. Each of these analyses was performed varying the level of phylogenetic inclusiveness. Bergmann's rule and the CVH were supported, but Rapoport's rule was not. More variance in the data for Bergmann's rule and the CVH was explained using species belonging to the L. boulengeri series rather than all species, and inclusion of multiple outgroups tended to obscure these macroecological patterns. Evidence for Bergmann's rule and the predicted patterns from the CVH remained after application of phylogenetic comparative methods, indicating a greater role of ecological processes rather than phylogeny in shaping the current species distributions of these lizards.

The relationship between morphology, escape behaviour and microhabitat occupation in the lizard cladeLiolaemus(Iguanidae: Tropidurinae*: Liolaemini)

Phenotypic differences among species are known to have functional consequences that in turn allow species to use different habitats. However, the role of behaviour in this ecomorphological paradigm is not well defined. We investigated the relationship between morphology, ecology and escape behaviour among 25 species of the lizard clade Liolaemus in a phylogenetic framework. We demonstrate that the relationship between morphology and characteristics of habitat structure shows little or no association, consistent with a previous study on this group. However, a significant relationship was found between morphology and escape behaviour with the distance a lizard moved from a potential predator correlated with body width, axilla-groin length, and pelvis width. A significant relationship between escape behaviour and habitat structure occupation was found; lizards that occupied tree trunks and open ground ran longer distances from predators and were found greater distances from shelter. Behavioural strategies used by these lizards in open habitats appear to have made unnecessary the evolution of limb morphology that has occurred in other lizards from other clades that are found in open settings. Understanding differences in patterns of ecomorphological relationships among clades is an important component for studying adaptive diversification.

A molecular phylogenetic study of ecological diversification in the Australian lizard genus Ctenophorus

We present phylogenetic analyses of the lizard genus Ctenophorus using 1,639 aligned positions of mitochondrial DNA sequences containing 799 parsimony-informative characters for samples of 22 species of Ctenophorus and 12 additional Australian agamid genera. Sequences from three protein-coding genes (ND1, ND2, and COI) and eight intervening tRNA genes are examined using both parsimony and maximum-likelihood analyses. Species of Ctenophorus form a monophyletic group with Rankinia adelaidensis, which we suggest placing in Ctenophorus. Ecological differentiation among species of Ctenophorus is most evident in the kinds of habitats used for shelter. Phylogenetic analyses suggest that the ancestral condition is to use burrows for shelter, and that habits of sheltering in rocks and shrubs/hummock grasses represent separately derived conditions. Ctenophorus appears to have undergone extensive cladogenesis approximately 10-12 million years ago, with all three major ecological modes being established at that time.

Evaluating trans-tethys migration: an example using acrodont lizard phylogenetics

A phylogenetic tree for acrodont lizards (Chamaeleonidae and Agamidae) is established based on 1434 bases (1041 informative) of aligned DNA positions from a 1685-1778 base pair region of the mitochondrial genome. Sequences from three protein-coding genes (ND1, ND2, and COI) are combined with sequences from eight intervening tRNA genes for samples of 70 acrodont taxa and two outgroups. Parsimony analysis of nucleotide sequences identifies eight major clades in the Acrodonta. Most agamid lizards are placed into three distinct clades. One clade is composed of all taxa occurring in Australia and New Guinea; Physignathus cocincinus from Southeast Asia is the sister taxon to the Australia-New Guinea clade. A second clade is composed of taxa occurring from Tibet and the Indian Subcontinent east through South and East Asia. A third clade is composed of taxa occurring from Africa east through Arabia and West Asia to Tibet and the Indian Subcontinent. These three clades contain all agamid lizards except Uromastyx, Leiolepis, and Hydrosaurus, which represent three additional clades of the Agamidae. The Chamaeleonidae forms another clade weakly supported as the sister taxon to the Agamidae. All eight clades of the Acrodonta contain members occurring on land masses derived from Gondwanaland. A hypothesis of agamid lizards rafting with Gondwanan plates is examined statistically. This hypothesis suggests that the African/West Asian clade is of African or Indian origin, and the South Asian clade is either of Indian or Southeast Asian origin. The shortest tree suggests a possible African origin for the former and an Indian origin for the latter, but this result is not statistically robust. The Australia-New Guinea clade rafted with the Australia-New Guinea plate and forms the sister group to a Southeast Asian taxon that occurs on plates that broke from northern Australia-New Guinea. Other acrodont taxa are inferred to be associated with the plates of Afro-Arabia and Madagascar (Chameleonidae), India (Uromastyx), or southeast Asia (Hydrosaurus and Leiolepis). Introduction of different biotic elements to Asia by way of separate Gondwanan plates may be a major theme of Asian biogeography. Three historical events may be responsible for the sharp faunal barrier between Southeast Asia and Australia-New Guinea, known as Wallace's line: (1) primary vicariance caused by plate separations; (2) secondary contact of Southeast Asian plates with Eurasia, leading to dispersal from Eurasia into Southeast Asia, and (3) dispersal of the Indian fauna (after collision of that subcontinent) to Southeast Asia. Acrodont lizards show the first and third of these biogeographic patterns and anguid lizards exhibit the second pattern. Modern faunal diversity may be influenced primarily by historical events such as tectonic collisions and land bridge connections, which are expected to promote episodic turnover of continental faunas by introducing new faunal elements into an area. Repeated tectonic collisions may be one of the most important phenomena promoting continental biodiversity. Phylogenetics is a powerful method for investigating these processes.

Evolution and phylogenetic information content of mitochondrial genomic structural features illustrated with acrodont lizards

DNA sequences from 195 squamate reptiles indicate that mitochondrial gene order is the most reliable phylogenetic character establishing monophyly of acrodont lizards and of the snake families Boidae, Colubridae, and Viperidae. Gene order shows no evidence of evolutionary parallelisms or reversals in these taxa. Derived secondary structures of mitochondrial tRNAs also prove to be useful phylogenetic characters showing no reversals. Parallelisms for secondary structures of tRNAs are restricted to deep lineages that are separated by at least 200 million years of independent evolution. Presence of a stem-and-loop structure between the genes encoding tRNA(Asn) and tRNA(Cys), where the replication origin for light-strand synthesis is typically located in vertebrate mitochondrial genomes, is found to undergo at least three and possibly as many as seven evolutionary shifts, most likely parallel losses. This character is therefore a less desirable phylogenetic marker than the other structural changes examined. Sequencing regions that contain multiple genes, including tRNA genes, may be preferable to the common practice of obtaining single-gene fragments for phylogenetic inference because it permits observation of major structural changes in the mitochondrial genome. Such characters may occasionally provide phylogenetic information on relatively short internal branches for which base substitutional changes are expected to be relatively uninformative.

Molecular phylogenetics, tRNA evolution, and historical biogeography in anguid lizards and related taxonomic families

Phylogenetic relationships among lizards of the families Anguidae, Anniellidae, Xenosauridae, and Shinisauridae are investigated using 2001 aligned bases of mitochondrial DNA sequence from the genes encoding ND1 (subunit one of NADH dehydrogenase), tRNA(Ile), tRNA(Gln), tRNA(Met), ND2, tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), tRNA(Tyr), and COI (subunit I of cytochrome c oxidase), plus the origin for light-strand replication (O(L)) between the tRNA(Asn) and the tRNA(Cys) genes. The aligned sequences contain 1013 phylogenetically informative characters. A well-resolved phylogenetic hypothesis is obtained. Because monophyly of the family Xenosauridae (Shinisaurus and Xenosaurus) is statistically rejected, we recommend placing Shinisaurus in a separate family, the Shinisauridae. The family Anniellidae and the anguid subfamilies Gerrhonotinae and Anguinae each form monophyletic groups receiving statistical support. The Diploglossinae*, which appears monophyletic, is retained as a metataxon (denoted with an asterisk) because its monophyly is statistically neither supported nor rejected. The family Anguidae appears monophyletic in analyses of the DNA sequence data, and statistical support for its monophyly is provided by reanalysis of previously published allozymic data. Anguid lizards appear to have had a northern origin in Laurasia. Taxa currently located on Gondwanan plates arrived there by dispersal from the north in two separate events, one from the West Indies to South America and another from a Laurasian plate to Morocco. Because basal anguine lineages are located in western Eurasia and Morocco, formation of the Atlantic Ocean (late Eocene) is implicated in the separation of the Anguinae from its North American sister taxon, the Gerrhonotinae. Subsequent dispersal of anguine lizards to East Asia and North America appears to have followed the Oligocene drying of the Turgai Sea. The alternative hypothesis, that anguine lizards originated in North America and dispersed to Asia via the Bering land bridge with subsequent colonization of Europe and Morocco, requires a phylogenetic tree seven steps longer than the most parsimonious hypothesis. North African, European, and West Asian anguines were isolated from others by the rapid uplift of Tibet in the late Oligocene to Miocene. Phylogenetic analysis of evolutionary changes in the gene encoding tRNA(Cys) suggests gradual reduction of dihydrouridine (D) stems by successive deletion of bases in some lineages. This evolutionary pattern contrasts with the one observed for parallel elimination of the D-stem in mitochondrial tRNAs of eight other reptile groups, in which replication slippage produces direct repeats. An unusual, enlarged TpsiC (T) stem is inferred for tRNA(Cys) in most species.

Molecular tests of phylogenetic taxonomies: a general procedure and example using four subfamilies of the lizard family Iguanidae

A general procedure is described for examining when results of molecular phylogenetic analyses warrant formal revision of taxonomies constructed using morphological characters. We illustrate this procedure with tests of monophyly for four subfamilies in the lizard family Iguanidae using 1561 aligned base positions (838 phylogenetically informative) of mitochondrial DNA sequences, representing coding regions for eight tRNAs, ND2, and portions of ND1 and COI. Ten new sequences ranging in length from 1732 to 1751 bases are compared with 12 previously reported sequences and 67 morphological characters (54 phylogenetically informative) from the literature. New morphological character states are provided for Sator. Phylogenies derived from the molecular and combined data are in agreement but both conflict with phylogenetic inferences from the morphological data alone. Strong support is found for the monophyly of the subfamilies Crotaphytinae and Phrynosomatinae. Monophyly of the Iguaninae is weakly supported in each analysis. All analyses suggest that the Tropidurinae is not monophyletic but the hypothesis of monophyly cannot be rejected. A phylogenetic taxonomy is proposed in which the Tropidurinae* is maintained as a metataxon (denoted with an asterisk), for which monophyly has not been demonstrated. Within the Phrynosomatinae, the close relationship of Sator and Sceloporus is questioned and an alternative hypothesis in which Sator is the sister taxon to a clade comprising Petrosaurus, Sceloporus, and Urosaurus is presented. Statistical tests of monophyly provide a powerful way to evaluate support for taxonomic groupings. Use of the metataxon prevents premature taxonomic rearrangements where support is lacking.

Phylogenetic relationships among Agamid lizards of the Laudakia caucasia species group: testing hypotheses of biogeographic fragmentation and an area cladogram for the Iranian Plateau

Phylogenetic relationships within the Laudakia caucasia species group on the Iranian Plateau were investigated using 1708 aligned bases of mitochondrial DNA sequence from the genes encoding ND1 (subunit one of NADH dehydrogenase), tRNAGln, tRNAIle, tRNAMet, ND2, tRNATrp, tRNAAla, tRNAAsn, tRNACys, tRNATyr, and COI (subunit I of cytochrome c oxidase). The aligned sequences contain 207 phylogenetically informative characters. Three hypotheses for historical fragmentation of Laudakia populations on the Iranian Plateau were tested. In two hypotheses, fragmentation of populations is suggested to have proceeded along continuous mountain belts that surround the Iranian Plateau. In another hypothesis, fragmentation is suggested to have resulted from a north-south split caused by uplifting of the Zagros Mountains in the late Miocene or early Pliocene [5-10 MYBP (million years before present)]. The shortest tree suggest the later hypothesis, and statistical tests reject the other two hypothesis. The phylogenetic tree is exceptional in that every branch is well supported. Geologic history provides dates for most branches of the tree. A plot of DNA substitutions against dates from geologic history refines the date for the north-south split across the Iranian Plateau to 9 MYBP (late Miocene). The rate of evolution for this segment of mtDNA is 0.65% (0.61-0.70%) change per lineage per million years. A hypothesis of area relationships for the biota of the Iranian Plateau is generated from the phylogenetic tree.

Phylogenetic relationships of toads in the Bufo bufo species group from the eastern escarpment of the Tibetan Plateau: a case of vicariance and dispersal

Phylogenetic relationships among Tibetan populations of the Bufo bufo species group are investigated using 1063 bases of mitochondrial DNA sequence from the genes encoding ND1 (subunit one of NADH dehydrogenase), tRNA(Ile), tRNA(Gln), tRNA(Met), and ND2. The aligned sequences contain 181 phylogenetically informative characters across all taxa sampled. Two hypotheses for colonization of the Tibetan Plateau are tested. A vicariant hypothesis predicts monophyly of populations from high elevations. A dispersalist hypothesis predicts monophyly of populations in each of two river drainages (Yangtze and Yellow rivers), which requires nonmonophyly of populations from high elevations. Both hypotheses are rejected in favor of a third hypothesis that combines elements of vicariance and dispersal. The most parsimonious phylogenetic tree places the high-elevation species, B. andrewsi, as the sister taxon to the other Asian Bufo populations; these high-elevation populations are postulated to have had a vicariant origin approximately 5 million years before present. The high-elevation population recognized as B. minshanicus is nested within low-elevation populations of B. gargarizans and is suggested to have dispersed onto the Tibetan Plateau more recently.

Tandem duplication via light-strand synthesis may provide a precursor for mitochondrial genomic rearrangement

A tandem duplication of the mitochondrial tRNA(Thr) and tRNA(Pro) genes in the amphisbaenian reptile Bipes biporus is the first case reported of a tandem duplication restricted to a single pair of tRNA genes in a vertebrate mitochondrial genome. Such duplications have been predicted, however, as intermediate steps in the evolution of observed mitochondrial genomic rearrangements through errors in light-strand replication. The tandem duplication reported here is evolutionarily associated with displacement of the origin for light-strand replication from its typical location in vertebrate mitochondrial genomes and loss of the dihydrouridine stem from the tRNA(Cys) gene; these factors implicate light-strand replicational errors in the tandem duplication of genic regions. Pseudogene formation in tandemly duplicated sequences appears to be an intermediate step in genomic rearrangement. However, formation of pseudogenes in the Bipes mitochondrial genome occurs in a pattern that precludes subsequent genomic rearrangement. Functional constraints placed on cleavage of mitochondrial transcripts by tRNA genes also may prevent mitochondrial genomic rearrangement.