Comparative phylogeography of sympatric sister species,Clevelandia ios and Eucyclogobius newberryi (Teleostei,Gobiidae), across the California Transition Zone

It is paradigmatic in marine species that greater dispersal ability often, but not always, results in greater gene flow and less population structure. Some of the exceptions may be attributable to studies confounded by comparison of species with dissimilar evolutionary histories, i.e. co-occurring species that are not closely related or species that are closely related but allopatric. Investigation of sympatric sister species, in contrast, should allow differences in phylogeographic structure to be attributed reliably to recently derived differences in dispersal ability. Here, using mitochondrial DNA control region sequence, we first confirm that Clevelandia ios and Eucyclogobius newberryi are sympatric sister taxa, then demonstrate considerably shallower phylogeographic structure in C. ios than in E. newberryi. This shallower phylogeographic structure is consistent with the higher dispersal ability of C. ios, which most likely results from the interaction of habitat and life-history differences between the species. We suggest that the paradigm will be investigated most rigorously by similar studies of other sympatric sister species, appended by thorough ecological studies, and by extending this sister-taxon approach to comparative phylogeographic studies of monophyletic clades of sympatric species.     

Climate-driven genetic divergence of limpets with different life histories across a southeast African marine biogeographic disjunction: different processes, same outcome

Genetic divergence among populations of marine broadcast spawners in the absence of past geological barriers presents an intriguing challenge to understanding speciation in the sea. To determine how differences in life history affect genetic divergence and demographic histories across incomplete dispersal barriers, we conducted a comparative phylogeographic study of three intertidal limpets (Siphonaria spp.) represented on either side of a biogeographic disjunction separating tropical and subtropical marine provinces in southeastern Africa. Using a combination of mitochondrial and nuclear sequence data, we identified two distinct evolutionary lineages each in both Siphonaria concinna (a planktonic disperser) and S. nigerrima (a direct developer), and panmixia in a second planktonic disperser, S. capensis. Although phylogeographic breaks were present in two species, how these became established differed depending on their life histories. In the direct developer, lack of gene flow following divergence, and demographic expansion from a small initial size in the species' subtropical population, point to a single colonisation event. In contrast, the evolutionary lineages of the planktonic disperser split into two genetic lineages with much larger initial population sizes and southward gene flow continued at least periodically, indicating that divergence in this species may have been driven by a combination of reduced larval dispersal and divergent selection. These findings help explain why the presence or absence of phylogeographic breaks often appears to be independent of species' dispersal potential.     

Phylogeographic analysis reveals a deep lineage split within North Atlantic Littorina saxatilis

Phylogeographic studies provide critical insight into the evolutionary histories of model organisms; yet, to date, range-wide data are lacking for the rough periwinkle Littorina saxatilis, a classic example of marine sympatric speciation. Here, we use mitochondrial DNA (mtDNA) sequence data to demonstrate that L. saxatilis is not monophyletic for this marker, but is composed of two distinct mtDNA lineages (I and II) that are shared with sister species Littorina arcana and Littorina compressa. Bayesian coalescent dating and phylogeographic patterns indicate that both L. saxatilis lineages originated in the eastern North Atlantic, around the British Isles, at approximately 0.64 Ma. Both lineages are now distributed broadly across the eastern, central and western North Atlantic, and show strong phylogeographic structure among regions. The Iberian Peninsula is genetically distinct, suggesting prolonged isolation from northeastern North Atlantic populations. Western North Atlantic populations of L. saxatilis lineages I and II predate the last glacial maximum and have been isolated from eastern North Atlantic populations since that time. This identification of two distinct, broadly distributed mtDNA lineages further complicates observed patterns of repeated incipient ecological speciation in L. saxatilis, because the sympatric origins of distinct ecotype pairs on eastern North Atlantic shores may be confounded by admixture of divergent lineages.     

Comparative phylogeography of Nearctic and Palearctic fishes

Combining phylogeographic data from mitochondrial DNA (mtDNA) of Nearctic and Palearctic freshwater and anadromous fishes, we used a comparative approach to assess the influence of historical events on evolutionary patterns and processes in regional fish faunas. Specifically, we (i) determined whether regional faunas differentially affected by Pleistocene glaciations show predictable differences in phylogeographic patterns; (ii) evaluated how processes of divergence and speciation have been influenced by such differential responses; and (iii) assessed the general contribution of phylogeographic studies to conservation issues. Comparisons among case studies revealed fundamental differences in phylogeographic patterns among regional faunas. Tree topologies were typically deeper for species from nonglaciated regions compared to northern species, whereas species with partially glaciated ranges were intermediate in their characteristics. Phylogeographic patterns were strikingly similar among southern species, whereas species in glaciated areas showed reduced concordance. The extent and locations of secondary contact among mtDNA lineages varied greatly among northern species, resulting in reduced intraspecific concordance of genetic markers for some northern species. Regression analysis of phylogeographic data for 42 species revealed significant latitudinal shifts in intraspecific genetic diversity. Both relative nucleotide diversity and estimates of evolutionary effective population size showed significant breakpoints matching the median latitude for the southern limit of the Pleistocene glaciations. Similarly, analysis of clade depth of phylogenetically distinct lineages vs. area occupied showed that evolutionary dispersal rates of species from glaciated and nonglaciated regions differed by two orders of magnitude. A negative relationship was also found between sequence divergence among sister species as a function of their median distributional latitude, indicating that recent bursts of speciation events have occurred in deglaciated habitats. Phylogeographic evidence for parallel evolution of sympatric northern species pairs in postglacial times suggested that differentiation of cospecific morphotypes may be driven by ecological release. Altogether, these results demonstrate that comparative phylogeography can be used to evaluate not only phylogeographic patterns but also evolutionary processes. As well as having significant implications for conservation programs, this approach enables new avenues of research for examining the regional, historical, and ecological factors involved in shaping intraspecific genetic diversity.     

A hitchhikers guide to the Galápagos: co-phylogeography of Galápagos mockingbirds and their parasites

Background

Parasites are evolutionary hitchhikers whose phylogenies often track the evolutionary history of their hosts. Incongruence in the evolutionary history of closely associated lineages can be explained through a variety of possible events including host switching and host independent speciation. However, in recently diverged lineages stochastic population processes, such as retention of ancestral polymorphism or secondary contact, can also explain discordant genealogies, even in fully co-speciating taxa. The relatively simple biogeographic arrangement of the Galápagos archipelago, compared with mainland biomes, provides a framework to identify stochastic and evolutionary informative components of genealogic data in these recently diverged organisms.

Results

Mitochondrial DNA sequences were obtained for four species of Galápagos mockingbirds and three sympatric species of ectoparasites - two louse and one mite species. These data were complemented with nuclear EF1α sequences in selected samples of parasites and with information from microsatellite loci in the mockingbirds. Mitochondrial sequence data revealed differences in population genetic diversity between all taxa and varying degrees of topological congruence between host and parasite lineages. A very low level of genetic variability and lack of congruence was found in one of the louse parasites, which was excluded from subsequent joint analysis of mitochondrial data. The reconciled multi-species tree obtained from the analysis is congruent with both the nuclear data and the geological history of the islands.

Conclusions

The gene genealogies of Galápagos mockingbirds and two of their ectoparasites show strong phylogeographic correlations, with instances of incongruence mostly explained by ancestral genetic polymorphism. A third parasite genealogy shows low levels of genetic diversity and little evidence of co-phylogeny with their hosts. These differences can mostly be explained by variation in life-history characteristics, primarily host specificity and dispersal capabilities. We show that pooling genetic data from organisms living in close ecological association reveals a more accurate phylogeographic history for these taxa. Our results have implications for the conservation and taxonomy of Galápagos mockingbirds and their parasites.     

Three divergent lineages within an Australian marsupial (Petrogale penicillata) suggest multiple major refugia for mesic taxa in southeast Australia

Mesic southeastern Australia represents the continent's ancestral biome and is highly biodiverse, yet its phylogeographic history remains poorly understood. Here, we examine mitochondrial DNA (mtDNA) control region and microsatellite diversity in the brush‐tailed rock‐wallaby (Petrogale penicillata; n = 279 from 31 sites), to assess historic evolutionary and biogeographic processes in southeastern Australia. Our results (mtDNA, microsatellites) confirmed three geographically discrete and genetically divergent lineages within brush‐tailed rock‐wallabies, whose divergence appears to date to the mid‐Pleistocene. These three lineages had been hypothesized previously but data were limited. While the Northern and Central lineages were separated by a known biogeographic barrier (Hunter Valley), the boundary between the Central and Southern lineages was not. We propose that during particularly cool glacial cycles, the high peaks of the Great Dividing Range and the narrow adjacent coastal plain resulted in a more significant north–south barrier for mesic taxa in southeastern Australia than has been previously appreciated. Similarly, located phylogeographic breaks in codistributed species highlight the importance of these regions in shaping the distribution of biodiversity in southeastern Australia and suggest the existence of three major refuge areas during the Pleistocene. Substructuring within the northern lineage also suggests the occurrence of multiple local refugia during some glacial cycles. Within the three major lineages, most brush‐tailed rock‐wallaby populations were locally highly structured, indicating limited dispersal by both sexes. The three identified lineages represent evolutionarily significant units and should be managed to maximize the retention of genetic diversity within this threatened species.     

The influence of life‐history strategy on genetic differentiation and lineage divergence in darters (Percidae: Etheostomatinae)

Recent studies determined that darters with specialized breeding strategies can exhibit deep lineage divergence over fine geographic scales without apparent physical barriers to gene flow. However, the extent to which intrinsic characteristics interact with extrinsic factors to influence population divergence and lineage diversification in darters is not well understood. This study employed comparative phylogeographic and population genetic methods to investigate the influence of life history on gene flow, dispersal ability, and lineage divergence in two sympatric sister darters with differing breeding strategies. Our results revealed highly disparate phylogeographic histories, patterns of genetic structure, and dispersal abilities between the two species suggesting that life history may contribute to lineage diversification in darters, especially by limiting dispersal among large river courses. Both species also showed striking differences in demographic history, indicating that extrinsic factors differentially affected each species during the Pleistocene. Collectively, our results indicate that intrinsic and extrinsic factors have influenced levels of gene flow among populations within both species examined. However, we suggest that life‐history strategy may play a more important role in lineage diversification in darters than previously appreciated, a finding that has potentially important implications for understanding diversification of the rich North American freshwater fish fauna.     

Phylogeographic patterns in the Australasian genus Chionohebe (Veronica s.l., Plantaginaceae) based on AFLP and chloroplast DNA sequences

The alpine genus Chionohebe is one of seven genera in the southern hemisphere Hebe complex. The main aims of this study were to infer the evolutionary relationships and assess phylogeographic patterns among the six species of Chionohebe, determine the origin of the two species with trans-Tasman distributions, and test species delimitations and specimen identifications based on morphology. Analyses of AFLP data recovered five major lineages within Chionohebe, some of which corresponded to species and varieties as currently circumscribed. Although the cushion chionohebes were strongly supported as monophyletic, the sole non-cushion species, C. densifolia, was sister to Parahebe trifida, and thus the AFLP data do not support a monophyletic Chionohebe as usually circumscribed. Strong north/south and west/east phylogeographic patterns were found among and within the main AFLP lineages in New Zealand. Analyses of chloroplast DNA (cpDNA) revealed eight haplotypes in Chionohebe, but these did not correspond to current taxonomy or geography due to widespread interspecific haplotype sharing. Based on both AFLP and cpDNA results, the two trans-Tasman species are shown to have originated in New Zealand and dispersed to Australia independently.     

Competitive exclusion: phylogeography’s ‘elephant in the room’?

Phylogeographic and evolutionary research programmes have successfully elucidated compelling genetic signatures of earth history. Particularly influential achievements include the demonstration of postglacial recolonization patterns for high-latitude taxa and phylogenetic demonstration of the 'progression rule' along oceanic island chains such as Hawaii. While both of these major biogeographic patterns clearly rely on rapid dispersal over long distances, their phylogeographic detection also apparently relies on the competitive exclusion of secondary dispersers. Such exclusion could occur either between or within species and might reflect fitness differences between lineages or, alternatively, neutral demographic processes (e.g. 'high-density blocking'). Regardless, such spatial genetic patterns would be rapidly eroded were it not for the failure of subsequent dispersers to contribute genetically to newly colonized populations. In addition to its role in revealing colonization patterns, competitive exclusion may also explain the maintenance of historic phylogeographic disjunctions long after the original physical barriers to dispersal have ceased to exist. Additionally, some of the most persuasive evidence of competitive exclusion comes from studies of anthropogenic extinction, where surviving lineages have subsequently expanded their ranges, apparently benefitting from the demise of their prehistoric sisters. Broadly, these biogeographic paradigms are united by the 'disconnect' between dispersal and colonization success, the latter being heavily influenced by inter- and intraspecific competition. Despite its apparent importance, such exclusion (especially within species) has received virtually no attention in the phylogeographic literature. Future studies should aim to test directly for the role of competitive exclusion in constraining the biogeography of highly dispersive taxa.     

Phylogeographic Structure in Penguin Ticks across an Ocean Basin Indicates Allopatric Divergence and Rare Trans-Oceanic Dispersal

The association of ticks (Acarina) and seabirds provides an intriguing system for assessing the influence of long-distance dispersal on the evolution of parasitic species. Recent research has focused on host-parasite evolutionary relationships and dispersal capacity of ticks parasitising flighted seabirds. Evolutionary research on the ticks of non-flighted seabirds is, in contrast, scarce. We conducted the first phylogeographic investigation of a hard tick species (Ixodes eudyptidis) that parasitises the Little Blue Penguin (Eudyptula minor). Using one nuclear (28S) and two mitochondrial (COI and 16S) markers, we assessed genetic diversity among several populations in Australia and a single population on the South Island of New Zealand. Our results reveal two deeply divergent lineages, possibly representing different species: one comprising all New Zealand samples and some from Australia, and the other representing all other samples from Australian sites. No significant population differentiation was observed among any Australian sites from within each major clade, even those separated by hundreds of kilometres of coastline. In contrast, the New Zealand population was significantly different to all samples from Australia. Our phylogenetic results suggest that the New Zealand and Australian populations are effectively isolated from each other; although rare long-distance dispersal events must occur, these are insufficient to maintain trans-Tasman gene flow. Despite the evidence for limited dispersal of penguin ticks between Australia and New Zealand, we found no evidence to suggest that ticks are unable to disperse shorter distances at sea with their hosts, with no pattern of population differentiation found among Australian sites. Our results suggest that terrestrial seabird parasites may be quite capable of short-distance movements, but only sporadic longer-distance (trans-oceanic) dispersal.     

Discordant patterns of genetic and phenotypic differentiation in five grasshopper species codistributed across a microreserve network

Conservation plans can be greatly improved when information on the evolutionary and demographic consequences of habitat fragmentation is available for several codistributed species. Here, we study spatial patterns of phenotypic and genetic variation among five grasshopper species that are codistributed across a network of microreserves but show remarkable differences in dispersal‐related morphology (body size and wing length), degree of habitat specialization and extent of fragmentation of their respective habitats in the study region. In particular, we tested the hypothesis that species with preferences for highly fragmented microhabitats show stronger genetic and phenotypic structure than codistributed generalist taxa inhabiting a continuous matrix of suitable habitat. We also hypothesized a higher resemblance of spatial patterns of genetic and phenotypic variability among species that have experienced a higher degree of habitat fragmentation due to their more similar responses to the parallel large‐scale destruction of their natural habitats. In partial agreement with our first hypothesis, we found that genetic structure, but not phenotypic differentiation, was higher in species linked to highly fragmented habitats. We did not find support for congruent patterns of phenotypic and genetic variability among any studied species, indicating that they show idiosyncratic evolutionary trajectories and distinctive demographic responses to habitat fragmentation across a common landscape. This suggests that conservation practices in networks of protected areas require detailed ecological and evolutionary information on target species to focus management efforts on those taxa that are more sensitive to the effects of habitat fragmentation.     

Contrasting patterns of phylogeographic relationships in sympatric sister species of ironclad beetles (Zopheridae: <Emphasis Type="Italic">Phloeodes</Emphasis> spp.) in California's Transverse Ranges

Background  

Comparative phylogeography of sympatric sibling species provides an opportunity to isolate the effects of geography and demographics on the evolutionary history of two lineages over the same, known time scale. In the current study, we investigated the phylogeographic structure of two zopherid beetle species, Phloeodes diabolicus and P. plicatus, where their ranges overlap in California's Transverse Ranges.     

Phylogeography of a widespread lizard complex reflects patterns of both geographic and ecological isolation

A primary challenge for modern phylogeography is understanding how ecology and geography, both contemporary and historical, shape the spatial distribution and evolutionary histories of species. Phylogeographic patterns are the result of many factors, including geology, climate, habitat, colonization history and lineage‐specific constraints. Assessing the relative influences of these factors is difficult because few species, regions and environments are sampled in enough detail to compare competing hypotheses rigorously and because a particular phylogeographic pattern can potentially result from different evolutionary scenarios. The silky anoles (Anolis sericeus complex) of Central America and Mexico are abundant and found in all types of lowland terrestrial habitat, offering an excellent opportunity to test the relative influences of the factors affecting diversification. Here, we performed a range‐wide statistical phylogeographic analysis on restriction site‐associated DNA (RAD) markers from silky anoles and compared the phylogeographic patterns we recovered to historical and contemporary environmental and topographic data. We constructed niche models to compare niche overlap between sister lineages and conducted coalescent simulations to characterize how the major lineages of silky anoles have diverged. Our results revealed that the mode of divergence for major lineage diversification events was geographic isolation, resulting in ecological divergence between lineages, followed by secondary contact. Moreover, comparisons of parapatric sister lineages suggest that ecological niche divergence contributed to isolation by environment in this system, reflecting the natural history differences among populations in divergent environments.     

Methods in comparative phylogeography, and their application to studying evolution in the north american aridlands

Historical biogeography and comparative phylogeography havemuch in common. Both seek to discover common historical patternsin the elements of biotas, although typically at different tiersof evolutionary history. Comparative phylogeography is basedon phylogeographic analyses of multiple taxa, usually widespreadspecies. By comparing the phylogeographic structures of numerouswidespread sympatric species, one can infer whether the currentfauna has been historically stable, as evidenced by the relativefrequency of geographically congruent reciprocally monophyleticgroups. Alternatively, if species distributions are ephemeralover evolutionary time, a mixture of phylogeographic structuresis expected. Coalescence analyses contribute information abouthistory irrespective of whether haplotype phylogenies are structuredor not. In the aridlands of North America, several isolatingevents are evident in the phylogeographic patterns of birds,mammals and herps. A mid-peninsular seaway in Baja California,dated at ca. one million years before present, had a pervasiveeffect, with 13 of 16 assayed species showing a concordant split.Hence, this community appears to have been a stable assemblageof species over the past one million years. In contrast, theavifauna of the Sonoran-Chihuahuan deserts consists of two specieswith a concordant split and three other species that are undifferentiatedacross both deserts. Hence, the species in this area have haddifferent histories. The Baja biota appears to resemble itsancestral configuration to a greater degree than the Sonoran-Chihuahuanone. A deeper evolutionary event separated taxa in Baja Californiafrom the eastern deserts, showing that the aridlands fauna wasaffected by events at different times resulting in overlaintiers of history.     

Shared geographic histories and dispersal contribute to congruent phylogenies between amphipods and their microsporidian parasites at regional and global scales

In parasites that strongly rely on a host for dispersal, geographic barriers that act on the host will simultaneously influence parasite distribution as well. If their association persists over macroevolutionary time it may result in congruent phylogenetic and phylogeographic patterns due to shared geographic histories. Here, we investigated the level of congruent evolutionary history at a regional and global scale in a highly specialised parasite taxon infecting hosts with limited dispersal abilities: the microsporidians Dictyocoela spp. and their amphipod hosts. Dictyocoela can be transmitted both vertically and horizontally and is the most common microsporidian genus occurring in amphipods in Eurasia. However, little is known about its distribution elsewhere. We started by conducting molecular screening to detect microsporidian parasites in endemic amphipod species in New Zealand; based on phylogenetic analyses, we identified nine species‐level microsporidian taxa including six belonging to Dictyocoela. With a distance‐based cophylogenetic analysis at the regional scale, we identified overall congruent phylogenies between Paracalliope, the most common New Zealand freshwater amphipod taxon, and their Dictyocoela parasites. Also, hosts and parasites showed similar phylogeographic patterns suggesting shared biogeographic histories. Similarly, at a global scale, phylogenies of amphipod hosts and their Dictyocoela parasites showed broadly congruent phylogenies. The observed patterns may have resulted from covicariance and/or codispersal, suggesting that the intimate association between amphipods and Dictyocoela may have persisted over macroevolutionary time. We highlight that shared biogeographic histories could play a role in the codiversification of hosts and parasites at a macroevolutionary scale.     

Comparative phylogeography of three codistributed stomatopods: origins and timing of regional lineage diversification in the Coral Triangle

The Indonesian-Australian Archipelago is the center of the world's marine biodiversity. Although many biogeographers have suggested that this region is a "center of origin," criticism of this theory has focused on the absence of processes promoting lineage diversification in the center. In this study we compare patterns of phylogeographic structure and gene flow in three codistributed, ecologically similar Indo-West Pacific stomatopod (mantis shrimp) species. All three taxa show evidence for limited gene flow across the Maluku Sea with deep genetic breaks between populations from Papua and Northern Indonesia, suggesting that limited water transport across the Maluku Sea may limit larval dispersal and gene flow across this region. All three taxa also show moderate to strong genetic structure between populations from Northern and Southern Indonesia, indicating limited gene flow across the Flores and Java Seas. Despite the similarities in phylogeographic structure, results indicate varied ages of the genetic discontinuities, ranging from the middle Pleistocene to the Pliocene. Concordance of genetic structure across multiple taxa combined with temporal discordance suggests that regional genetic structures have arisen from the action of common physical processes operating over extended time periods. The presence in all three species of both intraspecific genetic structure as well as deeply divergent lineages that likely represent cryptic species suggests that these processes may promote lineage diversification within the Indonesian-Australian Archipelago, providing a potential mechanism for the center of origin. Efforts to conserve biodiversity in the Coral Triangle should work to preserve both existing biodiversity as well as the processes creating the biodiversity.     

Cryptic lineage diversity,body size divergence,and sympatry in a species complex of Australian lizards (Gehyra)

Understanding the joint evolutionary and ecological underpinnings of sympatry among close relatives remains a key challenge in biology. This problem can be addressed through joint phylogenomic and phenotypic analysis of complexes of closely related lineages within, and across, species and hence representing the speciation continuum. For a complex of tropical geckos from northern Australia—Gehyra nana and close relatives—we combine mtDNA phylogeography, exon‐capture sequencing, and morphological data to resolve independently evolving lineages and infer their divergence history and patterns of morphological evolution. Gehyra nana is found to include nine divergent lineages and is paraphyletic with four other species from the Kimberley region of north‐west Australia. Across these 13 taxa, 12 of which are restricted to rocky habitats, several lineages overlap geographically, including on the diverse Kimberley islands. Morphological evolution is dominated by body size shifts, and both body size and shape have evolved gradually across the group. However, larger body size shifts are observed among overlapping taxa than among closely related parapatric lineages of G. nana, and sympatric lineages are more divergent than expected at random. Whether elevated body size differences among sympatric lineages are due to ecological sorting or character displacement remains to be determined.     

Systematics of the lizard family pygopodidae with implications for the diversification of Australian temperate biotas

We conducted a phylogenetic study of pygopodid lizards, a group of 38 species endemic to Australia and New Guinea, with two major goals: to reconstruct a taxonomically complete and robustly supported phylogeny for the group and to use this information to gain insights into the tempo, mode, and timing of the pygopodid radiation. Phylogenetic analyses of mitochondrial DNA (mtDNA), nuclear DNA (nDNA), and previously published morphological data using parsimony, maximum likelihood, and Bayesian methods on the independent and combined three data sets yielded trees with similar and largely stable ingroup topologies. However, relationships among the six most inclusive and unambiguously supported clades (Aprasia, Delma, Lialis, Ophidiocephalus, Pletholax, and Pygopus) varied depending on data set analyzed. We used parametric bootstrapping to help us understand which of the three-branch schemes linking these six taxa was most plausible given our data. We conclude based on our results that the arrangement ((((Delma, Lialis)Pygopus)Pletholax)(Aprasia, Ophidiocephalus)) represents the best hypothesis of intergeneric relationships. A second major problem to arise in our study concerned the inability of our two outgroup taxa (Diplodactylus) to root trees properly; three different rooting locations were suggested depending upon analysis. This long-branch attraction problem was so severe that the outgroup branch also interfered with estimation of ingroup relationships. We therefore used the molecular clock method to root the pygopodid tree. Results of two independent molecular clock analyses (mtDNA and nDNA) converged upon the same root location (branch leading to Delma). We are confident that we have found the correct root because the possibility of our clock estimates agreeing by chance alone is remote given that there are 65 possible root locations (branches) on the pygopodid tree (approximately 1 in 65 odds). Our analysis also indicated that Delma fraseri is not monophyletic, a result supported by a parametric bootstrapping test. We elevated the Western Australian race, Delma f. petersoni, to species status (i.e., Delma petersoni) because hybridization and incomplete lineage sorting could be ruled out as potential causes of this paraphyletic gene tree and because D. grayii is broadly sympatric with its sister species D. fraseri. Climate changes over the past 23 million years, which transformed Australia from a wet, green continent to one that is largely dry and brown, have been suspected as playing a major role in the diversification of Australia's temperate biotas. Our phylogenetic analyses of pygopodid speciation and biogeography revealed four important findings consistent with this climate change diversification model: (1) our fossil-calibrated phylogeny shows that although some extant pygopodid lineages predate the onset of aridification, 28 of 33 pygopodid species included in our study seem to have originated in the last 23 million years; (2) relative cladogenesis tests suggest that several major clades underwent higher than expected rates of speciation; (3) our findings support earlier studies showing that speciation of mesic-adapted biotas in the southeastern and southwestern corners of Australia largely occurred within each of these regions between 12 and 23 million years ago as opposed to repeated dispersal between these regions; and (4) we have identified for the first time the existence of several pairs of sympatric sister species of lizards living in arid and semiarid ecosystems. These sympatric sister species seem to be younger than allopatric or parapatric sister-species pairs, which is not consistent with previous beliefs.