Disjunct distributions during glacial and interglacial periods in mountain butterflies: Erebia epiphron as an example

The opposed forces of differentiation and admixture across glacial cycles in the butterfly Aglais urticae

Glacial cycles lead to periodic population interbreeding and isolation in warm-adapted species, which impact genetic structure and evolution. However, the effects of these processes on highly mobile and more cold-tolerant species are not well understood. This study aims to shed light on the phylogeographic history of Aglais urticae, a butterfly species with considerable dispersal ability, and a wide Palearctic distribution reaching the Arctic. Through the analysis of genomic data, four main genetic lineages are identified: European, Sierra Nevada, Sicily/Calabria/Peloponnese, and Eastern. The results indicate that the Sardo-Corsican endemic taxon ichnusa is a distinct species. The split between the relict lineages in southern Europe and the main European lineage is estimated to have happened 400-450 thousand years ago, with admixture observed during the Quaternary glacial cycles, and still ongoing, albeit to a much smaller extent. These results suggest that these lineages may be better treated as subspecific parapatric taxa. Ecological niche modelling supported the existence of both Mediterranean and extra-Mediterranean refugia during the glacial periods, with the main one located on the Atlantic coast. Nevertheless, gene flow between populations was possible, indicating that both differentiation and admixture have acted continuously across glacial cycles in this cold-tolerant butterfly, generally balancing each other but producing differentiated lineages in the southern peninsulas. We conclude that the population dynamics and the processes shaping the population genetic structure of cold-adapted species during the Quaternary ice ages may be different than those classically accepted for warm-adapted species.

The molecular population structure of Swertia perennis (Gentianaceae) in Central Europe

Phylogeographic analysis of Swertia perennis, a typical European subalpine springtime species, revealed the existence of at least five major phylogenetic lineages. A large phylogeographic separation exists among these geographical regions, which confirms the existence of glacial refugia in the Pyrenees, but also in the Eastern and Central Alps. The results obtained from the analyzes indicate the existence of three major differences between the populations of the Alps and the Pyrenees, but also between the populations of the Alps and those of other geographical regions (Carpathians, southern Czech Republic, Sudetes and northern Poland). The studied populations from the Black Forest and from north-eastern and southern Poland are a relic of the former wider distribution of these (periglacial) genetic groups within Swertia perennis. Our results also confirm the existence of biogeographical links between the Carpathians and the Hercynian Range and the Alps. Certainly, there was an exchange of genes between populations located in the eastern Alps, the Carpathians and the Hercynian ranges (Czech Republic, Jeseníky, Sudetes, Ore Mountains). This confirms previous results of comparative studies on the genetic diversity of populations of other vascular plant species.

Genomic evidence for three distinct species in the Erebia manto complex in Central Europe (Lepidoptera, Nymphalidae)

A problem to implement conservation strategies is that in many cases recognized taxa are in fact complexes of several cryptic species. Failure to properly delineate species may lead to misplaced priorities or to inadequate conservation measures. One such species complex is the yellow-spotted ringlet Erebia manto, which comprises several phenotypically distinct lineages, whose degree of genomic isolation has so far not been assessed. Some of these lineages are geographically restricted and thus possibly represent distinct units with conservation priorities. Using several thousand nuclear genomic markers, we evaluated to which degree the bubastis lineage from the Alps and the vogesiaca lineage from the Vosges, are genetically isolated from the widespread manto lineage. Our results suggest that both lineages are genetically as strongly differentiated from manto as other taxonomically well separated sibling species in this genus from each other, supporting a delineation of bubastis and vogesiaca as independent species. Given the restricted and isolated range of vogesiaca as well as the disjunct distribution of bubastis, our findings have significant implication for future conservation efforts on these formerly cryptic species and highlight the need to investigate the genomic identity within species complexes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-023-01501-w.

Secondary contact rather than coexistence-Erebia butterflies in the Alps

Secondary contact zones are ideal systems to study the processes that govern the evolution of reproductive barriers, especially at advanced stages of the speciation process. An increase in reproductive isolation resulting from selection against maladaptive hybrids is thought to contribute to reproductive barrier buildup in secondary contact zones. Although such processes have been invoked for many systems, it remains unclear to which extent they influence contact zone dynamics in nature. Here, we study a very narrow contact zone between the butterfly species Erebia cassioides and Erebia tyndarus in the Swiss Alps. We quantified phenotypic traits related to wing shape and reproduction as well as ecology to compare the degree of intra- and interspecific differentiation. Even though only very few first-generation hybrids occur, we find no strong indications for current reinforcing selection, suggesting that if reinforcement occurred in our system, it likely operated in the past. Additionally, we show that both species differ less in their ecological niche at the contact zone than elsewhere, which could explain why coexistence between these butterflies may currently not be possible.

Large- and small-scale geographic structures affecting genetic patterns across populations of an Alpine butterfly

Understanding factors influencing patterns of genetic diversity and the population genetic structure of species is of particular importance in the current era of global climate change and habitat loss. These factors include the evolutionary history of a species as well as heterogeneity in the environment it occupies, which in turn can change across time. Most studies investigating spatio-temporal genetic patterns have focused on patterns across wide geographic areas rather than local variation, but the latter can nevertheless be important particularly in topographically complex areas. Here, we consider these issues in the Sooty Copper butterfly (Lycaena tityrus) from the European Alps, using genome-wide SNPs identified through RADseq. We found strong genetic differentiation within the Alps with four genetic clusters, indicating western, central, and eastern refuges, and a strong reduction of genetic diversity from west to east. This reduction in diversity may suggest that the southwestern refuge was the largest one in comparison to other refuges. Also, the high genetic diversity in the west may result from (a) admixture of different western refuges, (b) more recent demographic changes, or (c) introgression of lowland L. tityrus populations. At small spatial scales, populations were structured by several landscape features and especially by high mountain ridges and large river valleys. We detected 36 outlier loci likely under altitudinal selection, including several loci related to membranes and cellular processes. We suggest that efforts to preserve alpine L. tityrus should focus on the genetically diverse populations in the western Alps, and that the dolomite populations should be treated as genetically distinct management units, since they appear to be currently more threatened than others. This study demonstrates the usefulness of SNP-based approaches for understanding patterns of genetic diversity, gene flow, and selection in a region that is expected to be particularly vulnerable to climate change.

Signature of mid-Pleistocene lineages in the European silver fir (Abies alba Mill.) at its geographic distribution margin

In a conservation and sustainable management perspective, we identify the ecological, climatic, and demographic factors responsible for the genetic diversity patterns of the European silver fir (Abies alba Mill.) at its southwestern range margin (Pyrenees Mountains, France, Europe). We sampled 45 populations throughout the French Pyrenees and eight neighboring reference populations in the Massif Central, Alps, and Corsica. We genotyped 1,620 individuals at three chloroplast and ten nuclear microsatellite loci. We analyzed within- and among-population genetic diversity using phylogeographic reconstructions, tests of isolation-by-distance, Bayesian population structure inference, modeling of demographic scenarios, and regression analyses of genetic variables with current and past environmental variables. Genetic diversity decreased from east to west suggesting isolation-by-distance from the Alps to the Pyrenees and from the Eastern to the Western Pyrenees. We identified two Pyrenean lineages that diverged from a third Alpine-Corsica-Massif Central lineage 0.8 to 1.1 M years ago and subsequently formed a secondary contact zone in the Central Pyrenees. Population sizes underwent contrasted changes, with a contraction in the west and an expansion in the east. Glacial climate affected the genetic composition of the populations, with the western genetic cluster only observed in locations corresponding to the coldest past climate and highest elevations. The eastern cluster was observed over a larger range of temperatures and elevations. All demographic events shaping the current spatial structure of genetic diversity took place during the Mid-Pleistocene Transition, long before the onset of the Holocene. The Western Pyrenees lineage may require additional conservation efforts, whereas the eastern lineage is well protected in in situ gene conservation units. Due to past climate oscillations and the likely emergence of independent refugia, east-west oriented mountain ranges may be important reservoir of genetic diversity in a context of past and ongoing climate change in Europe.

Prevalence and relationship of endosymbiotic Wolbachia in the butterfly genus Erebia

BACKGROUND

Wolbachia is an endosymbiont common to most invertebrates, which can have significant evolutionary implications for its host species by acting as a barrier to gene flow. Despite the importance of Wolbachia, still little is known about its prevalence and diversification pattern among closely related host species. Wolbachia strains may phylogenetically coevolve with their hosts, unless horizontal host-switches are particularly common. We address  these issues in the genus Erebia, one of the most diverse Palearctic butterfly genera.

RESULTS

We sequenced the Wolbachia genome from a strain infecting Erebia cassioides and showed that it belongs to the Wolbachia supergroup B, capable of infecting arthropods from different taxonomic orders. The prevalence of Wolbachia across 13 closely related Erebia host species based on extensive population-level genetic data revealed that multiple Wolbachia strains jointly infect all investigated taxa, but with varying prevalence. Finally, the phylogenetic relationships of Wolbachia strains are in some cases significantly associated to that of their hosts, especially among the most closely related Erebia species, demonstrating mixed evidence for phylogenetic coevolution.

CONCLUSIONS

Closely related host species can be infected by closely related Wolbachia strains, evidencing some phylogenetic coevolution, but the actual pattern of infection more often reflects historical or contemporary geographic proximity among host species. Multiple processes, including survival in distinct glacial refugia, recent host shifts in sympatry, and a loss of Wolbachia during postglacial range expansion seem to have jointly shaped the complex interactions between Wolbachia evolution and the diversification of its host among our studied Erebia species.

Reconstruction of forest dynamics in the Western Palaearctic based on phylogeographic analysis of the ringlet butterfly Erebia aethiops

Glacial refugia are centers of high biodiversity. Therefore, knowledge on their locations and reactions of associated populations and landscapes to climatic changes is crucial for conservation management. We here investigated the biogeography of a butterfly species linked to open forest habitats. Using mitochondrial and nuclear markers in combination with Bayesian simulations, we analyzed the location and age of potential glacial refugia of the species. We identified five putative refugia in Europe. Considering the ecological needs of our study species, tree density within these refugial areas, in contrast to earlier assumptions, must have exceeded the level of individually scattered trees. Our results also provide evidence that especially the refuge areas in the Carpathians were previously underestimated regarding their age: the refugia in the Southern Carpathians presented suitable conditions throughout several glacial cycles, probably since the Mindel or Riss cycles. Additionally, our analyses provided support for a forest refugium near the Tatra Mountains persisting the last glacial maximum. Our results underline the usefulness of this and probably other butterfly species as indicators of forest refugia.

Past, current, and potential future distributions of unique genetic diversity in a cold-adapted mountain butterfly

AIM

Climatic changes throughout the Pleistocene have strongly modified species distributions. We examine how these range shifts have affected the genetic diversity of a montane butterfly species and whether the genetic diversity in the extant populations is threatened by future climate change.

LOCATION

Europe.

TAXON

Erebia epiphron Lepidoptera: Nymphalidae.

METHODS

We analyzed mtDNA to map current genetic diversity and differentiation of E. epiphron across Europe to identify population refugia and postglacial range shifts. We used species distribution modeling (SDM) to hindcast distributions over the last 21,000 years to identify source locations of extant populations and to project distributions into the future (2070) to predict potential losses in genetic diversity.

RESULTS

We found substantial genetic diversity unique to specific regions within Europe (total number of haplotypes = 31, number of unique haplotypes = 27, H d = 0.9). Genetic data and SDM hindcasting suggest long-term separation and survival of discrete populations. Particularly, high rates of unique diversity in postglacially colonized sites in England (H d = 0.64) suggest this population was colonized from a now extinct cryptic refugium. Under future climate change, SDMs predict loss of climate suitability for E. epiphron, particularly at lower elevations (<1,000 meters above sea level) equating to 1 to 12 unique haplotypes being at risk under climate scenarios projecting 1°C and 2-3°C increases respectfully in global temperature by 2070.

MAIN CONCLUSIONS

Our results suggest that historical range expansion and retraction processes by a cold-adapted mountain species caused diversification between populations, resulting in unique genetic diversity which may be at risk if distributions of cold-adapted species shrink in future. Assisted colonizations of individuals from at-risk populations into climatically suitable unoccupied habitat might help conserve unique genetic diversity, and translocations into remaining populations might increase their genetic diversity and hence their ability to adapt to future climate change.

Multiple glacial refugia and contemporary dispersal shape the genetic structure of an endemic amphibian from the Pyrenees

Historical factors (colonization scenarios, demographic oscillations) and contemporary processes (population connectivity, current population size) largely contribute to shaping species' present-day genetic diversity and structure. In this study, we use a combination of mitochondrial and nuclear DNA markers to understand the role of Quaternary climatic oscillations and present-day gene flow dynamics in determining the genetic diversity and structure of the newt Calotriton asper (Al. Dugès, 1852), endemic to the Pyrenees. Mitochondrial DNA did not show a clear phylogeographic pattern and presented low levels of variation. In contrast, microsatellites revealed five major genetic lineages with admixture patterns at their boundaries. Approximate Bayesian computation analyses and linear models indicated that the five lineages likely underwent separate evolutionary histories and can be tracked back to distinct glacial refugia. Lineage differentiation started around the Last Glacial Maximum at three focal areas (western, central and eastern Pyrenees) and extended through the end of the Last Glacial Period in the central Pyrenees, where it led to the formation of two more lineages. Our data revealed no evidence of recent dispersal between lineages, whereas borders likely represent zones of secondary contact following expansion from multiple refugia. Finally, we did not find genetic evidence of sex-biased dispersal. This work highlights the importance of integrating past evolutionary processes and present-day gene flow and dispersal dynamics, together with multilocus approaches, to gain insights into what shaped the current genetic attributes of amphibians living in montane habitats.

Dispersal and adaptation strategies of the high mountain butterfly Boloria pales in the Romanian Carpathians

Background

Habitat quality is one main trigger for the persistence of butterflies. The effects of the influencing biotic and abiotic factors may be enhanced by the challenging conditions in high-alpine environments. To better our knowledge in this field, we performed a mark-release-recapture study with Boloria pales in the Southern Carpathians.

Methods

We analysed population structure, movement and foraging behaviour to investigate special adaptations to the alpine environment and to reveal differences between sexes. We compared these aspects in one sector with and one sector without grazing to address the effects of grazing intensity on habitat quality.

Results

We observed "soft" protandry, in which only a small number of males appeared before females, and an extended emergence of individuals over the observed flight period, dividing the population's age structure into three phases; both observations are considered adaptations to high mountain environments. Although both sexes were mostly sedentary, movement differences between them were obvious. Males flew larger distances than females and were more flight-active. This might explain the dimorphism in foraging behaviour: males preferred nectar sources of Asteraceae, females Caprifoliaceae. Transition from the grazed to the ungrazed sector was only observed for males and not for females, but the population density was higher and the flight distances of the individuals were significantly longer on the grazed sector compared with the ungrazed one.

Conclusion

Soft protandry, an extended emergence of the individuals and an adapted behavioural dimorphism between sexes render to represent a good adaptation of B. pales to the harsh environmental conditions of high mountain ecosystems. However, land-use intensity apparently has severe influence on population densities and movement behaviour. To protect B. pales and other high-alpine species from the negative consequences of overgrazing, areas without or just light grazing are needed.

Sexual dimorphism in the alpine butterflies Boloria pales and Boloria napaea: differences in movement and foraging behavior (Lepidoptera: Nymphalidae)

Sexual dimorphism is a widespread phenomenon in Lepidoptera. It is reflected in differences in life history, behavior and morphology. Analyses of differences in behavior are mostly difficult and time-consuming, especially in high mountain ecosystems. To enhance our knowledge on sexual dimorphisms of alpine butterflies, we performed a mark-release-recapture study on 2 species common in the Alps: Boloria pales and Boloria napaea. We analysed movement and foraging behavior to investigate differences between sexes. Both sexes were mostly sedentary and the movement distances of males and females similar. However, obvious differences in dispersal behavior between the sexes were found in the movement patterns. Three different patterns were distinguished. Most males showed intensive flight activity, but mostly flew only in a limited part of the entire habitat (i.e., their individual home range) searching for females, whereas females were less flight active and flew only to find places for oviposition or feeding. The third pattern, where individuals flew larger distances, was only observed in a small number of males, which always returned to their home range. Nearly all feeding was observed on Asteraceae. However, males preferred the genera Leontodon and Crepis, while females preferred Leontodon and Carduus. Apart from this sexual difference in foraging, individuals of both sexes were found to be more or less specialised on nectar sources. Flight activity was generally greater in males than females. Therefore, we think that sex-specific requirements in nectar ingredients exist, that is, sugar for the intensive flight activity of males and amino acids for egg production of females.

Jailed in the mountains: Genetic diversity and structure of an endemic newt species across the Pyrenees

The Pyrenees represent a natural laboratory for biogeographic, evolutionary and ecological research of mountain fauna as a result of the high variety of habitats and the profound effect of the glacial and interglacial periods. There is a paucity of studies providing a detailed insight into genetic processes and better knowledge on the patterns of genetic diversity and how they are maintained under high altitude conditions. This is of particular interest when considering the course of past climate conditions and glaciations in a species which is considered site tenacious, with long generation times. Here we analyzed the genetic patterns of diversity and structure of the endemic Pyrenean brook newt (Calotriton asper) along its distribution range, with special emphasis on the distinct habitat types (caves, streams, and lakes), and the altitudinal and geographical ranges, using a total set of 900 individuals from 44 different localities across the Pyrenean mountain range genotyped for 19 microsatellite loci. We found evidence for a negative longitudinal and positive altitudinal gradient of genetic diversity in C. asper populations. The fact that genetic diversity was markedly higher westwards is in accordance with other Pyrenean species. However, the impact of altitudinal gradient on the genetic diversity seems to differ from other species, and mostly from other amphibians. We found that lower altitudes can act as a barrier probably because the lowlands do not provide a suitable habitat for C. asper. Regarding the distinct habitat types, caves had significantly lower values of genetic diversity compared to streams or lakes. The mean F_ST value was relatively high (0.304) with maximum values as high as 0.771, suggesting a highly structured total population. Indeed, populations were grouped into five subclusters, the eastern populations (cluster 1) remained grouped into two subclusters and the central-western Pyrenees (cluster 2) into three subclusters. The increase of isolation with geographical distance is consistent with the population structure detected. In conclusion, C. asper seems to be adapted to high altitude mountain habitats, and its genetic diversity is higher in the western Pyrenees. In terms of conservation priority, we consider more relevant the populations that represent a reservoir of genetic diversity.

Specialized or opportunistic-how does the high mountain endemic butterfly Erebia nivalis survive in its extreme habitats?

High mountain ecosystems are a challenge for the survival of animal and plant species, which have to evolve specific adaptations to cope with the prevailing extreme conditions. The strategies to survive may reach from opportunistic to highly adapted traits. One species successfully surviving under these conditions is the here studied butterfly Erebia nivalis. In a mark-release-recapture study performed in the Hohe Tauern National Park (Austria) from 22 July to 26 August 2013, we marked 1386 individuals and recaptured 342 of these. For each capture event, we recorded the exact point of capture and various other traits (wing conditions, behavior, nectar sources). The population showed a partial protandrous demography with the minority of males emerging prior to the females, but the majority being synchronized with them. Males and females differed significantly in their behavior with males being more flight active and females nectaring and resting more. Both sexes showed preferences for the same plant species as nectar sources, but this specialization apparently is the result of a rapid individual adaptation to the locally available flowers. Estimates of the realized dispersal distances predicted a comparatively high amount of long-distance flights, especially for females. Therefore, the adaptation of Erebia nivalis to the unpredictable high mountain conditions might be a mixture of opportunism and specialized traits.

How Do Cold-Adapted Plants Respond to Climatic Cycles? Interglacial Expansion Explains Current Distribution and Genomic Diversity in Primula farinosa L

Understanding the effects of past climatic fluctuations on the distribution and population-size dynamics of cold-adapted species is essential for predicting their responses to ongoing global climate change. In spite of the heterogeneity of cold-adapted species, two main contrasting hypotheses have been proposed to explain their responses to Late Quaternary glacial cycles, namely, the interglacial contraction versus the interglacial expansion hypotheses. Here, we use the cold-adapted plant Primula farinosa to test two demographic models under each of the two alternative hypotheses and a fifth, null model. We first approximate the time and extent of demographic contractions and expansions during the Late Quaternary by projecting species distribution models across the last 72 ka. We also generate genome-wide sequence data using a Reduced Representation Library approach to reconstruct the spatial structure, genetic diversity, and phylogenetic relationships of lineages within P. farinosa. Finally, by integrating the results of climatic and genomic analyses in an Approximate Bayesian Computation framework, we propose the most likely model for the extent and direction of population-size changes in $P$. farinosa through the Late Quaternary. Our results support the interglacial expansion of $P$. farinosa, differing from the prevailing paradigm that the observed distribution of cold-adapted species currently fragmented in high altitude and latitude regions reflects the consequences of postglacial contraction processes.

Quaternary history, population genetic structure and diversity of the cold-adapted Alpine newt Ichthyosaura alpestris in peninsular Italy

Mediterranean peninsulas are major biodiversity hotspots, and cold-adapted species are an important component of this biodiversity. However, cold-adapted species contributed surprisingly little to our knowledge of the intimate links between Quaternary environmental changes, species' responses to these changes, and current patterns of intraspecific biodiversity. Here, we investigated the genetic structure and evolutionary history of a cold-adapted amphibian, the Alpine newt Ichthyosaura alpestris, within the Italian peninsula. Nuclear and mitochondrial markers consistently identified three distinct genetic lineages, whose divergence dates to the Early Pleistocene (1.9 and 0.8 million years ago). Our results show that the Italian peninsula provided multiple Pleistocene refugia to this cold-adapted species, and suggest that allopatric fragmentation followed by secondary admixture have been key events in the formation of its current pattern of genetic diversity. Indeed, estimates of population genetic diversity clearly identified contact populations as those achieving the highest levels of diversity. Such concordance among cold-adapted and temperate species in terms of processes triggering the formation of regional patterns of genetic diversity provides strong support for the hypothesis that gene exchange between divergent lineages, rather than long-term stability of refugial populations, has been the main step toward the formation of hotspots of intraspecific biodiversity.

Phylogeography of social polymorphism in a boreo-montane ant

Background

The disjunct distribution of several Palearctic species has been widely shaped by the changes in climatic conditions during the Quaternary. The observed genetic differentiation or reproductive isolation between extant populations may be the outcome of their contemporary geographic separation or reproductive incompatibility due to differences in phenotypic traits which have evolved in isolated refugia. In the boreal ant Leptothorax acervorum, colonies from central and peripheral populations differ in social structure: colonies from Central and Northern Europe may contain several equally reproductive queens (facultative polygyny), while in colonies from peripheral populations in Spain only one the most dominant of several queens lays eggs (functional monogyny). By reconstructing the specie’s evolutionary and demographic history in Southwestern Europe we examine whether variation in social organization is associated with restricted gene flow between the two social forms.

Results

We show that multi-queen colonies from all so far known inner Iberian populations of L. acervorum are functionally monogynous, whereas multi-queen colonies from all Pyrenean populations are polygynous, like those from other previously studied areas in Europe. Our analyses revealed complex spatial-genetic structure, but no association between spatial-genetic structure and social organization in SW-Europe. The population in the western Pyrenees diverged most strongly from other Iberian populations. Moreover, microsatellite data suggest the occurrence of recent bottlenecks in Pyrenean and inner Iberian populations.

Conclusions

Our study shows a lack of reproductive isolation between the two social forms in SW-Europe. This in turn suggests that demographic and spatial patterns in genetic variation as well as the distribution of social phenotypes are better explained by co-variation with climatic, ecological, and historical factors. Moreover, we for the first time show the existence of substantial spatial-genetic structure in L. acervorum, suggesting the existence of multiple refugia in SW-Europe, including two extra-Mediterranean refugia in France. While gene flow among inner Iberian refugia may have been larger during the late glacial, extra-Mediterranean refugia in southern France may have contributed to the post-glacial recolonization of W-Europe.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0711-3) contains supplementary material, which is available to authorized users.

Phylogeny, phylogeography, and evolution in the Mediterranean region: News from a freshwater mussel (Potomida, Unionida)

The Potomida genus (Bivalvia, Unionida) has a Circum-Mediterranean distribution and like other freshwater mussel species, its populations have suffered dramatic declines. Although this genus is currently considered as monotypic, it has a long history of taxonomic revisions and presently many aspects of its systematics and evolutionary history are unclear. We sampled a total of 323 individuals from 39 different sites across the Potomida genus distribution, and sequenced two mitochondrial (16S rDNA and Cytochrome c Oxidase Subunit I) and one nuclear (28S rDNA) genes to clarify its phylogeny and phylogeographic history. Our results show that the genus includes two well-supported clades, one comprising solely the western Mediterranean species Potomida littoralis, and the other including two eastern Mediterranean species, the Greek endemic P. acarnanica and the Anatolian and Middle Eastern P. semirugata. We suggest that Potomida started radiating during the upper Miocene, and that both vicariance and dispersal events shaped the diversification and distribution of the genus along the Mediterranean region. P. littoralis is further divided in two mitochondrial lineages, one restricted to Europe and the other occurring mostly in North Africa. Moreover, some European basins present both lineages in sympatry. The conservation status of the three recognized species should be reevaluated, particularly P. acarnanica, since it is restricted to two Greek river basins presenting a high risk of extinction. Overall, our results clarify some important gaps in knowledge concerning the phylogeny, phylogeography and evolution of the Potomida genus in the Mediterranean region with important taxonomical, ecological and conservational implications.

Three in One--Multiple Faunal Elements within an Endangered European Butterfly Species

Ice ages within Europe forced many species to retreat to refugia, of which three major biogeographic basic types can be distinguished: "Mediterranean", "Continental" and "Alpine / Arctic" species. However, this classification often fails to explain the complex phylogeography of European species with a wide range of latitudinal and altitudinal distribution. Hence, we tested for the possibility that all three mentioned faunal elements are represented within one species. Our data was obtained by scoring 1,307 Euphydryas aurinia individuals (46 European locations) for 17 allozyme loci, and sequencing a subset of 492 individuals (21 sites) for a 626 base pairs COI fragment. Genetic diversity indices, F statistics, hierarchical analyses of molecular variance, individual-based clustering, and networks were used to explore the phylogeographic patterns. The COI fragment represented 18 haplotypes showing a strong geographic structure. All but one allozyme loci analysed were polymorphic with a mean F_ST of 0.20, supporting a pronounced among population structure. Interpretation of both genetic marker systems, using several analytical tools, calls for the recognition of twelve genetic groups. These analyses consistently distinguished different groups in Iberia (2), Italy, Provence, Alps (3), Slovenia, Carpathian Basin, the lowlands of West and Central Europe as well as Estonia, often with considerable additional substructures. The genetic data strongly support the hypothesis that E. aurinia survived the last glaciation in Mediterranean, extra-Mediterranean and perialpine refugia. It is thus a rare example of a model organism that combines attributes of faunal elements from all three of these sources. The observed differences between allozymes and mtDNA most likely result from recent introgression of mtDNA into nuclear allozyme groups. Our results indicate discrepancies with the morphologically-based subspecies models, underlining the need to revise the current taxonomy.

Diversification of the cold-adapted butterfly genus Oeneis related to Holarctic biogeography and climatic niche shifts

Both geographical and ecological speciation interact during the evolution of a clade, but the relative contribution of these processes is rarely assessed for cold-dwelling biota. Here, we investigate the role of biogeography and the evolution of ecological traits on the diversification of the Holarctic arcto-alpine butterfly genus Oeneis (Lepidoptera: Satyrinae). We reconstructed the molecular phylogeny of the genus based on one mitochondrial (COI) and three nuclear (GAPDH, RpS5, wingless) genes. We inferred the biogeographical scenario and the ancestral state reconstructions of climatic and habitat requirements. Within the genus, we detected five main species groups corresponding to the taxonomic division and further paraphyletic position of Neominois (syn. n.). Next, we transferred O. aktashi from the hora to the polixenes species group on the bases of molecular relationships. We found that the genus originated in the dry grasslands of the mountains of Central Asia and dispersed over the Beringian Land Bridges to North America several times independently. Holarctic mountains, in particular the Asian Altai Mts. and Sayan Mts., host the oldest lineages and most of the species diversity. Arctic species are more recent, with Pliocene or Pleistocene origin. We detected a strong phylogenetic signal for the climatic niche, where one lineage diversified towards colder conditions. Altogether, our results indicate that both dispersal across geographical areas and occupation of distinct climatic niches promoted the diversification of the Oeneis genus.

Effects of recent and past climatic shifts on the genetic structure of the high mountain yellow-spotted ringlet butterfly Erebia manto (Lepidoptera, Satyrinae): a conservation problem

Mountain species have evolved important genetic differentiation due to past climatic fluctuations. The genetic uniqueness of many of these lineages is now at risk due to global warming. Here, we analyse allozyme polymorphisms of 1306 individuals (36 populations) of the mountain butterfly Erebia manto and perform Species Distribution Models (SDMs). As a consensus of analyses, we obtained six most likely genetic clusters: (i) Pyrenees with Massif Central; (ii) Vosges; (iii-v) Alps including the Slovakian Carpathians; (vi) southern Carpathians. The Vosges population showed the strongest genetic split from all other populations, being almost as strong as the split between E. manto and its sister species Erebia eriphyle. The distinctiveness of the Pyrenees-Massif Central group and of the southern Carpathians group from all other groups is also quite high. All three groups are assumed to have survived more than one full glacial-interglacial cycle close to their current distributions with up-hill and down-slope shifts conforming climatic conditions. In contrast with these well-differentiated groups, the three groups present in the Alps and the Slovakian Carpathians show a much shallower genetic structure and thus also should be of a more recent origin. As predicted by our SDM projections, rising temperatures will strongly impact the distribution of E. manto. While the populations in the Alps are predicted to shrink, the survival of the three lineages present here should not be at risk. The situation of the three other lineages is quite different. All models predict the extinction of the Vosges lineage in the wake of global warming, and also the southern Carpathians and Pyrenees-Massif Central lineages might be at high risk to disappear. Thus, albeit global warming will therefore be unlikely to threaten E. manto as a species, an important proportion of the species' intraspecific differentiation and thus uniqueness might be lost.

Microscale vicariance and diversification of Western Balkan caddisflies linked to karstification

The karst areas in the Dinaric region of the Western Balkan Peninsula are a hotspot of freshwater biodiversity. Many investigators have examined diversification of the subterranean freshwater fauna in these karst systems. However, diversification of surface-water fauna remains largely unexplored. We assessed local and regional diversification of surface-water species in karst systems and asked whether patterns of population differentiation could be explained by dispersal-diversification processes or allopatric diversification following karst-related microscale vicariance. We analyzed mitochondrial cytochrome c oxidase subunit I (mtCOI) sequence data of 4 caddisfly species (genus Drusus) in a phylogeographic framework to assess local and regional population genetic structure and Pliocene/Pleistocene history. We used BEAST software to assess the timing of intraspecific diversification of the target species. We compared climate envelopes of the study species and projected climatically suitable areas during the last glacial maximum (LGM) to assess differences in the species climatic niches and infer potential LGM refugia. The haplotype distribution of the 4 species (324 individuals from 32 populations) was characterized by strong genetic differentiation with few haplotypes shared among populations (16%) and deep divergence among populations of the 3 endemic species, even at local scales. Divergence among local populations of endemics often exceeded divergence among regional and continental clades of the widespread D. discolor. Major divergences among regional populations dated to 2.0 to 0.5 Mya. Species distribution model projections and genetic structure suggest that the endemic species persisted in situ and diversified locally throughout multiple Pleistocene climate cycles. The pattern for D. discolor was different and consistent with multiple invasions into the region. Patterns of population genetic structure and diversification were similar for the 3 regional endemic Drusus species and consistent with microscale vicariance after the onset of intensified karstification in the Dinaric region. Karstification may induce microscale vicariance of running surface-water habitats and probably promotes allopatric fragmentation of stream insects at small spatial scales.

Out of the alps: the biogeography of a disjunctly distributed mountain butterfly, the almond-eyed ringlet Erebia alberganus (Lepidoptera, Satyrinae)

Many studies on the biogeography of thermophilic and arctic-alpine species were performed during the past. Only little is known about species with intermediate characteristics. We analyzed the molecular biogeography of the butterfly Erebia alberganus (30 populations, representing 1106 individuals), sampled over the Alps, Apennines (Italy), and the Stara Planina (Bulgaria) using allozyme electrophoresis (17 loci). Genetic analyses revealed 3 major splits, with the strongest between the Stara Planina populations and all other populations, and a weaker split between the Alps and the Apennines. Individuals from the Apennines were genetically nested within the Alps group. The Alps cluster was segregated into 3 groups: the Southwestern, Western/Central, and Eastern Alps. The genetic diversities were highest for the Alps populations and significantly lower in the 2 isolates (Apennines, Stara Planina). The remarkable genetic split between Stara Planina and all other populations and the genetic distinctiveness of the former cluster might be interpreted as an ancient colonization event of this Balkan mountain range. The Apennines populations derive from a more recent expansion out of the Southwestern Alps. After surviving the Würm ice age most probably in the central Apennines, accompanied by genetic modification of some of these populations, northward expansion might have started from the western parts of the central Apennines reaching the northern Apennines during the early postglacial. The subtle genetic differentiation found among the Alps populations probably reflects 3 geographically disjunct Würm glacial centers located at the western slopes of the Southwestern Alps, at the southern slopes of the Central Alps, and in the Southeastern Alps.

Multiple glacial refugia of the low-dispersal ground beetle Carabus irregularis: molecular data support predictions of species distribution models

Classical glacial refugia such as the southern European peninsulas were important for species survival during glacial periods and acted as sources of post-glacial colonisation processes. Only recently, some studies have provided evidence for glacial refugia north of the southern European peninsulas. In the present study, we combined species distribution models (SDMs) with phylogeographic analyses (using mitochondrial DNA = mtDNA) to investigate if the cold-adapted, stenotopic and flightless ground beetle species, Carabus irregularis, survived the Last Glacial Maximum (LGM) in classical and/or other refugia. SDMs (for both a western European and for a Carpathian subgroup) were calculated with MAXENT on the basis of 645 species records to predict current and past distribution patterns. Two mtDNA loci (CO1 and ND5, concatenated sequence length: 1785 bp) were analyzed from 91 C. irregularis specimens to reconstruct the phylogeography of Central and eastern European populations and to estimate divergence times of the given lineages. Strong intra-specific genetic differentiation (inter-clade ΦST values ranged from 0.92 to 0.99) implied long-term isolation of major clades and subsclades. The high divergence between the nominate subspecies and the Carpathian subspecies C. i. montandoni points to two independent species rather than subspecies (K-2P distance 0.042 ± 0.004; supposed divergence of the maternal lineages dated back 1.6 to 2.5 million years BP) differing not only morphologically but also genetically and ecologically from each other. The SDMs also inferred classical as well as other refugia for C. irregularis, especially north of the Alps, in southeastern Europe and in the Carpathians. The coincidences between the results of both methods confirm the assumption of multiple glacial refugia for the studied species and the usefulness of combining methodological approaches for the understanding of the history of low-dispersal insect species.

Allopatry as a gordian knot for taxonomists: patterns of DNA barcode divergence in arctic-alpine lepidoptera

Many cold adapted species occur in both montane settings and in the subarctic. Their disjunct distributions create taxonomic complexity because there is no standardized method to establish whether their allopatric populations represent single or different species. This study employs DNA barcoding to gain new perspectives on the levels and patterns of sequence divergence among populations of 122 arctic-alpine species of Lepidoptera from the Alps, Fennoscandia and North America. It reveals intraspecific variability in the barcode region ranging from 0.00-10.08%. Eleven supposedly different species pairs or groups show close genetic similarity, suggesting possible synonymy in many cases. However, a total of 33 species show evidence of cryptic diversity as evidenced by the presence of lineages with over 2% maximum barcode divergence in Europe, in North America or between the two continents. Our study also reveals cases where taxonomic names have been used inconsistently between regions and exposes misidentifications. Overall, DNA barcodes have great potential to both increase taxonomic resolution and to make decisions concerning the taxonomic status of allopatric populations more objective.

Genetic consequences of interglacial isolation in a steppe bird

In response to climate changes that have occurred during Pleistocene glacial cycles, taxa associated to steppe vegetation might have followed a pattern of historical evolution in which isolation and fragmentation of populations occurred during the short interglacials and expansion events occurred during the long glacial periods, in contrast to the pattern described for temperate species. Here, we use molecular genetic data to evaluate this idea in a steppe bird with Palaearctic distribution, the little bustard (Tetrax tetrax). Overall, extremely low genetic diversity and differentiation was observed among eight little bustard populations distributed in Spain and France. Mismatch distribution analyses showed that most little bustard populations expanded during cooling periods previous to, and just after, the last interglacial period (127,000-111,000 years before present), when steppe habitats were widespread across Europe. Coalescent-based methods suggested that glacial expansions have resulted in substantial admixture in Western Europe due to the existence of different interglacial refugia. Our results are consistent with a model of evolution and genetic consequences of Pleistocene cycles with low between-population genetic differentiation as a result of short-term isolation periods during interglacials and long-term exchange during glacial periods.

Genetic population structure of the alpine species Rhododendron pseudochrysanthum sensu lato (Ericaceae) inferred from chloroplast and nuclear DNA

Background

A complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence. We examined the phylogeography and the evolutionary history of the Rhododendron pseudochrysanthum s. l.

Results

Systematic inconsistency was detected between gene genealogies of the cpDNA and nrDNA. Rooted at R. hyperythrum and R. formosana, both trees lacked reciprocal monophyly for all members of the complex. For R. pseudochrysanthum s.l., the spatial distribution of the cpDNA had a noteworthy pattern showing high genetic differentiation (F_ST = 0.56-0.72) between populations in the Yushan Mountain Range and populations of the other mountain ranges.

Conclusion

Both incomplete lineage sorting and interspecific hybridization/introgression may have contributed to the lack of monophyly among R. hyperythrum, R. formosana and R. pseudochrysanthum s.l. Independent colonizations, plus low capabilities of seed dispersal in current environments, may have resulted in the genetic differentiation between populations of different mountain ranges. At the population level, the populations of Central, and Sheishan Mountains may have undergone postglacial demographic expansion, while populations of the Yushan Mountain Range are likely to have remained stable ever since the colonization. In contrast, the single population of the Alishan Mountain Range with a fixed cpDNA haplotype may have experienced bottleneck/founder's events.

Climate change and evolutionary adaptations at species' range margins

During recent climate warming, many insect species have shifted their ranges to higher latitudes and altitudes. These expansions mirror those that occurred after the Last Glacial Maximum when species expanded from their ice age refugia. Postglacial range expansions have resulted in clines in genetic diversity across present-day distributions, with a reduction in genetic diversity observed in a wide range of insect taxa as one moves from the historical distribution core to the current range margin. Evolutionary increases in dispersal at expanding range boundaries are commonly observed in virtually all insects that have been studied, suggesting a positive feedback between range expansion and the evolution of traits that accelerate range expansion. The ubiquity of this phenomenon suggests that it is likely to be an important determinant of range changes. A better understanding of the extent and speed of adaptation will be crucial to the responses of biodiversity and ecosystems to climate change.