Quantifying Ecological, Morphological, and Genetic Variation to Delimit Species in the Coast Horned Lizard Species Complex (Phrynosoma)

Lineage separation and divergence form a temporally extended process whereby populations may diverge genetically, morphologically, or ecologically, and these contingent properties of species provide the operational criteria necessary for species delimitation. We inferred the historical process of lineage formation in the coast horned lizard (Phrynosoma coronatum) species complex by evaluating a diversity of operational species criteria, including divergence in mtDNA (98 specimens; 2,781 bp) and nuclear loci (RAG-1, 1,054 bp; BDNF 529 bp), ecological niches (11 bioclimatic variables; 285 unique localities), and cranial horn shapes (493 specimens; 16 landmarks). A phylogenetic analysis of mtDNA recovers 5 phylogeographic groups arranged latitudinally along the Baja California Peninsula and in California. The 2 southern phylogeographic groups exhibit concordance between genetic, morphological, and ecological divergence; however, differentiation is weak or absent at more recent levels defined by phylogeographic breaks in California. Interpreting these operational species criteria together suggests that there are 3 ecologically divergent and morphologically diagnosable species within the P. coronatum complex. Our 3-species taxonomic hypothesis invokes a deep coalescence event when fitting the mtDNA genealogy into the species tree, which is not unexpected for populations that have diverged recently. Although the hypothesis that the 3 phylogeographic groups distributed across California each represent distinctive species is not supported by all of the operational species criteria evaluated in this study, the conservation status of the imperiled populations represented by these genealogical units remains critical.