Development of molecular and morphological markers to improve species-specific monitoring and systematics of Northeast Atlantic and Mediterranean skates (Rajiformes)

The Northeast Atlantic and Mediterranean skates (Rajidae) showed remarkable species diversity but with high morphological and ecological conservatism. Since skates are particularly vulnerable to the bottom trawl fishery, species-specific demographic surveys as well as studies defining life history and evolutionary traits are important in prioritising conservation programs. However, the identification of juveniles and adults of some species may be difficult using referenced guidelines and identification keys. Therefore, we attempt to develop markers for species identification through the parallel analysis of a 16S rDNA gene sequence and of several morphological characters on 135 individuals collected by trawl surveys in the Adriatic Sea and putatively assigned to six taxa. Species-specific haplotypes were defined for Raja miraletus, Raja montagui, Dipturus oxyrinchus, since a solid accordance between species boundaries and well-differentiated haplotypes was observed. Comparative analysis of 16S rDNA sequences allowed the identification of three juvenile specimens of Leucoraja circularis, a species that rarely occurs in the Adriatic Sea. On the contrary, morphological traits and haplotype distribution were largely discordant in Raja asterias and R. clavata. While all putative R. clavata individuals showed a unique haplotype (H-CLA), only 8 of 30 putative R. asterias individuals possessed a second weakly divergent haplotype (H-AST). The remaining 22 R. asterias carried the H-CLA. The multivariate analyses of morphometric and meristic characters in putative R. clavata and R. asterias revealed the clustering of individuals regardless of haplotypes. However, a bimodal distribution of R. asterias and R. clavata samples would suggest that two separated taxa might exist, both sharing the two 16S rDNA haplotypes. The haplotype distribution appeared to be significantly correlated only to the standardised disc length/total length (DL/TL) variation. Three alternative explanations may support this scenario: (i) an incomplete lineage sorting process in two morphologically yet distinct taxa; (ii) a recent hybridisation between the two taxa; (iii) the two taxa are morphologically plastic species and all considered morphological characters may be misleading in discriminating between them at all maturity stages, except for the DL/TL. However, further analyses on larger data sets and using molecular key markers (i.e. nuclear genes) will be needed to definitely resolve the status of these taxa. Molecular relationships among rajid taxa are largely consistent with systematics based on internal and external anatomical features. This multidisciplinary study contributed to defining the pattern of species diversity and abundance of rajids in the Adriatic Sea. © 2003 Elsevier Science B.V. All rights reserved.