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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Molecular phylogeny, biogeography, and host plant shifts in the bee genus Melitta (Hymenoptera: Anthophila).

New molecular studies suggested that the family Melittidae is either a paraphyletic group from which all the other bees are derived, or the sister clade to all other existing bees. Studying the historical biogeography and evolution of each major lineage within this group is a key step to understand the origin and early radiation of bees. Melitta is the largest genus of melittid bees, for which a robust molecular phylogeny and a biogeographic analysis are still lacking. Here, we derive a phylogenetic hypothesis from the sequences of seven independent DNA fragments of mitochondrial and nuclear origin. This phylogenetic hypothesis is then used to infer the evolution of the species range and of the host-plant shifts in Melitta. Our results confirmed the monophyly of Melitta, but did not recover all previously defined clades within the genus. We propose new taxa by splitting the genus in three subgenera (including two new subgenera described in the Appendix: Afromelitta subgen. nov., Plesiomelitta subgen. nov.) and describe two new species: Melitta avontuurensis sp. n. and M. richtersveldensis sp. n. Regarding the evolution of host-plant use, our analysis suggests that all species currently specialized on one plant family originated from an ancestor that was specialized on Fabaceae plants. The inferred biogeographic history for the genus supported an African origin. In concordance with previous studies identifying Africa as the geographic origin for many clades of bees, our data bring new evidence for an African origin of melittid bees.

Large-scale pattern of genetic differentiation within African rainforest trees: insights on the roles of ecological gradients and past climate changes on the evolution of Erythrophleum spp (Fabaceae).

BACKGROUND: The evolutionary events that have shaped biodiversity patterns in the African rainforests are still poorly documented. Past forest fragmentation and ecological gradients have been advocated as important drivers of genetic differentiation but their respective roles remain unclear. Using nuclear microsatellites (nSSRs) and chloroplast non-coding sequences (pDNA), we characterised the spatial genetic structure of Erythrophleum (Fabaceae) forest trees in West and Central Africa (Guinea Region, GR). This widespread genus displays a wide ecological amplitude and taxonomists recognize two forest tree species, E. ivorense and E. suaveolens, which are difficult to distinguish in the field and often confused. RESULTS: Bayesian-clustering applied on nSSRs of a blind sample of 648 specimens identified three major gene pools showing no or very limited introgression. They present parapatric distributions correlated to rainfall gradients and forest types. One gene pool is restricted to coastal evergreen forests and corresponds to E. ivorense; a second one is found in gallery forests from the dry forest zone of West Africa and North-West Cameroon and corresponds to West-African E. suaveolens; the third gene pool occurs in semi-evergreen forests and corresponds to Central African E. suaveolens. These gene pools have mostly unique pDNA haplotypes but they do not form reciprocally monophyletic clades. Nevertheless, pDNA molecular dating indicates that the divergence between E. ivorense and Central African E. suaveolens predates the Pleistocene. Further Bayesian-clustering applied within each major gene pool identified diffuse genetic discontinuities (minor gene pools displaying substantial introgression) at a latitude between 0 and 2°N in Central Africa for both species, and at a longitude between 5° and 8°E for E. ivorense. Moreover, we detected evidence of past population declines which are consistent with historical habitat fragmentation induced by Pleistocene climate changes. CONCLUSIONS: Overall, deep genetic differentiation (major gene pools) follows ecological gradients that may be at the origin of speciation, while diffuse differentiation (minor gene pools) are tentatively interpreted as the signature of past forest fragmentation induced by past climate changes.

Climate change and the spread of vector-borne diseases: using approximate Bayesian computation to compare invasion scenarios for the bluetongue virus vector Culicoides imicola in Italy.

Bluetongue (BT) is a commonly cited example of a disease with a distribution believed to have recently expanded in response to global warming. The BT virus is transmitted to ruminants by biting midges of the genus Culicoides, and it has been hypothesized that the emergence of BT in Mediterranean Europe during the last two decades is a consequence of the recent colonization of the region by Culicoides imicola and linked to climate change. To better understand the mechanism responsible for the northward spread of BT, we tested the hypothesis of a recent colonization of Italy by C. imicola, by obtaining samples from more than 60 localities across Italy, Corsica, Southern France, and Northern Africa (the hypothesized source point for the recent invasion of C. imicola), and by genotyping them with 10 newly identified microsatellite loci. The patterns of genetic variation within and among the sampled populations were characterized and used in a rigorous approximate Bayesian computation framework to compare three competing historical hypotheses related to the arrival and establishment of C. imicola in Italy. The hypothesis of an ancient presence of the insect vector was strongly favoured by this analysis, with an associated P ≥ 99%, suggesting that causes other than the northward range expansion of C. imicola may have supported the emergence of BT in southern Europe. Overall, this study illustrates the potential of molecular genetic markers for exploring the assumed link between climate change and the spread of diseases.

Pre-adaptations and the evolution of pollination by sexual deception: Cope's rule of specialization revisited.

Pollination by sexual deception is arguably one of the most unusual liaisons linking plants and insects, and perhaps the most illustrative example of extreme floral specialization in angiosperms. While considerable progress has been made in understanding the floral traits involved in sexual deception, less is known about how this remarkable mimicry system might have arisen, the role of pre-adaptations in promoting its evolution and its extent as a pollination mechanism outside the few groups of plants (primarily orchids) where it has been described to date. In the Euro-Mediterranean region, pollination by sexual deception is traditionally considered to be the hallmark of the orchid genus Ophrys. Here, we introduce two new cases outside of Ophrys, in plant groups dominated by generalized, shelter-mimicking species. On the basis of phylogenetic reconstructions of ancestral pollination strategies, we provide evidence for independent and bidirectional evolutionary transitions between generalized (shelter mimicry) and specialized (sexual deception) pollination strategies in three groups of flowering plants, and suggest that pseudocopulation has evolved from pre-adaptations (floral colours, shapes and odour bouquets) that selectively attract male pollinators through shelter mimicry. These findings, along with comparative analyses of floral traits (colours and scents), shed light on particular phenotypic changes that might have fuelled the parallel evolution of these extraordinary pollination strategies. Collectively, our results provide the first substantive insights into how pollination sexual deception might have evolved in the Euro-Mediterranean region, and demonstrate that even the most extreme cases of pollinator specialization can reverse to more generalized interactions, breaking 'Cope's rule of specialization'.

Trees and/or networks to display intraspecific DNA sequence variation?

Phylogenetic trees and networks are both used in the scientific literature to display DNA sequence variation at the intraspecific level. Should we rather use trees or networks? I argue that the process of inferring the most parsimonious genealogical relationships among a set of DNA sequences should be dissociated from the problem of displaying this information in a graph. A network graph is probably more appropriate than a strict consensus tree if many alternative, equally most parsimonious, genealogies are to be included. Within the maximum parsimony framework, current phylogenetic inference and network-building algorithms are both unable to guarantee the finding of all most parsimonious (MP) connections. In fact, each approach can find MP connections that the other does not. Although it should be possible to improve at least the maximum parsimony approach, current implementations of these algorithms are such that it is advisable to use both approaches to increase the probability of finding all possible MP connections among a set of DNA sequences.

Conflicting mitochondrial and nuclear phylogeographic signals and evolution of host-plant shifts in the boreo-montane leaf beetle Chrysomela lapponica.

We conducted a phylogeographic study on the cold-adapted leaf beetle Chrysomela lapponica, that feeds on willow or birch, by sampling several populations throughout most of the geographic distribution of the species, and by sequencing for each individual one mitochondrial and two nuclear DNA fragments. Patterns of DNA sequence variation from the mitochondrial and nuclear loci, as displayed in the median-joining networks, appear to display contradicting historical signal: a deep genealogical divergence is observed with the mitochondrial genome between the Alpine population and all other populations found in the Euro-Siberian distribution of the species, that is completely absent with both nuclear loci. We use coalescence simulations of DNA sequence evolution to test the hypothesis that this apparent conflict is compatible with a neutral model of sequence evolution (i.e., to check whether the stochastic nature of the coalescence process can explain these patterns). Because the simulations show that this is highly unlikely, we consider two alternative hypotheses: (1) introgression of the mitochondrial genome of another species and (2) the effect of natural selection. Although introgression is the most plausible explanation, we fail to identify the source species of the introgressed mitochondrial genome among all known species closely related to C. lapponica. We therefore suggest that the putative introgression event is ancient and the source species is either extinct or currently outside the geographic range of C. lapponica explored in this study. The observed DNA sequence variation also suggests that a host-plant shift from willow to birch has occurred recently and independently in each of the three birch-feeding populations. This emphasizes further the relative ease with which these beetles can escape their ancestral host-plant specialization on willow, but shows at the same time that host-plant shifts are highly constrained, as they only occur between willow and birch.

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.