Morphology, nuclear SNPs and mate selection reveal that COI barcoding overestimates species diversity in a Mediterranean freshwater amphipod by an order of magnitude

DNA sequence information has revealed many morphologically cryptic species worldwide. For animals, DNA-based assessments of species diversity usually rely on the mitochondrial cytochrome c oxidase subunit I (COI) gene. However, a growing amount of evidence indicate that mitochondrial markers alone can lead to misleading species diversity estimates due to mito-nuclear discordance. Therefore, reports of putative species based solely on mitochondrial DNA should be verified by other methods, especially in cases where COI sequences are identical for different morphospecies or where divergence within the same morphospecies is high. Freshwater amphipods are particularly interesting in this context because numerous putative cryptic species have been reported. Here, we investigated the species status of the numerous mitochondrial molecular operational taxonomic units (MOTUs) found within Echinogammarus sicilianus. We used an integrative approach combining DNA barcoding with mate selection observations, detailed morphometrics and genome-wide double digest restriction site-associated DNA sequencing (ddRAD-seq). Within a relatively small sampling area, we detected twelve COI MOTUs (divergence = 1.8–20.3%), co-occurring in syntopy at two-thirds of the investigated sites. We found that pair formation was random and there was extensive nuclear gene flow among the ten MOTUs co-occurring within the same river stretch. The four most common MOTUs were also indistinguishable with respect to functional morphology. Therefore, the evidence best fits the hypothesis of a single, yet genetically diverse, species within the main river system. The only two MOTUs sampled outside the focal area were genetically distinct at the nuclear level and may represent distinct species. Our study reveals that COI-based species delimitation can significantly overestimate species diversity, highlighting the importance of integrative taxonomy for species validation, especially in hyperdiverse complexes with syntopically occurring mitochondrial MOTUs.