Phenotypic regionalization of the vertebral column in the thorny skate Amblyraja radiata: Stability and variation.

Regionalization of the vertebral column occurred early during vertebrate evolution and has been extensively investigated in mammals. However, less data are available on vertebral regions of crown gnathostomes. This is particularly true for batoids (skates, sawfishes, guitarfishes, and rays) whose vertebral column has long been considered to be composed of the same two regions as in teleost fishes despite the presence of a synarcual. However, the numerous vertebral units in chondrichthyans may display a more complex regionalization pattern than previously assumed and the intraspecific variation of such pattern deserves a thorough investigation. In this study, we use micro-computed tomography (µCT) scans of vertebral columns of a growth series of thorny skates Amblyraja radiata to provide the first fine-scale morphological description of vertebral units in a batoids species. We further investigate axial regionalization using a replicable clustering analysis on presence/absence of vertebral elements to decipher the regionalization of the vertebral column of A. radiata. We identify four vertebral regions in this species. The two anteriormost regions, named synarcual and thoracic, may undergo strong developmental or functional constraints because they display stable patterns of shapes and numbers of vertebral units across all growth stages. The third region, named hemal transitional, is characterized by high inter-individual morphological variation and displays a transition between the monospondylous (one centrum per somite) to diplospondylous (two centra per somite) conditions. The posteriormost region, named caudal, is subdivided into three sub-regions with shapes changing gradually along the anteroposterior axis. These regionalized patterns are discussed in light of ecological habits of skates.