DNA barcode reference library for Iberian butterflies enables a continental-scale preview of potential cryptic diversity

Performance of DNA metabarcoding, standard barcoding and morphological approaches in the identification of insect biodiversity

For two decades, DNA barcoding and, more recently, DNA metabarcoding have been used for molecular species identification and estimating biodiversity. Despite their growing use, few studies have systematically evaluated these methods. This study aims to evaluate the efficacy of barcoding methods in identifying species and estimating biodiversity, by assessing their consistency with traditional morphological identification and evaluating how assignment consistency is influenced by taxonomic group, sequence similarity thresholds and geographic distance. We first analysed 951 insect specimens across three taxonomic groups: butterflies, bumblebees and parasitic wasps, using both morphological taxonomy and single-specimen COI DNA barcoding. An additional 25,047 butterfly specimens were identified by COI DNA metabarcoding. Finally, we performed a systematic review of 99 studies to assess average consistency between insect species identity assigned via morphology and COI barcoding and to examine the distribution of research effort. Species assignment consistency was influenced by taxonomic group, sequence similarity thresholds and geographic distance. An average assignment consistency of 49% was found across taxonomic groups, with parasitic wasps displaying lower consistency due to taxonomic impediment. The number of missing matches doubled with a 100% sequence similarity threshold and COI intraspecific variation increased with geographic distance. Metabarcoding results aligned well with morphological biodiversity estimates and a strong positive correlation between sequence reads and species abundance was found. The systematic review revealed an 89% average consistency and also indicated taxonomic and geographic biases in research effort. Together, our findings demonstrate that while problems persist, barcoding approaches offer robust alternatives to traditional taxonomy for biodiversity assessment.

The last days of Aporia crataegi (L.) in Britain: Evaluating genomic erosion in an extirpated butterfly

Current rates of habitat degradation and climate change are causing unprecedented declines in global biodiversity. Studies on vertebrates highlight how conservation genomics can be effective in identifying and managing threatened populations, but it is unclear how vertebrate-derived metrics of genomic erosion translate to invertebrates, with their markedly different population sizes and life histories. The Black-veined White butterfly (Aporia crataegi) was extirpated from Britain in the 1920s. Here, we sequenced historical DNA from 17 specimens collected between 1854 and 1924 to reconstruct demography and compare levels of genomic erosion between extirpated British and extant European mainland populations. We contrast these results using modern samples of the Common Blue butterfly (Polyommatus icarus); a species with relatively stable demographic trends in Great Britain. We provide evidence for bottlenecks in both these species around the period of post-glacial colonization of the British Isles. Our results reveal different demographic histories and Ne for both species, consistent with their fates in Britain, likely driven by differences in life history, ecology and genome size. Despite a difference, by an order of magnitude, in historical effective population sizes (Ne), reduction in genome-wide heterozygosity in A. crataegi was comparable to that in P. icarus. Symptomatic of A. crataegi's disappearance were marked increases in runs-of-homozygosity (RoH), potentially indicative of recent inbreeding, and accumulation of putatively mildly and weakly deleterious variants. Our results provide a rare glimpse of genomic erosion in a regionally extinct insect and support the potential use of genomic erosion metrics in identifying invertebrate populations or species in decline.

Phylogeography of Two Enigmatic Sulphur Butterflies, Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897 (Lepidoptera, Pieridae), with Relations to Wolbachia Infection

The genus Colias Fabricius, 1807 includes numerous taxa and forms with uncertain status and taxonomic position. Among such taxa are Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897, interpreted in the literature either as conspecific forms, as subspecies of different but morphologically somewhat similar Colias species or as distinct species-level taxa. Based on mitochondrial and nuclear DNA markers, we reconstructed a phylogeographic pattern of the taxa in question. We recover and include in our analysis DNA barcodes of the century-old type specimens, the lectotype of C. tamerlana deposited in the Natural History Museum (Museum für Naturkunde), Berlin, Germany (ZMHU) and the paralectotype of C. tamerlana and the lectotype of C. mongola deposited in the Zoological Institute, Russian Academy of Sciences, St. Petersburg, Russia (ZISP). Our analysis grouped all specimens within four (HP_I-HP_IV) deeply divergent but geographically poorly structured clades which did not support nonconspecifity of C. mongola-C. tamerlana. We also show that all studied females of the widely distributed haplogroup HP_II were infected with a single Wolbachia strain belonging to the supergroup B, while the males of this haplogroup, as well as all other investigated specimens of both sexes, were not infected. Our data highlight the relevance of large-scale sampling dataset analysis and the need for testing for Wolbachia infection to avoid erroneous phylogenetic reconstructions and species misidentification.

Integrative analysis reveals the divergence and speciation between sister Sooty Copper butterflies Lycaena bleusei and L. tityrus

The comparison of closely related taxa is cornerstone in biology, as understanding mechanisms leading up to differentiation in relation to extant shared characters are powerful tools in interpreting the evolutionary process. Hotspots of biodiversity such as the west-Mediterranean, where many lineages meet are ideal grounds to study these processes. We set to explore the interesting example of Sooty Copper butterflies: widespread Eurasian Lycaena tityrus (Poda, 1761) comes into contact in Iberia with closely related and local endemic, L. bleusei (Oberthür, 1884), which hasn't always been considered a distinct species. An integrative analysis was designed, combining the use of extensive molecular data (five genes), geometric morphometrics analyses, verified and up-to-date distribution data, and environmental niche modelling, aimed at deciphering their true relationship, their placement within European Lycaena and trace their evolutionary history. We revealed several levels of differentiation: L. bleusei and L. tityrus appear to be reciprocally monophyletic independent gene-pools, distinct in all genes analysed, having mutually diverged 4.8 Ma ago. L. tityrus but not L. bleusei, further displays a genetic structure compatible with several glacial refugia, where populations assignable to infraspecific taxa surface. Conversely, L. bleusei shows a loss in mtDNA diversity in relation to nuDNA. Morphological analyses differentiate both species according to size and shape but also discriminate strong seasonal and sexual traits and a geographical phenotype segregation in L. tityrus. Finally, updated distribution and its modelling for current and glacial timeframes reveal both species respond differently to environmental variables, defining a mostly parapatric distribution and an overlapping belt where sympatry was recovered. During the last glacial maximum, a wider expansion in L. bleusei distribution explains current isolated populations. Our study highlights the importance of gathering several lines of evidence when deciphering the relationships between closely related populations in the fringe of cryptic species realm.

Diversification linked to larval host plant in the butterfly Eumedonia eumedon

It is widely accepted that the relationship between phytophagous insects and their host plants influences insect diversification. However, studies addressed at documenting host-associated genetic differentiation (HAD) and the mechanisms that drive reproductive isolation in host-associated lineages (or host races) are still scarce relative to insect diversity. To uncover further evidence on the HAD processes in Lepidoptera, we investigated the genetic structure of the geranium argus butterfly (Eumedonia eumedon) and tested for isolation by ecology (IBE) vs. isolation by distance (IBD). Genomic data revealed an array of host races (three of them in the same mountain range, the Cantabrian Mountains, northern Iberia) at apparently distinct levels of reproductive isolation. We found a pattern of IBE mediated by HAD at both local and European scales, in which genetic differentiation between populations and individuals correlated significantly with the taxonomic relatedness of the host plants. IBD was significant only when considered at the wider European scale. We hypothesize that, locally, HAD between Geranium-feeding populations was caused (at least partially) by allochrony, that is via adaptation of adult flight time to the flowering period of each host plant species. Nevertheless, the potential reproductive isolation between populations using Erodium and populations using Geranium cannot be explained by allochrony or IBD, and other mechanisms are expected to be at play.

Molecular Phylogeny and Taxonomy of the Butterfly Subtribe Scolitantidina with Special Focus on the Genera Pseudophilotes, Glaucopsyche and Iolana (Lepidoptera, Lycaenidae)

The Palearctic blue butterfly genus Pseudophilotes Beuret, 1958 is not homogenous regarding the morphology of its genital structures. For this reason, some of its species have been considered to be representatives of other genera of the subtribe Scolitantidina (subfamily Polyommatinae). Here, we address these taxonomic problems by analyzing the phylogenetic relationships between the genera, subgenera, and species of this subtribe inferred via the analysis of five nuclear and two mitochondrial DNA sequences. We demonstrate that the enigmatic Asian species P. panope (Eversmann, 1851) belongs to the genus Pseudophilotes but not to Praephilotes Forster, 1938 or Palaeophilotes Forster, 1938 and does not represent the independent genus Inderskia Korshunov, 2000, as hypothesized previously. We synonymize P. svetlana Yakovlev, 2003 (syn. nov.) and P. marina Zhdanko, 2004 (syn. nov.) with P. panope. We demonstrate a deep genetic divergence between lineages that were previously considered as subspecies of the single species Iolana iolas (Ochsenheimer, 1816). As a result, we confirm the multispecies concept of the genus Iolana Bethune-Baker, 1914. We show that the Holarctic genus Glaucopsyche can be divided into four subgenera: Glaucopsyche Scudder, 1872 (=Shijimiaeoides Beuret, 1958), Apelles Hemming, 1931, Bajluana Korshunov and Ivonin, 1990, and Phaedrotes Scudder, 1876.

The genome sequence of the scarce swallowtail, Iphiclides podalirius

The scarce swallowtail, Iphiclides podalirius (Linnaeus, 1758), is a species of butterfly in the family Papilionidae. Here, we present a chromosome-level genome assembly for Iphiclides podalirius as well as gene and transposable element annotations. We investigate how the density of genomic features differs between the 30 Iphiclides podalirius chromosomes. We find that shorter chromosomes have higher heterozygosity at four-fold-degenerate sites and a greater density of transposable elements. While the first result is an expected consequence of differences in recombination rate, the second suggests a counter-intuitive relationship between recombination and transposable element evolution. This high-quality genome assembly, the first for any species in the tribe Leptocircini, will be a valuable resource for population genomics in the genus Iphiclides and comparative genomics more generally.

One Species, Hundreds of Subspecies? New Insight into the Intraspecific Classification of the Old World Swallowtail (Papilio machaon Linnaeus, 1758) Based on Two Mitochondrial DNA Markers

The Old World swallowtail Papilio machaon Linnaeus, 1758 is one of the most well-known and most characteristic members of the family Papilionidae. Over the past two centuries, the butterfly has been the subject of many studies. P. machaon is characterised by a tendency to change the wing colour pattern. In turn, due to the great interest of collectors and amateur entomologists, these studies have been converted into the description of over 100 colour forms, aberrations and subspecies. In this study, mitochondrial DNA (mtDNA), 16S rDNA and cytochrome b sequences were used to examine the correlation between the intraspecific classification and genetic structure of P. machaon. The study used 87 specimens from 59 different localities covering the geographic distribution of this species in the Palaearctic. The phylogenetic relationships within and between the Old World swallowtail subspecies showed that the intraspecific classification proposed by various authors does not correlate with the variability in mitochondrial DNA sequences. In addition, populations occurring at the species distribution borders in the Palaearctic Region (i.e., Japan, Kamchatka, Morocco and Sakhalin) are genetically distinct from other species.

Butterfly–parasitoid–hostplant interactions in Western Palaearctic Hesperiidae: a DNA barcoding reference library

The study of ecological interactions between plants, phytophagous insects and their natural enemies is an essential but challenging component for understanding ecosystem dynamics. Molecular methods such as DNA barcoding can help elucidate these interactions. In this study, we employed DNA barcoding to establish hostplant and parasitoid interactions with hesperiid butterflies, using a complete reference library for Hesperiidae of continental Europe and north-western Africa (53 species, 100% of those recorded) based on 2934 sequences from 38 countries. A total of 233 hostplant and parasitoid interactions are presented, some recovered by DNA barcoding larval remains or parasitoid cocoons. Combining DNA barcode results with other lines of evidence allowed 94% species-level identification for Hesperiidae, but success was lower for parasitoids, in part due to unresolved taxonomy. Potential cases of cryptic diversity, both in Hesperiidae and Microgastrinae, are discussed. We briefly analyse the resulting interaction networks. Future DNA barcoding initiatives in this region should focus attention on north-western Africa and on parasitoids, because in these cases barcode reference libraries and taxonomy are less well developed.

Integrative taxonomy reveals cryptic diversity in North American Lasius ants, and an overlooked introduced species

Biological invasions are a grave threat to ecosystems. The black garden ant (Lasius niger) is a pest species in Europe. Current literature states that L. niger occupies a disjunct native distribution in the Holarctic, however, based on recent work, we re-evaluate this distribution. The native range of L. niger is reconsidered based on phylogenetic relationships (nine mitochondrial and nuclear markers, 5670 bp), DNA-barcoding (98 Holarctic specimens), morphometry (88 Holarctic specimens, 19 different measurements) and subjective assessment of phenotype. The potential spread of this species is estimated using ecological niche modeling. Lasius niger is more closely related to other Palearctic species than to the Nearctic ants known under this name. The latter are described as a distinct species, L. ponderosae sp. nov. However, DNA-barcoding discovered established populations of L. niger in metropolitan areas in Canada (Vancouver and Halifax). We describe a morphometrical method to delineate L. ponderosae sp. nov. and L. niger. MtDNA diversity and divergence is high within L. ponderosae sp. nov., but low within L. niger. More than 1,000,000 km² are suitable as a habitat for L. niger in North America. This case emphasizes the critical role of integrative taxonomy to detect cryptic species and identify potential biological invasions in their nascent stages.

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

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

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

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

Combining DNA metabarcoding and ecological networks to inform conservation biocontrol by small vertebrate predators

In multifunctional landscapes, diverse communities of flying vertebrate predators provide vital services of insect pest control. In such landscapes, conservation biocontrol should benefit service-providing species to enhance the flow, stability and resilience of pest control services supporting the production of food and fiber. However, this would require identifying key service providers, which may be challenging when multiple predators interact with multiple pests. Here we provide a framework to identify the functional role of individual species to pest control in multifunctional landscapes. First, we used DNA metabarcoding to provide detailed data on pest species predation by diverse predator communities. Then, these data were fed into an extensive network analysis, in which information relevant for conservation biocontrol is gained from parameters describing network structure (e.g., modularity) and species roles in such network (e.g., centrality, specialization). We applied our framework to a Mediterranean landscape, where 19 bat species were found to feed on 132 insect pest species. Metabarcoding data revealed potentially important bats that consumed insect pest species in high frequency and/or diversity. Network analysis showed a modular structure, indicating sets of bat species that are required to regulate specific sets of insect pests. A few generalist bats had particularly important roles, either at network or module levels. Extinction simulations highlighted six bats, including species of conservation concern, which were sufficient to ensure that over three-quarters of the pest species had at least one bat predator. Combining DNA metabarcoding and ecological network analysis provides a valuable framework to identify individual species within diverse predator communities that might have a disproportionate contribution to pest control services in multifunctional landscapes. These species can be regarded as candidate targets for conservation biocontrol, although additional information is needed to evaluate their actual effectiveness in pest regulation.

Cytochrome c oxidase subunit I barcode species delineation methods imply critically underestimated diversity in ‘common’ Hermeuptychia butterflies (Lepidoptera: Nymphalidae: Satyrinae)

DNA ‘barcoding’ has contributed greatly to resolving species limits in rapidly diverging tropical insect groups. Here, we explored species diversity in the widespread, cryptic Neotropical butterfly genus Hermeuptychia by generating 601 new sequences of the cytochrome c oxidase I (COI) barcode, tripling available information. Our dataset focused in particular on Ecuador, a biogeographically pivotal country, with additional sequences from Brazil, Costa Rica, French Guiana, Guatemala, Mexico, Panama and Peru. We examined the performance of two phylogeny-based approaches for objectively delineating species, Generalized Mixed Yule Coalescent (GMYC) and Poisson Tree Processes (PTP), as well as a clustering-based approach, Automatic Barcode Gap Discovery (ABGD), on the combined dataset, including our data and 302 published sequences. In general, GMYC and PTP tended to cluster or split likely species as assessed from morphology, depending on model settings, whereas ABGD was less sensitive and produced a more plausible classification. Numerous sequences formed well-supported clades, putative species, that were unrepresented in previously published datasets. The average diversity across all methods was 45 species, in contrast to the 11 species recognized in the current taxonomy. The resulting COI dataset, in combination with ongoing genomic and morphological research, should significantly clarify our understanding of Hermeuptychia species diversity.

Comparative mitogenomics of Drosophilidae and the evolution of the Zygothrica genus group (Diptera, Drosophilidae)

The Zygothrica genus group of Drosophilidae encompasses more than 437 species and five genera. Although knowledge regarding its diversity has increased, uncertainties about its monophyly and position within Drosophilidae remain. Genomic approaches have been widely used to address different phylogenetic questions and analyses involving the mitogenome have revealed a cost-efficient tool to these studies. Thus, this work aims to characterize mitogenomes of three species of the Zygothrica genus group (from the Hirtodrosophila, Paraliodrosophila and Zygothrica genera), while comparing them with orthologous sequences from other 23 Drosophilidae species and addressing their phylogenetic position. General content concerning gene order and overlap, nucleotide composition, start and stop codon, codon usage and tRNA structures were compared, and phylogenetic trees were constructed under different datasets. The complete mitogenomes characterized for H. subflavohalterata affinis H002 and P. antennta present the PanCrustacea gene order with 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, 13 protein coding genes and an A+T rich region with two T-stretched elements. Some peculiarities such as the almost complete overlap of genes tRNAH/ND4, tRNAF/ND5 and tRNAS2/ND1 are reported for different Drosophilidae species. Non-canonical secondary structures were encountered for tRNAS1 and tRNAY, revealing patterns that apply at different phylogenetic scales. According to the best depiction of the mitogenomes evolutionary history, the three Neotropical species of the Zygothrica genus group encompass a monophyletic lineage sister to Zaprionus, composing with this genus a clade that is sister to the Drosophila subgenus.

Phylogeny, species delimitation and biogeography of the endemic Palaearctic tribe Tomarini (Lepidoptera: Lycaenidae)

The tribe Tomarini is represented by the sole genus Tomares, comprising about eight species distributed from the western Mediterranean to Central Asia. We carried out a multilocus phylogenetic and a biogeographical analysis to test the taxonomy of the genus by several molecular species delimitation methods and reveal patterns shaping the current distribution of Tomares. The phylogenetic analysis based on four molecular markers recovered the monophyly of the genus and recovered two deep-branching lineages: an African clade and an Asian clade. Species delimitation analyses suggested six or ten putative species depending on the method applied. The haplotype network analysis of the Tomares nogelii clade revealed no phylogeographical and taxonomic structure. We consider the taxon Tomares nesimachus (syn. nov.) a synonym of T. nogelii and reinstate Tomares callimachus dentata stat. rev. for populations from south-eastern Turkey. Tomares originated between the early Oligocene and the early Miocene, most probably in south-west Asia. The split of the most recent common ancestor of Tomares occurred between the middle-late Miocene and middle-late Pliocene, probably as a response to increasing aridification and habitat fragmentation. Differentiation of the Asian clade took place in south-west Asia during the Pliocene and Pleistocene and coincided temporally with the evolution of Tomares host plants of the genus Astragalus (Fabaceae).

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.

Genetic variation in neotropical butterflies is associated with sampling scale, species distributions, and historical forest dynamics

Previous studies of butterfly diversification in the Neotropics have focused on Amazonia and the tropical Andes, while southern regions of the continent have received little attention. To address the gap in knowledge about the Lepidoptera of temperate South America, we analysed over 3000 specimens representing nearly 500 species from Argentina for a segment of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Representing 42% of the country's butterfly fauna, collections targeted species from the Atlantic and Andean forests, and biodiversity hotspots that were previously connected but are now isolated. We assessed COI effectiveness for species discrimination and identification and how its performance was affected by geographic distances and taxon coverage. COI data also allowed to study patterns of genetic variation across Argentina, particularly between populations in the Atlantic and Andean forests. Our results show that COI discriminates species well, but that identification success is reduced on average by ~20% as spatial and taxonomic coverage rises. We also found that levels of genetic variation are associated with species' spatial distribution type, a pattern which might reflect differences in their dispersal and colonization abilities. In particular, intraspecific distance between populations in the Atlantic and Andean forests was significantly higher in species with disjunct distributions than in those with a continuous range. All splits between lineages in these forests dated to the Pleistocene, but divergence dates varied considerably, suggesting that historical connections between the Atlantic and Andean forests have differentially affected their shared butterfly fauna. Our study supports the fact that large-scale assessments of mitochondrial DNA variation are a powerful tool for evolutionary studies.

DNA barcodes combined with multi-locus data of representative taxa can generate reliable higher-level phylogenies

Taxa are frequently labeled incertae sedis when their placement is debated at ranks above the species level, such as their subgeneric, generic, or subtribal placement. This is a pervasive problem in groups with complex systematics due to difficulties in identifying suitable synapomorphies. In this study, we propose combining DNA barcodes with a multilocus backbone phylogeny in order to assign taxa to genus or other higher-level categories. This sampling strategy generates molecular matrices containing large amounts of missing data that are not distributed randomly: barcodes are sampled for all representatives, and additional markers are sampled only for a small percentage. We investigate the effects of the degree and randomness of missing data on phylogenetic accuracy using simulations for up to 100 markers in 1000-tips trees, as well as a real case: the subtribe Polyommatina (Lepidoptera: Lycaenidae), a large group including numerous species with unresolved taxonomy. Our simulation tests show that when a strategic and representative selection of species for higher-level categories has been made for multigene sequencing (approximately one per simulated genus), the addition of this multigene backbone DNA data for as few as 5–10% of the specimens in the total data set can produce high-quality phylogenies, comparable to those resulting from 100% multigene sampling. In contrast, trees based exclusively on barcodes performed poorly. This approach was applied to a 1365-specimen data set of Polyommatina (including ca. 80% of described species), with nearly 8% of representative species included in the multigene backbone and the remaining 92% included only by mitochondrial COI barcodes, a phylogeny was generated that highlighted potential misplacements, unrecognized major clades, and placement for incertae sedis taxa. We use this information to make systematic rearrangements within Polyommatina, and to describe two new genera. Finally, we propose a systematic workflow to assess higher-level taxonomy in hyperdiverse groups. This research identifies an additional, enhanced value of DNA barcodes for improvements in higher-level systematics using large data sets. [Birabiro; DNA barcoding; incertae sedis; Kipepeo; Lycaenidae; missing data; phylogenomic; phylogeny; Polyommatina; supermatrix; systematics; taxonomy]

The efficacy of DNA barcoding in the classification, genetic differentiation, and biodiversity assessment of benthic macroinvertebrates

Macroinvertebrates have been recognized as key ecological indicators of aquatic environment and are the most commonly used approaches for water quality assessment. However, species identification of macroinvertebrates (especially of aquatic insects) proves to be very difficult due to the lack of taxonomic expertise in some regions and can become time-consuming. In this study, we evaluated the feasibility of DNA barcoding for the classification of benthic macroinvertebrates and investigated the genetic differentiation in seven orders (Insecta: Ephemeroptera, Plecoptera, Trichoptera, Diptera, Hemiptera, Coleoptera, and Odonata) from four large transboundary rivers of northwest China and further explored its potential application to biodiversity assessment. A total of 1,144 COI sequences, belonging to 176 species, 112 genera, and 53 families were obtained and analyzed. The barcoding gap analysis showed that COI gene fragment yielded significant intra- and interspecific divergences and obvious barcoding gaps. NJ phylogenetic trees showed that all species group into monophyletic species clusters whether from the same population or not, except two species (Polypedilum. laetum and Polypedilum. bullum). The distance-based (ABGD) and tree-based (PTP and MPTP) methods were utilized for grouping specimens into Operational Taxonomic Units (OTUs) and delimiting species. The ABGD, PTP, and MPTP analysis were divided into 177 (p = .0599), 197, and 195 OTUs, respectively. The BIN analysis generated 186 different BINs. Overall, our study showed that DNA barcoding offers an effective framework for macroinvertebrate species identification and sheds new light on the biodiversity assessment of local macroinvertebrates. Also, the construction of DNA barcode reference library of benthic macroinvertebrates in Eurasian transboundary rivers provides a solid backup for bioassessment studies of freshwater habitats using modern high-throughput technologies in the near future.

A DNA barcode library for the butterflies of North America

Although the butterflies of North America have received considerable taxonomic attention, overlooked species and instances of hybridization continue to be revealed. The present study assembles a DNA barcode reference library for this fauna to identify groups whose patterns of sequence variation suggest the need for further taxonomic study. Based on 14,626 records from 814 species, DNA barcodes were obtained for 96% of the fauna. The maximum intraspecific distance averaged 1/4 the minimum distance to the nearest neighbor, producing a barcode gap in 76% of the species. Most species (80%) were monophyletic, the others were para- or polyphyletic. Although 15% of currently recognized species shared barcodes, the incidence of such taxa was far higher in regions exposed to Pleistocene glaciations than in those that were ice-free. Nearly 10% of species displayed high intraspecific variation (>2.5%), suggesting the need for further investigation to assess potential cryptic diversity. Aside from aiding the identification of all life stages of North American butterflies, the reference library has provided new perspectives on the incidence of both cryptic and potentially over-split species, setting the stage for future studies that can further explore the evolutionary dynamics of this group.

High resolution DNA barcode library for European butterflies reveals continental patterns of mitochondrial genetic diversity

The study of global biodiversity will greatly benefit from access to comprehensive DNA barcode libraries at continental scale, but such datasets are still very rare. Here, we assemble the first high-resolution reference library for European butterflies that provides 97% taxon coverage (459 species) and 22,306 COI sequences. We estimate that we captured 62% of the total haplotype diversity and show that most species possess a few very common haplotypes and many rare ones. Specimens in the dataset have an average 95.3% probability of being correctly identified. Mitochondrial diversity displayed elevated haplotype richness in southern European refugia, establishing the generality of this key biogeographic pattern for an entire taxonomic group. Fifteen percent of the species are involved in barcode sharing, but two thirds of these cases may reflect the need for further taxonomic research. This dataset provides a unique resource for conservation and for studying evolutionary processes, cryptic species, phylogeography, and ecology.

Factors affecting the efficiency of molecular species delimitation in a species-rich insect family

In the context of global biodiversity loss, molecular species delimitation approaches can be very useful for accelerating species discovery through DNA taxonomy and inventory through DNA metabarcoding. In this study, the effect of some intrinsic factors on the efficiency of various single-marker species delimitation methods (fixed and variable nucleotide distance thresholds, ABGD, ASAP, GMYC, mPTP) was tested on more than 90 empirical data sets, derived from a set of 7,237 COI sequences attributed to 542 leaf beetles species (Coleoptera: Chrysomelidae). The considered factors were: (i) the number of haplotypes per species (as a proxy for genetic diversity), (ii) the geographic distance among conspecific collection localities (as a proxy of sampling width), (iii) the difficulty related to morphological identification of species, and (iv) the taxonomic rank. Distance-based methods, with on average more than 70% of match with morphological identification, outperformed those relying on phylogenetic trees, with less than 59%. A high number of haplotypes per species was found to have a negative effect on delimitation efficiency, whereas large geographic distances within species had a positive effect. All methods delimitations (except for GMYC) were significantly affected by the presence of species that are difficult to be identified, decreasing their efficiency. Finally, the only method influenced by the taxonomic rank of the data set was GMYC, showing lower efficiency in data sets at the genus than at higher levels. The observed biases we highlighted affecting efficiency could be accounted for when developing input data sets for species delimitation analyses to obtain a more reliable representation of biological diversity.

Evaluating DNA Barcoding for Species Identification and Discovery in European Gracillariid Moths

Gracillariidae is the most species-rich leaf-mining moth family with over 2,000 described species worldwide. In Europe, there are 263 valid named species recognized, many of which are difficult to identify using morphology only. Here we explore the use of DNA barcodes as a tool for identification and species discovery in European gracillariids. We present a barcode library including 6,791 COI sequences representing 242 of the 263 (92%) resident species. Our results indicate high congruence between morphology and barcodes with 91.3% (221/242) of European species forming monophyletic clades that can be identified accurately using barcodes alone. The remaining 8.7% represent cases of non-monophyly making their identification uncertain using barcodes. Species discrimination based on the Barcode Index Number system (BIN) was successful for 93% of species with 7% of species sharing BINs. We discovered as many as 21 undescribed candidate species, of which six were confirmed from an integrative approach; the other 15 require additional material and study to confirm preliminary evidence. Most of these new candidate species are found in mountainous regions of Mediterranean countries, the South-Eastern Alps and the Balkans, with nine candidate species found only on islands. In addition, 13 species were classified as deep conspecific lineages, comprising a total of 27 BINs with no intraspecific morphological differences found, and no known ecological differentiation. Double-digest restriction-site associated DNA sequencing (ddRAD) analysis showed strong mitonuclear discrepancy in four out of five species studied. This discordance is not explained by Wolbachia-mediated genetic sweeps. Finally, 26 species were classified as “unassessed species splits” containing 71 BINs and some involving geographical isolation or ecological specialization that will require further study to test whether they represent new cryptic species.

Use of species delimitation approaches to tackle the cryptic diversity of an assemblage of high Andean butterflies (Lepidoptera: Papilionoidea)

Cryptic biological diversity has generated ambiguity in taxonomic and evolutionary studies. Single-locus methods and other approaches for species delimitation are useful for addressing this challenge, enabling the practical processing of large numbers of samples for identification and inventory purposes. This study analyzed an assemblage of high Andean butterflies using DNA barcoding and compared the identifications based on the current morphological taxonomy with three methods of species delimitation (automatic barcode gap discovery, generalized mixed Yule coalescent model, and Poisson tree processes). Sixteen potential cryptic species were recognized using these three methods, representing a net richness increase of 11.3% in the assemblage. A well-studied taxon of the genus Vanessa, which has a wide geographical distribution, appeared with the potential cryptic species that had a higher genetic differentiation at the local level than at the continental level. The analyses were useful for identifying the potential cryptic species in Pedaliodes and Forsterinaria complexes, which also show differentiation along altitudinal and latitudinal gradients. This genetic assessment of an entire assemblage of high Andean butterflies (Papilionoidea) provides baseline information for future research in a region characterized by high rates of endemism and population isolation.

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".

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

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

DNA barcoding and species delimitation of butterflies (Lepidoptera) from Nigeria

Accurate identification of species is a prerequisite for successful biodiversity management and further genetic studies. Species identification techniques often require both morphological diagnostics and molecular tools, such as DNA barcoding, for correct identification. In particular, the use of the subunit I of the mitochondrial cytochrome c oxidase (COI) gene for DNA barcoding has proven useful in species identification for insects. However, to date, no studies have been carried out on the DNA barcoding of Nigerian butterflies. We evaluated the utility of DNA barcoding applied for the first time to 735 butterfly specimens from southern Nigeria. In total, 699 DNA barcodes, resulting in a record of 116 species belonging to 57 genera, were generated. Our study sample comprised 807 DNA barcodes based on sequences generated from our current study and 108 others retrieved from BOLD. Different molecular analyses, including genetic distance-based evaluation (Neighbor-Joining, Maximum Likelihood and Bayesian trees) and species delimitation tests (TaxonDNA, Automated Barcode Gap Discovery, General Mixed Yule-Coalescent, and Bayesian Poisson Tree Processes) were performed to accurately identify and delineate species. The genetic distance-based analyses resulted in 163 well-separated clusters consisting of 147 described and 16 unidentified species. Our findings indicate that about 90.20% of the butterfly species were explicitly discriminated using DNA barcodes. Also, our field collections reported the first country records of ten butterfly species-Acraea serena, Amauris cf. dannfelti, Aterica galena extensa, Axione tjoane rubescens, Charaxes galleyanus, Papilio lormieri lormeri, Pentila alba, Precis actia, Precis tugela, and Tagiades flesus. Further, DNA barcodes revealed a high mitochondrial intraspecific divergence of more than 3% in Bicyclus vulgaris vulgaris and Colotis evagore. Furthermore, our result revealed an overall high haplotype (gene) diversity (0.9764), suggesting that DNA barcoding can provide information at a population level for Nigerian butterflies. The present study confirms the efficiency of DNA barcoding for identifying butterflies from Nigeria. To gain a better understanding of regional variation in DNA barcodes of this biogeographically complex area, future work should expand the DNA barcode reference library to include all butterfly species from Nigeria as well as surrounding countries. Also, further studies, involving relevant genetic and eco-morphological datasets, are required to understand processes governing mitochondrial intraspecific divergences reported in some species complexes.

Rapid colour shift by reproductive character displacement in Cupido butterflies

Reproductive character displacement occurs when competition for successful breeding imposes a divergent selection on the interacting species, causing a divergence of reproductive traits. Here, we show that a disputed butterfly taxon is actually a case of male wing colour shift, apparently produced by reproductive character displacement. Using double digest restriction-site associated DNA sequencing and mitochondrial DNA sequencing we studied four butterfly taxa of the subgenus Cupido (Lepidoptera: Lycaenidae): Cupido minimus and the taxon carswelli, both characterized by brown males and females, plus C. lorquinii and C. osiris, both with blue males and brown females. Unexpectedly, taxa carswelli and C. lorquinii were close to indistinguishable based on our genomic and mitochondrial data, despite displaying strikingly different male coloration. In addition, we report and analysed a brown male within the C. lorquinii range, which demonstrates that the brown morph occurs at very low frequency in C. lorquinii. Such evidence strongly suggests that carswelli is conspecific with C. lorquinii and represents populations with a fixed male brown colour morph. Considering that these brown populations occur in sympatry with or very close to the blue C. osiris, and that the blue C. lorquinii populations never do, we propose that the taxon carswelli could have lost the blue colour due to reproductive character displacement with C. osiris. Since male colour is important for conspecific recognition during courtship, we hypothesize that the observed colour shift may eventually trigger incipient speciation between blue and brown populations. Male colour seems to be an evolutionarily labile character in the Polyommatinae, and the mechanism described here might be at work in the wide diversification of this subfamily of butterflies.

Molecular and morphometric divergence of four species of butterflies (Nymphalidae and Pieridae) from the Western Himalaya, India

Morphometric and molecular divergence among four butterfly species of the families Nymphalidae and Pieridae from the western Himalaya region were investigated using molecular tools, traditional morphometric measures and a truss network system. The considered species were Danaus chrysippus, Vanessa cardui, Pieris brassicae and Pieris canidia. Traditional taxonomy is sometimes unable to discriminate cryptic species or species that have close morphological features. Although taxonomists carefully examine external body features to differentiate the species; however, there is a risk for misidentification during a visual assessment of cryptic species. Therefore, we aimed to use the truss network system of 14 morphological landmarks interconnected to yield 90 variables about molecular taxonomy. Principal component analysis (PCA), discriminant function analysis (DFA) and cluster analysis (CA) were employed to determine morphometric variations. In the traditional analysis, 79, 68, 16 and 5 characters out of 90 were found significant (p < 0.05) for D. chrysippus, V. cardui, P. brassicae and P. canidia, respectively. One to seven principal components were extracted through PCA; they explained 87.5-100% of the total variance in samples. Notably, DFA correctly classified 100% of the original grouped cases and 100% of the cross-validated grouped cases. However, the variations were not the same for the two different methods (truss and traditional) employed for the analysis. We correctly identified all the species; the interspecies sequence divergence was between 0.1034 and 0.1398, and the intra-species sequence divergence range was 0.0001 to 0.0128 using the Cytochrome c oxidase subunit-I (COI) gene. The present study provides useful information about the application and complementary role of traditional with truss morphometric analysis for the precise identification and classification of the selected species.

Continental-scale patterns of hyper-cryptic diversity within the freshwater model taxon Gammarus fossarum (Crustacea, Amphipoda)

Traditional morphological diagnoses of taxonomic status remain widely used while an increasing number of studies show that one morphospecies might hide cryptic diversity, i.e. lineages with unexpectedly high molecular divergence. This hidden diversity can reach even tens of lineages, i.e. hyper cryptic diversity. Even well-studied model-organisms may exhibit overlooked cryptic diversity. Such is the case of the freshwater crustacean amphipod model taxon Gammarus fossarum. It is extensively used in both applied and basic types of research, including biodiversity assessments, ecotoxicology and evolutionary ecology. Based on COI barcodes of 4926 individuals from 498 sampling sites in 19 European countries, the present paper shows (1) hyper cryptic diversity, ranging from 84 to 152 Molecular Operational Taxonomic Units, (2) ancient diversification starting already 26 Mya in the Oligocene, and (3) high level of lineage syntopy. Even if hyper cryptic diversity was already documented in G. fossarum, the present study increases its extent fourfold, providing a first continental-scale insight into its geographical distribution and establishes several diversification hotspots, notably south-eastern and central Europe. The challenges of recording hyper cryptic diversity in the future are also discussed.

Longer is Not Always Better: Optimizing Barcode Length for Large-Scale Species Discovery and Identification

New techniques for the species-level sorting of millions of specimens are needed in order to accelerate species discovery, determine how many species live on earth, and develop efficient biomonitoring techniques. These sorting methods should be reliable, scalable, and cost-effective, as well as being largely insensitive to low-quality genomic DNA, given that this is usually all that can be obtained from museum specimens. Mini-barcodes seem to satisfy these criteria, but it is unclear how well they perform for species-level sorting when compared with full-length barcodes. This is here tested based on 20 empirical data sets covering ca. 30,000 specimens (5500 species) and six clade-specific data sets from GenBank covering ca. 98,000 specimens ($>$20,000 species). All specimens in these data sets had full-length barcodes and had been sorted to species-level based on morphology. Mini-barcodes of different lengths and positions were obtained in silico from full-length barcodes using a sliding window approach (three windows: 100 bp, 200 bp, and 300 bp) and by excising nine mini-barcodes with established primers (length: 94-407 bp). We then tested whether barcode length and/or position reduces species-level congruence between morphospecies and molecular operational taxonomic units (mOTUs) that were obtained using three different species delimitation techniques (Poisson Tree Process, Automatic Barcode Gap Discovery, and Objective Clustering). Surprisingly, we find no significant differences in performance for both species- or specimen-level identification between full-length and mini-barcodes as long as they are of moderate length ($>$200 bp). Only very short mini-barcodes (<200 bp) perform poorly, especially when they are located near the 5$^\prime$ end of the Folmer region. The mean congruence between morphospecies and mOTUs was ca. 75% for barcodes $>$200 bp and the congruent mOTUs contain ca. 75% of all specimens. Most conflict is caused by ca. 10% of the specimens that can be identified and should be targeted for re-examination in order to efficiently resolve conflict. Our study suggests that large-scale species discovery, identification, and metabarcoding can utilize mini-barcodes without any demonstrable loss of information compared to full-length barcodes. [DNA barcoding; metabarcoding; mini-barcodes; species discovery.].

DNA Barcoding and geographical scale effect: The problems of undersampling genetic diversity hotspots

DNA barcoding identification needs a good characterization of intraspecific genetic divergence to establish the limits between species. Yet, the number of barcodes per species is many times low and geographically restricted. A poor coverage of the species distribution range may hamper identification, especially when undersampled areas host genetically distinct lineages. If so, the genetic distance between some query sequences and reference barcodes may exceed the maximum intraspecific threshold for unequivocal species assignation. Taking a group of Quercus herbivores (moths) in Europe as model system, we found that the number of DNA barcodes from southern Europe is proportionally very low in the Barcoding of Life Data Systems. This geographical bias complicates the identification of southern query sequences, due to their high intraspecific genetic distance with respect to barcodes from higher latitudes. Pairwise intraspecific genetic divergence increased along with spatial distance, but was higher when at least one of the sampling sites was in southern Europe. Accordingly, GMYC (General Mixed Yule Coalescent) single-threshold model retrieved clusters constituted exclusively by Iberian haplotypes, some of which could correspond to cryptic species. The number of putative species retrieved was more reliable than that of multiple-threshold GMYC but very similar to results from ABGD and jMOTU. Our results support GMYC as a key resource for species delimitation within poorly inventoried biogeographic regions in Europe, where historical factors (e.g., glaciations) have promoted genetic diversity and singularity. Future European DNA barcoding initiatives should be preferentially performed along latitudinal gradients, with special focus on southern peninsulas.

A complete time-calibrated multi-gene phylogeny of the European butterflies

With the aim of supporting ecological analyses in butterflies, the third most species-rich superfamily of Lepidoptera, this paper presents the first time-calibrated phylogeny of all 496 extant butterfly species in Europe, including 18 very localised endemics for which no public DNA sequences had been available previously. It is based on a concatenated alignment of the mitochondrial gene COI and up to eleven nuclear gene fragments, using Bayesian inferences of phylogeny. To avoid analytical biases that could result from our region-focussed sampling, our European tree was grafted upon a global genus-level backbone butterfly phylogeny for analyses. In addition to a consensus tree, the posterior distribution of trees and the fully concatenated alignment are provided for future analyses. Altogether a complete phylogenetic framework of European butterflies for use by the ecological and evolutionary communities is presented.

DNA barcoding and species delimitation of the Old World tooth-carps, family Aphaniidae Hoedeman, 1949 (Teleostei: Cyprinodontiformes)

The fishes, which have currently named Aphanius Nardo, 1827 are the relict of the ancient ichthyofauna of the Tethys Sea. For a long time since 1827, the genus name has been subjected to revision by several researchers using mainly morphological features. Until recently, no comprehensive single- or multi-locus DNA barcoding study has been conducted on whole members of the family Aphaniidae. In the present study, by applying four conceptually different molecular species delimitation methods, including one distance-based method, one network-based and two topology-based methods, we examined a single-locus DNA barcode library (COI) diversity for the 268 sequences within the family Aphaniidae from the Old World (57 sequences are new in the present study and 211 sequences were retrieved from NCBI database). The molecular analyses revealed a clearer picture of intra-family relationships and allowed us to clarify the generic names, and also describe a new genus for the family Aphaniidae. Results supported distinction of three major clades related to three genera within this family: i) the first clade includes the A. mento group which are placed in a new genus, Paraphanius gen. nov., found in the Orontes (= Asi) and Tigris-Euphrates River drainage, the Levant in coastal waters and the Dead Sea basin, western Jordan, and in southern Turkey in the Mediterranean basins as well as in central Turkey. This clade positioned at the base of the phylogenetic tree, (ii) the second clade contains the A. dispar-like brackish water tooth-carps which are transferred to the genus Aphaniops Hoedeman, 1951 (type species, Lebias dispar), distributed in the coastal waters around the Red Sea and the Persian Gulf basins; and (iii) the third clade, the genus Aphanius Nardo, 1827 (type species Aphanius nanus = A. fasciatus) contains all the inland and inland-related tooth-carps, which are mainly distributed in the inland waters in Turkey and Iran and also in the inland-related drainages around the Mediterranean basin.

DNA barcode library for European Gelechiidae (Lepidoptera) suggests greatly underestimated species diversity

For the first time, a nearly complete barcode library for European Gelechiidae is provided. DNA barcode sequences (COI gene – cytochrome c oxidase 1) from 751 out of 865 nominal species, belonging to 105 genera, were successfully recovered. A total of 741 species represented by specimens with sequences ≥ 500bp and an additional ten species represented by specimens with shorter sequences were used to produce 53 NJ trees. Intraspecific barcode divergence averaged only 0.54% whereas distance to the Nearest-Neighbour species averaged 5.58%. Of these, 710 species possessed unique DNA barcodes, but 31 species could not be reliably discriminated because of barcode sharing or partial barcode overlap. Species discrimination based on the Barcode Index System (BIN) was successful for 668 out of 723 species which clustered from minimum one to maximum 22 unique BINs. Fifty-five species shared a BIN with up to four species and identification from DNA barcode data is uncertain. Finally, 65 clusters with a unique BIN remained unidentified to species level. These putative taxa, as well as 114 nominal species with more than one BIN, suggest the presence of considerable cryptic diversity, cases which should be examined in future revisionary studies.

Testing the Effectiveness of DNA Barcoding for Biodiversity Assessment of Moths from Nigeria

Comprehensive biodiversity assessment of moths in Nigeria rely greatly on accurate species identification. While most of the Nigerian moths are identified effortlessly using their morphological traits, some taxa are morphologically indistinguishable, which makes it difficult for taxon diagnosis. We investigated the efficiency of the DNA barcode, a fragment of the mitochondrial Cytochrome C oxidase subunit I, as a tool for the identification of Nigerian moths. We barcoded 152 individuals comprising 18 morphospecies collected from one of the remaining and threatened rainforest blocks of Nigeria – the Cross River National Park. Phenetic neighbor-joining tree and phylogenetic Maximum Likelihood approach were employed for the molecular-based species identification. Results showed that DNA barcodes enabled species-level identification of most of the individuals collected from the Park. Additionally, DNA barcoding unraveled the presence of at least six potential new and yet undescribed species—Amnemopsyche sp., Arctia sp., Deinypena sp., Hodebertia sp., Otroeda sp., and Palpita sp. The phylogenetic Maximum Likelihood using the combined dataset of all the newly assembled sequences from Nigeria showed that all species formed unique clades. The phylogenetic analyses provided evidence of population divergence in Euchromia lethe, Nyctemera leuconoe, and Deinypena lacista. This study thus illustrates the efficacy of DNA barcoding for species identification and discovery of potential new species, which demonstrates its relevance in biodiversity documentation of Nigerian moths. Future work should, therefore, extend to the creation of an exhaustive DNA barcode reference library comprising all species of moths from Nigeria to have a comprehensive insight on the diversity of moths in the country. Finally, we propose integrated taxonomic methods that would combine morphological, ecological, and molecular data in the identification and diversity studies of moths in Nigeria.

Twenty-six new species of Hoploscopa (Lepidoptera, Crambidae) from South-East Asia revealed by morphology and DNA barcoding

Hoploscopa Meyrick (Lepidoptera: Crambidae) is a fern-feeding genus found in montane areas of South-East Asia and Melanesia, eastwards up to the Samoan Islands. It includes sixteen described species, with at least 70 further undescribed species known from scientific collections. An iterative approach including morphological and molecular characters was used in order to explore the diversity of Hoploscopa. The hitherto described species are revised, and descriptions authored by T. Léger and M. Nuss are provided for an additional 26 new species: H.agtuuganonensissp. nov., H.albipunctasp. nov., H.albomaculatasp. nov., H.anacanthasp. nov., H.boletasp. nov., H.cynodontasp. nov., H.danaoensissp. nov., H.gombongisp. nov., H.gracilissp. nov., H.ignitamaculaesp. nov., H.isarogensissp. nov., H.jubatasp. nov., H.kelamasp. nov., H.kinabaluensissp. nov., H.mallyisp. nov., H.marijoweissaesp. nov., H.matheaesp. nov., H.niveofasciasp. nov., H.pangrangoensissp. nov., H.parvimaculasp. nov., H.pseudometacrossasp. nov., H.sepanggisp. nov., H.sumatrensissp. nov., H.titikasp. nov., H.tonsepisp. nov., H.ypsilonsp. nov. Using a protocol specific for the amplification of DNA from old museum specimens, we recovered 101 COI barcodes for all but one of the newly described species, with 76 being barcode compliant (>487 bp). Species delimitation analyses suggest cryptic diversity, with six cases reflecting allopatric divergence, and two further cases found in sympatry.

The conundrum of species delimitation: a genomic perspective on a mitogenetically super-variable butterfly

The Palaearctic butterfly Melitaea didyma stands out as one of the most striking cases of intraspecific genetic differentiation detected in Lepidoptera: 11 partially sympatric mitochondrial lineages have been reported, displaying levels of divergence of up to 7.4%. To better understand the evolutionary processes underlying the diversity observed in mtDNA, we compared mtDNA and genome-wide SNP data using double-digest restriction site-associated DNA sequencing (ddRADseq) results from 93 specimens of M. didyma ranging from Morocco to eastern Kazakhstan. We found that, between ddRADseq and mtDNA results, there is a match only in populations that probably remained allopatric for long periods of time. Other mtDNA lineages may have resulted from introgression events and were probably affected by Wolbachia infection. The five main ddRADseq clades supported by STRUCTURE were parapatric or allopatric and showed high pairwise F_ST values, but some were also estimated to display various levels of gene flow. Melitaea didyma represents one of the first cases of deep mtDNA splits among European butterflies assessed by a genome-wide DNA analysis and reveals that the interpretation of patterns remains challenging even when a high amount of genomic data is available. These findings actualize the ongoing debate of species delimitation in allopatry, an issue probably of relevance to a significant proportion of global biodiversity.

The determinants of genetic diversity in butterflies

Under the neutral theory, genetic diversity is expected to increase with population size. While comparative analyses have consistently failed to find strong relationships between census population size and genetic diversity, a recent study across animals identified a strong correlation between propagule size and genetic diversity, suggesting that r-strategists that produce many small offspring, have greater long-term population sizes. Here we compare genome-wide genetic diversity across 38 species of European butterflies (Papilionoidea), a group that shows little variation in reproductive strategy. We show that genetic diversity across butterflies varies over an order of magnitude and that this variation cannot be explained by differences in current abundance, propagule size, host or geographic range. Instead, neutral genetic diversity is negatively correlated with body size and positively with the length of the genetic map. This suggests that genetic diversity is determined both by differences in long-term population size and the effect of selection on linked sites.

A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris

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.

Genetic Variation and Differentiation of Hylesia metabus (Lepidoptera: Saturniidae): Moths of Public Health Importance in French Guiana and in Venezuela

Hylesia moths impact human health in South America, inducing epidemic outbreaks of lepidopterism, a puriginous dermatitis caused by the urticating properties of females' abdominal setae. The classification of the Hylesia genus is complex, owing to its high diversity in Amazonia, high intraspecific morphological variance, and lack of interspecific diagnostic traits which may hide cryptic species. Outbreaks of Hylesia metabus have been considered responsible for the intense outbreaks of lepidopterism in Venezuela and French Guiana since the C20, however, little is known about genetic variability throughout the species range, which is instrumental for establishing control strategies on H. metabus. Seven microsatellites and mitochondrial gene markers were analyzed from Hylesia moths collected from two major lepidopterism outbreak South American regions. The mitochondrial gene sequences contained significant genetic variation, revealing a single, widespread, polymorphic species with distinct clusters, possibly corresponding to populations evolving in isolation. The microsatellite markers validated the mitochondrial results, and suggest the presence of three populations: one in Venezuela, and two in French Guiana. All moths sampled during outbreak events in French Guiana were assigned to a single coastal population. The causes and implications of this finding require further research.

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

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