Joshua D. Stewart, Andrew D. Foote, Helen Hipperson, J. Kaden, Helen Senn, Alessandro Ponzo, Lauren R. Peel, Guy Stevens, Jehad Salah, Jane Hosegood, Rob Ogden, Sabine P. Wintner, Marloes Poortvliet, de Bruyn M, Daniel Fernando, Ramon Bonfil, Mohammed Abudaya, Gary R. Carvalho, Emily Humble, Simon Creer, Muhammad Moazzam, Stephen Pollett, Kim Bassos-Hull, and Rima W. Jabado
Practical biodiversity conservation relies on delineation of biologically meaningful units, particularly with respect to global conventions and regulatory frameworks. Traditional approaches have typically relied on morphological observation, resulting in artificially broad delineations and non-optimal species units for conservation. More recently, species delimitation methods have been revolutionised with High-Throughput Sequencing approaches, allowing study of diversity within species radiations using genome-wide data. The highly mobile elasmobranchs, manta and devil rays (Mobulaspp.), are threatened globally by targeted and bycatch fishing pressures resulting in recent protection under several global conventions. However, a lack of global data, morphological similarities, a succession of recent taxonomic changes and ineffectual traceability measures combine to impede development and implementation of a coherent and enforceable conservation strategy. Here, we generate genome-wide Single Nucleotide Polymorphism (SNP) data from among the most globally and taxonomically representative set of mobulid tissues. The resulting phylogeny and delimitation of species units represents the most comprehensive assessment of mobulid diversity with molecular data to date. We find a mismatch between current species classifications, and optimal species units for effective conservation. Specifically, we find robust evidence for an undescribed species of manta ray in the Gulf of Mexico and show that species recently synonymised are reproductively isolated. Further resolution is achieved at the population level, where cryptic diversity is detected in geographically distinct populations, and indicates potential for future traceability work determining regional location of catch. We estimate the optimal species tree and uncover substantial incomplete lineage sorting, where standing variation in extinct ancestral populations is identified as a driver of phylogenetic uncertainty, with further conservation implications. Our study provides a framework for molecular genetic species delimitation that is relevant to wide-ranging taxa of conservation concern, and highlights the potential for genomic data to support effective management, conservation and law enforcement strategies.