Biodiversity Data Journal 11: e96438 OO) doi: 10.3897/BDJ.11.e96438 open access Research Article Molecular Weevil Identification Project: A thoroughly curated barcode release of 1300 Western Palearctic weevil species (Coleoptera, Curculionoidea) André Schitte?, Peter E. Stiben§, Jonas J. Astrin? + Leibniz Institute for the Analysis of Biodiversity Change, Museum Koenig, Bonn, Germany § Curculio Institute, M6nchengladbach, Germany Corresponding author: Academic editor: Paulo Borges Received: 17 Oct 2022 | Accepted: 08 Dec 2022 | Published: 24 Jan 2023 Citation: Schitte A, StUben PE, Astrin JJ (2023) Molecular Weevil Identification Project: A thoroughly curated barcode release of 1300 Western Palearctic weevil species (Coleoptera, Curculionoidea). Biodiversity Data Journal 11: e96438. https://doi.org/10.3897/BDJ.11.e96438 Abstract The Molecular Weevil Identification project (MWI) studies the systematics of Western Palearctic weevils (Superfamily Curculionoidea) in an integrative taxonomic approach of DNA barcoding, morphology and ecology. This barcode release provides almost 3600 curated CO1 sequences linked to morphological vouchers in about 1300 weevil species. The dataset is presented in statistical distance tables and as a Neighbour-Joining tree. Bayesian Inference trees are computed for the subfamilies Cryptorhynchinae, Apioninae and Ceutorhynchinae. Altogether, 18 unresolved taxonomic issues are discussed. A new barcode primer set is presented. Finally, we establish group-specific genetic distances for many weevil genera to serve as a tool in species delineation. These values are statistically based on distances between "good species" and their congeners. With this morphologically calibrated approach, we could resolve most alpha-taxonomic questions within the MWI project. © Schutte A et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 2 Schutte A et al Keywords DNA barcoding, integrative taxonomy, thresholds, Cryptorhynchinae, Curculionidae, Apionidae, Western Palearctic, Europe, Canary Islands Introduction With 400,000 described species, beetles (Coleoptera) constitute the most diverse animal order (Slipinski et al. 2011, Bouchard et al. 2017). Amongst them, weevils (Superfamily Curculionoidea) form one of the most species-rich taxa, with 51,000 known species worldwide (Oberprieler et al. 2007). Exactly 15,407 weevil species are listed in the most recent catalogue covering the entire Palearctic realm (Alonso-Zarazaga et al. 2017) and about 3,500 species in the Western Palearctic (L6bl and Smetana 2011, L6bl and Smetana 2013). Weevils have a global distribution. Their larvae predominantly develop inside various plant parts, while adults mostly feed on leaves or roots. Many species are highly specialised; others feed on a wide range of plants (Zwolfer and Herbst 1988, Oberprieler et al. 2007, Letsch et al. 2018). Weevils play an important ecological role. Some species are pests in agriculture or forestry, for example, the large pine weevil Hylobius abietis (Leather et al. 1999), the rice weevil Sitophilus oryzae or the maize weevil Sitophilus zeamais (Wu and Yan 2018). Many Otiorhynchus species are greenhouse or horticultural pests (Thiem 1932, Sprick 2009). Taxonomic identification is easy for some common weevil species, but for many others, it is challenging and requires genital preparation and considerable taxonomic expertise. The taxonomic impediment (Godfray 2002, Godfray and Knapp 2004, Wheeler et al. 2004) implies that the number of experts able to identify organisms to species level is constantly decreasing (Irfanullah 2006), resulting - amongst other drawbacks - in inaccurate biodiversity assessments (Giangrande 2003). To compensate for this deficiency in times of the biodiversity crisis, two DNA-based approaches were simultaneously proposed (reviewed in Meier et al. (2006), Hansen et al. (2007), Teletchea (2010)): The DNA taxonomy concept by Tautz (Tautz et al. 2002, Tautz et al. 2003) proposed to utilise DNA sequences of several predefined standard genes as the scaffold for taxonomy, but not necessarily linked to the Linnaean binominal system. Hebert (Hebert et al. 2003a, Hebert et al. 2003b) envisioned relatively short DNA barcodes from a single gene as a universal system for re-identification purposes, ideally linked to current Linnaean names. Since then, DNA barcoding (Sequencing the 5'-half of the COI gene, Folmer et al. 1994, for animals) has been widely adopted by the scientific community and is the most commonly-used molecular marker in animals (Waugh 2007, Blaxter 2016), also demonstrated by over 11,000 barcoding-related publications by November 2022 (VWOS, Web of Science search for the term "DNA barcod*" in title/abstract/keywords). The Barcode of Life Database (BOLD, Ratnasingham and Hebert 2007), currently contains barcodes from almost 10 million specimens. Genetic distances can be measured as a proportion of different nucleotide positions in percent. An important prerequisite for DNA barcoding is thatinterspecific genetic distances Molecular Weevil Identification Project: A thoroughly curated barcode release ... 3 vary significantly for at least a large majority of cases from intraspecific genetic distances. This concept is often referred to as the barcoding gap (Meyer and Paulay 2005). Its existence was first stated for birds and various arthropod taxa (Hebert et al. 2003a, Hebert et al. 2003b, Hebert et al. 2004a, Hebert et al. 2004b, Barrett and Hebert 2005, Hajibabaei et al. 2006a, Hajibabaei et al. 2006b). The existence of the barcoding gap was thought to be an artefact of insufficient sampling across taxa (Meyer and Paulay 2005, Wiemers and Fiedler 2007, Bergsten et al. 2012). Over the years, it became apparent that some datasets showed more or less pronounced barcoding gaps, while others did not (depending on sampling, geographic region, taxon biology, degree of morphological crypsis, state of taxonomic revision of the group under study etc.). Underlying morphological misidentifications pose a major problem to DNA barcoding datasets and reference collections. Unfortunately, specimen misidentification is common in literature, in collections and particularly widespread in public sequence databases (Pentinsaari et al. 2020) and might reach up to 56% for taxa difficult to identify (Shea et al. 2011). The use of obsolete taxonomic names can lead to similar problematic effects (Mulcahy et al. 2022), especially in the absence of material vouchers (Pleijel et al. 2008, Astrin et al. 2013). The Molecular Weevil Identification project (MW) presented here strives to avoid pitfalls that arise in DNA barcoding studies, when not backed up by an extensive voucher collection. MWI created a reference database of high-quality DNA barcodes from scratch. Almost 1300 Western Palearctic weevil species have been barcoded, based on rigorous vouchering routines and project criteria: DNA was extracted non-invasively from specimens, then mounted as morphological vouchers for the dry collection, accessible at a public natural history collection (Leibniz Institute for the Analysis of Biodiversity Change, Museum Koenig, Bonn, Germany). These morphological vouchers are accompanied by stored DNA extracts and tissue samples in a dedicated biobank at the same institute. The laboratory infrastructure used in MWI was that of the German Barcode of Life (GBOL) project (Geiger et al. 2016). In several taxa, type localities were revisited to sample for the MWI project. Only experienced researchers from the European coleopterists' association Curculio Institute collected and identified the specimens morphologically (see specimen data table in Suppl. material 1). In this barcode release, we did not add any publicly available sequences from GenBank, BOLD or other third parties to the dataset for quality control reasons. For new species described during MWI, the DNA barcodes were mostly generated from paratypes (collected at the same location as the holotype). If no paratypes were available, the holotype was used for barcode generation. Specimens were recovered during the lysis step to allow future validation of the initial identification. Embedding the DNA barcoding method within an integrative taxonomic approach (VVill et al. 2005, Padial and De La Riva 2010) often helps to reveal cryptic diversity or synonyms (Hebert et al. 2004b, Pons et al. 2006, Kerr et al. 2007). Over several years of preparation for the present barcode release, dozens of alpha-taxonomic changes have been carried out, most of which began as conspicuous molecular findings and were then corroborated morphologically (often including the study of type material) and ecologically in a taxonomic feedback loop (Page et al. 2005). Until now, 157 taxonomic changes, including 80 new 4 Schutte A et al species descriptions, were based on MWI sequences. Most taxonomic changes refer to the subfamilies Cryptorhynchinae, Apioninae and Ceutorhynchinae, for example, Schutte et al. (2013), Stuben et al. (2013a), Stuben and Schutte (2013), Stuben (2014b), Stuben et al. (2015), Stuben and Bayer (2015), Stuben et al. (2016a), Stuben (2017b), Stuben (2018b), StUben (2018c), Stuben (2018a), Stuben and Schutte (2018). In practice, most of these initial conspicuous molecular findings consisted of simple genetic distances that were considerably higher (in relation to other intraspecific comparisons) or lower (in relation to other interspecific comparisons) than expected. Previous research shows that genetic distance values very often coincide with species limits, but vary widely by taxon and geographic setting, for example, 2.7% for a species delineation threshold for North American birds (Hebert et al. 2004a), 4% for North American spiders (Barrett and Hebert 2005) or 2% to 14% for Madagascan water beetles (Monaghan et al. 2005). The relevant question in this context is: How does one know which values to expect? Considerable discussion has gone into the topic of genetic thresholds as criteria for species delimitation (10x rule in Hebert et al. (2004a), Monaghan et al. (2005)). We agree that such threshold values cannot be used to delineate species (Wiemers and Fiedler 2007 , Hubert et al. 2010) as a subset under the argument that DNA barcodes alone are generally a poor criterion to describe species (Ahrens et al. 2021, Zamani et al. 2022). However, based on our experience, genetic thresholds can be profitably used as a heuristic criterion to highlight cryptic or problematic taxa in the vast majority of cases (new species, synonyms, species complexes), to be confirmed or refuted by morphological analysis. This study presents a validated, taxonomically thoroughly curated barcode release with almost 3600 sequences, the most extensive Western Palearctic weevil barcode dataset until now, covering ca. 1300 weevil species. Based on this data and considering different distribution patterns, we also give group-specific, statistically derived heuristic hints regarding which genetic distance values fall within the typical interspecific range. Cryptorhynchinae, Ceutorhynchinae and Apioninae are represented in our dataset with a high species coverage for the sampled region. Therefore, we specify such values only for genera of those three subfamilies, providing minimum and average p-distance values. By sharing these data, we hope to accelerate specimen re-identification within the discussed taxa and aid in prefiltering future cases for thorough integrative alpha-taxonomic investigation. Material and Methods We analysed 3573 mitochondrial CO1 sequences of the DNA barcoding region (Hebert et al. 2003a, based on Folmer et al. (1994)), representing 1296 valid species (1391 taxa if infraspecific epithets and taxon qualifiers such as "cf." and "sp." are also counted). The dataset contains 2017 sequences newly released in this study, see Suppl. material 1 (GenBank accession numbers followed by a "new" tag). See Suppl. material 2 for the Molecular Weevil Identification Project: A thoroughly curated barcode release ... 5 entire dataset in FASTA format. Both supplements can also be downloaded externally (DOI: 10.5281/zenodo.7430106). Sampling: collecting locations The geographic origin of the collected weevils is as follows: 2510 specimens (70% of the dataset) were collected throughout continental Europe including the Mediterranean islands; 889 specimens (25% of the dataset) were collected on the Macaronesian islands including the Canaries, Azores, Madeira Archipelago with Desertas Islands and Savage Islands (IIhas Selvagens); 164 specimens (5% of the dataset) were collected in continental North Africa, mostly Morocco and Tunisia. Collecting locations of the sampled specimens are plotted on two ArcGIS map baselayers with GPS Visualizer, see Fig. 1 (continental) and Fig. 2 (Atlantic islands). ~ - K e a? sigs ee ee agi q- coe a ee ~ . as? ce eS * . ee? ee im: F. cote 4 a “ * ed NS . . = ‘ : 2 rae ” uy hd . e 2a 3 ey wf "4 J a 2 pri > te “ . o* ee* . 7 tee? : % eo « ’ pre aes Pe ‘+ . 200k} " Figure 1. EES] Continental collecting locations. Black dots = previously-released MWI sequences, red dots = newly published with this study. The most frequently collected Curculionidae subfamilies were: Cryptorhynchinae (1190 sequences, 278 species, subspecies not differentiated in the species count), Entiminae (576 sequences, 269 species), Ceutorhynchinae (537 sequences, 203 species), Curculioninae (356 sequences, 168 species) and Apioninae (349 sequences, 115 species). A complete overview of sequences per subfamily is shown in Fig. 3. The number of specimens collected per species is illustrated in Fig. 4: 16% of sequences are singletons, the median is six specimens per species. Azores (northern) La Palma ode ‘“y & f# La Gomera ae Poed El Hierro Canary Islands Figure 2. EES] Collecting locations of Schitte A et al Maderia & Porto Santo Desertas J Azores (southern) Lanzarote Selvagens Tenerife ne aove*e +4 Fuerteventura Gran Canaria Atlantic islands. Black dots = previously-released MW sequences, red dots = newly published with this study. Scolytinae (Fam. Curculionidae) Nemonychinae (Fam. Nemonychidae) Brachycerinae (Fam. Brachyceridae) Otiorrhynchinae (Fam. Curculionidae) omer ta ae | Cimberidinae (Fam. Nemonychidae Apoderinae (Fam. Atte! wtadee\ Urodontinae (Fam. Anthribidae Urodontinae (Fam. Urodontidae’ Raymondionyminae (Fam. Raymondionym. Attelabinae (Fam. Attelabidae Conoderinae (Fam. Curculionidae Cyclominae (Fam. Curculionidae Anthribinae (Fam. Anthribidae Rhynchophorinae (Fam. Drysohthorides| Oxycoryninae (Fam. Belidae) Nanophyinae (Fam. Nanophyidae) Phytonominae (Fam. ‘Nonopide} Corimaliinae (Fam. Nanophyidae Mesoptiliinae (Fam. Curculionidae Erirhininae (Fam. Erirhinidae) Molytinae (Fam. Curculionidae) Rhynchitinae (Fam. Rhynchitidae) Bagoinae (Fam. Curculionidae) Hyperinae (Fam. Curculionidae) Cossoninae (Fam. Curculionidae) Lixinae (Fam. Curculionidae) Apioninae (Fam. Apionidae Curculioninae (Fam. Curculionidae Ceutorhynchinae (Fam. Curculionidae) Entiminae (Fam. Curculionidae) Cryptorhynchinae (Fam. Curculionidae) subfamily (family) Figure 3. EES] Hi NNNN ll > Lal ~~ | = ND aa w & > ws ~ br me are ~ a we 8 “ 10 specimens count per family 8 8 10000 Number of specimens per subfamily. Molecular Weevil Identification Project: A thoroughly curated barcode release ... 7 560 272 200 171 98 100 Tr 28 19 = 43 en Se ee ee ee, eS oe: a Bf ° Bite 20a ta oes ee ee ‘te Bi eae - ee ig 8 9 10 11 12 13 14 1S 16 17 18 19 20 23 number of species with x specimens number of specimen Figure 4. EES] Number of specimens sampled per species (x-axis). Number above the bars show the specimen count (y-axis). Subspecies level not considered, unclear identifications are not counted in this figure (58x cf., 18x sp.). For 560 species, only one specimen was sampled. Sampling: methodology Most specimens were collected directly into 96% non-denatured ethanol without killing agents. In some cases, we sequenced previously-collected dried specimens, usually not more than five years old. Field data, voucher numbers and GenBank accession numbers for all specimens are provided in Suppl. material 1. The pinned specimen vouchers (dry), DNA vouchers and, where available, tissue vouchers (frozen, same population as the sequenced individual) are deposited at the Coleoptera collection respectively at the Biobank of the Leibniz Institute for the Analysis of Biodiversity Change, Museum Koenig, Bonn, Germany (ZFMK). For most specimens, the pinned specimen voucher was the source of the non-destructively isolated DNA. In previous extractions, the DNA derived from the frozen tissue samples with equal emphasis on morphological integrity (MWI samples with specimen ID lower than 1823-PST; see Suppl. material 1 for naming scheme). Laboratory processing The laboratory routine for Cryptorhynchinae is described in Astrin et al. (2012). The laboratory routine for 91 sequences (specimen id 2906-PST to 3023-PST) is described in Stuben and Kramp (2019). The laboratory routine for all other samples is as follows: Genomic DNA was extracted from different parts of the beetle or non-destructively (for sample ID 1823-PST and higher) from whole specimens. Tissue lysis was performed at 56° Celsius overnight. For DNA extraction, a BioSprint 96 magnetic bead extractor was used with the corresponding kits, following the manufacturer's protocol for 200 ul elution volume 8 Schutte A et al (Qiagen: Hilden, Germany). We amplified the 5'-end of the CO1 (Cytochrome c oxidase subunit 1) gene with degenerate primers (Table 1): reaction volume of 20 ul; 2.5 ul of undiluted DNA template; Multiplex PCR Master Mix (Qiagen) with 16 pmol primer concentration for each primer (1.6 pul of 10 pmol/ul primer). The standard PCR product retrieved for weevils is 658 nucleotides (nt/bp) in length. Table 1. PCR primer sets. LCO1490-JJ & HCO2198-JJ (Astrin and Stuben 2008) were used for all samples first, success rate was about 95% in both directions. If PCR or sequencing failed, reactions were repeated with a different primer set, either LCO1490-JJ2 & HCO2198-JJ2 (higher degeneracy, Astrin et al. (2016)) or repeated with the newly-developed primer set LCO1490-MWI & HCO2198- MWI. Primer Name _ 5'-3' Read Direction Reference Optimised for LCO1490-MWI_ ACWAAYCATAARRAYATYGG this study (new) Apioninae & Ceutorhynchinae HCO2198- | TADACTTCDGGRTGDCCRAARAATCA this study (new) to tla MWI LCO1490-JJ. | CHACWAAYCATAAAGATATYGG Astrin and Stiiben Cryptorhynchinae (2008) HCO2198-JJ | AWACTTCVGGRTGVCCAAARAATCA _ Astrin and Stiiben (2008) LCO1490-JJ2 CHACWAAYCAYAARGAYATY GG Astrin et al. (2016) universal (arthropods) HCO2198-JJ2_ANACTTCNGGRTGNCCAAARAATCA _ Asirin et al. (2016) Thermal cycling was performed on GeneAmp PCR System 2700 instruments (Life Technologies, Carlsbad, USA) as follows: hot start Taq activation: 15 min at 95°C; first cycle set ("touch down" with 15 repeats): 35 s denaturation at 94°C, 90 s annealing at 55°C (-1°C/cycle) and 90 s extension at 72°C. Second cycle set (25 repeats): 35 s denaturation at 94°C, 90 s annealing at 40°C and 90 s extension at 72°C; final elongation: 10 min at 72°C. Amplicons were purified with the ExoSAP-IT kit (USB Corporation, Cleveland, Ohio) and sequenced bidirectionnally using the PCR primers (Table 1) at BGI Genomics (Shenzhen, China) or Macrogen (Amsterdam, The Netherlands) facility. Data analyses Contig assembly and trimming of primer regions were performed in Geneious Pro 6.1.8 (Kearse et al. 2012). The sequences were screened for: 1) pseudogenes (NUMTs, Song et al. (2008)) by inspection of the reading frame for stop-codons and 2) endosymbionts by visual inspection of each taxon position in the NJ tree and NCBI BLAST (Madden 2002). In case of inconsistency, the sequence was excluded, followed by re-amplification or re- extraction and re-amplification. Molecular Weevil Identification Project: A thoroughly curated barcode release ... 9 The dataset contains 3573 weevil barcodes, of which 3302 sequences cover the full barcode length (658 bp), 272 sequences are shorter. Two sequences (GU987885, MG229813) barely failed to reach the 500 bp minimum required by BOLD (Milton et al. 2013) and were kept in the dataset. Alignment. DNA sequences were aligned with the Muscle (Edgar 2004) plug-in in Geneious using default parameters (Drummond et al. 2012) and visually inspected for misaligned ends. The alignment is provided in Suppl. material 2. Neighbour-Joining tree. The Neighbour-Joining (NJ, Saitou and Nei (1987)) tree is based on the nucleotide sequence alignment of the entire dataset of 3573 weevil sequences plus One outgroup species (Chrysomelidae, GenBank FJ867810). The NJ tree was created in Geneious Pro 6.1.8 (Kearse et al. 2012) with JC69 nucleotide substitution model (Jukes and Cantor 1969). The tree is provided in Suppl. material 3. Bayesian Inference. Phylogenetic trees, based on Bayesian Inference, were reconstructed for three sub-datasets: 1) MrBayes sub-dataset for Cryptorhynchinae + Cossoninae: 1311 sequences in total, 1190 sequences from Cryptorhynchinae, 120 additional sequences from Cossoninae plus one outgroup species (Anthribidae, GenBank FJ867818). 2) MrBayes sub-dataset for Apioninae + Nanophyinae + Attelabidae: 367 sequences in total, 349 sequences from Apioninae, 5 additional from Attelabidae, 12 additional from Nanophyinae plus one outgroup species (Cryptorhynchus lapathi, D-0354-lap, GenBank EU286523). 3) MrBayes sub-dataset for Ceutorhynchinae: 537 Ceutorhynchinae sequences plus one outgroup (Cryptorhynchus lapathi, D-0354-lap, GenBank EU286523). Based on the Bayesian information criterion value (BIC, Schwarz (1978)), calculated with jModelTest 0.1.1 (Posada 2008), we applied the GTR+I+G substitution model (Lanave et al. 1984) for all Bayesian analyses. We ran MrBayes (Ronquist and Huelsenbeck 2003) MPI version 3.2.7 multiprocessor version with eight cores in two independent replicates, each with one cold chain and three chains of different temperatures (standard setting). The genetic code for metazoan mitochondrial DNA (metmt) was defined. The third codon position of the GO1 gene was unlinked in shape, revmat, statefreq and pinvar. The analyses ran for 20 million generations, sampling 20,000 trees. Negative log-likelihood score stabilisation was checked in a separate visualisation in Microsoft Excel 2013. Accordingly, we retained 19,900 trees after discarding the burn-in data, of which a 50%- majority rule consensus tree with posterior probabilities was built. Geneious was used for the graphical display of the tree. The trees are provided in Suppl. material 5. DiStats statistics (p-distance calculation) The Perl script DiStats (Astrin et al. 2016) simplifies the processing and statistical inspection of DNA barcode datasets. Amongst other functions, it calculates intraspecific 10 Schtitte A et al and interspecific genetic distances for a given nucleotide alignment. The interspecific p- distance values per genus and distribution (island, continental) are provided for the genera of the three weevil subfamilies in the focus of this study. See Suppl. material 6 for an in- depth description of input data selection, raw output files, DiStats results and data compilation. The following references are used exclusively in the suppl. material: Stuben and Behne (2010), Morris (2011), Kratky (2015), Morris and Barclay (2015), Russell and Velazquez de Castro (2015), Stuben (2017a) and Sprick (2019). The description below is the shortened version. Confidence groups. For DiStats analysis, only Cryptorhynchinae, Ceutorhynchinae and Apioninae species are taken into account (datasets of the best-sampled subfamilies). The sequences of each species are assigned to one of three confidence groups: ° Confidence group 1 (reference species / "good species"): Contains taxa, which were morphologically clear in the past; one synonym allowed for Cryptorhynchinae, three synonyms allowed for Ceutorhynchinae and Apioninae, otherwise moved to confidence group 2; ° confidence group 2 (congener dataset): Contains valid taxa which created some or many synonyms or subspecies, not evaluated as reference species (not a "good species"), but available as congeners in the dataset; taxa, which were morphologically difficult or ambiguous to identify; ° confidence group 3: (omitted species or specimens): Contains problematic taxa like species complexes or potentially new species, those were excluded from DiStats analyses. Only sequences / specimens from confidence groups 1 and 2 were used in DiStats statistics. Distribution groups. The reference species ("good species)" were also assigned to one out of four geographical distribution groups (Table 2) to assess the effect of different geographical distribution sizes with respect to interspecific genetic distances. The external supplement (DOI: 10.5281/zenodo.7430565) contains the geographical distribution maps used for the estimation of the maximum distribution range of each reference species. Table 2. Geographical distribution groups defined for the genera lists. Distribution group ISL (island) C1 (endemic) C2 (medium) C3 (large) Cryptorhynchinae island(s) up to 50 km 50 to 500 km 500 km and above Apioninae island(s) up to 50 km 50 to 2,000 km 2,000 km and above Ceutorhynchinae island(s) up to 50 km 50 to 2,000 km 2,000 km and above Interspecific distances per genus and distribution. VWWe examined the p-distances from the reference species to its closest congeners. Only for the reference species, the distance values to each closest congener were used to create genus lists with minimum and Molecular Weevil Identification Project: A thoroughly curated barcode release ... 11 average interspecific distance values per geographic distribution group. The closest congener can be another reference species or a taxon assigned to confidence group 2 (congener dataset). We never used the p-distances from taxa of confidence group 2; those were kept only to increase the amount of congeners in the DiStats dataset. ASAP analysis The programme ‘Assemble Species by Automatic Partitioning’ (ASAP, Puillandre et al. (2021)) estimates the species amount in a barcode dataset and suggests 10 p-distance threshold values for species delineation. ASAP is the successor of the programme ‘Automated Barcode Gap Discovery' (ABGD, Puillandre et al. (2012)). We used the three sub-datasets created for DiStats with ASAP (Cryptorhynchinae, Apioninae, Ceutorhynchinae). Those sub-datasets only contain taxa from confidence group 1 ("good species") and confidence group 2 ("congener dataset"), see DiStats above. We compare the ASAP-calculated genospecies (MOTUs, Floyd et al. (2002), Blaxter (2004)) with the morphological identification and counted wrongly assigned species for each dataset. We have two questions: 1) how reliable is a species delineation, based on a single threshold per dataset and 2) is the best fitting threshold suggested by the programme the best possible one to match the morphological species identifications? See Suppl. material 7 for further details about the data assembly. Results There remain 18 apparent contradictions between morphological identification and molecular results, see NJ tree in Suppl. material 3. These have been deliberately included in the dataset and are discussed in the results of taxonomy chapter in Suppl. material 4. Most of these cases await synonymisation or constitute very young species or species complexes challenging to disentangle with the CO1 gene. The following references are used exclusively in the suppl. material 4: Dieckmann (1979), Freude et al. (1981), Freude et al. (1983), Stuben (1994), Funk and Omland (2003), Bahr et al. (2008), Skuhrovec (2009), Stuben and Astrin (2010), Stuben et al. (2012), Stuben et al. (2013b), Stuben (2014a), Schutte and StUben (2015), Stuben et al. (2016b). Neighbour-Joining tree The NJ tree with the complete MWI dataset is shown in Suppl. material 3. At genus level, the NJ tree shows a very high congruence with the initial morphological identifications. Even higher taxa are mostly recovered as monophyletic and often cluster in a very similar way as reconstructed in Bayesian analysis, although NJ is neither a phylogenetic method nor is the mitochondrial CO1 gene alone considered suitable to resolve the relationships of higher taxa due to genetic saturation (Arbogast et al. 2002, Hebert and Gregory 2005). Additionally, at the species level, the neighbour-joining clustering algorithm delivers results that are consistently concordant with the Bayesian Inference. 12 Schutte A et al Misidentified specimens are easy to spot in trees when embedded into a matrix of congeneric sequences — misidentified singletons are much more difficult to detect. Beyond misidentified specimens, conflicts can be caused by cryptic species or unresolved synonyms. Several of such inconsistencies have been clarified by taxonomists of the Curculio Institute over the last years, especially in Cryptorhynchinae, Ceutorhynchinae and Apioninae (see !ntroduction), thus delivering a cleaner picture for this barcode release. Bayesian trees The Bayesian consensus trees focusing on three groups within the dataset are provided in Suppl. material 5: Cryptorhynchinae and Cossoninae with 1311 sequences; Apioninae, Nanophyinae and Attelabidae with 367 sequences; Ceutorhynchinae with 538 sequences. The Bayesian posterior probabilities mostly show full or at least high (> 90) support in between species. The phylogenetic trees show substantially more polytomies than the phenetic NJ tree. Nevertheless, taxon placements with regard to the closest related species in the dataset mostly coincide between both methods or have marginal deviations. Thus, the Bayesian tree overall confirms the morphological species identifications and also the naming contradictions, based on unresolved taxonomic issues in the same way as the NJ tree. DiStats analysis (p-distance values) The DiStats statistics are presented in Table 3 (Cryptorhynchinae), Table 4 (Apioninae) and Table 5 (Ceutorhynchinae). For each genus the minimum and the average value of the smallest available distance value to the closest congener is provided. The results per genus are separated into the four distribution groups defined in Table 2, referring to the size of the distribution area: island distribution (ISL), continental endemic (C1), medium distribution (C2) and large distribution (C3). For genera of the subfamily Cryptorhynchinae, the average distance value between the closest available congener (often sister species) ranges from 3.8% (Si/vacalles) to 19.9% ( Torneuma). For genera of Apioninae, the average distance between the closest related congener ranges from 1.7% (Taeniapion) to 18.2% (Pseudoperapion). For genera of Ceutorhynchinae, the average distance between the closest related congener ranges from 6.4% (Hesperorrhynchus) to 17.8% (Scleropterus). For some genera, the smallest p-distance value is significantly lower than the average one. This is often caused by just a single specimen within the genus. For example, in Exapion, the average distance between species is 7.0%, while the lowest value between two species is 2.3%. The closest conge Accept ner pair in this case is Exapion compactum vs. Exapion uliciperda. All taxa and their closest congeners are listed in the spreadsheets in tab "DiStats_results" (Suppl. material 6). Molecular Weevil Identification Project: A thoroughly curated barcode release ... 13 Table 3. Summarised DiStats results for genera of Cryptorhynchinae. Numbers indicate uncorrected p- distance values (genetic distances) expressed in percent. Two values are given per genus and distribution range: 1. minimum distance to the closest congener within all species in the dataset and 2. average distance to the closest congener within all species in the dataset. Abbreviations: min. = minimum, dist. = distance. Cryptorhynchinae island(s)/ island(s)/ endemic endemic medium medium large (> _ large (> : : (50 km) (50km) (50-500 (50-500 500km) 500 km) archipel _—_archipel km) km) Distribution group ISL ISL C1 C1 C2 C2 C3 C3 Genus min. dist. min. dist. min. dist. min. dist. min. dist. min. dist. min. dist. min. dist. to closest to closest toclosest toclosest to closest toclosest to closest to closest congener congener congener congener congener congener congener congener Acalles 8.4 9.0 3.2 7.0 6.2 10.6 7.1 11.7 Acallocrates 13.0 13.3 Acallorneuma 5.8 7.3 3.0 6.8 8.8 8.8 Aeoniacalles 8.8 9.1 Calacalles 3.2 6.0 14.3 14.3 Canariacalles 6.4 6.4 Caucasusacalles 15.3 15.3 Coloracalles 12.2 12.2 Dendroacalles 3.7 7.3 Dichromacalles 7.3 7.3 12.3 13.4 12.9 13.7 Echinodera 3.3 9.6 5.8 11.6 5.8 11.4 6.1 10.8 Echiumacalles 6.8 6.8 Elliptacalles 7.1 7.1 Ficusacalles 6.5 6.5 Kyklioacalles 8.8 8.8 4.3 7.8 4.9 9.0 7.1 10.2 Lauriacalles 10.2 10.2 Madeiracalles 1.8 8.5 Montanacalles 13.7 13.7 Onyxacalles 7.6 8.8 4.0 7.1 4.0 7.1 Pseudodichromacalles 6.4 7.7 Silvacalles 0.9 3.8 14 Cryptorhynchinae Distribution group Genus Sonchiacalles Torneuma Table 4. archipel _—_archipel ISL Schtitte A et al island(s)/ island(s)/ endemic endemic (50km) (50 km) ISL C1 C1 min. dist. min. dist. min. dist. min. dist. to closest to closest to closest to closest medium (50-500 km) C2 min. dist. congener congener congener congener congener 8.3 5.8 8.8 10.2 16.4 16.9 medium (50-500 km) C2 min. dist. to closest to closest congener large (> 500 km) C3 min. dist. to closest congener large (> 500 km) C3 min. dist. to closest congener Summarised DiStats results for genera of Apioninae. Numbers indicate uncorrected p-distance values (genetic distances) expressed in percent. Two values are given per genus and distribution range: 1. minimum distance to the closest congener within all species in the dataset and 2. average distance to the closest congener within all species in the dataset. Abbreviations: min. = minimum, dist. = distance. Apioninae Distribution group Genus Aizobius Alocentron Apion Aspidapion Catapion Ceratapion Cistapion Cyanapion Diplapion Eutrichapion Exapion Hemitrichapion Holotrichapion island(s)/ archipel ISL min. dist. to closest congener 44 10.8 8.1 5.0 island(s)/ endemic endemic archipe! (59Km) (50 km) ISL C1 C1 min. dist. min. dist. min. dist. to closest toclosest to closest congener congener congener 44 10.8 8.1 8.7 medium (50-2000 km) C2 min. dist. to closest congener 10.5 medium (50-2000 km) C2 min. dist. to closest congener 10.5 large (> 2000 km) C3 min. dist. to closest congener 10.8 10.5 6.3 4.9 8.5 4.4 15.8 10.9 3.0 11.0 2.3 12.6 7.5 large (> 2000 km) C3 min. dist. to closest congener 10.8 10.5 8.6 4.9 11.2 11.0 13.5 10.3 Molecular Weevil Identification Project: A thoroughly curated barcode release ... Apioninae island(s)/ _ island(s)/ archipel archipel Distribution ISL ISL group Genus min. dist. min. dist. to closest to closest congener congener Ischnopterapion Ixapion Kalcapion 4.3 4.7 Lepidapion 2.7 2.7 Loborhynchapion Malvapion Omphalapion Onychapion Oryxolaemus Oxystoma Perapion Phrissotrichum 8.2 8.2 Protapion Protopirapion Pseudapion Pseudaplemonus Pseudoperapion Pseudoprotapion Pseudostenapion Rhopalapion Stenopterapion Synapion Taeniapion 4.2 7.0 Taphrotopium Trichopterapion medium (50-2000 km) C2 min. dist. to closest congener 13.2 2.1 7.8 12.6 medium (50-2000 km) C2 min. dist. to closest congener 13.2 2.1 7.8 12.6 large (> 2000 km) C3 min. dist. to closest congener 13.3 12.6 4.1 14.3 15.5 18.2 14.6 17.6 11.6 12.6 13.2 1.7 13.4 16.4 15 large (> 2000 km) C3 min. dist. to closest congener 13.3 12.6 4.1 16.4 16 Table 5. Schtitte A et al Summarised DiStats results for genera of Ceutorhynchinae. Numbers indicate uncorrected p- distance values (genetic distances) expressed in percent. Two values are given per genus and distribution range: 1. minimum distance to the closest congener within all species in the dataset and 2. average distance to the closest congener within all species in the dataset. Abbreviations: min. = minimum, dist. = distance Ceutorhynchinae_island(s)/ archipel Distribution group = ISL Genus min. dist. to closest congener Aphytobius Auleutes Barioxyonyx Brachiodontus Ceutorhynchus 5.5 Coeliodinus Datonychidius Drupenatus Eucoeliodes Eubrychius Glocianus Hadroplontus Hesperorrhynchus 5.0 Homorosoma Marmaropus Mesoxyonyx Micrelus Microplontus Mogulones/ 7.8 Datonychus Mogulonoides Neoglocianus island(s)/ endemic endemic archipel (Km) (50 km) ISL C1 C1 min. dist. min. dist. min. dist. to closest toclosest to closest congener congener congener 9.8 6.4 8.7 medium (50-2000 km) C2 min. dist. to closest congener 7.6 14.1 15.7 6.8 15.4 14.6 14.6 11.6 medium (50-2000 km) C2 min. dist. to closest congener 7.6 14.9 15.8 9.8 15.4 14.6 14.6 13.1 large (> 2000 km) C3 min. dist. to closest congener 7.6 15.7 4.1 13.7 13.5 12.1 13.2 10.9 14.3 15.9 11.9 12.3 6.1 13.4 10.7 large (> 2000 km) C3 min. dist. to closest congener 7.6 15.7 9.1 13.7 13.5 12.1 13.3 10.9 14.3 15.9 13.4 10.7 Molecular Weevil Identification Project: A thoroughly curated barcode release ... 17 Ceutorhynchinae island(s)/ island(s)/ endemic endemic medium medium large (> large (> archipel | archipel (50 km) (50 km) (50-2000 (50-2000 2000km) 2000 km) km) km) Distribution group = ISL ISL C1 C1 C2 C2 C3 C3 Genus min. dist. min. dist. min. dist. min. dist. min. dist. min. dist. min. dist. min. dist. to closest toclosest toclosest toclosest toclosest toclosest toclosest to closest congener congener congener congener congener congener congener congener Neophytobius 14.9 14.9 Oprohinus 15.1 15.1 Oreorrhynchaeus 12.5 12.5 Parethelcus 8.6 8.6 Paroxyonyx 9.7 11.3 9.7 10.9 Pelenomus 11.5 13.4 Perioxyonyx 16.0 16.0 Phrydiuchus 12.2 12.2 10.5 11.5 Poophagus 14.4 14.4 Prisistus 15.7 15.7 12.9 12.9 Ranunculiphilus 12.8 12.8 Rhinoncus 8.4 8.4 7.0 10.1 Scleropterus 17.8 17.8 Scleropteridius 15.9 15.9 Sirocalodes 10.6 10.6 10.6 10.6 Thamiocolus 13.1 13.2 12.8 13.8 Trichosirocalus 12.0 12.0 12.6 13.6 Zacladus 12.2 12.2 ASAP analysis See Table 6 for summarised results of the ASAP calculation and subsequent evaluation of concordance between calculated MOTUs and morphospecies; see Suppl. material 7 for data assembly and full length results. Based on the programme's calculation, the ASAP- score is not a reliable identifier for the confidence level of the suggested threshold (ASAP- score: "the lower, the better"), at least not for the three sub-datasets: Cryptorhynchinae partition 1 with the lowest ASAP-score of "8.5" shows 19% wrongly assigned taxa, while partition 6 with an ASAP-score of "17" shows 16% wrongly assigned taxa (3% less). The ASAP-score does not deliver helpful information here. Apioninae partition 1 with the lowest ASAP-score of "3" shows 14% of wrongly assigned taxa, while partition 6 with an 18 Schutte A et al ASAP-score of "9.5" shows just 7% of wrongly assigned taxa. Thus, the partition/threshold with the higher ASAP-score is providing a much better threshold for that subfamily. Ceutorhynchinae partition 1, with the lowest ASAP-score of "1", shows 6% wrongly assigned taxa, while Partition 10, with an ASAP-score of "18.5", also shows 6% wrongly assigned taxa. The suggested p-distance threshold value for partition 1 is 5.1%, while it is 5.5% for partition 10. The threshold values only show minor differences and a similar amount of wrongly assigned taxa are to be expected. Still, one time, the ASAP-score is the lowest (highest confidence) and the other time, the highest (least confidence). On the contrary, partition 3, with a relatively high ASAP-score of "11" and a threshold value of 5.0% (nearly the same as partition 1), shows the lowest percentage of wrongly assigned taxa (5% error rate). Thus, the partition/threshold with the higher ASAP-score can provide a better one. Table 6. Left side of table: summarised ASAP results, right side of table: evaluation of concordance between MOTUs and morphospecies. For each subfamily dataset, 10 different thresholds ("ASAP partitions") and derived MOTUs are calculated by ASAP. The evaluation of concordance provides the deviations between MOTUs and morphospecies for each given threshold; wrongly assigned MOTUs are given in absolute numbers and in percent. Marked tables point to the threshold which fits best to each subfamily dataset (lowest number of deviation between MOTUs and morphospecies). ASAP ASAP ASAP ASAP evaluation of evaluation of evaluation of results results results results concordance concordance concordance ASAP MOTUs~ Threshold ASAP- no of wrongly no of wrongly % of wrongly Partition [%] score assigned taxa assigned seqs. assigned seqs. Cryptorhynchinae sub-dataset (contains 265 morphospecies, 1106 sequences) 1 236 7.3 8.5 74 214 19% 2 241 7.2 9.0 73 214 19% 3 315 3.8 9.5 84 190 17% 4 348 2.4 15.0 104 206 19% 5 639 0.4 16.5 373 480 43% 6 251 6.8 17.0 68 181 16% 7 325 3.4 24.5 86 193 17% 8 316 3.7 25.5 83 186 17% 9 302 4.3 29.0 84 205 19% 10 329 3.1 29.5 89 189 17% Apioninae sub-dataset (contains 114 morphospecies, 342 sequences) 1 95 6.0 3.0 19 47 14% Molecular Weevil Identification Project: A thoroughly curated barcode release ... 19 ASAP ASAP ASAP ASAP evaluation of evaluation of evaluation of ecu? results sulls rps ults concordance concordance concordance ASAP MOTUs~ Threshold ASAP- no of wrongly no of wrongly % of wrongly Partition [%] score assigned taxa assigned seqs. assigned seqs. 2 93 7.3 4.5 21 52 15% 3 92 7.5 5.0 22 52 15% 4 87 8.2 5:5 28 61 18% 5 94 6.8 7.5 17 41 12% 6 111 3.8 9.5 13 25 7% 7 88 8.1 11.0 24 57 17% 8 129 1.9 14.5 24 29 8% 9 116 3.0 14.5 16 28 8% 10 112 3.3 15.0 13 27 8% Ceutorhynchinae sub-dataset (contains 199 morphospecies, 491 sequences) 1 204 5.1 1.0 17 28 6% 2 191 6.9 7.0 19 38 8% 3 206 5.0 11.0 15 24 5% 4 186 7.7 11.5 19 40 8% 5 235 2.2 13.0 37 59 12% 6 190 7.1 13.5 18 37 8% 7 178 8.5 14.0 23 64 13% 8 183 7.8 16.0 20 42 9% AS) 178 8.6 16.5 23 46 9% 10 203 5.5 18.5 18 31 6% Discussion The present DNA barcode release provides results for almost 1300 Western Palearctic weevil taxa. This dataset's strength lies in its thorough validation of specimens, including the actual nomenclatorial resolution of many cases of previous taxonomic conflicts (in preceding publications within the MWI project). The correct identifications are mirrored in a high consistency between morphological identifications and molecular results. The ambiguous cases where molecular and morphological evidence could not be reconciled are discussed (Suppl. material 4). These conflicts mainly have their basis in pending 20 Schtitte A et al synonymisations or are caused by species complexes that cannot be resolved via DNA barcoding. In most cases, initial discrepancies between morphological identification and molecular results could be resolved by confirming or falsifying initial identification. In the latter case, additional resampling was needed and, in some situations, holotype comparisons, partly leading to taxonomic changes. DiStats statistics. The most densely sampled subfamilies in the dataset often show genus-specific distances between species; see Suppl. material 6. Within each subfamily in focus, there are genera with small mean interspecific distances and others with high mean interspecific distances. For example, the Cryptorhynchinae dataset shows an average distance of 3.8% in the genus Si/vacalles, while in the genus Torneuma, the average lies at 10.2% (both with island distribution). In Apioninae, the interspecific average ranges from 1.7% (Taeniapion) to 18.2% (Pseudoperapion; both taxa have a large distribution). In Ceutorhynchinae, it ranges from 6.4% (Hesperorrhynchus, island distribution) to 17.8% (Scleropterus, medium distribution, continental). Contrary to our initial expectations, it is clear that there does not exist a single threshold per subfamily that would characterise usual species limits. We also expected the species’ distribution scales to be correlated in some form with the average genetic distances between species, but this is not the case either. For example, in the genus Echinodera, the average distances for the mainland distribution groups are 11.6% (small, endemic), 11.4% (medium), 10.8% (large), showing a slightly decreasing tendency from small to large geographical distribution. For the genus Kyklioacalles, this tendency is reverted: 7.8% (small, endemic), 9% (medium), 10.2% (large). By summarising the statistical findings, it can be concluded that applying a single general genetic threshold for species delineation leads to mismatches between morphospecies and MOTUs, either false positives (oversplits) or false negatives (lumps). These mismatches are also clearly demonstrated in the ASAP results (see Table 6 and Suppl. material 7), with 10 widely-varying thresholds. No matter whether an increased or decreased threshold is applied, there remain significant deviations between morphospecies and calculated MOTUs, although no questionable taxa are included in the ASAP sub-datasets. Targeting alpha-taxonomic questions with a single threshold approach likely leads to unsatisfactory error rates between 5% and 43% (see ASAP results in Table 6, right column). ASAP or other single-threshold approaches are a convenient option to estimate species richness in widely-unknown biota or when there is no option to resort to using morphological information. Additionally, within a rough biodiversity assessment (e.g. metabarcoding), a small taxonomic error rate might not distort the final result. However, incorrect identifications can subsequently be incorporated into further studies. In the worst case, long-term environmental programmes could generate error cascades which can have a negative impact on environmental management and conservation (Bortolus 2008, Isaac et al. 2004). Molecular Weevil Identification Project: A thoroughly curated barcode release ... 21 It is known that undersampling leads to artificially increased interspecific genetic distances, creating deeper splits in trees and wider barcoding gaps (Moritz and Cicero 2004, Morando et al. 2003). Undersampling leads likewise to higher interspecific genetic distances in the DiStats statistics (Tables 3, 4, 5). This study's species coverage of the subfamilies Apioninae and Ceutorhynchinae is far from complete. The dataset contains 115 (29%) of roughly 400 binomial Apioninae taxa and 204 (51%) of roughly 400 binomial Ceutorhynchinae taxa of the Western Palearctic (L6bl and Smetana 2011, LOobl and Smetana 2013). By adding further sequences to the dataset, we assume the minimum average p-distance values will decrease significantly for some genera. Unlike in Apioninae and Ceutorhynchinae, the genera-specific distances of the Cryptorhynchinae dataset are not strongly affected by undersampling. We have covered 278 of 384 (72%) currently known binomial Cryptorhynchinae taxa (Stuben 2018b). From our observation, if the sister species is missing in the dataset, the remaining congeners are marginally more distant (higher p-distance values). For example, the adelphotaxon to Acalles granulimaculosus is Acalles pilula, both taxa are included in the dataset and they show a genetic distance of 11.3%. By removing Acalles pilula from the dataset, the closest congener is Acalles globulipennis, with a genetic distance of 12.3%. Thus, the minimum interspecific genetic distance for this taxon would increase by 1% in the DiStats output data table and the calculated average minimum distance for the genus Acalles from 9% to 9.2% in the island species compilation (Table 3). Besides missing several species, sampling usually could not cover the entire geographic distribution of most continental taxa in our dataset. Complete sampling within each species’ full geographic distribution range would likely reveal higher intraspecific distances than we can observe in the trees. Thus, we have not focused on intraspecific variation for the time being. Bergsten et al. (2012) found the intraspecific variation in Agabini diving beetles to be significantly correlated with the geographic scale of sampling: up to 70 individuals were required to sample 95% of the intraspecific variation. Future collecting of weevils on the Western Palearctic mainland should bear this in mind and should strive to fill the mentioned gaps. The Canary Islands were extensively sampled. The dataset usually contains at least one specimen per taxon from each island (for multi- island distributions) or several collecting spots per island (for endemic/single island distributions). Most species occurring in the Canaries do not occur on the mainland. Many mainland species, however, especially in Apioninae and Ceutorhynchinae, occur far beyond the Western Palearctic. Species distribution maps for the three subfamilies in focus are provided as external supplement under DOI 10.5281/zenodo.7430565. The many different examples from the literature (Coyne and Orr 1997, Sasa et al. 1998, Presgraves 2002, Mendelson 2003, Hebert et al. 2003a, Hebert et al. 2003b, Hebert et al. 2004a, Hebert et al. 2004b, Barrett and Hebert 2005, Monaghan et al. 2005, Vences et al. 2005, Zigler et al. 2005, Hickerson et al. 2006, Lefebure et al. 2006, Mikkelsen et al. 2007, Wiemers and Fiedler 2007, Hubert et al. 2008, Hundsdoerfer et al. 2009, Robinson et al. 2009, Acs et al. 2010, Hubert et al. 2010, Candek and Kuntner 2015, Spasojevic et al. 2016) clearly demonstrate that different groups have different genetic variabilities (see Suppl. material 8 for a summary of the previously mentioned references). Useful thresholds 22 Schutte A et al to delineate animal species boundaries using CO1 barcodes are often found between 4% and 15% genetic distance. However, almost every reference has deviant taxa with larger or smaller values. Besides subsampling, a taxon's age, modes of evolution and reproduction and its current and historical distribution (range size, climate, isolation barriers) can shape such distances. No single threshold would hold for all taxa, but knowing the average and minimum genetic distances between species within a specific (taxonomically limited group of interest) can be a vital supplementary tool in resolving taxonomic issues. For the weevil subfamilies we studied, we found that the taxonomic level for defining a common barcode threshold cannot meaningfully be established above the genus. Applying a threshold on genus level usually allows comparisons of species sharing a similar distribution pattern, evolutionary age and/or occupying similar ecological niches. The thresholds we propose to use as heuristic support tools for future research in these groups are well-calibrated morphologically and based on "good species". These values can provide a reference for future alpha-taxonomic weevil research consistent with the definitions or understanding of existing species. Genetic distances are easy to measure, but prompt the question of how much distance is needed to delineate species. On the morphological side, each weevil group has its own set of particular characteristics used for morphological identification, which specialists have agreed upon over time, often somewhat subjectively. Morphological variations (intraspecific and interspecific) have been the basis of discussions in taxonomy ever since and are crucial to study prior to a taxonomic change. Those morphological characters used for identification and delineation are mostly based on the consensus principle of the scientific community. Essential characters in one group, bristle length, for example, might not play any role in another group, where perhaps the colouration pattern on the elytra or the protrusion of the eyes may constitute the central diagnostic characters. Based on many years of experience, a specialist will Know about those morphological characters in his/her studied group. Known morphological variability within species is well factored in when examining differential characters for a new species description. Comparable situations arise when NJ trees or their underlying genetic distances are discussed (molecular intraspecific and interspecific variation), for example, when a single species appears in two neighbouring clusters. Based on a solely molecular point of view, those clusters might be separated by sufficient distance to infer the existence of a new species. Still, the morphologist may know from experience that those two clusters belong to a geographic variation. A - mostly historical - quantitative approach for species delineation was morphological phenetics or "numerical taxonomy" (Sneath and Sokal 1973). In weevils, for example, comparative results were published to delineate /schnopterapion modestum and Ischnopterapion plumbeomicans with dozens of minimum, maximum and average length measurements from each body part of the two weevil species (Ehret 1991). However, phenetics is unable to recover evolutionary relationships as it does not differentiate between homology and homoplasy. Thus, it has been substituted with computational methods which can deliver an approximation of phylogeny (VWVagele 2005). New species descriptions, based solely on DNA barcoding, have been carried out or at least suggested for cryptic species or species complexes soon after DNA barcoding was Molecular Weevil Identification Project: A thoroughly curated barcode release ... 23 established (Brower 2010, Cook et al. 2010, J6rger and Schrédl 2013, Doerder 2018). This strategy has recently been used for hyperdiverse taxa from the tropics as well (Butcher et al. 2012, Riedel et al. 2013a, Riedel et al. 2013b, Meierotto et al. 2019, Sharkey et al. 2021 ). Providing only a DNA barcode (Doerder 2018) or a consensus sequence (Sharkey et al. 2021) without thoroughly investigating the (group-specific) interspecific distances means postponing the (molecular) differential diagnosis. Discriminatory positions as differential diagnosis (proposed by Cook et al. (2010), applied by Jorger and Schrddl (2013) and Meierotto et al. (2019)) involve some uncertainty. They might become invalid as more specimens are collected, covering additional intraspecific variation from the same species or from closely-related species. Nevertheless, these approaches may constitute a way to capture the biodiversity of underdescribed taxa quickly. Riedel et al. (2013b) state: "A combination of digital imaging and molecular techniques allows the reduction of formal species descriptions to brief but highly accurate diagnoses. Although none of these tools is novel in itself, the progressive element is their combination and streamlining to produce a large number of usable species descriptions," provided the identifier has access to a sequencing facility and sufficient morphological knowledge to seek and find differential characters on the pictures provided. Yet, excluding morphology is not commonly accepted in the scientific community (Pante et al. 2015, Ahrens et al. 2021, Zamani et al. 2022), even not for protists (Warren et al. 2017). Using a universal genetic distance threshold to compellingly delineate species would decouple taxonomy from the previously-established systematics. It would lead into the direction of a "parallel taxonomy", which has been rejected 20 years ago (see Sperling (2003) on Tautz' DNA taxonomy concept in Tautz et al. (2002) and Tautz et al. (2003)). Relying solely or predominantly on DNA barcodes for species descriptions promises a turbo taxonomy (Butcher et al. 2012) or fast-track taxonomy (Riedel et al. 2013a). It seems appealing when morphology reaches its limits or performance increase is in focus (Fernandez-Triana 2022). In the long run, a trade-off between molecular quick-wins and morphological expertise may occur, for example, how to examine type specimens. Waiving morphological diagnoses in taxonomically challenging cases will most likely supersede conventional species descriptions soon if precautions are not being taken. The taxonomic inflation issue was addressed before DNA barcoding was introduced. Concerns were based on the practice of raising taxa from subspecies to species level, thus resulting in a change of the species concept rather than new species discoveries (Isaac et al. 2004, Rylands and Mittermeier 2014). This issue gains special relevance in light of DNA-based approaches. Even moderate differences in genetic distances between clusters of individuals can reach high statistical significance (Galtier 2019, Vences 2020), further prompting taxonomists to rank such clusters as subspecies or species (Hey 2009). However, we propose subspecies should not be seen as anything other than heuristic guides. They should prompt the community to consult an integrative array of methods, such as morphological, molecular, biogeographic, ecological or ethological, prior to elevating subspecies to species status. This would make decisions on species status much more sustainable. 24 Schutte A et al Adding nuclear markers in combination with phylogenetic models like multi-species coalescent model improve the accuracy of species delineation drastically (Eberle et al. 2020, Miralles et al. 2020, Dietz et al. 2021). They can easily uncover mitochondrial introgression, a possible downside of extrachromosal inheritance of the CO1 barcoding gene. Within the dataset of this study, we encountered two conspicuous cases of introgression: Hesperorrhynchus linaeotesselatus (Stuben and Kratky 2016) and Cionus griseus (Stuben and Behne 2015, Stuben et al. 2021, see further information in Suppl. material 4). Still, simply adding genes cannot assist in developing a molecular species concept (Galtier 2019) and does not answer questions regarding genetic thresholds in species delineation, since speciation circumstances (e.g. divergence times, population size) and sampling depth affect the dataset in the same way. Barcoding the holotype (as non-destructively as possible) is the gold standard in molecular taxonomy. This has been done for some MWI holotypes, for example, Madeiracalles beelzebubi Stuben & Kratky, 2018 or Torneuma alexi Stuben, 2018 (both described in Stuben (2018b)). Only a holotype provides an objective link to its Linnean binomen. Type material is often difficult to access. Most type-holding institutions still offer on-site inspections, which requires travelling. Often they only refer to photos of the specimens available on their web page. Loaning is still offered in some cases, but comes with waiting periods (sometimes even several years, based on the second author’s experience). To send out type specimens is time-consuming and exposes them to the risk of getting lost. If barcodes of most holotypes were openly available, requesting shipping of type material could be omitted in many cases. Using a paratype specimen to retrieve the DNA barcode is the second best choice if the collecting location matches the holotype's and sympatric occurrence is unlikely, for example, Torneuma korwitzi (Stuben and Schutte 2015). Surprisingly, the latter is rarely discussed in literature. Most insects were described before the advent of molecular tools. Barcoding of historical type material (hDNA) is not a new idea, but still controversial due to often invasive processing of the most valuable collection specimens (Townson et al. 1999, Mayer et al. 2021, Raxworthy and Smith 2021). However, we should be aware that retroactive barcoding of type material facilitates robust and sustainable knowledge gain in taxonomy to solve existing and future research questions (Strutzenberger et al. 2012, Prosser et al. 2015, Speidel et al. 2015, Hausmann et al. 2016, Scherz et al. 2020, Raxworthy and Smith 2021, Roycroft et al. 2022, Mulcahy et al. 2022). A final option to gain reliable DNA barcodes after the new species has been described, is recollecting specimens from the type locality (Jorger and Schrédl 2013, Bell et al. 2020), which was also one focus of the MWI project. During the past 250 years, almost every taxonomic change was based on morphological characters, continuously re-evaluating the underlying morphological characters. Hence for weevils, we can assume this "cleanup process" has built a strong foundation of valid morphological characters in most cases. We suggest preserving the already established and globally-accepted Linnean understanding of species as taxonomic backbone. This will ease the progression from a morphology-based past into a strong molecular-based future taxonomy, which will be compatible with the past. The risk of a disjunct parallel taxonomy would be decreased and the potential taxonomic inflation restrained to a minimum. The morphologically calibrated genus-specific distance values, based on "good Molecular Weevil Identification Project: A thoroughly curated barcode release ... 25 species" (Tables 3, 4, 5), constitute a somewhat reliable direction for species of the molecularly well-sampled subfamily Cryptorhynchinae and an initial approximation for Apioninae and Ceutorhynchinae. In general, the results of the molecular dataset (CO1) should be utilised in an integrative taxonomy approach, i.e. discussed with morphological and ecological aspects, geographical distribution and evolutionary age of the taxon in focus. Here, we should address some pitfalls to prevent future inflationary species descriptions: 1. Ignoring the minimum interspecific distances of the sister species. The interspecific genetic distances for weevils are mostly group-specific. A group can be a genus (e.g. Jorneuma or Silvacalles), but it can also be a subgenus (e.g. subgenus Euphorbioacalles of the genus Dendroacalles) or even a species complex (e.g. Acalles maraoensis complex). The interspecific distances of the sister species should be considered. If no sister species pair is available in the dataset, the closest congeners can be taken for an approximation. Otherwise, newly-collected specimens originating from a different population might be potentially classified as new species. Even small distances can create a split in a tree and might justify a new species description at first glance. If the interspecific distance of the potential new species falls below the previously known minimal one, the researcher should be cautious not to describe a synonym. A description can still be carried out if strong reasons justify the new species (Stuben and Schutte 2015, Garcia et al. 2019), for example, the young age of the new species or clear differences in morphologic characters under strong selective pressure, for example, in the aedeagi. 2. Single sequences per taxon or population. Using a single sequence per taxon drastically increases the risk of wrong conclusions when applied to alpha-taxonomic questions because intraspecific variation is not shown, but can be high for some taxa. In addition to the increased risk of misidentifications in singletons, not including intermediate specimens (of the same species) can create an artificial split in a tree which could be misinterpreted as a newly-discovered species, especially if the analysed individual was collected far from the previously-known sequence. If a species has a disjunct distribution, providing just one sequence from each population increases the likelihood of producing a synonym. This risk especially applies to islands. Artificial deep splits can be produced if the intraspecific distances within a population coincide with or even exceed the interspecific distances. On dataset compilation, the full sampling depth should be used. Using a single sequence from each population instead of all available sequences means leaving out all intermediate specimens belonging to the same species. The intraspecific distances then present themselves as an artificial deep split. The latter might be the case for some Laparocerus taxa described recently (Faria et al. 2016, Machado et al. 2017). For a more detailes explanation see Stuben (2022). 3. Gaps in existing sequence databases. A large genetic distance to the closest congener in a sequence database is not proof of having discovered a new species. Often, no reference sequences of the sequenced species have been previously deposited. Subsequently, a misinterpretation of the interspecific genetic distance to the closest 26 Schutte A et al database match, for example, 15% to the closest deposited one, can lead to describing a synonym, particularly if the sister species' type material is not consulted. Although potentially new species can be discovered very quickly with DNA barcoding (Riedel et al. 2013a, Meierotto et al. 2019, Sharkey et al. 2021), they must be taxonomically secured in the same way as traditionally done, at least for well-revised taxa: either by comparing all closely-related species morphologically to the potential new species and/or (especially where the former option is not an option) or by consulting barcodes from all closely-related species. Unfortunately, a described species lacking this validation can be laborious to refute. Most of the work (type comparison, sequencing) has to be carried out by a third party if the original author failed to do so. Unfortunately, it can be assumed that there will be a relatively large number of quickly described species in the future, only based on barcodes and not backed up by holotype comparisons. This situation can arise in island biota, for which one quickly tends to assume endemism: a candidate species is discovered by molecular means and described without holotype comparison to the continental fauna (e.g. in Garcia et al. (2022)). In this context, the candidate species carries the risk of being a synonym, because the species might have been described already from the mainland. The wrong conclusion "it must be a new species" is quickly made if no public sequence is available. Vice versa, new descriptions should always be supported by molecular data to prevent describing a synonym — also by other researchers at a later point. Sometimes, formerly established unique morphological traits used for species description can turn out to be misleading diagnostic characters after molecular data become available. The cryptorhynchine species Calacalles agana Stuben, 2010 (StUben 2010) can serve as an example: the author described the species without molecular support and synonymised it several years later (StUben 2015). Another case with conflicting results between morphology and DNA barcoding is Aeoniacalles aeonii bodegensis (see Suppl. material 9, StUben (2005), Stuben and Germann (2005), Stuben and Astrin (2011)). Following the biological species concept (Mayr 1942), one could consider conducting cross-breeding experiments (StUben 2005) for some generations prior to a new species description. External supplementary material Data Type: geographical distribution maps. Brief description: the ZIP file contains 613 distribution maps from Western Palearctic weevil taxa. The distribution maps showing Europe originate from the Curculio Institute's website (www.curci.de). Additional information on distribution range and known synonyms were based on the information from the L6bl catalogues (L6bIl and Smetana 2011, Lo6bl and Smetana 2013). The maximum distribution range of each species was measured in km with Google Earth's ruler function. Download via Zenodo DOI: 10.5281/zenodo.7430565 (368.1 MB). Data Type: Material Table and CO1 sequences. Brief description: alternative download source for the material table and the CO1 sequences used in this study. Download via Zenodo DOI: 10.5281/zenodo.7430106 (3.9 MB). Molecular Weevil Identification Project: A thoroughly curated barcode release ... 27 Acknowledgements This study was based on the genuine dedication of many Curculio Institute members. They took a significant role in laborious collecting, identifying and meticulously discussing the taxonomic results, especially Peter Sprick (Germany), Jiri Kratky (Czech Republic), Lutz Behne (Germany) and Christoph Bayer (Germany). We are extremely grateful to the laboratory technicians from the Museum Koenig who carried out the main part of the lab work and undertook every possible effort to get results from almost all samples, namely Hannah Petersen, Christina Blume and Laura von der Mark. We highly appreciate the cooperation with Katja Kramp and Eva Kleibusch from the Senckenberg German Entomological Institute (SDEI), who provided additional 91 sequences to this DNA barcode release. We would also like to thank all the authorities who issued collection permits, especially from the Canary Islands and the Azores. The publishing fees were partly covered by the Leibniz association's open access publishing fund. Grant title The samples collected in Germany were processed in the German Barcode of Life-Project (GBOL), supported by the German Federal Ministry of Education and Research, Berlin, Germany (BMBF FKZ 01L11101 & FKZ 01L11501A). Hosting institution Leibniz Institute for the Analysis of Biodiversity Change, Museum Koenig, Adenauerallee 160, 53113 Bonn, Germany. Curculio Institute - Center for Studies on western Palearctic Curculionoidea, Curculio- Institut e.V. (CURCI), Hauweg 62, 41066 Monchengladbach, Germany. References ° Acs Z, Challis RJ, Bihari P, Blaxter M, Hayward A, Melika G, Csoka G, Penzes Z, Pujade-Villar J, Nieves-Aldrey JL, Schonrogge K, Stone GN (2010) Phylogeny and DNA barcoding of inquiline Oak gallwasps (Hymenoptera: Cynipidae) of the Western Palaearctic. Molecular Phylogenetics and Evolution 55 (1): 210-225. https://doi.org/ 10.1016/j.ympev.2009.12.004 28 Schtitte A et al Ahrens D, Ahyong ST, Ballerio A, Barclay MV, Eberle J, Espeland M, Huber BA, Mengual X, Pacheco TL, Peters RS, Rulik B, Vaz-De-Mello F, Wesener T, Krell FT (2021) Is it time to describe new species without diagnoses? - A comment on Sharkey et al. (2021). Zootaxa 5027 (2): 151-159. https://doi.org/10.11646/zootaxa.5027.2.1 Alonso-Zarazaga MA, Barrios H, Borovec R, Bouchard P, Caldara R, Colonnelli E, Gultekin L, Hlavac P, Korotyaev BA, Lyal C, Machado A, Meregalli M, Pierotti H, Ren L, Sanchez-Ruiz M, Sforzi A, Silfverberg H, Skuhrovec J, Tryzna M, Yunakov N (2017) Cooperative catalogue of Palaearctic Coleoptera Curculionoidea. 1st Edition, Vol. 8. Monografias electronicas S.E.A., Zaragoza, Spain, 729 pp. URL: https://zoobank.org/ References/911EF526-33F0-4970-8EC2-A7F5AC1E1D3D Arbogast B, Edwards S, Wakeley J, Beerli P, Slowinski J (2002) Estimating divergence times from molecular data on phylogenetic and population genetic timescales. Annual Review of Ecology and Systematics 33 (1): 707-740. https://doi.org/10.1146/ annurev.ecolsys.33.010802.150500 Astrin JJ, Stuben PE (2008) Phylogeny in cryptic weevils: Molecules, morphology and new genera of Western Palaearctic Cryptorhynchinae (Coleoptera: Curculionidae). Invertebrate Systematics 22: 503-522. https://doi.org/10.1071/ISO7057 Astrin JJ, Stuben PE, Misof B, Wagele JW, Gimnich F, Raupach MJ, Ahrens D (2012) Exploring diversity in cryptorhynchine weevils (Coleoptera) using distance-, character- and tree-based species delineation. Molecular Phylogenetics and Evolution 63 (1): 1-14. https://doi.org/10.1016/j.ympev.2011.11.018 Astrin JJ, Zhou X, Misof B (2013) The importance of biobanking in molecular taxonomy, with proposed definitions for vouchers in a molecular context. Zookeys 365: 67-70. https://doi.org/10.3897/zookeys.365.5875 Astrin JJ, Hofer H, Spelda J, Holstein J, Bayer S, Hendrich L, Huber BA, Kielhorn KH, Krammer HJ, Lemke M, Monje JC, Moriniere J, Rulik B, Petersen M, Janssen H, Muster C (2016) Towards a DNA barcode reference database for spiders and harvestmen of Germany. PLOS One 11 (9): 1-24. https://doi.org/10.1371/journal.pone.0162624 Bahr F, Behne L, Sprick P, Stuben PE, Bayer C (2008) Digital-weevil-determination for Curculionoidea of West Palaearctic. Transalpina: Sitona (Entiminae: Sitonini). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 7 (86): 21-27. URL: https:/Awww.curci.de/?beitrag=86 Barrett RH, Hebert PN (2005) Identifying spiders through DNA barcodes. Canadian Journal of Zoology 83 (3): 481-491. https://doi.org/10.1139/z05-024 Bell R, Mulcahy D, Gotte S, Maley A, Mendoza C, Steffensen G, Barron li J, Hyman O, Flint W, Wynn A, McDiarmid R, McLeod D (2020) The type locality project: Collecting genomic-quality, Topotypic vouchers and training the next generation of specimen- based researchers. Systematics and Biodiversity 18 (6): 557-572. https://doi.org/ 10.1080/14772000.2020.1769224 Bergsten J, Bilton DT, Fujisawa T, Elliott M, Monaghan MT, Balke M, Hendrich L, Geijer J, Herrmann J, Foster GN, Ribera |, Nilsson AN, Barraclough TG, Vogler AP (2012) The effect of geographical scale of sampling on DNA barcoding. Systematic Biology 61 (5): 851-69. https://doi.org/10.1093/sysbio/sys037 Blaxter M (2016) Imagining Sisyphus happy: DNA barcoding and the unnamed majority. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 371 (1702): 1-5. https://doi.org/10.1098/rstb.2015.0329 Molecular Weevil Identification Project: A thoroughly curated barcode release ... 29 Blaxter ML (2004) The promise of a DNA taxonomy. Philosophical Transactions of the Royal Society B-Biological Sciences 359 (1444): 669-679. https://doi.org/10.1098/rstb. 2003.1447 Bortolus A (2008) Error Cascades in the biological sciences: The unwanted consequences of using bad taxonomy in ecology. Ambio 37 (2): 114-118. https://doi.org/ 10.1579/0044-7447(2008)37[114:ecitbs]2.0.co;2 Bouchard P, Smith A, Douglas H, Gimmel M, Brunke A, Kanda K (2017) Biodiversity of Coleoptera: Science and society (Chapter 11). In: Foottit RG, Adler PH (Eds) Insect Biodiversity. 2nd Edition, Vol. 1. John Wiley & Sons, New Jersey, USA, 80 pp. [ISBN 9781118945537]. https://doi.org/10.1002/9781118945568.ch11 Brower AVZ (2010) Alleviating the taxonomic impediment of DNA barcoding and setting a bad precedent: Names for ten species of Astraptes fulgerator (Lepidoptera: Hesperiidae: Eudaminae) with DNA-based diagnoses. Systematics and Biodiversity 8 (4): 485-491. https://doi.org/10.1080/14772000.2010.534512 Butcher BA, Smith MA, Sharkey MJ, Quicke DLJ (2012) A turbo-taxonomic study of Thai Aleiodes (Aleiodes) and Aleiodes (Arcaleiodes) (Hymenoptera: Braconidae: Rogadinae) based largely on COI barcoded specimens, with rapid descriptions of 179 new species. Zootaxa (3457)1-232. https://doi.org/10.11646/zootaxa.3457.1.1 Candek K, Kuntner M (2015) DNA barcoding gap: Reliable species identification over morphological and geographical scales. Molecular Ecology Resources 15 (2): 268-77. https://doi.org/10.1111/1755-0998.12304 Cook LG, Edwards RD, Crisp MD, Hardy NB (2010) Need morphology always be required for new species descriptions? Invertebrate Systematics 24 (3). https://doi.org/ 10.1071/1S10011 Coyne JA, Orr HA (1997) "Patterns of speciation in Drosophila" Revisited. Evolution 51 (1): 295-303. https://doi.org/10.2307/2410984 Dieckmann L (1979) Neue palaarktische Arten aus der Unterfamilie Ceutorhynchinae (Coleoptera, Curculionidae). Reichenbachia, Staatliches Museum fir Naturkunde in Dresden 17 (7): 49-56. Dietz L, Eberle J, Mayer C, Kukowka S, Bohacz C, Baur H, Espeland M, Huber BA, Hutter C, Mengual X, Peters RS, Vences M, Wesener T, Willmott K, Misof B, Niehuis O, Ahrens D (2021) Standardized nuclear markers advance metazoan taxonomy. bioRxiv (preprint) https://doi.org/10.1101/2021.05.07.443120 Doerder FP (2018) Barcodes reveal 48 new species of Tetrahymena, Dexiostoma, and Glaucoma: Phylogeny, ecology, and biogeography of new and established species. The Journal of Eukaryotic Microbiology 66 (1): 182-208. https://doi.org/10.1111/jeu.12642 Drummond AJ, Ashton B, Buxton S, Cheung M, Cooper A, Duran C, Field M, Heled J, Kearse M, Markowitz S, Moir R, Stones-Havas S, Sturrock S, Thierer T, Wilson A (2012) Geneious v5.6. htips://www.geneious.com. Accessed on: 2012-12-31. Eberle J, Ahrens D, Mayer C, Niehuis O, Misof B (2020) A plea for standardized nuclear markers in metazoan DNA taxonomy. Trends in Ecology & Evolution 35 (4): 336-345. https://doi.org/10.1016/j.tree.2019.12.003 Edgar RC (2004) MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research 32 (5): 1792-1797. https://doi.org/10.1093/nar/ gkh340 30 Schtitte A et al Ehret JM (1991) Etude du groupe d'Apion /oti avec la description d'une espéce et d'une sous espece nouvelle (Coleoptera: Curculionidae, Apioninae). Nouvelle Revue d'Entomologie. Nouvelle Serié, Paris 8 (2): 159-180. [In French]. Faria CMA, Machado A, Amorim IR, Gage MJG, Borges PAV, Emerson BC (2016) Evidence for multiple founding lineages and genetic admixture in the evolution of species within an oceanic island weevil (Coleoptera, Curculionidae) super-radiation. Journal of Biogeography 43 (1): 178-191. https://doi.org/10.1111/jbi.12606 Fernandez-Triana JL (2022) Turbo taxonomy approaches: Lessons from the past and recommendations for the future based on the experience with Braconidae (Hymenoptera) parasitoid wasps. ZooKeys 1087: 199-220. https://doi.org/10.3897/ zookeys.1087.76720 Floyd R, Abebe E, Papert A, Blaxter M (2002) Molecular barcodes for soil nematode identification. Molecular Ecology 11 (4): 839-850. https://doi.org/10.1046/j.1365-294x. 2002.01485.x Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R (1994) DNA primers for amplification of mitochondrial cytochrome c oxidase subunit | from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology 3 (5): 294-299. URL: https:// pubmed.ncbi.nlm.nih.gov/7881515/ Freude H, Harde KW, Lohse GA (1981) Die Kafer Mitteleuropas. Band 10. Curculionidae |. 1st, 10. Goecke & Evers, Krefeld, Germany, 310 pp. [ISBN 9783827406842] Freude H, Harde KW, Lohse GA (1983) Die Kafer Mitteleuropas. Band 11. Curculionidae Il. 1st, 11. Goecke & Evers, Krefeld, Germany, 344 pp. [ISBN 9783872630315] Funk D, Omland K (2003) Species-level paraphyly and polyphyly: Frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annual Review of Ecology, Evolution, and Systematics 34 (1): 397-423. https://doi.org/10.1146/ annurev.ecolsys.34.011802.132421 Galtier N (2019) Delineating species in the speciation continuum: A proposal. Evolutionary Applications 12 (4): 657-663. https://doi.org/10.1111/eva.12748 Garcia R, Andujar C, Oromi P, Emerson B, Lopez H (2019) The discovery of Barretonus (Curculionidae: Cossoninae) in the Canary Islands: barcoding, morphology and description of new species. Acta Entomologica Musei Nationalis Pragae 59 (2): 443-452. https://doi.org/10.2478/aemnp-2019-0033 Garcia R, Lopez H, Aguiar A (2022) Ceutorhynchus tenoensis n. sp. de Tenerife (Coleoptera, Curculionidae, Ceutorhynchinae). Revista de la Academia Canaria de Ciencias 33: 133-144. Geiger MF, Astrin JJ, Borsch T, Burkhardt U, Grobe P, Hand R, Hausmann A, Hohberg K, Krogmann L, Lutz M, Monje C, Misof B, Moriniere J, Muller K, Pietsch S, Quandt D, Rulik B, Scholler M, Traunspurger W, Haszprunar G, Wagele W (2016) How to tackle the molecular species inventory for an industrialized nation-lessons from the first phase of the German Barcode of Life initiative GBOL (2012-2015). Genome 59 (9): 661-70. https://doi.org/10.1139/gen-2015-0185 Giangrande A (2003) Biodiversity, conservation, and the ‘taxonomic impediment’. Aquatic Conservation: Marine and Freshwater Ecosystems 13 (5): 451-459. https:// doi.org/10.1002/aqc.584 Molecular Weevil Identification Project: A thoroughly curated barcode release ... 31 Godfray HC (2002) Challenges for taxonomy. Nature 417 (6884): 17-19. https://doi.org/ 10.1038/417017a Godfray HCJ, Knapp S (2004) Introduction (taxonomy for the twenty-first century). Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences 359 (1444): 559-569. https://doi.org/10.1098/rstb.2003.1457 Hajibabaei M, Janzen DH, Burns JM, Hallwachs W, Hebert PD (2006a) DNA barcodes distinguish species of tropical Lepidoptera. Proceedings of the National Academy of Sciences of the USA 103 (4): 968-971. https://doi.org/10.1073/pnas.0510466103 Hajibabaei M, Smith MA, Janzen D, Rodriguez J, Whitfield J, Hebert PN (2006b) A minimalist barcode can identify a soecimen whose DNA is degraded. Molecular Ecology Notes 6 (4): 959-964. https://doi.org/10.1111/j.1471-8286.2006.01470.x Hansen H, Bakke TA, Bachmann L (2007) DNA taxonomy and barcoding of monogenean parasites: lessons from Gyrodactylus. Trends in Parasitology 23 (8): 363-367. https://doi.org/10.1016/j.pt.2007.06.007 Hausmann A, Miller SE, Holloway JD, deWaard JR, Pollock D, Prosser SW, Hebert PD (2016) Calibrating the taxonomy of a megadiverse insect family: 3000 DNA barcodes from geometrid type specimens (Lepidoptera, Geometridae). Genome 59 (9): 671-684. https://doi.org/10.1139/gen-2015-0197 Hebert PD, Ratnasingham S, deWaard JR (2003a) Barcoding animal life: Cytochrome c oxidase subunit 1 divergences among closely related species. Proceedings of the Royal Society of London, Series B: Biological Sciences 270: 96-99. https://doi.org/10.1098/ rsbl.2003.0025 Hebert PD, Cywinska A, Ball SL, deWaard JR (2003b) Biological identifications through DNA barcodes. Proceedings of the Royal Society of London, Series B: Biological Sciences 270 (1512): 313-321. https://doi.org/10.1098/rspb.2002.2218 Hebert PD, Stoeckle MY, Zemlak TS, Francis CM (2004a) Identification of birds through DNA barcodes. PLOS Biology 2 (10): 1657-1633 (e312). https://doi.org/10.1371/ journal.pbio.0020312 Hebert PD, Penton EH, Burns JM, Janzen DH, Hallwachs W (2004b) Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proceedings of the National Academy of Sciences of the USA 101 (41): 14812-14817. https://doi.org/10.1073/pnas.0406166101 Hebert PD, Gregory TR (2005) The promise of DNA barcoding for taxonomy. Systematic Biology 54 (5): 852-859. https://doi.org/10.1080/10635150500354886 Hey J (2009) On the arbitrary identification of real species. In: Butlin R, Brider J, Schluter D (Eds) Speciation and patterns of diversity. 1st Edition. Cambridge University Press, Cambridge, UK. [ISBN 9780521883184]. https://doi.org/10.1017/ CBO9780511815683.003 Hickerson MJ, Meyer CP, Moritz C (2006) DNA barcoding will often fail to discover new animal species over broad parameter space. Systematic Biology 55 (5): 729-739. https://doi.org/10.1080/10635150600969898 Hubert N, Hanner R, Holm E, Mandrak N, Taylor E, Burridge M, Watkinson D, Dumont P, Curry R, Bentzen P, Zhang J, April J, Bernatchez L (2008) Identifying Canadian freshwater fishes through DNA barcodes. PLOS One 3 (6): 1-8 (€2490). https://doi.org/ 10.1371/journal.pone.0002490 Hubert N, Delrieu-Trottin E, lrisson JO, Meyer C, Planes S (2010) Identifying coral reef fish larvae through DNA barcoding: a test case with the families Acanthuridae and 32 Schtitte A et al Holocentridae. Molecular Phylogenetics and Evolution 55 (3): 1195-1203. https:// doi.org/10.1016/j.ympev.2010.02.023 Hundsdoerfer AK, Rubinoff D, Attie M, Wink M, Kitching IJ (2009) A revised molecular phylogeny of the globally distributed hawkmoth genus Hyles (Lepidoptera: Sphingidae), based on mitochondrial and nuclear DNA sequences. Molecular Phylogenetics and Evolution 52 (3): 852-65. https://doi.org/10.1016/j.ympev.2009.05.023 Irfanullah H (2006) Algal taxonomy in limnology: An example of the declining trend of taxonomic studies? Hydrobiologia 559 (1): 1-9. https://doi.org/10.1007/ $10750-005-9202-z Isaac NJ, Mallet J, Mace GM (2004) Taxonomic inflation: Its influence on macroecology and conservation. Trends in Ecology & Evolution 19 (9): 464-469. https://doi.org/ 10.1016/j.tree.2004.06.004 Jorger K, Schrédl M (2013) How to describe a cryptic species? Practical challenges of molecular taxonomy. Frontiers in Zoology 10 (59): 1-27. https://doi.org/ 10.1186/1742-9994-10-59 Jukes TH, Cantor CR (1969) Evolution of protein molecules (chapter 24). In: Munro HN (Ed.) Mammalian Protein Metabolism. 1st Edition. Academic Press, New York. https:// doi.org/10.1016/B978-1-4832-3211-9.50009-7 Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Mentjies P, Drummond A (2012) Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28 (12): 1647-1649. https:// doi.org/10.1093/bioinformatics/bts199 Kerr KC, Stoeckle MY, Dove CJ, Weigt LA, Francis CM, Hebert PD (2007) Comprehensive DNA barcode coverage of North American birds. Molecular Ecology Notes 7 (4): 535-543. https://doi.org/10.1111/.1471-8286.2007.01670.x Kratky J (2015) Aphytobius veronicae (Frivaldszky, 1884) (Ceutorhynchinae, Hypurini) — species status revised. SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 16 (245): 1-9. URL: https://Awww.curci.de/?beitrag=245 Lanave C, Preparata G, Saccone C, Serio G (1984) Anew method for calculating evolutionary substitution rates. Journal of Molecular Evolution 10 (1): 86-3. https:// doi.org/10.1007/bf02101990 Leather SR, Day KR, Salisbury AN (1999) The biology and ecology of the large pine weevil, Hylobius abietis (Coleoptera: Curculionidae): A problem of dispersal? Bulletin of Entomological Research 89 (1): 3-16. https://doi.org/10.1017/S0007485399000024 Lefebure T, Douady CJ, Gouy M, Gibert J (2006) Relationship between morphological taxonomy and molecular divergence within Crustacea: proposal of a molecular threshold to help species delimitation. Molecular Phylogenetics and Evolution 40 (2): 435-447. https://doi.org/10.1016/j.ympev.2006.03.014 Letsch H, Gottsberger B, Metzl C, Astrin J, Friedman ALL, McKenna DD, Fiedler K (2018) Climate and host-plant associations shaped the evolution of ceutorhynch weevils throughout the Cenozoic. Evolution; International Journal of Organic Evolution 72 (9): 1815-1828. https://doi.org/10.1111/evo.13520 L6bI L, Smetana A (2011) Catalogue of the Coleoptera. Curculionoidea |. 1st Edition, Vol. 7. Apollo Books, Stenstrup, 373 pp. [ISBN 9788788757934] L6bI L, Smetana A (2013) Catalogue of the Coleoptera. Curculionoidea II. 1st Edition, Vol. 8. Brill, Leiden & Boston, 700 pp. [ISBN 9789004252066] Molecular Weevil Identification Project: A thoroughly curated barcode release ... 33 Machado A, Rodriguez-Exposito E, Lopez M, Hernandez M (2017) Phylogenetic analysis of the genus Laparocerus, with comments on colonisation and diversification in Macaronesia (Coleoptera, Curculionidae, Entiminae). Zookeys 651: 1-77. https:// doi.org/10.3897/zookeys.651.10097 Madden T (2002) The BLAST sequence analysis tool (chapter 13). In: McEntyre J, Ostell J (Eds) The NCBI handbook. 2nd Edition. National Center for Biotechnology Information (US), Bethesda, Maryland, 13 pp. URL: http:/Awww.ncbi.nim.nih.gov/books/ NBK21097/ Mayer C, Dietz L, Call E, Kukowka S, Martin S, Espeland M (2021) Adding leaves to the Lepidoptera tree: Capturing hundreds of nuclear genes from old museum specimens. Systematic Entomology 46 (3): 1-23. https://doi.org/10.1111/syen.12481 Mayr E (1942) Systematics and the origin of species from the viewpoint of a zoologist. 1st Edition. Columbia University Press, New York, 334 pp. Meierotto S, Sharkey MJ, Janzen DH, Hallwachs W, Hebert PDN, Chapman EG, Smith MA (2019) A revolutionary protocol to describe understudied hyperdiverse taxa and overcome the taxonomic impediment. Deutsche Entomologische Zeitschrift 66 (2): 119-145. https://doi.org/10.3897/dez.66.34683 Meier R, Shiyang K, Vaidya G, Ng PL (2006) DNA barcoding and taxonomy in Diptera: A tale of high intraspecific variability and low identification success. Systematic Biology 55 (5): 715-728. https://doi.org/10.1080/10635150600969864 Mendelson TC (2003) Evidence of intermediate and asymmetrical behavioral isolation between orangethroat and orangebelly darters (Teleostei: Percidae). American Midland Naturalist 150 (2): 343-347. https://doi.org/ 10.1674/0003-0031 (2003)150[0343:EOIAAB]2.0.CO;2 Meyer CP, Paulay G (2005) DNA barcoding: Error rates based on comprehensive sampling. PLOS Biology 3 (12): 2229-2238 (e422). htips://doi.org/10.1371/journal.pbio. 0030422 Mikkelsen N, Schander C, Willassen E (2007) Local scale DNA barcoding of bivalves (Mollusca): A case study. Zoologica Scripta 36: 455-463. httos://doi.org/10.1111/j. 1463-6409.2006.00289.x Milton M, Pierossi P, Ratnasingham S (2013) Barcode of Life Data Systems Handbook, version 3.6. BOLD Systems, Biodiversity Institute of Ontario, Guelph, Ontario, Canada. URL: https://v3.boldsystems.org/libhtml_v3/static/BOLD Handbook Oct2013.pdf Miralles A, Geniez P, Beddek M, Aranda DM, Brito JC, Leblois R, Crochet PA (2020) Morphology and multilocus phylogeny of the spiny-footed lizard (Acanthodactylus erythrurus) complex reveal two new mountain species from the Moroccan Atlas. Zootaxa 4747 (2): 302-326. https://doi.org/10.11646/zootaxa.4747.2.4 Monaghan M, Balke M, Gregory TR, Vogler A (2005) DNA-based species delineation in tropical beetles using mitochondrial and nuclear markers. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 360: 1925-1933. https:// doi.org/10.1098/rstb.2005.1724 Morando M, Avila LJ, Sites JWJ (2003) Sampling strategies for delimiting species: Genes, individuals, and populations in the Liolaemus elongatus-kriegi complex (Squamata: Liolaemidae) in Andean-Patagonian South America. Systematic Biology 52 (2): 159-85. https://doi.org/10.1080/10635150390192717 Moritz C, Cicero C (2004) DNA barcoding: Promise and pitfalls. PLOS Biology 2 (10): 1529-1531 (e354). https://doi.org/10.1371/journal.pbio.0020354 34 Schutte A et al Morris M, Barclay M (2015) On the nomenclature of Microplontus melanostigma (Marsham, 1802), the weevil formerly referred to as Microplontus rugulosus (Herbst, 1795) (Curculionidae) and its synonyms, with lectotype and neotype designations. The Coleopterist 24: 133-139. Morris MG (2011) The Apionidae (Coleoptera) of the Canary Islands, with particular reference to the contribution of T. Vernon Wollaston. Acta Entomologica Musei Nationalis Pragae 51: 157-182. URL: http:/Awww.aemnp.eu/PDF/51_ 1/51 1 157.pdf Mulcahy DG, Ibanez R, Jaramillo CA, Crawford AJ, Ray JM, Gotte SW, Jacobs JF, Wynn AH, Gonzalez-Porter GP, McDiarmid RW, Crombie RI, Zug GR, de Queiroz K (2022) DNA barcoding of the National Museum of Natural History reptile tissue holdings raises concerns about the use of natural history collections and the responsibilities of scientists in the molecular age. PLOS One 17 (3): 1-23 (e€0264930). https://doi.org/ 10.1371/journal.pone.0264930 Oberprieler RG, Marvaldi AE, Anderson RS (2007) Weevils, weevils, weevils everywhere. Zootaxa (1668)491-520. https://doi.org/10.11646/zootaxa.1668.1.24 Padial JM, De La Riva | (2010) A response to recent proposals for integrative taxonomy. Biological Journal of the Linnean Society 101 (3): 747-756. https://doi.org/10.1111/j. 1095-8312.2010.01528.x Page TJ, Choy SC, Hughes JM (2005) The taxonomic feedback loop: Symbiosis of morphology and molecules. Biology Letters 1 (2): 139-142. https://doi.org/10.1098/rsbl. 2005.0298 Pante E, Schoelinck C, Puillandre N (2015) From integrative taxonomy to species description: One step beyond. Systematic Biology 64 (1): 152-160. https://doi.org/ 10.1093/sysbio/syu083 Pentinsaari M, Ratnasingham S, Miller SE, Hebert PDN (2020) BOLD and GenBank revisited — do identification errors arise in the lab or in the sequence libraries? PLOS One 15 (4): 1-10 (e€0231814). https://doi.org/10.1371/journal.pone.0231814 Pleijel F, Jondelius U, Norlinder E, Nygren A, Oxelman B, Schander C, Sundberg P, Thollesson M (2008) Phylogenies without roots? A plea for the use of vouchers in molecular phylogenetic studies. Molecular Phylogenetics and Evolution 48 (1): 369-371. https://doi.org/10.1016/j.ympev.2008.03.024 Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun §S, Sumlin WD, Vogler AP (2006) Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Systematic Biology 55 (4): 595-609. https://doi.org/ 10.1080/10635150600852011 Posada D (2008) jModelTest: phylogenetic model averaging. Molecular Biology and Evolution 25 (7): 1253-1256. https://doi.org/10.1093/molbev/msn083 Presgraves DC (2002) Patterns of postzygotic isolation in Lepidoptera. Evolution 56 (6): 1168-1183. https://doi.org/10.1111/j.0014-3820.2002.tb01430.x Prosser S, deWaard J, Miller S, Hebert PDN (2015) DNA barcodes from century-old type specimens using next-generation sequencing. Genome 58 (5): 267-268. https:// doi.org/10.1111/1755-0998.12474 Puillandre N, Lambert A, Broillet S, Achaz G (2012) ABGD, automatic barcode gap discovery for primary species delimitation. Molecular Ecology 21 (8): 1864-1877. hitps:// doi.org/10.1111/j.1365-294X.2011.05239.x Molecular Weevil Identification Project: A thoroughly curated barcode release ... 35 Puillandre N, Brouillet S, Achaz G (2021) ASAP: Assemble Species by Automatic Partitioning. Molecular Ecology Resources 21 (2): 609-620. https://doi.org/ 10.1111/1755-0998.13281 Ratnasingham S, Hebert PD (2007) BOLD: The Barcode of Life Data System (http:// www.barcodinglife.org). Molecular Ecology Notes 7 (3): 355-364. hitps://doi.org/10.1111/ j.1471-8286.2007.01678.x Raxworthy CJ, Smith BT (2021) Mining museums for historical DNA: Advances and challenges in museomics. Trends in Ecology & Evolution 36 (11): 1049-1060. https:// doi.org/10.1016/j.tree.2021.07.009 Riedel A, Sagata K, Suhardjono YR, Tanzler R, Balke M (2013a) Integrative taxonomy on the fast track - towards more sustainability in biodiversity research. Frontiers in Zoology 10 (15): 1-9. https://doi.org/10.1186/1742-9994-10-15 Riedel A, Sagata K, Surbakti S, Rene T, Michael B (2013b) One hundred and one new species of Trigonopterus weevils from New Guinea. Zookeys 280: 1-150. https://doi.org/ 10.3897/zookeys.280.3906 Robinson EA, Blagoev G, Hebert PN, Adamowicz S (2009) Prospects for using DNA barcoding to identify spiders in species-rich genera. Zookeys 16: 27-46. https://doi.org/ 10.3897/zookeys.16.239 Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19 (12): 1572-1574. https://doi.org/10.1093/ bioinformatics/btq180 Roycroft E, Moritz C, Rowe KC, Moussalli A, Eldridge MD, Portela Miguez R, Piggott MP, Potter S (2022) Sequence capture from historical museum specimens: maximizing value for population and phylogenomic studies. Frontiers in Ecology and Evolution 10 (931644): 1-12. https://doi.org/10.3389/fevo.2022.931644 Russell MI, Velazquez de Castro AJ (2015) A revision of the genus Pseudoprotapion Ehret, 1990 in the Iberian Peninsula, with description of a new species. Boletin de la Sociedad Entomologica Aragonesa 57: 1-18. URL: http://www.sea-socios.com/ Boletines/Boletin57.html Rylands A, Mittermeier R (2014) Primate taxonomy: Species and conservation. Evolutionary Anthropology: Issues, News, and Reviews 23 (1): 8-10. https://doi.org/ 10.1002/evan.21387 Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular Biology and Evolution 4 (4): 406-425. https://doi.org/ 10.1093/oxfordjournals.molbev.a040454 Sasa MM, Chippindale PT, Johnson NA (1998) Patterns of postzygotic isolation in frogs. Evolution 52 (6): 1811-1820. https://doi.org/10.1111/j.1558-5646.1998 .tb02258.x Scherz MD, Rasolonjatovo SM, Kohler J, Rancilhac L, Rakotoarison A, Raselimanana AP, Ohler A, Preick M, Hofreiter M, Glaw F, Vences M (2020) ‘Barcode fishing’ for archival DNA from historical type material overcomes taxonomic hurdles, enabling the description of a new frog species. Scientific Reports 10 (19109): 1-17. https://doi.org/ 10.1038/s41598-020-75431-9 Schitte A, Stuben PE, Sprick P (2013) The molecular weevil identification project (Coleoptera: Curculionoidea), part | — a contribution to integrative taxonomy and phylogenetic systematics. SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 14 (211): 1-77. URL: https://curci.de/institute/index.php?beitrag=211 36 Schutte A et al Schitte A, Stuben PE (2015) Molecular systematics and morphological identification of the cryptic species of the genus Acal/les Schoenherr, 1825, with descriptions of new species (Coleoptera: Curculionidae: Cryptorhynchinae). Zootaxa 3915 (1): 1-51. https:// doi.org/10.11646/zootaxa.3915.1.1 Schwarz G (1978) Estimating the dimension of a model. The Annals of Statistics 6 (2): 461-464. https://doi.org/10.1214/aos/11 76344136 Sharkey MJ, Janzen DH, Hallwachs W, Chapman EG, Smith MA, Dapkey T, Brown A, Ratnasingham S, Naik S, Manjunath R, Perez K, Milton M, Hebert P, Shaw SR, Kittel RN, Solis MA, Metz MA, Goldstein PZ, Brown JW, Quicke DLJ, van Achterberg C, Brown BV, Burns JM (2021) Minimalist revision and description of 403 new species in 11 subfamilies of Costa Rican braconid parasitoid wasps, including host records for 219 species. Zookeys 1013: 1-665. https://doi.org/10.3897/zookeys.1013.55600 Shea CP, Peterson JT, Wisniewski JM, Johnson NA (2011) Misidentification of freshwater mussel species (Bivalvia:Unionidae): contributing factors, management implications, and potential solutions. Journal of the North American Benthological Society 30 (2): 446-458. https://doi.org/10.1899/10-073.1 Skuhrovec J (2009) Digital-weevil-determination for Curculionoidea of West Palaearctic. Transalpina: Hypera / Limobius / Metadonus (Hyperinae: Hyperini). SNUDEBILLER: Studies on Taxonomy, biology and Ecology of Curculionoidea 10 (121): 39-47. URL: http:/Awww.curci.de/?article=247 Slipinski SA, Leschen RAB, Lawrence JF (2011) Order Coleoptera Linnaeus, 1758. In: Zhang, Z.-Q. (Ed.) animal biodiversity: an outline of higher-level classification and survey of taxonomic richness. Zootaxa 3148 (1): 203-208. https://doi.org/10.11646/ zootaxa.3148.1.39 Sneath PH, Sokal RR (1973) Numerical taxonomy: The principles and practice of numerical classification. 1st Edition. Freeman, W.H. & Co Ltd, San Francisco. [ISBN 9780716706977] Song H, Buhay JE, Whiting MF, Crandall KA (2008) Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proceedings of the National Academy of Sciences of the USA 105 (36): 13486-13491. https://doi.org/10.1073/pnas.0803076105 Spasojevic T, Kropf C, Nentwig W, Lasut L (2016) Combining morphology, DNA sequences, and morphometrics: Revising closely related species in the orb-weaving spider genus Araniella (Araneae, Araneidae). Zootaxa 4111 (4): 448-470. https://doi.org/ 10.11646/zootaxa.4111.4.6 Speidel W, Hausmann A, Muller GC, Kravchenko V, Mooser J, Witt TJ, Khallaayoune K, Prosser S, Hebert PD (2015) Taxonomy 2.0: Sequencing of old type specimens supports the description of two new species of the Lasiocampa decolorata group from Morocco (Lepidoptera, Lasiocampidae). Zootaxa 3999 (3): 401-412. https://doi.org/ 10.11646/zootaxa.3999.3.5 Sperling FA (2003) DNA barcoding: Deus ex machina. Newsletter of the Biological Survey of Canada (Terrestrial Arthropods) 22 (1): 50-53. Sprick P (2009) Monitoring von Russelkafern in Baumschulen, Staudengartnereien und Hopfengarten — Ergebnisse des ersten Untersuchungsjahres (2008). Deutsche Gesellschaft fur Allgemeine und Angewandte Entomologie 17. Sprick P (2019) Supplement to the Apionidae keys I-V (DWD: digital-weevil- determination for Curculionoidea of the West Palaearctic: Transalpina): distribution, host Molecular Weevil Identification Project: A thoroughly curated barcode release ... 37 plants, biology. / Nachtrag zu den Apionidae-Schlusseln I-V (DWD: Digital-Weevil- Determination der westpalaarktischen Curculionoidea: Transalpina): Verbreitung, Wirtspflanzen, Biologie. SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 20 (280): 1-29. URL: https://www.curci.de/?beitrag=280 Strutzenberger P, Brehm G, Fiedler K (2012) DNA barcode sequencing from old type specimens as a tool in taxonomy: A case study in the diverse genus Eo/s (Lepidoptera: Geometridae). PLOS One 7 (11): e49710. https://doi.org/10.1371/journal.pone.0049710 Stuben PE (1994) Das Trachyphloeus bifoveolatus-angustisetulus-Problem und die intermediaren Exemplare des Rheinlandes: Gibt es eine morphologische L6sung? (Coleoptera, Curculionidae). Mitteilungen der Arbeitsgemeinschaft Rheinischer Koleopterologen (Bonn) 4 (3): 111-129. Stuben PE (2005) "Basar Taxonomie"? - Ein erfolgreiches Kreuzungsexperiment zu Acalles aeonii Wollaston, 1864 (Coleoptera: Curculionidae). Weevil News 31: 1-13. URL: https:/Awww.curci.de/?beitrag=144 Stuben PE, Germann C (2005) Neue Erkenntnisse zur Taxonomie, Biologie und Okologie der Cryptorhynchinae von den Makaronesischen Inseln. 1. Beitrag: Kanaren / Tenerife (Coleoptera: Curculionidae: Cryptorhynchinae). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 6 (61): 37-83. URL: https:// www.curci.de/?beitrag=61 Stuben PE (2010) Calacalles (Nanoacalles) agana sp.n. (Coleoptera: Curculionidae: Cryptorhynchinae). Weevil News 58: 1-2. URL: https:/Awww.curci.de/?beitrag=167 Stuben PE, Astrin JJ (2010) Molecular phylogeny in endemic weevils: Revision of the genera of Macaronesian Cryptorhynchinae (Coleoptera: Curculionidae). Zoological Journal of the Linnean Society 160 (1): 40-87. https://doi.org/10.1111/j. 1096-3642.2009.00609.x Stuben PE, Behne L (2010) Hemitrichapion gomerense sp.n. (Coleoptera: Curculionoidea: Apionidae). Weevil News 53: 1-3. URL: https://www.curci.de/? beitrag=162 Stuben PE, Astrin JJ (2011) Aeoniacalles tabladoensis sp.n. and Aeoniacalles aeonii bodegensis (Stuben 2000) resyn. from the Canary Islands (Coleoptera: Curculionidae: Cryptorhynchinae). Weevil News 68: 1-5. URL: https:/Awww.curci.de/?beitrag=186 Stuben PE, Sprick P, Miller G, Bayer C, Behne L, Kratky J (2012) Digital-weevil- determination for Curculionoidea of West Palaearctic: Transalpina: Ceutorhynchinae (1. part). (Mononychini, Phytobiini, Hypurini, Cnemogonini, Scleropterini & Amalini). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 13 (191): 18-33. URL: https:/Awww.curci.de/?beitrag=191 Stuben PE, Schutte A (2013) Die Wiederentdeckung von Acalles xerampelinus Wollaston, 1864, und die Stellung im phylogenetischen System der kanarischen Cryptorhynchinae (Curculionidae). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 14 (218): 1-14. URL: https:/Awww.curci.de/?beitrag=218 Stuben PE, Schutte A, Astrin JJ (2013a) Molecular phylogeny of the weevil genus Dichromacalles Stuben (Curculionidae: Cryptorhynchinae) and description of a new species. Zotaxa 3718 (2): 101-127. https://doi.org/10.11646/zootaxa.3718.2.1 Stiiben PE, Miller G, Kratky J, Bayer C, Behne L, Sprick P (2013b) Digital-weevil- determination der westpalaarktischen Curculionoidea: Transalpina: Ceutorhynchinae (2. Part). (Ceuthorhynchini: Amalorrhynchus, Drupenatus, Poophagus, Coeliodes, Pseudocoeliodes, Coeliodinus, Eucoeliodes, Neoxyonyx, Thamiocolus, Micrelus, 38 Schutte A et al Zacladus, Phrydiuchus, Stenocarus, Nedyus, Ceutorhynchus: Marklissus). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 14 (210): 1-23. URL: https:/Awww.curci.de/?beitrag=210 Stuben PE (2014a) Die Curculionoidea (Coleoptera) Tenerifes. SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 15 (226): 1-5. URL: https:// www.curci.de/?beitrag=226 Stuben PE (2014b) New nomenclatural and taxonomic acts, and comments (2014). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 15 (231): 1-5. URL: https:/Avww.curci.de/?beitrag=231 Stuben PE (2015) Zwei neue Calacalles (Nanoacalles)-Arten von La Palma (Kanarische Inseln) (Coleoptera: Curculionidae: Cryptorhynchinae). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 16 (243): 1-10. URL: http:// www.curci.de/?beitrag=243 Stuben PE, Bayer C (2015) New nomenclatural and taxonomic acts, and comments (2015). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 16 (246): 1-8. URL: http:/Awww.curci.de/?beitrag=246 Stuben PE, Behne L (2015) Die Curculionoidea (Coleoptera) La Palmas. SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 16 (242): 1-86. URL: https:/Awww.curci.de/?beitrag=242 Stuben PE, Schutte A (2015) Torneuma korwitzi sp. n. from Madeira (Coleoptera: Curculionidae: Cryptorhynchinae). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 16 (36): 1-8. URL: https://www.curci.de/?beitrag=236 Stuben PE, Schutte A, Bayer C, Astrin JJ (2015) The molecular weevil identification project (Coleoptera: Curculionoidea), part Il — towards an integrative taxonomy. SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 16 (237): 1-294. URL: https:/Awww.curci.de/?beitrag=237 Stuben PE, Kratky J (2016) Hesperorrhynchus glutinosus sp.n. and H. palmensis sp.n. — with a review of all Hesperorrhynchus species from the Macaronesian Islands (Coleoptera: Curculionidae: Ceutorhynchinae) — an contribution to integrative taxonomy. SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 17 (247): 1-8. URL: https:/Avww.curci.de/?beitrag=247 Stuben PE, Schutte A, Lopez H, Astrin JJ (2016a) Molecular and morphological systematics of soil-inhabiting Cryptorhynchinae of the genus Acallorneuma and the tribe Torneumatini (Coleoptera: Curculionidae), with description of two new species. Contributions to Entomology 66 (2): 169-199. https://doi.org/10.21248/contrib.entomol. 66.2.169-199 Stuben PE, Bayer C, Bahr F, Sprick P, Behne L (2016b) Digital-weevil-determination for Curculionoidea of the West Palearctic: Transalpina: Anthonomini. SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 17 (249): 1-11. URL: https:/Awww.curci.de/?beitrag=249 Stuben PE (2017a) Die Curculionoidea (Coleoptera) von den Inseln Madeiras und den Selvagens. SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 18 (261): 1-92. URL: https://www.curci.de/?beitrag=261 Stuben PE (2017b) New nomenclatural and taxonomic acts, and comments (2016/2017). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 18 (265): 1-16. URL: https:/Awww.curci.de/?beitrag=265 Molecular Weevil Identification Project: A thoroughly curated barcode release ... 39 Stuben PE (2018a) Beschreibung neuer Arten aus der Gattung Caulotrupis Wollaston, 1854 (Coleoptera: Curculionidae: Cossoninae) von den Makaronesischen Inseln. Contributions to Entomology 68 (1): 83-96. https://doi.org/10.21248/contrib.entomol. 68.1.083-096 Stuben PE (2018b) The Cryptorhynchinae of the Western Palearctic. - Die Cryptorhynchinae der Westpalaarktis (Coleoptera: Curculionidae). 1st edition. Curculio Institute, Ménchengladbach, 518 pp. URL: hitps:/Avww.curci.de/?beitrag=281 [ISBN 9783000580383] Stuben PE (2018c) Kritische Anmerkungen zu den neuen Curculionoidea-Katalogen - nomenklatorische und taxonomische Anderungen (2018) (Coleoptera, Curculionoidea). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 19 (275): 1-13. URL: https:/Awww.curci.de/?beitrag=275 Stuben PE, Schitte A (2018) Revision of the genus Caulotrupis from Macaronesian Islands (Coleoptera: Curculionidae: Cossoninae). SNUDEBILLER: Studies on Taxonomy, Biology and Ecology of Curculionoidea 19 (272): 1-38. URL: hitps:// www.curci.de/?beitrag=272 Stuben PE, Kramp K (2019) Neue Echinodera aus Griechenland — Beitrag zur integrativen Taxonomie (Coleoptera: Curculionidae: Cryptorhynchinae). Contributions to Entomology 69 (2): 319-330. https://doi.org/10.21248/contrib.entomol.69.2.319-330 Stuben PE, Schutte A, Astrin JJ (2021) Barcoding and interspecific relationships of Macaronesian weevils (Coleoptera: Curculionoidea). Contributions to Entomology 71 (1): 127-135. https://doi.org/10.21248/contrib.entomol.71.1.127-135 Stuben PE (2022) Weevils of Macaronesia - Canary Islands, Madeira, Azores (Coleoptera: Curculionoidea). 1st Edition. Curculio Institute, Ménchengladbach, Germany, 784 pp. [ISBN 9783000684166] Tautz D, Arctander P, Minelli A, Thomas R, Vogler A (2002) DNA points the way ahead in taxonomy. Nature 418 (6897): 479-479. https://doi.org/10.1038/418479a Tautz D, Arctander P, Minelli A, Thomas R, Vogler A (2003) A plea for DNA taxonomy. Trends in Ecology & Evolution 18 (2): 70-74. https://doi.org/10.1016/ $0169-5347(02)00041-1 Teletchea F (2010) After 7 years and 1000 citations: Comparative assessment of the DNA barcoding and the DNA taxonomy proposals for taxonomists and non-taxonomists. Mitochondrial DNA 21 (6): 206-26. https://doi.org/10.3109/19401736.2010.532212 Thiem H (1932) Der gefurchte Dickmaulriler (Otiorrhynchus sulcatus F.) als Gewachshaus- und Freilandschadling. Die Gartenbauwissenschaft 6: 519-540. Townson H, Harbach RE, Callan TA (1999) DNA identification of museum specimens of the Anopheles gambiae complex: An evaluation of PCR as a tool for resolving the formal taxonomy of sibling species complexes. Systematic Entomology 24 (1): 95-100. https://doi.org/10.1046/j.1365-3113.1999.00084.x Vences M, Thomas M, Bonett RM, Vieites DR (2005) Deciphering amphibian diversity through DNA barcoding: Chances and challenges. Philosophical transactions of the Royal Society of London. Series B, Biological Sciences 360 (1462): 1859-1868. hitps:// doi.org/10.1098/rstb.2005.1717 Vences M (2020) The promise of next-generation taxonomy. Megataxa 001 (1): 35-38. https://doi.org/10.11646/megataxa.1.1.6 Wagele JW (2005) Foundations of phylogenetic systematics. 1st Edition. Pfeil, Munchen, Germany, 365 pp. [ISBN 9783899370560] 40 Schtitte A et al Warren A, Patterson DJ, Dunthorn M, Clamp JC, Achilles-Day UEM, Aescht E, Al-Farraj SA, Al-Quraishy S, Al-Rasheid K, Carr M, Day JG, Dellinger M, El-Serehy HA, Fan Y, Gao F, Gao S, Gong J, Gupta R, Hu X, Kamra K, Langlois G, Lin X, Lipscomb D, Lobban CS, Luporini P, Lynn DH, Ma H, Macek M, Mackenzie-Dodds J, Makhija S, Mansergh RI, Martin-Cereceda M, McMiller N, Montagnes DJS, Nikolaeva S, Ong'ondo GO, Perez-Uz B, Purushothaman J, Quintela-Alonso P, Rotterova J, Santoferrara L, Shao C, Shen Z, Shi X, Song W, Stoeck T, La Terza A, Vallesi A, Wang M, Weisse T, Wiackowski K, Wu L, Xu K, Yi Z, Zufall R, Agatha S (2017) Beyond the “code”: A guide to the description and documentation of biodiversity in ciliated protists (Alveolata, Ciliophora). The Journal of Eukaryotic Microbiology 64 (4): 539-554. https://doi.org/ 10.1111/jeu.12391 Waugh J (2007) DNA barcoding in animal species: Progress, potential and pitfalls. Bioessays 29 (2): 188-97. https://doi.org/10.1002/bies.20529 Wheeler QD, Raven PH, Wilson EO (2004) Taxonomy: Impediment or expedient? Science 303 (5656): 285-285. https://doi.org/10.1126/science.303.5656.285 Wiemers M, Fiedler K (2007) Does the DNA barcoding gap exist? - a case study in blue butterflies (Lepidoptera: Lycaenidae). Frontiers in Zoology 4 (8): 1-16. https://doi.org/ 10.1186/1742-9994-4-8 Will KW, Mishler BD, Wheeler QD (2005) The perils of DNA barcoding and the need for integrative taxonomy. Systematic Biology 54 (5): 844-851. https://doi.org/ 10.1080/10635150500354878 Wu F, Yan XP (2018) Distribution of the related weevil species Sitophilus oryzae and S. zeamais (Coleoptera: Curculionidae) in farmer stored grains of China. Journal of Economic Entomology 111 (3): 1461-1468. httos://doi.org/10.1093/jee/toy061 Zamani A, Fric ZF, Gante H, Hopkins T, Orfinger A, Scherz M, Bartonova AS, Pos DD (2022) DNA barcodes on their own are not enough to describe a species. Systematic Entomology 47 (3): 385-389. https://doi.org/10.1111/syen.12538 Zigler KS, McCartney MA, Levitan DR, Lessios HA (2005) Sea urchin bindin divergence predicts gamete compatibility. Evolution 59 (11): 2399-2404. https://doi.org/ 10.1554/05-098.1 Zwolfer V, Herbst J (1988) Praadaptation, Wirtskreiserweiterung and Parallel- Cladogenese in der Evolution von phytophagen Insekten *1, *2. Journal of Zoological Systematics and Evolutionary Research 26 (5): 320-340. https://doi.org/10.1111/j. 1439-0469.1988.tb00321.x Supplementary materials Suppl. material 1: Material Table EE Authors: Schutte A, Stuben PE, Astrin JJ Data type: spreadsheet with collecting data, voucher numbers and GenBank acc. numbers Brief description: The material table (248 pages in PDF format) contains information about each specimen's collecting spot, GPS position, collector, identifier, voucher numbers (DNA and tissue) and GenBank accession numbers. Download file (3.23 MB) Molecular Weevil Identification Project: A thoroughly curated barcode release ... 41 Suppl. material 2: CO1 Sequences FEI Authors: Schutte A, Astrin JJ Data type: CO1 sequences, DNA barcodes, nucleotide alignment, list with GenBank accession numbers Brief description: This study's complete DNA barcode dataset is provided as nucleotide alignment and unaligned sequences in *.fasta file format. The *.fasta files can be opened with any text editor. A read_me.txt file is included explaining the naming scheme. A list of GenBank accession numbers only is included. Download file (640.59 kb) Suppl. material 3: NJ Tree El Authors: Schitte A, Astrin JJ Data type: MW neighbour-joining tree Brief description: The ZIP file contains the neighbour-joining tree in four different file formats (newick, nex, png, svg). The tree is based on the complete DNA barcode dataset published in this study. Download file (12.38 MB) Suppl. material 4: Results of Taxonomy EE Authors: Schitte A, Stuben PE Data type: morphological identification versus molecular results Brief description: Contradictions between morphological identifications and molecular results are discussed (PDF file). Download file (104.23 kb) Suppl. material 5: Bayesian Trees EE Authors: Schutte A Data type: Bayesian Inference calculation, 50% majority rule consensus trees, Cryptorhynchinae, Apioninae, Ceutorhynchinae Brief description: The ZIP file contains the 50% majority rule consensus trees for three sub- datasets of the study, calculated with MrBayes: 1) Cryptorhynchinae + Cossoninae; 2) Apioninae + Nanophyinae + Attelabidae; 3) Ceutorhynchinae. Each tree is provided in three different file formats (NEX, PNG, SVG). The searchable SVG file can be opened with any internet browser. Download file (6.99 MB) 42 Schitte A et al Suppl. material 6: DiStats statistics FE Authors: Schitte A Data type: DiStats method description in detail (PDF), DiStats input data (nucleotide alignments) and output data (TXT), spreadsheets (XLSX) with data compilation Brief description: The supplement contains an in-depth description of the used method (which taxon was used as a reference taxon to create the morphologically calibrated p-distance statistics for Cryptorhynchinae, Apioninae and Ceutorhynchinae). Besides the input data to the DiStats scripts and unformatted output data, the data compilation leading to the output tables is described and available in three Excel spreadsheets. Download file (628.69 kb) Suppl. material 7: ASAP analyses EE Authors: Schutte A Data type: nucleotide alignments (FASTA), raw data output (HTML, TXT), Excel spreadsheets (XLSX) Brief description: Detailed description of the ASAP program and data compilation leading to the final summarised results table. Download file (2.49 MB) Suppl. material 8: Thresholds in various taxa EJ Authors: Schitte A Data type: text Brief description: Excursus about various taxa and their thresholds for species delineation. Download file (61.12 kb) Suppl. material 9: Aeoniacalles aeonii bodegensis [EI Authors: Schitte A, Stuben PE, Astrin, JJ Data type: text Brief description: Additional information to the unsolved taxonomic status of Aeoniacalles aeonii bodegensis (Stuben, 2000). Download file (45.36 kb)