The genetic legacy of aridification: Climate cycling fostered lizard diversification in Australian montane refugia and left low-lying deserts genetically depauperate

Dragons in the tropics - Phylogeography and speciation in Diporiphora lizards and common geographic breaks in co-distributed taxa

Co-distributed taxa can respond both similarly or differently to the same climatic and geological events, resulting in a range of phylogeographic patterns across the region. Using a nested approach on a taxonomically diverse yet morphologically conservative group of agamid lizards, we first aimed to evaluate more precisely the extent of phylogeographic structuring within the genus. Then, focusing on four lineages within the more widespread species, we assessed the impact of biogeographic barriers on phylogeographic structuring and demographic history of species, comparing to patterns previously observed in co-distributed taxa. These species occur in the Australian Monsoonal Tropics, a vast tropical savanna system with high richness and endemism associated with environmental heterogeneity and past climate fluctuations. The employment of genomic data helped to determine the relationships between specific taxa that were previously difficult to place. We found a local influence of biogeographic and climatic breaks on population dynamics, analogous to other species. We detected high levels of population structure in the West Kimberley and Arnhem Plateau, which are already known for high endemism. However, we also highlighted unique lineages in areas that have been overlooked until recently, in the South Kimberley and West Top End. Climatic and geographical features in the Arnhem Plateau act as a soft barrier between populations in the east and west regions of the Top End. These observations reflect patterns observed for other vertebrates across this rich biome, indicating how climatic variation, species' ecology, and landscape features interact to shape regional diversity and endemism.

Estimates of heterozygosity from single nucleotide polymorphism markers are context-dependent and often wrong

Genetic diversity is frequently described using heterozygosity, particularly in a conservation context. Often, it is estimated using single nucleotide polymorphisms (SNPs); however, it has been shown that heterozygosity values calculated from SNPs can be biased by both study design and filtering parameters. Though solutions have been proposed to address these issues, our own work has found them to be inadequate in some circumstances. Here, we aimed to improve the reliability and comparability of heterozygosity estimates, specifically by investigating how sample size and missing data thresholds influenced the calculation of autosomal heterozygosity (heterozygosity calculated from across the genome, i.e. fixed and variable sites). We also explored how the standard practice of tri- and tetra-allelic site exclusion could bias heterozygosity estimates and influence eventual conclusions relating to genetic diversity. Across three distinct taxa (a frog, Litoria rubella; a tree, Eucalyptus microcarpa; and a grasshopper, Keyacris scurra), we found heterozygosity estimates to be meaningfully affected by sample size and missing data thresholds, partly due to the exclusion of tri- and tetra-allelic sites. These biases were inconsistent both between species and populations, with more diverse populations tending to have their estimates more severely affected, thus having potential to dramatically alter interpretations of genetic diversity. We propose a modified framework for calculating heterozygosity that reduces bias and improves the utility of heterozygosity as a measure of genetic diversity, whilst also highlighting the need for existing population genetic pipelines to be adjusted such that tri- and tetra-allelic sites be included in calculations.

A revision of the gudgeon genus Hypseleotris (Gobiiformes: Gobioidei: Eleotridae) of northwest Australia, describing three new species and synonymizing the genus Kimberleyeleotris

Species within the northwest Australian clade of Hypseleotris (six species) and the genus Kimberleyeleotris (two species) are reviewed following the recording of new populations in the region and a molecular study of the group that identified three undescribed candidate species. Based on the analysis of extensive morphological and nuclear and mitochondrial molecular datasets, Kimberleyeleotris is here formally synonymised with Hypseleotris. Furthermore, three species from the Kimberley region, Western Australia, are described to science: Hypseleotris maranda sp. nov., Hypseleotris wunduwala sp. nov., and Hypseleotris garawudjirri sp. nov. The presence of, or number of scales across the head and body, the pattern of sensory papillae on the head, fin ray counts, dorsal and anal fin colouration (particularly in breeding males), and body depth, can be used to distinguish the members of the northwest Australia lineage. Furthermore, the newly described species were genetically separated from all northwest Australian congeners by K2P distances ranging from 7.8-11.3% based on the CO1 gene, and 7.7-16.3 % based on the entire mitochondrial genome. Two of the new species, H. maranda sp. nov. and H. wunduwala sp. nov., have extremely narrow ranges being found in single sub-catchments of the Roe and King Edward Rivers respectively. On the other hand, H. garawudjirri sp. nov. is moderately widespread, being found across the Charnley, Calder, and Sale rivers. While the conservation risk to H. maranda sp. nov. and H. wunduwala sp. nov. is inherently high due to their small range, there are currently no obvious local threatening processes to either of these species given their remote locations that are little impacted by human activities.

Plio–Pleistocene vicariance across arid Australia in the ‘Spiny Knob-tailed Geckos’ (

Across Australia’s monsoon tropics and vast arid zone isolated regions or ‘islands’ of upland or rocky habitat are home to disjunct populations of many taxa of plants and animals. Comparative analyses of lineages that occur across these habitat islands provide opportunities to understand when and how environmental change drove isolation and diversification across arid Australia. Here we present an analysis of mitochondrial genetic diversity across disjunct populations of geckos in the Nephrurus asper group. Dating analyses suggest that disjunct and genetically divergent populations spanning the northern half of Australia diverged through the Plio–Pleistocene. Based on the timing of divergence and current habitat associations we hypothesise that species in this lineage were isolated by the expansion of unsuitable arid-zone habitats from the late Pliocene onwards. Across most areas, these barriers appear to be sandy or stony deserts. However, in eastern Australia genetically divergent populations are separated by grassland on flat vertisol-dominated soils (‘blacksoils’), suggesting that these habitats also expanded during the late Pliocene aridification. Finally, we show that western Queensland populations formerly referred to N. asper are genetically divergent and diagnosable on the basis of colour pattern and, herein, recognise these populations as a distinct species. https://zoobank.org/urn:lsid:zoobank.org:pub:9508CAAA-D014-452D-A3DA-325851615FA7

Cytonuclear discordance, reticulation and cryptic diversity in one of North America's most common frogs

Complicated phylogenetic histories benefit from diverse sources of inference. Pseudacris crucifer (spring peeper) spans most of eastern North America and comprises six mtDNA lineages that form multiple contact zones. The putative Miocene or early Pliocene origins of the oldest lineages within Pseudacris crucifer imply sufficient time for species-level divergence. To understand why this species appears unified while congeners have radiated, we analyze and compare male advertisement calls, mitochondrial, and nuclear markers and speak to the complex processes that have potentially influenced its contemporary patterns. We find extensive geographic and topological mitonuclear discordance, with three nuclear lineages containing 6 more-structured mtDNA lineages, and nuclear introgression at some contact zones. Male advertisement call differentiation is incongruent with the genetic structure as only one lineage appears differentiated. Occupying the Interior Highlands of the central United States, this Western lineage also has the most concordant mitochondrial and nuclear geographic patterns. Based on our findings we suggest that the antiquity of common ancestors was not as important as the maintenance of allopatry in the divergence in P. crucifer genetic lineages. We use multiple lines of evidence to generate hypotheses of isolation, reticulation, and discordance within this species and to expand our understanding of the early stages of speciation.

Multi-Species Phylogeography of Arid-Zone Sminthopsinae (Marsupialia: Dasyuridae) Reveals Evidence of Refugia and Population Expansion in Response to Quaternary Change

Historical population contraction and expansion events associated with Pleistocene climate change are important drivers of intraspecific population structure in Australian arid-zone species. We compared phylogeographic patterns among arid-adapted Dasyuridae (Sminthopsis and Planigale) with close phylogenetic relationships and similar ecological roles to investigate the drivers of phylogeographic structuring and the importance of historical refugia. We generated haplotype networks for two mitochondrial (control region and cytochrome b) and one nuclear (omega-globin) gene from samples distributed across each species range. We used Φ_ST to test for a genetic population structure associated with the four Pilbara subregions, and we used expansion statistics and Bayesian coalescent skyline analysis to test for signals of historical population expansion and the timing of such events. Significant population structure associated with the Pilbara and subregions was detected in the mitochondrial data for most species, but not with the nuclear data. Evidence of population expansion was detected for all species, and it likely began during the mid-late Pleistocene. The timing of population expansion suggests that these species responded favorably to the increased availability of arid habitats during the mid-late Pleistocene, which is when previously patchy habitats became more widespread. We interpret our results to indicate that the Pilbara region could have acted as a refugium for small dasyurids.

Recent range expansion in Australian hummock grasses (Triodia) inferred using genotyping-by-sequencing

The Australian arid zone (AAZ) has undergone aridification and the formation of vast sandy deserts since the mid-Miocene. Studies on AAZ organisms, particularly animals, have shown patterns of mesic ancestry, persistence in rocky refugia and range expansions in arid lineages. There has been limited molecular investigation of plants in the AAZ, particularly of taxa that arrived in Australia after the onset of aridification. Here we investigate populations of the widespread AAZ grass Triodia basedowii to determine whether there is evidence for a recent range expansion, and if so, its source and direction. We also undertake a dating analysis for the species complex to which T. basedowii belongs, in order to place its diversification in relation to changes in AAZ climate and landscapes. We analyse a genomic single nucleotide polymorphism data set from 17 populations of T. basedowii in a recently developed approach for detecting the signal and likely origin of a range expansion. We also use alignments from existing and newly sequenced plastomes from across Poaceae for analysis in BEAST to construct fossil-calibrated phylogenies. Across a range of sampling parameters and outgroups, we detected a consistent signal of westward expansion for T. basedowii, originating in central or eastern Australia. Divergence time estimation indicates that Triodia began to diversify in the late Miocene (crown 7.0-8.8 million years (Ma)), and the T. basedowii complex began to radiate during the Pleistocene (crown 1.4-2.0 Ma). This evidence for range expansion in an arid-adapted plant is consistent with similar patterns in AAZ animals and likely reflects a general response to the opening of new habitat during aridification. Radiation of the T. basedowii complex through the Pleistocene has been associated with preferences for different substrates, providing an explanation why only one lineage is widespread across sandy deserts.

On and off the rocks: persistence and ecological diversification in a tropical Australian lizard radiation

Background

Congruent patterns in the distribution of biodiversity between regions or habitats suggest that key factors such as climatic and topographic variation may predictably shape evolutionary processes. In a number of tropical and arid biomes, genetic analyses are revealing deeper and more localised lineage diversity in rocky ranges than surrounding habitats. Two potential drivers of localised endemism in rocky areas are refugial persistence through climatic change, or ecological diversification and specialisation. Here we examine how patterns of lineage and phenotypic diversity differ across two broad habitat types (rocky ranges and open woodlands) in a small radiation of gecko lizards in the genus Gehyra (the australis group) from the Australian Monsoonal Tropics biome.

Results

Using a suite of approaches for delineating evolutionarily independent lineages, we find between 26 and 41 putative evolutionary units in the australis group (versus eight species currently recognised). Rocky ranges are home to a greater number of lineages that are also relatively more restricted in distribution, while lineages in open woodland habitats are fewer, more widely distributed, and, in one case, show evidence of range expansion. We infer at least two shifts out of rocky ranges and into surrounding woodlands. Phenotypic divergence between rocky ranges specialist and more generalist taxa is detected, but no convergent evolutionary regimes linked to ecology are inferred.

Conclusions

In climatically unstable biomes such as savannahs, rocky ranges have functioned as zones of persistence, generators of diversity and a source of colonists for surrounding areas. Phenotypic divergence can also be linked to the use of differing habitat types, however, the extent to which ecological specialisation is a primary driver or secondary outcome of localised diversification remains uncertain.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1408-1) contains supplementary material, which is available to authorized users.

How the Aridification of Australia Structured the Biogeography and Influenced the Diversification of a Large Lineage of Australian Cicadas

Over the last 30 million years, Australia's landscape has undergone dramatic cooling and drying due to the establishment of the Antarctic Circumpolar Current and change in global CO$_{2}$ levels. Studies have shown that many Australian organisms went extinct during these major cooling events, while others experienced adaptive radiations and increases in diversification rates as a result of exploiting new niches in the arid zone. Despite the many studies on diversification and biogeography in Australia, few have been continent-wide and none have focused on a group of organisms adapted to feeding on plants. We studied 162 species of cicadas in the Australian Pauropsalta complex, a large generic lineage within the tribe Cicadettini. We asked whether there were changes in the diversification rate of Pauropsalta over time and if so: 1) which clades were associated with the rate change? 2) did timing of rate shifts correspond to known periods of dramatic historical climate change, 3) did increases in diversification rate along select lineages correspond to adaptive radiations with movement into the arid zone? To address these questions, we estimated a molecular phylogeny of the Pauropsalta complex using ${\sim}$5300 bp of nucleotide sequence data distributed among five loci (one mtDNA locus and four nDNA loci). We found that this large group of cicadas did not diversify at a constant rate as they spread through Australia; instead the signature of decreasing diversification rate changed roughly around the time of the expansion of the east Antarctic ice sheets ${\sim}$16 Ma and the glaciation of the northern hemisphere ${\sim}$3 Ma. Unlike other Australian taxa, the Pauropsalta complex did not explosively radiate in response to an early invasion of the arid zone. Instead multiple groups invaded the arid zone and experienced rates of diversification similar to mesic-distributed taxa. We found evidence for relictual groups, located in pre-Mesozoic habitat, that have not diversified and continue to reside on mesic hosts in isolated "habitat islands". Future work should focus on groups of similar ages with similar distribution patterns to determine whether this tempo and pattern of diversification and biogeography is consistent with evidence from other phytophagous insects.

Evolutionary and natural history of the turtle frog, Myobatrachus gouldii, a bizarre myobatrachid frog in the southwestern Australian biodiversity hotspot

Southwest Australia (SWA) is a global biodiversity hotspot and a centre of diversity and endemism for the Australo-Papuan myobatrachid frogs. Myobatrachus gouldii (the turtle frog) has a highly derived morphology associated with its forward burrowing behaviour, largely subterranean habit, and unusual mode of reproduction. Its sister genera Metacrinia and Arenophryne have restricted distributions in Western Australia with significant phylogeographic structure, leading to the recent description of a new species in the latter. In contrast, Myobatrachus is distributed widely throughout SWA over multiple climatic zones, but little is known of its population structure, geographic variation in morphology, or reproduction. We generated molecular and morphological data to test for genetic and morphological variation, and to assess whether substrate specialisation in this species may have led to phylogeographic structuring similar to that of other plant and animal taxa in SWA. We assembled sequence data for one mitochondrial and four nuclear DNA loci (3628 base pairs) for 42 turtle frogs sampled throughout their range. Likelihood phylogenetic analyses revealed shallow phylogeographic structure in the mtDNA locus (up to 3.3% genetic distance) and little variation in three of the four nDNA loci. The mtDNA haplotype network suggests five geographically allopatric groups, with no shared haplotypes between regions. These geographic patterns are congruent with several other SWA species, with genetic groups restricted to major hydrological divisions, the Swan Coastal Plain, and the Darling Scarp. The geographically structured genetic groups showed no evidence of significant morphological differentiation (242 individuals), and there was little sexual size dimorphism, but subtle differences in reproductive traits suggest more opportunistic breeding in lower rainfall zones. Call data were compared to sister genera Metacrinia and Arenophryne and found to be highly conservative across the three genera. Like many taxa in SWA, topographic variation and Plio-Pleistocene arid fluctuations likely were historic drivers of diversification in M. gouldii.

Genetic and morphological analyses indicate that the Australian endemic scorpion Urodacus yaschenkoi (Scorpiones: Urodacidae) is a species complex

BACKGROUND

Australian scorpions have received far less attention from researchers than their overseas counterparts. Here we provide the first insight into the molecular variation and evolutionary history of the endemic Australian scorpion Urodacus yaschenkoi. Also known as the inland robust scorpion, it is widely distributed throughout arid zones of the continent and is emerging as a model organism in biomedical research due to the chemical nature of its venom.

METHODS

We employed Bayesian Inference (BI) methods for the phylogenetic reconstructions and divergence dating among lineages, using unique haplotype sequences from two mitochondrial loci (COXI, 16S) and one nuclear locus (28S). We also implemented two DNA taxonomy approaches (GMYC and PTP/dPTP) to evaluate the presence of cryptic species. Linear Discriminant Analysis was used to test whether the linear combination of 21 variables (ratios of morphological measurements) can predict individual's membership to a putative species.

RESULTS

Genetic and morphological data suggest that U. yaschenkoi is a species complex. High statistical support for the monophyly of several divergent lineages was found both at the mitochondrial loci and at a nuclear locus. The extent of mitochondrial divergence between these lineages exceeds estimates of interspecific divergence reported for other scorpion groups. The GMYC model and the PTP/bPTP approach identified major lineages and several sub-lineages as putative species. Ratios of several traits that approximate body shape had a strong predictive power (83-100%) in discriminating two major molecular lineages. A time-calibrated phylogeny dates the early divergence at the onset of continental-wide aridification in late Miocene and Pliocene, with finer-scale phylogeographic patterns emerging during the Pleistocene. This structuring dynamics is congruent with the diversification history of other fauna of the Australian arid zones.

DISCUSSION

Our results indicate that the taxonomic status of U. yaschenkoi requires revision, and we provide recommendations for such future efforts. A complex evolutionary history and extensive diversity highlights the importance of conserving U. yaschenkoi populations from different Australian arid zones in order to preserve patterns of endemism and evolutionary potential.

Phylogeographic evidence for evolutionary refugia in the Gulf sandstone ranges of northern Australia

Locating and protecting climate change refugia is important to conserving biodiversity with accelerating climate change. Comparative phylogeographic analysis provides an effective tool for locating such refugia, as long-term retention of one or more populations within a refugial landscape will generate unique genetic lineages. The ranges of the western Gulf region of northern Australia are thought to represent a significant arid-zone refugium, in which case low-dispersal organisms should have strong phylogeographic structure across the region. To test for this, we conducted extensive sampling of three species of Gehyra geckos and analysed diversity for mitochondrial DNA and eight nuclear loci. These analyses revealed congruent and high phylogeographic diversity, especially, but not exclusively, in rock-restricted species. This finding, and other recent phylogeographic evidence, demonstrates that these topographically variable landforms have enabled persistence of ecologically diverse vertebrate species through the climate changes of the late Pleistocene. Identification of this relatively under-protected region as a significant climate change refugium points to the need to expand protected areas in this region and to invest in ecological management across existing National Parks and Indigenous Protected Areas.

Contrasting Influences of Geographic Range and Distribution of Populations on Patterns of Genetic Diversity in Two Sympatric Pilbara Acacias

The influence of geographic range on species persistence has long been of interest and there is a need for a better understanding of the genetic consequences for species with restricted distributions, particularly with the increasing rate of global species extinctions. However, the genetic effects of restricted range are often confounded by the impacts of population distribution. We compared chloroplast and nuclear genetic diversity and differentiation in two acacias, the restricted, patchily distributed Acacia atkinsiana and the widespread, semi-continuously distributed A. ancistrocarpa. Lower intra-population diversity and higher differentiation between populations were seen in A. atkinsiana compared to its widespread congener, A. ancistrocarpa. There was little evidence of geographical influences on population genetic structure in A. ancistrocarpa whereas A. atkinsiana exhibited nuclear genetic structure with isolation by distance, differentiation of near-coastal populations from those in the ranges, and differentiation of peripheral populations from those in the centre of the distribution. These results are consistent with expectations of the effect of geographic range and population distribution on genetic diversity, but indicate that distribution of populations rather than geographic range has influenced the observed genetic structure. The contrasting patterns observed here demonstrate that conservation approaches for species management and ecological restoration need to consider the distribution of populations in geographically restricted species.

Young relicts and old relicts: a novel palaeoendemic vertebrate from the Australian Central Uplands

Climatic change, and in particular aridification, has played a dominant role in shaping Southern Hemisphere biotas since the mid-Neogene. In Australia, ancient and geologically stable ranges within the vast arid zone have functioned as refugia for populations of mesic taxa extirpated from surrounding areas, yet the extent to which relicts may be linked to major aridification events before or after the Pliocene has not been examined in detail. Here we use molecular phylogenetic and morphological data to show that isolated populations of saxicoline geckos in the genus Oedura from the Australian Central Uplands, formerly confounded as a single taxon, actually comprise two divergent species with contrasting histories of isolation. The recently resurrected Oedura cincta has close relatives occurring elsewhere in the Australian arid biomes with estimated divergence dates concentrated in the early Pliocene. A new taxon (described herein) diverged from all extant Oedura much earlier, well before the end of the Miocene. A review of data for Central Uplands endemic vertebrates shows that for most (including Oedura cincta), gene flow with other parts of Australia probably occurred until at least the very late Miocene or Pliocene. There are, however, a small number of palaeoendemic taxa-often ecologically specialized forms-that show evidence of having persisted since earlier intensification of aridity in the late Miocene.

Snails in the desert: Assessing the mitochondrial and morphological diversity and the influence of aestivation behavior on lineage differentiation in the Australian endemic Granulomelon Iredale, 1933 (Stylommatophora: Camaenidae)

Progressive aridification since the mid-Miocene has had a significant influence on the evolution of the biota in the arid zone of central Australia. Especially moisture sensitive groups, such as snails, are often restricted to topographically complex areas, which have acted as refugia in an otherwise inhospitable environment. This historical fragmentation is deemed to be a potent agent of allopatric lineage diversification. Camaenid land snails are amongst only a few terrestrial gastropods that have managed to survive in the arid zone probably due to their ability to escape desiccation through aestivation. Here, we present the first study of the mitochondrial lineage differentiation in an endemic land snail genus from the Australian 'Red Centre', Granulomelon Iredale, 1933. Exposing significant incongruence between mtDNA phylogeny and morphology-based taxonomy, we completely revise the species and genus level taxonomy of this camaenid group. We demonstrate that this genus contains three species, G. grandituberculatum, G. adcockianum and G. squamulosum, which have so far been assigned to different genera: Granulomelon Iredale, 1933 (junior synonym: Baccalena Iredale, 1937), Basedowena Iredale, 1937 and Pleuroxia Ancey, 1887. Two of these species are widespread comprising multiple divergent mitochondrial lineages. Based on a molecular clock estimate, these lineages diverged approximately during the mid-Pleistocene, a period of particularly severe aridification. The phylogeographic patterns are consistent with an isolation-by-distance model in one species but not the other. We suggest that these differences can be attributed to their distinctive aestivation behavior.

Surviving in mountain climate refugia: new insights from the genetic diversity and structure of the relict shrub Myrtus nivellei (Myrtaceae) in the Sahara Desert

The identification of past glacial refugia has become a key topic for conservation under environmental change, since they contribute importantly to shaping current patterns of biodiversity. However, little attention has been paid so far to interglacial refugia despite their key role for the survival of relict species currently occurring in climate refugia. Here, we focus on the genetic consequences of range contraction on the relict populations of the evergreen shrub Myrtus nivellei, endemic in the Saharan mountains since at least the end of the last Green Sahara period, around 5.5 ka B.P. Multilocus genotypes (nuclear microsatellites and AFLP) were obtained from 215 individuals collected from 23 wadis (temporary rivers) in the three main mountain ranges in southern Algeria (the Hoggar, Tassili n'Ajjer and Tassili n'Immidir ranges). Identical genotypes were found in several plants growing far apart within the same wadis, a pattern taken as evidence of clonality. Multivariate analyses and Bayesian clustering revealed that genetic diversity was mainly structured among the mountain ranges, while low isolation by distance was observed within each mountain range. The range contraction induced by the last episode of aridification has likely increased the genetic isolation of the populations of M. nivellei, without greatly affecting the genetic diversity of the species as a whole. The pattern of genetic diversity observed here suggests that high connectivity may have prevailed during humid periods, which is consistent with recent paleoenvironmental reconstructions.

Diversification and Demography of the Oriental Garden Lizard (Calotes versicolor) on Hainan Island and the Adjacent Mainland

The Oriental garden lizard (Calotes versicolor) is one of the few non-gekkonid lizards that are geographically widespread in the tropics. We investigated its population dynamics on Hainan Island and the adjacent mainland of China and Vietnam, focusing on the impact of cyclic upheaval and submergence of land bridges during the Pleistocene. Our Bayesian phylogenetic analysis reveals two mitochondrial lineages, A and B, which are estimated to have coalesced about 0.26 million years ago (95% credibility interval: 0.05–0.61 million years ago). Lineage A contains individuals mainly from central and southern Wuzhi Mountain on Hainan Island, whereas lineage B mainly comprises individuals from other sites on the island plus the adjacent mainland. The estimated coalescence times within lineages A (0.05 million years ago) and B (0.13 million years ago) fall within a period of cyclical land-bridge formation and disappearance in the Pleistocene. A spatial analysis of molecular variance identified two distinct population groupings: I, primarily containing lineage A, and II, mainly consisting of lineage B. However, haplotypes from lineages A and B occur sympatrically, suggesting that gene flow is ongoing. Neither Wuzhi Mountain nor Qiongzhou Strait and Gulf of Tonkin act as barriers to gene flow among C. versicolor populations. Analyses of the data using mismatch distributions and extended Bayesian skyline plots provide evidence of a relatively stable population size through time for Group I, and moderate population expansions and contractions during the end of the Pleistocene for Group II. We conclude that the phylogeographical patterns of C. versicolor are the combined product of Pleistocene sea-level oscillations and nonphysical barriers to gene flow.