Phylogeography and palaeodistribution modelling in the Patagonian steppe: the case of Mulinum spinosum (Apiaceae)

Aim An integrative study of the endemic, yet ubiquitous, Patagonian shrub Mulinum spinosum (Apiaceae) was performed: (1) to assess the historical processes that influenced its geographical pattern of genetic variation; (2) to test hypotheses of its survival in situ or in glacial refugia during glacial cycles; and (3) to model its extant and palaeoclimatic distributions to assess support for the phylogeographical patterns recovered. Location Chilean and Argentinian Andean region and Patagonian steppe. Methods Chloroplast DNA sequences, trnH–psbA, trnS–trnG and 3′trnV–ndhC, were obtained for 314 individuals of M. spinosum from 71 populations. The haplotype data matrix was analysed using nested clade analysis (NCA) to construct a network. Analysis of molecular variance (AMOVA), spatial analysis of molecular variance (SAMOVA) and neutrality tests were also used to test for genetic structure and range expansion in the species. The present potential geographical distribution of M. spinosum was modelled and projected onto a Last Glacial Maximum (LGM) model. Results Amongst the 29 haplotypes observed, one was widely distributed, but most were restricted to either northern or southern regions. The populations with highest haplotype diversity were found in southern Patagonia, the high Andean region, and northern Patagonia. AMOVA and SAMOVA showed latitudinal structure for Argentinian populations. NCA implied patterns of restricted gene flow or dispersal but with some long‐distance dispersal and also long‐distance colonization and/or past fragmentation. Neutrality tests did not support range expansions. The current distribution model was a fairly good representation of the extant geographical distribution of the species, and the distribution model for the LGM did not show important shifts of the extant range to lower latitudes, except for a shift towards the palaeoseashore. Main conclusions Based on agreement amongst phylogeographical patterns, distribution of genetic variability, equivocal evidence of putative refugia and palaeodistribution modelling, it is probable that glaciations did not greatly affect the distribution of Mulinum spinosum. Our results are consistent with the in situ survival hypothesis, and not with the latitudinal migration of plant communities to avoid adverse climate conditions during Pleistocene glaciations. It is possible that populations of northern Patagonia may have been isolated from the southern ones by the Chubut and Deseado basins.     

Past tree range dynamics in the Iberian Peninsula inferred through phylogeography and palaeodistribution modelling: A review

The forests in the Iberian Peninsula have been strongly influenced by past climatic changes, but reconstructing their historical distributions and dynamics is very difficult due to the complex climatic characteristics and relief of the region. Research disciplines such as phylogeography and species distribution modelling can describe the past range dynamics of individual tree species in relatively great detail and help elucidate how these species have reacted to climatic changes. Here we review phylogeographical and modelling studies from species representative of the major Iberian forest types and attempt to extract general trends from the diversity of individual species histories in the Peninsula. To date most studies focus on geographical range dynamics during the Pleistocene, but an increasing body of evidence shows that some species have also retained the genetic imprints of much more ancient processes. Many widespread Iberian species show a deep and often remarkably clear-cut divide between populations from the Mediterranean and from the Atlantic regions of the Iberian Peninsula, suggesting that both areas have independently sustained viable populations over extended periods. In fact, phylogeographical studies commonly find that species had several glacial refugia across the Iberian Peninsula. On the other hand, distribution models help identifying further suitable areas that could have sustained so far undetected refugia. Such studies are of interest for species conservation, because refugium populations are high-priority targets due to their long-term persistence and unique evolutionary trajectory. Overall, we conclude that palaeoecology, phylogeography and species distribution modelling have a great potential to inform each other because of their complementary perspectives and results. A true integration of these approaches is therefore fundamental for further progress in our understanding of past Iberian environments and the organisms they harboured.     

Phylogeography of the widespread African puff adder (Bitis arietans) reveals multiple Pleistocene refugia in southern Africa

Evidence from numerous Pan‐African savannah mammals indicates that open‐habitat refugia existed in Africa during the Pleistocene, isolated by expanding tropical forests during warm and humid interglacial periods. However, comparative data from other taxonomic groups are currently lacking. We present a phylogeographic investigation of the African puff adder (Bitis arietans), a snake that occurs in open‐habitat formations throughout sub‐Saharan Africa. Multiple parapatric mitochondrial clades occur across the current distribution of B. arietans, including a widespread southern African clade that is subdivided into four separate clades. We investigated the historical processes responsible for generating these phylogeographic patterns in southern Africa using species distribution modelling and genetic approaches. Our results show that interior regions of South Africa became largely inhospitable for B. arietans during glacial maxima, whereas coastal and more northerly areas remained habitable. This corresponds well with the locations of refugia inferred from mitochondrial data using a continuous phylogeographic diffusion model. Analysis of data from five anonymous nuclear loci revealed broadly similar patterns to mtDNA. Secondary admixture was detected between previously isolated refugial populations. In some cases, this is limited to individuals occurring near mitochondrial clade contact zones, but in other cases, more extensive admixture is evident. Overall, our study reveals a complex history of refugial isolation and secondary expansion for puff adders and a mosaic of isolated refugia in southern Africa. We also identify key differences between the processes that drove isolation in B. arietans and those hypothesized for sympatric savannah mammals.     

Phylogeography of Potentilla fruticosa (Rosaceae) in the Qinghai-Tibetan Plateau revisited: a reappraisal and new insights

Background: Most empirical studies to date have suggested that alpine plants in the central Qinghai-Tibetan Plateau (QTP) originated from the south-eastern QTP. However, previous phylogeographcial analysis of Potentilla fruticosa suggested that it diversified in the central QTP, which has also been assumed to have been a glacial refugium, and expanded to the north-eastern regions during the mid-Holocene period.

Aim: We reappraise this claim in the light of information acquired from sampling, haplotyping and analysing additional populations.

Methods: We sequenced the chloroplast trnT-L intergenic spacer of 346 individuals collected from 52 populations. We reconstructed phylogenetic relationships among haplotypes through Bayesian and Network analysis, and used several methods for demographic analysis.

Results: Some common and highly divergent haplotypes were distributed in the central, western and south-eastern QTP. Mismatch and other population genetic analyses applied to our data suggested that a distinct range expansion had occurred between approximately 30,000 and 325,000 years ago (kya).

Conclusions: Besides the possible postglacial expansion in the north-eastern QTP, our results further suggested that the range of the species radically expanded across the whole Plateau before the last glacial maximum (around 18 kya) and the south-eastern, western and central QTP regions had together provided important refugia during recent glacial stages. These findings, in contrast to the previous conclusions, highlight the importance of adequate sampling strategies in phylogeographical studies.     

Phylogeography and ecological niche modelling of the New Zealand stick insect Clitarchus hookeri (White) support survival in multiple coastal refugia

Aim Increasing our understanding of the effects of the Last Glacial Maximum (LGM) and determining the location of refugia requires studies on widely distributed species with dense sampling of populations. We have reconstructed the biogeographic history of Clitarchus hookeri (White), a widespread species of New Zealand stick insect that exhibits geographic parthenogenesis, using phylogeographic analysis and ecological niche modelling. Location New Zealand. Methods We used DNA sequence data from the mitochondrial cytochrome c oxidase subunit I gene to reconstruct phylogenetic relationships among haplotypes from C. hookeri and two undescribed Clitarchus species. We also used distribution data from our own field surveys and museum records to reconstruct the geographic distribution of C. hookeri during the present and the LGM, using ecological niche modelling. Results The ecological niche models showed that the geographic distribution of C. hookeri has expanded dramatically since the LGM. Our model predicted large areas of suitable LGM habitat in upper North Island, and small patches along the east coast of South Island. The phylogeographic analysis shows that populations in the northern half of North Island contain much higher levels of genetic variation than those from southern North Island and South Island, and is congruent with the ecological niche model. The distribution of bisexual populations is also non-random, with males completely absent from South Island and very rare in southern North Island. Main conclusions During the LGM C. hookeri was most likely restricted to several refugia in upper North Island and one or more smaller refugia along the east coast of South Island. The unisexual populations predominate in post-glacial landscapes and are clearly favoured in the recolonization of such areas. Our study exemplifies the utility of integrating ecological niche modelling and phylogeographic analysis.     

Phylogeography of the North American red fox: vicariance in Pleistocene forest refugia

Fossil, archaeological, and morphometric data suggest that indigenous red foxes in North America were derived from vicariance in two disjunct refugia during the last glaciation: one in Beringia and one in the contiguous USA. To test this hypothesis, we conducted a phylogeographical analysis of the North American red fox within its presettlement range. We sequenced portions of the mitochondrial cytochrome b (354 bp) gene and D-loop (342 bp) from 220 historical red fox specimens. Phylogenetic analysis of the cytochrome b gene produced two clades that diverged c . 400 000 years before present ( bp ): a Holarctic and a Nearctic clade. D-loop analyses of the Nearctic clade indicated three distinct subclades (≥ 99% Bayesian posterior probability); two that were more recently derived (rho estimate c . 20 000 bp ) and were restricted to the southwestern mountains and the eastern portion of North America, and one that was older (rho estimate c . 45 000 bp ) and more widespread in North America. Populations that migrated north from the southern refugium following deglaciation were derived from the colonization of North America during or prior to the Illinoian glaciation (300 000–130 000 bp ), whereas populations that migrated south from the northern refugium represent a more recent colonization event during the Wisconsin glaciation (100 000–10 000 bp ). Our findings indicate that Nearctic clade red foxes are phylogenetically distinct from their Holarctic counterparts, and reflect long-term isolation in two disjunct forest refugia during the Pleistocene. The montane lineage, which includes endangered populations, may be ecologically and evolutionarily distinct.     

Phylogeography of Schizopygopsis stoliczkai (Cyprinidae) in Northwest Tibetan Plateau area

Schizopygopsis stoliczkai (Cyprinidae, subfamily Schizothoracinae) is one of the major freshwater fishes endemic to the northwestern margin of the Tibetan Plateau. In the current study, we used mitochondrial DNA markers cytochrome b (Cyt b) and 16S rRNA (16S), as well as the nuclear marker, the second intron of the nuclear beta‐actin gene (Act2), to uncover the phylogeography of S. stoliczkai. In total, we obtained 74 haplotypes from 403 mitochondrial concatenated sequences. The mtDNA markers depict the phylogenetic structures of S. stoliczkai, which consist of clade North and clade South. The split time of the two clades is dated back to 4.27 Mya (95% HPD = 1.96–8.20 Mya). The estimated split time is earlier than the beginning of the ice age of Pleistocene (2.60 Mya), suggesting that the northwestern area of the Tibetan Plateau probably contain at least two glacial refugia for S. stoliczkai. SAMOVA supports the formation of four groups: (i) the Karakash River group; (ii) The Lake Pangong group; (iii) the Shiquan River group; (iv) the Southern Basin group. Clade North included Karakash River, Lake Pangong, and Shiquan River groups, while seven populations of clade South share the haplotypes. Genetic diversity, star‐like network, BSP analysis, as well as negative neutrality tests indicate recent expansions events of S. stoliczkai. Conclusively, our results illustrate the phylogeography of S. stoliczkai, implying the Shiquan River is presumably the main refuge for S. stoliczkai.     

Phylogeography and demography of Guianan harlequin toads (Atelopus): diversification within a refuge

We investigated the genetic structure of populations of Guianan harlequin toads (genus Atelopus) and their evolutionary affinities to extra-Guianan congeners. Phylogenetic analysis of mitochondrial cytochrome b (cyt b) and NADH dehydrogenase subunit 2 (ND2) gene sequences produced well-supported clades largely corresponding to the four recognized taxa in the Guianas (Atelopus spumarius hoogmoedi, Atelopus spumarius barbotini, Atelopus franciscus, and Atelopus flavescens). Our findings suggest that the Guianan A. spumarius represent distinct evolutionary lineages that merit distinction from Amazonian conspecifics, and that the status of A. flavescens and A. franciscus is somewhat less clear. Approximately 69% of the observed genetic variation is accounted for by differences between these four recognized taxa. Coalescent-based estimates of gene flow between taxa suggest that these lineages are largely isolated from one another. Negligible rates of migration between populations and significant divergence within such close proximity suggests that although the region inhabited by these taxa is almost entirely undisturbed, significant habitat heterogeneity exists as to have produced a remarkable diversification of Atelopus within the eastern Guiana Shield. These results contradict the commonly held view of the Guiana Shield as a 'refuge' whose stability during late Tertiary and Quaternary climatic fluctuations served as a biotic reservoir. Instead, we provide evidence that climatic fluctuations during this time had a diversifying effect within the Guianan region.     

Repeated Range Expansion and Glacial Endurance of Potentilla glabra (Rosaceae) in the Qinghai-Tibetan Plateau

To date, little is still known about how alpine species occurring in the Qinghai-Tibetan Plateau (QTP) responded to past climatic oscillations. Here, by using variations of the chloroplast trnJ-L, we examined the genetic distribution pattern of 101 individuals of Potentilla glabra, comprising both the interior QTP and the plateau edge. Phylogenetic and network analyses of 31 recovered haplotypes identified three tentative clades (A, B and C). Analysis of molecular variance (amova) revealed that most of the genetic variability was found within populations (0.693), while differentiations between populations were obviously distinct (Fst = 0.307). Two independent range expansions within clades A and B occurring at approximately 316 and 201 thousand years ago (kya) were recovered from the hierarchical mismatch analysis, and these two expansions were also confirmed by Fu's Fs values and 'g' tests. However, distant distributions of clade C and private haplotypes from clades A and B suggest that they had survived the Last Glacial Maximum (LGM) and previous glaciers in situ since their origins. Our findings based on available limited samples support that multiple refugia of a few cold-enduring species had been maintained in the QTP platform during LGM and/or previous glacial stages.     

Species diversification and phylogeographical patterns of birds in response to the uplift of the Qinghai-Tibet Plateau and Quaternary glaciations

The Qinghai-Tibet Plateau （QTP） is well known for being of great importance in the evolution of montane species due to its unique geological history and landform configuration, climate complexity, and diversified habitats. The effect of environmental changes since the Quaternary on species diversification, population genetic structure, and demography under environmental change can be studied using phylogenetic and phylogeographieal approaches. Birds are the most well-studied group of all terrestrial vertebrates with regard to their response to climatic changes over time. Herein, we briefly review the species diversification of birds in response to the uplift of the QTP, focusing on summarizing the different phylogeographical patterns of birds on the Plateau, its southeastern margin, and the Eastern Himalayas and the reasons underlying these patterns. Speciation was found to be closely related to the uplift of the QTP, with different patterns of intraspecific processes： （1） no divergence within a single refuge was identified in a restricted semi-continuous area of the eastern margin of the Plateau; （2） two divergent lineages with separated refugia were located at the south-eastern and north-eastern margins of the plateau; and （3） multiple divergent lineages within subregions were found in the Eastern Himalayas. Glacial movements and induced climate change are considered to be key factors in shaping these different patterns. The species distributed mainly in the heavily ice-covered platform regions of the Plateau experienced population expansion following the retreat of the extensive glaciations, whereas the species distributed on the ice-free edges of the plateau maintained their population size at a stable level. Demographic stresses on the edge species might have been mitigated by the milder climate in comparison to their platform-distributed counterparts. Various behavioral and ecological characteristics, including dispersal capacity, habitat preference, and elevation specificity, along with evolutionary history might have helped to shape these different phylogeographical patterns [Current Zoology 60 （2）： 149-161, 2014].     

Phylogeography of the endangered Eryngium alpinum L. (Apiaceae) in the European Alps

We studied the phylogeography of Eryngium alpinum by sequencing two intergenic chloroplast spacers, trnH-psbA and trnS-trnG (1322 bp). The sampling design included 36 populations and 397 individuals spanning the entire distribution range of the species, from France to Bosnia. Twenty-one haplotypes were characterized and polymorphism was observed both within and among populations. Population differentiation was strong (F(ST) = 0.92) and largely explained by the distinction of five geographic regions: Southwestern, Western, Middle, Eastern Alps and Balkans (F(CT) = 0.62). Moreover, N(ST) was significantly higher than G(ST) (P < 0.05), showing the existence of a phylogeographic pattern. Six major lineages were recognized using samova and median-joining networks. One lineage, highly divergent from the other ones, was only found in the Balkans and probably persisted in situ during last glaciations. All other lineages might have survived in a Southwestern refugium (Mercantour) and colonized the entire Alpine arc (Southwestern, Western, Middle and Eastern Alps) through repeated colonization events at different time periods. This is the first empirical study suggesting Southern refugia for calcareous Alpine plants, although the existence of a secondary refugium in northern Italy/Austria is also suspected. We also observed recent haplotype diversification, especially in the Southwestern Alps.     

Phylogeography of the Percichthyidae (Pisces) in Patagonia: roles of orogeny, glaciation, and volcanism

We used molecular evidence to examine the roles that vicariance mechanisms (mountain-building and drainage changes during the Pleistocene) have played in producing phylogeographical structure within and among South American fish species of the temperate perch family Percichthyidae. The percichthyids include two South American genera, Percichthys and Percilia, each containing several species, all of which are endemic to southern Argentina and Chile (Patagonia). Maximum-likelihood phylogenies constructed using mitochondrial DNA (mtDNA) control region haplotypes and nuclear GnRH3-2 intron allele sequences support the current taxonomy at the genus level (both Percichthys and Percilia form strongly supported, monophyletic clades) but indicate that species-level designations need revision. Phylogeographical patterns at the mtDNA support the hypothesis that the Andes have been a major barrier to gene flow. Most species diversity occurs in watersheds to the west of the Andes, together with some ancient divergences among conspecific populations. In contrast, only one species (Percichthys trucha) is found east of the Andes, and little to no phylogeographical structure occurs among populations in this region. Mismatch analyses of mtDNA sequences suggest that eastern populations last went through a major bottleneck c. 188 000 bp, a date consistent with the onset of the penultimate and largest Pleistocene glaciation in Patagonia. We suggest that eastern populations have undergone repeated founder-flush events as a consequence of glacial cycles, and that the shallow phylogeny is due to mixing during recolonization periods. The area of greater diversity west of the Andes lies outside the northern limit of the glaciers. mtDNA mismatch analysis of the genus Percilia which is restricted to this area suggests a long-established population at equilibrium. We conclude that patterns of genetic diversity in these South American genera have been primarily influenced by barriers to gene flow (Andean orogeny, and to a lesser extent, isolation in river drainages), and by glacial cycles, which have resulted in population contraction, re-arrangement of some watersheds, and the temporary breakdown of dispersal barriers among eastern river systems.     

Phylogeography of the widespread plant Ailanthus altissima (Simaroubaceae) in China indicated by three chloroplast DNA regions

Ailanthus altissima (Mill.) Swingle, a temperate tree species, has a wide distribution in China. To infer its refugia and patterns of migration during past climatic changes in China, genetic variations among different populations were studied. Gene sequences of three chloroplast DNA spacer regions, psbA-trnH, trnL-trnF, andtrnD-trnT, were obtained from 440 individuals of 44 populations. The distribution of haplotype and the relationships among them were investigated by haplotype network. In addition, the genetic diversity of the sampled regions was inferred, and the biogeographic history was also reconstructed. Twelve haplotypes were identified, among which, five were unique. The phylogenetic analysis and geographical distribution of haplotypes indicate that multiple glacial refugia existed in mainland China during the Quaternary oscillations. Due to the combined effects of contiguous range expansion and allopatric fragmentation, significant genetic structure was not found in this study. Based on biogeographic and demographic analysis, three main dispersal routes were identified for the major haplotypes, whereas others were more likely localized demographic expansion.     

Phylogeography of the vairone (Leuciscus souffia,Risso 1826) in Central Europe

The vairone Leuciscus souffia is a cyprinid fish that inhabits river systems in and around the Alps. The complete mitochondrial DNA control region (945 bp) was sequenced in 295 vairone from 22 populations in Central Europe. A total of 51 haplotypes were identified with a maximum pairwise distance between haplotypes of 5.6%. Phylogenetic analyses revealed two major clades in L. souffia, an ‘Italian’ clade, and an ‘Alpine’ clade. Two hybrid zones exist, in the Mediterranean Alps and in the Soca basin. The position of the sister species of L. souffia, L. turskyi, to the ‘Alpine’ and the ‘Italian’ clade could not be resolved unambiguously. However, a linearized tree analysis indicated that L. turskyi represents a third lineage, that originated at the same time as the ‘Alpine’ and the ‘Italian’ clades of L. souffia. In the ‘Alpine’ clade two groups were resolved, a subclade with haplotypes from the Rhône and Var basins and a cluster with haplotypes from the Danube and Rhine systems. Our data suggest a long history of the vairone in Central Europe, predating Pleistocene glacial cycles. Two main refugia during glaciations must have existed, one in Italy and another one most probably in the Danube system. However, age estimates based on molecular clock calibrations suggest the survival of ‘Alpine’ haplotypes in several drainages during the last glaciation cycles. The Rhine system was only recently colonized.     

Multiple periglacial refugia in the Patagonian steppe and post‐glacial colonization of the Andes: the phylogeography of Calceolaria polyrhiza

Aim We perform a phylogeographical study of an endemic Patagonian herbaceous plant to assess whether geographical patterns of genetic variation correspond to in situ Pleistocene survival or to glacial retreat and post‐glacial expansion. We also seek to determine the locations of potential glacial refugia and post‐glacial colonization routes. Location Southern Andes and Patagonian steppe. Methods We used Calceolaria polyrhiza, a widely distributed Patagonian herbaceous plant that occurs mainly in the understorey of Nothofagus rain forests and in the arid Patagonian steppe, as our model system. The chloroplast intergenic spacer trnH–psbA was sequenced for 590 individuals from 68 populations. Sequence data were analysed using phylogenetic (maximum parsimony, maximum likelihood and Bayesian inference) and population genetic (spatial analyses of molecular variance, mismatch distributions and neutrality tests) methods. Nested clade phylogeographic analyses, and divergence time estimates using a calibrated molecular clock, were also conducted. Results A total of 27 haplotypes identified in the present study clustered into four primary genealogical lineages, revealing three significant latitudinal phylogeographical breaks. The two high Andean lineages probably split first, during the late Miocene, and the Patagonian lineage split around 4 Ma, coincident with the establishment of the Patagonian steppe. Within each haplogroup, major diversification occurred in the Pleistocene. The Patagonian groups show a pattern consistent with a rapid post‐glacial expansion and colonization of the Andean flanks, achieved independently by four lineages. The highest haplotype diversity was found along a longitudinal transect that is remarkably congruent with the limit of the ice‐sheet extension during the Greatest Patagonian Glaciation. A north‐east expansion is evident, which is probably associated with the ‘Arid Diagonal’ fluctuations. Main conclusions Glacial climate fluctuations had a substantial impact on the diversification, distribution and demography of the study species. A scenario of multiple periglacial Pleistocene refugia and subsequent multiple recolonization routes, from eastern Patagonia to the Andean flanks, may explain the phylogeographical patterns observed. However, current genetic structure also preserves the imprints of older events that probably occurred in the Miocene and Pliocene, providing evidence that multiple processes, operating at different spatial and temporal scales, have moulded biodiversity in Patagonia.     

Imprints of multiple glacial refugia in the Pyrenees revealed by phylogeography and palaeodistribution modelling of an endemic spider

Mediterranean mountain ranges harbour highly endemic biota in islandlike habitats. Their topographic diversity offered the opportunity for mountain species to persist in refugial areas during episodes of major climatic change. We investigate the role of Quaternary climatic oscillations in shaping the demographic history and distribution ranges in the spider Harpactocrates ravastellus, endemic to the Pyrenees. Gene trees and multispecies coalescent analyses on mitochondrial and nuclear DNA sequences unveiled two distinct lineages with a hybrid zone around the northwestern area of the Catalan Pyrenees. The lineages were further supported by morphological differences. Climatic niche‐based species distribution models (SDMs) identified two lowland refugia at the western and eastern extremes of the mountain range, which would suggest secondary contact following postglacial expansion of populations from both refugia. Neutrality test and approximate Bayesian computation (ABC) analyses indicated that several local populations underwent severe bottlenecks followed by population expansions, which in combination with the deep population differentiation provided evidence for population survival during glacial periods in microrefugia across the mountain range, in addition to the main Atlantic and Mediterranean (western and eastern) refugia. This study sheds light on the complexities of Quaternary climatic oscillations in building up genetic diversity and local endemicity in the southern Europe mountain ranges.