The venom-gland transcriptome of the eastern diamondback rattlesnake (Crotalus adamanteus)

Background

Snake venoms have significant impacts on human populations through the morbidity and mortality associated with snakebites and as sources of drugs, drug leads, and physiological research tools. Genes expressed by venom-gland tissue, including those encoding toxic proteins, have therefore been sequenced but only with relatively sparse coverage resulting from the low-throughput sequencing approaches available. High-throughput approaches based on 454 pyrosequencing have recently been applied to the study of snake venoms to give the most complete characterizations to date of the genes expressed in active venom glands, but such approaches are costly and still provide a far-from-complete characterization of the genes expressed during venom production.

Results

We describe the de novo assembly and analysis of the venom-gland transcriptome of an eastern diamondback rattlesnake (Crotalus adamanteus) based on 95,643,958 pairs of quality-filtered, 100-base-pair Illumina reads. We identified 123 unique, full-length toxin-coding sequences, which cluster into 78 groups with less than 1% nucleotide divergence, and 2,879 unique, full-length nontoxin coding sequences. The toxin sequences accounted for 35.4% of the total reads, and the nontoxin sequences for an additional 27.5%. The most highly expressed toxin was a small myotoxin related to crotamine, which accounted for 5.9% of the total reads. Snake-venom metalloproteinases accounted for the highest percentage of reads mapping to a toxin class (24.4%), followed by C-type lectins (22.2%) and serine proteinases (20.0%). The most diverse toxin classes were the C-type lectins (21 clusters), the snake-venom metalloproteinases (16 clusters), and the serine proteinases (14 clusters). The high-abundance nontoxin transcripts were predominantly those involved in protein folding and translation, consistent with the protein-secretory function of the tissue.

Conclusions

We have provided the most complete characterization of the genes expressed in an active snake venom gland to date, producing insights into snakebite pathology and guidance for snakebite treatment for the largest rattlesnake species and arguably the most dangerous snake native to the United States of America, C. adamanteus. We have more than doubled the number of sequenced toxins for this species and created extensive genomic resources for snakes based entirely on de novo assembly of Illumina sequence data.

An empirical test of the mutational landscape model of adaptation using a single-stranded DNA virus

The primary impediment to formulating a general theory for adaptive evolution has been the unknown distribution of fitness effects for new beneficial mutations. By applying extreme value theory, Gillespie circumvented this issue in his mutational landscape model for the adaptation of DNA sequences, and Orr recently extended Gillespie's model, generating testable predictions regarding the course of adaptive evolution. Here we provide the first empirical examination of this model, using a single-stranded DNA bacteriophage related to phiX174, and find that our data are consistent with Orr's predictions, provided that the model is adjusted to incorporate mutation bias. Orr's work suggests that there may be generalities in adaptive molecular evolution that transcend the biological details of a system, but we show that for the model to be useful as a predictive or inferential tool, some adjustments for the biology of the system will be necessary.

The venom-gland transcriptome of the eastern coral snake (Micrurus fulvius) reveals high venom complexity in the intragenomic evolution of venoms

BACKGROUND

Snake venom is shaped by the ecology and evolution of venomous species, and signals of positive selection in toxins have been consistently documented, reflecting the role of venoms as an ecologically critical phenotype. New World coral snakes (Elapidae) are represented by three genera and over 120 species and subspecies that are capable of causing significant human morbidity and mortality, yet coral-snake venom composition is poorly understood in comparison to that of Old World elapids. High-throughput sequencing is capable of identifying thousands of loci, while providing characterizations of expression patterns and the molecular evolutionary forces acting within the venom gland.

RESULTS

We describe the de novo assembly and analysis of the venom-gland transcriptome of the eastern coral snake (Micrurus fulvius). We identified 1,950 nontoxin transcripts and 116 toxin transcripts. These transcripts accounted for 57.1% of the total reads, with toxins accounting for 45.8% of the total reads. Phospholipases A(2) and three-finger toxins dominated expression, accounting for 86.0% of the toxin reads. A total of 15 toxin families were identified, revealing venom complexity previously unknown from New World coral snakes. Toxins exhibited high levels of heterozygosity relative to nontoxins, and overdominance may favor gene duplication leading to the fixation of advantageous alleles. Phospholipase A(2) expression was uniformly distributed throughout the class while three-finger toxin expression was dominated by a handful of transcripts, and phylogenetic analyses indicate that toxin divergence may have occurred following speciation. Positive selection was detected in three of the four most diverse toxin classes, suggesting that venom diversification is driven by recurrent directional selection.

CONCLUSIONS

We describe the most complete characterization of an elapid venom gland to date. Toxin gene duplication may be driven by heterozygote advantage, as the frequency of polymorphic toxin loci was significantly higher than that of nontoxins. Diversification among toxins appeared to follow speciation reflecting species-specific adaptation, and this divergence may be directly related to dietary shifts and is suggestive of a coevolutionary arms race.

A high-throughput venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) and evidence for pervasive positive selection across toxin classes

Despite causing considerable human mortality and morbidity, animal toxins represent a valuable source of pharmacologically active macromolecules, a unique system for studying molecular adaptation, and a powerful framework for examining structure-function relationships in proteins. Snake venoms are particularly useful in the latter regard as they consist primarily of a moderate number of proteins and peptides that have been found to belong to just a handful of protein families. As these proteins and peptides are produced in dedicated glands, transcriptome sequencing has proven to be an effective approach to identifying the expressed toxin genes. We generated a venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) using Roche 454 sequencing technology. In the current work, we focus on transcripts encoding toxins. We identified 40 unique toxin transcripts, 30 of which have full-length coding sequences, and 10 have only partial coding sequences. These toxins account for 24% of the total sequencing reads. We found toxins from 11 previously described families of snake-venom toxins and have discovered two putative, previously undescribed toxin classes. The most diverse and highly expressed toxin classes in the C. adamanteus venom-gland transcriptome are the serine proteinases, metalloproteinases, and C-type lectins. The serine proteinases are the most abundant class, accounting for 35% of the toxin sequencing reads. Metalloproteinases are the most diverse; 11 different forms have been identified. Using our sequences and those available in public databases, we detected positive selection in seven of the eight toxin families for which sufficient sequences were available for the analysis. We find that the vast majority of the genes that contribute directly to this vertebrate trait show evidence for a role for positive selection in their evolutionary history.

Epistasis between beneficial mutations and the phenotype-to-fitness Map for a ssDNA virus

Epistatic interactions between genes and individual mutations are major determinants of the evolutionary properties of genetic systems and have therefore been well documented, but few quantitative data exist on epistatic interactions between beneficial mutations, presumably because such mutations are so much rarer than deleterious ones. We explored epistasis for beneficial mutations by constructing genotypes with pairs of mutations that had been previously identified as beneficial to the ssDNA bacteriophage ID11 and by measuring the effects of these mutations alone and in combination. We constructed 18 of the 36 possible double mutants for the nine available beneficial mutations. We found that epistatic interactions between beneficial mutations were all antagonistic—the effects of the double mutations were less than the sums of the effects of their component single mutations. We found a number of cases of decompensatory interactions, an extreme form of antagonistic epistasis in which the second mutation is actually deleterious in the presence of the first. In the vast majority of cases, recombination uniting two beneficial mutations into the same genome would not be favored by selection, as the recombinant could not outcompete its constituent single mutations. In an attempt to understand these results, we developed a simple model in which the phenotypic effects of mutations are completely additive and epistatic interactions arise as a result of the form of the phenotype-to-fitness mapping. We found that a model with an intermediate phenotypic optimum and additive phenotypic effects provided a good explanation for our data and the observed patterns of epistatic interactions.

Epistasis, the extent to which the effects of a mutation depend on its genetic context, can have profound effects on the evolutionary process and strongly affects our understanding of the prevalence of sexual reproduction. It has been investigated in a diverse array of organisms but almost exclusively for deleterious mutations. Interactions between beneficial mutations can impede adaptation, and we therefore investigated epistasis between beneficial mutations by constructing 18 bacteriophage genomes, each with two mutations that had been previously identified as beneficial, and measuring their fitnesses. We found universal evidence for epistasis—every pair of mutations conferred fitness lower than that expected from the single mutations alone. In many cases, a beneficial mutation became deleterious when in combination with another, and in fact, only one pair out of 18 could be shown to confer significantly greater fitness than its constituent mutations alone. To explain these results, we developed a model of the relationship between phenotype and fitness that posits an intermediate phenotypic optimum and assumes no epistasis at the phenotypic level. This model fit our data well and showed that the patterns we observed could result because mutants have phenotypes that overshoot the optimum.

Horizontal gene transfer and the evolution of microvirid coliphage genomes

Bacteriophage genomic evolution has been largely characterized by rampant, promiscuous horizontal gene transfer involving both homologous and nonhomologous source DNA. This pattern has emerged through study of the tailed double-stranded DNA (dsDNA) phages and is based upon a sparse sampling of the enormous diversity of these phages. The single-stranded DNA phages of the family Microviridae, including phiX174, appear to evolve through qualitatively different mechanisms, possibly as result of their strictly lytic lifestyle and small genome size. However, this apparent difference could reflect merely a dearth of relevant data. We sought to characterize the forces that contributed to the molecular evolution of the Microviridae and to examine the genetic structure of this single family of bacteriophage by sequencing the genomes of microvirid phage isolated on a single bacterial host. Microvirids comprised 3.5% of the detectable phage in our environmental samples, and sequencing yielded 42 new microvirid genomes. Phylogenetic analysis of the genes contained in these and five previously described microvirid phages identified three distinct clades and revealed at least two horizontal transfer events between clades. All members of one clade have a block of five putative genes that are not present in any member of the other two clades. Our data indicate that horizontal transfer does contribute to the evolution of the microvirids but is both quantitatively and qualitatively different from what has been observed for the dsDNA phages.

Beneficial fitness effects are not exponential for two viruses

The distribution of fitness effects for beneficial mutations is of paramount importance in determining the outcome of adaptation. It is generally assumed that fitness effects of beneficial mutations follow an exponential distribution, for example, in theoretical treatments of quantitative genetics, clonal interference, experimental evolution, and the adaptation of DNA sequences. This assumption has been justified by the statistical theory of extreme values, because the fitnesses conferred by beneficial mutations should represent samples from the extreme right tail of the fitness distribution. Yet in extreme value theory, there are three different limiting forms for right tails of distributions, and the exponential describes only those of distributions in the Gumbel domain of attraction. Using beneficial mutations from two viruses, we show for the first time that the Gumbel domain can be rejected in favor of a distribution with a right-truncated tail, thus providing evidence for an upper bound on fitness effects. Our data also violate the common assumption that small-effect beneficial mutations greatly outnumber those of large effect, as they are consistent with a uniform distribution of beneficial effects.

Quantity, Not Quality: Rapid Adaptation in a Polygenic Trait Proceeded Exclusively through Expression Differentiation

A trait's genomic architecture can affect the rate and mechanism of adaptation, and although many ecologically-important traits are polygenic, most studies connecting genotype, phenotype, and fitness in natural populations have focused on traits with relatively simple genetic bases. To understand the genetic basis of polygenic adaptation, we must integrate genomics, phenotypic data, ecology, and fitness effects for a genetically tractable, polygenic trait; snake venoms provide such a system for studying polygenic adaptation because of their genetic tractability and vital ecological role in feeding and defense. We used a venom transcriptome-proteome map, quantitative proteomics, genomics, and fitness assays in sympatric prey to construct a genotype-phenotype-fitness map for the venoms of an island-mainland pair of rattlesnake populations. Reciprocal fitness experiments demonstrated that each population was locally adapted to sympatric prey. We identified significant expression differentiation with little to no coding-sequence variation across populations, demonstrating that expression differentiation was exclusively the genetic basis of polygenic adaptation. Previous research on the genetics of adaptation, however, has largely been biased toward investigating protein-coding regions because of the complexity of gene regulation. Our results showed that biases at the molecular level can be in the opposite direction, highlighting the need for more systematic comparisons of different molecular mechanisms underlying rapid, adaptive evolution in polygenic traits.

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.

RNA-seq and high-definition mass spectrometry reveal the complex and divergent venoms of two rear-fanged colubrid snakes

Background

Largely because of their direct, negative impacts on human health, the venoms of front-fanged snakes of the families Viperidae and Elapidae have been extensively characterized proteomically, transcriptomically, and pharmacologically. However, relatively little is known about the molecular complexity and evolution of the venoms of rear-fanged colubrid snakes, which are, with a few notable exceptions, regarded as harmless to humans. Many of these snakes have venoms with major effects on their preferred prey, and their venoms are probably as critical to their survival as those of front-fanged elapids and viperids.

Results

We sequenced the venom-gland transcriptomes from a specimen of Hypsiglena (Desert Night Snake; family Colubridae, subfamily Dipsadinae) and of Boiga irregularis (Brown Treesnake; family Colubridae, subfamily Colubrinae) and verified the transcriptomic results proteomically by means of high-definition mass spectrometry. We identified nearly 3,000 nontoxin genes for each species. For B. irregularis, we found 108 putative toxin transcripts in 46 clusters with <1% nucleotide divergence, and for Hypsiglena we identified 79 toxin sequences that were grouped into 33 clusters. Comparisons of the venoms revealed divergent venom types, with Hypsiglena possessing a viper-like venom dominated by metalloproteinases, and B. irregularis having a more elapid-like venom, consisting primarily of three-finger toxins.

Conclusions

Despite the difficulty of procuring venom from rear-fanged species, we were able to complete all analyses from a single specimen of each species without pooling venom samples or glands, demonstrating the power of high-definition transcriptomic and proteomic approaches. We found a high level of divergence in the venom types of two colubrids. These two venoms reflected the hemorrhagic/neurotoxic venom dichotomy that broadly characterizes the difference in venom strategies between elapids and viperids.

The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network.

Phenotypic integration in the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus)

Selection can vary geographically across environments and temporally over the lifetime of an individual. Unlike geographic contexts, where different selective regimes can act on different alleles, age-specific selection is constrained to act on the same genome by altering age-specific expression. Snake venoms are exceptional traits for studying ontogeny because toxin expression variation directly changes the phenotype; relative amounts of venom components determine, in part, venom efficacy. Phenotypic integration is the dependent relationship between different traits that collectively produce a complex phenotype and, in venomous snakes, may include traits as diverse as venom, head shape and fang length. We examined the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus) across environments and over the lifetime of individuals and used a genotype-phenotype map approach, protein expression data and morphological data to demonstrate that: (i) ontogenetic effects explained more of the variation in toxin expression variation than geographic effects, (ii) both juveniles and adults varied geographically, (iii) toxin expression variation was a result of directional selection and (iv) different venom phenotypes covaried with morphological traits also associated with feeding in temporal (ontogenetic) and geographic (functional) contexts. These data are the first to demonstrate, to our knowledge, phenotypic integration between multiple morphological characters and a biochemical phenotype across populations and age classes. We identified copy number variation as the mechanism driving the difference in the venom phenotype associated with these morphological differences, and the parallel mitochondrial, venom and morphological divergence between northern and southern clades suggests that each clade may warrant classification as a separate evolutionarily significant unit.

Testing the extreme value domain of attraction for distributions of beneficial fitness effects

In modeling evolutionary genetics, it is often assumed that mutational effects are assigned according to a continuous probability distribution, and multiple distributions have been used with varying degrees of justification. For mutations with beneficial effects, the distribution currently favored is the exponential distribution, in part because it can be justified in terms of extreme value theory, since beneficial mutations should have fitnesses in the extreme right tail of the fitness distribution. While the appeal to extreme value theory seems justified, the exponential distribution is but one of three possible limiting forms for tail distributions, with the other two loosely corresponding to distributions with right-truncated tails and those with heavy tails. We describe a likelihood-ratio framework for analyzing the fitness effects of beneficial mutations, focusing on testing the null hypothesis that the distribution is exponential. We also describe how to account for missing the smallest-effect mutations, which are often difficult to identify experimentally. This technique makes it possible to apply the test to gain-of-function mutations, where the ancestral genotype is unable to grow under the selective conditions. We also describe how to pool data across experiments, since we expect few possible beneficial mutations in any particular experiment.

Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (Crotalus cerastes) lineages reveal little differential expression despite individual variation

Changes in gene expression can rapidly influence adaptive traits in the early stages of lineage diversification. Venom is an adaptive trait comprised of numerous toxins used for prey capture and defense. Snake venoms can vary widely between conspecific populations, but the influence of lineage diversification on such compositional differences are unknown. To explore venom differentiation in the early stages of lineage diversification, we used RNA-seq and mass spectrometry to characterize Sidewinder Rattlesnake (Crotalus cerastes) venom. We generated the first venom-gland transcriptomes and complementary venom proteomes for eight individuals collected across the United States and tested for expression differences across life history traits and between subspecific, mitochondrial, and phylotranscriptomic hypotheses. Sidewinder venom was comprised primarily of hemorrhagic toxins, with few cases of differential expression attributable to life history or lineage hypotheses. However, phylotranscriptomic lineage comparisons more than doubled instances of significant expression differences compared to all other factors. Nevertheless, only 6.4% of toxins were differentially expressed overall, suggesting that shallow divergence has not led to major changes in Sidewinder venom composition. Our results demonstrate the need for consensus venom-gland transcriptomes based on multiple individuals and highlight the potential for discrepancies in differential expression between different phylogenetic hypotheses.

Properties of adaptive walks on uncorrelated landscapes under strong selection and weak mutation

We examine properties of adaptive walks on uncorrelated (i.e. random) fitness landscapes starting from moderately fit genotypes under strong selection weak mutation. As an extension of Orr's model for a single step in an adaptive walk under these conditions, we show that the fitness rank of the dominant genotype in a population after the fixation of a beneficial mutation is, on average, (i+6)/4, where i is the fitness rank of the starting genotype. This accounts for the change in rank due to acquiring a new set of single-mutation neighbors after fixing a new allele through natural selection. Under this scenario, adaptive walks can be modeled as a simple Markov chain on the space of possible fitness ranks with an absorbing state at i = 1, from which no beneficial mutations are accessible. We find that these walks are typically short and are often completed in a single step when starting from a moderately fit genotype. As in Orr's original model, these results are insensitive to both the distribution of fitness effects and most biological details of the system under consideration.

The genetics of venom ontogeny in the eastern diamondback rattlesnake (Crotalus adamanteus)

The same selective forces that give rise to rapid inter- and intraspecific divergence in snake venoms can also favor differences in venoms across life-history stages. Ontogenetic changes in venom composition are well known and widespread in snakes but have not been investigated to the level of unambiguously identifying the specific loci involved. The eastern diamondback rattlesnake was previously shown to undergo an ontogenetic shift in venom composition at sexual maturity, and this shift accounted for more venom variation than geography. To characterize the genetics underlying the ontogenetic venom compositional change in C. adamanteus, we sequenced adult/juvenile pairs of venom-gland transcriptomes from five populations previously shown to have different adult venom compositions. We identified a total of 59 putative toxin transcripts for C. adamanteus, and 12 of these were involved in the ontogenetic change. Three toxins were downregulated, and nine were upregulated in adults relative to juveniles. Adults and juveniles expressed similar total levels of snake-venom metalloproteinases but differed substantially in their featured paralogs, and adults expressed higher levels of Bradykinin-potentiating and C-type natriuretic peptides, nerve growth factor, and specific paralogs of phospholipases A₂ and snake venom serine proteinases. Juvenile venom was more toxic to mice, indicating that the expression differences resulted in a phenotypically, and therefore potentially ecologically, significant difference in venom function. We also showed that adult and juvenile venom-gland transcriptomes for a species with known ontogenetic venom variation were equally effective at individually providing a full characterization of the venom genes of a species but that any particular individual was likely to lack several toxins in their transcriptome. A full characterization of a species’ venom-gene complement therefore requires sequencing more than one individual, although the ages of the individuals are unimportant.

Mutation-Driven Parallel Evolution during Viral Adaptation

Convergent evolution has been demonstrated across all levels of biological organization, from parallel nucleotide substitutions to convergent evolution of complex phenotypes, but whether instances of convergence are the result of selection repeatedly finding the same optimal solution to a recurring problem or are the product of mutational biases remains unsettled. We generated 20 replicate lineages allowed to fix a single mutation from each of four bacteriophage genotypes under identical selective regimes to test for parallel changes within and across genotypes at the levels of mutational effect distributions and gene, protein, amino acid, and nucleotide changes. All four genotypes shared a distribution of beneficial mutational effects best approximated by a distribution with a finite upper bound. Parallel adaptation was high at the protein, gene, amino acid, and nucleotide levels, both within and among phage genotypes, with the most common first-step mutation in each background fixing on an average in 7 of 20 replicates and half of the substitutions in two of the four genotypes occurring at shared sites. Remarkably, the mutation of largest beneficial effect that fixed for each genotype was never the most common, as would be expected if parallelism were driven by selection. In fact, the mutation of smallest benefit for each genotype fixed in a total of 7 of 80 lineages, equally as often as the mutation of largest benefit, leading us to conclude that adaptation was largely mutation-driven, such that mutational biases led to frequent parallel fixation of mutations of suboptimal effect.

Phenotypic Variation in Mojave Rattlesnake (Crotalus scutulatus) Venom Is Driven by Four Toxin Families

Phenotypic diversity generated through altered gene expression is a primary mechanism facilitating evolutionary response in natural systems. By linking the phenotype to genotype through transcriptomics, it is possible to determine what changes are occurring at the molecular level. High phenotypic diversity has been documented in rattlesnake venom, which is under strong selection due to its role in prey acquisition and defense. Rattlesnake venom can be characterized by the presence (Type A) or absence (Type B) of a type of neurotoxic phospholipase A2 (PLA2), such as Mojave toxin, that increases venom toxicity. Mojave rattlesnakes (Crotalus scutulatus), represent this diversity as both venom types are found within this species and within a single panmictic population in the Sonoran Desert. We used comparative venom gland transcriptomics of nine specimens of C. scutulatus from this region to test whether expression differences explain diversity within and between venom types. Type A individuals expressed significantly fewer toxins than Type B individuals owing to the diversity of C-type lectins (CTLs) and snake venom metalloproteinases (SVMPs) found in Type B animals. As expected, both subunits of Mojave toxin were exclusively found in Type A individuals but we found high diversity in four additional PLA2s that was not associated with a venom type. Myotoxin a expression and toxin number variation was not associated with venom type, and myotoxin a had the highest range of expression of any toxin class. Our study represents the most comprehensive transcriptomic profile of the venom type dichotomy in rattlesnakes and C. scutulatus. Even intra-specifically, Mojave rattlesnakes showcase the diversity of snake venoms and illustrate that variation within venom types blurs the distinction of the venom dichotomy.

Early significant ontogenetic changes in snake venoms

Snake venom plays a critical role in food acquisition, digestion, and defense. Venoms are known to change throughout the life of some snake species, but nothing is known about the venom composition of hatchling/neonate snakes prior to and just after their first shedding cycle, despite this being a critical time in the life of the snake. Using a cohort of Crotalus horridus and two cohorts of Crotalus adamanteus, we showed for the first time that snakes undergo significant changes in venom composition after the postnatal shedding event. The number of changes among cohorts ranged widely and there was wide variation in the direction of protein regulation, which appeared to be on a locus-specific level rather than protein-family level. These significant venom composition changes that take place in the first few weeks of life most likely play critical roles in venom economy and resource conservation and may partially explain the rare, post-birth maternal care found in some venomous species.

Functional characterizations of venom phenotypes in the eastern diamondback rattlesnake (Crotalus adamanteus) and evidence for expression-driven divergence in toxic activities among populations

Phenotypes frequently vary across and within species. The connection between specific phenotypic effects and function, however, is less understood despite being essential to our understanding of the adaptive process. Snake venoms are ideal for identifying functionally important phenotypic variation because venom variation is common, and venoms can be functionally characterized through simple assays and toxicity measurements. Previous work with the eastern diamondback rattlesnake (Crotalus adamanteus) used multivariate statistical approaches to identify six unique venom phenotypes. We functionally characterized hemolytic, gelatinase, fibrinogenolytic, and coagulant activity for all six phenotypes, as well as one additional venom, to determine if the statistically significant differences in toxin expression levels previously documented corresponded to differences in venom activity. In general, statistical differences in toxin expression predicted the identified functional differences, or lack thereof, in toxic activity, demonstrating that the statistical approach used to characterize C. adamanteus venoms was a fair representation of biologically meaningful differences. Minor differences in activity not accounted for by the statistical model may be the result of amino-acid differences and/or post-translational modifications, but overall we were able to link variation in protein expression levels to variation in function as predicted by multivariate statistical approaches.