Post-glacial recolonization of the Great Lakes region by the common gartersnake (Thamnophis sirtalis) inferred from mtDNA sequences

Pleistocene events played an important role in the differentiation of North American vertebrate populations. Michigan, in particular, and the Great Lakes region, in general, were greatly influenced by the last glaciation. While several hypotheses regarding the recolonization of this region have been advanced, none have been strongly supported. We generated 148 complete ND2 mitochondrial DNA (mtDNA) sequences from common gartersnake (Thamnophis sirtalis) populations throughout the Great Lakes region to evaluate phylogeographic patterns and population structure and to determine whether the distribution of haplotypic variants is related to the post-Pleistocene retreat of the Wisconsinan glacier. The common gartersnake was utilized, as it is believed to have been one of the primary vertebrate invaders of the Great Lakes region following the most recent period of glacial retreat and because it has been a model species for a variety of evolutionary, ecological, behavioral, and physiological studies. Several genetically distinct evolutionary lineages were supported by both genealogical and molecular population genetic analyses, although to different degrees. The geographic distribution of the majority of these lineages is interpreted as reflecting post-glacial recolonization dynamics during the late Pleistocene. These findings generally support previous hypotheses of range expansion in this region.     

Discordant temporal and geographic patterns in maternal lineages of eastern North American frogs, Rana catesbeiana (Ranidae) and Pseudacris crucifer (Hylidae)

Whether the Pleistocene has had a disproportionate impact on the recent diversification of temperate species, or played a lesser role in a more protracted process, has been a prominent evolutionary debate for the past decade. We used cytochrome b sequences to reconstruct the evolutionary histories of two widely co-distributed, and ecologically divergent frogs (Rana catesbeiana and Pseudacris crucifer) to examine the role of the Pleistocene in structuring these species. Results for R. catesbeiana reflect a pattern of allopatric fragmentation, likely in Coastal Plain refugia on either side of the Mississippi River dating to the mid to early Pleistocene. In contrast, P. crucifer contains numerous divergent lineages, including one west of the Mississippi River in the Interior Highlands, and in the east, multiple lineages that likely expanded from a number of southern Appalachian refugia with lineage sundering originating in the late Pliocene. Large-scale phylogeographic comparisons between these and other eastern North American species reflect both congruent and independent patterns of diversification, possibly reflecting the relative importance of dispersal ability and habitat associations. Although intra-lineage diversification has been structured by repeated Pleistocene glaciations, lineage sundering likely dates at least to the Pliocene in most (but not all) northern temperate amphibian and reptile species studied to date.     

When recent and evolutionary histories meet: deciphering temporal events from contemporary patterns of mtDNA from fishers (Martes pennanti) in north‐eastern North America

The current spatial distribution of genetic lineages across a region should reflect the complex interplay of both historical and contemporary processes. Postglacial expansion and recolonization in the distant past, in combination with more recent events with anthropogenic effects such as habitat fragmentation and overexploitation, can help shape the pattern of genetic structure observed in contemporary populations. In this study, we characterize the spatial distribution of mtDNA lineages for fisher (Martes pennanti) in north‐eastern North America. The history of fishers in this region is well understood and thus provides an opportunity to interpret patterns of genetic structure in the light of known historical (e.g. recolonization from glacial refugia) and contemporary events (e.g. reintroductions, fragmentation and natural recolonization). Our results indicate that fishers likely recolonized north‐eastern North America from a single Pleistocene refugium. Three genetically distinct remnant populations persisted through the population declines of the 1800s and served as sources for multiple reintroductions and natural recolonizations that have restored the fisher throughout north‐eastern North America. However, the spatial genetic structure of genetic lineages across the region still reflects the three remnant populations.     

Molecular evolution of the trnL(UAA) intron in bryophytes

Structure, variability, and molecular evolution of the trnL(UAA) intron in bryophytes (mosses and liverworts) is analyzed based on more than 1000 sequences representing all classes, including comparisons of lengths and GC-contents, sequence similarities, evolutionary rates and ti/tv ratios of the major lineages and selected genera. Secondary structure analyses of the more variable stem-loop regions facilitated recognition of sequence repeats and minute inversions that often occurred independently in non-related lineages, thus supporting alignment construction and homology assessment. The most length-variable stem-loop region P8 does not share a common evolutionary history across all major bryophyte lineages. Independent nucleotide additions such as internally repeated sequence segments resulted in non-homologous P8 sequences that cannot be folded into a common P8 secondary structure, neither for all bryophytes nor for liverworts or mosses. To address evolutionary patterns, separate analyses of P6/P8 and the remaining intron (core) have to be performed, as overall values of the complete intron are misleading. It is argued that a transition bias observed above the genus level in the core structure is caused by structural constraints, not by its higher GC-content in comparison to the more AT-rich P6 and P8. Compensating base pair changes detected in highly conserved elements are often characteristic of the major bryophyte lineages (classes). Sequence divergence and evolutionary rates are generally higher in liverworts than in mosses, resulting in ambiguous alignments of P6 and P8 even within classes. In mosses, trends towards length reduction of P8 and lower evolutionary rates of the intron are observed. Average intraspecific variation is less than 1%, corresponding to 2-3 mutations in the complete intron.     

Phylogenetic relationships of the North American chorus frogs (Pseudacris: Hylidae)

We examined phylogenetic relationships of the North American chorus frogs (Pseudacris: Hylidae) from 38 populations using 2.4 kb of 12S and 16S mtDNA to elucidate species relationships and examine congruence of previous phylogenetic hypotheses. Parsimony, maximum likelihood, and Bayesian phylogenies are consistent and reveal four strongly supported clades within Pseudacris: (1). A West Coast Clade containing regilla and cadaverina, (2). a Fat Frog Clade including ornata, streckeri, and illinoensis, (3). a Crucifer Clade consisting of crucifer and ocularis, and (4). a Trilling Frog Clade containing all other Pseudacris. Explicit hypothesis testing using parametric bootstrapping indicates that previous phylogenetic hypotheses are rejected by our sequence dataset. Within the Trilling Frog Clade, brimleyi and brachyphona form the sister group to the Nigrita Clade: nigrita, feriarum, triseriata, kalmi, clarkii, and maculata. The Nigrita Clade shows geographic division into three clades: (1). populations of maculata and triseriata west of the Mississippi River and Canadian populations, (2). southeastern US populations of feriarum and nigrita, and (3). northeastern US populations of feriarum, kalmi, and triseriata. We find that subspecific epithets for crucifer (crucifer and bartramiana) and nigrita (nigrita and verrucosa) are uninformative, therefore we discourage recognition of these subspecies. Pseudacris regilla, cadaverina, ocularis, and crucifer are maintained in Pseudacris.     

Comparative phylogeography of two coastal polychaete tubeworms in the Northeast Atlantic supports shared history and vicariant events

The historic processes which have led to the present-day patterns of genetic structure in the marine coastal fauna of the Northeast Atlantic are still poorly understood. While tectonic uplifts and changes in sea level may have caused large-scale vicariance, warmer conditions during glacial maxima may have allowed pockets of diversity to persist to a much wider extent than in the Northwestern Atlantic. The large-scale geographic distribution of deeply divergent lineages of the coastal polychaete tubeworms Pectinaria koreni (two clades) and Owenia fusiformis (three clades) were compared using a fragment of the mitochondrial cytochrome oxidase I gene (mtCOI). All lineages were present along the biogeographic transition zone on the north coast of Brittany (France) and we found evidence pointing towards congruence in the timing of cladogenic events between Pectinaria sp. (P. auricoma/P. belgica and P. koreni) and Owenia sp., suggesting a shared history of vicariant events. More conserved 16SrRNA sequences obtained from four species of Pectinariidae together with mtCOI sequences of P. koreni seem consistent with an initial establishment of pectinariids in the north, and a southward colonization of the Northeast Atlantic. Phylogeographic patterns in O. fusiformis were also consistent with a north/south pattern of lineage splitting and congruent levels of divergence were detected between lineages of both species. We observed signatures of both persistence in small northern glacial refugia, and of northwards range expansion from regions situated closer to the Mediterranean. However, whether the recolonization of the Northeast Atlantic by both species actually reflects separate interglacial periods is unclear with regards to the lack of molecular clock calibration in coastal polychaete species.     

Evolutionary dynamics of the Pgk1 gene in the polyploid genus Kengyilia (Triticeae: Poaceae) and its diploid relatives

The level and pattern of nucleotide variation in duplicate gene provide important information on the evolutionary history of polyploids and divergent process between homoeologous loci within lineages. Kengyilia is a group of allohexaploid species with the StYP genomic constitutions in the wheat tribe. To investigate the evolutionary dynamics of the Pgk1 gene in Kengyilia and its diploid relatives, three copies of Pgk1 homoeologues were isolated from all sampled hexaploid Kengyilia species and analyzed with the Pgk1 sequences from 47 diploid taxa representing 18 basic genomes in Triticeae. Sequence diversity patterns and genealogical analysis suggested that (1) Kengyilia species from the Central Asia and the Qinghai-Tibetan plateau have independent origins with geographically differentiated P genome donors and diverged levels of nucleotide diversity at Pgk1 locus; (2) a relatively long-time sweep event has allowed the Pgk1 gene within Agropyron to adapt to cold climate triggered by the recent uplifts of the Qinghai-Tibetan Plateau; (3) sweep event and population expansion might result in the difference in the d(N)/d(S) value of the Pgk1 gene in allopatric Agropyron populations, and this difference may be genetically transmitted to Kengyilia lineages via independent polyploidization events; (4) an 83 bp MITE element insertion has shaped the Pgk1 loci in the P genome lineage with different geographical regions; (5) the St and P genomes in Kengyilia were donated by Pseudoroegneria and Agropyron, respectively, and the Y genome is closely related to the Xp genome of Peridictyon sanctum. The interplay of evolutionary forces involving diverged natural selection, population expansion, and transposable events in geographically differentiated P genome donors could attribute to geographical differentiation of Kengyilia species via independent origins.     

Fine-scale spatial partitioning of genetic variation and evolutionary contestability in the invasive estuarine mussel Xenostrobus securis

Initial studies of Australian populations of the widely invasive mussel Xenostrobus securis raised the possibility of non-random spatial partitioning within estuaries of mitochondrial DNA (mtDNA) clades that have dispersed intercontinentally and those that have not. A fine-scale phylogeographic investigation was made to examine this possibility using cytochrome c oxidase DNA sequences and data from five microsatellite loci developed here. X. securis in the study region was found to comprise multiple robustly supported mtDNA clades, many of which were genetically very distinct. Frequently, the clades comprised numerous haplotypes that have narrow ranges and may have evolved in situ within drainages or nearby. Microsatellite data reveal significant intra- and inter-drainage differentiation but at lower levels than mtDNA. Variation in mtDNA appears to reflect not only recolonization after infrequent local extinction within drainages but also the low probability of migratory haplotypes successfully establishing during evolutionary contests, with resident haplotypes having selective or demographic advantages derived from local occupation. The patterns of mtDNA haplotype distribution suggest that central New South Wales is the source of the internationally invasive lineages of X. securis.     

Evolutionary history of the scorpionfly Dicerapanorpa magna (Mecoptera,Panorpidae)

Climate changes can have fundamental impacts on the distributional patterns of montane species, and range shifts frequently lead to allopatric divergence followed by the establishment of secondary contact zones. Many European and North American organisms have retreated to southern refugia during glacial periods and colonized northward during postglacial periods, but little is known about the evolutionary response of cold‐adapted insects to Pleistocene climate changes in eastern Asia. The scorpionfly Dicerapanorpa magna (Chou), with cold temperate habitat preference and weak dispersal ability, provides a good model system to explore how climate changes have influenced the distribution and divergence of cold‐adapted insects in eastern Asia. This study reconstructed the demographic dynamics and evolutionary history of D. magna with phylogeographic approaches, and predicted the species’ suitable areas under the Last Glacial Maximum (LGM) and current scenarios with the ecological niche modelling analysis. The mitochondrial cytochrome c oxidase subunit I resolved three phylogenetic lineages in D. magna dating back to Pleistocene, corresponding well with the geographically isolated Qinling, Bashan and Minshan Mountains. The ecological niche modelling recovered the suitable habitats for D. magna were the Qinling and Bashan Mountains under LGM and current conditions. The three lineages of D. magna might be in a process of incipient speciation, and likely derived their current distribution from separate glacial origins, followed by vicariance and divergence.     

Pleistocene evolutionary history of the Clouded Apollo (Parnassius mnemosyne): genetic signatures of climate cycles and a 'time-dependent' mitochondrial substitution rate

Genetic data are currently providing a large amount of new information on past distribution of species and are contributing to a new vision of Pleistocene ice ages. Nonetheless, an increasing number of studies on the 'time dependency' of mutation rates suggest that date assessments for evolutionary events of the Pleistocene might be overestimated. We analysed mitochondrial (mt) DNA (COI) sequence variation in 225 Parnassius mnemosyne individuals sampled across central and eastern Europe in order to assess (i) the existence of genetic signatures of Pleistocene climate shifts; and (ii) the timescale of demographic and evolutionary events. Our analyses reveal a phylogeographical pattern markedly influenced by the Pleistocene/Holocene climate shifts. Eastern Alpine and Balkan populations display comparatively high mtDNA diversity, suggesting multiple glacial refugia. On the other hand, three widely distributed and spatially segregated lineages occupy most of northern and eastern Europe, indicating postglacial recolonization from different refugial areas. We show that a conventional 'phylogenetic' substitution rate cannot account for the present distribution of genetic variation in this species, and we combine phylogeographical pattern and palaeoecological information in order to determine a suitable intraspecific rate through a Bayesian coalescent approach. We argue that our calibrated 'time-dependent' rate (0.096 substitutions/ million years), offers the most convincing time frame for the evolutionary events inferred from sequence data. When scaled by the new rate, estimates of divergence between Balkan and Alpine lineages point to c. 19 000 years before present (last glacial maximum), and parameters of demographic expansion for northern lineages are consistent with postglacial warming (5-11 000 years before present).     

Biogeographic history of an Australian freshwater shrimp, Paratya australiensis (Atyidae): the role life history transition in phylogeographic diversification

The widespread distribution of the freshwater shrimp Paratya australiensis in eastern Australia suggests that populations of this species have been connected in the past. Amphidromy is ancestral in these shrimps, although many extant populations are known to be restricted to freshwater habitats. In this study, we used a fragment of the cytochrome c oxidase I mitochondrial DNA (mtDNA) gene to examine diversity within P. australiensis and to assess the relative importance of amphidromy in its evolutionary history. We hypothesized that if transitions from an amphidromous to a freshwater life history were important, then we would find a number of divergent lineages restricted to single or groups of nearby drainages. Alternatively, if amphidromy was maintained within the species historically, we expected to find lineages distributed over many drainages. We assumed that the only way for divergence to occur within amphidromous lineages was if dispersal was limited to between nearby estuaries, which, during arid periods in the earth's history, became isolated from one another. We found nine highly divergent mtDNA lineages, estimated to have diverged from one another in the late Miocene/early Pliocene, when the climate was more arid than at present. Despite this, the geographic distribution of lineages and haplotypes within lineages did not support the notion of a stepping-stone model of dispersal between estuaries. We conclude that the extensive divergence has most likely arisen through a number of independent amphidromy-freshwater life history transitions, rather than via historical isolation of amphidromy populations. We also found evidence for extensive movement between coastal and inland drainages, supporting the notion that secondary contact between lineages may have occurred as a result of drainage rearrangements. Finally, our data indicate that P. australiensis is likely a complex of cryptic species, some of which are widely distributed, and others geographically restricted.     

Fossil evidence and phylogeography of temperate species: 'glacial refugia' and post-glacial recolonization

We present a short synthesis of the Pleistocene distribution dynamics and phylogeographic recolonization hypotheses for two temperate European mammal species, the red deer ( Cervus elaphus ) and the roe deer ( Capreolus capreolus ), for which high-resolution patterns of fossil evidence and genetic data sets are available. Such data are critical to an understanding of the role of hypothesized glacial refugia. Both species show a similar pattern: a relatively wide distribution in the southern part of Central Europe 60,000–25,000 years ago, and a strong restriction to areas in southern Europe for nearly 10,000 years during the Last Glacial Maximum (LGM) and the early Late Glacial (25,000–14,700 years ago). With the beginning of Greenland Interstadial 1 (Bølling/Allerød warming, c. 14,700–11,600 years ago) a sudden range expansion into Central Europe is visible, but the colonization of most of Central Europe, including the northern European Lowlands, only began in the early Holocene. In a European context, regions where the species were distributed during the LGM and early Late Glacial are most relevant as potential origins of recolonization processes, because during these c. 10,000 years distribution ranges were smaller than at any other time in the Late Quaternary. As far as the present distribution of temperate species and their genetic lineages is concerned, so-called 'cryptic refugia' are important only if the species are actually confirmed there during the LGM, as otherwise they could not possibly have contributed to the recolonization that eventually resulted in the present distribution ranges.     

Mitochondrial gene diversity in the common vole Microtus arvalis shaped by historical divergence and local adaptations

The phylogeography of the common vole (Microtus arvalis) was examined by analysing mitochondrial DNA (mtDNA) sequence variation in 1044 base pairs (bp) of the cytochrome b (cytb) gene and in 322 bp of the control region (ctr) among 106 individuals from 58 locations. The geographical distribution of four previously recognized cytb evolutionary lineages in Europe was refined and a new lineage was found in southern Germany. All lineages were distributed allopatrically, except in one sample that was probably located in a contact zone. The occurrence of several lineages in the Alps is in keeping with their recent recolonization from distinct sources. The translation of 84 cytb DNA sequences produced 33 distinct proteins with relationships that differed from those of the DNA haplotypes, suggesting that the mtDNA lineages did not diverge in response to selection. In comparison with M. agrestis, a neutrality test detected no overall evidence for selection in the cytb gene, but a closer examination of a structural model showed that evolutionarily conserved and functionally important positions were often affected. A new phylogeographical test of random accumulation of nonsynonymous mutations generated significant results in three lineages. We therefore conclude that the molecular diversity of cytb in M. arvalis is overall the result of the demographic history of the populations, but that there have been several episodes of local adaptation to peculiar environments.     

WHEN COLD IS BETTER: CLIMATE-DRIVEN ELEVATION SHIFTS YIELD COMPLEX PATTERNS OF DIVERSIFICATION AND DEMOGRAPHY IN AN ALPINE SPECIALIST (AMERICAN PIKA, OCHOTONA PRINCEPS)

The genetic consequences of climate-driven range fluctuation during the Pleistocene have been well studied for temperate species, but cold-adapted (e.g., alpine, arctic) species that may have responded uniquely to past climatic events have received less attention. In particular, we have no a priori expectation for long-term evolutionary consequences of elevation shifts into and out of sky islands by species adapted to alpine habitats. Here, we examined the influence of elevation shifts on genetic differentiation and historical demography in an alpine specialist, the American pika ( Ochotona princeps ). Pika populations are divided into five genetic lineages that evolved in association with separate mountain systems, rather than lineages that reflect individual sky islands. This suggests a role for glacial-period elevation shifts in promoting gene flow among high-elevation populations and maintaining regional cohesion of genetic lineages. We detected a signature of recent demographic decline in all lineages, consistent with the expectation that Holocene climate warming has driven range retraction in southern lineages, but unexpected for northern populations that presumably represent postglacial expansion. An ecological niche model of past and future pika distributions highlights the influence of climate on species range and indicates that the distribution of genetic diversity may change dramatically with continued climate warming.     

Postglacial expansion and not human influence best explains the population structure in the endangered kea (Nestor notabilis)

Inferring past demography is a central question in evolutionary and conservation biology. It is, however, sometimes challenging to infer the processes that shaped the current patterns of genetic variation in endangered species. Population substructuring can occur as a result of survival in several isolated refugia and subsequent recolonization processes or via genetic drift following a population decline. The kea (Nestor notabilis) is an endemic parrot widely distributed in the mountains of the South Island of New Zealand that has gone through a major human‐induced population decline during the 1860s–1970s. The aims of this study were to understand the glacial and postglacial history of kea and to determine whether the recent population decline played a role in the shaping of the current genetic variation. We examined the distribution of genetic variation, differentiation and admixture in kea using 17 microsatellites and the mitochondrial control region. Mitochondrial data showed a shallow phylogeny and a genetic distinction between the North and South of the range consistent with the three genetic clusters identified with microsatellite data. Both marker types indicated an increase in genetic isolation by geographic distance. Approximate Bayesian Computation supported a scenario of postglacial divergence from a single ancestral glacial refugium, suggesting that the contemporary genetic structure has resulted from recolonization processes rather than from a recent population decline. The recent evolutionary origin of this genetic structure suggests that each genetic cluster does not need to be considered as independent conservation units.     

The role of topography in structuring the demographic history of the pine processionary moth,Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

Aim We investigated the Quaternary history of the pine processionary moth, Thaumetopoea pityocampa, an oligophagous insect currently expanding its range. We tested the potential role played by mountain ranges during the post‐glacial recolonization of western Europe. Location Western Europe, with a focus on the Pyrenees, Massif Central and western Alps. Methods Maternal genetic structure was investigated using a fragment of the mitochondrial cytochrome c oxidase subunit I (COI) gene. We analysed 412 individuals from 61 locations and performed maximum likelihood and maximum parsimony phylogenetic analyses and hierarchical analysis of molecular variance, and we investigated signs of past expansion. Results A strong phylogeographic pattern was found, with two deeply divergent clades. Surprisingly, these clades were not separated by the Pyrenees but rather were distributed from western to central Iberia and from eastern Iberia to the Italian Peninsula, respectively. This latter group consisted of three shallowly divergent lineages that exhibited strong geographic structure and independent population expansions. The three identified lineages occurred: (1) on both sides of the Pyrenean range, with more genetically diverse populations in the east, (2) from eastern Iberia to western France, with a higher genetic diversity in the south, and (3) from the western Massif Central to Italy. Admixture areas were found at the foot of the Pyrenees and Massif Central. Main conclusions The identified genetic lineages were geographically structured, but surprisingly the unsuitable high‐elevation areas of the main mountainous ranges were not responsible for the spatial separation of genetic groups. Rather than acting as barriers to dispersal, mountains appear to have served as refugia during the Pleistocene glaciations, and current distributions largely reflect expansion from these bottlenecked refugial populations. The western and central Iberian clade did not contribute to the northward post‐glacial recolonization of Europe, yet its northern limit does not correspond to the Pyrenees. The different contributions of the identified refugia to post‐glacial expansion might be explained by differences in host plant species richness. For example, the Pyrenean lineage could have been trapped elevationally by tracking montane pines, while the eastern Iberian lineage could have expanded latitudinally by tracking thermophilic lowland pine species.     

Establishing a perimeter position: speciation around the Indian Ocean Basin

Historical biological interactions among peripheral landmasses on the periphery of the Indian Ocean Basin (IOB) are generally poorly understood. While interactions based on early Gondwanan vicariance have been used to explain present day lineage distributions, several recent studies have instead inferred dispersal across the IOB. This inference is often advanced because lineages under study have species inhabiting IOB islands. Here we examine the roles of continental vicariance vs. trans-IOB dispersal in the distribution of an avian genus found around the perimeter of the IOB. A molecular phylogeny does reveal evidence of a relationship that would require the inference of trans-IOB dispersal between eastern Africa and Sri Lanka. However, molecular clock data, ancestral area analyses and paleoclimatic reconstructions suggest that vicariance related to tropical forest expansion and retraction is more likely to have facilitated African-Asian interchange, with an initial colonization of Africa from Asia quickly followed by a recolonization of Asia. Subsequent dispersal from Asia to Sri Lanka and islands east of the Sunda Shelf are inferred; these latter islands were colonized in a stepping-stone fashion that culminated in colonization of the Sunda Shelf, and a recolonization of mainland Asia. We propose that circum-IOB distributions, which post-date early Gondwanan breakup, are most likely the result of continent-based vicariant events, particularly those events related to large-scale habitat alterations, and not trans-IOB dispersals.     

Rapid postglacial diversification and long‐term stasis within the songbird genus Junco: phylogeographic and phylogenomic evidence

Natural systems composed of closely related taxa that vary in the degree of phenotypic divergence and geographic isolation provide an opportunity to investigate the rate of phenotypic diversification and the relative roles of selection and drift in driving lineage formation. The genus Junco (Aves: Emberizidae) of North America includes parapatric northern forms that are markedly divergent in plumage pattern and colour, in contrast to geographically isolated southern populations in remote areas that show moderate phenotypic divergence. Here, we quantify patterns of phenotypic divergence in morphology and plumage colour and use mitochondrial DNA genes, a nuclear intron, and genomewide SNPs to reconstruct the demographic and evolutionary history of the genus to infer relative rates of evolutionary divergence among lineages. We found that geographically isolated populations have evolved independently for hundreds of thousands of years despite little differentiation in phenotype, in sharp contrast to phenotypically diverse northern forms, which have diversified within the last few thousand years as a result of the rapid postglacial recolonization of North America. SNP data resolved young northern lineages into reciprocally monophyletic lineages, indicating low rates of gene flow even among closely related parapatric forms, and suggesting a role for strong genetic drift or multifarious selection acting on multiple loci in driving lineage divergence. Juncos represent a compelling example of speciation in action, where the combined effects of historical and selective factors have produced one of the fastest cases of speciation known in vertebrates.     

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.     

Molecular phylogeography of Jerdon’s pitviper (Protobothrops jerdonii): importance of the uplift of the Tibetan plateau

Aim To investigate the population history and demographics of Jerdon’s pitviper, Protobothrops jerdonii, and elucidate how the unique physical conditions and heterogeneous mountain environments resulting from the uplift of the Tibetan Plateau shaped the genetic diversity and evolutionary history of the species. Location China and Vietnam. Methods We sequenced and analysed a total of 1752 base pairs from two mitochondrial genes, cytochrome b (cyt b) and NADH dehydrogenase subunit (ND4), for 81 specimens sampled from 27 localities across the species’ range, and a total of 464 base pairs from two nuclear genes for 28 representative samples from all mitochondrial DNA lineages. Based on these data, we constructed the genealogical relationships and estimated the divergence times of the mitochondrial DNA clades. Results The mitochondrial DNA results revealed the existence of five distinct, strongly supported and geographically structured DNA lineages within populations of P. jerdonii that are paraphyletic with respect to Protobothrops xiangchengensis. Estimation of divergence dates suggested that P. jerdonii possibly evolved in the western Hengduan Mountains region c. 6.6 Ma in the late Miocene. Nuclear DNA data did not provide sufficient resolution to distinguish the mitochondrial DNA lineages. Main conclusions Based on the present‐day distribution and intraspecific genealogy, the evolutionary history of P. jerdonii can be explained by a pattern of dispersal followed by vicariance. All lines of evidence suggest that historical biogeographical factors, particularly the north–south orientation of the higher mountains, as well as low‐elevation areas in western China, had the greatest influence on the population structure, lineage formation and species distribution of this snake. However, highly heterogeneous habitats and glacial cycles appear to have affected patterns of intraspecific differentiation. While our mitochondrial data provide evidence for clear phylogeographical structure, our small sampling of nuclear genes does not, suggesting that nuclear markers may not have had sufficient time to coalesce to match patterns observed in the mitochondrial data.