Understanding genomic landscapes of differentiation in round-tailed horned lizards ( Phrynosoma modestum )

Population divergence is promoted by divergent selection and inhibited by gene flow, but the mechanisms of and relationship between these two processes remain poorly understood. Developing a well-informed hypothesis of the selective pressures underlying divergence in a natural population requires a thorough understanding of both species structure and demographic history. In this study, we assess whole-genome sequences of round-tailed horned lizards (Phrynosoma modestum) from throughout the species range and combine phylogenetic analyses with genomic landscape scans to understand how current genetic diversity has been influenced by demographic histories and evolutionary pressures. Maximum likelihood phylogenetic analysis supports two lineages within the species, corresponding to a North/South population divide that developed around 7 Ma and displays little migration. However, intermediate genealogical divergence index values between the two lineages ultimately leave us unable to recommend a full taxonomic distinction. Genome-wide scans of population genetic statistics identified islands of divergence exhibiting differentiation patterns linked to models of reproductive isolation and within-population selection. Significantly negative values of Tajima's D and positive selection statistics in these islands offer support for selection acting on P. modestum, but patterns may also stem from recent population expansions. We posit that selection within populations has played a large role in shaping genomic divergence across the species' range. Taken together, our results provide perspective into how variable selective pressures shape the genomics of two divergent populations currently maintaining species integrity, despite significant signatures of geographic structure and divergence.