A re‐examination of the molecular phylogeny and biogeography of the genus Schizothorax (Teleostei: Cyprinidae) through enhanced sampling,with emphasis on the species in the Yunnan–Guizhou Plateau,China

The mitochondrial DNA cytochrome b sequences of 36 Schizothorax species from 51 localities in the Yunnan–Guizhou Plateau (YGP) and its adjacent areas were analysed. Maximum parsimony, Maximum likelihood and Bayesian phylogenetic analyses were performed to examine the relationships of Schizothorax species. A hypothesis of the phylogenetic relationships of the species is given. A relaxed molecular clock based on Bayesian evolutionary analysis was used to tentatively calculate the divergence times of Schizothorax. Samples from the YGP were tentatively grouped into three geographically distributed clades: the Tsangpo‐Irrawaddy, the Mekong‐Salween and the Trans‐Jinsha River (including Jinsha, Red, Nanpan and Beipan Rivers). Calibration of the molecular clock revealed that two geological periods, the late Miocene about 10 million years before present (Myr BP) and the Pliocene (4.0 Myr BP), were important times in the vicariant speciation of Schizothorax. The phylogenetic history of the species is congruent with events caused by the uplift of the Tibetan Plateau and the YGP. The divergence of Schizothorax species in YGP began in the Pliocene. Our phylogenetic trees did not support the hypothesis that the paleo Jinsha River was drained through the Yangtze River‐Jianchuan Lake‐Erhai Lake to the Red River. Schizothorax in the Beipan River were derived from the Jinsha River.     

Out of the forest: past and present range expansion of a parthenogenetic weevil pest,or how to colonize the world successfully

Previous research revealed complex diversification patterns in the parthenogenetic weevil Naupactus cervinus. To understand the origin of clonal diversity and successful spreading of this weevil, we investigated its geographic origin and possible dispersal routes and whether parthenogens can persist in habitats under unsuitable environmental conditions. This study is based on samples taken throughout a broad area of the species’ range. We used both mitochondrial and nuclear markers and applied phylogenetic and network analyses to infer possible relationships between haplotypes. Bayesian phylogeographic analyses and ecological niche modeling were used to investigate the processes that shaped genetic diversity and enabled the colonization of new geographic areas. Southeastern Brazil emerges as the original distribution area of N. cervinus. We detected two range expansions, one along natural corridors during the Pleistocene and the other in countries outside South America during recent times. Isolation due to climate shifts during the early Pleistocene led to diversification in two divergent clades, which probably survived in different refugia of the Paranaense Forest and the Paraná River delta. The origin of the clonal diversity was probably a complex process including mutational diversification, hybridization, and secondary colonization. The establishment of N. cervinus in areas outside its native range may indicate adaptation to drier and cooler conditions. Parthenogenesis would be advantageous for the colonization of new environments by preventing the breakup of successful gene combinations. As in other insect pests, the present distribution of N. cervinus results from both its evolutionary history and its recent history related to human activities.     

Allopatric divergence and secondary contact without genetic admixture for Arichanna perimelaina (Lepidoptera: Geometridae), an alpine moth endemic to the Hengduan Mountains

Mountain systems, especially at high altitudes, are an excellent model for determining the mechanisms underlying high species diversity and endemism. Herein, we elucidate the evolutionary history of the alpine moth Arichanna perimelaina (Wehrli), which is endemic to the Hengduan Mountains (HM) region in southwest China, based on three mitochondrial genes and two nuclear genes. Our results revealed six deeply divergent clades that corresponded to populations in different mountain systems in the HM region. Bayesian divergence time estimations suggested a mid‐ to late Pleistocene genetic divergence. The results also showed that the Mt Yulong (YL) region was a refugium and valley corridors established by glaciation during the Pleistocene allowed populations on the separate mountains to migrate. The reproductive isolation among the different clades on contact zone in the YL region may be associated with the asynchronous mating rhythms and/or the divergent mate recognition caused by the ecological source of divergent selection. Allopatric divergence associated with complex topographies and climatic oscillations, regional dispersal via valley corridors and the suitable refugium of the YL region shaped the genetic divergence and distribution pattern of A. perimelaina in the HM region. These findings highlight the essential role of complex terrain and climatic fluctuations in shaping the unique phylogeographic history of a narrow alpine moth, and provide insights into the mechanisms underlying high species richness and endemism in the HM region.     

Population structure of the Indonesian giant tiger shrimp Penaeus monodon: a window into evolutionary similarities between paralogous mitochondrial DNA sequences and their genomes

Here we used both microsatellites and mtCR (mitochondrial DNA control region) sequences as genetic markers to examine the genetic diversity and population structure of Penaeus monodon shrimp from six Indonesian regions. The microsatellite data showed that shrimp from the Indian and the Pacific Ocean were genetically distinct from each other. It has been reported previously that P. monodon mtCR sequences from the Indo‐Pacific group into two major paralogous clades of unclear origin. Here we show that the population structure inferred from mtCR sequences matches the microsatellite‐based population structure for one of these clades. This is consistent with the notion that this mtCR clade shares evolutionary history with nuclear DNA and may thus represent nuclear mitochondrial pseudogenes (Numts).     

Pliocene intraspecific divergence and Plio‐Pleistocene range expansions within Picea likiangensis (Lijiang spruce), a dominant forest tree of the Qinghai‐Tibet Plateau

A knowledge of intraspecific divergence and range dynamics of dominant forest trees in response to past geological and climate change is of major importance to an understanding of their recent evolution and demography. Such knowledge is informative of how forests were affected by environmental factors in the past and may provide pointers to their response to future environmental change. However, genetic signatures of such historical events are often weak at individual loci due to large effective population sizes and long generation times of forest trees. This problem can be overcome by analysing genetic variation across multiple loci. We used this approach to examine intraspecific divergence and past range dynamics in the conifer Picea likiangensis, a dominant tree of forests occurring in eastern and southern areas of the Qinghai‐Tibet Plateau (QTP). We sequenced 13 nuclear loci, two mitochondrial DNA regions and three plastid (chloroplast) DNA regions in 177 individuals sampled from 22 natural populations of this species, and tested the hypothesis that its evolutionary history was markedly affected by Pliocene QTP uplifts and Quaternary climatic oscillations. Consistent with the taxonomic delimitation of the three morphologically divergent varieties examined, all individuals clustered into three genetic groups with intervariety admixture detected in regions of geographical overlap. Divergence between varieties was estimated to have occurred within the Pliocene and ecological niche modelling based on 20 ecological variables suggested that niche differentiation was high. Furthermore, modelling of population‐genetic data indicated that two of the varieties (var. rubescens and var. linzhiensis) expanded their population sizes after the largest Quaternary glaciation in the QTP, while expansion of the third variety (var. likiangensis) began prior to this, probably following the Pliocene QTP uplift. These findings point to the importance of geological and climatic changes during the Pliocene and Pleistocene as causes of intraspecific diversification and range shifts of dominant tree species in the QTP biodiversity hot spot region.     

Phylogeny and taxonomy of the genus Ilyodon Eigenmann, 1907 (Teleostei: Goodeidae), based on mitochondrial and nuclear DNA sequences

Taxonomy of the live‐bearing fish of the genus Ilyodon Eigenmann, 1907 (Goodeidae), in Mexico, is controversial, with morphology and mitochondrial genetic analyses in disagreement about the number of valid species. The present study accumulated a comprehensive DNA sequences dataset of 98 individuals of all Ilyodon species and mitochondrial and nuclear loci to reconstruct the evolutionary history of the genus. Phylogenetic inference produced five clades, one with two sub‐clades, and one clade including three recognized species. Genetic distances in mitochondrial genes (cytb: 0.5%–2.1%; coxI: 0.5%–1.1% and d‐loop: 2.3%–10.2%) were relatively high among main clades, while, as expected, nuclear genes showed low variation (0.0%–0.2%), with geographic concordance and few shared haplotypes among river basins. High genetic structure was observed among clades and within basins. Our genetic analyses, applying the priority principle, suggest the recognition only of Ilyodon whitei and Ilyodon furcidens, with I. cortesae relegated to an invalid species, the populations of which belong to I. whitei.     

High levels of population subdivision in a morphologically conserved Mediterranean toad (Alytes cisternasii) result from recent,multiple refugia: evidence from mtDNA,microsatellites and nuclear genealogies

Pleistocene glaciations often resulted in differentiation of taxa in southern European peninsulas, producing the high levels of endemism characteristic of these regions (e.g. the Iberian Peninsula). Despite their small ranges, endemic species often exhibit high levels of intraspecific differentiation as a result of a complex evolutionary history dominated by successive cycles of fragmentation, expansion and subsequent admixture of populations. Most evidence so far has come from the study of species with an Atlantic distribution in northwestern Iberia, and taxa restricted to Mediterranean‐type habitats remain poorly studied. The Iberian Midwife toad (Alytes cisternasii) is a morphologically conserved species endemic to southwestern and central Iberia and a typical inhabitant of Mediterranean habitats. Applying highly variable genetic markers from both mitochondrial and nuclear genomes to samples collected across the species’ range, we found evidence of high population subdivision within A. cisternasii. Mitochondrial haplotypes and microsatellites show geographically concordant patterns of genetic diversity, suggesting population fragmentation into several refugia during Pleistocene glaciations followed by subsequent events of geographical and demographic expansions with secondary contact. In addition, the absence of variation at the nuclear β‐fibint7 and Ppp3caint4 gene fragments suggests that populations of A. cisternasii have been recurrently affected by episodes of extinction and recolonization, and that documented patterns of population subdivision are the outcome of recent and multiple refugia. We discuss the evolutionary history of the species with particular interest in the increasing relevance of Mediterranean refugia for the survival of genetically differentiated populations during the Pleistocene glaciations as revealed by studies in co‐distributed taxa.     

History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: Aconitum gymnandrum (Ranunculaceae)

How Quaternary climatic oscillations affected range distributions and intraspecific divergence of alpine plants on the Qinghai‐Tibetan Plateau (QTP) remains largely unknown. Here, we report a survey of chloroplast DNA (cpDNA) and nuclear ribosomal internal transcribed spacer (ITS) DNA variation aimed at exploring the phylogeographical history of the QTP alpine endemic Aconitum gymnandrum. We sequenced three cpDNA fragments (rpl20–rps12 intergenic spacer, the trnV intron and psbA‐trnH spacer) and also the nuclear (ITS) region in 245 individuals from 23 populations sampled throughout the species’ range. Two distinct lineages, with eastern and western geographical distributions respectively, were identified from a phylogenetic analysis of ITS sequence variation. Based on a fast substitution rate, these were estimated to have diverged from each other in the early Pleistocene approximately 1.45 Ma. The analysis of cpDNA variation identified nine chlorotypes that clustered into two major clades that were broadly congruent in geographical distribution with the two ITS lineages. The east–west split of cpDNA divergence was supported by an amova which partitioned approximately half of the total variance between these two groups of populations. Analysis of the spatial distribution of chlorotypes showed that each clade was subdivided into two groups of populations such that a total of four population groups existed in the species. It is suggested that these different groups derive from four independent glacial refugia that existed during the Last Glacial Maximum (LGM), and that three of these refugia were located at high altitude on the QTP platform itself at that time. Coalescent simulation of chlorotype genealogies supported both an early Pleistocene origin of the two main cpDNA clades and also the ‘four‐refugia’ hypothesis during the LGM. Two previous phylogeographical studies of QTP alpine plants indicated that such plants retreated to refugia at the eastern/south‐eastern plateau edge during the LGM and/or previous glacial maxima. However, the results for A. gymnandrum suggest that at least some of these cold‐tolerant species may have also survived centrally on the QTP platform throughout the Quaternary.     

Contrasting responses to Pleistocene climate changes: a case study of two sister species Allium cyathophorum and A. spicata (Amaryllidaceae) distributed in the eastern and western Qinghai–Tibet Plateau

It has been hypothesized that species occurring in the eastern and the western Qinghai–Tibet Plateau (QTP) responded differently to climate changes during the Pleistocene. Here, we test this hypothesis by phylogeographic analysis of two sister species, Allium cyathophorum and A. spicata. We sequenced two chloroplast DNA (cpDNA) fragments (accD‐psaI and the rpl16 intron) of 150 individuals, and the nuclear (ITS) region of 114 individuals, from 19 populations throughout the distributional ranges of these species. The divergence between the two species was dated at 779 ‐ 714 thousand years before the present and was likely initiated by the most major glaciation in the QTP. Analysis of chlorotype diversity showed that A. spicata, the species occurring in the western QTP, contains much lower genetic diversity (0.25) than A. cyathophorum (0.93), which is distributed in the eastern QTP. Moreover, multiple independent tests suggested that the A. spicata population had expanded recently, while no such expansion was detected in A. cyathophorum, indicating a contrasting pattern of responses to Pleistocene climate changes. These findings highlight the importance of geographical topography in determining how species responded to the climate changes that took place in the QTP during the Pleistocene.     

High genetic variation and phylogeographic relations among Palearctic fairy shrimp populations reflect persistence in multiple southern refugia during Pleistocene ice ages and postglacial colonisation

相似文献   

River islands,refugia and genetic structuring in the endemic brown frog Rana kukunoris (Anura,Ranidae) of the Qinghai‐Tibetan Plateau

Frequently, Pleistocene climatic cycling has been found to be the diver of genetic structuring in populations, even in areas that did not have continental ice sheets, such as on the Qinghai‐Tibetan Plateau (QTP). Typically, species distributed on the plateau have been hypothesized to re‐treat to south‐eastern refugia, especially during the Last Glacial Maximum (LGM). We evaluated sequence variation in the mitochondrial DNA gene Cytb and the nuclear DNA gene RAG‐1 in Rana kukunoris, a species endemic to the QTP. Two major lineages, N and S, were identified, and lineage N was further subdivided into N1 and N2. The geographical distribution and genealogical divergences supported the hypothesis of multiple refugia. However, major lineages and sublineages diverged prior to the LGM. Demographical expansion was detected only in lineage S and sublineage N2. Sublineage N1 might have survived several glacial cycles in situ and did not expand after the LGM because of the absence of suitable habitat; it survived in river islands. Genetic analysis and environment modelling suggested that the north‐eastern edge of QTP contained a major refugium for R. kukunoris. From here, lineage S dispersed southwards after the LGM. Two microrefugia in northern Qilian Mountains greatly contributed to current level of intraspecific genetic diversity. These results were found to have important implications for the habitat conservation in Northwest China.     

The Pleistocene demography of an alpine juniper of the Qinghai‐Tibetan Plateau: tabula rasa,cryptic refugia or something else?

Aim Numerous palaeoecological and genetic studies have shown that different tree species responded in very different ways to Pleistocene climatic oscillations. Some were forced into small refugia far from their current range, while others were able to survive in small refugia close to, or even within, their current natural range. In this study we examine the Pleistocene demography of a juniper species (Juniperus przewalskii, Cupressaceae) from the Qinghai‐Tibetan Plateau. Location The Qinghai‐Tibetan Plateau (QTP). Methods Eight nuclear loci were sequenced in 141 individuals from 20 natural populations distributed across the entire natural range of J. przewalskii, and coalescent analysis was used to test demographic hypotheses. Results The overall nucleotide diversity in the sample was low (π_sil = 0.0029), with few rare alleles and pronounced population genetic structure (F_ST = 0.181). We detected a division previously found using chloroplast DNA markers: all segregating sites in populations from the central part of the QTP appear to be a subset of those found around the edge of the plateau, confirming the relatively young age of the former. In contrast to the middle Pleistocene bottlenecks detected in boreal tree species, the coalescent‐based analyses failed to reject the standard neutral model for the juniper species considered here. Main conclusions Juniperus przewalskii did not undergo marked changes in population sizes during the Pleistocene, although this species seems to have experienced recent, post‐glacial expansion. This finding is largely consistent with the limited number of previous studies on conifer species of the QTP, but contradicts findings of studies on boreal species. These findings have wide implications for understanding plant species’ responses to past climatic oscillations on the high‐elevation QTP.     

Pleistocene climate changes shaped the divergence and demography of Asian populations of the great tit Parus major: evidence from phylogeographic analysis and ecological niche models

Different scales and frequencies of glaciations developed in Europe and Asia during the Pleistocene. Because species’ responses to climate change are influenced by interactive factors including ecology and local topography, the pattern and tempo of species diversification may vary significantly across regions. The great tit Parus major is a widespread Eurasian passerine with a range that encircles the central Asian desert and high‐altitude areas of the Tibetan Plateau. A number of genetic studies have assessed the effect of paleo‐climate changes on the distribution of the European population. However, none have comprehensively addressed how paleo‐climate change affected the distribution of the great tit in China, an apparent hotspot of P. major subspecific diversity. Here, we describe likely paleo‐climatic effects on P. major populations in China based on a combination of phylogeography and ecological niche models (ENMs). We sequenced three mitochondrial DNA markers from 28 populations (213 individuals), and downloaded 112 sequences from outside its Chinese range. As the first step in clarifying the intra‐specific relationships among haplotypes, we attempted to clarify the divergence and demography of populations in China. Phylogeographic analysis revealed that P. major is comprised of five highly divergent clades with geographic breaks corresponding to steep mountains and dry deserts. A previously undescribed monophyletic clade with high genetic diversity, stable niches and a long and independent evolutionary history was detected in the mountainous areas of southwest China. The estimated times at which these clades diverged was traced back to the Early‐Middle Pleistocene (2.19–0.61 mya). Contrary to the post‐LGM (the Last Glacial Maximum) expansion of European populations, demographic history indicates that Asian populations expanded before the LGM after which they remained relatively stable or grew slowly through the LGM. ENMs support this conclusion and predict a similar distribution in the present and the LGM. Our genetic and ecological results demonstrate that Pleistocene climate changes shaped the divergence and demography of P. major in China.     

Diversification in wild populations of the model organism Anolis carolinensis: A genome‐wide phylogeographic investigation

The green anole (Anolis carolinensis) is a lizard widespread throughout the southeastern United States and is a model organism for the study of reproductive behavior, physiology, neural biology, and genomics. Previous phylogeographic studies of A. carolinensis using mitochondrial DNA and small numbers of nuclear loci identified conflicting and poorly supported relationships among geographically structured clades; these inconsistencies preclude confident use of A. carolinensis evolutionary history in association with morphological, physiological, or reproductive biology studies among sampling localities and necessitate increased effort to resolve evolutionary relationships among natural populations. Here, we used anchored hybrid enrichment of hundreds of genetic markers across the genome of A. carolinensis and identified five strongly supported phylogeographic groups. Using multiple analyses, we produced a fully resolved species tree, investigated relative support for each lineage across all gene trees, and identified mito‐nuclear discordance when comparing our results to previous studies. We found fixed differences in only one clade—southern Florida restricted to the Everglades region—while most polymorphisms were shared between lineages. The southern Florida group likely diverged from other populations during the Pliocene, with all other diversification during the Pleistocene. Multiple lines of support, including phylogenetic relationships, a latitudinal gradient in genetic diversity, and relatively more stable long‐term population sizes in southern phylogeographic groups, indicate that diversification in A. carolinensis occurred northward from southern Florida.     

Evolutionary history of the Maltese wall lizard Podarcis filfolensis: insights on the ‘Expansion–Contraction’ model of Pleistocene biogeography

The expansion–contraction (EC) model predicts demographic and range contraction of temperate species during Pleistocene glaciations as a consequence of climate‐related habitat changes, and provides a paradigm for explaining the high intraspecific diversity found in refugia in terms of long‐term demographic stability. However, recent evidence has revealed a weak predictive power of this model for terrestrial species in insular and coastal settings. We investigated the Pleistocene EC dynamics and their evolutionary consequences on temperate species using the Maltese archipelago and its endemic lizard Podarcis filfolensis as a model system. The evolutionary and demographic history of P. filfolensis as inferred from mitochondrial and nuclear sequences data does not conform to the EC model predictions, supporting (i) demographic and spatial stability or expansion, rather than contraction, of the northern and southern lineages during the last glacial period; and (ii) a major role for allopatric differentiation primed by sea‐level dynamics, rather than prolonged demographic stability, in the formation of the observed genetic diversity. When combined with evidence from other Mediterranean refugia, this study shows how the incorporation of Pleistocene sea‐level variations in the EC model accounts for a reverse demographic and range response of insular and coastal temperate biotas relative to continental ones. Furthermore, this cross‐archipelago pattern in which allopatric diversity is formed and shaped by EC cycles resembles that seen between isolated populations within mainland refugia and suggests that the EC model, originally developed to explain population fluctuations into and out‐of refugia, may be appropriate for describing the demographic and evolutionary dynamics driving the high genetic diversity observed in these areas.     

Phylogeography of the genus Dasiphora (Rosaceae) in the Qinghai‐Tibetan Plateau: divergence blurred by expansion

Plateau uprisings and climatic oscillations are considered to have caused extensive allopatric divergences that account for the rich species diversity of the Qinghai‐Tibetan Plateau (QTP). However, secondary contact during range shifts in the Quaternary glacial cycles or inter‐uplift stages may have restored the gene flow between species and so counteracted these divergences, particularly in rapidly‐adapting dominant elements. We tested this hypothesis by determining the phylogeographical history of Dasiphora (Rosaceae), a genus of two species that are widely distributed on the QTP and co‐exist in numerous localities. We sequenced two chloroplast DNA fragments (rbcL, trnT‐L) for 559 individuals from 87 populations. Bayesian methods were used to identify phylogenetic relationships and to estimate divergence times. Demographic histories were inferred using neutrality tests, mismatch distribution analysis, and coalescent simulation. A total of 112 haplotypes that clustered into three major groups were identified. The formation of these groups and their subgroups was dated to between the Pliocene and the late Pleistocene. In addition, we found that some groups underwent multiple extensive expansions. Species‐specific haplotypes were identified for each species, although these haplotypes phylogenetically intermixed. These results suggest that recent plateau uplifts and climatic oscillations might have caused the deep divergences observed within this genus. However, later range expansions probably blurred these divergences and possible species boundaries. Our results shed new light on the complex evolutionary history of the QTP alpine plants. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 777–788.     

DEEP DIVERSIFICATION AND LONG‐TERM PERSISTENCE IN THE SOUTH AMERICAN ‘DRY DIAGONAL’: INTEGRATING CONTINENT‐WIDE PHYLOGEOGRAPHY AND DISTRIBUTION MODELING OF GECKOS

The relative influence of Neogene geomorphological events and Quaternary climatic changes as causal mechanisms on Neotropical diversification remains largely speculative, as most divergence timing inferences are based on a single locus and have limited taxonomic or geographic sampling. To investigate these influences, we use a multilocus (two mitochondrial and 11 nuclear genes) range‐wide sampling of Phyllopezus pollicaris, a gecko complex widely distributed across the poorly studied South American ‘dry diagonal’ biomes. Our approach couples traditional and model‐based phylogeography with geospatial methods, and demonstrates Miocene diversification and limited influence of Pleistocene climatic fluctuations on P. pollicaris. Phylogeographic structure and distribution models highlight that persistence across multiple isolated regions shaped the diversification of this species complex. Approximate Bayesian computation supports hypotheses of allopatric and ecological/sympatric speciation between lineages that largely coincide with genetic clusters associated with Chaco, Cerrado, and Caatinga, standing for complex diversification between the ‘dry diagonal’ biomes. We recover extremely high genetic diversity and suggest that eight well‐supported clades may be valid species, with direct implications for taxonomy and conservation assessments. These patterns exemplify how low‐vagility species complexes, characterized by strong genetic structure and pre‐Pleistocene divergence histories, represent ideal radiations to investigate broad biogeographic histories of associated biomes.     

Phylogeography of the Petunia integrifolia complex in southern Brazil

In this study, we evaluated the genetic diversity of the Petunia integrifolia species group using a phylogeographical approach, and attempted to understand better its diversification and taxonomy. Plants from five morphological groups were collected, covering a large part of the geographical distribution of most of the species. Two major clades were found in the phylogenetic tree, and an additional lineage, corresponding to P. inflata, was found in the haplotype network obtained for plastid markers. All three lineages are clearly delimited geographically, but, with the exception of P. inflata, the morphological groups were not genetically distinct. Our results suggest that a population expansion after a size reduction resulted in the establishment of two distinct and allopatric groups c. 0.5 Mya, one group occurring in a geologically ancient area, and the other occurring in areas that were under the influence of a series of marine transgressions/regressions at the end of the Pleistocene. These two clades are evolutionarily significant units with significantly different allele frequencies in their nuclear genome and reciprocal monophyly in maternal, uniparentally inherited markers. All our results suggest that the morphology‐based taxonomy in this group does not reflect its evolutionary history, and revision of its species limits should incorporate the distribution of the genetic diversity. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 174 , 199–213.     

Evolutionary drivers of phylogeographical diversity in the highlands of Mexico: a case study of the Crotalus triseriatus species group of montane rattlesnakes

Aim To assess the genealogical relationships of widespread montane rattlesnakes in the Crotalus triseriatus species group and to clarify the role of Late Neogene mountain building and Pleistocene pine–oak forest fragmentation in driving the diversification of Mexican highland taxa. Location Highlands of mainland Mexico and the south‐western United States (Texas, New Mexico, and Arizona). Methods A synthesis of inferences was used to address several associated questions about the biogeography of the Mexican highlands and the evolutionary drivers of phylogeographical diversity in co‐distributed taxa. We combined extensive range‐wide sampling (130 individuals representing five putative species) and mixed‐model phylogenetic analyses of 2408 base pairs of mitochondrial DNA to estimate genealogical relationships and divergence times within the C. triseriatus species group. We then assessed the tempo of diversification using a maximum likelihood framework based on the birth–death process. Estimated times of divergences provided a probabilistic temporal component and questioned whether diversification rates have remained constant or varied over time. Finally, we looked for phylogeographical patterns in other co‐distributed taxa. Results We identified eight major lineages within the C. triseriatus group, and inferred strong correspondence between maternal and geographic history within most lineages. At least one cryptic species was detected. Relationships among lineages were generally congruent with previous molecular studies, with differences largely attributable to our expanded taxonomic and geographic sampling. Estimated divergences between most major lineages occurred in the Late Miocene and Pliocene. Phylogeographical structure within each lineage appeared to have been generated primarily during the Pleistocene. Although the scale of genetic diversity recognized affected estimated rates of diversification, rates appeared to have been constant through time. Main conclusions The biogeographical history of the C. triseriatus group implies a dynamic history for the highlands of Mexico. The Neogene formation of the Transvolcanic Belt appears responsible for structuring geographic diversity among major lineages. Pleistocene glacial–interglacial climatic cycles and resultant expansions and contractions of the Mexican pine–oak forest appear to have driven widespread divergences within lineages. Climatic change, paired with the complex topography of Mexico, probably produced a myriad of species‐specific responses in co‐distributed Mexican highland taxa. The high degree of genetic differentiation recovered in our study and others suggests that the Mexican highlands may contain considerably more diversity than currently recognized.