Italian Odonates in the Pandora's Box: A Comprehensive DNA Barcoding Inventory Shows Taxonomic Warnings at the Holarctic Scale

The Odonata are considered among the most endangered freshwater faunal taxa. Their DNA-based monitoring relies on validated reference datasets that are often lacking or do not cover important biogeographical centres of diversification. This study presents the results of a DNA barcoding campaign on Odonata, based on the standard 658 bp 5 end region of the mitochondrial COI gene, involving the collection of 812 specimens (409 of which barcoded) from peninsular Italy and its main islands (328 localities), belonging to all the 88 species (31 Zygoptera and 57 Anisoptera) known from the country. Additional BOLD and GenBank data from Holarctic samples expanded the dataset to 1294 DNA barcodes. A multi-approach species delimitation analysis involving two distance (OT and ABGD) and four tree-based (PTP, MPTP, GMYC, bGMYC) methods were used to explore these data. Of the 88 investigated morphospecies, 75 (85%) unequivocally corresponded to distinct Molecular Operational Units, whereas the remaining ones were classified as warnings (i.e., showing a mismatch between morphospecies assignment and DNA-based species delimitation). These results are in contrast with other DNA barcoding studies on Odonata showing up to 95% of identification success. The species causing warnings were grouped in three categories depending on if they showed low, high, or mixed genetic divergence patterns. The analysis of haplotype networks revealed unexpected intraspecific complexity at the Italian, Palearctic, and Holarctic scale, possibly indicating the occurrence of cryptic species. Overall, this study provides new insights into the taxonomy of odonates and a valuable basis for future DNA and eDNA-based monitoring studies.