Extra-Mediterranean glacial refugia in a Mediterranean faunal element: the phylogeography of the chalk-hill blue Polyommatus coridon (Lepidoptera, Lycaenidae)

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.

From DNA barcodes to ecology: Meta-analysis of central European beetles reveal link with species ecology but also to data pattern and gaps

DNA barcoding has been used worldwide to identify biological specimens and to delimit species. It represents a cost-effective, fast, and efficient way to assess biodiversity with help of the public Barcode of Life Database (BOLD) accounting for more than 236,000 animal species and more than 10 million barcode sequences. Here, we performed a meta-analysis of available barcode data of central European Coleoptera to detect intraspecific genetic patterns among ecological groups in relation to geographic distance with the aim to investigate a possible link between infraspecific variation and species ecology. We collected information regarding feeding style, body size, as well as habitat and biotope preferences. Mantel tests and two variants of Procrustes analysis, both involving the Principal Coordinates Neighborhood Matrices (PCNM) approach, were applied on genetic and geographic distance matrices. However, significance levels were too low to further use the outcome for further trait investigation: these were in mean for all ecological guilds only 7.5, 9.4, or 15.6% for PCNM + PCA, NMDS + PCA, and Mantel test, respectively, or at best 28% for a single guild. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. Results suggest that increased numbers of species, sampling localities, and specimens for a chosen area of interest may give new insights to explore barcode data and species ecology for the scope of conservation on a larger scale. We performed a meta-analysis of available barcode data of central European beetles to detect intraspecific genetic patterns among ecological groups in relation to geographic distance, regarding feeding style, body size, as well as habitat and biotope preferences. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. However, significance levels were too low to further use the outcome for further trait investigation.

Linking large-scale genetic structure of three Argynnini butterfly species to geography and environment

Understanding which factors and processes are associated with genetic differentiation within and among species remains a major goal in evolutionary biology. To explore differences and similarities in genetic structure and its association with geographical and climatic factors in sympatric sister species, we conducted a large‐scale (>32° latitude and >36° longitude) comparative phylogeographical study on three Argynnini butterfly species (Speyeria aglaja, Fabriciana adippe and F. niobe) that have similar life histories, but differ in ecological generalism and dispersal abilities. Analyses of nuclear (ddRAD‐sequencing derived SNP markers) and mitochondrial (COI sequences) data revealed differences between species in genetic structure and how genetic differentiation was associated with climatic factors (temperature, solar radiation, precipitation, wind speed). Geographical proximity accounted for much of the variation in nuclear and mitochondrial structure and evolutionary relationships in F. adippe and F. niobe, but only explained the pattern observed in the nuclear data in S. aglaja, for which mitonuclear discordance was documented. In all species, Iberian and Balkan individuals formed genetic clusters, suggesting isolation in glacial refugia and limited postglacial expansion. Solar radiation and precipitation were associated with the genetic structure on a regional scale in all species, but the specific combinations of environmental and geographical factors linked to variation within species were unique, pointing to species‐specific responses to common environments. Our findings show that the species share similar colonization histories, and that the same ecological factors, such as niche breadth and dispersal capacity, covary with genetic differentiation within these species to some extent, thereby highlighting the importance of comparative phylogeographical studies in sympatric sister species.

Geographical disjunction and environmental conditions drive intraspecific differentiation in the chalk-hill blue butterfly

Drivers of evolution are often related to geographical isolation and/or diverging environmental conditions. Spatial variation in neutral genetic markers mostly reflects past geographical isolation, i.e. long-lasting allopatry, whereas morphology is often driven by local environmental conditions, resulting in more rapid evolution. In Europe, most thermophilic species persisted during the past glacial periods in geographically disjunct refugia, representing long-lasting isolates, frequently with diverging environmental conditions. This situation has driven the evolution of intraspecific signatures in species. Here, we analysed wing shape and wing pigmentation of the chalk-hill blue butterfly, Polyommatus coridon, across its entire distribution range restricted to the western Palaearctic. In addition, we compiled abiotic environmental parameters for each sampling site. Wing colour patterns differentiated a western and an eastern lineage. These lineages might represent two main Pleistocene refugia and differentiation centres, one located on the Italian Peninsula and the other in the Balkan region. The two lineages showed evidence of hybridization across Central Europe, from the Alps and across Germany. The intraspecific differentiation was strongest in the width of the brown band on the outer margin of the wings. The morphological structures obtained are in line with genetic signatures found in previous studies, but the latter are more fine-grained. Current environmental conditions, such as mean temperatures, were only marginally correlated with colour patterns. Our study underlines that Pleistocene range shifts, often resulting in allopatric isolation, shape intraspecific phenotypic structures within species; that pigmentation responds in a more sensitive manner to spatial disjunction than wing shape; and that morphometric and genetic structures in P. coridon provide concordant patterns and thus support identical biogeographical conclusions.

The effects of drift and selection on latitudinal genetic variation in Scandinavian common toads (Bufo bufo) following postglacial recolonisation

Clinal variation is paramount for understanding the factors shaping genetic diversity in space and time. During the last glacial maximum, northern Europe was covered by glacial ice that rendered the region uninhabitable for most taxa. Different evolutionary processes during and after the recolonisation of this area from different glacial refugia have affected the genetic landscape of the present day European flora and fauna. In this study, we focus on the common toad (Bufo bufo) in Sweden and present evidence suggesting that these processes have resulted in two separate lineages of common toad, which colonised Sweden from two directions. Using ddRAD sequencing data for demographic modelling, structure analyses, and analysis of molecular variance (AMOVA), we provide evidence of a contact zone located between Uppland and Västerbotten in central Sweden. Genetic diversity was significantly higher in southern Sweden compared to the north, in accordance with a pattern of decreased genetic diversity with increasing distance from glacial refugia. Candidate genes under putative selection are identified through outlier detection and gene-environment association methods. We provide evidence of divergent selection related to stress response and developmental processes in these candidate genes. The colonisation of Sweden by two separate lineages may have implications for how future conservation efforts should be directed by identifying management units and putative local adaptations.

Wolbachia affects mitochondrial population structure in two systems of closely related Palaearctic blue butterflies

The bacterium Wolbachia infects many insect species and spreads by diverse vertical and horizontal means. As co-inherited organisms, these bacteria often cause problems in mitochondrial phylogeny inference. The phylogenetic relationships of many closely related Palaearctic blue butterflies (Lepidoptera: Lycaenidae: Polyommatinae) are ambiguous. We considered the patterns of Wolbachia infection and mitochondrial diversity in two systems: Aricia agestis/Aricia artaxerxes and the Pseudophilotes baton species complex. We sampled butterflies across their distribution ranges and sequenced one butterfly mitochondrial gene and two Wolbachia genes. Both butterfly systems had uninfected and infected populations, and harboured several Wolbachia strains. Wolbachia was highly prevalent in A. artaxerxes and the host's mitochondrial structure was shallow, in contrast to A. agestis. Similar bacterial alleles infected both Aricia species from nearby sites, pointing to a possible horizontal transfer. Mitochondrial history of the P. baton species complex mirrored its Wolbachia infection and not the taxonomical division. Pseudophilotes baton and P. vicrama formed a hybrid zone in Europe. Wolbachia could obscure mitochondrial history, but knowledge on the infection helps us to understand the observed patterns. Testing for Wolbachia should be routine in mitochondrial DNA studies.

Ancient divergence, a crisis of salt and another of ice shaped the evolution of the west Mediterranean butterfly Euchloe tagis

The Mediterranean region is an extremely complex hotspot where, since the Miocene, extensive geological, habitat and climate changes have taken place, alternating between warm and cold periods. These phenomena have taken a toll on the genetic composition of species, and surviving lineages have often adapted locally and diverged to the point of (complete) speciation. To study these phenomena, in this study we used one of the most enigmatic butterflies, the Portuguese dappled white, Euchloe tagis, a west Mediterranean endemic with fragmented, morphologically differentiated populations whose status have long been disputed. Even its affiliations with other Anthocharidini are largely unresolved. We used mitochondrial and nuclear markers under a phylogenetic and phylogeographical framework to evaluate its placement among relatives and population differentiation, reconstructing its evolutionary history. We found that this species had a Miocene origin ~15 Mya and was nearest to Euchloe s.s. and Elphinstonia. Its populations showed high genetic diversity but all coalesced to 5.3 Mya, when European and all but one African population diverged. Our multiple findings concerning the evolution of E. tagis through a changing, narrow habitable area might provide a more general perspective on how species survive within this hotspot of paramount importance.

Rapid climate change results in long-lasting spatial homogenization of phylogenetic diversity

Scientific understanding of biodiversity dynamics, resulting from past climate oscillations and projections of future changes in biodiversity, has advanced over the past decade. Little is known about how these responses, past or future, are spatially connected. Analyzing the spatial variability in biodiversity provides insight into how climate change affects the accumulation of diversity across space. Here, we evaluate the spatial variation of phylogenetic diversity of European seed plants among neighboring sites and assess the effects of past rapid climate changes during the Quaternary on these patterns. Our work shows a marked homogenization in phylogenetic diversity across Central and Northern Europe linked to high climate change velocity and large distances to refugia. Our results suggest that the future projected loss in evolutionary heritage may be even more dramatic, as homogenization in response to rapid climate change has occurred among sites across large landscapes, leaving a legacy that has lasted for millennia.

How past climate change has affected biodiversity over large spatial scales remains underexplored. Here, the authors find marked homogenization in flowering plant phylogenetic diversity across Central and Northern Europe linked to rapid climate change and large distances to glacial refugia.

Spatio-temporal formation of the genetic diversity in the Mediterranean dwelling lichen during the Neogene and Quaternary epochs

Genetic patterns of lichenized fungi often display a mosaic-like and difficult to interpret structure blurring their evolutionary history. The genetic diversity and phylogeographic pattern of a mycobiont of the predominantly Mediterranean dwelling lichen Solenopsora candicans were investigated on the base of extensive sampling (361 individuals, 77 populations) across its entire distribution range. We tested whether the genetic pattern of S. candicans mirrors paleoclimatic and paleogeological events in the Mediterranean and adjacent regions. The divergence time estimates indicated a Tertiary origin for S. candicans, with formation of intraspecific diversity initiated in the Late Miocene. The distribution of the most divergent haplotypes, mostly of a pre-Pleistocene origin, was restricted to the eastern or western extremities of the Mediterranean exhibiting Kiermack disjunction. The population genetic diversity analyses indicated multiple diversity centres and refugia for S. candicans across the entire Mediterranean Basin. While the south Mediterranean regions harboured both the Tertiary and Quaternary born diversity, conforming to the 'cumulative refugia' paradigm, the Apennine and Balkan Peninsulas in the north hosted mostly younger Pleistocene haplotypes and lineages. The recent population expansion of S. candicans might have occurred in the middle Pleistocene with a population burst in the Apennine and Balkan peninsulas. The presence of unique haplotypes in Central Europe indicates the existence of extra-Mediterranean microrefugia. This study presents the first comprehensive lichen phylogeography from the Mediterranean region and simultaneously reports for the first time the glacial survival of a warm-adapted lichen in the temperate zone.

Polyploid evolution: The ultimate way to grasp the nettle

Polyploidy is one of the major forces of plant evolution and widespread mixed-ploidy species offer an opportunity to evaluate its significance. We therefore selected the cosmopolitan species Urtica dioica (stinging nettle), examined its cytogeography and pattern of absolute genome size, and assessed correlations with bioclimatic and ecogeographic data (latitude, longitude, elevation). We evaluated variation in ploidy level using an extensive dataset of 7012 samples from 1317 populations covering most of the species' distribution area. The widespread tetraploid cytotype (87%) was strongly prevalent over diploids (13%). A subsequent analysis of absolute genome size proved a uniform Cx-value of core U. dioica (except for U. d. subsp. cypria) whereas other closely related species, namely U. bianorii, U. kioviensis and U. simensis, differed significantly. We detected a positive correlation between relative genome size and longitude and latitude in the complete dataset of European populations and a positive correlation between relative genome size and longitude in a reduced dataset of diploid accessions (the complete dataset of diploids excluding U. d. subsp. kurdistanica). In addition, our data indicate an affinity of most diploids to natural and near-natural habitats and that the tetraploid cytotype and a small part of diploids (population from the Po river basin in northern Italy) tend to inhabit synanthropic sites. To sum up, the pattern of ploidy variation revealed by our study is in many aspects unique to the stinging nettle, being most likely first of all driven by the greater ecological plasticity and invasiveness of the tetraploid cytotype.

Red deer in Iberia: Molecular ecological studies in a southern refugium and inferences on European postglacial colonization history

The red deer (Cervus elaphus) is a widespread wild ungulate in Europe that has suffered strong anthropogenic impacts over their distribution during the last centuries, but also at the present time, due its economic importance as a game species. Here we focus on the evolutionary history of the red deer in Iberia, one of the three main southern refugial areas for temperate species in Europe, and addressed the hypothesis of a cryptic refugia at higher latitudes during the Last Glacial Maximum (LGM). A total of 911 individuals were sampled, genotyped for 34 microsatellites specifically developed for red deer and sequenced for a fragment of 670 bp of the mitochondrial (mtDNA) D-loop. The results were combined with published mtDNA sequences, and integrated with species distribution models and historical European paleo-distribution data, in order to further examine the alternative glacial refugial models and the influence of cryptic refugia on European postglacial colonization history. Clear genetic differentiation between Iberian and European contemporary populations was observed at nuclear and mtDNA levels, despite the mtDNA haplotypes central to the phylogenetic network are present across western Europe (including Iberia) suggesting a panmictic population in the past. Species distribution models, fossil records and genetic data support a timing of divergence between Iberian and European populations that overlap with the LGM. A notable population structure was also found within the Iberian Peninsula, although several populations displayed high levels of admixture as a consequence of recent red deer translocations. Five D-loop sub-lineages were found in Iberia that belong to the Western European mtDNA lineage, while there were four main clusters based on analysis of nuclear markers. Regarding glacial refugial models, our findings provide detailed support for the hypothesis that red deer may have persisted in cryptic northern refugia in western Europe during the LGM, most likely in southern France, southern Ireland, or in a region between them (continental shelf), and these regions were the source of individuals during the European re-colonization. This evidence heightens the importance of conserving the high mitochondrial and nuclear diversity currently observed in Iberian populations.

Extra-Mediterranean glacial refuges in barred and common grass snakes (Natrix helvetica, N. natrix)

Extra-Mediterranean glacial refugia of thermophilic biota, in particular in northern latitudes, are controversial. In the present study we provide genetic evidence for extra-Mediterranean refugia in two species of grass snake. The refuge of a widely distributed western European lineage of the barred grass snake (Natrix helvetica) was most likely located in southern France, outside the classical refuges in the southern European peninsulas. One genetic lineage of the common grass snake (N. natrix), distributed in Scandinavia, Central Europe and the Balkan Peninsula, had two distinct glacial refuges. We show that one was located in the southern Balkan Peninsula. However, Central Europe and Scandinavia were not colonized from there, but from a second refuge in Central Europe. This refuge was located in between the northern ice sheet and the Alpine glaciers of the last glaciation and most likely in a permafrost region. Another co-distributed genetic lineage of N. natrix, now massively hybridizing with the aforementioned lineage, survived the last glaciation in a structured refuge in the southern Balkan Peninsula, according to the idea of 'refugia-within-refugia'. It reached Central Europe only very recently. This study reports for the first time the glacial survival of a thermophilic egg-laying reptile species in Central Europe.