Trait differentiation and modular toxin expression in palm-pitvipers

Background Modularity is the tendency for systems to organize into semi-independent units and can be a key to the evolution and diversification of complex biological systems. Snake venoms are highly variable modular systems that exhibit extreme diversification even across very short time scales. One well-studied venom phenotype dichotomy is a trade-off between neurotoxicity versus hemotoxicity that occurs through the high expression of a heterodimeric neurotoxic phospholipase A2 (PLA2) or snake venom metalloproteinases (SVMPs). We tested whether the variation in these venom phenotypes could occur via variation in regulatory sub-modules through comparative venom gland transcriptomics of representative Black-Speckled Palm-Pitvipers (Bothriechis nigroviridis) and Talamancan Palm-Pitvipers (B. nubestris). Results We assembled 1517 coding sequences, including 43 toxins for B. nigroviridis and 1787 coding sequences including 42 toxins for B. nubestris. The venom gland transcriptomes were extremely divergent between these two species with one B. nigroviridis exhibiting a primarily neurotoxic pattern of expression, both B. nubestris expressing primarily hemorrhagic toxins, and a second B. nigroviridis exhibiting a mixed expression phenotype. Weighted gene coexpression analyses identified six submodules of transcript expression variation, one of which was highly associated with SVMPs and a second which contained both subunits of the neurotoxic PLA2 complex. The sub-module association of these toxins suggest common regulatory pathways underlie the variation in their expression and is consistent with known patterns of inheritance of similar haplotypes in other species. We also find evidence that module associated toxin families show fewer gene duplications and transcript losses between species, but module association did not appear to affect sequence diversification. Conclusion Sub-modular regulation of expression likely contributes to the diversification of venom phenotypes within and among species and underscores the role of modularity in facilitating rapid evolution of complex traits

Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (\u3ci\u3eCrotalus cerastes\u3c/i\u3e) 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

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

The Genetics of Adaptation for Eight Microvirid Bacteriophages

Theories of adaptive molecular evolution have recently experienced significant expansion, and their predictions and assumptions have begun to be subjected to rigorous empirical testing. However, these theories focus largely on predicting the first event in adaptive evolution, the fixation of a single beneficial mutation. To address long-term adaptation it is necessary to include new assumptions, but empirical data are needed for guidance. To empirically characterize the general properties of adaptive walks, eight recently isolated relatives of the single-stranded DNA (ssDNA) bacteriophage φX174 (family Microviridae) were adapted to identical selective conditions. Three of the eight genotypes were adapted in replicate, for a total of 11 adaptive walks. We measured fitness improvement and identified the genetic changes underlying the observed adaptation. Nearly all phages were evolvable; nine of the 11 lineages showed a significant increase in fitness. However, fitness plateaued quickly, and adaptation was achieved through only three substitutions on average. Parallel evolution was rampant, both across replicates of the same genotype as well as across different genotypes, yet adaptation of replicates never proceeded through the exact same set of mutations. Despite this, final fitnesses did not vary significantly among replicates. Final fitnesses did vary significantly across genotypes but not across phylogenetic groupings of genotypes. A positive correlation was found between the number of substitutions in an adaptive walk and the magnitude of fitness improvement, but no correlation was found between starting and ending fitness. These results provide an empirical framework for future adaptation theory

Evidence for divergent patterns of local selection driving venom variation in Mojave Rattlesnakes (\u3ci\u3eCrotalus scutulatus\u3c/i\u3e)

Snake venoms represent an enriched system for investigating the evolutionary processes that lead to complex and dynamic trophic adaptations. It has long been hypothesized that natural selection may drive geographic variation in venom composition, yet previous studies have lacked the population genetic context to examine these patterns. We leverage range-wide sampling of Mojave Rattlesnakes (Crotalus scutulatus) and use a combination of venom, morphological, phylogenetic, population genetic, and environmental data to characterize the striking dichotomy of neurotoxic (Type A) and hemorrhagic (Type B) venoms throughout the range of this species. We find that three of the four previously identified major lineages within C. scutulatus possess a combination of Type A, Type B, and a ‘mixed’ Type A + B venom phenotypes, and that fixation of the two main venom phenotypes occurs on a more fine geographic scale than previously appreciated. We also find that Type A + B individuals occur in regions of inferred introgression, and that this mixed phenotype is comparatively rare. Our results support strong directional local selection leading to fixation of alternative venom phenotypes on a fine geographic scale, and are inconsistent with balancing selection to maintain both phenotypes within a single population. Our comparisons to biotic and abiotic factors further indicate that venom phenotype correlates with fang morphology and climatic variables. We hypothesize that links to fang morphology may be indicative of co-evolution of venom and other trophic adaptations, and that climatic variables may be linked to prey distributions and/or physiology, which in turn impose selection pressures on snake venoms

Identification and Characterization of Novel Proteins from Arizona Bark Scorpion Venom That Inhibit Nav1.8, a Voltage-Gated Sodium Channel Regulator of Pain Signaling

The voltage-gated sodium channel Nav1.8 is linked to neuropathic and inflammatory pain, highlighting the potential to serve as a drug target. However, the biophysical mechanisms that regulate Nav1.8 activation and inactivation gating are not completely understood. Progress has been hindered by a lack of biochemical tools for examining Nav1.8 gating mechanisms. Arizona bark scorpion (Centruroides sculpturatus) venom proteins inhibit Nav1.8 and block pain in grasshopper mice (Onychomys torridus). These proteins provide tools for examining Nav1.8 structure-activity relationships. To identify proteins that inhibit Nav1.8 activity, venom samples were fractioned using liquid chromatography (reversed-phase and ion exchange). A recombinant Nav1.8 clone expressed in ND7/23 cells was used to identify subfractions that inhibited Nav1.8 Na+ current. Mass-spectrometry-based bottom-up proteomic analyses identified unique peptides from inhibitory subfractions. A search of the peptides against the AZ bark scorpion venom gland transcriptome revealed four novel proteins between 40 and 60% conserved with venom proteins from scorpions in four genera (Centruroides, Parabuthus, Androctonus, and Tityus). Ranging from 63 to 82 amino acids, each primary structure includes eight cysteines and a "CXCE" motif, where X = an aromatic residue (tryptophan, tyrosine, or phenylalanine). Electrophysiology data demonstrated that the inhibitory effects of bioactive subfractions can be removed by hyperpolarizing the channels, suggesting that proteins may function as gating modifiers as opposed to pore blockers

Report from the First Snake Genomics and Integrative Biology Meeting

This report summarizes the proceedings of the 1st Snake Genomics and Integrative Biology Meeting held in Vail, CO USA, 5-8 October 2011. The meeting had over twenty registered participants, and was conducted as a single session of presentations. Goals of the meeting included coordination of genomic data collection and fostering collaborative interactions among researchers using snakes as model systems.Organismic and Evolutionary Biolog

Report from the First Snake Genomics and Integrative Biology Meeting

This report summarizes the proceedings of the 1st Snake Genomics and Integrative Biology Meeting held in Vail, CO USA, 5-8 October 2011. The meeting had over twenty registered participants, and was conducted as a single session of presentations. Goals of the meeting included coordination of genomic data collection and fostering collaborative interactions among researchers using snakes as model systems

Distributed under Creative Commons CC-BY 4.0 Viral recombination blurs taxonomic lines: examination of single-stranded DNA viruses in a wastewater treatment plant

ABSTRACT Understanding the structure and dynamics of microbial communities, especially those of economic concern, is of paramount importance to maintaining healthy and efficient microbial communities at agricultural sites and large industrial cultures, including bioprocessors. Wastewater treatment plants are large bioprocessors which receive water from multiple sources, becoming reservoirs for the collection of many viral families that infect a broad range of hosts. To examine this complex collection of viruses, full-length genomes of circular ssDNA viruses were isolated from a wastewater treatment facility using a combination of sucrose-gradient size selection and rolling-circle amplification and sequenced on an Illumina MiSeq. Single-stranded DNA viruses are among the least understood groups of microbial pathogens due to genomic biases and culturing difficulties, particularly compared to the larger, more often studied dsDNA viruses. However, the group contains several notable well-studied examples, including agricultural pathogens which infect both livestock and crops (Circoviridae and Geminiviridae), and model organisms for genetics and evolution studies (Microviridae). Examination of the collected viral DNA provided evidence for 83 unique genotypic groupings, which were genetically dissimilar to known viral types and exhibited broad diversity within the community. Furthermore, although these genomes express similarities to known viral families, such as Circoviridae, Geminiviridae, and Microviridae, many are so divergent that they may represent new taxonomic groups. This study demonstrated the efficacy of the protocol for separating bacteria and large viruses from the sought after ssDNA viruses and the ability to use this protocol to obtain an in-depth analysis of the diversity within this group