Species radiation in the Alps: multiple range shifts caused diversification in Ringlet butterflies in the European high mountains

Museomics help resolving the phylogeny of snowfinches (Aves, Passeridae, Montifringilla and allies)

Historical specimens from museum collections provide a valuable source of material also from remote areas or regions of conflict that are not easily accessible to scientists today. With this study, we are providing a taxon-complete phylogeny of snowfinches using historical DNA from whole skins of an endemic species from Afghanistan, the Afghan snowfinch, Pyrgilauda theresae. To resolve the strong conflict between previous phylogenetic hypotheses, we generated novel mitogenome sequences for selected taxa and genome-wide SNP data using ddRAD sequencing for all extant snowfinch species endemic to the Qinghai-Tibet Plateau (QTP) and for an extended intraspecific sampling of the sole Central and Western Palearctic snowfinch species (Montifringilla nivalis). Our phylogenetic reconstructions unanimously refuted the previously suggested paraphyly of genus Pyrgilauda. Misplacement of one species-level taxon (Onychostruthus tazcanowskii) in previous snowfinch phylogenies was undoubtedly inferred from chimeric mitogenomes that included heterospecific sequence information. Furthermore, comparison of novel and previously generated sequence data showed that the presumed sister-group relationship between M. nivalis and the QTP endemic M. henrici was suggested based on flawed taxonomy. Our phylogenetic reconstructions based on genome-wide SNP data and on mitogenomes were largely congruent and supported reciprocal monophyly of genera Montifringilla and Pyrgilauda with monotypic Onychostruthus being sister to the latter. The Afghan endemic P. theresae likely originated from a rather ancient Pliocene out-of-Tibet dispersal probably from a common ancestor with P. ruficollis. Our extended trans-Palearctic sampling for the white-winged snowfinch, M. nivalis, confirmed strong lineage divergence between an Asian and a European clade dated to 1.5 - 2.7 million years ago (mya). Genome-wide SNP data suggested subtle divergence among European samples from the Alps and from the Cantabrian mountains.

A macroevolutionary role for chromosomal fusion and fission in Erebia butterflies

The impact of large-scale chromosomal rearrangements, such as fusions and fissions, on speciation is a long-standing conundrum. We assessed whether bursts of change in chromosome numbers resulting from chromosomal fusion or fission are related to increased speciation rates in Erebia, one of the most species-rich and karyotypically variable butterfly groups. We established a genome-based phylogeny and used state-dependent birth-death models to infer trajectories of karyotype evolution. We demonstrated that rates of anagenetic chromosomal changes (i.e., along phylogenetic branches) exceed cladogenetic changes (i.e., at speciation events), but, when cladogenetic changes occur, they are mostly associated with chromosomal fissions rather than fusions. We found that the relative importance of fusion and fission differs among Erebia clades of different ages and that especially in younger, more karyotypically diverse clades, speciation is more frequently associated with cladogenetic chromosomal changes. Overall, our results imply that chromosomal fusions and fissions have contrasting macroevolutionary roles and that large-scale chromosomal rearrangements are associated with bursts of species diversification.

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.

Climate, immigration and speciation shape terrestrial and aquatic biodiversity in the European Alps

Quaternary climate fluctuations can affect speciation in regional biodiversity assembly in two non-mutually exclusive ways: a glacial species pump, where isolation in glacial refugia accelerates allopatric speciation, and adaptive radiation in underused adaptive zones during ice-free periods. We detected biogeographic and genetic signatures associated with both mechanisms in the assembly of the biota of the European Alps. Age distributions of endemic and widespread species within aquatic and terrestrial taxa (amphipods, fishes, amphibians, butterflies and flowering plants) revealed that endemic fish evolved only in lakes, are highly sympatric, and mainly of Holocene age, consistent with adaptive radiation. Endemic amphipods are ancient, suggesting preglacial radiation with limited range expansion and local Pleistocene survival, perhaps facilitated by a groundwater-dwelling lifestyle. Terrestrial endemics are mostly of Pleistocene age and are thus more consistent with the glacial species pump. The lack of evidence for Holocene adaptive radiation in the terrestrial biome is consistent with faster recolonization through range expansion of these taxa after glacial retreats. More stable and less seasonal ecological conditions in lakes during the Holocene may also have contributed to Holocene speciation in lakes. The high proportion of young, endemic species makes the Alpine biota vulnerable to climate change, but the mechanisms and consequences of species loss will likely differ between biomes because of their distinct evolutionary histories.

Pronounced mito-nuclear discordance and various Wolbachia infections in the water ringlet Erebia pronoe have resulted in a complex phylogeographic structure

Several morphological and mitochondrial lineages of the alpine ringlet butterfly species Erebia pronoe have been described, indicating a complex phylogenetic structure. However, the existing data were insufficient and allow neither a reconstruction of the biogeographic history, nor an assessment of the genetic lineages. Therefore, we analysed mitochondrial (COI, NDI) and nuclear (EF1α, RPS5) gene sequences and compared them with sequences from the sister species Erebia melas. Additionally, we combined this information with morphometric data of the male genitalia and the infection patterns with Wolbachia strains, based on a WSP analysis. We obtained a distinct phylogeographic structure within the E. pronoe-melas complex with eight well-distinguishable geographic groups, but also a remarkable mito-nuclear discordance. The mito-nuclear discordance in E. melas and E. pronoe glottis can be explained by different ages of Wolbachia infections with different Wolbachia strains, associated selective sweeps, and hybridisation inhibition. Additionally, we found indications for incipient speciation of E. pronoe glottis in the Pyrenees and a pronounced range dynamic within and among the other high mountain systems of Europe. Our results emphasize the importance of combined approaches in reconstructing biogeographic patterns and evaluating phylogeographic splits.

Genomics Reveal Admixture and Unexpected Patterns of Diversity in a Parapatric Pair of Butterflies

We studied the evolutionary relationship of two widely distributed parapatric butterfly species, Melitaea athalia and Melitaea celadussa, using the ddRAD sequencing approach, as well as genital morphology and mtDNA data. M. athalia was retrieved as paraphyletic with respect to M. celadussa. Several cases of mito-nuclear discordance and morpho-genetic mismatch were found in the contact zone. A strongly diverged and marginally sympatric clade of M. athalia from the Balkans was revealed. An in-depth analysis of genomic structure detected high levels of admixture between M. athalia and M. celadussa at the contact zone, though not reaching the Balkan clade. The demographic modelling of populations supported the intermediate genetic make-up of European M. athalia populations with regards to M. celadussa and the Balkan clade. However, the dissimilarity matrix of genotype data (PCoA) suggested the Balkan lineage having a genetic component that is unrelated to the athalia-celadussa group. Although narrowly sympatric, almost no signs of gene flow were found between the main M. athalia group and the Balkan clade. We propose two possible scenarios on the historical evolution of our model taxa and the role of the last glacial maximum in shaping their current distribution. Finally, we discuss the complexities regarding the taxonomic delimitation of parapatric taxa.

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.

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.

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.

Overlooked cryptic diversity in Muschampia (Lepidoptera: Hesperiidae) adds two species to the European butterfly fauna

Cryptic species represent a challenge for documenting global biodiversity. Even in well-studied groups, such as European butterflies, the application of integrative approaches has allowed the recognition of an unexpected number of cryptic taxa. Here, we combine the analysis of mitochondrial (cytochrome c oxidase I, COI) and nuclear (internal transcribed spacer 2, ITS2) markers with geometric morphometrics of the male genitalia to study diversity within the butterfly Muschampia proto. The nuclear marker reveals three well-supported and deeply diverged lineages, which are also detected based on mitochondrial DNA, although the latter recovers one of them as paraphyletic with poor support. These lineages also present distinct male genital characters, which allow blind assignment of > 97% of specimens when applying a jackknife procedure. We conclude that M. proto comprises three cryptic species that started to differentiate ~2 Mya: M. proto, distributed in northern Africa, the Iberian Peninsula and southern France; Muschampia alta comb. & stat. nov., occurring in southern Italy and the Balkan Peninsula; and Muschampia proteides, present in the easternmost part of Europe, the Near East and Iran. This discovery adds two new species to the European butterfly fauna and highlights the necessity to continue investigating potential cryptic diversity.

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.

Secondary contact zones of closely-related Erebia butterflies overlap with narrow phenotypic and parasitic clines

Zones of secondary contact between closely related taxa are a common legacy of the Quaternary ice ages. Despite their abundance, the factors that keep species apart and prevent hybridization are often unknown. Here, we study a very narrow contact zone between three closely related butterfly species of the Erebia tyndarus species complex. Using genomic data, we first determined whether gene flow occurs and then assessed whether it might be hampered by differences in chromosome number between some species. We found interspecific gene flow between sibling species that differ in karyotype by one chromosome. Conversely, only F1 hybrids occurred between two species that have the same karyotype, forming a steep genomic cline. In a second step, we fitted clines to phenotypic, ecological and parasitic data to identify the factors associated with the genetic cline. We found clines for phenotypic data and the prevalence of the endosymbiont parasite Wolbachia to overlap with the genetic cline, suggesting that they might be drivers for separating the two species. Overall, our results highlight that some gene flow is possible between closely related species despite different chromosome numbers, but that other barriers restrict such gene flow.

Recent diversification of Chrysoritis butterflies in the South African Cape (Lepidoptera: Lycaenidae)

Although best known for its extraordinary radiations of endemic plant species, the South African fynbos is home to a great diversity of phytophagous insects, including butterflies in the genus Chrysoritis (Lepidoptera: Lycaenidae). These butterflies are remarkably uniform morphologically; nevertheless, they comprise 43 currently accepted species and 68 currently valid taxonomic names. While many species have highly restricted, dot-like distributions, others are widespread. Here, we investigate the phylogenetic and biogeographic history underlying their diversification by analyzing molecular markers from 406 representatives of all described species throughout their respective ranges. We recover monophyletic clades for both C. chrysaor and C. thysbe species-groups, and identify a set of lineages that fall between them. The estimated age of divergence for the genus is 32 Mya, and we document significantly rapid diversification of the thysbe species-group in the Pleistocene (~2 Mya). Using ancestral geographic range reconstruction, we show that West Fynbos is the most likely region of origin for the radiation of the thysbe species-group. The colonization of this region occurred 9 Mya and appears to have been followed by a long period of relative stasis before a recent increase in diversification. Thus, the thysbe radiation does not appear to have resulted from the colonization of new biogeographic areas. Rather, the impact of species interactions (with ants and plants), the appearance of key innovations, and/or the opening of new ecological niche space in the region might explain the sudden burst of speciation that occurred in this group 2 Mya. The biogeographic model suggests two different diversification processes with few historical cross-colonisations, one in eastern South Africa for the C. chrysaor group and the other in western South Africa for the remaining taxa. Distributional range assessments and ecological niche models for each species show important niche overlap, and in a few cases, complete overlap. However, these shared traits are not explained by phylogenetic history. Chrysoritis taxa frequently fly in sympatry and gene tree reticulation appears to be widespread at the species level, suggesting that several episodes of range shifts might have led to secondary sympatries, allowing limited gene flow that challenges species delimitation efforts. In addition, the unusually high diversification rate for the thysbe clade of 1.35 [0.91-1.81] lineages per million years also suggests the possibility of taxonomic oversplitting. The phylogeny presented here provides a framework for a taxonomic revision of the genus. We highlight cases of potential synonymy both in allopatry and sympatry, and stress the importance of dedicated studies to assess potential pre- and post-zygotic barriers giving rise to species delimitations of the thysbe group.

Hide-and-seek with hoverflies: Merodon aureus – a species, a complex or a subgroup?

In order to disentangle the currently confused interpretations and nomenclature of Merodon aureus and M. aeneus, we have reviewed all existing type material and species names known to us as assigned synonyms of these taxa. We resolve M. aeneus as being a junior synonym of M. aureus. We designate a lectotype for M. aureus and a neotype for M. aeneus. Additionally, we provide evidence that M. aureus, together with two newly discovered taxa (M. calidus sp. nov. and M. ortus sp. nov.), represent a complex of cryptic species named the M. aureus species complex. This complex, together with the M. unicolor species complex and the species M. pumilus, is part of the M. aureus subgroup. The M. unicolor species complex comprises two cryptic species: M. unicolor and M. albidus sp. nov. The new species are described by applying an integrative taxonomic approach using several data types (COI and 28S rRNA genes, geometric morphometry of the wings, ecological and distributional data). Based on the COI gene sequence analysis and distributional data, the pupa previously described as an immature stage of the species M. aureus is redefined as an immature stage of the new species M. calidus. Speciation within the M. aureus subgroup is discussed in the context of the phylogeographic history in the studied region.

An updated checklist of the European Butterflies (Lepidoptera, Papilionoidea)

This paper presents an updated checklist of the butterflies of Europe, together with their original name combinations, and their occurrence status in each European country. According to this checklist, 496 species of the superfamily Papilionoidea occur in Europe. Changes in comparison with the last version (2.6.2) of Fauna Europaea are discussed. Compared to that version, 16 species are new additions, either due to cryptic species most of which have been discovered by molecular methods (13 cases) or due to discoveries of Asian species on the eastern border of the European territory in the Ural mountains (three cases). On the other hand, nine species had to be removed from the list, because they either do not occur in Europe or lost their species status due to new evidence. In addition, three species names had to be changed and 30 species changed their combination due to new evidence on phylogenetic relationships. Furthermore, minor corrections were applied to some authors’ names and years of publication. Finally, the namePolyommatusottomanusLefèbvre, 1831, which is threatened by its senior synonymLycaenalegeriFreyer, 1830, is declared anomen protectum, thereby conserving its name in the current combinationLycaenaottomana.

M. van Swaay, Vila R, Vliegenthart A, Wahlberg N, Verovnik R. An updated checklist of the European Butterflies (Lepidoptera, Papilionoidea)

This paper presents an updated checklist of the butterflies of Europe, together with their original name combinations, and their occurrence status in each European country. According to this checklist, 496 species of the superfamily Papilionoidea occur in Europe. Changes in comparison with the last version (2.6.2) of Fauna Europaea are discussed. Compared to that version, 16 species are new additions, either due to cryptic species most of which have been discovered by molecular methods (13 cases) or due to discoveries of Asian species on the eastern border of the European territory in the Ural mountains (three cases). On the other hand, nine species had to be removed from the list, because they either do not occur in Europe or lost their species status due to new evidence. In addition, three species names had to be changed and 30 species changed their combination due to new evidence on phylogenetic relationships. Furthermore, minor corrections were applied to some authors' names and years of publication. Finally, the name Polyommatusottomanus Lefèbvre, 1831, which is threatened by its senior synonym Lycaenalegeri Freyer, 1830, is declared a nomen protectum, thereby conserving its name in the current combination Lycaenaottomana.

Scientific Opinion addressing the state of the science on risk assessment of plant protection products for in-soil organisms

Following a request from EFSA, the Panel on Plant Protection Products and their Residues developed an opinion on the science behind the risk assessment of plant protection products for in-soil organisms. The current risk assessment scheme is reviewed, taking into account new regulatory frameworks and scientific developments. Proposals are made for specific protection goals for in-soil organisms being key drivers for relevant ecosystem services in agricultural landscapes such as nutrient cycling, soil structure, pest control and biodiversity. Considering the time-scales and biological processes related to the dispersal of the majority of in-soil organisms compared to terrestrial non-target arthropods living above soil, the Panel proposes that in-soil environmental risk assessments are made at in- and off-field scale considering field boundary levels. A new testing strategy which takes into account the relevant exposure routes for in-soil organisms and the potential direct and indirect effects is proposed. In order to address species recovery and long-term impacts of PPPs, the use of population models is also proposed.