Diverse Gene Regulatory Mechanisms Alter Rattlesnake Venom Gene Expression at Fine Evolutionary Scales

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

A snapshot of Bothrops jararaca snake venom gland subcellular proteome

Snake venom protein synthesis undergoes finely regulated processes in the specialized secretory epithelium within the venom gland. Such processes occur within a defined period in the cell and at specific cellular locations. Thus, the determination of subcellular proteomes allows the characterization of protein groups for which the site may be relevant to their biological roles, thereby allowing the deconvolution of complex biological circuits into functional information. In this regard, we performed subcellular fractionation of proteins from B. jararaca venom gland, focusing on nuclear proteins since this cellular compartment comprises key effectors that shape gene expression. Our results provided a snapshot of B. jararaca's subcellular venom gland proteome and pointed to a 'conserved' proteome core among different life stages (newborn and adult) and between sexes (adult male and female). Overall, the top 15 highly abundant proteins identified in B. jararaca venom glands mirrored the panel of highly expressed genes in human salivary glands. Therefore, the expression profile observed for such a protein set could be considered a conserved core signature of salivary gland secretory epithelium. Moreover, the newborn venom gland displayed a unique expression signature of transcription factors involved in regulating transcription and biosynthetic processes and may mirror biological constraints of the ontogenetic development of B. jararaca, contributing to venom proteome diversity.

Snake venom gene expression is coordinated by novel regulatory architecture and the integration of multiple co-opted vertebrate pathways

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis-regulatory sequences (enhancers and promoters), trans-regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis-regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade.

Venom production and secretion in reptiles

The venom glands of reptiles, particularly those of front-fanged advanced snakes, must satisfy conflicting biological demands: rapid synthesis of potentially labile and highly toxic proteins, storage in the gland lumen for long periods, stabilization of the stored secretions, immediate activation of toxins upon deployment and protection of the animal from the toxic effects of its own venom. This dynamic system could serve as a model for the study of a variety of different phenomena involving exocrine gland activation, protein synthesis, stabilization of protein products and secretory mechanisms. However, these studies have been hampered by a lack of a long-term model that can be propagated in the lab (as opposed to whole-animal studies). Numerous attempts have been made to extend the lifetime of venom gland secretory cells, but only recently has an organoid model been shown to have the requisite qualities of recapitulation of the native system, self-propagation and long-term viability (>1 year). A tractable model is now available for myriad cell- and molecular-level studies of venom glands, protein synthesis and secretion. However, venom glands of reptiles are not identical, and many differ very extensively in overall architecture, microanatomy and protein products produced. This Review summarizes the similarities among and differences between venom glands of helodermatid lizards and of rear-fanged and front-fanged snakes, highlighting those areas that are well understood and identifying areas where future studies can fill in significant gaps in knowledge of these ancient, yet fascinating systems.

Contrasting patterns of venom regeneration in a centipede (Scolopendra viridis) and a scorpion (Centruroides hentzi)

As biochemical traits with clear fitness consequences, venoms serve a critical ecological role for the animals that produce them. Understanding how venoms are maintained and regenerated after use will, therefore, provide valuable insight into the ecology of venomous animals. Furthermore, most studies on venomous organisms often require removing animals from the wild and waiting extended periods of time between venom extractions. Uncovering the patterns of venom regeneration across different species will likely lead to the development of more efficient venom extraction protocols, reducing both experimental time and the number of animals required. Using reversed-phase high-performance liquid chromatography, we identified asynchronous regeneration of venom protein component abundances in the centipede Scolopendra viridis but found no evidence for asynchronous venom regeneration in the scorpion Centruroides hentzi. We also observed high levels of intraspecific venom variation in C. hentzi, emphasizing the importance of testing for intraspecific venom variation in studies evaluating the synchronicity of venom regeneration. Although the regeneration of relative venom protein component abundances is an asynchronous process in S. viridis, we provide evidence that the presence-absence of major venom components is not an asynchronous process and suggest that studies relying on just the presence/absence of individual proteins (e.g. bioprospecting, drug discovery) could use catch-and-release methods of venom extraction to reduce the number of animals removed from the wild.

Experimental Manipulation of Corticosterone Does Not Affect Venom Composition or Functional Activity in Free-Ranging Rattlesnakes

AbstractVenom is an integral feeding trait in many animal species. Although venom often varies ontogenetically, little is known about the proximate physiological mediators of venom variation within individuals. The glucocorticoid hormone corticosterone (CORT) can alter the transcription and activation of proteins, including homologues of snake venom components such as snake venom metalloproteinases (SVMPs) and phospholipase A₂ (PLA₂). CORT is endogenously produced by snakes, varies seasonally and also in response to stress, and is a candidate endogenous mediator of changes in venom composition and functional activity. Here, we tested the hypothesis that CORT induces changes in snake venom by sampling the venom of wild adult rattlesnakes before and after they were treated with either empty (control) or CORT-filled (treatment) Silastic implants. We measured longitudinal changes in whole-venom composition, whole-venom total protein content, and enzymatic activity of SVMP and PLA₂ components of venom. We also assessed the within-individual repeatability of venom components. Despite successfully elevating plasma CORT in the treatment group, we found no effect of CORT treatment or average plasma CORT level on any venom variables measured. Except for total protein content, venom components were highly repeatable within individuals ([Formula: see text]). Our results indicate that the effects of CORT, a hormone commonly associated with stress and metabolic functions, in adult rattlesnake venom are negligible. Our findings bode well for venom researchers and biomedical applications that rely on the consistency of venoms produced from potentially stressed individuals and provide an experimental framework for future studies of proximate mediators of venom variation across an individual's life span.

Toxinology provides multidirectional and multidimensional opportunities: A personal perspective

In nature, toxins have evolved as weapons to capture and subdue the prey or to counter predators or competitors. When they are inadvertently injected into humans, they cause symptoms ranging from mild discomfort to debilitation and death. Toxinology is the science of studying venoms and toxins that are produced by a wide variety of organisms. In the past, the structure, function and mechanisms of most abundant and/or most toxic components were characterized to understand and to develop strategies to neutralize their toxicity. With recent technical advances, we are able to evaluate and determine the toxin profiles using transcriptomes of venom glands and proteomes of tiny amounts of venom. Enormous amounts of data from these studies have opened tremendous opportunities in many directions of basic and applied research. The lower costs for profiling venoms will further fuel the expansion of toxin database, which in turn will provide greater exciting and bright opportunities in toxin research.

Venom Systems as Models for Studying the Origin and Regulation of Evolutionary Novelties

A central goal in biology is to determine the ways in which evolution repeats itself. One of the most remarkable examples in nature of convergent evolutionary novelty is animal venom. Across diverse animal phyla, various specialized organs and anatomical structures have evolved from disparate developmental tissues to perform the same function, that is, produce and deliver a cocktail of potent molecules to subdue prey or predators. Venomous organisms therefore offer unique opportunities to investigate the evolutionary processes of convergence of key adaptive traits, and the molecular mechanisms underlying the emergence of novel genes, cells, and tissues. Indeed, some venomous species have already proven to be highly amenable as models for developmental studies, and recent work with venom gland organoids provides manipulatable systems for directly testing important evolutionary questions. Here, we provide a synthesis of the current knowledge that could serve as a starting point for the establishment of venom systems as new models for evolutionary and molecular biology. In particular, we highlight the potential of various venomous species for the study of cell differentiation and cell identity, and the regulatory dynamics of rapidly evolving, highly expressed, tissue-specific, gene paralogs. We hope that this review will encourage researchers to look beyond traditional study organisms and consider venom systems as useful tools to explore evolutionary novelties.

The Primary Duct of Bothrops jararaca Glandular Apparatus Secretes Toxins

Despite numerous studies concerning morphology and venom production and secretion in the main venom gland (and some data on the accessory gland) of the venom glandular apparatus of Viperidae snakes, the primary duct has been overlooked. We characterized the primary duct of the Bothrops jararaca snake by morphological analysis, immunohistochemistry and proteomics. The duct has a pseudostratified epithelium with secretory columnar cells with vesicles of various electrondensities, as well as mitochondria-rich, dark, basal, and horizontal cells. Morphological analysis, at different periods after venom extraction, showed that the primary duct has a long cycle of synthesis and secretion, as do the main venom and accessory glands; however, the duct has a mixed mode venom storage, both in the lumen and in secretory vesicles. Mouse anti-B. jararaca venom serum strongly stained the primary duct's epithelium. Subsequent proteomic analysis revealed the synthesis of venom toxins-mainly C-type lectin/C-type lectin-like proteins. We propose that the primary duct's toxin synthesis products complement the final venom bolus. Finally, we hypothesize that the primary duct and the accessory gland (components of the venom glandular apparatus) are part of the evolutionary path from a salivary gland towards the main venom gland.

Venom regeneration in the centipede Scolopendra polymorpha: evidence for asynchronous venom component synthesis

Venom regeneration comprises a vital process in animals that rely on venom for prey capture and defense. Venom regeneration in scolopendromorph centipedes likely influences their ability to subdue prey and defend themselves, and may influence the quantity and quality of venom extracted by researchers investigating the venom's biochemistry. We investigated venom volume and total protein regeneration during the 14-day period subsequent to venom extraction in the North American centipede Scolopendra polymorpha. We further tested the hypothesis that venom protein components, separated by reversed-phase fast protein liquid chromatography (RP-FPLC), undergo asynchronous (non-parallel) synthesis. During the first 48 h, volume and protein mass increased linearly. Protein regeneration lagged behind volume regeneration, with 65–86% of venom volume and 29–47% of protein mass regenerated during the first 2 days. No additional regeneration occurred over the subsequent 12 days, and neither volume nor protein mass reached initial levels 7 months later (93% and 76%, respectively). Centipede body length was negatively associated with rate of venom regeneration. Analysis of chromatograms of individual venom samples revealed that 5 of 10 chromatographic regions and 12 of 28 peaks demonstrated changes in percent of total peak area (i.e., percent of total protein) among milking intervals, indicating that venom proteins are regenerated asynchronously. Moreover, specimens from Arizona and California differed in relative amounts of some venom components. The considerable regeneration of venom occurring within the first 48 h, despite the reduced protein content, suggests that predatory and defensive capacities are minimally constrained by the timing of venom replacement.

Venom-related transcripts from Bothrops jararaca tissues provide novel molecular insights into the production and evolution of snake venom

Attempts to reconstruct the evolutionary history of snake toxins in the context of their co-option to the venom gland rarely account for nonvenom snake genes that are paralogous to toxins, and which therefore represent important connectors to ancestral genes. In order to reevaluate this process, we conducted a comparative transcriptomic survey on body tissues from a venomous snake. A nonredundant set of 33,000 unigenes (assembled transcripts of reference genes) was independently assembled from six organs of the medically important viperid snake Bothrops jararaca, providing a reference list of 82 full-length toxins from the venom gland and specific products from other tissues, such as pancreatic digestive enzymes. Unigenes were then screened for nontoxin transcripts paralogous to toxins revealing 1) low level coexpression of approximately 20% of toxin genes (e.g., bradykinin-potentiating peptide, C-type lectin, snake venom metalloproteinase, snake venom nerve growth factor) in body tissues, 2) the identity of the closest paralogs to toxin genes in eight classes of toxins, 3) the location and level of paralog expression, indicating that, in general, co-expression occurs in a higher number of tissues and at lower levels than observed for toxin genes, and 4) strong evidence of a toxin gene reverting back to selective expression in a body tissue. In addition, our differential gene expression analyses identify specific cellular processes that make the venom gland a highly specialized secretory tissue. Our results demonstrate that the evolution and production of venom in snakes is a complex process that can only be understood in the context of comparative data from other snake tissues, including the identification of genes paralogous to venom toxins.