Phylogenetic relationships of phrynosomatid lizards based on nuclear and mitochondrial data,and a revised phylogeny for Sceloporus

Phrynosomatid lizards are among the most common and diverse groups of reptiles in western North America, Mexico, and Central America. Phrynosomatidae includes 136 species in 10 genera. Phrynosomatids are used as model systems in many research programs in evolution and ecology, and much of this research has been undertaken in a comparative phylogenetic framework. However, relationships among many phrynosomatid genera are poorly supported and in conflict between recent studies. Further, previous studies based on mitochondrial DNA sequences suggested that the most species-rich genus (Sceloporus) is possibly paraphyletic with respect to as many as four other genera (Petrosaurus, Sator, Urosaurus, and Uta). Here, we collect new sequence data from five nuclear genes and combine them with published data from one additional nuclear gene and five mitochondrial gene regions. We compare trees from nuclear and mitochondrial data from 37 phrynosomatid taxa, including a “species tree” (from BEST) for the nuclear data. We also present a phylogeny for 122 phrynosomatid species based on maximum likelihood analysis of the combined data, which provides a strongly-supported hypothesis for relationships among most phrynosomatid genera and includes most phrynosomatid species. Our results strongly support the monophyly of Sceloporus (including Sator) and many of the relationships within it. We present a new classification for phrynosomatid lizards and the genus Sceloporus, and offer a new tree with branch lengths for use in comparative studies.     

Delimiting species using DNA and morphological variation and discordant species limits in spiny lizards (Sceloporus)

Haplotype phylogenies based on DNA sequence data are increasingly being used to test traditional species-level taxonomies based on morphology. However, few studies have critically compared species limits based on morphological and DNA data, and the methods used to delimit species using either type of data are only rarely explained. In this paper, we review three approaches for species delimitation (tree-based with DNA data and tree-based and character-based with morphological data) and propose explicit protocols for each. We then compare species limits inferred from these approaches, using morphological and mtDNA data for the Yarrow's spiny lizard (Sceloporus jarrovii), a traditionally polytypic species from the southwestern United States and Mexico. All three approaches support division of S. jarrovii into five species, but only two species are the same among the three approaches. We find the greatest support for the five species that are delimited based on mtDNA data, and we argue that mtDNA data may have important (and previously unappreciated) advantages for species delimitation. Because different data and approaches can disagree so extensively, our results demonstrate that the methodology of species delimitation is a critical issue in systematics.     

Phylogenetic relationships of horned lizards (Phrynosoma) based on nuclear and mitochondrial data: evidence for a misleading mitochondrial gene tree

It has proven remarkably difficult to obtain a well-resolved and strongly supported phylogeny for horned lizards (Phrynosoma) because of incongruence between morphological and mitochondrial DNA sequence data. We infer the phylogenetic relationships among all 17 extant Phrynosoma species using >5.1 kb of mtDNA (12S rRNA, 16S rRNA, ND1, ND2, ND4, Cyt b, and associated tRNA genes), and >2.2kb from three nuclear genes (RAG-1, BDNF, and GAPD) for most taxa. We conduct separate and combined phylogenetic analyses of these data using maximum parsimony, maximum likelihood, and Bayesian methods. The phylogenetic relationships inferred from the mtDNA data are congruent with previous mtDNA analyses based on fewer characters and provide strong support for most branches. However, we detected strong incongruence between the mtDNA and nuclear data using comparisons of branch support and Shimodaira-Hasegawa tests, with the (P. platyrhinos+P. goodei) clade identified as the primary source of this conflict. Our analysis of a P. mcalliixP. goodei hybrid suggests that this incongruence is caused by reticulation via introgressive hybridization. Our preferred phylogeny based on an analysis of the combined data (excluding the introgressed mtDNA data) provides a new framework for interpreting character evolution and biogeography within Phrynosoma. In the context of this improved phylogeny we propose a phylogenetic taxonomy highlighting four clades: (1) Tapaja, containing the viviparous short-horned lizards P. ditmarsi, P. hernandesi, P. douglasii, and P. orbiculare; (2) Anota, containing species with prominent cranial horns (P. solare, P. mcallii, and the P. coronatum group); (3) Doliosaurus, containing three species lacking antipredator blood-squirting (P. modestum, P. platyrhinos, and P. goodei); and (4) Brevicauda, containing two viviparous species with extremely short tails that lack blood-squirting (P. braconnieri and P. taurus).     

Species trees for the tree swallows (Genus Tachycineta): An alternative phylogenetic hypothesis to the mitochondrial gene tree

The New World swallow genus Tachycineta comprises nine species that collectively have a wide geographic distribution and remarkable variation both within- and among-species in ecologically important traits. Existing phylogenetic hypotheses for Tachycineta are based on mitochondrial DNA sequences, thus they provide estimates of a single gene tree. In this study we sequenced multiple individuals from each species at 16 nuclear intron loci. We used gene concatenated approaches (Bayesian and maximum likelihood) as well as coalescent-based species tree inference to reconstruct phylogenetic relationships of the genus. We examined the concordance and conflict between the nuclear and mitochondrial trees and between concatenated and coalescent-based inferences. Our results provide an alternative phylogenetic hypothesis to the existing mitochondrial DNA estimate of phylogeny. This new hypothesis provides a more accurate framework in which to explore trait evolution and examine the evolution of the mitochondrial genome in this group.     

Prevalence of the nematode Spauligodon giganticus (Oxyurida: Pharyngodonidae) in neonatal Yarrow's spiny lizards, Sceloporus jarrovii (Sauria: Iguanidae).

Examination of the digestive tracts of 31 field-collected neonatal Sceloporus jarrovii Cope, 1875, estimated to be between 1 day and 2 wk of age revealed a prevalence of 74% (23/31) infection for eggs and/or larvae of the nematode Spauligodon giganticus (Read and Amrein, 1953) Skrjabin, Schikhobalova, and Lagodovskaja, 1960, and a prevalence of 3% (1/31) for the cestode Oochoristica sp. Twenty-two captive-born neonatal Sceloporus jarrovii, no older than 1 day revealed a prevalence of 14% for eggs of Spauligodon giganticus. Substrate licking, a behavior evident soon after birth, apparently is responsible for this almost immediate infection by S. jarrovii.     

Differences in the thermal physiology of adult Yarrow's spiny lizards (Sceloporus jarrovii) in relation to sex and body size

Sexual size dimorphism (SSD) is often assumed to reflect the phenotypic consequences of differential selection operating on each sex. Species that exhibit SSD may also show intersexual differences in other traits, including field‐active body temperatures, preferred temperatures, and locomotor performance. For these traits, differences may be correlated with differences in body size or reflect sex‐specific trait optima. Male and female Yarrow's spiny lizards, Sceloporus jarrovii, in a population in southeastern Arizona exhibit a difference in body temperature that is unrelated to variation in body size. The observed sexual variation in body temperature may reflect divergence in thermal physiology between the sexes. To test this hypothesis, we measured the preferred body temperatures of male and female lizards when recently fed and fasted. We also estimated the thermal sensitivity of stamina at seven body temperatures. Variation in these traits provided an opportunity to determine whether body size or sex‐specific variation unrelated to size shaped their thermal physiology. Female lizards, but not males, preferred a lower body temperature when fasted, and this pattern was unrelated to body size. Larger individuals exhibited greater stamina, but we detected no significant effect of sex on the shape or height of the thermal performance curves. The thermal preference of males and females in a thermal gradient exceeded the optimal temperature for performance in both sexes. Our findings suggest that differences in thermal physiology are both sex‐ and size‐based and that peak performance at low body temperatures may be adaptive given the reproductive cycles of this viviparous species. We consider the implications of our findings for the persistence of S. jarrovii and other montane ectotherms in the face of climate warming.     

Phylogenetic relationships in the spoon tarsus subgroup of Hawaiian Drosophila: conflict and concordance between gene trees

The Hawaiian Drosophilidae contains approximately 1000 species, placed in species groups and subgroups based largely on secondary sexual modifications to wings, forelegs and mouthparts. Members of the spoon tarsus subgroup possess a cup-shaped structure on the foretarsi of males. Eight of the twelve species in this subgroup are found only on the Big Island of Hawaii, suggesting that they have diverged within the past 600,000 years. This rapid diversification has made determining the relationships within this group difficult to infer. We use 13 genes, including nine rapidly evolving nuclear loci, to estimate relationships within the spoon tarsus species, as well as to test the monophyly of this subgroup. A variety of analytical approaches are used, including individual and concatenated analyses, Bayesian estimation of species trees and Bayesian untangling of concordance knots. We find widespread agreement between phylogenetic estimates derived from different methods, although some incongruence is present. Notably, our analyses suggest that the spoon tarsus subgroup, as currently defined, is not monophyletic.     

Phylogeny, divergence times and species limits of spiny lizards (Sceloporus magister species group) in western North American deserts and Baja California

The broad distribution of the Sceloporus magister species group (squamata: phrynosomatidae) throughout western North America provides an appropriate model for testing biogeographical hypotheses explaining the timing and origins of diversity across mainland deserts and the Baja California Peninsula. We inferred concordant phylogenetic trees describing the higher-level relationships within the magister group using 1.6 kb of mitochondrial DNA (mtDNA) and 1.7 kb of nuclear DNA data. These data provide strong support for the parallel divergence of lineages endemic to the Baja California Peninsula (S. zosteromus and the orcutti complex) in the form of two sequential divergence events at the base of the magister group phylogeny. A relaxed phylogenetic analysis of the mtDNA data using one fossil and one biogeographical constraint provides a chronology of these divergence events and evidence that further diversification within the Baja California clades occurred simultaneously, although patterns of geographical variation and speciation between clades differ. We resolved four major phylogeographical clades within S. magister that (i) provide a novel phylogenetic placement of the Chihuahuan Desert populations sister to the Mojave Desert; (ii) illustrate a mixed history for the Colorado Plateau that includes Mojave and Sonoran Desert components; and (iii) identify an area of overlap between the Mojave and Sonoran Desert clades near Yuma, Arizona. Estimates of bidirectional migration rates among populations of S. magister using four nuclear loci support strong asymmetries in gene flow among the major mtDNA clades. Based on the nonexclusivity of mtDNA haplotypes, nuclear gene flow among populations and wide zones of phenotypic intergradation, S. magister appears to represent a single geographically variable and widespread species.     

Locomotor performance of hatchling fence lizards (Sceloporus occidentalis): Quantitative genetics and morphometric correlates

Summary We examined heritabilities and correlations among measures of locomotor performance (speed, stamina) and among possible morphometric determinants of performance (hindlimb span, tail length) in families of hatchling lizards (Sceloporus occidentalis). We were particularly interested in determining whether these traits were heritable and thus might potentially respond genetically to selection. Moreover, we wished to determine whether speed and stamina are negatively genetically correlated, as suggested bya priori physiological and empirical considerations. All four traits appeared to be significantly heritable. Broadsense heritabilities were 0.33–0.36 for speed, 0.35–0.36 for stamina, 0.45–0.51 for hindlimb span, and 0.46–0.47 for tail length. Contrary to expectations, speed and stamina were not negatively genetically correlated. Hindlimb span and tail length, however, were negatively genetically correlated (but not phenotypically correlated). Hindlimb span and stamina were positively phenotypically correlated. Thus, for example, selection for longer hindlimb span could potentially result in shorter tails, contrary to evolutionary predictions based only on phenotypic correlations.     

Timing major conflict between mitochondrial and nuclear genes in species relationships of Polygonia butterflies (Nymphalidae: Nymphalini)

Background  

Major conflict between mitochondrial and nuclear genes in estimating species relationships is an increasingly common finding in animals. Usually this is attributed to incomplete lineage sorting, but recently the possibility has been raised that hybridization is important in generating such phylogenetic patterns. Just how widespread ancient and/or recent hybridization is in animals and how it affects estimates of species relationships is still not well-known.     

The interplay between natural and sexual selection in the evolution of sexual size dimorphism in Sceloporus lizards (Squamata: Phrynosomatidae)

Sexual size dimorphism (SSD) evolves because body size is usually related to reproductive success through different pathways in females and males. Female body size is strongly correlated with fecundity, while in males, body size is correlated with mating success. In many lizard species, males are larger than females, whereas in others, females are the larger sex, suggesting that selection on fecundity has been stronger than sexual selection on males. As placental development or egg retention requires more space within the abdominal cavity, it has been suggested that females of viviparous lizards have larger abdomens or body size than their oviparous relatives. Thus, it would be expected that females of viviparous species attain larger sizes than their oviparous relatives, generating more biased patterns of SSD. We test these predictions using lizards of the genus Sceloporus. After controlling for phylogenetic effects, our results confirm a strong relationship between female body size and fecundity, suggesting that selection for higher fecundity has had a main role in the evolution of female body size. However, oviparous and viviparous females exhibit similar sizes and allometric relationships. Even though there is a strong effect of body size on female fecundity, once phylogenetic effects are considered, we find that the slope of male on female body size is significantly larger than one, providing evidence of greater evolutionary divergence of male body size. These results suggest that the relative impact of sexual selection acting on males has been stronger than fecundity selection acting on females within Sceloporus lizards.     

Delimiting species: comparing methods for Mendelian characters using lizards of the Sceloporus grammicus (Squamata: Phrynosomatidae) complex

Species form the fundamental units of analysis in many areas of biology and, therefore, rigorous delimitation of this unit is important to a broad array of researchers. Recently, many new empirical methods have been proposed to delimit species in nature, and a large literature exists on the theoretical merit and superiority of each method. However, few empirical studies actually compare the results of these methods applied in the same study system. We used a large allozyme and chromosome dataset to apply a number of genetic-distance, character-based, and tree-based methods to a well-studied, data-rich system: the Sceloporus grammicus lizard complex of central Mexico. We hypothesized species boundaries under a general lineage or evolutionary species conceptual framework in an a priori fashion using mapped restriction-site data (mitochondrial DNA and nuclear rDNA), allozymes, and morphology. We then compared the ability of different methods to recover the "hypothesized evolutionary species" (HES). Highton's genetic-distance method and a tree-based method consistently recovered all four HES, although sometimes with weak support. With two exceptions, other methods recovered the same HES, but additional groups were weakly delimited and nested within the HES. Given the apparent recent divergence of some of the chromosome races and distinct populations in this complex, these are encouraging results. We emphasize the value of specifying testable criteria as clearly as possible and testing these with methods that make use of different properties of a single dataset.     

Molecular phylogeny of Allomyces macrogynus: Congruency between nuclear ribosomal RNA- and mitochondrial protein-based trees

We have sequenced the nuclear and mitochondrial small subunit rRNA genes (rns) and the mitochondrial genes coding for subunits 1 and 3 of the cytochrome oxidase (cox1 and cox3, respectively) of the chytridiomycete Allomyces macrogynus. Phylogenetic trees inferred from the derived COX1 and COX3 proteins and the nuclear rns sequences show with good bootstrap support that A. macrogynus is an early diverging fungus. The trees inferred from mitochondrial rns sequences do not yield a topology that is supported by bootstrap analysis. The similarity and the relative robustness of the nuclear rns and the mitochondrial protein-derived phylogenetic trees suggest that protein sequences are of higher value than rRNA sequences for reconstructing mitochondrial evolution. In addition, our trees support a monophyletic origin of mitochondria for the range of analyzed eukaryotes. Correspondence to: B. Franz Lang     

An experimental manipulation of food and water: growth and time-space utilization of hatchling lizards (Sceloporus undulatus)

Summary At a site in western Nebraska, we studied the effects of supplemental food and water on growth, survival, and habitat utilization of hatchling iguanid lizards, Sceloporus undulatus. Hatchlings were randomly assigned to one of four experimental groups and received either (1) supplemental food only, (2) supplemental water only, (3) supplemental food and water, or (4) neither supplemental food nor water (control). Neither supplemental food nor supplemental water affected growth rate in length, growth rate in mass, or size-specific body mass (i.e. robustness) at the end of the treatment period. Over the 45 day treatment pierod, there were no significant differences among treatment groups in survival, home-range size, or habitat utilization. There were, however, significant differences in the index of activity. Hatchlings in the treatment groups receiving supplemental water were sighted more often than hatchlings in the control group. We conclude that, during our study, neither food nor water limited growth but that water limited activity.     

Phylogenetic relationships of the Dactyloa clade of Anolis lizards based on nuclear and mitochondrial DNA sequence data

The Dactyloa clade, one of two major subgroups of mainland Anolis lizards, is distributed from Costa Rica to Peru, including the Amazon region and the southern Lesser Antilles. We estimated the phylogenetic relationships within Dactyloa based on mitochondrial (ND2, five transfer-RNAs, COI) and nuclear (RAG1) gene regions using likelihood and Bayesian methods under different partition strategies. In addition, we tested the monophyly of five previously recognized groups within Dactyloa. The data strongly support the monophyly of Dactyloa and five major clades: eastern, latifrons, Phenacosaurus, roquet and western, each of which exhibits a coherent geographic range. Relationships among the five major clades are less clear: support for basal nodes within Dactyloa is weak and some contradictory relationships are supported by different datasets and/or phylogenetic methods. Of the previously recognized subgroups within Dactyloa, only the roquet series consistently passed the topology tests applied. The monophyly of the aequatorialis, latifrons (as traditionally circumscribed) and punctatus series was strongly rejected, and the monophyly of Phenacosaurus (as traditionally circumscribed) yielded mixed results. The results of the phylogenetic analyses suggest the need for a revised taxonomy and have implications for the biogeography and tempo of the Dactyloa radiation.     

Phylogeny of Passerida (Aves: Passeriformes) based on nuclear and mitochondrial sequence data

Passerida is a monophyletic group of oscine passerines that includes almost 3500 species (about 36%) of all bird species in the world. The current understanding of higher-level relationships within Passerida is based on DNA-DNA hybridizations [C.G. Sibley, J.E. Ahlquist, Phylogeny and Classification of Birds, 1990, Yale University Press, New Haven, CT]. Our results are based on analyses of 3130 aligned nucleotide sequence data obtained from 48 ingroup and 13 outgroup genera. Three nuclear genes were sequenced: c-myc (498-510 bp), RAG-1 (930 bp), and myoglobin (693-722 bp), as well one mitochondrial gene; cytochrome b (879 bp). The data were analysed by parsimony, maximum-likelihood, and Bayesian inference. The African rockfowl and rockjumper are found to constitute the deepest branch within Passerida, but relationships among the other taxa are poorly resolved--only four major clades receive statistical support. One clade corresponds to Passeroidea of [C.G. Sibley, B.L. Monroe, Distribution and Taxonomy of Birds of the World, 1990, Yale University Press, New Haven, CT] and includes, e.g., flowerpeckers, sunbirds, accentors, weavers, estrilds, wagtails, finches, and sparrows. Starlings, mockingbirds, thrushes, Old World flycatchers, and dippers also group together in a clade corresponding to Muscicapoidea of Sibley and Monroe [op. cit.]. Monophyly of their Sylvioidea could not be corroborated--these taxa falls either into a clade with wrens, gnatcatchers, and nuthatches, or one with, e.g., warblers, bulbuls, babblers, and white-eyes. The tits, penduline tits, and waxwings belong to Passerida but have no close relatives among the taxa studied herein.     

Phylogeny of Thaumastodermatidae (Gastrotricha: Macrodasyida) inferred from nuclear and mitochondrial sequence data

Background

Phylogenetic relationships within Gastrotricha are poorly known. Attempts to shed light on this subject using morphological traits have led to hypotheses lacking satisfactory statistical support; it seemed therefore that a different approach was needed.

Methodology/Principal Findings

In this paper we attempt to elucidate the relationships within the taxonomically vast family Thaumastodermatidae (Macrodasyida) using molecular sequence data. The study includes representatives of all the extant genera of the family and for the first time uses a multi-gene approach to infer evolutionary liaisons within Gastrotricha. The final data set comprises sequences of three genes (18S, 28S rDNA and COI mtDNA) from 41 species, including 29 thaumastodermatids, 11 non-thaumastodermatid macrodasyidans and a single chaetonotidan. Molecular data was analyzed as a combined set of 3 genes and as individual genes, using Bayesian and maximum likelihood approaches. Two different outgroups were used: Xenotrichula intermedia (Chaetonotida) and members of the putative basal Dactylopodola (Macrodasyida). Thaumastodermatidae and all other sampled macrodasyidan families were found monophyletic except for Cephalodasyidae. Within Thaumastodermatidae Diplodasyinae and Thaumastodermatinae are monophyletic and so are most genera. Oregodasys turns out to be the most basal group within Thaumastodermatinae in analyses of the concatenated data set as well as in analyses of the nuclear genes. Thaumastoderma appears as the sister taxon to the remaining species. Surprisingly, Tetranchyroderma is non-monophyletic in our analyses as one group of species clusters with Ptychostomella while another appears as the sister group of Pseudostomella.

Conclusions/Significance

Results in general agree with the current classification; however, a revision of the more derived thaumastodermatid taxa seems necessary. We also found that the ostensible COI sequences from several species do not conform to the general invertebrate or any other published mitochondrial genetic code; they may be mitochondrially derived nuclear genes (numts), or one or more modifications of the mitochondrial genetic code within Gastrotricha.