Of glaciers and refugia: a decade of study sheds new light on the phylogeography of northwestern North America

DNA fingerprinting in botany: past, present, future

Almost three decades ago Alec Jeffreys published his seminal Nature papers on the use of minisatellite probes for DNA fingerprinting of humans (Jeffreys and colleagues Nature 1985, 314:67-73 and Nature 1985, 316:76-79). The new technology was soon adopted for many other organisms including plants, and when Hilde Nybom, Kurt Weising and Alec Jeffreys first met at the very First International Conference on DNA Fingerprinting in Berne, Switzerland, in 1990, everybody was enthusiastic about the novel method that allowed us for the first time to discriminate between humans, animals, plants and fungi on the individual level using DNA markers. A newsletter coined "Fingerprint News" was launched, T-shirts were sold, and the proceedings of the Berne conference filled a first book on "DNA fingerprinting: approaches and applications". Four more conferences were about to follow, one on each continent, and Alec Jeffreys of course was invited to all of them. Since these early days, methodologies have undergone a rapid evolution and diversification. A multitude of techniques have been developed, optimized, and eventually abandoned when novel and more efficient and/or more reliable methods appeared. Despite some overlap between the lifetimes of the different technologies, three phases can be defined that coincide with major technological advances. Whereas the first phase of DNA fingerprinting ("the past") was dominated by restriction fragment analysis in conjunction with Southern blot hybridization, the advent of the PCR in the late 1980s gave way to the development of PCR-based single- or multi-locus profiling techniques in the second phase. Given that many routine applications of plant DNA fingerprinting still rely on PCR-based markers, we here refer to these methods as "DNA fingerprinting in the present", and include numerous examples in the present review. The beginning of the third phase actually dates back to 2005, when several novel, highly parallel DNA sequencing strategies were developed that increased the throughput over current Sanger sequencing technology 1000-fold and more. High-speed DNA sequencing was soon also exploited for DNA fingerprinting in plants, either in terms of facilitated marker development, or directly in the sense of "genotyping-by-sequencing". Whereas these novel approaches are applied at an ever increasing rate also in non-model species, they are still far from routine, and we therefore treat them here as "DNA fingerprinting in the future".

First report of Taenia arctos (Cestoda: Taeniidae) from grizzly (Ursus arctos horribilis) and black bears (Ursus americanus) in North America

The cestode Taenia arctos was found at necropsy in the small intestine of a grizzly (Ursus arctos horribilis) and a black bear (Ursus americanus) from Kananaskis Country in southwestern Alberta, Canada. The autolysis of the tapeworm specimens precluded detailed morphological characterization of the parasites but molecular analysis based on mitochondrial DNA cytochrome c oxidase subunit 1 gene confirmed their identity as T. arctos. This is the first report of T. arctos from definitive hosts in North America. Its detection in Canadian grizzly and black bears further supports the Holarctic distribution of this tapeworm species and its specificity for ursids as final hosts. Previously, T. arctos was unambiguously described at its adult stage in brown bears (Ursus arctos arctos) from Finland, and as larval stages in Eurasian elk (Alces alces) from Finland and moose (Alces americanus) from Alaska, USA. Given the morphological similarity between T. arctos and other Taenia species, the present study underlines the potential for misidentification of tapeworm taxa in previous parasitological reports from bears and moose across North America. The biogeographical history of both definitive and intermediate hosts in the Holarctic suggests an ancient interaction between U. arctos, Alces spp., and T. arctos, and a relatively recent host-switching event in U. americanus.

Phylogeography and post-glacial recolonization in wolverines (Gulo gulo) from across their circumpolar distribution

Interglacial-glacial cycles of the Quaternary are widely recognized in shaping phylogeographic structure. Patterns from cold adapted species can be especially informative - in particular, uncovering additional glacial refugia, identifying likely recolonization patterns, and increasing our understanding of species' responses to climate change. We investigated phylogenetic structure of the wolverine, a wide-ranging cold adapted carnivore, using a 318 bp of the mitochondrial DNA control region for 983 wolverines (n=209 this study, n=774 from GenBank) from across their full Holarctic distribution. Bayesian phylogenetic tree reconstruction and the distribution of observed pairwise haplotype differences (mismatch distribution) provided evidence of a single rapid population expansion across the wolverine's Holarctic range. Even though molecular evidence corroborated a single refugium, significant subdivisions of population genetic structure (0.01< ΦST <0.99, P<0.05) were detected. Pairwise ΦST estimates separated Scandinavia from Russia and Mongolia, and identified five main divisions within North America - the Central Arctic, a western region, an eastern region consisting of Ontario and Quebec/Labrador, Manitoba, and California. These data are in contrast to the nearly panmictic structure observed in northwestern North America using nuclear microsatellites, but largely support the nuclear DNA separation of contemporary Manitoba and Ontario wolverines from northern populations. Historic samples (c. 1900) from the functionally extirpated eastern population of Quebec/Labrador displayed genetic similarities to contemporary Ontario wolverines. To understand these divergence patterns, four hypotheses were tested using Approximate Bayesian Computation (ABC). The most supported hypothesis was a single Beringia incursion during the last glacial maximum that established the northwestern population, followed by a west-to-east colonization during the Holocene. This pattern is suggestive of colonization occurring in accordance with glacial retreat, and supports expansion from a single refugium. These data are significant relative to current discussions on the conservation status of this species across its range.

Deep genetic divergence between disjunct Refugia in the Arctic-Alpine King's Crown, Rhodiola integrifolia (Crassulaceae)

Despite the strength of climatic variability at high latitudes and upper elevations, we still do not fully understand how plants in North America that are distributed between Arctic and alpine areas responded to the environmental changes of the Quaternary. To address this question, we set out to resolve the evolutionary history of the King's Crown, Rhodiola integrifolia using multi-locus population genetic and phylogenetic analyses in combination with ecological niche modeling. Our population genetic analyses of multiple anonymous nuclear loci revealed two major clades within R. integrifolia that diverged from each other ~ 700 kya: one occurring in Beringia to the north (including members of subspecies leedyi and part of subspecies integrifolia), and the other restricted to the Southern Rocky Mountain refugium in the south (including individuals of subspecies neomexicana and part of subspecies integrifolia). Ecological niche models corroborate our hypothesized locations of refugial areas inferred from our phylogeographic analyses and revealed some environmental differences between the regions inhabited by its two subclades. Our study underscores the role of geographic isolation in promoting genetic divergence and the evolution of endemic subspecies in R. integrifolia. Furthermore, our phylogenetic analyses of the plastid spacer region trnL-F demonstrate that among the native North American species, R. integrifolia and R. rhodantha are more closely related to one another than either is to R. rosea. An understanding of these historic processes lies at the heart of making informed management decisions regarding this and other Arctic-alpine species of concern in this increasingly threatened biome.

Genetic roadmap of the Arctic: plant dispersal highways, traffic barriers and capitals of diversity

We provide the first comparative multispecies analysis of spatial genetic structure and diversity in the circumpolar Arctic using a common strategy for sampling and genetic analyses. We aimed to identify and explain potential general patterns of genetic discontinuity/connectivity and diversity, and to compare our findings with previously published hypotheses. We collected and analyzed 7707 samples of 17 widespread arctic-alpine plant species for amplified fragment length polymorphisms (AFLPs). Genetic structure, diversity and distinctiveness were analyzed for each species, and extrapolated to cover the geographic range of each species. The resulting maps were overlaid to produce metamaps. The Arctic and Atlantic Oceans, the Greenlandic ice cap, the Urals, and lowland areas between southern mountain ranges and the Arctic were the strongest barriers against gene flow. Diversity was highest in Beringia and gradually decreased into formerly glaciated areas. The highest degrees of distinctiveness were observed in Siberia. We conclude that large-scale general patterns exist in the Arctic, shaped by the Pleistocene glaciations combined with long-standing physical barriers against gene flow. Beringia served as both refugium and source for interglacial (re)colonization, whereas areas further west in Siberia served as refugia, but less as sources for (re)colonization.

The diversity and biogeography of the Coleoptera of Churchill: insights from DNA barcoding

BACKGROUND

Coleoptera is the most diverse order of insects (>300,000 described species), but its richness diminishes at increasing latitudes (e.g., ca. 7400 species recorded in Canada), particularly of phytophagous and detritivorous species. However, incomplete sampling of northern habitats and a lack of taxonomic study of some families limits our understanding of biodiversity patterns in the Coleoptera. We conducted an intensive biodiversity survey from 2006-2010 at Churchill, Manitoba, Canada in order to quantify beetle species diversity in this model region, and to prepare a barcode library of beetles for sub-arctic biodiversity and ecological research. We employed DNA barcoding to provide estimates of provisional species diversity, including for families currently lacking taxonomic expertise, and to examine the guild structure, habitat distribution, and biogeography of beetles in the Churchill region.

RESULTS

We obtained DNA barcodes from 3203 specimens representing 302 species or provisional species (the latter quantitatively defined on the basis of Molecular Operational Taxonomic Units, MOTUs) in 31 families of Coleoptera. Of the 184 taxa identified to the level of a Linnaean species name, 170 (92.4%) corresponded to a single MOTU, four (2.2%) represented closely related sibling species pairs within a single MOTU, and ten (5.4%) were divided into two or more MOTUs suggestive of cryptic species. The most diverse families were the Dytiscidae (63 spp.), Staphylinidae (54 spp.), and Carabidae (52 spp.), although the accumulation curve for Staphylinidae suggests that considerable additional diversity remains to be sampled in this family. Most of the species present are predatory, with phytophagous, mycophagous, and saprophagous guilds being represented by fewer species. Most named species of Carabidae and Dytiscidae showed a significant bias toward open habitats (wet or dry). Forest habitats, particularly dry boreal forest, although limited in extent in the region, were undersampled.

CONCLUSIONS

We present an updated species list for this region as well as a species-level DNA barcode reference library. This resource will facilitate future work, such as biomonitoring and the study of the ecology and distribution of larvae.

Phylogeographical patterns of an alpine plant, Rhodiola dumulosa (Crassulaceae), inferred from chloroplast DNA sequences

The phylogeographical patterns of Rhodiola dumulosa, an alpine plant species restrictedly growing in the crevices of rock piles, were investigated based on 4 fragments of the chloroplast genome. To cover the full distribution of R. dumulosa in China, 19 populations from 3 major disjunct distribution areas (northern, central, and northwestern China) were sampled. A total of 5881bp (after alignment) of chloroplast DNA (cpDNA) from 100 individuals were sequenced. The combined cpDNA data set yielded 36 haplotypes. The total genetic diversity of R. dumulosa was remarkably high (H(T) = 0.981). The interpopulation genetic differentiation was significantly large (F(ST) = 0.8537, P < 0.001), possibly due to the long-term isolation of the natural populations. N(ST) was significantly larger than G(ST) (P < 0.001), indicating the presence of phylogeographical structure among the R. dumulosa populations. We propose 2 migration steps to explain the current distribution of R. dumulosa in China. First, this species migrated from refugia in the Qinghai-Tibetan Plateau to northern areas via the intervening highlands when temperatures increased; second, the highland populations migrated toward the mountaintops when temperatures increased further because R. dumulosa is adapted to cold environments. During the second migration step, the common ancestral haplotypes may have been gradually lost.

Evaluating the significance of paleophylogeographic species distribution models in reconstructing quaternary range-shifts of nearctic chelonians

The climatic cycles of the Quaternary, during which global mean annual temperatures have regularly changed by 5-10°C, provide a special opportunity for studying the rate, magnitude, and effects of geographic responses to changing climates. During the Quaternary, high- and mid-latitude species were extirpated from regions that were covered by ice or otherwise became unsuitable, persisting in refugial retreats where the environment was compatible with their tolerances. In this study we combine modern geographic range data, phylogeny, Pleistocene paleoclimatic models, and isotopic records of changes in global mean annual temperature, to produce a temporally continuous model of geographic changes in potential habitat for 59 species of North American turtles over the past 320 Ka (three full glacial-interglacial cycles). These paleophylogeographic models indicate the areas where past climates were compatible with the modern ranges of the species and serve as hypotheses for how their geographic ranges would have changed in response to Quaternary climate cycles. We test these hypotheses against physiological, genetic, taxonomic and fossil evidence, and we then use them to measure the effects of Quaternary climate cycles on species distributions. Patterns of range expansion, contraction, and fragmentation in the models are strongly congruent with (i) phylogeographic differentiation; (ii) morphological variation; (iii) physiological tolerances; and (iv) intraspecific genetic variability. Modern species with significant interspecific differentiation have geographic ranges that strongly fluctuated and repeatedly fragmented throughout the Quaternary. Modern species with low genetic diversity have geographic distributions that were highly variable and at times exceedingly small in the past. Our results reveal the potential for paleophylogeographic models to (i) reconstruct past geographic range modifications, (ii) identify geographic processes that result in genetic bottlenecks; and (iii) predict threats due to anthropogenic climate change in the future.

Evolutionary consequences of putative intra-and interspecific hybridization in agaric fungi

Agaric fungi of the southern Appalachian Mountains including Great Smoky Mountains National Park are often heterozygous for the rDNA internal transcribed spacer region (ITS) with >42% of collections showing some heterozygosity for indels and/or base-pair substitutions. For these collections, intra-individual haplotype divergence is typically less than 2%, but for 3% of these collections intra-individual haplotype divergence exceeds that figure. We hypothesize that high intra-individual haplotype divergence is due to hybridization between agaric fungi with divergent haplotypes, possibly migrants from geographically isolated glacial refugia. Four species with relatively high haplotype divergence were examined: Armillaria mellea, Amanita citrina f. lavendula, Gymnopus dichrous and the Hygrocybe flavescens/chlorophana complex. The ITS region was sequenced, haplotypes of heterozygotes were resolved through cloning, and phylogenetic analyses were used to determine the outcome of hybridization events. Within Armillaria mellea and Amanita citrina f. lavendula, we found evidence of interbreeding and recombination. Within G. dichrous and H. flavescens/chlorophana, hybrids were identified but there was no evidence for F2 or higher progeny in natural populations suggesting that the hybrid fruitbodies might be an evolutionary dead end and that the genetically divergent Mendelian populations from which they were derived are, in fact, different species. The association between ITS haplotype divergence of less than 5% (Armillaria mellea = 2.6% excluding gaps; Amanita citrina f. lavendula = 3.3%) with the presence of putative recombinants and greater than 5% (Gymnopus dichrous = 5.7%; Hygrocybe flavescens/chlorophana = 14.1%) with apparent failure of F1 hybrids to produce F2 or higher progeny in populations may suggest a correlation between genetic distance and reproductive isolation.

Genetic and morphological structure of a spruce hybrid (Picea sitchensis x P. glauca) zone along a climatic gradient

PREMISE OF THE STUDY

Historic colonization and contemporary evolutionary processes contribute to patterns of genetic variation and differentiation among populations. However, separating the respective influences of these processes remains a challenge, particularly for natural hybrid zones, where standing genetic variation may result from evolutionary processes both preceding and following contact, influencing the evolutionary trajectory of hybrid populations. Where adaptation to novel environments may be facilitated by interspecific hybridization, teasing apart these processes will have practical implications for forest management in changing environments.

METHODS

We evaluated the neutral genetic architecture of the Picea sitchensis (Sitka spruce) × P. glauca (white spruce) hybrid zone along the Nass and Skeena river valleys in northwestern British Columbia using chloroplast, mitochondrial, and nuclear microsatellite markers, in combination with cone morphological traits.

KEY RESULTS

Sitka spruce mitotype "capture", evidenced by this species dominating the maternal lineage, is consistent with earlier colonization of the region by Sitka spruce. This "capture" differs from the spatial distribution of chloroplast haplotypes, indicating pollen dispersal and its contribution to geographic structure. Genetic ancestry, based on nuclear markers, was strongly influenced by climate and geography. Highly parallel results for replicate transects along environmental gradients provide support for the bounded hybrid superiority model of hybrid zone maintenance. •

CONCLUSIONS

This broad-scale analysis of neutral genetic structure indicates the importance of historic and contemporary gene flow, environmental selection, and their interaction in shaping neutral genetic variation within this hybrid zone, informative to seed transfer development and reforestation for future climates.

Mitochondrial DNA haplotype distribution patterns in Pinus ponderosa (Pinaceae): range-wide evolutionary history and implications for conservation

PREMISE OF THE STUDY

Ponderosa pine (Pinus ponderosa Douglas ex P. Lawson & C. Lawson) exhibits complicated patterns of morphological and genetic variation across its range in western North America. This study aims to clarify P. ponderosa evolutionary history and phylogeography using a highly polymorphic mitochondrial DNA marker, with results offering insights into how geographical and climatological processes drove the modern evolutionary structure of tree species in the region.

METHODS

We amplified the mtDNA nad1 second intron minisatellite region for 3,100 trees representing 104 populations, and sequenced all length variants. We estimated population-level haplotypic diversity and determined diversity partitioning among varieties, races and populations. After aligning sequences of minisatellite repeat motifs, we evaluated evolutionary relationships among haplotypes.

KEY RESULTS

The geographical structuring of the 10 haplotypes corresponded with division between Pacific and Rocky Mountain varieties. Pacific haplotypes clustered with high bootstrap support, and appear to have descended from Rocky Mountain haplotypes. A greater proportion of diversity was partitioned between Rocky Mountain races than between Pacific races. Areas of highest haplotypic diversity were the southern Sierra Nevada mountain range in California, northwestern California, and southern Nevada.

CONCLUSIONS

Pinus ponderosa haplotype distribution patterns suggest a complex phylogeographic history not revealed by other genetic and morphological data, or by the sparse paleoecological record. The results appear consistent with long-term divergence between the Pacific and Rocky Mountain varieties, along with more recent divergences not well-associated with race. Pleistocene refugia may have existed in areas of high haplotypic diversity, as well as the Great Basin, Southwestern United States/northern Mexico, and the High Plains.

On the evolutionary and biogeographic history of Saxifraga sect. Trachyphyllum (Gaud.) Koch (Saxifragaceae Juss.)

Arctic-alpine plants in the genus Saxifraga L. (Saxifragaceae Juss.) provide an excellent system for investigating the process of diversification in northern regions. Yet, sect. Trachyphyllum (Gaud.) Koch, which is comprised of about 8 to 26 species, has still not been explored by molecular systematists even though taxonomists concur that the section needs to be thoroughly re-examined. Our goals were to use chloroplast trnL-F and nuclear ITS DNA sequence data to circumscribe the section phylogenetically, test models of geographically-based population divergence, and assess the utility of morphological characters in estimating evolutionary relationships. To do so, we sequenced both genetic markers for 19 taxa within the section. The phylogenetic inferences of sect. Trachyphyllum using maximum likelihood and Bayesian analyses showed that the section is polyphyletic, with S. aspera L. and S bryoides L. falling outside the main clade. In addition, the analyses supported several taxonomic re-classifications to prior names. We used two approaches to test biogeographic hypotheses: i) a coalescent approach in Mesquite to test the fit of our reconstructed gene trees to geographically-based models of population divergence and ii) a maximum likelihood inference in Lagrange. These tests uncovered strong support for an origin of the clade in the Southern Rocky Mountains of North America followed by dispersal and divergence episodes across refugia. Finally we adopted a stochastic character mapping approach in SIMMAP to investigate the utility of morphological characters in estimating evolutionary relationships among taxa. We found that few morphological characters were phylogenetically informative and many were misleading. Our molecular analyses provide a foundation for the diversity and evolutionary relationships within sect. Trachyphyllum and hypotheses for better understanding the patterns and processes of divergence in this section, other saxifrages, and plants inhabiting the North Pacific Rim.

Model selection as a tool for phylogeographic inference: an example from the willow Salix melanopsis

Phylogeographic inference has typically relied on analyses of data from one or a few genes to provide estimates of demography and population histories. While much has been learned from these studies, all phylogeographic analysis is conditioned on the data, and thus, inferences derived from data that represent a small sample of the genome are unavoidably tenuous. Here, we demonstrate one approach for moving beyond classic phylogeographic research. We use sequence capture probes and Illumina sequencing to generate data from >400 loci in order to infer the phylogeographic history of Salix melanopsis, a riparian willow with a disjunct distribution in coastal and the inland Pacific Northwest. We evaluate a priori phylogeographic hypotheses using coalescent models for parameter estimation, and the results support earlier findings that identified post-Pleistocene dispersal as the cause of the disjunction in S. melanopsis. We also conduct a series of model selection exercises using IMa2, Migrate-n and ∂a∂i. The resulting ranking of models indicates that refugial dynamics were complex, with multiple regions in the inland regions serving as the source for postglacial colonization. Our results demonstrate that new sources of data and new approaches to data analysis can rejuvenate phylogeographic research by allowing for the identification of complex models that enable researchers to both identify and estimate the most relevant parameters for a given system.

When east meets west: population structure of a high-latitude resident species, the boreal chickadee (Poecile hudsonicus)

The population genetic structure of northern boreal species has been strongly influenced both by the Quaternary glaciations and the presence of contemporary barriers, such as mountain ranges and rivers. We used a combination of mitochondrial DNA (mtDNA), nuclear microsatellites and spatial distribution modelling to study the population genetic structure of the boreal chickadee (Poecile hudsonicus), a resident passerine, and to investigate whether historical or contemporary barriers have influenced this northern species. MtDNA data showed evidence of eastern and western groups, with secondary admixture occurring in central Canada. This suggests that the boreal chickadee probably persisted in multiple glacial refugia, one in Beringia and at least one in the east. Palaeo-distribution modelling identified suitable habitat in Beringia (Alaska), Atlantic Canada and the southern United States, and correspond to divergence dates of 60-96 kya. Pairwise FST values for both mtDNA and microsatellites were significant for all comparisons involving Newfoundland, though mtDNA data suggest a more recent separation. Furthermore, unlike mtDNA data, nuclear data support population connectivity among the continental populations, possibly due to male-biased dispersal. Although both are significant, the isolation-by-distance signal is much stronger for mtDNA (r(2)=0.51) than for microsatellites (r(2)=0.05), supporting the hypothesis of male-biased dispersal. The population structure of the boreal chickadee was influenced by isolation in multiple refugia and contemporary barriers. In addition to geographical distance, physical barriers such as the Strait of Belle Isle and northern mountains in Alaska are restricting gene flow, whereas the Rocky Mountains in the west are a porous barrier.

Phylogeography and postglacial expansion of the endangered semi-aquatic mammal Galemys pyrenaicus

Background

Species with strict ecological requirements may provide new insights into the forces that shaped the geographic variation of genetic diversity. The Pyrenean desman, Galemys pyrenaicus, is a small semi-aquatic mammal that inhabits clean streams of the northern half of the Iberian Peninsula and is endangered in most of its geographic range, but its genetic structure is currently unknown. While the stringent ecological demands derived from its aquatic habitat might have caused a partition of the genetic diversity among river basins, Pleistocene glaciations would have generated a genetic pattern related to glacial refugia.

Results

To study the relative importance of historical and ecological factors in the genetic structure of G. pyrenaicus, we used mitochondrial and intronic sequences of specimens covering most of the species range. We show, first, that the Pyrenean desman has very low levels of genetic diversity compared to other mammals. In addition, phylogenetic and dating analyses of the mitochondrial sequences reveal a strong phylogeographic structure of a Middle Pleistocene origin, suggesting that the main lineages arose during periods of glacial isolation. Furthermore, both the spatial distribution of nuclear and mitochondrial diversity and the results of species distribution modeling suggest the existence of a major glacial refugium in the northwestern part of the Iberian Peninsula. Finally, the main mitochondrial lineages show a striking parapatric distribution without any apparent exchange of mitochondrial haplotypes between the lineages that came into secondary contact (although with certain permeability to nuclear genes), indicating incomplete mixing after the post-glacial recolonization. On the other hand, when we analyzed the partition of the genetic diversity among river basins, the Pyrenean desman showed a lower than expected genetic differentiation among main rivers.

Conclusions

The analysis of mitochondrial and intronic markers in G. pyrenaicus showed the predominant effects of Pleistocene glaciations on the genetic structure of this species, while the distribution of the genetic diversity was not greatly influenced by the main river systems. These results and, particularly, the discovery of a marked phylogeographic structure, may have important implications for the conservation of the Pyrenean desman.

Phylogeographic evidence for a link of species divergence of Ephedra in the Qinghai-Tibetan Plateau and adjacent regions to the Miocene Asian aridification

The Qinghai-Tibetan Plateau (QTP) has become one of the hotspots for phylogeographical studies due to its high species diversity. However, most previous studies have focused on the effects of the Quaternary glaciations on phylogeographical structures and the locations of glacial refugia, and little is known about the effects of the aridization of interior Asia on plant population structure and speciation. Here the chloroplast DNA (cpDNA) trnT-trnF and trnS-trnfM sequences were used to investigate the differentiation and phylogeographical history of 14 Ephedra species from the QTP and northern China, based on a sampling of 107 populations. The phylogeographical analysis, together with phylogenetic reconstruction based on combined four cpDNA fragments (rbcL, rpl16, rps4, and trnS-trnfM), supports three main lineages (eastern QTP, southern QTP, and northern China) of these Ephedra species. Divergence of each lineage could be dated to the Middle or Late Miocene, and was very likely linked to the uplift of the QTP and the Asian aridification, given the high drought and/or cold tolerance of Ephedra. Most of the Ephedra species had low intraspecific variation and lacked a strong phylogeographical structure, which could be partially attributed to clonal reproduction and a relatively recent origin. In addition, ten of the detected 25 cpDNA haplotypes are shared among species, suggesting that a wide sampling of species is helpful to investigate the origin of observed haplotypes and make reliable phylogeographical inference. Moreover, the systematic positions of some Ephedra species are discussed.

Impact of climate changes from Middle Miocene onwards on evolutionary diversification in Eurasia: insights from the mesobuthid scorpions

The aridification from Middle Miocene onwards has transformed the Asian interior into an arid environment, and the Pleistocene glacial-interglacial oscillations exerted further ecological impact. Therefore, both aridification and glaciation would have considerably influenced the evolution of many mid-latitude species in temperate Asia. Here, we tested this perspective by a phylogeographic study of the mesobuthid scorpions across temperate Asia using one mitochondrial and three nuclear genes. Concordant mitochondrial and nuclear gene trees were obtained, which are consistent with species tree inferred using a Bayesian approach. The age of the most recent common ancestor (MRCA) of all the studied scorpions was estimated to be 12.49 Ma (late Middle Miocene); Mesobuthus eupeus diverged from the clade composing Mesobuthus caucasicus and Mesobuthus martensii in early Late Miocene (10.21 Ma); M. martensii diverged from M. caucasicus at 5.53 Ma in Late Miocene. The estimated MRCA ages of M. martensii and the Chinese lineage of M. eupeus were 2.37 and 0.68 Ma, respectively. Central Asia was identified as the ancestral area for the lineage leading to M. martensii and M. caucasicus and the Chinese lineage of M. eupeus. The ancestral habitat of the genus Mesobuthus is likely to have been characterized by an arid environment; a shift towards more humid habitat occurred in the MRCA of M. martensii and a lineage of M. caucasicus, finally leading to the adaptation of M. martensii to humid environment. Our data strongly support the idea that the stepwise intensified aridifications from Mid-Miocene onwards drove the diversification of mesobuthid scorpions, and suggest that M. martensii and M. eupeus observed today in China originated from an ancestral lineage distributed in Central Asia. Both the colonization and the ensuing evolution of these species in East Asia appear to have been further moulded by Quaternary glaciations.

Phylogeographical analysis of mtDNA data indicates postglacial expansion from multiple glacial refugia in woodland caribou (Rangifer tarandus caribou)

Glacial refugia considerably shaped the phylogeographical structure of species and may influence intra-specific morphological, genetic, and adaptive differentiation. However, the impact of the Quaternary ice ages on the phylogeographical structure of North American temperate mammalian species is not well-studied. Here, we surveyed ~1600 individuals of the widely distributed woodland caribou (Rangifer tarandus caribou) using mtDNA control region sequences to investigate if glacial refugia contributed to the phylogeographical structure in this subspecies. Phylogenetic tree reconstruction, a median-joining network, and mismatch distributions supported postglacial expansions of woodland caribou from three glacial refugia dating back to 13544-22005 years. These three lineages consisted almost exclusively of woodland caribou mtDNA haplotypes, indicating that phylogeographical structure was mainly shaped by postglacial expansions. The putative centres of these lineages are geographically separated; indicating disconnected glacial refugia in the Rocky Mountains, east of the Mississippi, and the Appalachian Mountains. This is in congruence with the fossil record that caribou were distributed in these areas during the Pleistocene. Our results suggest that the last glacial maximum substantially shaped the phylogeographical structure of this large mammalian North American species that will be affected by climatic change. Therefore, the presented results will be essential for future conservation planning in woodland caribou.

Phylogeographical patterns shed light on evolutionary process in South America

The South American continent is composed of several biogeographical regions harbouring the highest biodiversity on the globe, encompassing five of the world's biodiversity 'hot spots'. Nonetheless, the patterns and processes responsible for shaping its astonishing species diversity are largely unknown. Here, we present a review of current South American phylogeographical knowledge based on published articles on this topic. An appraisal of the literature reveals emerging phylogeographical patterns in the biota of South America. The striking phylogeographical divergence observed among organism lineages in South American studies is suggestive of high levels of undocumented species diversity. The interplay between Pleistocene climatic oscillations and Pliocene/Miocene orogenic events has contributed to shaping the current diversity and distribution of modern lineages in both the tropical and temperate regions of South America. Although older divergence times were observed for a range of species, most herpetofauna underwent an intraspecific lineage split much earlier than other organisms. The geographical ranges of species associated with forest habitats were reduced mainly during glacial cycles, whereas species associated with open vegetation domains have shown variable responses to climatic oscillations. The results suggest a highly complex mosaic of phylogeographical patterns in South America. We suggest future research directions to promote a better understanding of the origin and maintenance of the South American biota.

Parasites in ungulates of Arctic North America and Greenland: a view of contemporary diversity, ecology, and impact in a world under change

Parasites play an important role in the structure and function of arctic ecosystems, systems that are currently experiencing an unprecedented rate of change due to various anthropogenic perturbations, including climate change. Ungulates such as muskoxen, caribou, moose and Dall's sheep are also important components of northern ecosystems and are a source of food and income, as well as a focus for maintenance of cultural traditions, for northerners. Parasites of ungulates can influence host health, population dynamics and the quality, quantity and safety of meat and other products of animal origin consumed by people. In this article, we provide a contemporary view of the diversity of nematode, cestode, trematode, protozoan and arthropod parasites of ungulates in arctic and subarctic North America and Greenland. We explore the intricate associations among host and parasite assemblages and identify key issues and gaps in knowledge that emerge in a regime of accelerating environmental transition.

Hybridization among Arctic white-headed gulls (Larus spp.) obscures the genetic legacy of the Pleistocene

We studied the influence of glacial oscillations on the genetic structure of seven species of white-headed gull that breed at high latitudes (Larus argentatus, L. canus, L. glaucescens, L. glaucoides, L. hyperboreus, L. schistisagus, and L. thayeri). We evaluated localities hypothesized as ice-free areas or glacial refugia in other Arctic vertebrates using molecular data from 11 microsatellite loci, mitochondrial DNA (mtDNA) control region, and six nuclear introns for 32 populations across the Holarctic. Moderate levels of genetic structure were observed for microsatellites (F_ST= 0.129), introns (Φ_ST= 0.185), and mtDNA control region (Φ_ST= 0.461), with among-group variation maximized when populations were grouped based on subspecific classification. Two haplotype and at least two allele groups were observed across all loci. However, no haplotype/allele group was composed solely of individuals of a single species, a pattern consistent with recent divergence. Furthermore, northernmost populations were not well differentiated and among-group variation was maximized when L. argentatus and L. hyberboreus populations were grouped by locality rather than species, indicating recent hybridization. Four populations are located in putative Pleistocene glacial refugia and had larger τ estimates than the other 28 populations. However, we were unable to substantiate these putative refugia using coalescent theory, as all populations had genetic signatures of stability based on mtDNA. The extent of haplotype and allele sharing among Arctic white-headed gull species is noteworthy. Studies of other Arctic taxa have generally revealed species-specific clusters as well as genetic structure within species, usually correlated with geography. Aspects of white-headed gull behavioral biology, such as colonization ability and propensity to hybridize, as well as their recent evolutionary history, have likely played a large role in the limited genetic structure observed.

The impact of Pleistocene climate change on an ancient arctic-alpine plant: multiple lineages of disparate history in Oxyria digyna

The ranges of arctic-alpine species have shifted extensively with Pleistocene climate changes and glaciations. Using sequence data from the trnH-psbA and trnT-trnL chloroplast DNA spacer regions, we investigated the phylogeography of the widespread, ancient (>3 million years) arctic-alpine plant Oxyria digyna (Polygonaceae). We identified 45 haplotypes and six highly divergent major lineages; estimated ages of these lineages (time to most recent common ancestor, T(MRCA)) ranged from ∼0.5 to 2.5 million years. One lineage is widespread in the arctic, a second is restricted to the southern Rocky Mountains of the western United States, and a third was found only in the Himalayan and Altai regions of Asia. Three other lineages are widespread in western North America, where they overlap extensively. The high genetic diversity and the presence of divergent major cpDNA lineages within Oxyria digyna reflect its age and suggest that it was widespread during much of its history. The distributions of individual lineages indicate repeated spread of Oxyria digyna through North America over multiple glacial cycles. During the Last Glacial Maximum it persisted in multiple refugia in western North America, including Beringia, south of the continental ice, and within the northern limits of the Cordilleran ice sheet. Our data contribute to a growing body of evidence that arctic-alpine species have migrated from different source regions over multiple glacial cycles and that cryptic refugia contributed to persistence through the Last Glacial Maximum.

One species, three Pleistocene evolutionary histories: phylogeography of the Italian crested newt, Triturus carnifex

Phylogeographic patterns of temperate species from the Mediterranean peninsulas have been investigated intensively. Nevertheless, as more phylogeographies become available, either unique patterns or new lines of concordance continue to emerge, providing new insights on the evolution of regional biotas. Here, we investigated the phylogeography and evolutionary history of the Italian crested newt, Triturus carnifex, through phylogenetic, molecular dating and population structure analyses of two mitochondrial gene fragments (ND2 and ND4; overall 1273 bp). We found three main mtDNA lineages having parapatric distribution and estimated divergence times between Late Pliocene and Early Pleistocene. One lineage (S) was widespread south of the northern Apennine chain and was further geographically structured into five sublineages, likely of Middle Pleistocene origin. The second lineage (C) was widespread throughout the Padano-Venetian plain and did not show a clear phylogeographic structure. The third lineage (N) was observed in only two populations located on western Croatia/Slovenia. Results of analysis of molecular variance suggested that partitioning populations according to the geographic distribution of these lineages and sublineages explains 76% of the observed genetic variation. The phylogeographic structure observed within T. carnifex and divergence time estimates among its lineages, suggest that responses to Pleistocene environmental changes in this single species have been as diverse as those found previously among several codistributed temperate species combined. Consistent with the landscape heterogeneity, physiographic features, and palaeogeographical evolution of its distribution range, these responses encompass multiple refugia along the Apennine chain, lowland refugia in large peri-coastal plains, and a 'cryptic' northern refugium.

Parallel and lineage-specific molecular adaptation to climate in boreal black spruce

In response to selective pressure, adaptation may follow different genetic pathways throughout the natural range of a species due to historical differentiation in standing genetic variation. Using 41 populations of black spruce (Picea mariana), the objectives of this study were to identify adaptive genetic polymorphisms related to temperature and precipitation variation across the transcontinental range of the species, and to evaluate the potential influence of historical events on their geographic distribution. Population structure was first inferred using 50 control nuclear markers. Then, 47 candidate gene SNPs identified in previous genome scans were tested for relationship with climatic factors using an F(ST) -based outlier method and regressions between allele frequencies and climatic variations. Two main intraspecific lineages related to glacial vicariance were detected at the transcontinental scale. Within-lineage analyses of allele frequencies allowed the identification of 23 candidate SNPs significantly related to precipitation and/or temperature variation, among which seven were common to both lineages, eight were specific to the eastern lineage and eight were specific to the western lineage. The implication of these candidate SNPs in adaptive processes was further supported by gene functional annotations. Multiple evidences indicated that the occurrence of lineage-specific adaptive SNPs was better explained by selection acting on historically differentiated gene pools rather than differential selection due to heterogeneity of interacting environmental factors and pleiotropic effects. Taken together, these findings suggest that standing genetic variation of potentially adaptive nature has been modified by historical events, hence affecting the outcome of recent selection and leading to different adaptive routes between intraspecific lineages.

Genetic signals of demographic expansion in Downy Woodpecker (Picoides pubescens) after the last North American glacial maximum

The glacial cycles of the Pleistocene have been recognized as important, large-scale historical processes that strongly influenced the demographic patterns and genetic structure of many species. Here we present evidence of a postglacial expansion for the Downy Woodpecker (Picoides pubescens), a common member of the forest bird communities in North America with a continental distribution. DNA sequences from the mitochondrial tRNA-Lys, and ATPase 6 and 8 genes, and microsatellite data from seven variable loci were combined with a species distribution model (SDM) to infer possible historical scenarios for this species after the last glacial maximum. Analyses of Downy Woodpeckers from 23 geographic areas suggested little differentiation, shallow genealogical relationships, and limited population structure across the species’ range. Microsatellites, which have higher resolution and are able to detect recent differences, revealed two geographic groups where populations along the eastern edge of the Rocky Mountains (Montana, Utah, Colorado, and southern Alberta) were genetically isolated from the rest of the sampled populations. Mitochondrial DNA, an important marker to detect historical patterns, recovered only one group. However, populations in Idaho and southeast BC contained high haplotype diversity and, in general were characterized by the absence of the most common mtDNA haplotype. The SDM suggested several areas in the southern US as containing suitable Downy Woodpecker habitat during the LGM. The lack of considerable geographic structure and the starburst haplotype network, combined with several population genetic tests, suggest a scenario of demographic expansion during the last part of Pleistocene and early Holocene.

Large-scale asymmetric introgression of cytoplasmic DNA reveals Holocene range displacement in a North American boreal pine complex

Jack pine (Pinus banksiana) and lodgepole pine (Pinus contorta var. latifolia) are two North American boreal hard pines that hybridize in their zone of contact in western Canada. The main objective of this study was to characterize their patterns of introgression resulting from past and recent gene flow, using cytoplasmic markers having maternal or paternal inheritance. Mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) diversity was assessed in allopatric populations of each species and in stands from the current zone of contact containing morphological hybrids. Cluster analyses were used to identify genetic discontinuities among groups of populations. A canonical analysis was also conducted to detect putative associations among cytoplasmic DNA variation, tree morphology, and site ecological features. MtDNA introgression was extensive and asymmetric: it was detected in P. banksiana populations from the hybrid zone and from allopatric areas, but not in P. contorta populations. Very weak cpDNA introgression was observed, and only in P. banksiana populations. The mtDNA introgression pattern indicated that central Canada was first colonized by migrants from a P. contorta glacial population located west of the Rocky Mountains, before being replaced by P. banksiana migrating westward during the Holocene. In contrast, extensive pollen gene flow would have erased the cpDNA traces of this ancient presence of P. contorta. Additional evidence for this process was provided by the results of canonical analysis, which indicated that the current cpDNA background of trees reflected recent pollen gene flow from the surrounding dominant species rather than historical events that took place during the postglacial colonization.

Comparative phylogeography of African savannah ungulates

The savannah biome of sub-Saharan Africa harbours the highest diversity of ungulates (hoofed mammals) on Earth. In this review, we compile population genetic data from 19 codistributed ungulate taxa of the savannah biome and find striking concordance in the phylogeographic structuring of species. Data from across taxa reveal distinct regional lineages, which reflect the survival and divergence of populations in isolated savannah refugia during the climatic oscillations of the Pleistocene. Data from taxa across trophic levels suggest distinct savannah refugia were present in West, East, Southern and South-West Africa. Furthermore, differing Pleistocene evolutionary biogeographic scenarios are proposed for East and Southern Africa, supported by palaeoclimatic data and the fossil record. Environmental instability in East Africa facilitated several spatial and temporal refugia and is reflected in the high inter- and intraspecific diversity of the region. In contrast, phylogeographic data suggest a stable, long-standing savannah refuge in the south.

Reticulate evolution on a global scale: a nuclear phylogeny for New World Dryopteris (Dryopteridaceae)

Reticulate, or non-bifurcating, evolution is now recognized as an important phenomenon shaping the histories of many organisms. It appears to be particularly common in plants, especially in ferns, which have relatively few barriers to intra- and interspecific hybridization. Reticulate evolutionary patterns have been recognized in many fern groups, though very few have been studied rigorously using modern molecular phylogenetic techniques in order to determine the causes of the reticulate patterns. In the current study, we examine patterns of branching and reticulate evolution in the genus Dryopteris, the woodferns. The North American members of this group have long been recognized as a classic example of reticulate evolution in plants, and we extend analysis of the genus to all 30 species in the New World, as well as numerous taxa from other regions. We employ sequence data from the plastid and nuclear genomes and use maximum parsimony (MP), maximum likelihood (ML), Bayesian inference (BI), and divergence time analyses to explore the relationships of New World Dryopteris to other regions and to reconstruct the timing and events which may have led to taxa displaying reticulate rather than strictly branching histories. We find evidence for reticulation among both the North and Central/South American groups of species, and our data support a classic hypothesis for reticulate evolution via allopolyploid speciation in the North America taxa, including an extinct diploid progenitor in this group. In the Central and South American species, we find evidence of extensive reticulation involving unknown ancestors from Asia, and we reject deep coalescent processes such as incomplete lineage sorting in favor of more recent intercontinental hybridization and chloroplast capture as an explanation for the origin of the Latin American reticulate taxa.

A continental-wide perspective: the genepool of nuclear encoded ribosomal DNA and single-copy gene sequences in North American Boechera (Brassicaceae)

74 of the currently accepted 111 taxa of the North American genus Boechera (Brassicaceae) were subject to pyhlogenetic reconstruction and network analysis. The dataset comprised 911 accessions for which ITS sequences were analyzed. Phylogenetic analyses yielded largely unresolved trees. Together with the network analysis confirming this result this can be interpreted as an indication for multiple, independent, and rapid diversification events. Network analyses were superimposed with datasets describing i) geographical distribution, ii) taxonomy, iii) reproductive mode, and iv) distribution history based on phylogeographic evidence. Our results provide first direct evidence for enormous reticulate evolution in the entire genus and give further insights into the evolutionary history of this complex genus on a continental scale. In addition two novel single-copy gene markers, orthologues of the Arabidopsis thaliana genes At2g25920 and At3g18900, were analyzed for subsets of taxa and confirmed the findings obtained through the ITS data.

Current and historical drivers of landscape genetic structure differ in core and peripheral salamander populations

With predicted decreases in genetic diversity and greater genetic differentiation at range peripheries relative to their cores, it can be difficult to distinguish between the roles of current disturbance versus historic processes in shaping contemporary genetic patterns. To address this problem, we test for differences in historic demography and landscape genetic structure of coastal giant salamanders (Dicamptodon tenebrosus) in two core regions (Washington State, United States) versus the species' northern peripheral region (British Columbia, Canada) where the species is listed as threatened. Coalescent-based demographic simulations were consistent with a pattern of post-glacial range expansion, with both ancestral and current estimates of effective population size being much larger within the core region relative to the periphery. However, contrary to predictions of recent human-induced population decline in the less genetically diverse peripheral region, there was no genetic signature of population size change. Effects of current demographic processes on genetic structure were evident using a resistance-based landscape genetics approach. Among core populations, genetic structure was best explained by length of the growing season and isolation by resistance (i.e. a ‘flat’ landscape), but at the periphery, topography (slope and elevation) had the greatest influence on genetic structure. Although reduced genetic variation at the range periphery of D. tenebrosus appears to be largely the result of biogeographical history rather than recent impacts, our analyses suggest that inherent landscape features act to alter dispersal pathways uniquely in different parts of the species' geographic range, with implications for habitat management.

Pleistocene and ecological effects on continental-scale genetic differentiation in the bobcat (Lynx rufus)

The potential for widespread, mobile species to exhibit genetic structure without clear geographic barriers is a topic of growing interest. Yet the patterns and mechanisms of structure--particularly over broad spatial scales--remain largely unexplored for these species. Bobcats occur across North America and possess many characteristics expected to promote gene flow. To test whether historical, topographic or ecological factors have influenced genetic differentiation in this species, we analysed 1 kb mtDNA sequence and 15 microsatellite loci from over 1700 samples collected across its range. The primary signature in both marker types involved a longitudinal cline with a sharp transition, or suture zone, occurring along the Great Plains. Thus, the data distinguished bobcats in the eastern USA from those in the western half, with no obvious physical barrier to gene flow. Demographic analyses supported a scenario of expansion from separate Pleistocene refugia, with the Great Plains representing a zone of secondary contact. Substructure within the two main lineages likely reflected founder effects, ecological factors, anthropogenic/topographic effects or a combination of these forces. Two prominent topographic features, the Mississippi River and Rocky Mountains, were not supported as significant genetic barriers. Ecological regions and environmental correlates explained a small but significant proportion of genetic variation. Overall, results implicate historical processes as the primary cause of broad-scale genetic differentiation, but contemporary forces seem to also play a role in promoting and maintaining structure. Despite the bobcat's mobility and broad niche, large-scale landscape changes have contributed to significant and complex patterns of genetic structure.

Phylogeny, divergence times, and historical biogeography of New World Dryopteris (Dryopteridaceae)

PREMISE OF THE STUDY

Dryopteris is a large, cosmopolitan fern genus ideal for addressing questions about diversification, biogeography, hybridization, and polyploidy, which have historically been understudied in ferns. We constructed a highly resolved, well-supported phylogeny for New World Dryopteris and used it to investigate biogeographic patterns and divergence times.

METHODS

We analyzed relationships among 97 species of Dryopteris, including taxa from all major biogeographic regions, with analyses based on 5699 aligned nucleotides from seven plastid loci. Phylogenetic analyses used maximum parsimony, maximum likelihood, and Bayesian inference. We conducted divergence time analyses using BEAST and biogeographic analyses using maximum parsimony, maximum likelihood, Bayesian, and S-DIVA approaches. We explored the monophyly of subgenera and sections in the most recent generic classification and of geographic groups of taxa using Templeton tests.

KEY RESULTS

The genus Dryopteris arose ca. 42 million years ago (Ma). Most of the Central and South American species form a well-supported clade which arose 32 Ma, but the remaining New World species are the result of multiple, independent dispersal and vicariance events involving Asia, Europe, and Africa over the last 15 Myr. We identified six long-distance dispersal events and three vicariance events in the immediate ancestry of New World species; reconstructions for another four lineages were ambiguous.

CONCLUSIONS

New World Dryopteris are not monophyletic; vicariance has dominated the history of the North American species, while long-distance dispersal prevails in the Central and South American species, a pattern not previously seen in plants.

Historical biogeography of the North American glacier ice worm, Mesenchytraeus solifugus (Annelida: Oligochaeta: Enchytraeidae)

North American ice worms are the largest glacially-obligate metazoans, inhabiting coastal, temperate glaciers between southcentral Alaska and Oregon. We have collected ice worm specimens from 10 new populations, completing a broad survey throughout their geographic range. Phylogenetic analyses of 87 individuals using fragments of nuclear 18S rRNA, and mitochondrial 12S rRNA and cyctochrome c oxidase subunit 1 (CO1) identified 18 CO1 haplotypes with divergence values up to ~10%. Phylogeographic interpretations suggest a St. Elias Range, Alaskan ancestry from an aquatic mesenchytraeid oligochaete during the early-Pliocene. A gradual, northward expansion by active dispersal from the central St. Elias clade characterizes a northern clade that is confined to Alaska (with one exception on Vancouver Island, British Columbia), while a distinct southern clade representing worms from British Columbia, Washington and Oregon was likely founded by a passive dispersal event originating from a northern ancestor. The geographic boundary between central and southern clades coincides with an ice worm distribution gap located in southern Alaska, which appears to have restricted active gene flow throughout the species' evolutionary history.