Regionalization of the vertebral column and its correlation with heart position in snakes: Implications for evolutionary pathways and morphological diversification.

Spinal regionalization has important implications for the evolution of vertebrate body plans. We determined the variation in the number and morphology of vertebrae across the vertebral column (i.e., vertebral formula) for 63 snake species representing 13 families using intracolumnar variation in vertebral shape. Vertebral counts were used to determine the position of the heart, pylorus, and left kidney for each species. Across all species we observed a conspicuous midthoracic transition in vertebral shape, indicating four developmental domains of the precloacal vertebral column (cervical, anterior thoracic, posterior thoracic, and lumbar). Using phylogenetic analyses, the boundary between the anterior and posterior thoracic vertebrae was correlated with heart position. No associations were found between shifts in morphology of the vertebral column and either the pylorus or left kidney. We observed that among taxa, the number of preapex and postapex vertebrae could change independently from one another and from changes in the total number of precloacal vertebrae. Ancestral state reconstruction of the preapex and postapex vertebrae illustrated several evolutionary pathways by which diversity in the vertebral column and heart position have been attained. In addition, no conspicuous pattern was observed among the heart, pylorus, or kidney indicating that their relative positions to each other evolve independently. We conclude that snakes exhibit four morphologically distinct regions of the vertebral column. We discuss the implications of the forebody and hindbody vertebral formula on the morphological diversification of snakes.

The best of both worlds? Rattlesnake hybrid zones generate complex combinations of divergent venom phenotypes that retain high toxicity.

Studying the consequences of hybridization between closely related species with divergent traits can reveal patterns of evolution that shape and maintain extreme trophic adaptations. Snake venoms are an excellent model system for examining the evolutionary and ecological patterns that underlie highly selected polymorphic traits. Here we investigate hybrid venom phenotypes that result from natural introgression between two rattlesnake species that express highly divergent venom phenotypes: Crotalus o. concolor and C. v. viridis. Though not yet documented, interbreeding between these species may lead to novel venom phenotypes with unique activities that break the typical trends of venom composition in rattlesnakes. The characteristics of these unusual phenotypes could unveil the roles of introgression in maintaining patterns of venom composition and variation, including the near ubiquitous dichotomy between neurotoxic or degradative venoms observed across rattlesnakes. We use RADseq data to infer patterns of gene flow and hybrid ancestry between these diverged lineages and link these genetic data with analyses of venom composition, biological activity, and whole animal model toxicity tests to understand the impacts of introgression on venom composition. We find that introgressed populations express admixed venom phenotypes that do not sacrifice biological activity (lethal toxicity) or overall abundance of dominant toxins compared to parental venoms. These hybridized venoms therefore do not represent a trade-off in functionality between the typical phenotypic extremes but instead represent a unique combination of characters whose expression appears limited to the hybrid zone.

When Food Fights Back: Skull Morphology and Feeding Behavior of Centipede-Eating Snakes.

Feeding is a complex process that involves an integrated response of multiple functional systems. Animals evolve phenotypic integration of complex morphological traits to covary and maximize performance of feeding behaviors. Specialization, such as feeding on dangerous prey, can further shape the integration of behavior and morphology as traits are expected to evolve and maintain function in parallel. Feeding on centipedes, with their powerful forcipules that pinch and inject venom, has evolved multiple times within snakes, including the genus Tantilla. However, the behavioral and morphological adaptations used to consume this dangerous prey are poorly understood. By studying snakes with varying degrees of dietary specialization, we can test the integration of diet, morphology, and behavior to better understand the evolution of consuming difficult prey. We studied the prey preference and feeding behavior of Tantilla using the flat-headed snake (T. gracilis) and the crowned snake (T. coronata), which differ in the percentage of centipedes in their diet. We then quantified cranial anatomy using geometric morphometric data from CT scans. To test prey preference, we offered multiple types of prey and recorded snake behavior. Both species of snakes showed interest in multiple prey types, but only struck or consumed centipedes. To subdue centipedes, crowned snakes used coiling and holding (envenomation) immediately after striking, while flat-headed snakes used the novel behavior of pausing and holding onto centipedes for a prolonged time prior to the completion of swallowing. Each skull element differed in shape after removing the effects of size, position, and orientation. The rear fang was larger in crowned snakes, but the mechanical advantage of the lower jaw was greater in flat-headed snakes. Our results suggest that the integration of behavioral and morphological adaptations is important for the success of subduing and consuming dangerous prey.

Evidence that genomic incompatibilities and other multilocus processes impact hybrid fitness in a rattlesnake hybrid zone.

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

The roles of balancing selection and recombination in the evolution of rattlesnake venom.

The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance.

Use of scalation landmarks in geometric morphometrics of squamate reptiles: a comment on homology.

Geometric morphometrics (GM) is a powerful analytical approach for evaluating phenotypic variation relevant to taxonomy and systematics, and as with any statistical methodology, requires adherence to fundamental assumptions for inferences to be strictly valid. An important consideration for GM is how landmark configurations, which represent sets of anatomical loci for evaluating shape variation through Cartesian coordinates, relate to underlying homology (Zelditch et al. 1995; Polly 2008). Perhaps more so than with traditional morphometrics, anatomical homology is a crucial assumption for GM because of the mathematical and biological interpretations associated with shape change depicted by deformation grids, such as the thin plate spline (Klingenberg 2008; Zelditch et al. 2012). GM approaches are often used to analyze shapes or outlines of structures, which are not necessarily related to common ancestry, and in this respect GM approaches that use linear semi-landmarks and related methods are particularly amenable to evaluating primary homology, or raw similarity between structures (De Pinna 1991; Palci Lee 2019). This relaxed interpretation of homology that focuses more on recognizable and repeatable landmarks is defensible so long as authors are clear regarding the purpose of the analyses and in defining their landmark configurations (Palci Lee 2019). Secondary homology, or similarity due to common ancestry, can also be represented with GM methods and is often assumed to be reflected in fixed Type 1 (juxtaposition of tissues) or Type 2 (self-evident geometry) landmarks (Bookstein 1991).

Snake Recombination Landscapes Are Concentrated in Functional Regions despite PRDM9.

Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions-a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions-a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages.

Geographic variation in morphology in the Mohave Rattlesnake (Crotalus scutulatus Kennicott 1861) (Serpentes: Viperidae): implications for species boundaries.

The Mohave Rattlesnake (Crotalus scutulatus) is a highly venomous pitviper inhabiting the arid interior deserts, grasslands, and savannas of western North America. Currently two subspecies are recognized: the Northern Mohave Rattlesnake (C. s. scutulatus) ranging from southern California to the southern Central Mexican Plateau, and the Huamantla Rattlesnake (C. s. salvini) from the region of Tlaxcala, Veracruz, and Puebla in south-central Mexico. Although recent studies have demonstrated extensive geographic variation in venom composition and cryptic genetic diversity in this species, no modern studies have focused on geographic variation in morphology. Here we analyzed a series of qualitative, meristic, and morphometric traits from 347 specimens of C. scutulatus and show that this species is phenotypically cohesive without discrete subgroups, and that morphology follows a continuous cline in primarily color pattern and meristic traits across the major axis of its expansive distribution. Interpreted in the context of previously published molecular evidence, our morphological analyses suggest that multiple episodes of isolation and secondary contact among metapopulations during the Pleistocene were sufficient to produce distinctive genetic populations, which have since experienced gene flow to produce clinal variation in phenotypes without discrete or diagnosable distinctions among these original populations. For taxonomic purposes, we recommend that C. scutulatus be retained as a single species, although it is possible that C. s. salvini, which is morphologically the most distinctive population, could represent a peripheral isolate in the initial stages of speciation.

The origins and evolution of chromosomes, dosage compensation, and mechanisms underlying venom regulation in snakes.

Here we use a chromosome-level genome assembly of a prairie rattlesnake (Crotalus viridis), together with Hi-C, RNA-seq, and whole-genome resequencing data, to study key features of genome biology and evolution in reptiles. We identify the rattlesnake Z Chromosome, including the recombining pseudoautosomal region, and find evidence for partial dosage compensation driven by an evolutionary accumulation of a female-biased up-regulation mechanism. Comparative analyses with other amniotes provide new insight into the origins, structure, and function of reptile microchromosomes, which we demonstrate have markedly different structure and function compared to macrochromosomes. Snake microchromosomes are also enriched for venom genes, which we show have evolved through multiple tandem duplication events in multiple gene families. By overlaying chromatin structure information and gene expression data, we find evidence for venom gene-specific chromatin contact domains and identify how chromatin structure guides precise expression of multiple venom gene families. Further, we find evidence for venom gland-specific transcription factor activity and characterize a complement of mechanisms underlying venom production and regulation. Our findings reveal novel and fundamental features of reptile genome biology, provide insight into the regulation of snake venom, and broadly highlight the biological insight enabled by chromosome-level genome assemblies.

Cryptic genetic diversity, population structure, and gene flow in the Mojave rattlesnake (Crotalus scutulatus).

The Mojave rattlesnake (Crotalus scutulatus) inhabits deserts and arid grasslands of the western United States and Mexico. Despite considerable interest in its highly toxic venom and the recognition of two subspecies, no molecular studies have characterized range-wide genetic diversity and population structure or tested species limits within C. scutulatus. We used mitochondrial DNA and thousands of nuclear loci from double-digest restriction site associated DNA sequencing to infer population genetic structure throughout the range of C. scutulatus, and to evaluate divergence times and gene flow between populations. We find strong support for several divergent mitochondrial and nuclear clades of C. scutulatus, including splits coincident with two major phylogeographic barriers: the Continental Divide and the elevational increase associated with the Central Mexican Plateau. We apply Bayesian clustering, phylogenetic inference, and coalescent-based species delimitation to our nuclear genetic data to test hypotheses of population structure. We also performed demographic analyses to test hypotheses relating to population divergence and gene flow. Collectively, our results support the existence of four distinct lineages within C. scutulatus, and genetically defined populations do not correspond with currently recognized subspecies ranges. Finally, we use approximate Bayesian computation to test hypotheses of divergence among multiple rattlesnake species groups distributed across the Continental Divide, and find evidence for co-divergence at this boundary during the mid-Pleistocene.

Minimum area thresholds for rattlesnakes and colubrid snakes on islands in the Gulf of California, Mexico.

We expand a framework for estimating minimum area thresholds to elaborate biogeographic patterns between two groups of snakes (rattlesnakes and colubrid snakes) on islands in the western Gulf of California, Mexico. The minimum area thresholds for supporting single species versus coexistence of two or more species relate to hypotheses of the relative importance of energetic efficiency and competitive interactions within groups, respectively. We used ordinal logistic regression probability functions to estimate minimum area thresholds after evaluating the influence of island area, isolation, and age on rattlesnake and colubrid occupancy patterns across 83 islands. Minimum area thresholds for islands supporting one species were nearly identical for rattlesnakes and colubrids (~1.7 km2), suggesting that selective tradeoffs for distinctive life history traits between rattlesnakes and colubrids did not result in any clear advantage of one life history strategy over the other on islands. However, the minimum area threshold for supporting two or more species of rattlesnakes (37.1 km2) was over five times greater than it was for supporting two or more species of colubrids (6.7 km2). The great differences between rattlesnakes and colubrids in minimum area required to support more than one species imply that for islands in the Gulf of California relative extinction risks are higher for coexistence of multiple species of rattlesnakes and that competition within and between species of rattlesnakes is likely much more intense than it is within and between species of colubrids.

Allometry of sexual size dimorphism in turtles: a comparison of mass and length data.

BACKGROUND: The macroevolutionary pattern of Rensch's Rule (positive allometry of sexual size dimorphism) has had mixed support in turtles. Using the largest carapace length dataset and only large-scale body mass dataset assembled for this group, we determine (a) whether turtles conform to Rensch's Rule at the order, suborder, and family levels, and (b) whether inferences regarding allometry of sexual size dimorphism differ based on choice of body size metric used for analyses. METHODS: We compiled databases of mean body mass and carapace length for males and females for as many populations and species of turtles as possible. We then determined scaling relationships between males and females for average body mass and straight carapace length using traditional and phylogenetic comparative methods. We also used regression analyses to evalutate sex-specific differences in the variance explained by carapace length on body mass. RESULTS: Using traditional (non-phylogenetic) analyses, body mass supports Rensch's Rule, whereas straight carapace length supports isometry. Using phylogenetic independent contrasts, both body mass and straight carapace length support Rensch's Rule with strong congruence between metrics. At the family level, support for Rensch's Rule is more frequent when mass is used and in phylogenetic comparative analyses. Turtles do not differ in slopes of sex-specific mass-to-length regressions and more variance in body size within each sex is explained by mass than by carapace length. DISCUSSION: Turtles display Rensch's Rule overall and within families of Cryptodires, but not within Pleurodire families. Mass and length are strongly congruent with respect to Rensch's Rule across turtles, and discrepancies are observed mostly at the family level (the level where Rensch's Rule is most often evaluated). At macroevolutionary scales, the purported advantages of length measurements over weight are not supported in turtles.

A reconnaissance analysis of groundwater quality in the Eagle Ford shale region reveals two distinct bromide/chloride populations.

The extraction of oil and natural gas from unconventional shale formations has prompted a series of investigations to examine the quality of the groundwater in the overlying aquifers. Here we present a reconnaissance analysis of groundwater quality in the Eagle Ford region of southern Texas. These data reveal two distinct sample populations that are differentiable by bromide/chloride ratios. Elevated levels of fluoride, nitrate, sulfate, various metal ions, and the detection of exotic volatile organic compounds highlight a high bromide group of samples, which is geographically clustered, while encompassing multiple hydrogeological strata. Samples with bromide/chloride ratios representative of connate water displayed elevated levels of total organic carbon, while revealing the detection of alcohols and chlorinated compounds. These findings suggest that groundwater quality in the Western Gulf Basin is, for the most part, controlled by a series of natural processes; however, there is also evidence of episodic contamination events potentially attributed to unconventional oil and gas development or other anthropogenic activities. Collectively, this characterization of natural groundwater constituents and exogenous compounds will guide targeted remediation efforts and provides insight for agricultural entities, industrial operators, and rural communities that rely on groundwater in southern Texas.

Point source attribution of ambient contamination events near unconventional oil and gas development.

We present an analysis of ambient benzene, toluene, and xylene isomers in the Eagle Ford shale region of southern Texas. In situ air quality measurements using membrane inlet mobile mass spectrometry revealed ambient benzene and toluene concentrations as high as 1000 and 5000 parts-per-billion, respectively, originating from specific sub-processes on unconventional oil and gas well pad sites. The detection of highly variant contamination events attributable to natural gas flaring units, condensate tanks, compressor units, and hydrogen sulfide scavengers indicates that mechanical inefficiencies, and not necessarily the inherent nature of the extraction process as a whole, result in the release of these compounds into the environment. This awareness of ongoing contamination events contributes to an enhanced knowledge of ambient volatile organic compounds on a regional scale. While these reconnaissance measurements on their own do not fully characterize the fluctuations of ambient BTEX concentrations that likely exist in the atmosphere of the Eagle Ford Shale region, they do suggest that contamination events from unconventional oil and gas development can be monitored, controlled, and reduced.

Temporal variation in groundwater quality in the Permian Basin of Texas, a region of increasing unconventional oil and gas development.

The recent expansion of natural gas and oil extraction using unconventional oil and gas development (UD) practices such as horizontal drilling and hydraulic fracturing has raised questions about the potential for environmental impacts. Prior research has focused on evaluations of air and water quality in particular regions without explicitly considering temporal variation; thus, little is known about the potential effects of UD activity on the environment over longer periods of time. Here, we present an assessment of private well water quality in an area of increasing UD activity over a period of 13months. We analyzed samples from 42 private water wells located in three contiguous counties on the Eastern Shelf of the Permian Basin in Texas. This area has experienced a rise in UD activity in the last few years, and we analyzed samples in four separate time points to assess variation in groundwater quality over time as UD activities increased. We monitored general water quality parameters as well as several compounds used in UD activities. We found that some constituents remained stable over time, but others experienced significant variation over the period of study. Notable findings include significant changes in total organic carbon and pH along with ephemeral detections of ethanol, bromide, and dichloromethane after the initial sampling phase. These data provide insight into the potentially transient nature of compounds associated with groundwater contamination in areas experiencing UD activity.

Limitations of climatic data for inferring species boundaries: insights from speckled rattlesnakes.

Phenotypes, DNA, and measures of ecological differences are widely used in species delimitation. Although rarely defined in such studies, ecological divergence is almost always approximated using multivariate climatic data associated with sets of specimens (i.e., the "climatic niche"); the justification for this approach is that species-specific climatic envelopes act as surrogates for physiological tolerances. Using identical statistical procedures, we evaluated the usefulness and validity of the climate-as-proxy assumption by comparing performance of genetic (nDNA SNPs and mitochondrial DNA), phenotypic, and climatic data for objective species delimitation in the speckled rattlesnake (Crotalus mitchellii) complex. Ordination and clustering patterns were largely congruent among intrinsic (heritable) traits (nDNA, mtDNA, phenotype), and discordance is explained by biological processes (e.g., ontogeny, hybridization). In contrast, climatic data did not produce biologically meaningful clusters that were congruent with any intrinsic dataset, but rather corresponded to regional differences in atmospheric circulation and climate, indicating an absence of inherent taxonomic signal in these data. Surrogating climate for physiological tolerances adds artificial weight to evidence of species boundaries, as these data are irrelevant for that purpose. Based on the evidence from congruent clustering of intrinsic datasets, we recommend that three subspecies of C. mitchellii be recognized as species: C. angelensis, C. mitchellii, and C. Pyrrhus.

A Comprehensive Analysis of Groundwater Quality in The Barnett Shale Region.

The exploration of unconventional shale energy reserves and the extensive use of hydraulic fracturing during well stimulation have raised concerns about the potential effects of unconventional oil and gas extraction (UOG) on the environment. Most accounts of groundwater contamination have focused primarily on the compositional analysis of dissolved gases to address whether UOG activities have had deleterious effects on overlying aquifers. Here, we present an analysis of 550 groundwater samples collected from private and public supply water wells drawing from aquifers overlying the Barnett shale formation of Texas. We detected multiple volatile organic carbon compounds throughout the region, including various alcohols, the BTEX family of compounds, and several chlorinated compounds. These data do not necessarily identify UOG activities as the source of contamination; however, they do provide a strong impetus for further monitoring and analysis of groundwater quality in this region as many of the compounds we detected are known to be associated with UOG techniques.

The phenotypic effects of spontaneous mutations in different environments.

Understanding the context dependence of mutation represents the current frontier of mutation research. In particular, understanding how traits vary in their abilities to accrue mutational variation and how the environment influences expression of mutant phenotypes yields insight into evolutionary processes. We conducted phenotypic assays in four environments using a set of Daphnia pulex mutation accumulation lines to examine the context dependence of mutation. Life-history traits accrued mutational variance faster than morphological traits when considered in individual environments. Across environments, the mutational variance in plasticity was also greater for life-history traits than for morphological traits, although this pattern was less robust. In addition, the expression of mutational variance depended on the environment, which resulted in changes in the rank order of genotype performance across environments in some cases. Such cryptic genetic variation resulting from mutation may maintain genetic diversity and allow for rapid adaptation in spatially or temporally variable environments.