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.

Hybridization fuelled diversification in Spialia butterflies.

The importance of hybridization and introgression is well documented in the evolution of plants but, in insects, their role is not fully understood. Given the fact that insects are the most diverse group of organisms, assessing the impact of reticulation events on their evolution may be key to comprehend the emergence of such remarkable diversity. Here, we used an insect model, the Spialia butterflies, to gather genomic evidence of hybridization as a promoter of novel diversity. By using double-digest RADseq (ddRADseq), we explored the phylogenetic relationships between Spialia orbifer, S. rosae and S. sertorius, and documented two independent events of interspecific gene flow. Our data support that the Iberian endemism S. rosae probably received genetic material from S. orbifer in both mitochondrial and nuclear DNA, which could have contributed to a shift in the ecological preferences of S. rosae. We also show that admixture between S. sertorius and S. orbifer probably occurred in Italy. As a result, the admixed Sicilian populations of S. orbifer are differentiated from the rest of populations both genetically and morphologically, and display signatures of reproductive character displacement in the male genitalia. Additionally, our analyses indicated that genetic material from S. orbifer is present in S. sertorius along the Italian Peninsula. Our findings add to the view that hybridization is a pervasive phenomenon in nature and in butterflies in particular, with important consequences for evolution due to the emergence of novel phenotypes.

Two ways to be endemic. Alps and Apennines are different functional refugia during climatic cycles.

Endemics co-occur because they evolved in situ and persist regionally or because they evolved ex situ and later dispersed to shared habitats, generating evolutionary or ecological endemicity centres, respectively. We investigate whether different endemicity centres can intertwine in the region ranging from Alps to Sicily, by studying their butterfly fauna. We gathered an extensive occurrence data set for butterflies of the study area (27,123 records, 269 species, in cells of 0.5 × 0.5 degrees of latitude-longitude). We applied molecular-based delimitation methods (GMYC model) to 26,557 cytochrome c oxidase subunit 1 (COI) sequences of Western Palearctic butterflies. We identified entities based on molecular delimitations and/or the checklist of European butterflies and objectively attributed occurrences to their most probable entity. We obtained a zoogeographic regionalisation based on the 69 endemics of the area. Using phylogenetic ANOVA we tested if endemics from different centres differ from each other and from nonendemics for key ecological traits and divergence time. Endemicity showed high incidence in the Alps and Southern Italy. The regionalisation separated the Alps from the Italian Peninsula and Sicily. The endemics of different centres showed a high turnover and differed in phylogenetic distances, phenology and distribution traits. Endemics are on average younger than nonendemics and the Peninsula-Sicily endemics also have lower variance in divergence than those from the Alps. The observed variation identifies Alpine endemics as paleoendemics, now occupying an ecological centre, and the Peninsula-Sicily ones as neoendemics, that diverged in the region since the Pleistocene. The results challenge the common view of the Alpine-Apennine area as a single "Italian refugium".

Evidence for multiple colonisations and Wolbachia infections shaping the genetic structure of the widespread butterfly Polyommatus icarus in the British Isles.

The paradigm of isolation in southern refugia during glacial periods followed by expansions during interglacials, producing limited genetic differentiation in northern areas, dominates European phylogeography. However, the existence of complex structured populations in formerly glaciated areas, and islands connected to mainland areas during glacial maxima, call for alternative explanations. We reconstructed the mtDNA phylogeography of the widespread Polyommatus Icarus butterfly with an emphasis on the formerly glaciated and connected British Isles. We found distinct geographical structuring of CO1 haplogroups, with an ancient lineage restricted to the marginal European areas, including Northern Scotland and Outer Hebrides. Population genomic analyses, using ddRADSeq genomic markers, also reveal substantial genetic structuring within Britain. However, there is negligble mito-nuclear concordance consistent with independent demographic histories of mitochondrial versus nuclear DNA. While mtDNA-Wolbachia associations in northern Britain could account for the geographic structuring of mtDNA across most of the British Isles, for nuclear DNA markers (derived from ddRADseq data) butterflies from France cluster between northern and southern British populations - an observation consistent with a scenario of multiple recolonisation. Taken together our results suggest that contemporary mtDNA structuring in the British Isles (and potentially elsewhere in Europe) largely results from Wolbachia infections, however, nuclear genomic structuring suggests a history of at least two distinct colonisations. This two-stage colonisation scenario has previously been put forth to explain genetic diversity and structuring in other British flora and fauna. Additionally, we also present preliminary evidence for potential Wolbachia-induced feminization in the Outer Hebrides.

The Pleistocene species pump past its prime: Evidence from European butterfly sister species.

The Pleistocene glacial cycles had a profound impact on the ranges and genetic make-up of organisms. While it is clear that the contact zones that have been described for many sister taxa are secondary and have formed in the current interglacial, it is unclear when the taxa involved began to diverge. Previous estimates based on small numbers of loci are unreliable given the stochasticity of genetic drift and the contrasting effects of incomplete lineage sorting and gene flow on gene divergence. Here, we use genome-wide transcriptome data to estimate divergence for 18 sister species pairs of European butterflies showing either sympatric or contact zone distributions. We find that in most cases, species divergence predates the mid-Pleistocene transition or even the entire Pleistocene period. We also show that although post-divergence gene flow is restricted to contact zone pairs, they are not systematically younger than sympatric pairs. This suggests that contact zones are not limited to the initial stages of the speciation process, but can involve notably old taxa. Finally, we show that mitochondrial divergence and nuclear divergence are only weakly correlated and mitochondrial divergence is higher for contact zone pairs.

How long is 3 km for a butterfly? Ecological constraints and functional traits explain high mitochondrial genetic diversity between Sicily and the Italian Peninsula.

Populations inhabiting Mediterranean islands often show contrasting genetic lineages, even on islands that were connected to the mainland during glacial maxima. This pattern is generated by forces acting in historical and contemporary times. Understanding these phenomena requires comparative studies integrating genetic structure, functional traits and dispersal constraints. Using as a model the butterfly species living across the Messina strait (3 km wide) separating Sicily from the Italian Peninsula, we aimed to unravel the mechanisms limiting the dispersal of matrilines and generating genetic differentiation across a narrow sea strait. We analysed the mitochondrial COI gene of 84 butterfly species out of 90 documented in Sicily and compared them with populations from the neighbouring southern Italian Peninsula (1,398 sequences) and from the entire Palaearctic region (8,093 sequences). For each species, we regressed 13 functional traits and 2 ecological constraints to dispersal (winds experienced at the strait and climatic suitability) against genetic differentiation between Sicily and Italian Peninsula to understand the factors limiting dispersal. More than a third of the species showed different haplogroups across the strait and most of them also represented endemic haplogroups for this island. One fifth of Sicilian populations (and 32.3% of endemic lineages) had their closest relatives in distant areas, instead of the neighbouring Italian Peninsula, which suggests high relictuality. Haplotype diversity was significantly explained by the length of the flight period, an intrinsic phenology trait, while genetic differentiation was explained by both intrinsic traits (wingspan and degree of generalism) and contemporary local constraints (winds experienced at the strait and climatic suitability). A relatively narrow sea strait can produce considerable differentiation among butterfly matrilines and this phenomenon showed a largely deterministic fingerprint. Because of unfavourable winds, populations of the less dispersive Sicilian butterflies tended to differentiate into endemic variants or to maintain relict populations. Understanding these phenomena required the integration of DNA sequences, species traits and physical constraints for a large taxon at continental scale. Future studies may reveal if the patterns here shown for mitochondrial DNA are also reflected in the nuclear genome or, alternatively, are the product of limited female dispersal.

Population genomics reveals evolution and variation of Saccharomyces cerevisiae in the human and insects gut.

The quest to discover the variety of ecological niches inhabited by Saccharomyces cerevisiae has led to research in areas as diverse as wineries, oak trees and insect guts. The discovery of fungal communities in the human gastrointestinal tract suggested the host's gut as a potential reservoir for yeast adaptation. Here, we report the existence of yeast populations associated with the human gut (HG) that differ from those isolated from other human body sites. Phylogenetic analysis on 12 microsatellite loci and 1715 combined CDSs from whole-genome sequencing revealed three subclusters of HG strains with further evidence of clonal colonization within the host's gut. The presence of such subclusters was supported by other genomic features, such as copy number variation, absence/introgressions of CDSs and relative polymorphism frequency. Functional analysis of CDSs specific of the different subclusters suggested possible alterations in cell wall composition and sporulation features. The phenotypic analysis combined with immunological profiling of these strains further showed that sporulation was related with strain-specific genomic characteristics in the immune recognition pattern. We conclude that both genetic and environmental factors involved in cell wall remodelling and sporulation are the main drivers of adaptation in S. cerevisiae populations in the human gut.

Home economics in an oak gall: behavioural and chemical immune strategies against a fungal pathogen in Temnothorax ant nests.

Nest architecture is a fundamental character shaping immune strategies of social insects. The arboreal ant Temnothorax unifasciatus nests in cavities such as oak galls where the entire colony lives in a unique small chamber. In these conditions, physiological and behavioural strategies likely prevail over compartmentalisation and are presumably tuned with colony size. We designed two experiments to study chemical and behavioural immune strategies against the entomopathogenic fungus Metarhizium anisopliae in colonies of different sizes. First, we compared spore germination and length of germinal tubes inside artificial nests, designed to impede the contact between the ants and the fungus, in colonies of different size. In the absence of direct contact, Temnothorax unifasciatus colonies inhibit fungal growth inside their nests, presumably through volatile compounds. The analysis revealed a positive correlation between fungistatic activity and colony size, indicating that workers of smaller colonies do not invest a higher per capita effort in producing such substances compared to larger colonies. Second, we performed a removal experiment of contaminated and non-contaminated items introduced inside the nests of colonies of different size. Small colonies challenged with contaminated fibres showed an increased removal of all the items (both contaminated and non-contaminated) compared to small colonies challenged with non-contaminated fibres only. Conversely, larger colonies moved items regardless of the presence of the spores inside the nest. Colony size qualitatively affected removal of waste items showing a pathogen elicited reaction in small colonies to optimise the reduced workforce, while the removal behaviour in larger colonies revealed to be expressed constitutively.

Social wasps are a Saccharomyces mating nest.

The reproductive ecology of Saccharomyces cerevisiae is still largely unknown. Recent evidence of interspecific hybridization, high levels of strain heterozygosity, and prion transmission suggest that outbreeding occurs frequently in yeasts. Nevertheless, the place where yeasts mate and recombine in the wild has not been identified. We found that the intestine of social wasps hosts highly outbred S. cerevisiae strains as well as a rare S. cerevisiae×S. paradoxus hybrid. We show that the intestine of Polistes dominula social wasps favors the mating of S. cerevisiae strains among themselves and with S. paradoxus cells by providing a succession of environmental conditions prompting cell sporulation and spores germination. In addition, we prove that heterospecific mating is the only option for European S. paradoxus strains to survive in the gut. Taken together, these findings unveil the best hidden secret of yeast ecology, introducing the insect gut as an environmental alcove in which crosses occur, maintaining and generating the diversity of the ascomycetes.

Trans-generational immunization in the acrobat ant Crematogaster scutellaris.

Trans-generational immunization is defined as the transmission of an enhanced resistance to a pathogen from parents to offspring. By using the host-parasite system of the ant Crematogaster scutellaris and the entomopathogenic fungus Metarhizium anisopliae, we describe this phenomenon for the first time in ants. We exposed four groups of hibernating queens to different treatments (i) a non-lethal dose of live conidiospores in Triton, (ii) a dose of heat-killed conidiospores in Triton, (iii) a control Triton solution, and (iv) a naive control. We exposed their first workers to a high dose of conidiospores and measured mortality rates. Workers produced by queens exposed to live conidiospores survived longer than those belonging to the other groups, while exposure to Triton and dead spores had no effect. Starved workers showed a significantly higher mortality. The treatments did not influence queen mortality, nor the number of offspring they produced at the emergence of the first worker, showing no evidence of immunization costs-at least for these parameters in the first year of colony development. We propose that trans-generational immunization represents an important component of social immunity that could affect colony success, particularly during the critical phase of claustral foundation.

Use of genetic, climatic, and microbiological data to inform reintroduction of a regionally extinct butterfly.

Species reintroductions are increasingly used as means of mitigating biodiversity loss. Besides habitat quality at the site targeted for reintroduction, the choice of source population can be critical for success. The butterfly Melanargia russiae (Esper´s marbled white) was extirpated from Hungary over 100 years ago, and a reintroduction program has recently been approved. We used museum specimens of this butterfly, mitochondrial DNA data (mtDNA), endosymbiont screening, and climatic-similarity analyses to determine which extant populations should be used for its reintroduction. The species displayed 2 main mtDNA lineages across its range: 1 restricted to Iberia and southern France (Iberian lineage) and another found throughout the rest of its range (Eurasian lineage). These 2 lineages possessed highly divergent wsp alleles of the bacterial endosymbiont Wolbachia. The century-old Hungarian specimens represented an endemic haplotype belonging to the Eurasian lineage, differing by one mutation from the Balkan and eastern European populations. The Hungarian populations of M. russiae occurred in areas with a colder and drier climate relative to most sites with extant known populations. Our results suggest the populations used for reintroduction to Hungary should belong to the Eurasian lineage, preferably from eastern Ukraine (genetically close and living in areas with the highest climatic similarity). Materials stored in museum collections can provide unique opportunities to document historical genetic diversity and help direct conservation.

Social wasp intestines host the local phenotypic variability of Saccharomyces cerevisiae strains.

Nowadays, the presence of Saccharomyces cerevisiae has been assessed in both wild and human-related environments. Social wasps have been shown to maintain and vector S. cerevisiae among different environments. The availability of strains isolated from wasp intestines represents a striking opportunity to assess whether the strains found in wasp intestines are characterized by peculiar traits. We analysed strains isolated from the intestines of social wasps and compared them with strains isolated from other sources, all collected in a restricted geographic area. We evaluated the production of volatile metabolites during grape must fermentation, the resistance to different stresses and the ability to exploit various carbon sources. Wasp strains, in addition to representing a wide range of S. cerevisiae genotypes, also represent large part of the phenotypes characterizing the sympatric set of yeast strains; their higher production of acetic acid and ethyl acetate could reflect improved ability to attract insects. Our findings suggest that the relationship between yeasts and wasps should be preserved, to safeguard not only the natural variance of this microorganism but also the interests of wine-makers, who could take advantage from the exploitation of their phenotypic variability. Copyright © 2016 John Wiley & Sons, Ltd.

Integrative analyses unveil speciation linked to host plant shift in Spialia butterflies.

Discovering cryptic species in well-studied areas and taxonomic groups can have profound implications in understanding eco-evolutionary processes and in nature conservation because such groups often involve research models and act as flagship taxa for nature management. In this study, we use an array of techniques to study the butterflies in the Spialia sertorius species group (Lepidoptera, Hesperiidae). The integration of genetic, chemical, cytogenetic, morphological, ecological and microbiological data indicates that the sertorius species complex includes at least five species that differentiated during the last three million years. As a result, we propose the restitution of the species status for two taxa often treated as subspecies, Spialia ali (Oberthür, 1881) stat. rest. and Spialia therapne (Rambur, 1832) stat. rest., and describe a new cryptic species Spialia rosae Hernández-Roldán, Dapporto, Dinca, Vicente & Vila sp. nov. Spialia sertorius (Hoffmannsegg, 1804) and S. rosae are sympatric and synmorphic, but show constant differences in mitochondrial DNA, chemical profiles and ecology, suggesting that S. rosae represents a case of ecological speciation involving larval host plant and altitudinal shift, and apparently associated with Wolbachia infection. This study exemplifies how a multidisciplinary approach can reveal elusive cases of hidden diversity.

Role of social wasps in Saccharomyces cerevisiae ecology and evolution.

Saccharomyces cerevisiae is one of the most important model organisms and has been a valuable asset to human civilization. However, despite its extensive use in the last 9,000 y, the existence of a seasonal cycle outside human-made environments has not yet been described. We demonstrate the role of social wasps as vector and natural reservoir of S. cerevisiae during all seasons. We provide experimental evidence that queens of social wasps overwintering as adults (Vespa crabro and Polistes spp.) can harbor yeast cells from autumn to spring and transmit them to their progeny. This result is mirrored by field surveys of the genetic variability of natural strains of yeast. Microsatellites and sequences of a selected set of loci able to recapitulate the yeast strain's evolutionary history were used to compare 17 environmental wasp isolates with a collection of strains from grapes from the same region and more than 230 strains representing worldwide yeast variation. The wasp isolates fall into subclusters representing the overall ecological and industrial yeast diversity of their geographic origin. Our findings indicate that wasps are a key environmental niche for the evolution of natural S. cerevisiae populations, the dispersion of yeast cells in the environment, and the maintenance of their diversity. The close relatedness of several wasp isolates with grape and wine isolates reflects the crucial role of human activities on yeast population structure, through clonal expansion and selection of specific strains during the biotransformation of fermented foods, followed by dispersal mediated by insects and other animals.

A genetic atlas for the butterflies of continental Canada and United States.

Multi-locus genetic data for phylogeographic studies is generally limited in geographic and taxonomic scope as most studies only examine a few related species. The strong adoption of DNA barcoding has generated large datasets of mtDNA COI sequences. This work examines the butterfly fauna of Canada and United States based on 13,236 COI barcode records derived from 619 species. It compiles i) geographic maps depicting the spatial distribution of haplotypes, ii) haplotype networks (minimum spanning trees), and iii) standard indices of genetic diversity such as nucleotide diversity (π), haplotype richness (H), and a measure of spatial genetic structure (GST). High intraspecific genetic diversity and marked spatial structure were observed in the northwestern and southern North America, as well as in proximity to mountain chains. While species generally displayed concordance between genetic diversity and spatial structure, some revealed incongruence between these two metrics. Interestingly, most species falling in this category shared their barcode sequences with one at least other species. Aside from revealing large-scale phylogeographic patterns and shedding light on the processes underlying these patterns, this work also exposed cases of potential synonymy and hybridization.

A forest butterfly in sahara desert oases: isolation does not matter.

Numerous studies addressing the impact of habitat fragmentation on genetic diversity have been performed. In this study, we analyze the effects of a seemingly nonpermeable matrix on the population structure of the forest-dwelling butterfly Pararge aegeria in geographically isolated oases at the northern margin of the Sahara desert using microsatellites, morphological characters, and species distribution modeling. Results from all analyses are mostly congruent and reveal 1) a split between European and North African populations, 2) rather low divergence between populations from the eastern and western part of North Africa (Morocco vs. Tunisia), 3) a lack of differentiation between the oasis and Atlas Mountain populations, 4) as well as among the oasis populations, and 5) no reduction of genetic variability in oasis populations. However, one exception to this general trend resulted from the analyses of wing shape; wings of butterflies from oases are more elongated compared with those from the other habitats. This pattern of phenotypic divergence may suggest a recent colonization of the oasis habitats by individuals, which might be accompanied by a rather dispersive behavior. Species distribution modeling suggests a fairly recent reexpansion of the species' climatic niche starting in the Holocene at about 6000 before present. The combined results indicate a rather recent colonization of the oases by highly mobile individuals from genetically diverse founder populations. The colonization was likely followed by the expansion and persistence of these founder populations under relatively stable environmental conditions. This, together with low rates of gene flow, likely prevented differentiation of populations via drift and led to the maintenance of high genetic diversity.

Assessing the aesthetic attractivity of European butterflies: A web-based survey protocol.

Aesthetic attractivity stands as an underestimated yet fundamental feature of species in conservation biology, significantly driving disproportionate protection efforts towards charismatic species. Despite the evidence, few attempts sought to precisely quantify the impact of aesthetic attractivity in defining priority of species for conservation actions (e.g. inclusion in International Union for Conservation of Nature red lists and protection lists). This study protocol describes the setting of an online test (available from April 2022 to April 2023 at www.unveiling.eu) designed to i) quantify the aesthetic attractivity to humans of the 496 European butterfly species and ii) identify which features (both in the perceived animal and in the perceiver) influence the aesthetic attractivity of a given butterfly species. The test is divided in 5 sections (personal data, ranking, single morphological features, emotional engagement, dispositional variables) aimed at profiling the relation each participant has with the species examined. In the long-term, evaluating butterflies' aesthetic attractivity could facilitate the critical assessment of current conservation strategies, such as the process of selection of flag and umbrella species by research institutions, environmental associations and Non Governative Organizations. This is expected to provide the much-needed evidence to set up unbiased biodiversity conservation strategies and counteract the selective anthropogenic pressure which favours the extinction of unattractive species, being no or less protected compared to charismatic species.

A combined genetic-morphometric analysis unravels the complex biogeographical history of Polyommatus icarus and Polyommatus celina common blue butterflies.

Widespread species have the potential to reveal large-scale biogeographical patterns, as well as responses to environmental changes possibly unique to habitat generalists. This study presents a continental-scale phylogeographical analysis of Polyommatus icarus, one of the most common Palaearctic butterflies, and the morphologically and ecologically similar Polyommatus celina, a recently discovered cryptic species. By combining data from mitochondrial [cytochrome c oxidase subunit I (COI)] and nuclear [internal transcribed spacer (ITS2)] molecular markers with geometric morphometrics, we document a complex phylogeographical history for the two species. Despite morphological similarities, the genetic divergence between these two species is high (more than 5% at COI) and they are not sister species. For the first time, we show that P. celina occurs not only in North Africa but also in Europe, where it inhabits several west Mediterranean islands, as well as large parts of Iberia, where it occurs in parapatry with P. icarus. The two species appear to completely exclude each other on islands, but we provide morphological and molecular evidence that introgression occurred in the Iberian Peninsula. We discovered strongly diverged lineages that seem to represent relict populations produced by past range expansions and contractions: Crete and Iberian isolates for P. icarus, Balearics-Sardinia and Sicily-Lipari for P. celina. This study shows that a combined genetic-morphometric approach can shed light on cryptic diversity while providing the necessary resolution to reconstruct a fine-scale phylogeographical history of species at both spatial and temporal levels.

Cuticular Hydrocarbon Profile of Parasitic Beetles, Aethina tumida (Coleoptera: Nitidulidae).

Cuticular hydrocarbons (CHCs) cover insects' bodies and play important roles in chemical communication, including nestmate recognition, for social insects. To enter colonies of a social host species, parasites may acquire host-specific CHCs or covertly maintain their own CHC profile by lowering its quantity. However, the chemical profile of small hive beetles (SHBs), Aethina tumida, which are parasites of honey bee, Apis mellifera, colonies, and other bee nests, is currently unknown. Here, adults of SHB and honey bee host workers were collected from the same field colonies and their CHC profiles were analysed using GC-MS. The chemical profiles of field-sampled SHBs were also compared with those of host-naive beetles reared in the laboratory. Laboratory-reared SHBs differed in their CHC profiles from field-sampled ones, which showed a more similar, but ten-fold lower, generic host CHC profile compared to host workers. While the data confirm colony-specific CHCs of honey bee workers, the profile of field-collected SHBs was not colony-specific. Adult SHBs often commute between different host colonies, thereby possibly preventing the acquisition of a colony-specific CHC profiles. An ester was exclusive to both groups of SHBs and might constitute an intraspecific recognition cue. Our data suggest that SHBs do not use any finely tuned chemical strategy to conceal their presence inside host colonies and instead probably rely on their hard exoskeleton and defence behaviours.

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.