DNA Barcoding Southwestern Atlantic Skates: A 20-Year Effort in Building a Species Identification Library

The skate fauna in the Southwest Atlantic Ocean (SWA; 34–55° S) is represented by ~32 species, many of which share external features that have led to misidentifications and deficient fishery statistics. The use of DNA barcoding to discriminate SWA skate species was explored after 20 years of surveys. COI sequences were subjected to distance-based neighbor-joining (NJ), maximum likelihood (ML), barcode index number (BIN), automatic barcode gap discovery (ABGD), and nucleotide diagnostic character (NDC) analyses. For widely distributed species, a haplotype network was built. Overall, 187 specimens and 31 egg cases from 26 skate species were barcoded. NJ and ML analyses showed that nearly all species exhibited unique barcodes or clusters of closely related haplotypes, except for Psammobatis normani/P. rudis and Dipturus trachyderma/D. argentinensis. The first pair was discriminated by NCD. BIN analysis recovered 17 groups, whereas ABGD recovered 23, better reflecting taxonomic diversity. In summary, 24 species were resolved by COI. Phylogeographic signals were observed for Amblyraja doellojuradoi and Zearaja brevicaudata. Compiling our results with data from BOLD, almost all the species occurring in the area possess barcodes, contributing to completing and curating the BOLD reference library, which constitutes an important tool for resolving taxonomic issues, tracing fishery products, and performing eDNA biomonitoring.