At the Gate of Mutualism: Identification of Genomic Traits Predisposing to Insect-Bacterial Symbiosis in Pathogenic Strains of the Aphid Symbiont Serratia symbiotica.

Mutualistic associations between insects and heritable bacterial symbionts are ubiquitous in nature. The aphid symbiont Serratia symbiotica is a valuable candidate for studying the evolution of bacterial symbiosis in insects because it includes a wide diversity of strains that reflect the diverse relationships in which bacteria can be engaged with insects, from pathogenic interactions to obligate intracellular mutualism. The recent discovery of culturable strains, which are hypothesized to resemble the ancestors of intracellular strains, provide an opportunity to study the mechanisms underlying bacterial symbiosis in its early stages. In this study, we analyzed the genomes of three of these culturable strains that are pathogenic to aphid hosts, and performed comparative genomic analyses including mutualistic host-dependent strains. All three genomes are larger than those of the host-restricted S. symbiotica strains described so far, and show significant enrichment in pseudogenes and mobile elements, suggesting that these three pathogenic strains are in the early stages of the adaptation to their host. Compared to their intracellular mutualistic relatives, the three strains harbor a greater diversity of genes coding for virulence factors and metabolic pathways, suggesting that they are likely adapted to infect new hosts and are a potential source of metabolic innovation for insects. The presence in their genomes of secondary metabolism gene clusters associated with the production of antimicrobial compounds and phytotoxins supports the hypothesis that S. symbiotia symbionts evolved from plant-associated strains and that plants may serve as intermediate hosts. Mutualistic associations between insects and bacteria are the result of independent transitions to endosymbiosis initiated by the acquisition of environmental progenitors. In this context, the genomes of free-living S. symbiotica strains provide a rare opportunity to study the inventory of genes held by bacterial associates of insects that are at the gateway to a host-dependent lifestyle.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Genome Sequencing from Post-Mortem Formalin-Fixed, Paraffin-Embedded Lung Tissues.

Implementation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) testing in the daily practice of pathology laboratories requires procedure adaptation to formalin-fixed, paraffin-embedded (FFPE) samples. So far, one study reported the feasibility of SARS-CoV-2 genome sequencing on FFPE tissues with only one contributory case of two. This study optimized SARS-CoV-2 genome sequencing using the Ion AmpliSeq SARS-CoV-2 Panel on 22 FFPE lung tissues from 16 deceased coronavirus disease 2019 (COVID-19) patients. SARS-CoV-2 was detected in all FFPE blocks using a real-time RT-qPCR targeting the E gene with crossing point (Cp) values ranging from 16.02 to 34.16. Sequencing was considered as contributory (i.e. with a uniformity >55%) for 17 FFPE blocks. Adapting the number of target amplification PCR cycles according to the RT-qPCR Cp values allowed optimization of the sequencing quality for the contributory blocks (i.e. 20 PCR cycles for blocks with a Cp value <28 and 25 PCR cycles for blocks with a Cp value between 28 and 30). Most blocks with a Cp value >30 were non-contributory. Comparison of matched frozen and FFPE tissues revealed discordance for only three FFPE blocks, all with a Cp value >28. Variant identification and clade classification was possible for 13 patients. This study validates SARS-CoV-2 genome sequencing on FFPE blocks and opens the possibility to explore correlation between virus genotype and histopathologic lesions.

Whole-genome sequencing reveals asymmetric introgression between two sister species of cold-resistant leaf beetles.

A large number of genetic variation studies have identified cases of mitochondrial genome introgression in animals, indicating that reproductive barriers among closely related species are often permeable. Because of its sheer size, the impact of hybridization on the evolution of the nuclear genome is more difficult to apprehend. Only a few studies have explored it recently thanks to recent progress in DNA sequencing and genome assembly. Here, we analysed whole-genome sequence variation among multiple individuals of two sister species of leaf beetles inside their hybrid zone, in which asymmetric mitochondrial genome introgression had previously been established. We used a machine learning approach based on computer simulations for training to identify regions of the nuclear genome that were introgressed. We inferred asymmetric introgression of ≈2% of the genome, in the same direction that was observed for the mitochondrial genome. Because a previous study based on a reduced-representation sequencing approach was not able to detect this introgression, we conclude that whole-genome sequencing is necessary when the fraction of the introgressed genome is small. We also analysed the whole-genome sequence of a hybrid individual, demonstrating that hybrids have the capacity to backcross with the species for which virtually no introgression was observed. Our data suggest that one species has recently invaded the range of the other and/or some alleles that where transferred from the invaded into the invading species could be under positive selection and may have favoured the adaptation of the invading species to the Alpine environment.

Genome Assembly of the Cold-Tolerant Leaf Beetle Gonioctena quinquepunctata, an Important Resource for Studying Its Evolution and Reproductive Barriers between Species.

Coleoptera is the most species-rich insect order, yet is currently underrepresented in genomic databases. An assembly was generated for ca. 1.7-Gb genome of the leaf beetle Gonioctena quinquepunctata by first assembling long-sequence reads (Oxford Nanopore; ± 27-fold coverage) and subsequently polishing the resulting assembly with short sequence reads (Illumina; ± 85-fold coverage). The unusually large size (most Coleoptera species are associated with a reported size below 1 Gb) was at least partially attributed to the presence of a large fraction of repeated elements (73.8%). The final assembly was characterized by an N50 length of 432 kb and a BUSCO score of 95.5%. The heterozygosity rate was ±0.6%. Automated genome annotation informed by RNA-Seq resulted in 40,568 predicted proteins, which is much larger than the typical range 17,000-23,000 predicted for other Coleoptera. However, no evidence of a genome duplication was detected. This new reference genome will contribute to our understanding of genetic variation in the Coleoptera. Among others, it will also allow exploring reproductive barriers between species, investigating introgression in the nuclear genome, and identifying genes involved in resistance to extreme climate conditions.

Evolutionary history of inquiline social parasitism in Plagiolepis ants.

Social parasitism, i.e. the parasitic dependence of a social species on another free-living social species, is one of the most intriguing phenomena in social insects. It has evolved to various levels, the most extreme form being inquiline social parasites which have lost the worker caste, and produce only male and female sexual offspring that are reared by the host worker force. The inquiline syndrome has been reported in 4 species within the ant genus Plagiolepis, in Europe. Whether inquiline social parasitism evolved once or multiple times within the genus remains however unknown. To address this question, we generated data for 5 inquiline social parasites - 3 species previously described and 2 unidentified species - and their free-living hosts from Europe, and we inferred their phylogenetic relationships. We tested Emery's rule, which predicts that inquiline social parasites and their hosts are close relatives. Our results show that inquiline parasitism evolved independently at least 5 times in the genus. Furthermore, we found that all inquilines were associated with one of the descendants of their most related free-living species, suggesting sympatric speciation is the main process leading to the emergence of the parasitic species, consistent with the stricter version of Emery's rule.

What do tropical cryptogams reveal? Strong genetic structure in Amazonian bryophytes.

Lowland tropical bryophytes have been perceived as excellent dispersers. In such groups, the inverse isolation hypothesis proposes that spatial genetic structure is erased beyond the limits of short-distance dispersal. Here, we determine the influence of environmental variation and geographic barriers on the spatial genetic structure of a widely dispersed and phylogenetically independent sample of Amazonian bryophytes. Single nucleotide polymorphism data were produced from a restriction site-associated DNA sequencing protocol for 10 species and analyzed through F-statistics and Mantel tests. Neither isolation-by-environment nor the impact of geographic barriers were recovered from the analyses. However, significant isolation-by-distance patterns were observed for 8 out of the 10 investigated species beyond the scale of short-distance dispersal (> 1 km), offering evidence contrary to the inverse isolation hypothesis. Despite a cadre of life-history traits and distributional patterns suggesting that tropical bryophytes are highly vagile, our analyses reveal spatial genetic structures comparable to those documented for angiosperms, whose diaspores are orders of magnitude larger. Dispersal limitation for tropical bryophytes flies in the face of traditional assumptions regarding their dispersal potential, and suggests that the plight of this component of cryptic biodiversity is more dire than previously considered in light of accelerated forest fragmentation in the Amazon.

Mitonuclear mismatch alters performance and reproductive success in naturally introgressed populations of a montane leaf beetle.

Coordination between nuclear and mitochondrial genomes is critical to metabolic processes underlying animals' ability to adapt to local environments, yet consequences of mitonuclear interactions have rarely been investigated in populations where individuals with divergent mitochondrial and nuclear genomes naturally interbreed. Genetic variation in the leaf beetle Chrysomela aeneicollis was assessed along a latitudinal thermal gradient in California's Sierra Nevada. Variation at mitochondrial cytochrome oxidase II (COII) and the nuclear gene phosphoglucose isomerase (PGI) shows concordance and was significantly greater along a 65 km transect than 10 other loci. STRUCTURE analyses using neutral loci identified a southern and northern subpopulation, which interbreed in the central drainage Bishop Creek. COII and PGI were used as indicators of mitochondrial and nuclear genetic variation in field and laboratory experiments conducted on beetles from this admixed population. Fecundity, larval development rate, running speed and male mating frequency were higher for beetles with geographically "matched" than "mismatched" mitonuclear genotypes. Effects of mitonuclear mismatch were largest for individuals with northern nuclear genotypes possessing southern mitochondria and were most pronounced after heat treatment or at high elevation. These findings suggest that mitonuclear incompatibility diminishes performance and reproductive success in nature, effects that could intensify at environmental extremes.

Unraveling heteroplasmy patterns with NOVOPlasty.

Heteroplasmy, the existence of multiple mitochondrial haplotypes within an individual, has been studied across different scientific fields. Mitochondrial genome polymorphisms have been linked to multiple severe disorders and are of interest to evolutionary studies and forensic science. Before the development of massive parallel sequencing (MPS), most studies of mitochondrial genome variation were limited to short fragments and to heteroplasmic variants associated with a relatively high frequency (>10%). By utilizing ultra-deep sequencing, it has now become possible to uncover previously undiscovered patterns of intra-individual polymorphisms. Despite these technological advances, it is still challenging to determine the origin of the observed intra-individual polymorphisms. We therefore developed a new method that not only detects intra-individual polymorphisms within mitochondrial and chloroplast genomes more accurately, but also looks for linkage among polymorphic sites by assembling the sequence around each detected polymorphic site. Our benchmark study shows that this method is capable of detecting heteroplasmy more accurately than any method previously available and is the first tool that is able to completely or partially reconstruct the sequence for each mitochondrial haplotype (allele). The method is implemented in our open source software NOVOPlasty that can be downloaded at https://github.com/ndierckx/NOVOPlasty.

Limited gene exchange between two sister species of leaf beetles within a hybrid zone in the Alps.

Analysing genomic variation within and between sister species is a first step towards understanding species boundaries. We focused on two sister species of cold-resistant leaf beetles, Gonioctena quinquepunctata and G. intermedia, whose ranges overlap in the Alps. A previous study of DNA sequence variation had revealed multiple instances of mitochondrial genome introgression in this region, suggesting recent hybridization between the two species. To evaluate the extent of gene exchange resulting from these hybridization events, we sampled individuals of both species inside and outside the hybrid zone and analysed genomic variation among them using RAD-seq markers. Individual levels of introgression in the nuclear genome were estimated first by defining species-specific SNPs (displaying a fixed difference between species) a priori and second by using model-based methods. Both types of analyses indicated little gene exchange, if any, between species at the level of the nuclear genome. Whereas the first method suggested slightly more gene flow, we argue that it has likely overestimated introgression in the phylogeographic context of this study. We conclude that strong intrinsic barriers prevent genetic exchange at the level of the nuclear genome between the two species. The apparent discrepancy observed between introgression occurring in the nuclear and mitochondrial genomes could be explained by selection acting in favour of the latter. Also, these results have consequences for the phylogeographic study of each species, since we can assume that genetic diversity in the overlapping portion of their ranges is not the product of introgression.

Molecular chaperoning helps safeguarding mitochondrial integrity and motor functions in the Sahara silver ant Cataglyphis bombycina.

The Sahara silver ant Cataglyphis bombycina is one of the world's most thermotolerant animals. Workers forage for heat-stricken arthropods during the hottest part of the day, when temperatures exceed 50 °C. However, the physiological adaptations needed to cope with such harsh conditions remain poorly studied in this desert species. Using transcriptomics, we screened for the most heat-responsive transcripts of C. bombycina with aim to better characterize the molecular mechanisms involved with macromolecular stability and cell survival to heat-stress. We identified 67 strongly and consistently expressed transcripts, and we show evidences of both evolutionary selection and specific heat-induction of mitochondrial-related molecular chaperones that have not been documented in Formicidae so far. This indicates clear focus of the silver ant's heat-shock response in preserving mitochondrial integrity and energy production. The joined induction of small heat-shock proteins likely depicts the higher requirement of this insect for proper motor function in response to extreme burst of heat-stresses. We discuss how those physiological adaptations may effectively help workers resist and survive the scorching heat and burning ground of the midday Sahara Desert.

Widespread co-occurrence of two distantly related mitochondrial genomes in individuals of the leaf beetle Gonioctena intermedia.

Mitochondrial genome heteroplasmy-the presence of more than one genomic variant in individuals-is considered only occasional in animals, and most often involves molecules differing only by a few recent mutations. Thanks to new sequencing technologies, a large number of DNA fragments from a single individual can now be sequenced and visualized separately, allowing new insights into intra-individual mitochondrial genome variation. Here, we report evidence from both (i) massive parallel sequencing (MPS) of genomic extracts and (ii) Sanger sequencing of PCR products, for the widespread co-occurrence of two distantly related (greater than 1% nucleotide divergence, excluding the control region) mitochondrial genomes in individuals of a natural population of the leaf beetle Gonioctena intermedia Sanger sequencing of PCR products using universal primers previously failed to identify heteroplasmy in this population. Its occurrence was detected with MPS data and may have important implications for evolutionary studies. It suggests the need to re-evaluate, using MPS techniques, the proportion of animal species displaying heteroplasmy.

Glacial survival of trophically linked boreal species in northern Europe.

Whether non-arctic species persisted in northern Europe during the Last Glacial Maximum (LGM) is highly debated. Until now, the debate has mostly focused on plants, with little consideration for other groups of organisms, e.g. the numerous plant-dependent insect species. Here, we study the late-Quaternary evolution of the European range of a boreo-montane leaf beetle, Gonioctena intermedia, which feeds exclusively on the boreal and temperate trees Prunus padus and Sorbus aucuparia Using species distribution models, we estimated the congruence between areas of past and present suitable climate for this beetle and its host plants. Then we derived historical hypotheses from the congruent range estimates, and evaluated their compatibility with observed DNA sequence variation at five independent loci. We investigated compatibility using computer simulations of population evolution under various coalescence models. We find strong evidence for range modifications in response to late-Quaternary climate changes, and support for the presence of populations of G. intermedia in northern Europe since the beginning of the last glaciation. The presence of a co-dependent insect in the region through the LGM provides new evidence supporting the glacial survival of cold-tolerant tree species in northern Europe.

NOVOPlasty: de novo assembly of organelle genomes from whole genome data.

The evolution in next-generation sequencing (NGS) technology has led to the development of many different assembly algorithms, but few of them focus on assembling the organelle genomes. These genomes are used in phylogenetic studies, food identification and are the most deposited eukaryotic genomes in GenBank. Producing organelle genome assembly from whole genome sequencing (WGS) data would be the most accurate and least laborious approach, but a tool specifically designed for this task is lacking. We developed a seed-and-extend algorithm that assembles organelle genomes from whole genome sequencing (WGS) data, starting from a related or distant single seed sequence. The algorithm has been tested on several new (Gonioctena intermedia and Avicennia marina) and public (Arabidopsis thaliana and Oryza sativa) whole genome Illumina data sets where it outperforms known assemblers in assembly accuracy and coverage. In our benchmark, NOVOPlasty assembled all tested circular genomes in less than 30 min with a maximum memory requirement of 16 GB and an accuracy over 99.99%. In conclusion, NOVOPlasty is the sole de novo assembler that provides a fast and straightforward extraction of the extranuclear genomes from WGS data in one circular high quality contig. The software is open source and can be downloaded at https://github.com/ndierckx/NOVOPlasty.

Mitochondrial DNA hyperdiversity and its potential causes in the marine periwinkle Melarhaphe neritoides (Mollusca: Gastropoda).

We report the presence of mitochondrial DNA (mtDNA) hyperdiversity in the marine periwinkle Melarhaphe neritoides (Linnaeus, 1758), the first such case among marine gastropods. Our dataset consisted of concatenated 16S-COI-Cytb gene fragments. We used Bayesian analyses to investigate three putative causes underlying genetic variation, and estimated the mtDNA mutation rate, possible signatures of selection and the effective population size of the species in the Azores archipelago. The mtDNA hyperdiversity in M. neritoides is characterized by extremely high haplotype diversity (Hd = 0.999 ± 0.001), high nucleotide diversity (π = 0.013 ± 0.001), and neutral nucleotide diversity above the threshold of 5% (πsyn = 0.0677). Haplotype richness is very high even at spatial scales as small as 100m2. Yet, mtDNA hyperdiversity does not affect the ability of DNA barcoding to identify M. neritoides. The mtDNA hyperdiversity in M. neritoides is best explained by the remarkably high mutation rate at the COI locus (µ = 5.82 × 10-5 per site per year or µ = 1.99 × 10-4 mutations per nucleotide site per generation), whereas the effective population size of this planktonic-dispersing species is surprisingly small (Ne = 5, 256; CI = 1,312-3,7495) probably due to the putative influence of selection. Comparison with COI nucleotide diversity values in other organisms suggests that mtDNA hyperdiversity may be more frequently linked to high µ values and that mtDNA hyperdiversity may be more common across other phyla than currently appreciated.

High migration rates shape the postglacial history of amphi-Atlantic bryophytes.

Paleontological evidence and current patterns of angiosperm species richness suggest that European biota experienced more severe bottlenecks than North American ones during the last glacial maximum. How well this pattern fits other plant species is less clear. Bryophytes offer a unique opportunity to contrast the impact of the last glacial maximum in North America and Europe because about 60% of the European bryoflora is shared with North America. Here, we use population genetic analyses based on approximate Bayesian computation on eight amphi-Atlantic species to test the hypothesis that North American populations were less impacted by the last glacial maximum, exhibiting higher levels of genetic diversity than European ones and ultimately serving as a refugium for the postglacial recolonization of Europe. In contrast with this hypothesis, the best-fit demographic model involved similar patterns of population size contractions, comparable levels of genetic diversity and balanced migration rates between European and North American populations. Our results thus suggest that bryophytes have experienced comparable demographic glacial histories on both sides of the Atlantic. Although a weak, but significant genetic structure was systematically recovered between European and North American populations, evidence for migration from and towards both continents suggests that amphi-Atlantic bryophyte population may function as a metapopulation network. Reconstructing the biogeographic history of either North American or European bryophyte populations therefore requires a large, trans-Atlantic geographic framework.

Comparative multilocus phylogeography of two Palaearctic spruce bark beetles: influence of contrasting ecological strategies on genetic variation.

While phylogeographic patterns of organisms are often interpreted through past environmental disturbances, mediated by climate changes, and geographic barriers, they may also be strongly influenced by species-specific traits. To investigate the impact of such traits, we focused on two Eurasian spruce bark beetles that share a similar geographic distribution, but differ in their ecology and reproduction. Ips typographus is an aggressive tree-killing species characterized by strong dispersal, whereas Dendroctonus micans is a discrete inbreeding species (sib mating is the rule), parasite of living trees and a poor disperser. We compared genetic variation between the two species over both beetles' entire range in Eurasia with five independent gene fragments, to evaluate whether their intrinsic differences could have an influence over their phylogeographic patterns. We highlighted widely divergent patterns of genetic variation for the two species and argue that the difference is indeed largely compatible with their contrasting dispersal strategies and modes of reproduction. In addition, genetic structure in I. typographus divides European populations in a northern and a southern group, as was previously observed for its host plant, and suggests past allopatric divergence. A long divergence time was estimated between East Asian and other populations of both species, indicating their long-standing presence in Eurasia, prior to the last glacial maximum. Finally, the strong population structure observed in D. micans for the mitochondrial locus provides insights into the recent colonization history of this species, from its native European range to regions where it was recently introduced.

Approximate Bayesian Computation Reveals the Crucial Role of Oceanic Islands for the Assembly of Continental Biodiversity.

The perceived low levels of genetic diversity, poor interspecific competitive and defensive ability, and loss of dispersal capacities of insular lineages have driven the view that oceanic islands are evolutionary dead ends. Focusing on the Atlantic bryophyte flora distributed across the archipelagos of the Azores, Madeira, the Canary Islands, Western Europe, and northwestern Africa, we used an integrative approach with species distribution modeling and population genetic analyses based on approximate Bayesian computation to determine whether this view applies to organisms with inherent high dispersal capacities. Genetic diversity was found to be higher in island than in continental populations, contributing to mounting evidence that, contrary to theoretical expectations, island populations are not necessarily genetically depauperate. Patterns of genetic variation among island and continental populations consistently fitted those simulated under a scenario of de novo foundation of continental populations from insular ancestors better than those expected if islands would represent a sink or a refugium of continental biodiversity. We, suggest that the northeastern Atlantic archipelagos have played a key role as a stepping stone for transoceanic migrants. Our results challenge the traditional notion that oceanic islands are the end of the colonization road and illustrate the significant role of oceanic islands as reservoirs of novel biodiversity for the assembly of continental floras.

Impact of past climatic changes and resource availability on the population demography of three food-specialist bees.

Past climate change is known to have strongly impacted current patterns of genetic variation of animals and plants in Europe. However, ecological factors also have the potential to influence demographic history and thus patterns of genetic variation. In this study, we investigated the impact of past climate, and also the potential impact of host plant species abundance, on intraspecific genetic variation in three codistributed and related specialized solitary bees of the genus Melitta with very similar life history traits and dispersal capacities. We sequenced five independent loci in samples collected from the three species. Our analyses revealed that the species associated with the most abundant host plant species (Melitta leporina) displays unusually high genetic variation, to an extent that is seldom reported in phylogeographic studies of animals and plants. This suggests a potential role of food resource abundance in determining current patterns of genetic variation in specialized herbivorous insects. Patterns of genetic variation in the two other species indicated lower overall levels of diversity, and that M. nigricans could have experienced a recent range expansion. Ecological niche modelling of the three Melitta species and their main host plant species suggested a strong reduction in range size during the last glacial maximum. Comparing observed sequence data with data simulated using spatially explicit models of coalescence suggests that M. leporina recovered a range and population size close to their current levels at the end of the last glaciation, and confirms recent range expansion as the most likely scenario for M. nigricans. Overall, this study illustrates that both demographic history and ecological factors may have contributed to shape current phylogeographic patterns.

Comparing phylogeographic hypotheses by simulating DNA sequences under a spatially explicit model of coalescence.

Computer simulations of genetic data are increasingly used to investigate the impact of complex historical scenarios on patterns of genetic variation. Yet, in most empirical studies, relatively large portions of species ranges are often treated as panmictic populations, ignoring the underlying spatial context. In some cases, however, a more accurate spatial model is required. We use a spatially explicit model of coalescence (easily constructed by overlaying a two-dimensional grid on maps displaying an estimate of past and current species ranges) to evaluate the potential of several summary statistics to differentiate three typical phylogeographic scenarios. We first explore the variation of each summary statistic within the boundaries of each phylogeographic scenario, and identify those that appear most promising for a comparison of historical scenarios and/or to infer historical parameters. We then combine a selected set of summary statistics in a single chi-square statistic and evaluate whether it can be used to differentiate past geographic fragmentation or range expansion from a simple scenario of isolation by distance. We also investigate the benefits of using a spatially explicit model by comparing its performance to alternative models that are less spatially explicit (lower geographic resolution). The results identify conditions in which each summary statistic is useful to infer the evolution of a species range, and allow us to validate our spatially explicit model of coalescence and our procedure to compare simulated and observed sequence data. We also provide a detailed description of the spatially explicit model of coalescence used, which is currently lacking.

Multi-locus DNA sequence variation in a complex of four leaf beetle species with parapatric distributions: mitochondrial and nuclear introgressions reveal recent hybridization.

While the importance for including multiple independent loci in phylogeographic studies is largely acknowledged, a majority of these still focus on a single species. We combine the study of multilocus DNA sequence variation (one mitochondrial and four unlinked nuclear fragments) at both the inter- and intraspecific levels to explore the evolutionary history of four closely related specialized herbivorous insects (Chrysomelidae, genus Gonioctena). Analyzing the sequences for samples collected across their entire European range allows us to (1) characterize the genetic boundaries among species, i.e. the degree of lineage sorting, (2) infer their phylogenetic relationships and (3) explore reproductive barriers among them in regions where their ranges overlap. For two sister species, we identify multiple independent cases of unidirectional transfer of genetic material (introgression) at both mitochondrial and nuclear loci, demonstrating recent hybridization between both species in the overlapping regions of their range. The highlighted pattern of genetic variation suggests that Gonioctena intermedia expanded its range into that of G. quinquepunctata, and that both species may experience mutual exclusion. Overall, this study illustrates that interpreting intraspecific genetic variation for the purpose of evolutionary inference without the broader context of the other closely related species could lead to erroneous conclusions.