Extreme habitats that emerged in the Pleistocene triggered divergence of weedy Youngia (Asteraceae) in Taiwan

Weeds with broad distributions and large morphological variation are challenging for systematists and evolutionarily intriguing because their intensive dispersal would likely prevent local morphological differentiation. Study on weeds will help to understand divergence in plants unlikely to be affected by geographical and ecological barriers. We studied Youngia japonica based on nrDNA and cpDNA sequences. This is a widespread native in Asia and invasive worldwide; nevertheless, three subspecies (japonica, longiflora, and formosana) and an undescribed variant occur in Taiwan. Bayesian and the most parsimonious phylogenies revealed that subspecies longiflora is a different linage and independently arrived in Taiwan during the Pleistocene via land connection to the Asian Continent. Bayesian time estimation suggested that Youngia in Taiwan diverged in the lower Pleistocene or more recently. Extreme habitats that emerged in the Pleistocene, i.e., cold mountain ranges for subspecies formosana and xeric, raised coral reefs for the undescribed Youngia variant probably had triggered the divergence. Components of Youngia in Taiwan are not monophyletic; a coalescent-based test suggested incomplete lineage sorting. Nevertheless, the samples within each taxon share unique morphological features suggesting a common gene pool and each taxon has different dominant ITS and/or cpDNA types; these conditions suggest ongoing process toward monophyly via coalescent processes and support the delimitation of intraspecific taxa.     

Ecological divergence among sympatric colour morphs in blood pythons, Python brongersmai

Blood pythons in northeastern Sumatra display a series of discrete colour morphs, even among hatchlings within a single clutch. The first step towards understanding the maintenance of this polymorphism is to test the null hypothesis that colour variation in this species has no major biological correlates. Data on >2,000 blood pythons killed for the commercial leather industry enabled us to test, and reject, this hypothesis. The four colour morphs differed significantly in most of the traits that we measured, including temporal and spatial abundances, sex ratios, age structures, mean adult body sizes, body shapes (tail length and body mass relative to snout-vent length), energy stores, numbers of gut parasites, prey types, feeding frequencies and clutch sizes. The causal basis for these associations remains unclear, but is likely to involve three processes: direct effects of colour, linkages between genes for colour and other traits, and correlated spatial heterogeneity in colour, morphology and ecology. The colour polymorphism may be maintained by frequency-dependent selection and genotype-specific habitat selection, because these sedentary ambush predators are under strong selection for effective camouflage to hide them from both predators and potential prey. In support of this hypothesis, similar colour polymorphisms have evolved independently in several other snake taxa that rely upon ambush predation. Received: 18 December 1997 / Accepted: 23 March 1998     

Discordant patterns of evolutionary differentiation in two Neotropical treefrogs

Comparative studies of codistributed taxa test the degree to which historical processes have shaped contemporary population structure. Discordant patterns of lineage divergence among taxa indicate that species differ in their response to common historical processes. The complex geologic landscape of the Isthmus of Central America provides an ideal setting to test the effects of vicariance and other biogeographic factors on population history. We compared divergence patterns between two codistributed Neotropical frogs ( Dendropsophus ebraccatus and Agalychnis callidryas ) that exhibit colour pattern polymorphisms among populations, and found significant differences between them in phenotypic and genetic divergence among populations. Colour pattern in D. ebraccatus did not vary with genetic or geographic distance, while colour pattern co-varied with patterns of gene flow in A. callidryas . In addition, we detected significant species differences in the phylogenetic history of populations, gene flow among them, and the extent to which historical diversification and recent gene flow have been restricted by five biogeographic barriers in Costa Rica and Panama. We inferred that alternate microevolutionary processes explain the unique patterns of diversification in each taxon. Our study underscores how differences in selective regimes and species-typical ecological and life-history traits maintain spatial patterns of diversification.     

Population genetics and conservation of the threatened southeastern beach mouse (Peromyscus polionotus niveiventris): subspecies and evolutionary units

We investigated genetic diversity within the southeastern beach mouse (SEBM-Peromyscus polionotus niveiventris) and also tested the hypothesis that the subspecies recognition of P.p. niveiventris, based on size and color differences, is congruent with this taxon representing a discrete evolutionary lineage. We used ten polymorphic microsatellite loci and mitochondrial cytochrome-b gene DNA sequences to investigate genetic diversity and population structure within the SEBM, and to determine the level of divergence between the SEBM and the nearest known inland subspecies of the oldfield mouse (Peromyscus polionotus rhoadsi). Moderate genetic distances were observed between the SEBM and the inland oldfield mouse based on microsatellite data, with F _ST values ranging from 0.11 to 0.22 between these taxa. Additionally, mitochondrial DNA haplotypes of the SEBM formed a distinct monophyletic group relative to haplotypes sampled from P. p. rhoadsi. Based on previous estimates of rates of mitochondrial DNA evolution in rodents, we inferred that Pleistocene sea-level fluctuations are likely responsible for the historical isolation of the SEBM lineage from mainland P. polionotus. Our data demonstrate the genetic distinctiveness of the SEBM, justifying the current subspecies designation for the SEBM and its continued protection under the United States Endangered Species Act. We classify the Cape Canaveral and Smyrna Dunes Park populations of SEBM as a single evolutionary significant unit. The two known extant allopatric populations of the SEBM showed some differentiation in microsatellite frequencies and were moderately reciprocally distinguishable based on assignment to distinct genetic clusters by a Bayesian admixture procedure. These results justify the classification of these two extant SEBM populations as distinct management units that should be independent targets of management and conservation attention.     

Allozyme variation and differentiation in African and Indian rhinoceroses

We studied variation at 25 to 31 allozymic loci in African and Asian rhinoceroses. Four taxa in three genera were examined: African Ceratotherium simum simum (northern white rhinoceros), C. s. cottoni (southern white rhinoceros), Diceros bicornis (black rhinoceros), and Rhinoceros unicornis (Indian rhinoceros). Extremely small amounts of intraspecific variation were observed in sample sizes of 2 to 10 presumably unrelated individuals per taxon: P = .00-.10, H = 0.00-0.02. We examined demographic bottlenecks and sampling errors as possible reasons for the low levels of detectable variation. The very small intraspecific genetic distance (D = 0.005) between the two living white rhinoceros subspecies is far less than the distance that has been reported for other mammal subspecies. The mean D value of 0.32 +/- 0.11 between the two African genera was also less than expected given the divergence time of greater than 7 million years suggested by the fossil record. Rhinoceroses may be evolving more slowly at the structural gene loci than are some other mammal groups. The estimate of D = 1.05 +/- 0.24 for the African-Indian split supports this idea, as the lineage diverged at least 26 million years ago. Our results contribute to the currently available scientific information on which management decisions aimed toward saving endangered rhinoceroses should be based.     

Plant performance and morpho-functional differentiation in response to edaphic variation in Iberian columbines: cues for range distribution?

Aims Experimental studies of the response to environmental variation of closely related taxa are needed to understand the mechanisms underlying phenotypic divergence, habitat segregation and range span within a radiating genus. We explored the magnitude of phenotypic differentiation and adaptive plasticity in relation to edaphic variation and its significance for habitat segregation and range span in Iberian columbines.Methods We performed a common garden experiment varying soil type (calcareous vs. siliceous) and depth (rocky-shallow vs. deep) with two pairs of widespread and narrowly distributed subspecies of Iberian columbines: Aquilegia vulgaris (subspp. vulgaris and nevadensis) and Aquilegia pyrenaica (subspp. pyrenaica and cazorlensis). We compared tolerance to edaphic variation, trait differentiation and adaptive plasticity of 10 morpho-functional traits between species and distribution ranges. Additionally, we obtained estimates of the competitive environment faced by each taxon in two to four populations per taxa.Important findings Results partially support hypotheses of higher competitive ability and tolerance to edaphic variation in widespread than in narrow endemic taxa. At the species level, the widely distributed taxon was the most tolerant to edaphic variation. Within species, no consistent pattern emerged since the population from the most widespread subspecies was the most tolerant in A. vulgaris but not in A. pyrenaica. Columbines were differentiated in many traits at species and range level. However, the pattern of differentiation does not fully support the hypothesis of higher specialization and stress tolerance in narrow endemics. Although plasticity was generally low, the results support the hypothesis of adaptive plasticity in widespread but not in restricted taxa at least at species level. Ecological differences (adaptive plasticity and competitive ability among others) may have contributed to phenotypic divergence and edaphic niche segregation, as well as to differences in range span among columbines.     

Conservation Implications of Rapid Shifts in the Trade of Wild African and Asian Pythons

Pythons are harvested for the international leather industry and pet trade. We analyzed the CITES export records (1999–2008) of the most intensively commercialized wild pythons (Python regius, Python sebae, Python reticulatus, Python molurus, Python curtus species complex) from African and Asian countries where reliable data on trade rates and ecology are available. Mean declared annual numbers of exported pythons were 30,000 in five African countries and 164,000 in Indonesia. Trade intensity tripled in Indonesia over the last decade, but declined in Africa. African international trade is exclusively associated with the pet market (mainly United States and Europe), whereas Asian pythons are sold mostly for luxury leather, albeit more recently also for the pet trade. A negative correlation between the annual numbers of pythons traded in Africa vs. Asia suggests a rapid and recent shift of the pressure exerted on wild populations in the two main exporting continents. We also found a strong effect of the currency exchange rate (i.e., U.S.$/€, the currencies used by the major importing countries) on African python exports: when the cost per African python increased, importers relied increasingly on Asian providers for pet trade. Overall, our data indicate that Asian pythons (especially P. reticulatus) might be threatened due to the rapidly increasing pressure, whereas the decreasing international trade in African pythons is likely more sustainable.     

Lineage Divergence and Historical Gene Flow in the Chinese Horseshoe Bat (Rhinolophus sinicus)

Closely related taxa living in sympatry provide good opportunities to investigate the origin of barriers to gene flow as well as the extent of reproductive isolation. The only two recognized subspecies of the Chinese rufous horseshoe bat Rhinolophus sinicus are characterized by unusual relative distributions in which R. s. septentrionalis is restricted to a small area within the much wider range of its sister taxon R. s. sinicus. To determine the history of lineage divergence and gene flow between these taxa, we applied phylogenetic, demographic and coalescent analyses to multi-locus datasets. MtDNA gene genealogies and microsatellite-based clustering together revealed three divergent lineages of sinicus, corresponding to Central China, East China and the offshore Hainan Island. However, the central lineage of sinicus showed a closer relationship with septentrionalis than with other lineages of R. s. sinicus, in contrary to morphological data. Paraphyly of sinicus could result from either past asymmetric mtDNA introgression between these two taxa, or could suggest septentrionalis evolved in situ from its more widespread sister subspecies. To test between these hypotheses, we applied coalescent-based phylogenetic reconstruction and Approximate Bayesian Computation (ABC). We found that septentrionalis is likely to be the ancestral taxon and therefore a recent origin of this subspecies can be ruled out. On the other hand, we found a clear signature of asymmetric mtDNA gene flow from septentrionalis into central populations of sinicus yet no nuclear gene flow, thus strongly pointing to historical mtDNA introgression. We suggest that the observed deeply divergent lineages within R. sinicus probably evolved in isolation in separate Pleistocene refugia, although their close phylogeographic correspondence with distinct eco-environmental zones suggests that divergent selection might also have promoted broad patterns of population genetic structure.     

Harvest Effects on Blood Pythons in North Sumatra

Wildlife populations exposed to intense harvesting often exhibit shifts in attributes of the commercial offtake (e.g., numbers, body sizes, body condition, sex ratios, size at maturation, reproductive frequency). To look for such changes, we examined >2,500 specimens of field-collected blood pythons (Python brongersmai) brought to 2 processing facilities in North Sumatra over 2 survey periods (1996–1997 and 2014–2015). Over the 18 years between our surveys, approximately 900,000 blood pythons were taken from Indonesia for the commercial trade; North Sumatra accounted for about 35% of that trade volume. Between survey periods, numbers of snakes brought to 1 processing facility decreased and mean body size increased, whereas at the second facility numbers and sizes of snakes remained the same. Overall, the pythons collected in 2015–2016 were thinner-bodied than in 1996–1997, with fewer immature animals and a reduced size at maturation and reproductive frequency among females. Clutch size relative to maternal size decreased also. These temporal shifts in python demographics likely reflect a response to harvesting, although ecosystem changes may have played a role. Compared to other squamate reptiles, blood pythons may be more vulnerable to over-harvesting because most snakes are intensively collected from discrete habitat patches (oil palm plantations), and the species can be targeted by specific hunting methods. To ensure sustainability of the harvest, we recommend minimum and maximum size limits for skins used in the commercial trade, and more intensive monitoring at processing facilities. © 2019 The Wildlife Society.     

Floral trait differentiation in Anacamptis coriophora: Phenotypic selection on scents,but not on colour

Current divergent selection may promote floral trait differentiation among conspecific populations in flowering plants. However, whether this applies to complex traits such as colour or scents has been little studied, even though these traits often vary within species. In this study, we compared floral colour and odour as well as selective pressures imposed upon these traits among seven populations belonging to three subspecies of the widespread, generalist orchid Anacamptis coriophora. Colour was characterized using calibrated photographs, and scents were sampled using dynamic headspace extraction and analysed using gas chromatography–mass spectrometry. We then quantified phenotypic selection exerted on these traits by regressing fruit set values on floral trait values. We showed that the three studied subspecies were characterized by different floral colour and odour, with one of the two predominant floral volatiles emitted by each subspecies being taxon‐specific. Plant size was positively correlated with fruit set in most populations, whereas we found no apparent link between floral colour and female reproductive success. We detected positive selection on several taxon‐specific compounds in A. coriophora subsp. fragrans, whereas no selection was found on floral volatiles of A. coriophora subsp. coriophora and A. coriophora subsp. martrinii. This study is one of the first to document variation in phenotypic selection exerted on floral scents among conspecific populations. Our results suggest that selection could contribute to ongoing chemical divergence among A. coriophora subspecies.     

A new cryptic species of pond turtle from southern Italy, the hottest spot in the range of the genus Emys (Reptilia, Testudines, Emydidae)

Geographic variation in the mtDNA haplotypes (cytochrome b gene) of 127 European pond turtles from Italy was investigated. Thirty‐eight of the Italian samples were also studied by nuclear fingerprinting (ISSR PCR) and compared with samples from other parts of the range representing all nine currently known mtDNA lineages of Emys orbicularis. Our genetic findings were compared against morphological data sets (measurements, colour pattern) for 109 adult turtles from southern Italy. Italy is displaying on a small geographical scale the most complicated variation known over the entire distributional area of Emys (North Africa over Europe and Asia Minor to the Caspian and Aral Seas). The Tyrrhenic coast of the Apennine Peninsula, the Mt. Pollino area and Basilicata are inhabited by Emys orbicularis galloitalica, a subspecies harbouring a distinct mtDNA lineage. The same lineage is also found in Sardinia. Along the Adriatic coast of Italy and on the Salentine Peninsula (Apulia, southern Italy), another morphologically distinctive subspecies (Emys orbicularis hellenica) occurs, which also bears a different mtDNA lineage. A higher diversity of mtDNA haplotypes in the south of the Apennine Peninsula suggests that the glacial refugia of E. o. galloitalica and E. o. hellenica were located here. A further refuge of E. o. hellenica probably existed in the southern Balkans. The west coasts of the Balkans and Corfu have probably been colonized from Italy and not from the geographically closer southern Balkanic refuge. In Sicily, a third mtDNA lineage is distributed, which is sister to all other known lineages of Emys. Morphologically, Sicilian pond turtles resemble E. o. galloitalica. However, nuclear fingerprinting revealed a clear distinctiveness of the Sicilian taxon, whereas no significant divergence was detected between representatives of the other eight mtDNA lineages of Emys. Furthermore, nuclear fingerprinting provided no evidence for current or past gene flow between the Sicilian taxon and the mainland subspecies of E. orbicularis. Therefore, Sicilian pond turtles are described here as a species new to science. Some populations in Calabria and on the Salentine Peninsula comprise individuals of different mtDNA lineages. We interpret this as a natural contact. However, we cannot exclude that these syntopic occurrences are the result of human activity. For example, in other parts of Italy, the natural distribution pattern of Emys is obscured by allochthonous turtles. This could also be true for southern Italy. The discovery of the complex taxonomic differentiation in southern Italy requires reconsidering conservation strategies.     

Distinctiveness in the face of gene flow: hybridization between specialist and generalist gartersnakes

Patterns of divergence and polymorphism across hybrid zones can provide important clues as to their origin and maintenance. Unimodal hybrid zones or hybrid swarms are composed predominantly of recombinant individuals whose genomes are patchworks of alleles derived from each parental lineage. In contrast, bimodal hybrid zones contain few identifiable hybrids; most individuals fall within distinct genetic clusters. Distinguishing between hybrid swarms and bimodal hybrid zones can be important for taxonomic and conservation decisions regarding the status and value of hybrid populations. In addition, the causes of bimodality are important in understanding the generation and maintenance of biological diversity. For example, are distinct clusters mostly reproductively isolated and co‐adapted gene complexes, or can distinctiveness be maintained by a few ‘genomic islands’ despite rampant gene flow across much of the genome? Here we focus on three patterns of distinctiveness in the face of gene flow between gartersnake taxa in the Great Lakes region of North America. Bimodality, the persistence of distinct clusters of genotypes, requires strong barriers to gene flow and supports recognition of distinct specialist (Thamnophis butleri) and generalist (Thamnophis radix) taxa. Concordance of DNA‐based clusters with morphometrics supports the hypothesis that trophic morphology is a key component of divergence. Finally, disparity in the level of differentiation across molecular markers (amplified fragment length polymorphisms) indicates that distinctiveness is maintained by strong selection on a few traits despite high gene flow currently or in the recent past.     

Mimetic colour pattern evolution in the highly polymorphic Bombus trifasciatus (Hymenoptera: Apidae) species complex and its comimics

Müllerian mimetic systems have uncovered some of the dynamic processes by which natural selection can drive the radiation of convergent and divergent phenotypes. We examined evolution involving Müllerian mimicry in bumble bees by documenting the distribution and evolution of colour patterns amongst three colour‐polymorphic lineages –Bombus trifasciatus Smith, Bombus haemorrhoidalis Smith, and Bombus breviceps Smith – that mimic each other across ～14 colour groups in South‐East Asia. Using mitochondrial DNA sequence data, we estimated relationships within each lineage to infer the processes that gave rise to the colour diversity and develop hypotheses on species recognition. We expanded on our assessment of species delineation in the B. trifasciatus lineage using three nuclear gene fragments and morphometrics. Comparison of colour patterns amongst georeferenced specimens showed considerable variation in the degree and geographical range of mimicry amongst mimicry groups. Phylogenetic estimates show high rates of colour pattern evolution, with colour variation often exceeding variation within the fast‐evolving mitochondrial genes. The molecular data, and to some degree the morphometric data, support unique histories for several taxa recognized previously within the B. trifasciatus lineage, which may include several species. Early vicariant events within the B. trifasciatus lineage are likely to have occurred ～2.2 Mya in the mountains of south‐west China. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 166 , 805–826.     

Python phylogenetics: inference from morphology and mitochondrial DNA

We used nucleotide sequences from four mitochondrial genes and structural features of the mitochondrial control region, combined with a revised, previously published, morphological data set to infer phylogenetic relationships among the pythons. We aimed to determine which of two competing hypotheses of relationships of the genera Aspidites and Python best explains the evolutionary and bioegeographical history of the family. All analyses of the combined data recover a set of relationships in which (1) the genus Python is paraphyletic with the two east Asian species, P. reticulatus and P. timoriensis , as the sister lineage to the seven Australo-Papuan python genera. We support recognition of a distinct genus for the P. reticulatus  +  P. timoriensis clade; (2) the remaining species of the genus Python form a clade which is the sister lineage to the remainder of the family; (3) the genus Aspidites is embedded among the Australo-Papuan genera. The seemingly primitive characteristics of Aspidites may be better interpreted as reversals or specializations that have accompanied a switch to burrowing in this genus. Resolution of the relationships among the Australo-Papuan lineages is weak, possibly because of rapid diversification early in the history of the radiation. We assessed the tempo of the Indo-Australian python radiation using a maximum likelihood framework based on the birth–death process. We find strong support for elevated speciation rates during the period when Australia collided with the proto-Indonesian archipelago. The data support an origin for pythons outside Australia, followed by a radiation into Australia during the mid-Tertiary.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 93 , 603–619.     

Niche divergence and lineage diversification among closely related Sistrurus rattlesnakes

Comparing niche divergence among closely related taxa can yield important insights into the ecological distinctiveness of genetically similar forms, and identify the processes that are responsible for diversification in such organisms. Here, we apply newly developed techniques for analysing niche divergence to assess how ecologically distinct a group of closely related rattlesnakes (Sistrurus sp.) are and to explore the role that niche divergence may have played in their diversification. We find that all taxa even the most recently evolved subspecies (approximately 100,000 years old) are now ecologically distinct, implying a role for ecology in the diversification process. Statistical analysis based on comparisons with null models show that niche divergence between forms is more common than niche conservation. Finally, there is nonlinear relationship between phylogenetic and niche divergence in this group whereby niche divergence develops more rapidly between recently diverged subspecies than more distantly related forms. Overall, our results argue that ecology may play an important role in the diversification process in these snakes.     

Ancestral polymorphisms in genetic markers obscure detection of evolutionarily distinct populations in the endangered Florida grasshopper sparrow (Ammodramus savannarum floridanus)

Genetic analyses of bird subspecies designated as conservation units can address whether they represent units with independent evolutionary histories and provide insights into the evolutionary processes that determine the degree to which they are genetically distinct. Here we use mitochondrial DNA control region sequence and six microsatellite DNA loci to examine phylogeographical structure and genetic differentiation among five North American grasshopper sparrow (Ammodramus savannarum) populations representing three subspecies, including a population of the endangered Florida subspecies (A. s. floridanus). This federally listed taxon is of particular interest because it differs phenotypically from other subspecies in plumage and behaviour and has also undergone a drastic decline in population size over the past century. Despite this designation, we observed no phylogeographical structure among populations in either marker: mtDNA haplotypes and microsatellite genotypes from floridanus samples did not form clades that were phylogenetically distinct from variants found in other subspecies. However, there was low but significant differentiation between Florida and all other populations combined in both mtDNA (FST = 0.069) and in one measure of microsatellite differentiation (theta = 0.016), while the non-Florida populations were not different from each other. Based on analyses of mtDNA variation using a coalescent-based model, the effective sizes of these populations are large (approximately 80,000 females) and they have only recently diverged from each other (< 26,000 ybp). These populations are probably far from genetic equilibrium and therefore the lack of phylogenetic distinctiveness of the floridanus subspecies and minimal genetic differentiation is due most probably to retained ancestral polymorphism. Finally, levels of variation in Florida were similar to other populations supporting the idea that the drastic reduction in population size which has occurred within the last 100 years has not yet had an impact on levels of variation in floridanus. We argue that despite the lack of phylogenetic distinctiveness of floridanus genotypes the observed genetic differentiation and previously documented phenotypic differences justify continued designation of this subspecies as a protected population segment.     

Strikingly variable divergence times inferred across an Amazonian butterfly 'suture zone'

'Suture zones' are areas where hybrid and contact zones of multiple taxa are clustered. Such zones have been regarded as strong evidence for allopatric divergence by proponents of the Pleistocene forest refugia theory, a vicariance hypothesis frequently used to explain diversification in the Amazon basin. A central prediction of the refugia and other vicariance theories is that the taxa should have a common history so that divergence times should be coincident among taxa. A suture zone for Ithomiinae butterflies near Tarapoto, NE Peru, was therefore studied to examine divergence times of taxa in contact across the zone. We sequenced 1619bp of the mitochondrial COI/COII region in 172 individuals of 31 species from across the suture zone. Inferred divergence times differed remarkably, with divergence between some pairs of widespread species (each of which may have two or more subspecies interacting in the zone, as in the genus Melinaea) being considerably less than that between hybridizing subspecies in other genera (for instance in Oleria). Our data therefore strongly refute a simple hypothesis of simultaneous vicariance and suggest that ongoing parapatric or other modes of differentiation in continuous forest may be important in driving diversification in Amazonia.     

Phylogeography of speciation: allopatric divergence and secondary contact between outcrossing and selfing Clarkia

The origins of hybrid zones between parapatric taxa have been of particular interest for understanding the evolution of reproductive isolation and the geographic context of species divergence. One challenge has been to distinguish between allopatric divergence (followed by secondary contact) versus primary intergradation (parapatric speciation) as alternative divergence histories. Here, we use complementary phylogeographic and population genetic analyses to investigate the recent divergence of two subspecies of Clarkia xantiana and the formation of a hybrid zone within the narrow region of sympatry. We tested alternative phylogeographic models of divergence using approximate Bayesian computation (ABC) and found strong support for a secondary contact model and little support for a model allowing for gene flow throughout the divergence process (i.e. primary intergradation). Two independent methods for inferring the ancestral geography of each subspecies, one based on probabilistic character state reconstructions and the other on palaeo-distribution modelling, also support a model of divergence in allopatry and range expansion leading to secondary contact. The membership of individuals to genetic clusters suggests geographic substructure within each taxon where allopatric and sympatric samples are primarily found in separate clusters. We also observed coincidence and concordance of genetic clines across three types of molecular markers, which suggests that there is a strong barrier to gene flow. Taken together, our results provide evidence for allopatric divergence followed by range expansion leading to secondary contact. The location of refugial populations and the directionality of range expansion are consistent with expectations based on climate change since the last glacial maximum. Our approach also illustrates the utility of combining phylogeographic hypothesis testing with species distribution modelling and fine-scale population genetic analyses for inferring the geography of the divergence process.