MITOCHONDRIAL INTROGRESSION AND INCOMPLETE LINEAGE SORTING THROUGH SPACE AND TIME: PHYLOGENETICS OF CROTAPHYTID LIZARDS

Reconstructing the evolutionary history of the radiation of the land snail genus Xerocrassa on Crete based on mitochondrial sequences and AFLP markers

BACKGROUND

A non-adaptive radiation triggered by sexual selection resulted in ten endemic land snail species of the genus Xerocrassa on Crete. Only five of these species and a more widespread species are monophyletic in a mitochondrial gene tree. The reconstruction of the evolutionary history of such closely related species can be complicated by incomplete lineage sorting, introgression or inadequate taxonomy. To distinguish between the reasons for the nonmonophyly of several species in the mitochondrial gene tree we analysed nuclear AFLP markers.

RESULTS

Whereas six of the eleven morphologically delimited Xerocrassa species from Crete are monophyletic in the mitochondrial gene tree, nine of these species are monophyletic in the tree based on AFLP markers. Only two morphologically delimited species could not be distinguished with the multilocus data and might have diverged very recently or might represent extreme forms of a single species. The nonmonophyly of X. rhithymna with respect to X. kydonia is probably the result of incomplete lineage sorting, because there is no evidence for admixture in the AFLP data and the mitochondrial haplotype groups of these species coalesce deeply. The same is true for the main haplotype groups of X. mesostena. The nonmonophyly of X. franciscoi might be the result of mitochondrial introgression, because the coalescences of the haplotypes of this species with some X. mesostena haplotypes are shallow and there is admixture with neighbouring X. mesostena.

CONCLUSION

The most likely causes for the nonmonophyly of species in the mitochondrial gene tree of the Xerocrassa radiation on Crete could be inferred using AFLP data by a combination of several criteria, namely the depth of the coalescences in the gene tree, the geographical distribution of shared genetic markers, and concordance with results of admixture analyses of nuclear multilocus markers. The strongly subdivided population structure increases the effective population size of land snail species and, thus, the likelihood of a long persistence of ancestral polymorphisms. Our study suggests that ancestral polymorphisms are a frequent cause for nonmonophyly of species with a strongly subdivided population structure in gene trees.

Rampant historical mitochondrial genome introgression between two species of green pond frogs, Pelophylax nigromaculatus and P. plancyi

BACKGROUND

Mitochondrial introgression may result in the mitochondrial genome of one species being replaced by that of another species without leaving any trace of past hybridization in its nuclear genome. Such introgression can confuse the species genealogy estimates and lead to absurd inferences of species history. We used a phylogenetic approach to explore the potential mitochondrial genome introgression event(s) between two closely related green pond frog species, Pelophylax nigromaculatus and P. plancyi.

RESULTS

DNA sequence data of one mitochondrial and two nuclear genes from an extensive sampling of the two species were collected, and the genealogies of the three genes were constructed and compared. While the two nuclear genes congruently showed mutual reciprocal monophyly of both species, the mitochondrial phylogeny separated a Korean P. nigromaculatus clade, a paraphyletic central China P. plancyi assemblage, and a large well-supported introgression clade. Within the introgression clade, the mitochondrial haplotypes of the two species were mixed together. This reticulated pattern can be most parsimoniously explained by an ancient mitochondrial introgression event from P. plancyi to P. nigromaculatus that occurred at least 1.36 MYA, followed by multiple recent introgression events from P. nigromaculatus back to P. plancyi within the last 0.63 MY. The re-constitution of previously co-adapted genomes in P. plancyi may be responsible for the recent rampant introgression events. The Korean P. nigromaculatus clade likely represents the only surviving "true" mitochondrial lineage of P. nigromaculatus, and the central China P. plancyi assemblage likely represents the "original" P. plancyi mitochondrial lineage. Refugia in the Korean Peninsula and central China may have played a significant role in preserving these ancient lineages.

CONCLUSIONS

The majority of individuals in the two species have either introgressed (P. nigromaculatus) or reclaimed (P. plancyi) mitochondrial genomes while no trace of past hybridization in their nuclear genomes was detected. Asymmetrical reproductive ability of hybrids and continuous backcrossing are likely responsible for the observed mitochondrial introgression. This case is unique in that it includes an ancient "forward" introgression and many recent "backward" introgressions, which re-constitutes the original nuclear and mitochondrial genomes of P. plancyi. This hybrid system provides an excellent opportunity to study cyto-nuclear interaction and co-adaptation.

Species delineation using Bayesian model-based assignment tests: a case study using Chinese toad-headed agamas (genus Phrynocephalus)

Background

Species are fundamental units in biology, yet much debate exists surrounding how we should delineate species in nature. Species discovery now requires the use of separate, corroborating datasets to quantify independently evolving lineages and test species criteria. However, the complexity of the speciation process has ushered in a need to infuse studies with new tools capable of aiding in species delineation. We suggest that model-based assignment tests are one such tool. This method circumvents constraints with traditional population genetic analyses and provides a novel means of describing cryptic and complex diversity in natural systems. Using toad-headed agamas of the Phrynocephalus vlangalii complex as a case study, we apply model-based assignment tests to microsatellite DNA data to test whether P. putjatia, a controversial species that closely resembles P. vlangalii morphologically, represents a valid species. Mitochondrial DNA and geographic data are also included to corroborate the assignment test results.

Results

Assignment tests revealed two distinct nuclear DNA clusters with 95% (230/243) of the individuals being assigned to one of the clusters with > 90% probability. The nuclear genomes of the two clusters remained distinct in sympatry, particularly at three syntopic sites, suggesting the existence of reproductive isolation between the identified clusters. In addition, a mitochondrial ND2 gene tree revealed two deeply diverged clades, which were largely congruent with the two nuclear DNA clusters, with a few exceptions. Historical mitochondrial introgression events between the two groups might explain the disagreement between the mitochondrial and nuclear DNA data. The nuclear DNA clusters and mitochondrial clades corresponded nicely to the hypothesized distributions of P. vlangalii and P. putjatia.

Conclusions

These results demonstrate that assignment tests based on microsatellite DNA data can be powerful tools for distinguishing closely related species and support the validity of P. putjatia. Assignment tests have the potential to play a significant role in elucidating biodiversity in the era of DNA data. Nonetheless, important limitations do exist and multiple independent datasets should be used to corroborate results from assignment programs.

Molecular patterns of differentiation in canyon treefrogs (Hyla arenicolor): evidence for introgressive hybridization with the Arizona treefrog (H. wrightorum) and correlations with advertisement call differences

Detection of genetic and behavioural diversity within morphologically similar species has led to the discovery of cryptic species complexes. We tested the hypothesis that US populations of the canyon treefrog (Hyla arenicolor) may consist of cryptic species by examining mate-attraction signals among three divergent clades defined by mtDNA. Using a multi-locus approach, we re-analysed phylogenetic relationships among the three clades and a closely related, but morphologically and behaviourally dissimilar species, the Arizona treefrog (H. wrightorum). We found evidence for introgression of H. wrightorum's mitochondrial genome into H. arenicolor. Additionally, the two-clade topology based on nuclear data is more congruent with patterns of call variation than the three-clade topology from the mitochondrial dataset. The magnitude of the call divergence is probably insufficient to promote isolation of the nuclear DNA-defined clades should they become sympatric, but further divergence in call properties significant in species identification could promote speciation in the future.

Combined mitochondrial 16S and 12S rDNA sequences: an effective genetic marker for inter-species phylogenetic analysis of zoonotic trematodes

The present study studied the genetic variation among Schistosoma japonicum isolates from different endemic regions in mainland China and examined the phylogenetic relationships of zoonotic trematodes using the combined mitochondrial 16S and 12S ribosomal DNA sequences. The fragments of 16S and 12S rDNA were amplified from 22 S. japonicum isolates, and sequenced, and the relevant sequences of other nine trematode species belonging to six genera in four families were downloaded from GenBank, and their phylogenetic relationships were re-constructed by unweighted pair-group method with arithmetic averages analyses using the combined 16S and 12S rDNA sequences, with Trichinella spiralis as outgroup. The results showed that the partial sequences of mitochondrial 16S and 12S rDNA of S. japonicum were 757 and 797 bp, respectively, and they were quite conserved among the S. japonicum isolates. Phylogenetic analysis revealed that the combined 16S and 12S rDNA sequences were not able to distinguish S. japonicum isolates in mountainous areas from those in lake/marshland areas in mainland China. However, the combined sequences could distinguish different species of zoonotic trematodes. Therefore, the combined mitochondrial 16S and 12S rDNA sequences provide an effective molecular marker for the inter-species phylogenetic analysis and differential identification of zoonotic trematodes.

Elucidation of cryptic diversity in a widespread nearctic treefrog reveals episodes of mitochondrial gene capture as frogs diversified across a dynamic landscape

We investigate the evolutionary history of the wide-ranging Nearctic treefrog Hyla arenicolor through the integration of extensive range-wide sampling, phylogenetic analyses of multilocus genetic data, and divergence dating. Previous phylogeographic studies of this frog documented a potential signature of introgressive hybridization from an ecologically and morphologically divergent sister species. Based on our Bayesian phylogenetic analyses of mitochondrial DNA, we inferred strong phylogeographic structure in H. arenicolor as indicated by seven well-supported clades, five of which correspond to well-defined biogeographic regions. Clades from the Balsas Basin and southwestern Central Mexican Plateau in Mexico, and the Grand Canyon of Arizona, group with the morphologically, behaviorally, and ecologically divergent mountain treefrogs in the H. eximia group, rendering H. arenicolor as paraphyletic. The phylogenetic position of at least two of these three H. arenicolor clades within the H. eximia group, however, is most likely the result of several episodes of introgressive hybridization and subsequent mitochondrial gene capture separated in time and space, as supported by evidence from the nuclear genes. Hyla arenicolor from the Balsas Basin appear to be deeply divergent from other H. arenicolor and represent a distinctly different species. Results suggests that introgressive hybridization events, both ancient and contemporary, coupled with late Neogene vicariance and Pleistocene climate-driven range shifts, have all played a role in the historical diversification of H. arenicolor.

Gene flow and species delimitation: a case study of two pine species with overlapping distributions in southeast China

Species delimitation detected by molecular markers is complicated by introgression and incomplete lineage sorting between species. Recent modeling suggests that fixed genetic differences between species are highly related to rates of intraspecific gene flow. However, it remains unclear whether such differences are due to high levels of intraspecific gene flow overriding the spread of introgressed alleles or favoring rapid lineage sorting between species. In pines, chloroplast (cp) and mitochondrial (mt) DNAs are normally paternally and maternally inherited, respectively, and thus their relative rates of intraspecific gene flow are expected to be high and low, respectively. In this study, we used two pine species with overlapping geographical distributions in southeast China, P. massoniana and P. hwangshanensis, as a model system to examine the association between organelle gene flow and variation within and between species. We found that cpDNA variation across these two pine species is more species specific than mtDNA variation and almost delimits taxonomic boundaries. The shared mt/cp DNA genetic variation between species shows no bias in regard to parapatric versus allopatric species' distributions. Our results therefore support the hypothesis that high intraspecific gene flow has accelerated cpDNA lineage sorting between these two pine species.

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.

Molecular phylogeographic analyses of the loach Oxynoemacheilus bureschi reveal post-glacial range extensions across the Balkans

Rivers on the Balkan Peninsula can be separated into ichthyofaunistic areas with different endemic fish species. The Vardar River contains a particularly large number of endemics, indicating its complete and long-term isolation from neighbouring river systems. One of the few species shared with other rivers is the loach species Oxynoemacheilus bureschi. In this study, the genetic analyses of 175 individuals of O. bureschi from 17 sites, covering the entire distribution of the species, including the Rivers Vardar (= Axios), Struma (= Strymon), Mesta (= Nestos) and Danube, were performed using one mitochondrial and one nuclear marker. Genetic differentiation among populations was in general low. Shared haplotypes were common and occurred even between distant localities and different river systems. This points to a high degree of gene flow among populations and rejects the hypothesis that the population in the Vardar River represents a relict from an early colonization of the Balkan Peninsula. In contrast, the results suggest that populations in the Vardar River, as well as those in the Danube River, are of recent origin, and a human-mediated introduction cannot be excluded. On the other hand, the populations in the Aggitis River, a left tributary of the lower Struma River, were clearly separated from the rest of the species and represent a long-term isolated lineage. Demographic analyses suggest a recent population expansion for O. bureschi, in which the population in the Aggitis River was not involved.

Cryptic failure of partitioned Bayesian phylogenetic analyses: lost in the land of long trees

Partitioned Bayesian phylogenetic analyses of routine genetic data sets, constructed using MrBayes (Ronquist and Huelsenbeck 2003), can become trapped in regions of parameter space characterized by unrealistically long trees and distorted partition rate multipliers. Such analyses commonly fail to reach stationarity during hundreds of millions of generations of sampling-many times longer than most published analyses. Some data sets are so prone to this problem that paired MrBayes runs begun from different starting trees repeatedly find the same incorrect long-tree solutions and consequently pass the most commonly employed tests of stationarity, including the average standard deviation of split frequencies (ASDSF) and the potential scale reduction factor (PSRF) statistics offered by MrBayes (Gelman and Rubin 1992). In these situations, failure to reach stationarity is recognizable only in light of prior knowledge of model parameters, such as the expectation that third-codon-position sites usually evolve fastest in protein-coding genes. The conditions that lead to the long-tree problem are frequently encountered in phylogenetic studies today, and I present 6 demonstration examples from the literature. Although the effects on tree length (TL) are often dramatic, effects on topology appear to be subtle. Susceptibility to the problem is sometimes predicted by the difference between the true TL and the starting TL. In some cases, the problems described here can be avoided or reduced by manipulation of the starting TL and/or by adjustments to the prior on branch lengths. In more difficult situations, accurate branch length estimation may not be possible with Bayesian methods because of dependence of the solution on the branch length prior.

Delimiting species boundaries within the Neotropical bamboo Otatea (Poaceae: Bambusoideae) using molecular, morphological and ecological data

Species delimitation is a task that has engaged taxonomists for more than two centuries. Recently, it has been demonstrated that molecular data and ecological niche modeling are useful in species delimitation. In this paper multiple data sets (molecular, morphological, ecological) were utilized to set limits for the species belonging to the Neotropical bamboo Otatea, because there is disagreement about species circumscriptions and also because the genus has an interesting distribution, with most of its populations in Mexico and a single disjunct population in Colombia. Molecular and morphological phylogenetic analyses recovered trees with conflicting topologies. Tree-based morphological and character-based analyses recognized the same entities. Ecological niche models and PCA/MANOVAS agreed with the recognition of the same entities that resulted from the morphological analyses. Morphological analyses retrieved clades supported by diagnostic characters and coherent geographical distributions. Based on these results seven entities should be recognized in Otatea, instead of the three previously described species.

Prokaryotic evolution and the tree of life are two different things

BACKGROUND

The concept of a tree of life is prevalent in the evolutionary literature. It stems from attempting to obtain a grand unified natural system that reflects a recurrent process of species and lineage splittings for all forms of life. Traditionally, the discipline of systematics operates in a similar hierarchy of bifurcating (sometimes multifurcating) categories. The assumption of a universal tree of life hinges upon the process of evolution being tree-like throughout all forms of life and all of biological time. In multicellular eukaryotes, the molecular mechanisms and species-level population genetics of variation do indeed mainly cause a tree-like structure over time. In prokaryotes, they do not. Prokaryotic evolution and the tree of life are two different things, and we need to treat them as such, rather than extrapolating from macroscopic life to prokaryotes. In the following we will consider this circumstance from philosophical, scientific, and epistemological perspectives, surmising that phylogeny opted for a single model as a holdover from the Modern Synthesis of evolution.

RESULTS

It was far easier to envision and defend the concept of a universal tree of life before we had data from genomes. But the belief that prokaryotes are related by such a tree has now become stronger than the data to support it. The monistic concept of a single universal tree of life appears, in the face of genome data, increasingly obsolete. This traditional model to describe evolution is no longer the most scientifically productive position to hold, because of the plurality of evolutionary patterns and mechanisms involved. Forcing a single bifurcating scheme onto prokaryotic evolution disregards the non-tree-like nature of natural variation among prokaryotes and accounts for only a minority of observations from genomes.

CONCLUSION

Prokaryotic evolution and the tree of life are two different things. Hence we will briefly set out alternative models to the tree of life to study their evolution. Ultimately, the plurality of evolutionary patterns and mechanisms involved, such as the discontinuity of the process of evolution across the prokaryote-eukaryote divide, summons forth a pluralistic approach to studying evolution.

REVIEWERS

This article was reviewed by Ford Doolittle, John Logsdon and Nicolas Galtier.

A 2000 km genetic wake yields evidence for northern glacial refugia and hybrid zone movement in a pair of songbirds

Hybrid zones are natural experiments that expose the forces maintaining species differences. But for cases where a trait of one of the hybridizing pair appears shifted into the range of the other, the underlying mechanism can be difficult to infer. For example, hybridization between hermit warbler (Dendroica occidentalis) and Townsend's warbler (Dendroica townsendi) is restricted to narrow hybrid zones in Washington and Oregon, yet hermit mtDNA can be found in phenotypically pure Townsend's populations up to 2000 km north along the Pacific coast. This could reflect introgression of selectively favoured hermit mitochondria north across the hybrid zones, or a neutral genetic wake left behind following southern zone movement. Hermit mitochondrial haplotypes in populations of coastal Townsend's exhibit relatively high genetic diversity and significant divergence from those found in populations of hermit warblers. This contradicts the predictions of selective introgression, but is consistent with a northern population of hermits diverging in a glacial refugium before being replaced by Townsend's via aggressive hybridization. Previous field studies showing Townsend's males to be competitively superior to hermit males support this scenario, and suggest that the extreme hybrid zone movement evidenced by the hermit mitochondrial wake represents an extinction in progress.

Frequent mitochondrial gene introgression among high elevation Tibetan megophryid frogs revealed by conflicting gene genealogies

Historical mitochondrial introgression causes differences between a species' mitochondrial gene genealogy and its nuclear gene genealogy, making tree-based species delineation ambiguous. Using sequence data from one mitochondrial gene (cytochrome b) and three nuclear genes (introns), we examined the evolutionary history of four high elevation Tibetan megophryid frog species, Scutiger boulengeri, Scutiger glandulatus, Scutiger mammatus and Scutiger tuberculatus. The three nuclear genes shared a similar history but the mitochondrial gene tree suggested a drastically different evolutionary scenario. The conflicts between them were explained by multiple episodes of mitochondrial introgression events via historical interspecific hybridization. 'Foreign' mitochondrial genomes might have been fixed in populations and extended through a large portion of the species' distribution. Some hybridization events were probably as old as 10 Myr, while others were recent. An F(1) hybrid was also identified. Historical hybridization events among the four species appeared to be persistent and were not restricted to the period of Pleistocene glaciation, as in several other well-studied cases. Furthermore, hybridization involved several species and occurred in multiple directions, and there was no indication of one mitochondrial genome being superior to others. In addition, incomplete lineage sorting resulting from budding speciation may have also explained some discrepancies between the mitochondrial DNA and nuclear gene trees. Combining all evidences, the former 'Scutiger mammatus' appeared to be two species, including a new species. With the availability of a wide range of highly variable nuclear gene markers, we recommend using a combination of mitochondrial gene and multiple nuclear genes to reveal a complete species history.

Porous genomes and species integrity in the branchiopod Artemia

Over the years, studies on interspecific hybridization have highlighted cases where gene exchange between taxa continues for a significant amount of time after speciation. The reasons for this lag of reproductive isolation relative to genetic isolation are largely unclear, and the question still remains whether the resulting hybrids represent novel biological (and taxonomic) diversity or merely an evolutionary liability. We provide strong indications in the branchiopod Artemia that hybrids between distantly related species may not be evolutionary inconsequential. Based on a global sampling of published and newly derived nuclear (ITS1) and mitochondrial (16S rRNA) sequence data from all representatives of the genus, we have identified natural hybrids between Artemia species (A. persimilisxA. franciscana, A. salinaxA. franciscana) separated by evolutionary interludes of tens of millions of years. Our combined analytical framework of cladistic and network methods provides evidence that hybridizations are the result of recent secondary contact following pronounced allopatric differentiation. The detection of mitochondrial introgression from A. persimilis to A. franciscana attests F(1) hybrid fertility. The reasons for this apparent unidirectionality of introgression are currently unknown but a likely explanation is provided based on morphometric divergence. We discuss the evolutionary implications of our results within the broader context of continental zooplankters.

Phylogenetic relationships among diploid species of Symphyotrichum (Asteraceae: Astereae) based on two nuclear markers, ITS and GAPDH

The mostly North American subtribe Symphyotrichinae (Asteraceae: Astereae) comprises Canadanthus, Ampelaster, Psilactis, Almutaster, and Symphyotrichum. Intergeneric and interspecific relationships within the subtribe have been investigated in the past, particularly by Nesom [Nesom, G.L., 1994. Review of the taxonomy of Aster sensu lato (Asteraceae: Astereae), emphasizing the new world species, Phytologia 77, 141-297] and Semple [Semple, J.C., 2005. Classification of Symphyotrichum. Available from: <http://www.jcsemple.uwaterloo.ca/Symphyotrichumclassification.htm/>], using morphological and cytological approaches. Symphyotrichum is the largest and most complex genus within the subtribe and includes four subgenera: Symphyotrichum (x=7, 8), Virgulus (x=4, 5), Astropolium (x=5), and Chapmaniani (x=7). In this study we used two nuclear markers, the nrDNA internal transcribed spacer (ITS) and the low-copy nuclear gene glyceraldehyde 3-phosphate dehydrogenase (GAPDH), to resolve intergeneric and interspecific relationships within the subtribe at the diploid level, and to determine whether our phylogenies validate the classifications of Nesom or Semple. Our results confirm the distinct generic status of Canadanthus and Ampelaster, whereas Psilactis and Almutaster form a polytomy with Symphyotrichum. Within Symphyotrichum, subg. Virgulus is monophyletic based on ITS and appears polyphyletic based on GAPDH. Neither the ITS nor the GAPDH analyses support a distinct status for subg. Astropolium, which groups within subg. Symphyotrichum. In general, interspecific relationships within Symphyotrichum are unresolved. Lack of resolution may be interpreted as a case of recent and rapid evolutionary radiation.

Pleistocene glacial refugia across the Appalachian Mountains and coastal plain in the millipede genus Narceus: evidence from population genetic, phylogeographic, and paleoclimatic data

Background

Species that are widespread throughout historically glaciated and currently non-glaciated areas provide excellent opportunities to investigate the role of Pleistocene climatic change on the distribution of North American biodiversity. Many studies indicate that northern animal populations exhibit low levels of genetic diversity over geographically widespread areas whereas southern populations exhibit relatively high levels. Recently, paleoclimatic data have been combined with niche-based distribution modeling to locate possible refugia during the Last Glacial Maximum. Using phylogeographic, population, and paleoclimatic data, we show that the distribution and mitochondrial data for the millipede genus Narceus are consistent with classical examples of Pleistocene refugia and subsequent post-glacial population expansion seen in other organismal groups.

Results

The phylogeographic structure of Narceus reveals a complex evolutionary history with signatures of multiple refugia in southeastern North America followed by two major northern expansions. Evidence for refugial populations were found in the southern Appalachian Mountains and in the coastal plain. The northern expansions appear to have radiated from two separate refugia, one from the Gulf Coastal Plain area and the other from the mid-Atlantic coastal region. Distributional models of Narceus during the Last Glacial Maximum show a dramatic reduction from the current distribution, with suitable ecological zones concentrated along the Gulf and Atlantic coastal plain. We found a strong correlation between these zones of ecological suitability inferred from our paleo-model with levels of genetic diversity derived from phylogenetic and population estimates of genetic structuring.

Conclusion

The signature of climatic change, during and after the Pleistocene, on the distribution of the millipede genus Narceus is evident in the genetic data presented. Niche-based historical distribution modeling strengthens the conclusions drawn from the genetic data and proves useful in identifying probable refugia. Such interdisciplinary biogeographic studies provide a comprehensive approach to understanding these processes that generate and maintain biodiversity as well as the framework necessary to explore questions regarding evolutionary diversification of taxa.