Multiple refugia and barriers explain the phylogeography of the Valais shrew, Sorex antinorii (Mammalia: Soricomorpha), Sorex antinorii (Mammalia: Soricomorpha)

The aim of the present study was to investigate the genetic structure of the Valais shrew (Sorex antinorii) by a combined phylogeographical and landscape genetic approach, and thereby to infer the locations of glacial refugia and establish the influence of geographical barriers. We sequenced part of the mitochondrial cytochrome b (cyt b) gene of 179 individuals of S. antinorii sampled across the entire species' range. Six specimens attributed to S. arunchi were included in the analysis. The phylogeographical pattern was assessed by Bayesian molecular phylogenetic reconstruction, population genetic analyses, and a species distribution modelling (SDM)-based hindcasting approach. We also used landscape genetics (including isolation-by-resistance) to infer the determinants of current intra-specific genetic structure. The phylogeographical analysis revealed shallow divergence among haplotypes and no clear substructure within S. antinorii. The starlike structure of the median-joining network is consistent with population expansion from a single refugium, probably located in the Apennines. Long branches observed on the same network also suggest that another refugium may have existed in the north-eastern part of Italy. This result is consistent with SDM, which also suggests several habitable areas for S. antinorii in the Italian peninsula during the LGM. Therefore S. antinorii appears to have occupied disconnected glacial refugia in the Italian peninsula, supporting previous data for other species showing multiple refugia within southern refugial areas. By coupling genetic analyses and SDM, we were able to infer how past climatic suitability contributed to genetic divergence of populations. The genetic differentiation shown in the present study does not support the specific status of S. arunchi.

Conservation phylogeography: does historical diversity contribute to regional vulnerability in European tree frogs (Hyla arborea)?

Documenting and preserving the genetic diversity of populations, which conditions their long-term survival, have become a major issue in conservation biology. The loss of diversity often documented in declining populations is usually assumed to result from human disturbances; however, historical biogeographic events, otherwise known to strongly impact diversity, are rarely considered in this context. We apply a multilocus phylogeographic study to investigate the late-Quaternary history of a tree frog (Hyla arborea) with declining populations in the northern and western part of its distribution range. Mitochondrial and nuclear polymorphisms reveal high genetic diversity in the Balkan Peninsula, with a spatial structure moulded by the last glaciations. While two of the main refugial lineages remained limited to the Balkans (Adriatic coast, southern Balkans), a third one expanded to recolonize Northern and Western Europe, loosing much of its diversity in the process. Our findings show that mobile and a priori homogeneous taxa may also display substructure within glacial refugia ('refugia within refugia') and emphasize the importance of the Balkans as a major European biodiversity centre. Moreover, the distribution of diversity roughly coincides with regional conservation situations, consistent with the idea that historically impoverished genetic diversity may interact with anthropogenic disturbances, and increase the vulnerability of populations. Phylogeographic models seem important to fully appreciate the risks of local declines and inform conservation strategies.

Multiscale phylogeography: glacial history and autopolyploidy of Biscutella laevigata L. (Brassicaceae) in the Western Alps

Résumé Les changements climatiques du Quaternaire ont eu une influence majeure sur la distribution et l'évolution des biota septentrionaux. Les Alpes offrent un cadre spatio-temporel bien étudié pour comprendre la réactivité de la flore et le potentiel d'adaptation d'une espèce végétale face aux changements climatiques. Certaines hypothèses postulent une diversification des espèces en raison de la disparition complète de la flore des Alpes et d'un isolement important des espèces dans des refuges méridionaux durant les dernières glaciations (Tabula Rasa). Une autre hypothèse stipule le maintien de poches de résistance pour la végétation au coeur des Alpes (Nunataks). Comme de nombreuses espèces végétales présentant un grand succès écologique semblent avoir réagi aux glaciations par la multiplication de leur génome (autopolyploïdie), leur étude en milieu naturel devrait permettre de comprendre les avantages inhérents à la polyploïdie. Biscutella laevigata est un modèle emblématique de biogéographie historique, diverses études ayant montré que des populations diploïdes sont actuellement isolées dans les zones restées déglacées durant le dernier maximum glaciaire, alors que des tétraploïdes ont recolonisé l'ensemble des zones alpines mises à nu par le retrait des glaciers. Si le contexte périglaciaire semble avoir favorisé ce jeune complexe autopolyploïde, les circonstances et les avantages de cette mutation génomique ne sont pas encore clairs. Y a-t-il eu de multiples événements de polyploïdisation ? Dans quelle mesure affecte(nt)il(s) la diversité génétique et le potentiel évolutif des polyploïdes ? Les polyploïdes ont-ils une grande flexibilité génomique, favorisant une radiation adaptative, ou doivent-ils leur succès à une grande plasticité écologique ? Cette étude aborde ces questions à différentes échelles spatiales et temporelles. L'échelle régionale des Alpes occidentales permet d'aborder les facteurs distaux (aspects historiques), alors que l'échelle locale cherche à appréhender les facteurs proximaux (mécanismes évolutifs). Dans les Alpes occidentales, des populations ont été densément échantillonnées et étudiées grâce à (1) leur cytotype, (2) leur appartenance taxonomique, (3) leur habitat et (4) des marqueurs moléculaires de l'ADN chloroplastique, en vue d'établir leurs affinités évolutives. Á l'échelle locale, deux systèmes de population ont été étudiés : l'un où les populations persistent en périphérie de l'aire de distribution et l'autre au niveau du front actif de colonisation, en marge altitudinale. Les résultats à l'échelle des Alpes occidentales révèlent les sites d'intérêt (refuges glaciaires, principales barrières et voies de recolonisation) pour une espèce représentative des pelouses alpines, ainsi que pour la biodiversité régionale. Les Préalpes ont joué un rôle important dans le maintien de populations à proximité immédiate des Alpes centrales et dans l'évolution du taxon, voire de la végétation. Il est aussi démontré que l'époque glaciaire a favorisé l'autopolyploïdie polytopique et la recolonisation des Alpes occidentales par des lignées distinctes qui s'hybrident au centre des Alpes, influençant fortement leur diversité génétique et leur potentiel évolutif. L'analyse de populations locales en situations contrastées à l'aide de marqueurs AFLP montre qu'au sein d'une lignée présentant une grande expansion, la diversité génétique est façonnée par des forces évolutives différentes selon le contexte écologique et historique. Les populations persistant présentent une dispersion des gènes restreinte, engendrant une diversité génétique assez faible, mais semblent adaptées aux conditions locales de l'environnement. À l'inverse, les populations colonisant la marge altitudinale sont influencées par les effets de fondation conjugués à une importante dispersion des gènes et, si ces processus impliquent une grande diversité génétique, ils engendrent une répartition aléatoire des génotypes dans l'environnement. Les autopolyploïdes apparaissent ainsi comme capables de persister face aux changements climatiques grâce à certaines facultés d'adaptation locale et de grandes capacités à maintenir une importante diversité génétique lors de la recolonisation post-glaciaire. Summary The extreme climate changes of the Quaternary have had a major influence on species distribution and evolution. The European Alps offer a great framework to investigate flora reactivity and the adaptive potential of species under changing climate. Some hypotheses postulate diversification due to vegetation removal and important isolation in southern refugia (Tabula Rasa), while others explain phylogeographic patterns by the survival of species in favourable Nunataks within the Alps. Since numerous species have successfully reacted to past climate changes by genome multiplication (autopolyploidy), studies of such taxa in natural conditions is likely to explain the ecological success and the advantages of autopolyploidy. Early cytogeographical surveys of Biscutella laevigata have shed light on the links between autopolyploidy and glaciations by indicating that diploids are now spatially isolated in never-glaciated areas, while autotetraploids have recolonised the zones covered by glaciers- during the last glacial maximum. A periglacial context apparently favoured this young autopolyploid complex but the circumstances and the advantages of this genomic mutation remain unclear. What is the glacial history of the B. laevigata autopolyploid complex? Are there multiple events of polyploidisation? To what extent do they affect the genetic diversity and the evolutionary potential of polyploids? Is recolonisation associated with adaptive processes? How does long-term persistence affect genetic diversity? The present study addresses these questions at different spatiotemporal scales. A regional survey at the Western Alps-scale tackles distal factors (evolutionary history), while local-scale studies explore proximal factors (evolutionary mechanisms). In the Western Alps, populations have been densely sampled and studied from the (1) cytotypic, (2) morphotaxonomic, (3) habitat point of views, as well as (4) plastid DNA molecular markers, in order to infer their relationships and establish the maternal lineages phylogeography. At the local scale, populations persisting at the rear edge and populations recolonising the attitudinal margin at the leading edge have been studied by AFLPs to show how genetic diversity is shaped by different evolutionary forces across the species range. The results at the regional scale document the glacial history of a widespread species, representative of alpine meadows, in a regional area of main interest (glacial refugia, main barriers and recolonisation routes) and points out to sites of interest for regional biodiversity. The external Alps have played a major role in the maintenance of populations near the central Alps during the Last Glacial Maximum and influenced the evolution of the species, and of vegetation. Polytopic autopolyploidy in different biogeographic districts is also demonstrated. The species has had an important and rapid radiation because recolonisation took place from different refugia. The subsequent recolonisation of the Western Alps was achieved by independent lineages that are presently admixing in the central Alps. The role of the Pennic summit line is underlined as a great barrier that was permeable only through certain favourable high-altitude passes. The central Alps are thus viewed as an important crossroad where genomes with different evolutionary histories are meeting and admixing. The AFLP analysis and comparison of local populations growing in contrasted ecological and historical situations indicate that populations persisting in the external Alps present restricted gene dispersal and low genetic diversity but seem in equilibrium with their environment. On the contrary, populations colonising the attitudinal margin are mainly influenced by founder effects together with great gene dispersal and genotypes have a nearly random distribution, suggesting that recolonisation is not associated with adaptive processes. Autopolyploids that locally persist against climate changes thus seem to present adaptive ability, while those that actively recolonise the Alps are successful because of their great capacity to maintain a high genetic diversity against founder effects during recolonisation.

Comparative phylogeography of mutualists and the effect of the host on the genetic structure of its partners.

Whether or not species participating in specialized and obligate interactions display similar and simultaneous demographic variations at the intraspecific level remains an open question in phylogeography. In the present study, we used the mutualistic nursery pollination occurring between the European globeflower Trollius europaeus and its specialized pollinators in the genus Chiastocheta as a case study. Explicitly, we investigated if the phylogeographies of the pollinating flies are significantly different from the expectation under a scenario of plant-insect congruence. Based on a large-scale sampling, we first used mitochondrial data to infer the phylogeographical histories of each fly species. Then, we defined phylogeographical scenarios of congruence with the plant history, and used maximum likelihood and Bayesian approaches to test for plant-insect phylogeographical congruence for the three Chiastocheta species. We show that the phylogeographical histories of the three fly species differ. Only Chiastocheta lophota and Chiastocheta dentifera display strong spatial genetic structures, which do not appear to be statistically different from those expected under scenarios of phylogeographical congruence with the plant. The results of the present study indicate that the fly species responded in independent and different ways to shared evolutionary forces, displaying varying levels of congruence with the plant genetic structure

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data.

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.

Molecular phylogeography of the nose-horned viper (Vipera ammodytes, Linnaeus (1758)): evidence for high genetic diversity and multiple refugia in the Balkan peninsula

The nose-horned viper (Vipera ammodytes) occurs in a large part of the south-eastern Europe and Asia Minor. Phylogenetic relationships were reconstructed for a total of 59 specimens using sequences from three mitochondrial regions (16S and cytochrome b genes, and control region, totalling 2308 bp). A considerable number of clades were observed within this species, showing a large genetic diversity within the Balkan peninsula. Splitting of the basal clades was evaluated to about 4 million years ago. Genetic results are in contradiction with presently accepted taxonomy based on morphological characters: V. a. gregorwallneri and V. a. ruffoi do not display any genetic difference compared with the nominotypic subspecies (V. a. ammodytes), involving that these subspecies can be regarded as synonyms. High genetic divergence in the central part of the Balkan peninsula is not concordant with low morphological differentiation. Finally, the extensive genetic diversity within the Balkan peninsula and the colonisation routes are discussed

Comparative phylogeography and postglacial colonization routes in Europe.

The Quaternary cold periods in Europe are thought to have heavily influenced the amount and distribution of intraspecific genetic variation in both animals and plants. The phylogeographies of 10 taxa, including mammals (Ursus arctos, Sorex spp., Crocidura suaveolens, Arvicola spp.), amphibians (Triturus spp.), arthropods (Chorthippus parallelus), and plants (Abies alba, Picea abies, Fagus sylvatica, Quercus spp.), were analysed to elucidate general trends across Europe. Only a small degree of congruence was found amongst the phylogeographies of the 10 taxa, but the likely postglacial colonization routes exhibit some similarities. A Brooks parsimony analysis produced an unrooted area phylogram, showing that: (i) the northern regions were colonized generally from the Iberic and Balkanic refugia; and (ii) the Italian lineages were often isolated due to the presence of the Alpine barrier. The comparison of colonization routes highlighted four main suture-zones where lineages from the different refugia meet. Some of the intraspecific genetic distances among lineages indicated a prequaternary divergence that cannot be connected to any particular cold period, but are probably related mainly to the date of arrival of each taxon in the European continent. As a consequence, molecular genetics so far appears to be of limited use in dating Quaternary events.

Integrating species distribution models (SDMs) and phylogeography for two species of Alpine Primula.

The major intention of the present study was to investigate whether an approach combining the use of niche-based palaeodistribution modeling and phylo-geography would support or modify hypotheses about the Quaternary distributional history derived from phylogeographic methods alone. Our study system comprised two closely related species of Alpine Primula. We used species distribution models based on the extant distribution of the species and last glacial maximum (LGM) climate models to predict the distribution of the two species during the LGM. Phylogeographic data were generated using amplified fragment length polymorphisms (AFLPs). In Primula hirsuta, models of past distribution and phylogeographic data are partly congruent and support the hypothesis of widespread nunatak survival in the Central Alps. Species distribution models (SDMs) allowed us to differentiate between alpine regions that harbor potential nunatak areas and regions that have been colonized from other areas. SDMs revealed that diversity is a good indicator for nunataks, while rarity is a good indicator for peripheral relict populations that were not source for the recolonization of the inner Alps. In P. daonensis, palaeo-distribution models and phylogeographic data are incongruent. Besides the uncertainty inherent to this type of modeling approach (e.g., relatively coarse 1-km grain size), disagreement of models and data may partly be caused by shifts of ecological niche in both species. Nevertheless, we demonstrate that the combination of palaeo-distribution modeling with phylogeographical approaches provides a more differentiated picture of the distributional history of species and partly supports (P. hirsuta) and partly modifies (P. daonensis and P. hirsuta) hypotheses of Quaternary distributional history. Some of the refugial area indicated by palaeodistribution models could not have been identified with phylogeographic data.

Phylogeography and recolonization of the Swiss Alps by the Valais shrew (Sorex antinorii), inferred with autosomal and sex-specific markers.

Using one male-inherited, one female-inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii) in the Swiss Alps. Bayesian analysis on autosomal microsatellites suggests a clear genetic differentiation between two groups of populations. This geographically based structure is consistent with two separate postglacial recolonization routes of the species into Switzerland from Italian refugia after the last Pleistocene glaciations. Sex-specific markers also confirm genetic structuring among western and eastern areas, since very few haplotypes for either Y chromosome or mtDNA genome are shared between the two regions. Overall, these results suggest that two already well-differentiated genetic lineages colonized the Swiss Alps and came into secondary contact in the Rhône Valley. Low level of admixture between the two lineages is likely explained by the mountainous landscape structure of lateral valleys orthogonal to the main Rhône valley.

Comparative phylogeography in a specific and obligate pollination antagonism.

In specific and obligate interactions the nature and abundance of a given species can have important effects on the survival and population dynamics of associated organisms. In a phylogeographic framework, we therefore expect that the fates of organisms interacting specifically are also tightly interrelated. Here we investigate such a scenario by analyzing the genetic structures of species interacting in an obligate plant-insect pollination lure-and-trap antagonism, involving Arum maculatum (Araceae) and its specific psychodid (Diptera) visitors Psychoda phalaenoides and Psycha grisescens. Because the interaction is asymmetric (i.e., only the plant depends on the insect), we expect the genetic structure of the plant to be related with the historical pollinator availability, yielding incongruent phylogeographic patterns between the interacting organisms.Using insect mtDNA sequences and plant AFLP genome fingerprinting, we inferred the large-scale phylogeographies of each species and the distribution of genetic diversities throughout the sampled range, and evaluated the congruence in their respective genetic structures using hierarchical analyses of molecular variances (AMOVA). Because the composition of pollinator species varies in Europe, we also examined its association with the spatial genetic structure of the plant.Our findings indicate that while the plant presents a spatially well-defined genetic structure, this is not the case in the insects. Patterns of genetic diversities also show dissimilar distributions among the three interacting species. Phylogeographic histories of the plant and its pollinating insects are thus not congruent, a result that would indicate that plant and insect lineages do not share the same glacial and postglacial histories. However, the genetic structure of the plant can, at least partially, be explained by the type of pollinators available at a regional scale. Differences in life-history traits of available pollinators might therefore have influenced the genetic structure of the plant, the dependent organism, in this antagonistic interaction.

Multilocus phylogeography of the common lizard Zootoca vivipara at the Ibero-Pyrenean suture zone reveals lowland barriers and high-elevation introgression.

BACKGROUND: The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure RESULTS: The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France. CONCLUSIONS: The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.

The phylogeography of an alpine leaf beetle: divergence within Oreina elongata spans several ice ages.

The genetic landscape of the European flora and fauna was shaped by the ebb and flow of populations with the shifting ice during Quaternary climate cycles. While this has been well demonstrated for lowland species, less is known about high altitude taxa. Here we analyze the phylogeography of the leaf beetle Oreina elongata from 20 populations across the Alps and Apennines. Three mitochondrial and one nuclear region were sequenced in 64 individuals. Within an mtDNA phylogeny, three of seven subspecies are monophyletic. The species is chemically defended and aposematic, with green and blue forms showing geographic variation and unexpected within-population polymorphism. These warning colors show pronounced east-west geographical structure in distribution, but the phylogeography suggests repeated origin and loss. Basal clades come from the central Alps. Ancestors of other clades probably survived across northern Italy and the northern Adriatic, before separation of eastern, southern and western populations and rapid spread through the western Alps. After reviewing calibrated gene-specific substitution rates in the literature, we use partitioned Bayesian coalescent analysis to date our phylogeography. The major clades diverged long before the last glacial maximum, suggesting that O. elongata persisted many glacial cycles within or at the edges of the Alps and Apennines. When analyzing additional barcoding pairwise distances, we find strong evidence to consider O. elongata as a species complex rather than a single species.

Phylogeography and dispersal in the velvet gecko (Oedura lesueurii), and potential implications for conservation of an endangered snake (Hoplocephalus bungaroides).

ABSTRACT: BACKGROUND: To conserve critically endangered predators, we also need to conserve the prey species upon which they depend. Velvet geckos (Oedura lesueurii) are a primary prey for the endangered broad-headed snake (Hoplocephalus bungaroides), which is restricted to sandstone habitats in southeastern Australia. We sequenced the ND2 gene from 179 velvet geckos, to clarify the lizards' phylogeographic history and landscape genetics. We also analysed 260 records from a longterm (3-year) capture-mark-recapture program at three sites, to evaluate dispersal rates of geckos as a function of locality, sex and body size. RESULTS: The genetic analyses revealed three ancient lineages in the north, south and centre of the species' current range. Estimates of gene flow suggest low dispersal rates, constrained by the availability of contiguous rocky habitat. Mark-recapture records confirm that these lizards are highly sedentary, with most animals moving < 30 m from their original capture site even over multi-year periods. CONCLUSION: The low vagility of these lizards suggests that they will be slow to colonise vacant habitat patches; and hence, efforts to restore degraded habitats for broad-headed snakes may need to include translocation of lizards.

Discordances between phylogenetic and morphological patterns in alpine leaf beetles attest to an intricate biogeographic history of lineages in postglacial Europe.

Pleistocene glacial and interglacial periods have moulded the evolutionary history of European cold-adapted organisms. The role of the different mountain massifs has, however, not been accurately investigated in the case of high-altitude insect species. Here, we focus on three closely related species of non-flying leaf beetles of the genus Oreina (Coleoptera, Chrysomelidae), which are often found in sympatry within the mountain ranges of Europe. After showing that the species concept as currently applied does not match barcoding results, we show, based on more than 700 sequences from one nuclear and three mitochondrial genes, the role of biogeography in shaping the phylogenetic hypothesis. Dating the phylogeny using an insect molecular clock, we show that the earliest lineages diverged more than 1 Mya and that the main shift in diversification rate occurred between 0.36 and 0.18 Mya. By using a probabilistic approach on the parsimony-based dispersal/vicariance framework (MP-DIVA) as well as a direct likelihood method of state change optimization, we show that the Alps acted as a cross-roads with multiple events of dispersal to and reinvasion from neighbouring mountains. However, the relative importance of vicariance vs. dispersal events on the process of rapid diversification remains difficult to evaluate because of a bias towards overestimation of vicariance in the DIVA algorithm. Parallels are drawn with recent studies of cold-adapted species, although our study reveals novel patterns in diversity and genetic links between European mountains, and highlights the importance of neglected regions, such as the Jura and the Balkanic range.