From Africa to Europe and back: refugia and range shifts cause high genetic differentiation in the Marbled White butterfly Melanargia galathea

From Western Asia to the Mediterranean Basin: Diversification of the Widespread Euphorbia nicaeensis Alliance (Euphorbiaceae)

The Mediterranean Basin is an important biodiversity hotspot and one of the richest areas in the world in terms of plant diversity. Its flora parallels in several aspects that of the Eurasian steppes and the adjacent Irano-Turanian floristic region. The Euphorbia nicaeensis alliance spans this immense area from the western Mediterranean to Central Asia. Using an array of complementary methods, ranging from phylogenomic and phylogenetic data through relative genome size (RGS) estimation to morphometry, we explored relationships and biogeographic connections among taxa of this group. We identified the main evolutionary lineages, which mostly correspond to described taxa. However, despite the use of highly resolving Restriction Site Associated DNA (RAD) sequencing data, relationships among the main lineages remain ambiguous. This is likely due to hybridisation, lineage sorting triggered by rapid range expansion, and polyploidisation. The phylogenomic data identified cryptic diversity in the Mediterranean, which is also correlated with RGS and, partly, also, morphological divergence, rendering the description of a new species necessary. Biogeographic analyses suggest that Western Asia is the source area for the colonisation of the Mediterranean by this plant group and highlight the important contribution of the Irano-Turanian region to the high diversity in the Mediterranean Basin. The diversification of the E. nicaeensis alliance in the Mediterranean was triggered by vicariance in isolated Pleistocene refugia, morphological adaptation to divergent ecological conditions, and, to a lesser extent, by polyploidisation.

Stability in the South, Turbulence Toward the North: Evolutionary History of Aurinia saxatilis (Brassicaceae) Revealed by Phylogenomic and Climatic Modelling Data

The Balkan Peninsula played an important role in the evolution of many Mediterranean plants and served as a major source for post-Pleistocene colonisation of central and northern Europe. Its complex geo-climatic history and environmental heterogeneity significantly influenced spatiotemporal diversification and resulted in intricate phylogeographic patterns. To explore the evolutionary dynamics and phylogeographic patterns within the widespread eastern Mediterranean and central European species Aurinia saxatilis, we used a combination of phylogenomic (restriction-site associated DNA sequencing, RADseq) and phylogenetic (sequences of the plastid marker ndhF) data as well as species distribution models generated for the present and the Last Glacial Maximum (LGM). The inferred phylogenies retrieved three main geographically distinct lineages. The southern lineage is restricted to the eastern Mediterranean, where it is distributed throughout the Aegean area, the southern Balkan Peninsula, and the southern Apennine Peninsula, and corresponds to the species main distribution area during the LGM. The eastern lineage extends from the eastern Balkan Peninsula over the Carpathians to central Europe, while the central lineage occupies the central Balkan Peninsula. Molecular dating places the divergence among all the three lineages to the early to middle Pleistocene, indicating their long-term independent evolutionary trajectories. Our data revealed an early divergence and stable in situ persistence of the southernmost, eastern Mediterranean lineage, whereas the mainland, south-east European lineages experienced more complex and turbulent evolutionary dynamics triggered by Pleistocene climatic oscillations. Our data also support the existence of multiple glacial refugia in southeast Europe and highlight the central Balkan Peninsula not only as a cradle of lineage diversifications but also as a source of lineage dispersal. Finally, the extant genetic variation within A. saxatilis is congruent with the taxonomic separation of peripatric A. saxatilis subsp. saxatilis and A. saxatilis subsp. orientalis, whereas the taxonomic status of A. saxatilis subsp. megalocarpa remains doubtful.

Climatic refugia and reduced extinction correlate with underdispersion in mammals and birds in Africa

Macroevolutionary patterns, often inferred from metrics of community relatedness, are often used to ascertain major evolutionary processes shaping communities. These patterns have been shown to be informative of biogeographic barriers, of habitat suitability and invasibility (especially with regard to environmental filtering), and of regions that function as evolutionary cradles (i.e., sources of diversification) or museums (i.e., regions of reduced extinction). Here, we analyzed continental datasets of mammal and bird distributions to identify primary drivers of community evolution on the African continent for mostly endothermic vertebrates. We find that underdispersion (i.e., relatively low phylogenetic diversity compared to species richness) closely correlates with specific ecoregions that have been identified as climatic refugia in the literature, regardless of whether these specific regions have been touted as cradles or museums. Using theoretical models of identical communities that differ only with respect to extinction rates, we find that even small suppressions of extinction rates can result in underdispersed communities, supporting the hypothesis that climatic stability can lead to underdispersion. We posit that large-scale patterns of under- and overdispersion between regions of similar species richness are more reflective of a particular region's extinction potential, and that the very nature of refugia can lead to underdispersion via the steady accumulation of species richness through diversification within the same ecoregion during climatic cycles. Thus, patterns of environmental filtering can be obfuscated by environments that coincide with biogeographic refugia, and considerations of regional biogeographic history are paramount for inferring macroevolutionary processes.

Combined climatic and anthropogenic stress threaten resilience of important wetland sites in an arid region

Climate change and anthropogenic perturbation threaten resilience of wetlands globally, particularly in regions where environmental conditions are already hot and dry, and human impacts are rapidly intensifying and expanding. Here we assess the vulnerability of Ramsar wetlands of six North African countries (Western Sahara, Morocco, Algeria, Tunisia, Libya, and Egypt) by asking three questions: (1) what are the recent anthropogenic changes that the wetlands experienced? (2) what are the projected future climatic changes? (3) how wetlands with different conservation priorities and globally threatened species are impacted by anthropogenic pressures? We used climatic data (historical and future projections) from WorldClim 2, drought index (SPEI), and human footprint index (HFI for 2000 and 2019) to estimate anthropogenic pressures, as well as waterbird conservation value (WCV: a metric indicating conservation priority of sites) and the breeding distribution of three threatened waterbird species (Aythya nyroca, Marmaronetta angustirostris, and Oxyura leucocephala) to understand how biodiversity is impacted by anthropogenic pressure. We found that temperature, precipitation, drought, and human footprint index (HFI) increased during earlier decades. Interestingly, areas with high HFI are projected to encounter lower warming but more severe drought. We also found that WCV was positively correlated with the magnitude of current HFI, indicating that sites of high conservation value for waterbirds encounter higher levels of anthropogenic pressure. The breeding range of the three threatened species of waterbirds showed a marked increase in HFI and is projected to experience a severe increase in temperature by 2081-2100, especially under the high emission scenario (SSP8.5) where environmental temperature becomes closer to the species critical maximum. Our results highlight the importance of integrating new conservation measures that increase the resilience of North African protected wetlands to reduce extinction risk to biodiversity.

Identification of winter moth (Operophtera brumata) refugia in North Africa and the Italian Peninsula during the last glacial maximum

Numerous studies have shown that the genetic diversity of species inhabiting temperate regions has been shaped by changes in their distributions during the Quaternary climatic oscillations. For some species, the genetic distinctness of isolated populations is maintained during secondary contact, while for others, admixture is frequently observed. For the winter moth (Operophtera brumata), an important defoliator of oak forests across Europe and northern Africa, we previously determined that contemporary populations correspond to genetic diversity obtained during the last glacial maximum (LGM) through the use of refugia in the Iberian and Aegean peninsulas, and to a lesser extent the Caucasus region. Missing from this sampling were populations from the Italian peninsula and from North Africa, both regions known to have played important roles as glacial refugia for other species. Therefore, we genotyped field-collected winter moth individuals from southern Italy and northwestern Tunisia-the latter a region where severe oak forest defoliation by winter moth has recently been reported-using polymorphic microsatellite. We reconstructed the genetic relationships of these populations in comparison to moths previously sampled from the Iberian and Aegean peninsulas, the Caucasus region, and western Europe using genetic distance, Bayesian clustering, and approximate Bayesian computation (ABC) methods. Our results indicate that both the southern Italian and the Tunisian populations are genetically distinct from other sampled populations, and likely originated in their respective refugium during the LGM after diverging from a population that eventually settled in the Iberian refugium. These suggest that winter moth populations persisted in at least five Mediterranean LGM refugia. Finally, we comment that outbreaks by winter moth in northwestern Tunisia are not the result of a recent introduction of a nonnative species, but rather are most likely due to land use or environmental changes.

A continent-scale study of the social structure and phylogeography of the bent-wing bat, Miniopterus schreibersii (Mammalia: Chiroptera), using new microsatellite data

Miniopterus schreibersii is a cave-dwelling bat species with a wide distribution in the western Palearctic spanning southern and central Europe, North Africa, Anatolia, the Middle East, and the Caucasus. We investigated the social structure and its effects on the genetic makeup of this species, using 10 nuclear microsatellite markers and a partial fragment of the mitochondrial cytochrome b gene. Samples were examined from the species' entire circum-Mediterranean range. Local structuring that was previously detected among populations of M. schreibersii using mitochondrial markers was not observed for microsatellite markers, indicating male-biased dispersal for the species. Some support was found for postglacial expansions in Europe, with Anatolia potentially acting as the primary refugium during the Last Glacial Maximum (LGM). However, support for this hypothesis is not as strong as that previously detected using mitochondrial DNA markers. This is likely due to the diminishing effect of male-mediated dispersal, replenishing the nuclear diversity faster than the mitochondrial diversity in regions that are relatively far from the glacial refugia.

Biogeographical patterns in the structural blue of male Polyommatus icarus butterflies

Color is a widely used communication channel in the living world for a variety of functions ranging from sexual communication to warning colors. A particularly rich spectrum of colors appears on the wings of many butterflies. The males of lycaenid butterflies often exhibit a conspicuous blue coloration generated by photonic nanoarchitectures on their dorsal wing surfaces. Using UV-VIS spectroscopy, we investigated the spatio-temporal variations of this coloration for Polyommatus icarus butterflies, considering an interval of more than 100 years and a geographical range spanning Europe (west) and Asia (east). The blue coloration in Hungary is very stable both within a year (three broods typical in Hungary) and within the period of 100 years (more than 300 generations). East-west geographic variation was investigated among 314 male P. icarus butterflies. In agreement with earlier genetic and morphometric studies, it was found that the western males are not divided in distinct lineages. Clear differences in coloration were found between the eastern and western groups, with a transition in the region of Turkey. These differences are tentatively attributed to bottleneck effects during past glaciations.

Comparative phylogeography of amphibians and reptiles in Algeria suggests common causes for the east-west phylogeographic breaks in the Maghreb

A series of phylogeographic studies in the Maghreb identified a repeated pattern of deep genetic divergence between an eastern (Tunisia) and western (Morocco) lineage for several taxa but lack of sampling in Algeria made it difficult to know if the range limits between the eastern and western lineages were shared among taxa or not. To address this question, we designed a comparative phylogeographic study using 8 reptile and 3 amphibian species with wide distribution in the Maghreb as models. We selected species where previous studies had identified an East-West phylogeographic divide and collected sampled in Algeria to 1) examine whether the simple East-West divergence pattern still holds after filling the sampling gap in Algeria or if more complex diversity patterns emerge; 2) if the E-W pattern still holds, test whether the limits between the E and W clades are shared between species, suggesting that common historical process caused the E-W divergences; 3) if E-W limits are shared between species, use information on the age of the divergence to identify possible geological or climatic events that could have triggered these E-W differentiations. We found that the E-W pattern was generally maintained after additional sampling in Algeria and identified two common disjunction areas, one around the Algeria-Morocco border, the other one in Kabylia (central Algeria), suggesting that common historical mechanisms caused the E-W divergences in the Maghreb. Our estimates for the times to most common recent ancestors to the E and W clades span a wide range between the Messinian salinity crisis and the Plio-Pleistocene limit (except for one older split), suggesting different origins for the initial divergences and subsequent preservation of the E and W lineages in common climatic refugia in the west and the east of the Maghreb.

Varying and unchanging whiteness on the wings of dusk-active and shade-inhabiting Carystoides escalantei butterflies

Whiteness, although frequently apparent on the wings, legs, antennae, or bodies of many species of moths and butterflies, along with other colors and shades, has often escaped our attention. Here, we investigate the nanostructure and microstructure of white spots on the wings of Carystoides escalantei, a dusk-active and shade-inhabiting Costa Rican rain forest butterfly (Hesperiidae). On both males and females, two types of whiteness occur: angle dependent (dull or bright) and angle independent, which differ in the microstructure, orientation, and associated properties of their scales. Some spots on the male wings are absent from the female wings. Whether the angle-dependent whiteness is bright or dull depends on the observation directions. The angle-dependent scales also show enhanced retro-reflection. We speculate that the biological functions and evolution of Carystoides spot patterns, scale structures, and their varying whiteness are adaptations to butterfly's low light habitat and to airflow experienced on the wing base vs. wing tip.

Postglacial recolonization shaped the genetic diversity of the winter moth (Operophtera brumata) in Europe

Changes in climate conditions, particularly during the Quaternary climatic oscillations, have long been recognized to be important for shaping patterns of species diversity. For species residing in the western Palearctic, two commonly observed genetic patterns resulting from these cycles are as follows: (1) that the numbers and distributions of genetic lineages correspond with the use of geographically distinct glacial refugia and (2) that southern populations are generally more diverse than northern populations (the “southern richness, northern purity” paradigm). To determine whether these patterns hold true for the widespread pest species the winter moth (Operophtera brumata), we genotyped 699 individual winter moths collected from 15 Eurasian countries with 24 polymorphic microsatellite loci. We find strong evidence for the presence of two major genetic clusters that diverged ~18 to ~22 ka, with evidence that secondary contact (i.e., hybridization) resumed ~ 5 ka along a well‐established hybrid zone in Central Europe. This pattern supports the hypothesis that contemporary populations descend from populations that resided in distinct glacial refugia. However, unlike many previous studies of postglacial recolonization, we found no evidence for the “southern richness, northern purity” paradigm. We also find evidence for ongoing gene flow between populations in adjacent Eurasian countries, suggesting that long‐distance dispersal plays an important part in shaping winter moth genetic diversity. In addition, we find that this gene flow is predominantly in a west‐to‐east direction, suggesting that recently debated reports of cyclical outbreaks of winter moth spreading from east to west across Europe are not the result of dispersal.

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.

Why do cryptic species tend not to co-occur? A case study on two cryptic pairs of butterflies

As cryptic diversity is being discovered, mostly thanks to advances in molecular techniques, it is becoming evident that many of these taxa display parapatric distributions in mainland and that they rarely coexist on islands. Genetic landscapes, haplotype networks and ecological niche modeling analyses were performed for two pairs of non-sister cryptic butterfly species, Aricia agestis-A. cramera and Polyommatus icarus-P. celina (Lycaenidae), to specifically assess non-coexistence on western Mediterranean islands, and to test potential causes producing such chequered distribution patterns. We show that the morphologically and ecologically equivalent pairs of species do not coexist on any of the studied islands, although nearly all islands are colonized by one of them. According to our models, the cryptic pairs displayed marked climatic preferences and 'precipitation during the driest quarter' was recovered as the most important climatic determinant. However, neither dispersal capacity, nor climatic or ecological factors fully explain the observed distributions across particular sea straits, and the existence of species interactions resulting in mutual exclusion is suggested as a necessary hypothesis. Given that the studied species are habitat generalists, feeding on virtually unlimited resources, we propose that reproductive interference, together with climatic preferences, sustain density-dependent mechanisms like "founder takes all" and impede coexistence on islands. Chequered distributions among cryptic taxa, both sister and non-sister, are common in butterflies, suggesting that the phenomenon revealed here could be important in determining biodiversity patterns.

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.

Spatial genetic structure patterns of phenotype-limited and boundary-limited expanding populations: a simulation study

Range expansions may create a unique spatial genetic pattern characterized by alternate genetically homogeneous domains and allele frequency clines. Previous attempts to model range expansions have mainly focused on the loss of genetic diversity during expansions. Using individual-based models, we examined spatial genetic patterns under two expansion scenarios, boundary-limited range expansions (BLRE) and phenotype-limited range expansions (PhLRE). Our simulation revealed that the genetic diversity within populations lost quickly during the range expansion, while the genetic difference accumulated between populations. Consequently, accompanying the expansions, the overall diversity featured a slow decrease. Specifically, during BLREs, high speed of boundary motion facilitated the maintenance of total genetic diversity and sharpened genetic clines. Very slight constraints on boundary motion of BLREs drastically narrowed the homogeneous domains and increased the allele frequency fluctuations from those levels exhibited by PhLREs. Even stronger constraints, however, surprisingly brought the width of homogeneous domains and the allele frequency fluctuations back to the normal levels of PhLREs. Furthermore, high migration rates maintained a higher total genetic diversity than low ones did during PhLREs. Whereas, the total genetic diversities during BLREs showed a contrary pattern: higher when migration was low than those when migration was high. Besides, the increase of migration rates helped maintain a greater number of homogeneous domains during PhLREs, but their effects on the number of homogeneous domains during BLREs were not monotonous. Previous studies have showed that the homogenous domains can merge to form a few broad domains as the expansion went on, leading to fewer homogeneous domains. Our simulations, meanwhile, revealed that the range expansions could also rebuild homogeneous domains from the clines during the range expansion. It is possible that that the number of homogeneous domains was determined by the interaction of merging and newly emerging homogeneous domains.

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.

Evolution and radiation in the scorpion Buthus elmoutaouakili Lourenco and Qi 2006 (Scorpiones: Buthidae) at the foothills of the Atlas Mountains (North Africa)

When low dispersal ability of an organism meets geographical barriers, the evolution of inter- and intraspecific differentiation is often facilitated. In the Atlas massif of North Africa, the genus Buthus splits into several species and diverges into numerous genetic lineages, often following the orographic structures of mountain systems. Such high mountain ranges often act as barriers for species with reduced mobility even on small spatial scales. To study the effect of orographic structures on organisms with low dispersal ability, we collected 61 individuals of the scorpion species Buthus elmoutaouakili at 18 locations around the southwestern foothills of the High Atlas and Antiatlas and in the Sousse valley (western Morocco). We analyzed intraspecific differentiation patterns within this geographically restricted area of about 100 × 50 km using 452 bp of the cytochrome oxidase I mitochondrial gene. We detected 5 distinct genetic lineages. In a second analysis, we added 61 previously published sequences from Buthus species from Europe and North Africa. Using a molecular clock approach, we detected old splits (4-5 Ma) separating the samples from 1) the western High Atlas and north of these mountains, 2) the Sousse valley and adjoining mountain areas, and 3) the southwestern Antiatlas. Further differentiation happened in the first 2 geographical groups about 3 Ma. Thus, the divergence time estimates based on a Bayesian approach support the onset of differentiation into these main clades along the Pliocene (5-2.3 Ma) when climatic oscillations started and a constant global cooling preceded the glacial-interglacial cycles of the Pleistocene. Further genetic splits into parapatric groups are detectable for the Sousse valley main group in the early Pleistocene. The climatic oscillations of the Pliocene and early Pleistocene might have caused repeated range shifts, expansions, and retractions leading to repeated vicariance, hereby producing the hierarchical structure of genetic differentiation in B. elmoutaouakili. A taxonomic revision, including morphological and molecular data, is needed to assess the status of each of these Buthus scorpion lineages.