DISCORDANT DISTRIBUTION OF POPULATIONS AND GENETIC VARIATION IN A SEA STAR WITH HIGH DISPERSAL POTENTIAL

Patiria miniata, a broadcast‐spawning sea star species with high dispersal potential, has a geographic range in the intertidal zone of the northeast Pacific Ocean from Alaska to California that is characterized by a large range gap in Washington and Oregon. We analyzed spatial genetic variation across the P. miniata range using multilocus sequence data (mtDNA, nuclear introns) and multilocus genotype data (microsatellites). We found a strong phylogeographic break at Queen Charlotte Sound in British Columbia that was not in the location predicted by the geographical distribution of the populations. However, this population genetic discontinuity does correspond to previously described phylogeographic breaks in other species. Northern populations from Alaska and Haida Gwaii were strongly differentiated from all southern populations from Vancouver Island and California. Populations from Vancouver Island and California were undifferentiated with evidence of high gene flow or very recent separation across the range disjunction between them. The surprising and discordant spatial distribution of populations and alleles suggests that historical vicariance (possibly caused by glaciations) and contemporary dispersal barriers (possibly caused by oceanographic conditions) both shape population genetic structure in this species.     

Mitochondrial phylogeographic structure of the white-browed piculet (Sasia ochracea): cryptic genetic differentiation and endemism in Indochina

Aim Our understanding of the geographic patterns of gene flow between populations of birds in the Indo‐Malayan faunal region is surprisingly poor compared with that in other parts of the world. A thorough knowledge of general patterns of phylogeographic structure is, however, of utmost importance for conservation purposes. Species with poor dispersal capabilities could serve as indicators of endemism and genetic isolation in the Indochinese subregion. From their morphology (tiny size, short tail, short and rounded wings), piculets of the genus Sasia are inferred to have poor dispersal capabilities, and thus form a suitable focal species. This study analysed the pattern of genetic variation within the White‐browed Piculet (Sasia ochracea). Location Southeast Asia, north of the Isthmus of Kra. Methods We sampled 43 individuals throughout the breeding range of S. ochracea. DNA was extracted both from fresh tissues (n = 15) and from toe pads from ancient museum skins (n = 28). We amplified a 801‐bp fragment of the mitochondrial ND2 gene to reconstruct the phylogeographic history of the White‐browed Piculet. The sequence data were analysed using Bayesian inference, statistical parsimony, and population genetics methods (analysis of molecular variance, mismatch distributions). We estimated the amount of ongoing gene flow between populations using the coalescent‐based method implemented in Mdiv . Results The analysis of molecular variance indicated that the current taxonomy does not adequately reflect the amount of genetic variation within S. ochracea, as the great majority of genetic variation was nested within the nominal subspecies, which is distributed from Nepal to southern Vietnam. Bayesian inference analyses and haplotype networks suggested the occurrence of five main lineages that are strongly correlated with geography. Our coalescent‐based analyses indicated a very limited amount of ongoing gene flow between these five lineages. Our dating analyses suggested that the genetic structuring probably occurred during the last 400,000 years. Main conclusions Our analyses revealed that S. ochracea is composed of at least five lineages: south Vietnam (South Annam and ‘Cochinchina’), India and Nepal, Myanmar and India, the remainder of Indochina, and probably southern Myanmar (Tenasserim). We strongly recommend that studies aiming to understand the phylogeographic structure within Indo‐Malayan species sample these areas.     

It's about time: divergence, demography, and the evolution of developmental modes in marine invertebrates

Differences in larval developmental mode are predicted to affect ecological and evolutionary processes ranging from gene flow and population bottlenecks to rates of population recovery from anthropogenic disturbance and capacity for local adaptation. The most powerful tests of these predictions use comparisons among species to ask how phylogeographic patterns are correlated with the evolution and loss of prolonged planktonic larval development. An important and largely untested assumption of these studies is that interspecific differences in population genetic structure are mainly caused by differences in dispersal and gene flow (rather than by differences in divergence times among populations or changes in effective population sizes), and that species with similar patterns of spatial genetic variation have similar underlying temporal demographic histories. Teasing apart these temporal and spatial patterns is important for understanding the causes and consequences of evolutionary changes in larval developmental mode. New analytical methods that use the coalescent history of allelic diversity can reveal these temporal patterns, test the strength of traditional population-genetic explanations for variation in spatial structure based on differences in dispersal, and identify strongly supported alternative explanations for spatial structure based on demographic history rather than on gene flow alone. We briefly review some of these recent analytical developments, and show their potential for refining ideas about the correspondence between the evolution of larval developmental mode, population demographic history, and spatial genetic variation.     

Late Pleistocene speciation of three closely related tree peonies endemic to the Qinling–Daba Mountains,a major glacial refugium in Central China

Determining the factors promoting speciation is a major task in ecological and evolutionary research and can be aided by phylogeographic analysis. The Qinling–Daba Mountains (QDM) located in central China form an important geographic barrier between southern subtropical and northern temperate regions, and exhibit complex topography, climatic, and ecological diversity. Surprisingly, few phylogeographic analyses and studies of plant speciation in this region have been conducted. To address this issue, we investigated the genetic divergence and evolutionary histories of three closely related tree peony species (Paeonia qiui, P. jishanensis, and P. rockii) endemic to the QDM. Forty populations of the three tree peony species were genotyped using 22 nuclear simple sequence repeat markers (nSSRs) and three chloroplast DNA sequences to assess genetic structure and phylogenetic relationships, supplemented by morphological characterization and ecological niche modeling (ENM). Morphological and molecular genetic analyses showed the three species to be clearly differentiated from each other. In addition, coalescent analyses using DIYABC conducted on nSSR variation indicated that the species diverged from each other in the late Pleistocene, while ecological niche modeling (ENM) suggested they occupied a larger area during the Last Glacial Maximum (LGM) than at present. The combined genetic evidence from nuclear and chloroplast DNA and the results of ENM indicate that each species persisted through the late Pleistocene in multiple refugia in the Qinling, Daba, and Taihang Mountains with divergence favored by restricted gene flow caused by geographic isolation, ecological divergence, and limited pollen and seed dispersal. Our study contributes to a growing understanding of the origin and population structure of tree peonies and provides insights into the high level of plant endemism present in the Qinling–Daba Mountains of Central China.     

msBayes: Pipeline for testing comparative phylogeographic histories using hierarchical approximate Bayesian computation

Background  

Although testing for simultaneous divergence (vicariance) across different population-pairs that span the same barrier to gene flow is of central importance to evolutionary biology, researchers often equate the gene tree and population/species tree thereby ignoring stochastic coalescent variance in their conclusions of temporal incongruence. In contrast to other available phylogeographic software packages, msBayes is the only one that analyses data from multiple species/population pairs under a hierarchical model.     

Comparative rangewide phylogeography of four endemic Taiwanese bat species

Phylogeographic reconstructions of codistributed taxa can help reveal the interplay between abiotic factors, such as altitude and climate, and species‐specific attributes, in shaping patterns of population genetic structure. Recent studies also demonstrate the value of both rangewide sampling and species distribution modelling (SDM) in comparative phylogeography. Here, we combine these approaches to study the population histories of four phylogenetically related forest‐dependent bat species. All are endemic to the mountainous island of Taiwan but show differences in their tolerance to altitude, with Murina gracilis considered to be a high‐altitude specialist, M. recondita and Kerivoula sp. low‐altitude specialists, and M. puta an altitudinal generalist. We tested the prediction that contrasting habitat preferences would impact on patterns of past and contemporary gene flow and found broad concordance between the results of population genetic analyses and species distribution models based on the Model for Interdisciplinary Research on Climate. Both lowland species showed evidence of genetic divergence between the east and west of the island, consistent with SDMs that indicated the Central Mountain Range (CMR) has presented a long‐term and continuous barrier to gene flow since before the Last Glacial Maximum. In contrast, Murina gracilis and M. puta showed lower degrees of historical isolation and genetic differentiation associated with the CMR, reflecting greater gene flow, possibly coupled with past population growth in M. puta. Together our results highlight the usefulness of combining distribution models with phylogeographic analyses to understand the drivers of genetic structure.     

Amplicon sequencing of 42 nuclear loci supports directional gene flow between South Pacific populations of a hydrothermal vent limpet

In the past few decades, population genetics and phylogeographic studies have improved our knowledge of connectivity and population demography in marine environments. Studies of deep‐sea hydrothermal vent populations have identified barriers to gene flow, hybrid zones, and demographic events, such as historical population expansions and contractions. These deep‐sea studies, however, used few loci, which limit the amount of information they provided for coalescent analysis and thus our ability to confidently test complex population dynamics scenarios. In this study, we investigated population structure, demographic history, and gene flow directionality among four Western Pacific hydrothermal vent populations of the vent limpet Lepetodrilus aff. schrolli. These vent sites are located in the Manus and Lau back‐arc basins, currently of great interest for deep‐sea mineral extraction. A total of 42 loci were sequenced from each individual using high‐throughput amplicon sequencing. Amplicon sequences were analyzed using both genetic variant clustering methods and evolutionary coalescent approaches. Like most previously investigated vent species in the South Pacific, L. aff. schrolli showed no genetic structure within basins but significant differentiation between basins. We inferred significant directional gene flow from Manus Basin to Lau Basin, with low to no gene flow in the opposite direction. This study is one of the very few marine population studies using >10 loci for coalescent analysis and serves as a guide for future marine population studies.     

Nuclear and mitochondrial phylogeography of the European fire-bellied toads Bombina bombina and Bombina variegata supports their independent histories

Exact location and number of glacial refugia still remain unclear for many European cold‐blooded terrestrial vertebrates. We performed a fine‐scaled multilocus phylogeographic analysis of two Bombina species combining mitochondrial variation of 950 toads from 385 sites and nuclear genes (Rag‐1, Ncx‐1) from a subset of samples to reconstruct their colonization and contemporary variation patterns. We identified the lowlands northwest of the Black Sea and the Carpathians to be important refugial areas for B. bombina and B. variegata, respectively. This result emphasizes the importance of Central European refugia for ectothermic terrestrial species, far north of the Mediterranean areas regarded as exclusive glacial refugia for the animals. Additional refugia for B. variegata have been located in the southern Apennines and Balkans. In contrast, no evidence for the importance of other east European plains as refugial regions has been found. The distribution of mtDNA and Ncx‐1 variation suggests the presence of local refugia near the Black Sea for B. bombina; however, coalescent simulations did not allow to distinguish whether one or two refugia were present in the region. Strong genetic drift apparently accompanied postglacial expansions reducing diversity in the colonization areas. Extended sampling, coupled with the multilocus isolation with migration analysis, revealed a limited and geographically restricted gene flow from the Balkan to Carpathian populations of B. variegata. However, despite proximity of inferred B. bombina and B. variegata refugia, gene exchange between them was not detected.     

Defining spatial and temporal patterns of phylogeographic structure in Madagascar's iguanid lizards (genus Oplurus)

Understanding the remarkably high species diversity and levels of endemism found among Madagascar’s flora and fauna has been the focus of many studies. One hypothesis that has received much attention proposes that Quaternary climate fluctuations spurred diversification. While spatial patterns of distribution and phylogenetic relationships can provide support for biogeographic predictions, temporal estimates of divergence are required to determine the fit of these geospatial patterns to climatic or biogeographic mechanisms. We use multilocus DNA sequence data to test whether divergence times among Malagasy iguanid lizards of the subfamily Oplurinae are compatible with a hypotheses of Pliocene–Pleistocene diversification. We estimate the oplurine species tree and associated divergence times under a relaxed‐clock model. In addition, we examine the phylogeographic structure and population divergence times within two sister species of Oplurus primarily distributed in the north‐west and south‐west of Madagascar (Oplurus cuvieri and Oplurus cyclurus, respectively). We find that divergence events among oplurine lineages occurred in the Oligocene and Miocene and are thus far older and incompatible with the hypothesis that recent climate fluctuations are related to current species diversity. However, the timing of intraspecific divergences and spatial patterns of population genetic structure within O. cuvieri and O. cyclurus suggest a role for both intrinsic barriers and recent climate fluctuations at population‐level divergences. Integrating information across spatial and temporal scales allows us to identify and better understand the mechanisms generating patterns diversity.     

Unravelling the evolutionary history and future prospects of endemic species restricted to former glacial refugia

The contemporary distribution and genetic composition of biodiversity bear a signature of species’ evolutionary histories and the effects of past climatic oscillations. For many European species, the Mediterranean peninsulas of Iberia, Italy and the Balkans acted as glacial refugia and the source of range recolonization, and as a result, they contain disproportionately high levels of diversity. As these areas are particularly threatened by future climate change, it is important to understand how past climatic changes affected their biodiversity. We use an integrated approach, combining markers with different evolutionary rates and combining phylogenetic analysis with approximate Bayesian computation and species distribution modelling across temporal scales. We relate phylogeographic processes to patterns of genetic variation in Myotis escalerai, a bat species endemic to the Iberian Peninsula. We found a distinct population structure at the mitochondrial level with a strong geographic signature, indicating lineage divergence into separate glacial refugia within the Iberian refugium. However, microsatellite markers suggest higher levels of gene flow resulting in more limited structure at recent time frames. The evolutionary history of M. escalerai was shaped by the effects of climatic oscillations and changes in forest cover and composition, while its future is threatened by climatically induced range contractions and the role of ecological barriers due to competition interactions in restricting its distribution. This study warns that Mediterranean peninsulas, which provided refuge for European biodiversity during past glaciation events, may become a trap for limited dispersal and ecologically limited endemic species under future climate change, resulting in loss of entire lineages.     

Contrasting support for alternative models of genomic variation based on microhabitat preference: species‐specific effects of climate change in alpine sedges

Deterministic processes may uniquely affect codistributed species’ phylogeographic patterns such that discordant genetic variation among taxa is predicted. Yet, explicitly testing expectations of genomic discordance in a statistical framework remains challenging. Here, we construct spatially and temporally dynamic models to investigate the hypothesized effect of microhabitat preferences on the permeability of glaciated regions to gene flow in two closely related montane species. Utilizing environmental niche models from the Last Glacial Maximum and the present to inform demographic models of changes in habitat suitability over time, we evaluate the relative probabilities of two alternative models using approximate Bayesian computation (ABC) in which glaciated regions are either (i) permeable or (ii) a barrier to gene flow. Results based on the fit of the empirical data to data sets simulated using a spatially explicit coalescent under alternative models indicate that genomic data are consistent with predictions about the hypothesized role of microhabitat in generating discordant patterns of genetic variation among the taxa. Specifically, a model in which glaciated areas acted as a barrier was much more probable based on patterns of genomic variation in Carex nova, a wet‐adapted species. However, in the dry‐adapted Carex chalciolepis, the permeable model was more probable, although the difference in the support of the models was small. This work highlights how statistical inferences can be used to distinguish deterministic processes that are expected to result in discordant genomic patterns among species, including species‐specific responses to climate change.     

Retrospective coalescent methods and the reconstruction of metapopulation histories in the sea

Phylogeographic analyses are a key interface between ecological and evolutionary ways of knowing because such analyses integrate the cumulative effects of demographic (ecological) processes over geological (evolutionary) time scales. Newly developed coalescent methods allow evolutionary ecologists to overcome some limitations associated with inferring population history from classic methods such as Wright’s F _ST. Here we briefly contrast classic and coalescent methods for looking backward in time through a population genetic lens, focusing on the key advantages of the isolation-with-migration (IM) class of coalescent methods for distinguishing ancient connectedness from actual recurrent contemporary gene flow as causes of genetic similarity or differentiation among populations. Making this critical distinction can lead to the discovery of otherwise obscured histories underlying conventional patterns of spatial variation. We illustrate the importance of these insights using analyses of Pacific fishes, snails, and sea stars in which population sizes and divergence times are more important than rates of contemporary gene flow as determinants of population genetic differentiation. We then extend the IM method to genetic data from two model metapopulation species (California abalone, Australian damselfish). The analyses show the potential use of non-equilibrium IM methods for differentiating among metapopulation models that make different predictions about population parameters and have different implications for the design of marine protected areas and other conservation goals. At face value, the results largely rule out classic metapopulation dynamics (dominated by extinction and colonization rather than connectivity via ongoing recurrent gene flow) but, at the same time, do not strongly support a modern marine metapopulation dynamic (ecologically significant connectivity between demes). However, the results also highlight the need for much more data (i.e., loci) sampled on different spatial scales in order to determine whether metapopulation dynamics might exist on smaller scales than are typically sampled by most phylogeographers and landscape geneticists.     

Pleistocene climatic cycling drives intra‐specific diversification in the intermediate horseshoe bat (Rhinolophus affinis) in Southern China

The repeated formation and loss of land‐bridges during the Pleistocene have had lasting impacts on population genetic structure. In the tropics, where island populations persisted through multiple glacial cycles, alternating periods of isolation and contact are expected to have driven population and taxonomic divergence. Here, we combine mitochondrial and nuclear sequence data with microsatellites to dissect the impact of Pleistocene climate change on intra‐specific diversification in the horseshoe bat Rhinolophus affinis. This taxon shows considerable morphological and acoustic variation: two parapatric subspecies (himalayanus and macrurus) occur on mainland China and a third (hainanus) on Hainan Island. Our phylogeographic reconstruction and coalescent analyses suggest the island subspecies formed from an ancestral population of himalayanus via two colonization events c. 800 000 years before present. R. a. hainanus then recolonized the mainland, forming macrurus and thus a secondary contact zone with himalayanus. Finally, macrurus recolonized Hainan following the LGM. We found that all three biological events corresponded to known periods of land‐bridge formation. Evidence of introgression was detected between macrurus and both its sister taxa, with geographical proximity rather than length of separation appearing to be the biggest determinant of subsequent genetic exchange. Our study highlights the important role of climate‐mediated sea level changes have had in shaping current processes and patterns of population structure and taxonomic diversification.     

Evolutionary history of bulldog bats (genus Noctilio): recent diversification and the role of the Caribbean in Neotropical biogeography

We present phylogenetic relationships and phylogeographic patterns of the two species of bulldog bats, genus Noctilio. Using a comprehensive sampling of 118 individuals throughout the species distribution, we investigated the distribution of molecular variation in one nuclear and two mitochondrial markers. Phylogenetic trees do not recover Noctilio albiventris as a monophyletic group and point to three similar‐age intraspecific genetic lineages, suggesting cryptic diversity in this taxon. These lineages correspond to the subspecies previously proposed, and are strongly associated with major river basins in South America. Analyses also suggest a very recent origin for the fishing bat Noctilio leporinus, which probably originates from N. albiventris, with a population expansion corresponding to its invasion in South America. Based on our analysis, the speciation event was dated in the Pleistocene epoch and seems to be associated with the variation of the sea level in the Caribbean islands. The present work indicates how phylogeographic studies support the identification of independent evolutionary lineages, driving new systematic/taxonomic investigations, while at the same time shed light on the role of the Caribbean in shaping Neotropical bat fauna diversity. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Genetic diversity of the imperilled bath sponge Spongia officinalis Linnaeus, 1759 across the Mediterranean Sea: patterns of population differentiation and implications for taxonomy and conservation

The Mediterranean bath sponge Spongia officinalis is an iconic species with high socio‐economic value and precarious future owing to unregulated harvesting, mortality incidents and lack of established knowledge regarding its ecology. This study aims to assess genetic diversity and population structure of the species at different geographical scales throughout its distribution. For this purpose, 11 locations in the Eastern Mediterranean (Aegean Sea), Western Mediterranean (Provence coast) and the Strait of Gibraltar were sampled; specimens were analysed using partial mitochondrial cytochrome oxidase subunit I (COI) sequences, along with a set of eight microsatellite loci. According to our results (i) no genetic differentiation exists among the acknowledged Mediterranean morphotypes, and hence, S. officinalis can be viewed as a single, morphologically variable species; (ii) a notable divergence was recorded in the Gibraltar region, indicating the possible existence of a cryptic species; (iii) restriction to gene flow was evidenced between the Aegean Sea and Provence giving two well‐defined regional clusters, thus suggesting the existence of a phylogeographic break between the two systems; (iv) low levels of genetic structure, not correlated to geographical distance, were observed inside geographical sectors, implying mechanisms (natural or anthropogenic) that enhance dispersal and gene flow have promoted population connectivity; (v) the genetic diversity of S. officinalis is maintained high in most studied locations despite pressure from harvesting and the influence of devastating epidemics. These findings provide a basis towards the effective conservation and management of the species.     

PERCHED AT THE MITO‐NUCLEAR CROSSROADS: DIVERGENT MITOCHONDRIAL LINEAGES CORRELATE WITH ENVIRONMENT IN THE FACE OF ONGOING NUCLEAR GENE FLOW IN AN AUSTRALIAN BIRD

Relationships among multilocus genetic variation, geography, and environment can reveal how evolutionary processes affect genomes. We examined the evolution of an Australian bird, the eastern yellow robin Eopsaltria australis, using mitochondrial (mtDNA) and nuclear (nDNA) genetic markers, and bioclimatic variables. In southeastern Australia, two divergent mtDNA lineages occur east and west of the Great Dividing Range, perpendicular to latitudinal nDNA structure. We evaluated alternative scenarios to explain this striking discordance in landscape genetic patterning. Stochastic mtDNA lineage sorting can be rejected because the mtDNA lineages are essentially distinct geographically for > 1500 km. Vicariance is unlikely: the Great Dividing Range is neither a current barrier nor was it at the Last Glacial Maximum according to species distribution modeling; nuclear gene flow inferred from coalescent analysis affirms this. Female philopatry contradicts known female‐biased dispersal. Contrasting mtDNA and nDNA demographies indicate their evolutionary histories are decoupled. Distance‐based redundancy analysis, in which environmental temperatures explain mtDNA variance above that explained by geographic position and isolation‐by‐distance, favors a nonneutral explanation for mitochondrial phylogeographic patterning. Thus, observed mito‐nuclear discordance accords with environmental selection on a female‐linked trait, such as mtDNA, mtDNA–nDNA interactions or genes on W‐chromosome, driving mitochondrial divergence in the presence of nuclear gene flow.     

Ancient differentiation in the single-island avian radiation of endemic Hispaniolan chat-tanagers (Aves: Calyptophilus)

The simple geographic structure of island systems often makes them tractable for studies of the patterns and processes of biological diversification. The Calyptophilus chat-tanagers of Hispaniola are of general evolutionary interest because their multiple lineages might have arisen on a single island, of conservation concern because several isolated populations are nearly extinct, and taxonomically ambiguous because they have been variously lumped or split into one to four species. To explore the context of diversification of the seven extant Calyptophilus populations, we conducted a multilocus coalescent analysis based on sequences of mitochondrial ND2 and three nuclear intron loci. We then compared patterns of phylogeographic genetic variation with the morphological differences that distinguish these populations. Mitochondrial haplotypes formed two reciprocally monophyletic groups separated by a large magnitude of nucleotide divergence. Intron structure largely paralleled the geographic grouping pattern of the mitochondrial DNA (mtDNA), but these groups were only reciprocally monophyletic at one of the three introns. Also, the magnitude of between-group divergence was much lower in the introns than mtDNA genealogies. Multilocus coalescent analyses inferred a nonzero divergence time between these two major geographic groups, but suggested that they have experienced a low level of gene flow. All four markers showed substantial allele sharing within each of the two groups, demonstrating that many now separated montane populations do not have long histories of isolation. Considered in concert, our multilocus phylogeographic reconstructions support the recognition of two species within the Calyptophilus complex, and raise the possibility that these taxa differentiated prior to the fusion of the two palaeo-islands that form present-day Hispaniola.     

Historical male‐mediated introgression in horseshoe bats revealed by multilocus DNA sequence data

Instances of hybridization between mammalian taxa in the wild are rarely documented. To test for introgression between sibling species of horseshoe bat (Rhinolophus yunanensis and R. pearsoni) and two subspecies of the latter (R. p. pearsoni and R. p. chinensis), we sequenced two mtDNA and two ncDNA markers in individuals sampled from multiple localities within their overlapping ranges. The interspecific mtDNA gene tree corresponded to the expected taxonomic divisions, and coalescent‐based analyses suggested divergence occurred around 4 MYA. However, these relationships strongly conflicted with those recovered from two independent nuclear gene trees, in which R. yunanensis clustered with R. p. pearsoni to the exclusion of R. p. chinensis. This geographically widespread discordance is best explained by large‐scale historical introgression of ncDNA from R. yunanensis to R. pearsoni by male‐mediated exchange in mixed species colonies during Pleistocene glacial periods, when ranges may have contracted and overlapped more than at present. Further species tree–gene tree conflicts were detected between R. p. pearsoni and R. p. chinensis, also indicating past and/or current introgression in their overlapping regions. However, here the patterns point to asymmetric mtDNA introgression without ncDNA introgression. Analyses of coalescence times indicate this exchange has occurred subsequent to the divergence of these subspecies from their common ancestor. Our work highlights the importance of using multiple data sets for reconstructing phylogeographic histories and resolving taxonomic relationships.