Assembly of the largest squamate reference genome to date: The western fence lizard, Sceloporus occidentalis.

Spiny lizards (genus Sceloporus) have long served as important systems for studies of behavior, thermal physiology, dietary ecology, vector biology, speciation, and biogeography. The western fence lizard, Sceloporus occidentalis, is found across most of the major biogeographical regions in the western United States and northern Baja California, Mexico, inhabiting a wide range of habitats, from grassland to chaparral to open woodlands. As small ectotherms, Sceloporus lizards are particularly vulnerable to climate change, and S. occidentalis has also become an important system for studying the impacts of land use change and urbanization on small vertebrates. Here, we report a new reference genome assembly for S. occidentalis, as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomics strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 608 scaffolds spanning 2,856 Mb, has a contig N50 of 18.9 Mb, a scaffold N50 of 98.4 Mb, and BUSCO completeness score of 98.1% based on the tetrapod gene set. This reference genome will be valuable for understanding ecological and evolutionary dynamics in S. occidentalis, the species status of the California endemic island fence lizard (S. becki), and the spectacular radiation of Sceloporus lizards.

Temporal climatic variability predicts thermal tolerance in two sympatric lizard species.

Thermal acclimatization, plastic shifts in thermal physiology in response to recent climatic conditions, is thought to be adaptive in highly seasonal environments where thermal variability is high but predictable. Thus, lizards from mid-latitude, desert environments should exhibit plasticity in their thermal tolerance limits, the upper (CTmax) and lower (CTmin) body temperatures they can withstand while maintaining physiological functioning, associated with changes in seasonal changes in climatic variation (i.e., when daily fluctuations in temperature are greater, lizards should have wider thermal tolerance breadths [CTmax-CTmin]). We measured the thermal tolerance limits of two Phrynosomatid lizard species, Uta stansburiana and Sceloporus tristichus, occurring in sympatry at three time points to test for temporal variation in thermal physiology in response to climatic variation. We found that lizards of both species measured during times when climatic variability was high had wider thermal tolerance breadths than lizards measured when climatic variability was lower. While CTmax was largely invariable, CTmin varied in response to minimum air temperature, driving the observed difference in thermal tolerance breadth among the sampling periods.

Phylogenomics of a rapid radiation: is chromosomal evolution linked to increased diversification in north american spiny lizards (Genus Sceloporus)?

BACKGROUND: Resolving the short phylogenetic branches that result from rapid evolutionary diversification often requires large numbers of loci. We collected targeted sequence capture data from 585 nuclear loci (541 ultraconserved elements and 44 protein-coding genes) to estimate the phylogenetic relationships among iguanian lizards in the North American genus Sceloporus. We tested for diversification rate shifts to determine if rapid radiation in the genus is correlated with chromosomal evolution. RESULTS: The phylogenomic trees that we obtained for Sceloporus using concatenation and coalescent-based species tree inference provide strong support for the monophyly and interrelationships among nearly all major groups. The diversification analysis supported one rate shift on the Sceloporus phylogeny approximately 20-25 million years ago that is associated with the doubling of the speciation rate from 0.06 species/million years (Ma) to 0.15 species/Ma. The posterior probability for this rate shift occurring on the branch leading to the Sceloporus species groups exhibiting increased chromosomal diversity is high (posterior probability = 0.997). CONCLUSIONS: Despite high levels of gene tree discordance, we were able to estimate a phylogenomic tree for Sceloporus that solves some of the taxonomic problems caused by previous analyses of fewer loci. The taxonomic changes that we propose using this new phylogenomic tree help clarify the number and composition of the major species groups in the genus. Our study provides new evidence for a putative link between chromosomal evolution and the rapid divergence and radiation of Sceloporus across North America.

Short Tree, Long Tree, Right Tree, Wrong Tree: New Acquisition Bias Corrections for Inferring SNP Phylogenies.

Single nucleotide polymorphisms (SNPs) are useful markers for phylogenetic studies owing in part to their ubiquity throughout the genome and ease of collection. Restriction site associated DNA sequencing (RADseq) methods are becoming increasingly popular for SNP data collection, but an assessment of the best practises for using these data in phylogenetics is lacking. We use computer simulations, and new double digest RADseq (ddRADseq) data for the lizard family Phrynosomatidae, to investigate the accuracy of RAD loci for phylogenetic inference. We compare the two primary ways RAD loci are used during phylogenetic analysis, including the analysis of full sequences (i.e., SNPs together with invariant sites), or the analysis of SNPs on their own after excluding invariant sites. We find that using full sequences rather than just SNPs is preferable from the perspectives of branch length and topological accuracy, but not of computational time. We introduce two new acquisition bias corrections for dealing with alignments composed exclusively of SNPs, a conditional likelihood method and a reconstituted DNA approach. The conditional likelihood method conditions on the presence of variable characters only (the number of invariant sites that are unsampled but known to exist is not considered), while the reconstituted DNA approach requires the user to specify the exact number of unsampled invariant sites prior to the analysis. Under simulation, branch length biases increase with the amount of missing data for both acquisition bias correction methods, but branch length accuracy is much improved in the reconstituted DNA approach compared to the conditional likelihood approach. Phylogenetic analyses of the empirical data using concatenation or a coalescent-based species tree approach provide strong support for many of the accepted relationships among phrynosomatid lizards, suggesting that RAD loci contain useful phylogenetic signal across a range of divergence times despite the presence of missing data. Phylogenetic analysis of RAD loci requires careful attention to model assumptions, especially if downstream analyses depend on branch lengths.

Repeated evolution of viviparity in phrynosomatid lizards constrained interspecific diversification in some life-history traits.

In vertebrates, viviparity has evolved independently multiple times, apparently increasing morphological diversification and speciation rates as a consequence. We tested whether the evolution of viviparity has also increased diversification of life-history traits by estimating evolutionary rates of lizards from the North American family Phrynosomatidae. Using modern phylogenetic comparative methods, we compared these rates between oviparous and viviparous species, and found no support for this hypothesis. Instead, we found higher evolutionary rates for oviparous species in some life-history traits. Our results suggest that the evolution of viviparity may have constrained rather than facilitated evolution of life histories.

The effect of parity on morphological evolution among phrynosomatid lizards.

The shift from egg laying to live-bearing is one of the most well-studied transitions in evolutionary biology. Few studies, however, have assessed the effect of this transition on morphological evolution. Here, we evaluated the effect of reproductive mode on the morphological evolution of 10 traits, among 108 species of phrynosomatid lizards. We assess whether the requirement for passing shelled eggs through the pelvic girdle has led to morphological constraints in oviparous species and whether long gestation times in viviparous species have led to constraints in locomotor morphology. We fit models to the data that vary both in their tempo (strength and rate of selection) and mode of evolution (Brownian or Ornstein-Uhlenbeck) and estimates of trait optima. We found that most traits are best fit by a generalized multipeak OU model, suggesting differing trait optima for viviparous vs. oviparous species. Additionally, rates (σ(2) ) of both pelvic girdle and forelimb trait evolution varied with parity; viviparous species had higher rates. Hindlimb traits, however, exhibited no difference in σ(2) between parity modes. In a functional context, our results suggest that the passage of shelled eggs constrains the morphology of the pelvic girdle, but we found no evidence of morphological constraint of the locomotor apparatus in viviparous species. Our results are consistent with recent lineage diversification analyses, leading to the conclusion that transitions to viviparity increase both lineage and morphological diversification.

Information out of the blue: phenotypic correlates of abdominal color patches in Sceloporus lizards.

Colorful ornaments are important visual signals for animal communication that can provide critical information about the quality of the signaler. In this study, we focused on different color characteristics of the abdominal patches of males of six lizard species from the genus Sceloporus. We addressed three main objectives. First, we examined if size, brightness, saturation, and conspicuousness of these ornaments are indicative of body size, condition, immune function, or levels of testosterone and corticosterone. Second, we evaluated if the distinct components of these abdominal patches (blue or green patches and black stripes) transmit similar information about the signaler, which would support the redundant signal hypothesis, or if these components are related to different phenotypic traits, which would support the multiple message hypothesis. Third, we compared the phenotypic correlates of these ornaments among our six species to understand the degree of conservatism in the signaling patterns or to find species-specific signals. Using data collected from males in natural conditions and a multi-model inference framework, we found that in most species the area of the patches and the brightness of the blue component are positively related to body size. Thus, these color characteristics are presumably indicative of the physical strength and competitive ability of males and these shared signals were likely inherited from a common ancestor. In half of the species, males in good body condition also exhibit relatively larger blue and black areas, suggesting that the expression of these ornaments is condition-dependent. Abdominal patches also provide information about immunocompetence of the males as indicated by different correlations between certain color characteristics and ectoparasite load, counts of heterophils, and the heterophil:lymphocyte ratio. Our findings reveal that area and brightness of the abdominal patches signal the size and body condition of males, whereas blue saturation and conspicuousness with respect to the surrounding substrate are indicative of immune condition, thus supporting the multiple message hypothesis. However, some of these correlations were not shared by all species and, hence, point to intriguing species-specific signals.

The complete mitochondrial genome of the horned lizard Phrynosoma blainvillii (Squamata: Phrynosomatidae) from California, USA.

Analysis of Phrynosoma blainvillii Gray from Marina, Monterey County, California, using 150 bp paired-end Illumina sequences (Illumina, San Diego, CA) resulted in the assembly of its complete mitogenome. The mitogenome is 16,946 bp in length and contains a putative origin of light strand replication (OL), control region, 22 tRNA, 2 rRNA, and 13 protein-coding genes. Its content and organization are similar to other Squamata. Phylogenetic analysis of P. blainvillii resolves it in a clade with P. sherbrookei Nieto-Montes de Oca, Arenas-Moreno, Beltrán-Sánchez & Leaché, sister in position to Uma notata Baird. Mitochondrial marker analysis of P. blainvillii from Marina shows that it belongs to a coastal Santa Lucia Mountain Range haplogroup that is distinct from other populations of P. blainvillii in California.

Two new species of Parapharyngodon parasites of Sceloporus pyrocephalus, with a key to the species found in Mexico (Nematoda, Pharyngodonidae).

Two new species of Parapharyngodon collected from the intestine of the Mexican boulder spiny lizard Sceloporus pyrocephalus are described. This study increases to 49 the number of valid species assigned to Parapharyngodon worldwide, 11 of them distributed in Mexico. Males of the two new species share the presence of four pairs of caudal papillae, an anterior echinate cloacal lip and the presence of lateral alae; however, both differ from each other in lateral alae extension and echinate cloacal anterior lip morphology. Females of both species have a prebulbar uterus and eggs shell punctuate with pores, characteristics shared with few other species of Parapharyngodon. Both new species differ from other congeneric species in the papillar arrangement, the anterior cloacal lip morphology, the lateral alae extension and total length/spicule ratio. A taxonomic key for the species of Parapharyngodon distributed in Mexico is provided.

The ultrastructure of spermatid development during spermiogenesis within the rosebelly lizard, Sceloporus variabilis (Reptilia, Squamata, Phrynosomatidae).

Several recent studies have mapped out the characters of spermiogenesis within several species of squamates. Many of these data have shown both conserved and possibly apomorphic morphological traits that could be important in future phylogenetic analysis within Reptilia. There, however, has not been a recent study that compares spermiogenesis and its similarities or differences between two species of reptile that reside in the same genus. Thus, the present analysis details the changes to spermiogenesis in Sceloporus variabilis and then compares spermatid morphologies to that of Sceloporus bicanthalis. Many of the morphological changes that the spermatids undergo in these two species are similar or conserved, which is similar to what has been reported in other squamates. There are six main character differences that can be observed during the development of the spermatids between these two sceloporid lizards. They include the presence (S. variabilis) or absence (S. bicanthalis) of a mitochondrial/endoplasmic reticulum complex near the Golgi apparatus during acrosome development, a shallow (S. variabilis) or deep (S. bicanthalis) nuclear indentation that accommodates the acrosomal vesicle, filamentous (S. variabilis) or granular (S. bicanthalis) chromatin condensation, no spiraling (S. variabilis) or spiraling (S. bicanthalis) of chromatin during condensation, absence (S. variabilis) or presence (S. bicanthalis) of the longitudinal manchette microtubules, and the lack of (S. variabilis) or presence (S. bicanthalis) of nuclear lacunae. This is the first study that compares spermiogenic ultrastructural characters between species within the same genus. The significance of the six character differences between two distantly related species within Sceloporus is still unknown, but these data do suggest that spermiogenesis might be a good model to study the hypothesis that spermatid ontogeny is species specific.