Phylogeny of snakes (Serpentes): combining morphological and molecular data in likelihood, Bayesian and parsimony analyses

The phylogeny of living and fossil snakes is assessed using likelihood and parsimony approaches and a dataset combining 263 morphological characters with mitochondrial (2693 bp) and nuclear (1092 bp) gene sequences. The 'no common mechanism' (NCMr) and 'Markovian' (Mkv) models were employed for the morphological partition in likelihood analyses; likelihood scores in the NCMr model were more closely correlated with parsimony tree lengths. Both models accorded relatively less weight to the molecular data than did parsimony, with the effect being milder in the NCMr model. Partitioned branch and likelihood support values indicate that the mtDNA and nuclear gene partitions agree more closely with each other than with morphology. Despite differences between data partitions in phylogenetic signal, analytic models, and relative weighting, the parsimony and likelihood analyses at retrieved the following widely accepted groups: scolecophidians, alethinophidians, cylindrophiines, macrostomatans (sensu lato) and caenophidians. Anilius alone emerged as the most basal alethinophidian; the combined analyses resulted in a novel and stable position of uropeltines and cylindrophiines as the second-most basal clade of alethinophidians. The limbed marine pachyophiids, along with Dinilysia and Wonambi, were always basal to all living snakes. Other results stable in all combined analyses include: Xenopeltis and Loxocemus were sister taxa (fide morphology) but clustered with pythonines (fide molecules), and Ungaliophis clustered with a boine-erycine clade (fide molecules). Tropidophis remains enigmatic; it emerges as a basal alethinophidian in the parsimony analyses (fide molecules) but a derived form in the likelihood analyses (fide morphology), largely due to the different relative weighting accorded to data partitions.