Editorial Type:
Article Category: Research Article
 | 
Online Publication Date: 20 Aug 2015

Phylogenomics of Horned Lizards (Genus: Phrynosoma) Using Targeted Sequence Capture Data

and
Page Range: 586 – 594
DOI: 10.1643/CH-15-248
Save
Download PDF

New genome sequencing techniques are enabling phylogenetic studies to scale-up from using a handful of loci to hundreds or thousands of loci from throughout the genome. In this study, we use targeted sequence capture (TSC) data from 540 ultraconserved elements and 44 protein-coding genes to estimate the phylogenetic relationships among all 17 species of horned lizards in the genus Phrynosoma. Previous molecular phylogenetic analyses of Phrynosoma based on a few nuclear genes, restriction site associated DNA (RAD) sequencing, or mitochondrial DNA (mtDNA) have produced conflicting relationships. Some of these conflicts are likely the result of rapid speciation at the start of Phrynosoma diversification, whereas other examples of gene tree discordance appear to be caused by active and residual traces of hybridization. Concatenation and coalescent-based species tree phylogenetic analyses of these new TSC data support the same topology, and a divergence dating analysis suggests that the Phrynosoma crown group is up to 30 million years old. The new phylogenomic tree supports the recognition of four main clades within Phrynosoma, including Anota (P. mcallii, P. solare, and the P. coronatum complex), Doliosaurus (P. modestum, P. goodei, and P. platyrhinos), Tapaja (P. ditmarsi, P. douglasii, P. hernandesi, and P. orbiculare), and Brevicauda (P. braconnieri, P. sherbrookei, and P. taurus). The phylogeny provides strong support for the relationships among all species of Phrynosoma and provides a robust new framework for conducting comparative analyses.

Copyright: © 2015 by the American Society of Ichthyologists and Herpetologists
Fig. 1.
Fig. 1.

Time-calibrated phylogeny for Phrynosoma and sand lizards based on a Bayesian relaxed-clock analysis of 584 nuclear loci. The dual error bars on each node show the 95% highest posterior density (HPD) for divergence times from the incomplete (top; 584 loci) and complete (bottom; 223 loci) data sets. The support values that correspond to each lettered branch are provided in Table 3.


Contributor Notes

Associate Editor: B. Stuart.

Received: 13 Feb 2015
Accepted: 07 May 2015
  • Download PDF