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 data sets 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 data set to 1,294 DNA barcodes. A multi-approach species delimitation analysis involving two distance (OT and ABGD) and four tree-based (PTP, MPTP, GMYC and bGMYC) methods was 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 into 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.