When taxonomy meets genomics: lessons from a common songbird

Taxonomy is being increasingly informed by genomics. Traditionally, taxonomy has relied extensively on phenotypic traits for the identification and delimitation of species, though with a growing influence from molecular phylogenetics in recent decades. Now, genomics opens up new and more powerful tools for analysing the evolutionary history and relatedness among species, as well as understanding the genetic basis for phenotypic traits and their role in reproductive isolation. New insights gained from genomics will therefore have major effects on taxonomic classifications and species delimitation. How a genomics approach can inform a flawed taxonomy is nicely exemplified by Mason & Taylor ( 2015 ) in this issue of Molecular Ecology. They studied redpolls, which comprise a genus (Acanthis) of fringillid finches with a wide distribution in the Holarctic region, and whose species taxonomy has been a matter of much controversy for decades (Fig.  1 ). Current authoritative checklists classify them into one, two or three species, and five or six subspecies, based largely on geographical differences in phenotypic traits. Previous studies, including a recent one of the subspecies on Iceland (Amouret et al. 2015 ), have found no evidence of differentiation between these taxa in conventional molecular markers. The lack of genetic structure has been interpreted as incomplete lineage sorting among rapidly evolving lineages. Now Mason & Taylor ( 2015 ), using a large data set of genomewide SNPs, verify that they all belong to a single gene pool with a common evolutionary history, and with little or no geographical structuring. They also show that phenotypic traits used in taxonomic classifications (plumage and bill morphology) are closely associated with polygenic patterns of gene expression, presumably driven by ecological selection on a few regulatory genes. Several lessons can be learned from this study. Perhaps the most important one for taxonomy is the risk of taxonomic inflation resulting from overemphasizing phenotypic traits under local adaptation and ignoring a lack of phylogenetic signal in molecular markers.     

Speciation and phylogeography in the cosmopolitan marine moon jelly, Aurelia sp

Background  

The cosmopolitan moon jelly Aurelia is characterized by high degrees of morphological and ecological plasticity, and subsequently by an unclear taxonomic status. The latter has been revised repeatedly over the last century, dividing the genus Aurelia in as many as 12 or as little as two species. We used molecular data and phenotypic traits to unravel speciation processes and phylogeographic patterns in Aurelia.     

Habitat‐driven diversification,hybridization and cryptic diversity in the Fork‐tailed Drongo (Passeriformes: Dicruridae: Dicrurus adsimilis)

Species complexes of widespread African vertebrates that include taxa distributed across different habitats are poorly understood in terms of their phylogenetic relationships, levels of genetic differentiation and diversification dynamics. The Fork‐tailed Drongo (Dicrurus adsimilis) species complex includes seven Afrotropical taxa with parapatric distributions, each inhabiting a particular bioregion. Various taxonomic hypotheses concerning the species limits of the Fork‐tailed Drongo have been suggested, based largely on mantle and upperpart coloration, but our understanding of diversity and diversification patterns remains incomplete. Especially given our lack of knowledge about how well these characters reflect taxonomy in a morphologically conservative group. Using a thorough sampling across Afrotropical bioregions, we suggest that the number of recognized species within the D. adsimilis superspecies complex has likely been underestimated and that mantle and upperpart coloration reflects local adaptation to different habitat structure, rather than phylogenetic relationships. Our results are consistent with recent phylogeographic studies of sub‐Saharan African vertebrates, indicating that widespread and often morphologically uniform species comprise several paraphyletic lineages, often with one or more of the lineages being closely related to phenotypically distinct forms inhabiting a different, yet geographically close, biome.     

Evolution along the Great Rift Valley: phenotypic and genetic differentiation of East African white‐eyes (Aves,Zosteropidae)

The moist and cool cloud forests of East Africa represent a network of isolated habitats that are separated by dry and warm lowland savannah, offering an opportunity to investigate how strikingly different selective regimes affect species diversification. Here, we used the passerine genus Zosterops (white‐eyes) from this region as our model system. Species of the genus occur in contrasting distribution settings, with geographical mountain isolation driving diversification, and savannah interconnectivity preventing differentiation. We analyze (1) patterns of phenotypic and genetic differentiation in high‐ and lowland species (different distribution settings), (2) investigate the potential effects of natural selection and temporal and spatial isolation (evolutionary drivers), and (3) critically review the taxonomy of this species complex. We found strong phenotypic and genetic differentiation among and within the three focal species, both in the highland species complex and in the lowland taxa. Altitude was a stronger predictor of phenotypic patterns than the current taxonomic classification. We found longitudinal and latitudinal phenotypic gradients for all three species. Furthermore, wing length and body weight were significantly correlated with altitude and habitat type in the highland species Z. poliogaster. Genetic and phenotypic divergence showed contrasting inter‐ and intraspecific structures. We suggest that the evolution of phenotypic characters is mainly driven by natural selection due to differences in the two macro‐habitats, cloud forest and savannah. In contrast, patterns of neutral genetic variation appear to be rather driven by geographical isolation of the respective mountain massifs. Populations of the Z. poliogaster complex, as well as Z. senegalensis and Z. abyssinicus, are not monophyletic based on microsatellite data and have higher levels of intraspecific differentiation compared to the currently accepted species.     

Comparative phylogeographies of six species of hinged terrapins (Pelusios spp.) reveal discordant patterns and unexpected differentiation in the P. castaneus/P. chapini complex and P. rhodesianus

Using up to 2117 bp of mitochondrial DNA and up to 2012 bp of nuclear DNA, we analysed phylogeographic differentiation of six widely distributed species of African hinged terrapins (Pelusios spp.) representing different habitat types. Two taxa each live in savannahs or in forests and mesic savannahs, respectively, and the remaining two species occur in intermediate habitats. The species living in forests and mesic savannahs do not enter dry savannahs, whereas the savannah species may occur in dry and wet savannahs and even in semi‐arid steppe regions. We found no obvious correlation between habitat type and phylogeographic pattern: one savannah species (P. rhodesianus) shows phylogeographic structure, i.e. pronounced genetic differences among geographically distinct populations, and the other (P. nanus) not. One species inhabiting forests and mesic savannahs (P. carinatus) has phylogeographic structure, the other (P. gabonensis) not. The same pattern is true for the two ecologically intermediate species, with phylogeographic structure present in P. castaneus and absent in P. chapini. Nuclear evidence suggests that the latter two taxa with abutting and partially overlapping ranges are distinct, while mtDNA is only weakly differentiated. Pelusios castaneus shows pronounced phylogeographic structure, which could reflect Pleistocene range interruptions correlated with the fluctuating forest cover in West and Central Africa. Our results do not support the recognition of an extinct subspecies of P. castaneus for the Seychelles. Pelusios carinatus contains two well supported clades, which are separated by the Congo River. This species is closely related to P. rhodesianus, a taxon consisting of two deeply divergent mitochondrial clades. One of these clades is paraphyletic with respect to P. carinatus, but the two clades of P. rhodesianus are not differentiated in the studied nuclear markers and, again, paraphyletic with respect to P. carinatus. Using mtDNA sequences from the type material of P. rhodesianus, we were able to allocate this name to one of the two clades. However, owing to the confusing relationships of P. rhodesianus and P. carinatus, we refrain from taxonomic decisions.     

Biogeography and host‐related factors trump parasite life history: limited congruence among the genetic structures of specific ectoparasitic lice and their rodent hosts

Parasites and hosts interact across both micro‐ and macroevolutionary scales where congruence among their phylogeographic and phylogenetic structures may be observed. Within southern Africa, the four‐striped mouse genus, Rhabdomys, is parasitized by the ectoparasitic sucking louse, Polyplax arvicanthis. Molecular data recently suggested the presence of two cryptic species within P. arvicanthis that are sympatrically distributed across the distributions of four putative Rhabdomys species. We tested the hypotheses of phylogeographic congruence and cophylogeny among the two parasite lineages and the four host taxa, utilizing mitochondrial and nuclear sequence data. Despite the documented host‐specificity of P. arvicanthis, limited phylogeographic correspondence and nonsignificant cophylogeny was observed. Instead, the parasite–host evolutionary history is characterized by limited codivergence and several duplication, sorting and host‐switching events. Despite the elevated mutational rates found for P. arvicanthis, the spatial genetic structure was not more pronounced in the parasite lineages compared with the hosts. These findings may be partly attributed to larger effective population sizes of the parasite lineages, the vagility and social behaviour of Rhabdomys, and the lack of host‐specificity observed in areas of host sympatry. Further, the patterns of genetic divergence within parasite and host lineages may also be largely attributed to historical biogeographic changes (expansion‐contraction cycles). It is thus evident that the association between P. arvicanthis and Rhabdomys has been shaped by the synergistic effects of parasite traits, host‐related factors and biogeography over evolutionary time.     

Phylogeographic history of Patagonian lizards of the Liolaemus elongatus complex (Iguania: Liolaemini) based on mitochondrial and nuclear DNA sequences

In this study, we present a phylogeographic analysis of a group of lizards distributed in north‐western Patagonia, the Liolaemus elongatus complex. We sequenced 581 individuals for one mitochondrial gene (cytochrome‐b), and for a subset, we sequenced another mitochondrial gene (12S rRNA) and two nuclear genes: kinesin family member 24 (KIF24) and the anonymous nuclear locus LDAB1D. We estimated gene trees, mitochondrial and nuclear haploytpe networks, standard molecular diversity indices, genetic distances between lineages and Bayesian skyline plots. Our results provide evidence for recognition of seven species previously described within the L. elongatus complex: Liolaemus antumalguen, Liolaemus chillanensis, Liolaemus carlosgarini, Liolaemus burmeisteri, Liolaemus smaug, Liolaemus elongatus and Liolaemus crandalli, but we did not find sufficient evidence to support Liolaemus choique, Liolaemus shitan or Liolaemus sp. 6 as distinct species. We identified four candidate species (Liolaemus sp. 1, Liolaemus sp. 2, Liolaemus sp. 3 and Liolaemus sp. 7), and we discuss evolutionary processes that may have contributed to the origin of these lineages and their taxonomic and conservation implications.     

Unraveling the diversification and systematic puzzle of the highly polymorphic Psammobates tentorius (Bell, 1828) complex (Reptilia: Testudinidae) through phylogenetic analyses and species delimitation approaches

The high level of phenotypic diversity in southern African tent tortoises (Psammobates tentorius complex) has for decades prevented systematists from developing a stable taxonomy for the group. Here, we used a comprehensive DNA sequence dataset (mtDNA: Cytb, ND4, ND4 adjacent tRNA-His, and tRNA-Ser, 12S, 16S; and nDNA: PRLR gene) of 455 specimens, and the latest phylogenetic and species delimitation analytical procedures, to unravel the long-standing P. tentorius complex systematic puzzle. Our results for mtDNA and nDNA were incongruent, with the poorly supported nDNA phylogeny differentiating the three recognized subspecies, and showing potential hybridization in some regions. In contrast, the concatenated mtDNA phylogeny identified seven operational taxonomic units, with strong support. Clades 1, 4, 5, and 7 corresponded to tortoises identified as P. t. tentorius, clade 3 to P. t. trimeni, and clades 2 and 6 to P. t. verroxii. Our analyses showed conflicting topologies for the placement of C6 (P. t. verroxii north of the Orange River), with stronger support for it being sister to C2 + C3 than to the other clades. Clades 1, 2, and 6 had significantly higher genetic diversity than clades 3, 4, 5, and 7, perhaps because these clades inhabit substantially larger areas. The potential for future cladogenic radiations seems high in C1 and C6, particularly in C6 for which the within-clade diversification level was highest. Further research involving microsatellite DNA, phylogeographic evaluations, and morphological variation among clades is crucial for understanding the adaptive radiation of the P. tentorius complex and for modifying their taxonomy.     

Phylogeny of European Pyrochroa (Coleoptera,Pyrochroidae) reveals cryptic taxa and different glacial histories

Only three saproxylic species of Pyrochroinae (Coleoptera: Pyrochroidae) are distributed in Europe, two of which belonging to Pyrochroa: P. coccinea and P. serraticornis. However, P. serraticornis is polytypic, for the presence of the endemic subspecies P. s. kiesenwetteri in southern Italy. Using both molecular and morphological data, we explored the phylogeny of the European Pyrochroa species. A multilocus (COI, CAD, 28S) phylogenetic analysis helped highlight different evolutionary histories for the two examined species. First, P. coccinea, distributed throughout Europe, showed a high differentiation among Italian and European populations. Furthermore, three different taxonomic entities were identified within P. serraticornis, among which the cryptic species Pyrochroa bifoveata sp. n. from central Europe is described and illustrated. A comprehensive identification key to the European Pyrochroinae is also provided. Our results also suggested an historical survival of P. coccinea and P. s. kiesenwetteri in glacial refugia in Italy, and a subsequent post-glacial spread of the former species throughout the Peninsula. In contrast, the current distribution of P. s. serraticornis likely originated from a post-glacial colonization of western European relict populations, while the survival of P. bifoveata plausibly occurred in more eastern glacial refugia (e.g. Carpathian or Balkan regions). Similarly, the European populations of P. coccinea could have originated from relict populations in glacial refugia out from the Italian Peninsula. More comprehensive data on the taxonomy, ecology and biogeography of Pyrochroa are needed to learn more about these species and to help preserve the European saproxylic fauna.     

Comparative phylogeography of codistributed aquatic insects in western North America: insights into dispersal and regional patterns of genetic structure

相似文献   

Genomic,ecological, and morphological approaches to investigating species limits: A case study in modern taxonomy from Tropical Eastern Pacific surgeonfishes

A wide variety of species are distinguished by slight color variations. However, molecular analyses have repeatedly demonstrated that coloration does not always correspond to distinct evolutionary histories between closely related groups, suggesting that this trait is labile and can be misleading for species identification. In the present study, we analyze the evolutionary history of sister species of Prionurus surgeonfishes in the Tropical Eastern Pacific (TEP), which are distinguished by the presence or absence of dark spots on their body. We examined the species limits in this system using comparative specimen‐based approaches, a mitochondrial gene (COI), more than 800 nuclear loci (Ultraconserved Elements), and abiotic niche comparisons. The results indicate there is a complete overlap of meristic counts and morphometric measurements between the two species. Further, we detected multiple individuals with intermediate spotting patterns suggesting that coloration is not diagnostic. Mitochondrial data recovered a single main haplotype shared between the species and all locations resulting in a complete lack of structure (Φ_ST = 0). Genomic analyses also suggest low levels of genetic differentiation (F_ST = 0.013), and no alternatively fixed SNPs were detected between the two phenotypes. Furthermore, niche comparisons could not reject niche equivalency or similarity between the species. These results suggest that these two phenotypes are conspecific and widely distributed in the TEP. Here, we recognize Prionurus punctatus Gill 1862 as a junior subjective synonym of P. laticlavius (Valenciennes 1846). The underlying causes of phenotypic variation in this species are unknown. However, this system gives insight into general evolutionary dynamics within the TEP.     

Nuclear and mitochondrial data reveal different evolutionary processes in the Lake Tanganyika cichlid genus <Emphasis Type="Italic">Tropheus</Emphasis>

Background  

Cichlid fishes are notorious for their wealth of intra- and interspecific colour pattern diversity. In Lake Tanganyika, the endemic genus Tropheus represents the most impressive example for geographic variation in the pattern and hue of integument colouration, but the taxonomy of the over 100 mostly allopatric colour morphs remains to a large degree unresolved. Previous studies of mitochondrial DNA sequence data revealed polyphyly of the six nominally described species and complex phylogeographic patterns influenced by lake level fluctuations and population admixture, and suggested the parallel evolution of similar colour patterns in divergent evolutionary lineages. A gene tree of a rapidly radiating group may be subject to incomplete and stochastic lineage sorting, and to overcome this problem we used multi-locus, nuclear AFLP data in comparison with mtDNA sequences to study diversification, migration and introgression in Tropheus colour morphs in Lake Tanganyika.     

Deep phylogeographic divergence of a migratory passerine in Sino‐Himalayan and Siberian forests: the Red‐flanked Bluetail (Tarsiger cyanurus) complex

Enormous mountainous forests in Sino‐Himalayans and Siberia harbor important avian biodiversity in the Northern Hemisphere. Numerous studies in last two decades have been contributed to systematics and taxonomy of passerines birds in these regions and have revealed various and complex phylogeographic patterns. A passerine species Red‐flanked Bluetail Tarsiger cyanurus provided a good system to manifest such evolutionary complexity. The subspecies T. c. cyanurus and T. c. rufilatus (or/and T. c. pallidior), divergent in morphology, acoustics, and migratory strategies are allopatric in Siberia and Sino‐Himalayan forests, respectively. The two taxa most likely deserve full species status but rigorous genetic analysis is missing. In this study, multilocus phylogeography based on mitochondrial DNA and Z‐linked DNA reveals that T. c. cyanurus and T. c. rufilatus are reciprocally monophyletic with significant statistical support and differ with a large number of diagnostic nucleotide sites resulting substantial genetic divergence. Our finding supports the proposed split of Tarsiger cyanurus s.l. that T. cyanurus and T. rufilatus should be treated as two full species. Whether “pallidior” is a subspecies or geographical form of T. rufilatus is still uncertain. Additionally, these two forest passerine species may have diverged 1.88 (3.25–1.30) Mya, which might be shaped by geographical vicariance due to grassland and desert steppe on the central Loess Plateau during the Pliocene. Taken together, this study and further suggests another independent example of North Palearctic–Sino‐Himalayan phylogeographic pattern in Palearctic birds.     

Phylogeography of Chinese cherry (Prunus pseudocerasus Lindl.) inferred from chloroplast and nuclear DNA: insights into evolutionary patterns and demographic history

Chinese cherry (Prunus pseudocerasus Lindl.) is a commercially valuable fruit crop in China. In order to obtain new insights into its evolutionary history and provide valuable recommendations for resource conservation, phylogeographic patterns of 26 natural populations (305 total individuals) from six geographic regions were analyzed using chloroplast and nuclear DNA fragments. Low levels of haplotype and nucleotide diversity were found in these populations, especially in landrace populations. It is likely that a combined effect of botanical characteristics impact the effective population size, such as inbreeding mating system, long life span, as well as vegetative reproduction. In addition, strong bottleneck effect caused by domestication, together with founder effect after dispersal and subsequent demographic expansion, might also accelerate the reduction of the genetic variation in landrace populations. Interestingly, populations from Longmen Mountain (LMM) and Daliangshan Mountain (DLSM) exhibited relatively higher levels of genetic diversity, inferring the two historical genetic diversity centers of the species. Moreover, moderate population subdivision was also detected by both chloroplast DNA (G_ST = 0.215; N_ST = 0.256) and nuclear DNA (G_ST = 0.146; N_ST = 0.342), respectively. We inferred that the episodes of efficient gene flow through seed dispersal, together with features of long generation cycle and inbreeding mating system, were likely the main contributors causing the observed phylogeographic patterns. Finally, factors that led to the present demographic patterns of populations from these regions and taxonomic varieties were also discussed.     

Phylogeography and diversification of an Amazonian understorey hummingbird: paraphyly and evidence for widespread cryptic speciation in the Plio‐Pleistocene

Straight‐billed Hermit Phaethornis bourcieri inhabits the understorey of upland terra firme forest throughout most of the Amazon basin. Currently, two allopatric taxa regarded as subspecies are recognized: P. b. bourcieri and P. b. major. However, the validity, interspecific limits and evolutionary history of these taxa are not yet fully elucidated. We use molecular characters to propose a phylogenetic hypothesis for populations and taxa grouped under Phaethornis bourcieri. Our results showed that P. bourcieri is part of the ‘Ametrornis’ clade, along with Phaethornis philippii and Phaethornis koepckeae, and that the subspecies major is more closely related to the latter two species than to populations grouped under nominate bourcieri. Our phylogenetic hypotheses recovered three main reciprocally monophyletic clades under nominate bourcieri separated by the lower Negro River and the Branco River or the Branco–Negro interfluve (clades B and C) and the upper Amazon (Solimões) or lower Marañon/Ucayali Rivers (clades C and D). Based on multi‐locus phylogeographic and population genetics approaches, we show that P. b. major is best treated as a separate species, and that P. b. bourcieri probably includes more than one evolutionary species, whose limits remain uncertain. The diversification of the ‘Ametrornis’ clade (P. bourcieri, P. philippii and P. koepckeae) is centred in the Amazon and appears to be closely linked to the formation of the modern Amazon drainage during the Plio‐Pleistocene.     

Phylogeny and systematics of Old World serotine bats (genus Eptesicus,Vespertilionidae, Chiroptera): an integrative approach

Integrative taxonomy aims to document biodiversity by incorporating all useful characters to increase confidence in hypotheses about phylogenetic relationships. In this study, we combine data obtained independently from morphology, two maternally inherited mtDNA genes and two biparentally inherited nuDNA genes to make phylogenetic and taxonomic hypotheses about the Palaearctic members of the bat genus Eptesicus (Vespertilionidae). This genus is distributed worldwide (except for Antarctica) and is highly diversified, presenting one of the most entangled taxonomic puzzles among all mammals. Our results support restoring the genus Rhyneptesicus and separating E. isabellinus and E. pachyomus from E. serotinus and E. ognevi and E. anatolicus from E. bottae. Differences in the phylogenetic hypotheses from mtDNA and nuDNA data suggest the occurrence within E. serotinus of evolutionary processes such as mtDNA capture and secondary contacts between partially differentiated ecomorphs. These two evolutionary processes deserve more in‐depth studies within the group.     

Mitochondrial phylogeny and phylogeography of East African squeaker catfishes (Siluriformes: Synodontis)

Background  

Squeaker catfishes (Pisces, Mochokidae, Synodontis) are widely distributed throughout Africa and inhabit a biogeographic range similar to that of the exceptionally diverse cichlid fishes, including the three East African Great Lakes and their surrounding rivers. Since squeaker catfishes also prefer the same types of habitats as many of the cichlid species, we hypothesized that the East African Synodontis species provide an excellent model group for comparative evolutionary and phylogeographic analyses.     

Structural diversity of nuptial pads in Phyllomedusinae (Amphibia: Anura: Hylidae)

We studied the morphological variation of the nuptial pads using light microscopy and scanning electron microscopy (SEM) in 26 species of phyllomedusines (Anura: Hylidae), representing the five currently recognized genera. All phyllomedusines have single nuptial pads with dark colored epidermal projections (EPs). Spine‐shaped EPs occur in Cruziohyla calcarifer, Phrynomedusa appendiculata and in one species of Phasmahyla. The other species have roundish EPs. The density of the EPs on the pad is variable. Species in the Phyllomedusa hypochondrialis Group have EPs with a density that varies between 764 ± 58/mm² and 923 ± 160/mm². In all other studied species (including the Phyllomedusa burmeisteri and Phyllomedusa perinesos groups, Phyllomedusa camba, Phyllomedusa boliviana, Phyllomedusa sauvagii, Phyllomedusa bicolor, and Phyllomedusa tomopterna) the density of EPs varies between 108 ± 20/mm² and 552 ± 97/mm². Pores were observed with SEM in C. calcarifer, Agalychnis lemur, Agalychnis moreletii, but its presence is confirmed through histological sections on several other species. Its visibility using SEM seems to be related with the level of separation between adjacent EPs. The pores in the four studied species of Agalychnis are shown with SEM and histological sections to have a characteristic epidermal rim, that is absent in the otherphyllomedusines. Unlike most previous reports on breeding glands, those of phyllomedusines are alcian blue positive, indicating the presence of acidic mucosubstances on its secretions. J. Morphol. 2012. © 2012 Wiley Periodicals, Inc.     

Phylogeographic inference using Bayesian model comparison across a fragmented chorus frog species complex

Fragmented species complexes provide an interesting system for investigating biogeographic history and the present distribution of genetic variation. Recent advances in sequencing technology and statistical phylogeography enable the collection and rigorous analysis of large multilocus data sets, but designing studies that produce meaningful phylogeographic inferences remains challenging. We implemented a Bayesian model comparison approach to investigate previous biogeographic hypotheses while simultaneously inferring the presence of genetic structure in a chorus frog species complex. The Illinois chorus frog (Pseudacris illinoensis), originally described as a subspecies of the broadly distributed Strecker's chorus frog (Pseudacris streckeri), occurs in small, disjunct regions associated with scarce sand prairie habitats that have been impacted by human development. We used high‐throughput sequencing to develop and collect a multitiered genetic data set comprised of three different marker types (23 anonymous nuclear sequence loci, four mitochondrial genes and 14 microsatellite loci) designed to address questions across different evolutionary timescales. Phylogenetic analyses uncovered a deep divergence between populations in the Edwards Plateau of central Texas and all other P. streckeri/P. illinoensis populations, but suggest the disjunct distribution of P. illinoensis occurred more recently. Our best‐supported migration model is consistent with the hypothesis that central Texas represented a refugium from which populations expanded via multiple routes. This model also indicates that disjunct northern and southern regions of P. illinoensis should be considered genetically distinct management units. Our study provides an evolutionary context for future studies and conservation efforts in P. illinoensis and demonstrates the utility of model‐based approaches for phylogeographic inference.