The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom.

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake.

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

Topology Testing and Demographic Modeling Illuminate a Novel Speciation Pathway in the Greater Caribbean Sea Following the Formation of the Isthmus of Panama.

Recent genomic analyses have highlighted the prevalence of speciation with gene flow in many taxa and have underscored the importance of accounting for these reticulate evolutionary processes when constructing species trees and generating parameter estimates. This is especially important for deepening our understanding of speciation in the sea where fast-moving ocean currents, expanses of deep water, and periodic episodes of sea level rise and fall act as soft and temporary allopatric barriers that facilitate both divergence and secondary contact. Under these conditions, gene flow is not expected to cease completely while contemporary distributions are expected to differ from historical ones. Here, we conduct range-wide sampling for Pederson's cleaner shrimp (Ancylomenes pedersoni), a species complex from the Greater Caribbean that contains three clearly delimited mitochondrial lineages with both allopatric and sympatric distributions. Using mtDNA barcodes and a genomic ddRADseq approach, we combine classic phylogenetic analyses with extensive topology testing and demographic modeling (10 site frequency replicates × 45 evolutionary models × 50 model simulations/replicate = 22,500 simulations) to test species boundaries and reconstruct the evolutionary history of what was expected to be a simple case study. Instead, our results indicate a history of allopatric divergence, secondary contact, introgression, and endemic hybrid speciation that we hypothesize was driven by the final closure of the Isthmus of Panama and the strengthening of the Gulf Stream Current ~3.5 Ma. The history of this species complex recovered by model-based methods that allow reticulation differs from that recovered by standard phylogenetic analyses and is unexpected given contemporary distributions. The geologically and biologically meaningful insights gained by our model selection analyses illuminate what is likely a novel pathway of species formation not previously documented that resulted from one of the most biogeographically significant events in Earth's history.

Genetic mechanisms and biological processes underlying host response to ophidiomycosis (snake fungal disease) inferred from tissue-specific transcriptome analyses.

Emerging infectious diseases in wildlife species caused by pathogenic fungi are of growing concern, yet crucial knowledge gaps remain for diseases with potentially large impacts. For example, there is detailed knowledge about host pathology and mechanisms underlying response for chytridiomycosis in amphibians and white-nose syndrome in bats, but such information is lacking for other more recently described fungal infections. One such disease is ophidiomycosis, caused by the fungus Ophidiomyces ophidiicola, which has been identified in many species of snakes, yet the biological mechanisms and molecular changes occurring during infection are unknown. To gain this information, we performed a controlled experimental infection in captive Prairie rattlesnakes (Crotalus viridis) with O. ophidiicola at two different temperatures: 20 and 26°C. We then compared liver, kidney, and skin transcriptomes to assess tissue-specific genetic responses to O. ophidiicola infection. Given previous histopathological studies and the fact that snakes are ectotherms, we expected highest fungal activity on skin and a significant impact of temperature on host response. Although we found fungal activity to be localized on skin, most of the differential gene expression occurred in internal tissues. Infected snakes at the lower temperature had the highest host mortality whereas two-thirds of the infected snakes at the higher temperature survived. Our results suggest that ophidiomycosis is likely a systemic disease with long-term effects on host response. Our analysis also identified candidate protein coding genes that are potentially involved in host response, providing genetic tools for studies of host response to ophidiomycosis in natural populations.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome.

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers.

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

Venom Gene Sequence Diversity and Expression Jointly Shape Diet Adaptation in Pitvipers.

Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA2s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms-such as higher absolute levels of expression-are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey.

Genetic characterization of potential venom resistance proteins in California ground squirrels (Otospermophilus beecheyi) using transcriptome analyses.

Understanding the molecular basis of adaptations in coevolving species requires identifying the genes that underlie reciprocally selected phenotypes, such as those involved in venom in snakes and resistance to the venom in their prey. In this regard, California ground squirrels (CGS; Otospermophilus beecheyi) are eaten by northern Pacific rattlesnakes (Crotalus oreganus oreganus), but individual squirrels may still show substantial resistance to venom and survive bites. A recent study using proteomics identified venom interactive proteins (VIPs) in the blood serum of CGS. These VIPs represent possible resistance proteins, but the sequences of genes encoding them are unknown despite the value of such data to molecular studies of coevolution. To address this issue, we analyzed a de novo assembled transcriptome from CGS liver tissue-where many plasma proteins are synthesized-and other tissues from this species. We then examined VIP sequences in terms of three characteristics that identify them as possible resistance proteins: evidence for positive selection, high liver expression, and nonsynonymous variation across CGS populations. Based on these characteristics, we identified five VIPs (i.e., α-2-macroglobulin, α-1-antitrypsin-like protein GS55-LT, apolipoprotein A-II, hibernation-associated plasma protein HP-20, and hibernation-associated plasma protein HP-27) as the most likely candidates for resistance proteins among VIPs identified to date. Four of these proteins have been previously implicated in conferring resistance to the venom in mammals, validating our approach. When combined with the detailed information available for rattlesnake venom proteins, these results set the stage for future work focused on understanding coevolutionary interactions at the molecular level between these species.

Low-coverage whole genome sequencing for highly accurate population assignment: Mapping migratory connectivity in the American Redstart (Setophaga ruticilla).

Understanding the geographic linkages among populations across the annual cycle is an essential component for understanding the ecology and evolution of migratory species and for facilitating their effective conservation. While genetic markers have been widely applied to describe migratory connections, the rapid development of new sequencing methods, such as low-coverage whole genome sequencing (lcWGS), provides new opportunities for improved estimates of migratory connectivity. Here, we use lcWGS to identify fine-scale population structure in a widespread songbird, the American Redstart (Setophaga ruticilla), and accurately assign individuals to genetically distinct breeding populations. Assignment of individuals from the nonbreeding range reveals population-specific patterns of varying migratory connectivity. By combining migratory connectivity results with demographic analysis of population abundance and trends, we consider full annual cycle conservation strategies for preserving numbers of individuals and genetic diversity. Notably, we highlight the importance of the Northern Temperate-Greater Antilles migratory population as containing the largest proportion of individuals in the species. Finally, we highlight valuable considerations for other population assignment studies aimed at using lcWGS. Our results have broad implications for improving our understanding of the ecology and evolution of migratory species through conservation genomics approaches.

Inferring population connectivity in eastern massasauga rattlesnakes (Sistrurus catenatus) using landscape genetics.

Assessing the environmental factors that influence the ability of a threatened species to move through a landscape can be used to identify conservation actions that connect isolated populations. However, direct observations of species' movement are often limited, making the development of alternate approaches necessary. Here we use landscape genetic analyses to assess the impact of landscape features on the movement of individuals between local populations of a threatened snake, the eastern massasauga rattlesnake (Sistrurus catenatus). We linked connectivity data with habitat information from two landscapes of similar size: a large region of unfragmented habitat and a previously studied fragmented landscape consisting of isolated patches of habitat. We used this analysis to identify features of the landscape where modification or acquisition would enhance population connectivity in the fragmented region. We found evidence that current connectivity was impacted by both contemporary land-cover features, especially roads, and inherent landscape features such as elevation. Next, we derived estimates of expected movement ability using a recently developed pedigree-based approach and least-cost paths through the unfragmented landscape. We then used our pedigree and resistance map to estimate resistance polygons of the potential extent for S. catenatus movement in the fragmented landscape. These polygons identify possible sites for future corridors connecting currently isolated populations in this landscape by linking the impact of future habitat modification or land acquisition to dispersal ability in this species. Overall, our study shows how modeling landscape resistance across differently fragmented landscapes can identify habitat features that affect contemporary movement in threatened species in fragmented landscapes and how this information can be used to guide mitigation actions whose goal is to connect isolated populations.

The Chemosensory Repertoire of the Eastern Diamondback Rattlesnake (Crotalus adamanteus) Reveals Complementary Genetics of Olfactory and Vomeronasal-Type Receptors.

Pitviper sensory perception incorporates diverse stimuli through the integration of trichromatic color vision, bifocal heat-sensing, and dual-system chemoperception. Chemoperception, or olfaction, is mediated by chemoreceptors in the olfactory bulb and the vomeronasal organ, but the true genomic complexity of the gene families and their relative contributions is unknown. A full genomic accounting of pitviper chemoperception directly complements our current understanding of their venoms by generating a more complete polyphenic representation of their predatory arsenal. To characterize the genetic repertoire of pitviper chemoperception, we analyzed a full-genome assembly for Crotalus adamanteus, the eastern diamondback rattlesnake. We identified hundreds of genes encoding both olfactory receptors (ORs; 362 full-length genes) and type-2 vomeronasal receptors (V2Rs; 430 full-length genes). Many chemoreceptor genes are organized into large tandem repeat arrays. Comparative analysis of V2R orthologs across squamates demonstrates how gene array expansion and contraction underlies the evolution of the chemoreceptor repertoire, which likely reflects shifts in life history traits. Chromosomal assignments of chemosensory genes identified sex chromosome specific chemoreceptor genes, providing gene candidates underlying observed sex-specific chemosensory-based behaviors. We detected widespread episodic evolution in the extracellular, ligand-binding domains of both ORs and V2Rs, suggesting the diversification of chemoreceptors is driven by transient periods of positive selection. We provide a robust genetic framework for studying pitviper chemosensory ecology and evolution.

Genomic signatures of inbreeding and mutation load in a threatened rattlesnake.

Theory predicts that threatened species living in small populations will experience high levels of inbreeding that will increase their genetic load, but recent work suggests that the impact of load may be minimized by purging resulting from long-term population bottlenecks. Empirical studies that examine this idea using genome-wide estimates of inbreeding and genetic load in threatened species are limited. Here we use individual genome resequencing data to compare levels of inbreeding, levels of genetic load (estimated as mutation load) and population history in threatened Eastern massasauga rattlesnakes (Sistrurus catenatus), which exist in small isolated populations, and closely related yet outbred Western massasauga rattlesnakes (Sistrurus tergeminus). In terms of inbreeding, S. catenatus genomes had a greater number of runs of homozygosity of varying sizes, indicating sustained inbreeding through repeated bottlenecks when compared to S. tergeminus. At the species level, outbred S. tergeminus had higher genome-wide levels of mutation load in the form of greater numbers of derived deleterious mutations compared to S. catenatus, presumably due to long-term purging of deleterious mutations in S. catenatus. In contrast, mutations that escaped species-level drift effects within S. catenatus populations were in general more frequent and more often found in homozygous genotypes than in S. tergeminus, suggesting a reduced efficiency of purifying selection in smaller S. catenatus populations for most mutations. Our results support an emerging idea that the historical demography of a threatened species has a significant impact on the type of genetic load present, which impacts implementation of conservation actions such as genetic rescue.

Drift, selection and adaptive variation in small populations of a threatened rattlesnake.

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long- and short-term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre-date the effects of recent drift, and that functional variation in these loci persists despite small short-term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of "drift debt," a nonequilibrium state where present-day levels of variation overestimate the amount of functional genetic diversity present in future populations.

The molecular basis of venom resistance in a rattlesnake-squirrel predator-prey system.

Understanding how interspecific interactions mould the molecular basis of adaptations in coevolving species is a long-sought goal of evolutionary biology. Venom in predators and venom resistance proteins in prey are coevolving molecular phenotypes, and while venoms are highly complex mixtures it is unclear if prey respond with equally complex resistance traits. Here, we use a novel molecular methodology based on protein affinity columns to capture and identify candidate blood serum resistance proteins ("venom interactive proteins" [VIPs]) in California Ground Squirrels (Otospermophilus beecheyi) that interact with venom proteins from their main predator, Northern Pacific Rattlesnakes (Crotalus o. oreganus). This assay showed that serum-based resistance is both population- and species-specific, with serum proteins from ground squirrels showing higher binding affinities for venom proteins of local snakes compared to allopatric individuals. Venom protein specificity assays identified numerous and diverse candidate prey resistance VIPs but also potential targets of venom in prey tissues. Many specific VIPs bind to multiple snake venom proteins and, conversely, single venom proteins bind multiple VIPs, demonstrating that a portion of the squirrel blood serum "resistome" involves broad-based inhibition of nonself proteins and suggests that resistance involves a toxin scavenging mechanism. Analyses of rates of evolution of VIP protein homologues in related mammals show that most of these proteins evolve under purifying selection possibly due to molecular constraints that limit the evolutionary responses of prey to rapidly evolving snake venom proteins. Our method represents a general approach to identify specific proteins involved in co-evolutionary interactions between species at the molecular level.

High levels of functional divergence in toxicity towards prey among the venoms of individual pigmy rattlesnakes.

Venom is a complex molecular phenotype that shows high levels of variation in expressed proteins between individuals within and between populations. However, the functional significance of this variation in terms of toxicity towards prey is largely unknown. Here, we assessed the relative toxicity of venom from individual pygmy rattlesnakes ( Sistrurus miliarius) on brown anoles ( Anolis sagrei) using a novel assay involving tests of fixed doses of venom from individual snakes on individual lizards. We found high levels of functional variation between individual venoms within populations with individual differences (nested within population) explaining 3.6 times more variation in toxicity than population differences. Our results suggest a previously unappreciated adaptive significance to within-population variation in venom. They argue that selective mechanisms that maintain variation within populations may be of equal or greater importance to divergent selection leading to local adaption between populations as evolutionary explanations of venom variation within species.

Local prey community composition and genetic distance predict venom divergence among populations of the northern Pacific rattlesnake (Crotalus oreganus).

Identifying the environmental correlates of divergence in functional traits between populations can provide insights into the evolutionary mechanisms that generate local adaptation. Here, we assess patterns of population differentiation in expressed venom proteins in Northern Pacific rattlesnakes (Crotalus oreganus) from 13 locations across California. We evaluate the relative importance of major biotic (prey species community composition), abiotic (temperature, precipitation and elevation) and genetic factors (genetic distance based on RAD-seq loci) as correlates of population divergence in venom phenotypes. We found that over half of the variation in venom composition is associated with among-population differentiation for genetic and environmental variables and that this variation occurred along axes defining previously observed functional trade-offs between venom proteins that have neurotoxic, myotoxic and hemorrhagic effects. Surprisingly, genetic differentiation among populations was the best predictor of venom divergence, accounting for 46% of overall variation, whereas differences in prey community composition and abiotic factors explained smaller amounts of variation (23% and 19%, respectively). The association between genetic differentiation and venom composition could be due to an isolation-by-distance effect or, more likely, an isolation-by-environment effect where selection against recent migrants is strong, producing a correlation between neutral genetic differentiation and venom differentiation. Our findings suggest that even coarse estimates of prey community composition can be useful in understanding the selection pressures acting on patterns of venom protein expression. Additionally, our results suggest that factors other than adaptation to spatial variation in prey need to be considered when explaining population divergence in venom.

Evaluating local adaptation of a complex phenotype: reciprocal tests of pigmy rattlesnake venoms on treefrog prey.

Theory predicts that predator-prey interactions can generate reciprocal selection pressures on species pairs, which can result in local adaptation, yet the presence and pattern of local adaptation is poorly studied in vertebrate predator-prey systems. Here, we used a reciprocal common garden (laboratory) experimental design involving comparisons between local and foreign populations to determine if local adaptation was present between a generalist predator-the pigmy rattlesnake (Sistrurus miliarius)-and a co-occurring prey-the squirrel treefrog (Hyla squirella). We conducted toxicity trials using snake venom from two populations separated by 340 km tested on prey from sympatric and allopatric populations, resulting in data from four venom origin-frog origin combinations. We assessed venom effectiveness using two measures (frog mortality at 24 h and time to frog death) and then used regression analyses to look for a signal of local adaptation with either measure. We found evidence for local adaptation for one measure (time to death), but not the other (frog mortality). We argue that in this system, the time to death of a prey item is a more ecologically relevant measure of venom effectiveness than is frog mortality at 24 h. Our results document an example of local adaptation between two interacting vertebrates using a whole-organism assay and a local versus foreign criteria and provide evidence that population-level variation in snake venom is adaptive.

Hybrid origin of European Vipers (Vipera magnifica and Vipera orlovi) from the Caucasus determined using genomic scale DNA markers.

BACKGROUND: Studying patterns of introgression can illuminate the role of hybridization in speciation, and help guide decisions relevant to the conservation of rare taxa. Vipera magnifica and Vipera orlovi are small vipers that have high conservation status due to their rarity and restricted distributions in an area of the Caucasus region where two other related species are present - V. kaznakovi and V. renardi. Despite numerous observations of hybridization between different species of small vipers, and the potential of a hybrid origin for V. magnifica and V. orlovi based on their distribution with respect to V. kaznakovi and V. renardi, hypotheses of a hybrid origin have not been formally tested. Here we generate genomic-scale data by performing next generation sequencing of double digest restriction-site associated DNA libraries, and use these multilocus data to test whether these two species are of hybrid origin. RESULTS: We generated over nine hundred loci for 38 specimens of six taxa, and analysed the dataset using Bayesian clustering and multivariate methods, as well as Patterson D-statistics, which can distinguish between incomplete lineage sorting and introgression as explanations for shared polymorphism. The results demonstrate a pattern of historical admixture in the two purported hybrids that is consistent with past gene flow from V. renardi into V. kaznakovi. The average admixture proportion in individuals was low (6.39 %) in the case of V. magnifica, but was higher in V. orlovi (19.02 %). We also show that the specific individual samples used in D-statistic tests can have a significant impact on inferences regarding the magnitude of introgression, suggesting the importance of including multiple individuals in these analyses. CONCLUSIONS: Our results support the conclusion that both V. orlovi and V. magnifica had formed through a hybridization event between V. kaznakovi and V. renardi. Given a low proportion of admixture and absence of clear ecological and morphological differences V. magnifica should be treated as a marginal population of V. kaznakovi. Further studies that include analyses of ecological segregation of V. orlovi from parental taxa and search for evolutionary consequences of hybridisation would clarify if V. orlovi is a distinct hybrid species. Until this we recommend preserving the current taxonomy and protection status of V. orlovi.

Coevolution of venom function and venom resistance in a rattlesnake predator and its squirrel prey.

Measuring local adaptation can provide insights into how coevolution occurs between predators and prey. Specifically, theory predicts that local adaptation in functionally matched traits of predators and prey will not be detected when coevolution is governed by escalating arms races, whereas it will be present when coevolution occurs through an alternate mechanism of phenotype matching. Here, we analyse local adaptation in venom activity and prey resistance across 12 populations of Northern Pacific rattlesnakes and California ground squirrels, an interaction that has often been described as an arms race. Assays of venom function and squirrel resistance show substantial geographical variation (influenced by site elevation) in both venom metalloproteinase activity and resistance factor effectiveness. We demonstrate local adaptation in the effectiveness of rattlesnake venom to overcoming present squirrel resistance, suggesting that phenotype matching plays a role in the coevolution of these molecular traits. Further, the predator was the locally adapted antagonist in this interaction, arguing that rattlesnakes are evolutionarily ahead of their squirrel prey. Phenotype matching needs to be considered as an important mechanism influencing coevolution between venomous animals and resistant prey.

Spread and persistence of influenza A viruses in waterfowl hosts in the North American Mississippi migratory flyway.

UNLABELLED: While geographic distance often restricts the spread of pathogens via hosts, this barrier may be compromised when host species are mobile. Migratory waterfowl in the order Anseriformes are important reservoir hosts for diverse populations of avian-origin influenza A viruses (AIVs) and are assumed to spread AIVs during their annual continental-scale migrations. However, support for this hypothesis is limited, and it is rarely tested using data from comprehensive surveillance efforts incorporating both the temporal and spatial aspects of host migratory patterns. We conducted intensive AIV surveillance of waterfowl using the North American Mississippi Migratory Flyway (MMF) over three autumn migratory seasons. Viral isolates (n = 297) from multiple host species were sequenced and analyzed for patterns of gene dispersal between northern staging and southern wintering locations. Using a phylogenetic and nucleotide identity framework, we observed a larger amount of gene dispersal within this flyway rather than between the other three longitudinally identified North American flyways. Across seasons, we observed patterns of regional persistence of diversity for each genomic segment, along with limited survival of dispersed AIV gene lineages. Reassortment increased with both time and distance, resulting in transient AIV constellations. This study shows that within the MMF, AIV gene flow favors spread along the migratory corridor within a season, and also that intensive surveillance during bird migration is important for identifying virus dispersal on time scales relevant to pandemic responsiveness. In addition, this study indicates that comprehensive monitoring programs to capture AIV diversity are critical for providing insight into AIV evolution and ecology in a major natural reservoir. IMPORTANCE: Migratory birds are a reservoir for antigenic and genetic diversity of influenza A viruses (AIVs) and are implicated in the spread of virus diversity that has contributed to previous pandemic events. Evidence for dispersal of avian-origin AIVs by migratory birds is rarely examined on temporal scales relevant to pandemic or panzootic threats. Therefore, characterizing AIV movement by hosts within a migratory season is important for implementing effective surveillance strategies. We conducted surveillance following birds along a major North American migratory route and observed that within a migratory season, AIVs rapidly reassorted and gene lineages were dispersed primarily within the migratory corridor. Patterns of regional persistence were observed across seasons for each gene segment. We show that dispersal of AIV gene lineages by migratory birds occurs quickly along migratory routes and that surveillance for AIVs threatening human and animal health should focus attention on these routes.