Regionalization of the vertebral column and its correlation with heart position in snakes: Implications for evolutionary pathways and morphological diversification

Spinal regionalization has important implications for the evolution of vertebrate body plans. We determined the variation in the number and morphology of vertebrae across the vertebral column (i.e., vertebral formula) for 63 snake species representing 13 families using intracolumnar variation in vertebral shape. Vertebral counts were used to determine the position of the heart, pylorus, and left kidney for each species. Across all species we observed a conspicuous midthoracic transition in vertebral shape, indicating four developmental domains of the precloacal vertebral column (cervical, anterior thoracic, posterior thoracic, and lumbar). Using phylogenetic analyses, the boundary between the anterior and posterior thoracic vertebrae was correlated with heart position. No associations were found between shifts in morphology of the vertebral column and either the pylorus or left kidney. We observed that among taxa, the number of preapex and postapex vertebrae could change independently from one another and from changes in the total number of precloacal vertebrae. Ancestral state reconstruction of the preapex and postapex vertebrae illustrated several evolutionary pathways by which diversity in the vertebral column and heart position have been attained. In addition, no conspicuous pattern was observed among the heart, pylorus, or kidney indicating that their relative positions to each other evolve independently. We conclude that snakes exhibit four morphologically distinct regions of the vertebral column. We discuss the implications of the forebody and hindbody vertebral formula on the morphological diversification of snakes.