RESUMEN
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.
Asunto(s)
Venenos de Crotálidos , Crotalus , Evolución Molecular , Animales , Crotalus/genética , Venenos de Crotálidos/genética , Redes Reguladoras de Genes , Regulación de la Expresión GénicaRESUMEN
Disintegrins, a family of snake venom protein, which are capable of modulating the activity of integrins that play a fundamental role in the regulation of many physiological and pathological processes. The main purpose of this study is to obtain the recombinant disintegrin (r-DI) and evaluate its biological activity. In this study, we explored a high-level expression prokaryotic system and purification strategy for r-DI. Then, r-DI was treated to assay effects on cell growth, migration, and invasion. The affinity for the interactions of r-DI with integrin was determined using Surface plasmon resonance (SPR) analyses. The r-DI can be expressed in Escherichia coli and purified by one-step chromatography. The r-DI can inhibit B16F10 cells proliferation, migration, and invasion. Also, we found that r-DI could interact with the integrin αIIbß3 (GPIIb/IIIa). The r-DI can be expressed, purified, characterized through functional assays, and can also maintain strong biological activities. Thus, this study showed potential therapeutic effects of r-DI for further functional and structural studies.
Asunto(s)
Desintegrinas , Escherichia coli , Proteínas Recombinantes , Escherichia coli/genética , Escherichia coli/metabolismo , Animales , Desintegrinas/química , Desintegrinas/genética , Desintegrinas/aislamiento & purificación , Desintegrinas/farmacología , Proteínas Recombinantes/genética , Proteínas Recombinantes/química , Proteínas Recombinantes/aislamiento & purificación , Proteínas Recombinantes/biosíntesis , Proteínas Recombinantes/metabolismo , Ratones , Viperidae/genética , Complejo GPIIb-IIIa de Glicoproteína Plaquetaria/genética , Complejo GPIIb-IIIa de Glicoproteína Plaquetaria/química , Complejo GPIIb-IIIa de Glicoproteína Plaquetaria/metabolismo , Línea Celular Tumoral , Expresión Génica , Movimiento Celular/efectos de los fármacos , Proliferación Celular/efectos de los fármacos , Venenos de Crotálidos/química , Venenos de Crotálidos/genética , Crotalinae , Serpientes VenenosasRESUMEN
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.
Asunto(s)
Venenos de Crotálidos , Serpientes Venenosas , Animales , Venenos de Crotálidos/genética , Venenos de Crotálidos/metabolismo , Epigenómica , Crotalus/genética , Crotalus/metabolismoRESUMEN
RNA purification and cDNA synthesis represents the starting point for molecular analyses of snake venom proteins-enzymes. Usually, the sacrifice of snakes is necessary for venom gland extraction to identify protein-coding transcripts; however, the venom can be used as a source of transcripts. Although there are methods for obtaining RNA from venom, no comparative analysis has been conducted in the Bothrops genus. In the present study, we compared four commercial methods for RNA purification and cDNA synthesis from venom (liquid, lyophilized, or long-term storage) of four clinically relevant species of Peruvian Bothrops. Our results show that the TRIzol method presents the highest yield of RNA purified from venom (59 ± 11 ng/100 µL or 10 mg). The SuperScript First-Strand Synthesis System kit produced high amounts of cDNA (3.2 ± 1.2 ng cDNA/ng RNA), and the highest value was from combination with the Dynabeads mRNA DIRECT kit (4.8 ± 2.0 ng cDNA/ng RNA). The utility of cDNA was demonstrated with the amplification of six relevant toxins: thrombin-like enzymes, P-I and P-III metalloproteinases, acid and basic phospholipases A2, and disintegrins. To our knowledge, this is the first comparative study of RNA purification and cDNA synthesis methodologies from Bothrops genus venom.
Asunto(s)
Bothrops , Venenos de Crotálidos , Animales , ADN Complementario/genética , Bothrops/genética , Perú , Relevancia Clínica , Venenos de Crotálidos/genética , Proteínas , ARNRESUMEN
Studying the consequences of hybridization between closely related species with divergent traits can reveal patterns of evolution that shape and maintain extreme trophic adaptations. Snake venoms are an excellent model system for examining the evolutionary and ecological patterns that underlie highly selected polymorphic traits. Here we investigate hybrid venom phenotypes that result from natural introgression between two rattlesnake species that express highly divergent venom phenotypes: Crotalus o. concolor and C. v. viridis. Though not yet documented, interbreeding between these species may lead to novel venom phenotypes with unique activities that break the typical trends of venom composition in rattlesnakes. The characteristics of these unusual phenotypes could unveil the roles of introgression in maintaining patterns of venom composition and variation, including the near ubiquitous dichotomy between neurotoxic or degradative venoms observed across rattlesnakes. We use RADseq data to infer patterns of gene flow and hybrid ancestry between these diverged lineages and link these genetic data with analyses of venom composition, biological activity, and whole animal model toxicity tests to understand the impacts of introgression on venom composition. We find that introgressed populations express admixed venom phenotypes that do not sacrifice biological activity (lethal toxicity) or overall abundance of dominant toxins compared to parental venoms. These hybridized venoms therefore do not represent a trade-off in functionality between the typical phenotypic extremes but instead represent a unique combination of characters whose expression appears limited to the hybrid zone.
Asunto(s)
Venenos de Crotálidos , Toxinas Biológicas , Animales , Crotalus/genética , Crotalus/metabolismo , Toxinas Biológicas/metabolismo , Venenos de Serpiente , Fenotipo , Venenos de Crotálidos/genética , Venenos de Crotálidos/toxicidadRESUMEN
BACKGROUND: Snake venoms are trophic adaptations that represent an ideal model to examine the evolutionary factors that shape polymorphic traits under strong natural selection. Venom compositional variation is substantial within and among venomous snake species. However, the forces shaping this phenotypic complexity, as well as the potential integrated roles of biotic and abiotic factors, have received little attention. Here, we investigate geographic variation in venom composition in a wide-ranging rattlesnake (Crotalus viridis viridis) and contextualize this variation by investigating dietary, phylogenetic, and environmental variables that covary with venom. RESULTS: Using shotgun proteomics, venom biochemical profiling, and lethality assays, we identify 2 distinct divergent phenotypes that characterize major axes of venom variation in this species: a myotoxin-rich phenotype and a snake venom metalloprotease (SVMP)-rich phenotype. We find that dietary availability and temperature-related abiotic factors are correlated with geographic trends in venom composition. CONCLUSIONS: Our findings highlight the potential for snake venoms to vary extensively within species, for this variation to be driven by biotic and abiotic factors, and for the importance of integrating biotic and abiotic variation for understanding complex trait evolution. Links between venom variation and variation in biotic and abiotic factors indicate that venom variation likely results from substantial geographic variation in selection regimes that determine the efficacy of venom phenotypes across populations and snake species. Our results highlight the cascading influence of abiotic factors on biotic factors that ultimately shape venom phenotype, providing evidence for a central role of local selection as a key driver of venom variation.
Asunto(s)
Venenos de Crotálidos , Crotalus , Animales , Crotalus/genética , Filogenia , Venenos de Serpiente/genética , Venenos de Serpiente/química , Fenotipo , Venenos de Crotálidos/genética , Venenos de Crotálidos/químicaRESUMEN
Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.
RESUMEN: El veneno de las serpientes puede variar entre y dentro de las especies. Mientras algunos grupos de viperidos del Nuevo Mundocomo las cascabeleshan sido bien estudiadas, muy poco se sabe acerca del veneno de las nauyacas de frío (Cerrophidion) que se encuentran en las zonas altas de Mesoamérica. Comparadas con las extensamente estudiadas cascabeles, que estan ampliamente distribuidas, las poblaciones de Cerrophidion, aisladas en montañas, pueden poseer trayectorias evolutivas y diferenciación en su veneno unicos. En el presente trabajo, describimos el transcriptoma de las glándulas de veneno de poblaciones de C. petlalcalensis, C. tzotzilorum, y C. godmani de México, y un individuo de C. sasai de Costa Rica. Exploramos la variación en la expresión de toxinas en Cerrophidion y la evolución en las secuencias geneticas en C. godmani específicamente. El transcriptoma de la glándula de veneno de Cerrophidion esta compuesto principalmente de Metaloproteinasas de Veneno de Serpiente, Fosfolipasas A[Formula: see text] (PLA[Formula: see text]s), y Serin Proteasas de Veneno de Serpiente. Cerrophidion petlalcalensis presenta poca variación intraespecífica; sin embargo, los transcriptomas de la glandula de veneno de C. godmani y C. tzotzilorum difieren significativamente entre poblaciones geográficamente aisladas. Curiosamente, la variación intraespecífica estuvo atribuida principalmente a la expresión de las toxinas ya que no encontramos señales de selección en las toxinas de C. godmani. Adicionalmente, encontramos miotoxinas similares a PLA[Formula: see text] en todas las especies excepto C. petlalcalensis, y PLA[Formula: see text]s similares a crotoxina en la población sureña de C. godmani. Nuestros resultados demuestran la presencia de variacion intraespecífica presente en el veneno de C. godmani y C. tzotzilorum. Las toxinas de Cerrophidion godmani muestran poca evidencia de selección direccional, y la variación en la secuencias de las toxinas es consistente con evolucion bajo un modelo de equilibrio de mutación-deriva. Algunos individuos de C. godmani de la población del sur potencialmente tienen un veneno neurotóxico dada la presencia de PLA[Formula: see text]s similares a la crotoxina, sin embargo, se necesita más evidencia para corroborar esta hipótesis.
Asunto(s)
Venenos de Crotálidos , Crotalinae , Crotoxina , Viperidae , Humanos , Animales , Crotalinae/genética , Crotalinae/metabolismo , Viperidae/metabolismo , Crotoxina/metabolismo , Venenos de Crotálidos/genética , Venenos de Crotálidos/metabolismo , Venenos de Crotálidos/toxicidad , Venenos de Serpiente/metabolismo , Poliésteres/metabolismoRESUMEN
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.
Asunto(s)
Adaptación Fisiológica , Venenos de Crotálidos , Sciuridae , Animales , Adaptación Fisiológica/genética , Venenos de Crotálidos/genética , Perfilación de la Expresión Génica , Sciuridae/genética , Resistencia a Medicamentos/genéticaRESUMEN
The complexity of snake venom composition reflects adaptation to the diversity of prey and may be driven at times by a coevolutionary arms race between snakes and venom-resistant prey. However, many snakes are also resistant to their own venom due to serum-borne inhibitors of venom toxins, which raises the question of how snake autoinhibitors maintain their efficacy as venom proteins evolve. To investigate this potential three-way arms race among venom, prey, and autoinhibitors, we have identified and traced the evolutionary origin of serum inhibitors of snake venom metalloproteinases (SVMPs) in the Western Diamondback rattlesnake Crotalus atrox which possesses the largest known battery of SVMP genes among crotalids examined. We found that C. atrox expresses five members of a Fetuin A-related metalloproteinase inhibitor family but that one family member, FETUA-3, is the major SVMP inhibitor that binds to approximately 20 different C. atrox SVMPs and inhibits activities of all three SVMP classes. We show that the fetua-3 gene arose deep within crotalid evolution before the origin of New World species but, surprisingly, fetua-3 belongs to a different paralog group than previously identified SVMP inhibitors in Asian and South American crotalids. Conversely, the C. atrox FETUA-2 ortholog of previously characterized crotalid SVMP inhibitors shows limited activity against C. atrox SVMPs. These results reveal that there has been a functional evolutionary shift in the major SVMP inhibitor in the C. atrox lineage as the SVMP family expanded and diversified in the Crotalus lineage. This broad-spectrum inhibitor may be of potential therapeutic interest.
Asunto(s)
Venenos de Crotálidos , Toxinas Biológicas , Animales , Crotalus/genética , Venenos de Crotálidos/genética , Venenos de Crotálidos/metabolismo , Metaloproteasas/genética , Metaloproteasas/metabolismo , Venenos de Serpiente/metabolismo , Toxinas Biológicas/metabolismoRESUMEN
The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance.
Asunto(s)
Venenos de Crotálidos , Animales , Venenos de Crotálidos/genética , Crotalus/genética , Genoma , Recombinación Genética , Venenos de Serpiente/genéticaRESUMEN
Crotamine, myotoxin a and homologs are short peptides that often comprise major fractions of rattlesnake venoms and have been extensively studied for their bioactive properties. These toxins are thought to be important for rapidly immobilizing mammalian prey and are implicated in serious, and sometimes fatal, responses to envenomation in humans. While high quality reference genomes for multiple venomous snakes are available, the loci that encode myotoxins have not been successfully assembled in any existing genome assembly. Here, we integrate new and existing genomic and transcriptomic data from the Prairie Rattlesnake (Crotalus viridis viridis) to reconstruct, characterize, and infer the chromosomal locations of myotoxin-encoding loci. We integrate long-read transcriptomics (Pacific Bioscience's Iso-Seq) and short-read RNA-seq to infer gene sequence diversity and characterize patterns of myotoxin and paralogous ß-defensin expression across multiple tissues. We also identify two long non-coding RNA sequences which both encode functional myotoxins, demonstrating a newly discovered source of venom coding sequence diversity. We also integrate long-range mate-pair chromatin contact data and linked-read sequencing to infer the structure and chromosomal locations of the three myotoxin-like loci. Further, we conclude that the venom-associated myotoxin is located on chromosome 1 and is adjacent to non-venom paralogs. Consistent with this locus contributing to venom composition, we find evidence that the promoter of this gene is selectively open in venom gland tissue and contains transcription factor binding sites implicated in broad trans-regulatory pathways that regulate snake venoms. This study provides the best genomic reconstruction of myotoxin loci to date and raises questions about the physiological roles and interplay between myotoxin and related genes, as well as the genomic origins of snake venom variation.
Asunto(s)
Venenos de Crotálidos , Crotalus/fisiología , Neurotoxinas , Animales , Secuencia de Bases , Venenos de Crotálidos/química , Venenos de Crotálidos/genética , Crotalus/genética , Variaciones en el Número de Copia de ADN , Genómica , Humanos , Mamíferos , Venenos de Serpiente/química , Venenos de Serpiente/genética , TranscriptomaRESUMEN
Opossums in the tribe Didelphini are resistant to pit viper venoms and are hypothesized to be coevolving with venomous snakes. Specifically, a protein involved in blood clotting (von Willebrand factor [vWF] which is targeted by snake venom C-type lectins [CTLs]) has been found to undergo rapid adaptive evolution in Didelphini. Several unique amino acid changes in vWF could explain their resistance; however, experimental evidence that these changes disrupt binding to venom CTLs was lacking. Furthermore, without explicit testing of ancestral phenotypes to reveal the mode of evolution, the assertion that this system represents an example of coevolution rather than noncoevolutionary adaptation remains unsupported. Using expressed vWF proteins and purified venom CTLs, we quantified binding affinity for vWF proteins from all resistant taxa, their venom-sensitive relatives, and their ancestors. We show that CTL-resistant vWF is present in opossums outside clade Didelphini and likely across a wider swath of opossums (family Didelphidae) than previously thought. Ancestral reconstruction and in vitro testing of vWF phenotypes in a clade of rapidly evolving opossums reveal a pattern consistent with trench warfare coevolution between opossums and their venomous snake prey.
Asunto(s)
Venenos de Crotálidos , Crotalinae , Animales , Venenos de Crotálidos/genética , Zarigüeyas/metabolismo , Venenos de Serpiente/metabolismo , Serpientes/metabolismo , Factor de von Willebrand/genética , Factor de von Willebrand/metabolismoRESUMEN
Interspecific differences in snake venom compositions can result from distinct regulatory mechanisms acting in each species. However, comparative analyses focusing on identifying regulatory elements and patterns that led to distinct venom composition are still scarce. Among venomous snakes, Bothrops cotiara and Bothrops fonsecai represent ideal models to complement our understanding of the regulatory mechanisms of venom production. These recently diverged species share a similar specialized diet, habitat, and natural history, but each presents a distinct venom phenotype. Here, we integrated data from the venom gland transcriptome and miRNome and the venom proteome of B. fonsecai and B. cotiara to better understand the regulatory mechanisms that may be acting to produce differing venom compositions. We detected not only the presence of similar toxin isoforms in both species but also distinct expression profiles of phospholipases A2 (PLA2) and some snake venom metalloproteinases (SVMPs) and snake venom serine proteinases (SVSPs) isoforms. We found evidence of modular expression regulation of several toxin isoforms implicated in venom divergence and observed correlated expression of several transcription factors. We did not find strong evidence for miRNAs shaping interspecific divergence of the venom phenotypes, but we identified a subset of toxin isoforms whose final expression may be fine-tuned by specific miRNAs. Sequence analysis on orthologous toxins showed a high rate of substitutions between PLA2s, which indicates that these toxins may be under strong positive selection or represent paralogous toxins in these species. Our results support other recent studies in suggesting that gene regulation is a principal mode of venom evolution across recent timescales, especially among species with conserved ecotypes.
Asunto(s)
Bothrops , Venenos de Crotálidos , MicroARNs , Toxinas Biológicas , Animales , Bothrops/genética , Bothrops/metabolismo , Brasil , Venenos de Crotálidos/genética , Venenos de Crotálidos/metabolismo , MicroARNs/metabolismo , Fosfolipasas A2/genética , Fosfolipasas A2/metabolismo , Venenos de Serpiente/metabolismo , Toxinas Biológicas/metabolismoRESUMEN
Proteins isolated from natural sources can be composed of a mixture of isoforms with similar physicochemical properties that coexist in the final steps of purification. Yet, even where unverified, the assumed sequence is enforced throughout the structural studies. Herein, we propose a novel perspective to address the usually neglected sequence heterogeneity of natural products by integrating biophysical, genetic and structural data in our program SEQUENCE SLIDER. The aim is to assess the evidence supporting chemical composition in structure determination. Locally, we interrogate the experimental map to establish which side chains are supported by the structural data, and the genetic information relating sequence conservation is integrated into this statistic. Hence, we build a constrained peptide database, containing most probable sequences to interpret mass spectrometry data (MS). In parallel, we perform MS de novo sequencing with genomic-based algorithms to detect point mutations. We calibrated SLIDER with Gallus gallus lysozyme, whose sequence is unequivocally established and numerous natural isoforms are reported. We used SLIDER to characterize a metalloproteinase and a phospholipase A2-like protein from the venom of Bothrops moojeni and a crotoxin from Crotalus durissus collilineatus. This integrated approach offers a more realistic structural descriptor to characterize macromolecules isolated from natural sources.
Asunto(s)
Mezclas Complejas/química , Isoformas de Proteínas/análisis , Programas Informáticos , Animales , Venenos de Crotálidos/química , Venenos de Crotálidos/genética , Crotalus/genética , Crotoxina/química , Crotoxina/genética , Fosfolipasas A2/químicaRESUMEN
Ecologically divergent selection can lead to the evolution of reproductive isolation through the process of ecological speciation, but the balance of responsible evolutionary forces is often obscured by an inadequate assessment of demographic history and the genetics of traits under selection. Snake venoms have emerged as a system for studying the genetic basis of adaptation because of their genetic tractability and contributions to fitness, and speciation in venomous snakes can be associated with ecological diversification such as dietary shifts and corresponding venom changes. Here, we explored the neurotoxic (type A)-hemotoxic (type B) venom dichotomy and the potential for ecological speciation among Timber Rattlesnake (Crotalus horridus) populations. Previous work identified the genetic basis of this phenotypic difference, enabling us to characterize the roles geography, history, ecology, selection, and chance play in determining when and why new species emerge or are absorbed. We identified significant genetic, proteomic, morphological, and ecological/environmental differences at smaller spatial scales, suggestive of incipient ecological speciation between type A and type B C. horridus. Range-wide analyses, however, rejected the reciprocal monophyly of venom type, indicative of varying intensities of introgression and a lack of reproductive isolation across the range. Given that we have now established the phenotypic distributions and ecological niche models of type A and B populations, genome-wide data are needed and capable of determining whether type A and type B C. horridus represent distinct, reproductively isolated lineages due to incipient ecological speciation or differentiated populations within a single species.
Asunto(s)
Venenos de Crotálidos/genética , Crotalus/genética , Introgresión Genética , Especiación Genética , AnimalesRESUMEN
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.
Asunto(s)
Bothrops/genética , Venenos de Crotálidos/genética , Evolución Molecular , Genoma/genética , Venenos de Serpiente/genética , Secuencia de Aminoácidos , Animales , Secuencia de Bases , Bothrops/clasificación , Venenos de Crotálidos/clasificación , Femenino , Perfilación de la Expresión Génica/métodos , Filogenia , Proteoma/metabolismo , Proteómica/métodos , RNA-Seq/métodos , Análisis de Secuencia de ADN/métodos , Venenos de Serpiente/clasificaciónRESUMEN
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.
Asunto(s)
Venenos de Crotálidos/genética , Crotalus/genética , Genoma/genética , Anotación de Secuencia Molecular , Animales , Regulación de la Expresión Génica/genética , Genotipo , Transcriptoma/genética , Secuenciación Completa del GenomaRESUMEN
Ontogenetic changes in venom composition have been described in Bothrops snakes, but only a few studies have attempted to identify the targeted paralogues or the molecular mechanisms involved in modifications of gene expression during ontogeny. In this study, we decoded B. jararacussu venom gland transcripts from six specimens of varying sizes and analyzed the variability in the composition of independent venom proteomes from 19 individuals. We identified 125 distinct putative toxin transcripts, and of these, 73 were detected in venom proteomes and only 10 were involved in the ontogenetic changes. Ontogenetic variability was linearly related to snake size and did not correspond to the maturation of the reproductive stage. Changes in the transcriptome were highly predictive of changes in the venom proteome. The basic myotoxic phospholipases A2 (PLA2s) were the most abundant components in larger snakes, while in venoms from smaller snakes, PIII-class SVMPs were the major components. The snake venom metalloproteinases (SVMPs) identified corresponded to novel sequences and conferred higher pro-coagulant and hemorrhagic functions to the venom of small snakes. The mechanisms modulating venom variability are predominantly related to transcriptional events and may consist of an advantage of higher hematotoxicity and more efficient predatory function in the venom from small snakes.
Asunto(s)
Tamaño Corporal/genética , Bothrops/genética , Venenos de Crotálidos/genética , Proteómica/métodos , Transcriptoma/genética , Animales , Venenos de Crotálidos/análisis , Venenos de Crotálidos/química , Femenino , Ontología de Genes , Masculino , Análisis de Secuencia de ADN/métodosRESUMEN
Ontogenetic shifts in venom occur in many snakes but establishing their nature as gradual or discrete processes required additional study. We profiled shifts in venom expression from the neonate to adult sizes of two rattlesnake species, the eastern diamondback and the timber rattlesnake. We used serial sampling and venom chromatographic profiling to test if ontogenetic change occurs gradually or discretely. We found evidence for gradual shifts in overall venom composition in six of eight snakes, which sometimes spanned more than two years. Most chromatographic peaks shift gradually, but one quarter shift in a discrete fashion. Analysis of published diet data showed gradual shifts in overall diet composition across the range of body sizes attained by our eight study animals, while the shifts in abundance of different prey classes varied in form from gradual to discrete. Testosterone concentrations were correlated with the change in venom protein composition, but the relationship is not strong enough to suggest causation. Venom research employing simple juvenile versus adult size thresholds may be failing to account for continuous variation in venom composition lifespan. Our results imply that venom shifts represent adaptive matches to dietary shifts and highlight venom for studies of alternative gene regulatory mechanisms.
Asunto(s)
Venenos de Crotálidos/metabolismo , Crotalus/metabolismo , Ecosistema , Proteínas de Reptiles/metabolismo , Testosterona/metabolismo , Factores de Edad , Animales , Tamaño Corporal , Venenos de Crotálidos/genética , Crotalus/genética , Crotalus/crecimiento & desarrollo , Dieta , Regulación del Desarrollo de la Expresión Génica , Proteínas de Reptiles/genéticaRESUMEN
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.