BACKGROUND:Snake venoms are trophic adaptations that represent an ideal model to examine the evolutionary factors that shape polymorphic traits under strong natural selection. Venom compositional variation is substantial within and among venomous snake species. However, the forces shaping this phenotypic complexity, as well as the potential integrated roles of biotic and abiotic factors, have received little attention. Here, we investigate geographic variation in venom composition in a wide-ranging rattlesnake (Crotalus viridis viridis) and contextualize this variation by investigating dietary, phylogenetic, and environmental variables that covary with venom.RESULTS:Using shotgun proteomics, venom biochemical profiling, and lethality assays, we identify 2 distinct divergent phenotypes that characterize major axes of venom variation in this species: a myotoxin-rich phenotype and a snake venom metalloprotease (SVMP)-rich phenotype. We find that dietary availability and temperature-related abiotic factors are correlated with geographic trends in venom composition.CONCLUSIONS:Our findings highlight the potential for snake venoms to vary extensively within species, for this variation to be driven by biotic and abiotic factors, and for the importance of integrating biotic and abiotic variation for understanding complex trait evolution. Links between venom variation and variation in biotic and abiotic factors indicate that venom variation likely results from substantial geographic variation in selection regimes that determine the efficacy of venom phenotypes across populations and snake species. Our results highlight the cascading influence of abiotic factors on biotic factors that ultimately shape venom phenotype, providing evidence for a central role of local selection as a key driver of venom variation.

The eastern Pinesnake (Pituophis melanoleucus) is found throughout eastern United States. Taxonomy in this group has been controversial with several conflicting species designations. Three subspecies of the eastern Pinesnake have prevailed in the literature to their geographic locations and scale coloration: the northern Pinesnake (P. m. melanoleucus), the Florida Pinesnake (P. m. mugitus), and the Black Pinesnake (P. m. lodingi). Within the region, there are several major barriers to dispersal, particularly major river drainage systems and human modification of the longleaf pine habitat. Consistently, a lack of phylogenetic resolution has plagued these taxa in prior studies. The goal of this study was to examine the taxonomic validity of the eastern Pinesnake complex using single nucleotide polymorphisms (SNPs) isolated from ultra-conserved elements (UCEs) in phylogenetic and population genetic approaches. Molecular species delimitation approaches indicated that the population of eastern Pinesnake exhibits population structure across its range that may rise to the level of being new species.

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.

Understanding the processes and mechanisms that promote lineage divergence is a central goal in evolutionary biology. For instance, studies investigating the spatial distribution of genomic variation often highlight biogeographic barriers underpinning geographic isolation or patterns of isolation by environment and isolation by distance that can also lead to lineage divergence. However, the patterns and processes that shape genomic variation and drive lineage divergence may be taxon‐specific, even across closely related taxa co‐occurring within the same biogeographic region. Here, we use molecular data in the form of ultraconserved elements (UCEs) to infer the evolutionary relationships and population genomic structure of the Eastern Pinesnake complex ( Pituophis melanoleucus )—a polytypic wide‐ranging species that occupies much of the Eastern Nearctic. In addition to inferring evolutionary relationships, population genomic structure and gene flow, we also test relationships between genomic diversity and putative barriers to dispersal, environmental variation, and geographic distance. We present results that reveal shallow population genomic structure and ongoing gene flow, despite an extensive geographic range that transcends geographic features found to reduce gene flow among many taxa, including other squamate reptiles within the Eastern Nearctic. Further, our results indicate that the observed genomic diversity is spatially distributed as a pattern of isolation by distance and suggest that the current subspecific taxonomy do not adhere to independent lineages, but, rather, show a significant amount of admixture across the entire P. melanoleucus range.

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.

Understanding the genetic underpinnings of schistosome host preferences is critical. Luo et al. recently identified genes associated with intermediate and definitive host-switching based on a new chromosome-level genome for Schistosoma japonicum, population genetic comparisons, and follow-up experiments. This represents a guide to fully map-selected schistosome genes using population genetics.

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 global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

Schistosomiasis is a neglected tropical disease caused by multiple parasitic Schistosoma species, and which impacts over 200 million people globally, mainly in low- and middle-income countries. Genomic surveillance to detect evidence for natural selection in schistosome populations represents an emerging and promising approach to identify and interpret schistosome responses to ongoing control efforts or other environmental factors. Here we review how genomic variation is used to detect selection, how these approaches have been applied to schistosomes, and how future studies to detect selection may be improved. We discuss the theory of genomic analyses to detect selection, identify experimental designs for such analyses, and review studies that have applied these approaches to schistosomes. We then consider the biological characteristics of schistosomes that are expected to respond to selection, particularly those that may be impacted by control programs. Examples include drug resistance, host specificity, and life history traits, and we review our current understanding of specific genes that underlie them in schistosomes. We also discuss how inherent features of schistosome reproduction and demography pose substantial challenges for effective identification of these traits and their genomic bases. We conclude by discussing how genomic surveillance for selection should be designed to improve understanding of schistosome biology, and how the parasite changes in response to selection.

Schistosomiasis persists in Asian regions despite aggressive elimination measures. To identify factors enabling continued parasite transmission, we performed reduced representation genome sequencing on Schistosoma japonicum miracidia collected across multiple years from transmission hotspots in Sichuan, China. We discovered strong geographic structure, suggesting that local, rather than imported, reservoirs are key sources of persistent infections in the region. At the village level, parasites collected after referral for praziquantel treatment are closely related to local pre-treatment populations. Schistosomes within villages are also highly related, suggesting that only a few parasites from a limited number of hosts drive re-infection. The close familial relationships among miracidia from different human hosts also implicate short transmission routes among humans. At the individual host level, genetic evidence indicates that multiple humans retained infections following referral for treatment. Our findings suggest that end-game schistosomiasis control measures should focus on completely extirpating local parasite reservoirs and confirming successful treatment of infected human hosts.

Male-biased mutation rates occur in a diverse array of organisms. The ratio of male-to-female mutation rate may have major ramifications for evolution across the genome, and for sex-linked genes in particular. In ZW species, the Z chromosome is carried by males two-thirds of the time, leading to the prediction that male-biased mutation rates will have a disproportionate effect on the evolution of Z-linked genes relative to autosomes and the W chromosome. Colubroid snakes (including colubrids, elapids, and viperids) have ZW sex determination, yet male-biased mutation rates have not been well studied in this group. Here we analyze a population genomic dataset from rattlesnakes to quantify genetic variation within and genetic divergence between species. We use a new method for unbiased estimation of population genetic summary statistics to compare variation between the Z chromosome and autosomes and to calculate net nucleotide differentiation between species. We find evidence for a 2.03-fold greater mutation rate in male rattlesnakes relative to females, corresponding to an average μZ/μA ratio of 1.1. Our results from snakes are quantitatively similar to birds, suggesting that male-biased mutation rates may be a common feature across vertebrate lineages with ZW sex determination.

Understanding the processes and mechanisms that promote lineage divergence is a central goal in evolutionary biology. For instance, studies investigating the spatial distribution of genomic variation often highlight biogeographic barriers underpinning geographic isolation, as well as patterns of isolation by environment and isolation by distance that can also lead to lineage divergence. However, the patterns and processes that shape genomic variation and drive lineage divergence may be taxa-specific, even across closely related taxa co-occurring within the same biogeographic region. Here, we use molecular data in the form of ultra-conserved elements (UCEs) to infer the evolutionary relationships and population genomic structure of the Eastern Pinesnake complex ( Pituophis melanoleucus ) – a polytypic wide-ranging species that occupies much of the Eastern Nearctic. In addition to inferring evolutionary relationships, population genomic structure, and gene flow, we also test relationships between genomic diversity and putative barriers to dispersal, environmental variation, and geographic distance. We present results that reveal shallow population genomic structure and ongoing gene flow, despite an extensive geographic range that transcends geographic features found to reduce gene flow among many taxa, including other squamate reptiles within the Eastern Nearctic. Further, our results indicate that the observed genomic diversity is spatially distributed as a pattern of isolation by distance and suggest that the current subspecific taxonomy do not adhere to independent lineages, but rather, show a significant amount of admixture across the entire P. melanoleucus range.

Genomic approaches hold great promise for resolving unanswered questions about transmission patterns and responses to control efforts for schistosomiasis and other neglected tropical diseases. However, the cost of generating genomic data and the challenges associated with obtaining sufficient DNA from individual schistosome larvae (miracidia) from mammalian hosts have limited the application of genomic data for studying schistosomes and other complex macroparasites. Here, we demonstrate the feasibility of utilizing whole genome amplification and sequencing (WGS) to analyze individual archival miracidia. As an example, we sequenced whole genomes of 22 miracidia from 11 human hosts representing two villages in rural Sichuan, China, and used these data to evaluate patterns of relatedness and genetic diversity. We also down-sampled our dataset to test how lower coverage sequencing could increase the cost effectiveness of WGS while maintaining power to accurately infer relatedness. Collectively, our results illustrate that population-level WGS datasets are attainable for individual miracidia and represent a powerful tool for ultimately providing insight into overall genetic diversity, parasite relatedness, and transmission patterns for better design and evaluation of disease control efforts.

Over the past decade, the increasing availability and decreasing cost of genome sequencing technologies have led to significant advancements in our understanding of the evolution, variation and function of reptile venoms. This includes multiple annotated genomes for venomous reptile species as well as countless public data sets of venom gene sequences, gene expression data, and data availability from a growing diversity of species. These genomic resources have provided new insight into reptile venoms, and the continuing development of additional genomic resources and related data sets holds promise for continued advancement of our understanding of reptile venom systems. Here, we present an overview of currently available genomic resources for the study of reptile venoms, discuss genomic approaches used to study reptile venoms, and highlight recent advances in our understanding of how venom evolves and functions in reptiles based on genomic data. We also discuss ongoing progress and future work that will be important for closing major gaps in our current understanding of the genomic basis of venom variation, evolution and regulation.

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.

The study of recently diverged lineages whose geographical ranges come into contact can provide insight into the early stages of speciation and the potential roles of reproductive isolation in generating and maintaining species. Such insight can also be important for understanding the strategies and challenges for delimiting species within recently diverged species complexes. Here, we use mitochondrial and nuclear genetic data to study population structure, gene flow and demographic history across a geographically widespread rattlesnake clade, the western rattlesnake species complex (Crotalus cerberus, Crotalus viridis, Crotalus oreganus and relatives), which contains multiple lineages with ranges that overlap geographically or contact one another. We find evidence that the evolutionary history of this group does not conform to a bifurcating tree model and that pervasive gene flow has broadly influenced patterns of present-day genetic diversity. Our results suggest that lineage diversity has been shaped largely by drift and divergent selection in isolation, followed by secondary contact, in which reproductive isolating mechanisms appear weak and insufficient to prevent introgression, even between anciently diverged lineages. The complexity of divergence and secondary contact with gene flow among lineages also provides new context for why delimiting species within this complex has been difficult and contentious historically.