Mitochondrial DNA (mtDNA) sequencing has led to an unprecedented rise in the identification of cryptic species. However, it is widely acknowledged that nuclear DNA (nuDNA) sequence data are also necessary to properly define species boundaries. Next generation sequencing techniques provide a wealth of nuclear genomic data, which can be used to ascertain both the evolutionary history and taxonomic status of putative cryptic species. Here, we focus on the intriguing case of the butterfly Thymelicus sylvestris (Lepidoptera: Hesperiidae). We identified six deeply diverged mitochondrial lineages; three distributed all across Europe and found in sympatry, suggesting a potential case of cryptic species. We then sequenced these six lineages using double‐digest restriction‐site associated DNA sequencing (ddRADseq). Nuclear genomic loci contradicted mtDNA patterns and genotypes generally clustered according to geography, i.e., a pattern expected under the assumption of postglacial recolonization from different refugia. Further analyses indicated that this strong mtDNA/nuDNA discrepancy cannot be explained by incomplete lineage sorting, sex‐biased asymmetries, NUMTs, natural selection, introgression or Wolbachia‐mediated genetic sweeps. We suggest that this mitonuclear discordance was caused by long periods of geographic isolation followed by range expansions, homogenizing the nuclear but not the mitochondrial genome. These results highlight T. sylvestris as a potential case of multiple despeciation and/or lineage fusion events. We finally argue, since mtDNA and nuDNA do not necessarily follow the same mechanisms of evolution, their respective evolutionary history reflects complementary aspects of past demographic and biogeographic events.

AimWe describe fine‐scale diversity patterns of the entire butterfly fauna occurring on the Tuscan Archipelago. By assessing the traits associated with population diversification, haplotype uniqueness and extinction, we aim to identify the factors determining the origin and maintenance of genetic diversity, and population vulnerability to environmental changes.LocationTuscan Archipelago, Sardinia, Tuscany (Italy) and Corsica (France).MethodsWe built a mtDNA dataset (1,303 COI sequences) for the 52 butterfly species reported in the Archipelago, also including specimens from neighbouring areas, and compiled data on 12 species traits and on the apparent extinction of species from the main islands. We calculated indices that measure genetic differentiation, and using phylogenetic regressions we evaluated the relationships between these indices and species traits. Finally, we inferred which traits are associated with disappearance of species on individual islands using phylogenetic regression.ResultsThe overall spatial pattern of genetic diversity corresponded with the proximity of the areas, but strong contrasts were also identified between geographically close areas. Together with the island endemics, several common and widespread species had a high genetic diversification among islands and mainland. Phylogenetic regressions revealed that smaller‐sized, more specialized species, with a preference for drier regions, displayed greater genetic structure and/or haplotype uniqueness. Species that disappeared from islands had a higher population diversification. Capraia has experienced a notable loss of diversity, which significantly affected species with shorter flight periods.Main conclusionsTuscan island butterflies are characterized by strong genetic contrasts and species differ in their contribution to the overall genetic diversity. By ranking the species for their contribution to genetic diversity and identifying the traits linked to the emergence and maintenance of diversity, we have developed a valuable tool for prioritizing populations as targets for monitoring and conservation action. The dataset constructed also represents a valuable resource for testing biogeographical hypotheses.

No abstract

Motivation: Butterflies represent a model in biology and a flagship group for invertebrate conservation. We provide four new resources for the Western Palaearctic butterflies: (1) an updated checklist comprising 552 species; (2) a curated dataset of 32,126 mitochondrial cytochrome c oxidase subunit I (COI) sequences for 532 species, including a de novo reference library for the Maghreb (Morocco, northern Algeria and Tunisia) and Macaronesia (Azores, Madeira and Canary Islands); (3) seven indexes of intraspecific genetic variation (IGV): observed and expected number of haplotypes, haplotype and nucleotide diversity, two fixation indexes and maximum p-distance; and (4) species-level maps illustrating the distribution of COI variability and haplotype networks. The updated checklist will be fundamental for any application dealing with butterfly diversity in the Western Palaearctic.The IGV indexes provide measures for genetic polymorphism and spatial structure and represent proxies for dispersal capacity. These resources will facilitate comparative studies of macrogenetics, foster integrative taxonomy and aid conservation strategies. Main types of variables contained:A complete species checklist in table format, 32,126 mitochondrial DNA barcodes provided with metadata (species membership, WGS84 coordinates and sequence length) and a book in PDF format, including the IGV atlas and indexes, are provided. Spatial location and grain:The checklist encompasses Europe up to the Urals in the east, north Macaronesia (the Azores, Madeira and the Canary Islands) and the Maghreb (Morocco, northern Algeria and Tunisia). COI sequences have been retained in the geographical interval of −31.3 to 67.5° of longitude and 27.5-71.2° of latitude.