1
|
Entiauspe-Neto OM, Nachtigall PG, Borges-Martins M, Junqueira-de-Azevedo ILM, Grazziotin FG. Highly conserved and extremely variable: The paradoxical pattern of toxin expression revealed by comparative venom-gland transcriptomics of Phalotris (Serpentes: Dipsadidae). Toxicon 2024; 244:107740. [PMID: 38705487 DOI: 10.1016/j.toxicon.2024.107740] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2024] [Revised: 04/11/2024] [Accepted: 04/27/2024] [Indexed: 05/07/2024]
Abstract
Although non-front fanged snakes account for almost two-thirds of snake diversity, most studies on venom composition and evolution focus exclusively on front-fanged species, which comprise most of the clinically relevant accidents. Comprehensive reports on venom composition of non-front fanged snakes are still scarce for several groups. In this study, we address such shortage of knowledge by providing new insights about the venom composition among species of Phalotris, a poorly studied Neotropical dipsadid genus. Phalotris are known for their specialized venom delivery system and toxic venoms, which can cause life-threatening accidents in humans. We evaluate the venom-gland transcriptome of Phalotris, comparing the following three South American species: P. reticulatus for the Araucaria Pine forests, P. lemniscatus for the Pampa grasslands, and P. mertensi for the Brazilian Cerrado. Our results indicate similar venom profiles, in which they share a high expression level of Kunitz-type inhibitors (KUNZ). On the other hand, comparative analyses revealed substantial differences in the expression levels of C-type lectins (CTL) and snake venom metalloproteinases (SVMP). The diverse set of SVMP and CTL isoforms shows signals of positive selection, and we also identified truncated forms of type III SVMPs, which resemble type II and type I SVMPs of viperids. Additionally, we identified a CNP precursor hosting a proline-rich region containing a BPP motif resembling those commonly detected in viperid venoms with hypotensive activity. Altogether, our results suggest an evolutionary history favoring high expression levels of few KUNZ isoforms in Phalotris venoms, contrasting with a highly diverse set of SVMP and CTL isoforms. Such diversity can be comparable with the venom variability observed in some viperids. Our findings highlight the extreme phenotypic diversity of non-front fanged snakes and the importance to allocate greater effort to study neglected groups of Colubroidea.
Collapse
Affiliation(s)
- Omar M Entiauspe-Neto
- Laboratório de Coleções Zoológicas, Instituto Butantan, 05503-900, Av. Vital Brazil, 1500, Butantã, São Paulo, SP, Brazil; Programa de Pós-graduação em Biologia Animal, Departamento de Zoologia, Instituto de Biociências, Universidade Federal do Rio Grande do Sul, Av. Bento Gonçalves, CEP 91501-970, Porto Alegre, RS, Brazil.
| | - Pedro G Nachtigall
- Laboratório de Toxinologia Aplicada, Instituto Butantan, 05503-900, Av. Vital Brazil, 1500, Butantã, São Paulo, SP, Brazil
| | - Márcio Borges-Martins
- Programa de Pós-graduação em Biologia Animal, Departamento de Zoologia, Instituto de Biociências, Universidade Federal do Rio Grande do Sul, Av. Bento Gonçalves, CEP 91501-970, Porto Alegre, RS, Brazil
| | | | - Felipe G Grazziotin
- Laboratório de Coleções Zoológicas, Instituto Butantan, 05503-900, Av. Vital Brazil, 1500, Butantã, São Paulo, SP, Brazil
| |
Collapse
|
2
|
Hogan MP, Holding ML, Nystrom GS, Colston TJ, Bartlett DA, Mason AJ, Ellsworth SA, Rautsaw RM, Lawrence KC, Strickland JL, He B, Fraser P, Margres MJ, Gilbert DM, Gibbs HL, Parkinson CL, Rokyta DR. The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom. Proc Natl Acad Sci U S A 2024; 121:e2313440121. [PMID: 38578985 PMCID: PMC11032440 DOI: 10.1073/pnas.2313440121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Accepted: 03/13/2024] [Indexed: 04/07/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Michael P. Hogan
- Department of Biological Science, Florida State University, Tallahassee, FL32306
| | - Matthew L. Holding
- Department of Biological Science, Florida State University, Tallahassee, FL32306
- Life Sciences Institute, University of Michigan, Ann Arbor, MI48109
| | - Gunnar S. Nystrom
- Department of Biological Science, Florida State University, Tallahassee, FL32306
| | - Timothy J. Colston
- Department of Biological Science, Florida State University, Tallahassee, FL32306
- Department of Biology, University of Puerto Rico at Mayagüez, Mayagüez, PR00681
| | - Daniel A. Bartlett
- Department of Biological Science, Florida State University, Tallahassee, FL32306
| | - Andrew J. Mason
- Department of Biological Sciences, Clemson University, Clemson, SC29634
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University, Columbus, OH43210
| | - Schyler A. Ellsworth
- Department of Biological Science, Florida State University, Tallahassee, FL32306
| | - Rhett M. Rautsaw
- Department of Biological Sciences, Clemson University, Clemson, SC29634
- Department of Integrative Biology, University of South Florida, Tampa, FL33620
- School of Biological Sciences, Washington State University, Pullman, WA99164
| | - Kylie C. Lawrence
- Department of Biological Science, Florida State University, Tallahassee, FL32306
| | - Jason L. Strickland
- Department of Biological Sciences, Clemson University, Clemson, SC29634
- Department of Biology, University of South Alabama, Mobile, AL36688
| | - Bing He
- Department of Biological Science, Florida State University, Tallahassee, FL32306
| | - Peter Fraser
- Department of Biological Science, Florida State University, Tallahassee, FL32306
| | - Mark J. Margres
- Department of Integrative Biology, University of South Florida, Tampa, FL33620
| | - David M. Gilbert
- Laboratory of Chromosome Replication and Epigenome Regulation, San Diego Biomedical Research Institute, San Diego, CA92121
| | - H. Lisle Gibbs
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University, Columbus, OH43210
| | - Christopher L. Parkinson
- Department of Biological Sciences, Clemson University, Clemson, SC29634
- Department of Forestry and Environmental Conservation, Clemson University, Clemson, SC29634
| | - Darin R. Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL32306
| |
Collapse
|
3
|
Freitas-de-Sousa LA, Colombini M, Souza VC, Silva JPC, Mota-da-Silva A, Almeida MRN, Machado RA, Fonseca WL, Sartim MA, Sachett J, Serrano SMT, Junqueira-de-Azevedo ILM, Grazziotin FG, Monteiro WM, Bernarde PS, Moura-da-Silva AM. Venom Composition of Neglected Bothropoid Snakes from the Amazon Rainforest: Ecological and Toxinological Implications. Toxins (Basel) 2024; 16:83. [PMID: 38393161 PMCID: PMC10891915 DOI: 10.3390/toxins16020083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 01/18/2024] [Accepted: 01/30/2024] [Indexed: 02/25/2024] Open
Abstract
Snake venoms have evolved in several families of Caenophidae, and their toxins have been assumed to be biochemical weapons with a role as a trophic adaptation. However, it remains unclear how venom contributes to the success of venomous species for adaptation to different environments. Here we compared the venoms from Bothrocophias hyoprora, Bothrops taeniatus, Bothrops bilineatus smaragdinus, Bothrops brazili, and Bothrops atrox collected in the Amazon Rainforest, aiming to understand the ecological and toxinological consequences of venom composition. Transcriptomic and proteomic analyses indicated that the venoms presented the same toxin groups characteristic from bothropoids, but with distinct isoforms with variable qualitative and quantitative abundances, contributing to distinct enzymatic and toxic effects. Despite the particularities of each venom, commercial Bothrops antivenom recognized the venom components and neutralized the lethality of all species. No clear features could be observed between venoms from arboreal and terrestrial habitats, nor in the dispersion of the species throughout the Amazon habitats, supporting the notion that venom composition may not shape the ecological or toxinological characteristics of these snake species and that other factors influence their foraging or dispersal in different ecological niches.
Collapse
Affiliation(s)
| | - Mônica Colombini
- Laboratório de Imunopatologia, Instituto Butantan, São Paulo 05503-900, SP, Brazil; (L.A.F.-d.-S.); (M.C.)
| | - Vinicius C. Souza
- Laboratório de Toxinologia Aplicada, Instituto Butantan, São Paulo 05503-900, SP, Brazil; (V.C.S.); (J.P.C.S.); (S.M.T.S.); (I.L.M.J.-d.-A.)
| | - Joanderson P. C. Silva
- Laboratório de Toxinologia Aplicada, Instituto Butantan, São Paulo 05503-900, SP, Brazil; (V.C.S.); (J.P.C.S.); (S.M.T.S.); (I.L.M.J.-d.-A.)
| | - Ageane Mota-da-Silva
- Instituto Federal do Acre, Campus de Cruzeiro do Sul, Cruzeiro do Sul 69980-000, AC, Brazil;
| | - Marllus R. N. Almeida
- Laboratório de Herpetologia, Universidade Federal do Acre, Campus Floresta, Cruzeiro do Sul 69895-000, AC, Brazil; (M.R.N.A.); (R.A.M.); (W.L.F.); (P.S.B.)
| | - Reginaldo A. Machado
- Laboratório de Herpetologia, Universidade Federal do Acre, Campus Floresta, Cruzeiro do Sul 69895-000, AC, Brazil; (M.R.N.A.); (R.A.M.); (W.L.F.); (P.S.B.)
| | - Wirven L. Fonseca
- Laboratório de Herpetologia, Universidade Federal do Acre, Campus Floresta, Cruzeiro do Sul 69895-000, AC, Brazil; (M.R.N.A.); (R.A.M.); (W.L.F.); (P.S.B.)
| | - Marco A. Sartim
- Instituto de Pesquisa Clínica Carlos Borborema, Fundação de Medicina Tropical Dr. Heitor Vieira Dourado, Manaus 69040-000, AM, Brazil; (M.A.S.); (J.S.); (W.M.M.)
| | - Jacqueline Sachett
- Instituto de Pesquisa Clínica Carlos Borborema, Fundação de Medicina Tropical Dr. Heitor Vieira Dourado, Manaus 69040-000, AM, Brazil; (M.A.S.); (J.S.); (W.M.M.)
| | - Solange M. T. Serrano
- Laboratório de Toxinologia Aplicada, Instituto Butantan, São Paulo 05503-900, SP, Brazil; (V.C.S.); (J.P.C.S.); (S.M.T.S.); (I.L.M.J.-d.-A.)
| | - Inácio L. M. Junqueira-de-Azevedo
- Laboratório de Toxinologia Aplicada, Instituto Butantan, São Paulo 05503-900, SP, Brazil; (V.C.S.); (J.P.C.S.); (S.M.T.S.); (I.L.M.J.-d.-A.)
| | - Felipe G. Grazziotin
- Laboratório de Coleções Zoológicas, Instituto Butantan, São Paulo 05503-900, SP, Brazil;
| | - Wuelton M. Monteiro
- Instituto de Pesquisa Clínica Carlos Borborema, Fundação de Medicina Tropical Dr. Heitor Vieira Dourado, Manaus 69040-000, AM, Brazil; (M.A.S.); (J.S.); (W.M.M.)
| | - Paulo S. Bernarde
- Laboratório de Herpetologia, Universidade Federal do Acre, Campus Floresta, Cruzeiro do Sul 69895-000, AC, Brazil; (M.R.N.A.); (R.A.M.); (W.L.F.); (P.S.B.)
| | - Ana M. Moura-da-Silva
- Laboratório de Imunopatologia, Instituto Butantan, São Paulo 05503-900, SP, Brazil; (L.A.F.-d.-S.); (M.C.)
| |
Collapse
|
4
|
Shin MK, Park HR, Hwang IW, Bu KB, Jang BY, Lee SH, Oh JW, Yoo JS, Sung JS. In Silico-Based Design of a Hybrid Peptide with Antimicrobial Activity against Multidrug-Resistant Pseudomonas aeruginosa Using a Spider Toxin Peptide. Toxins (Basel) 2023; 15:668. [PMID: 38133172 PMCID: PMC10747792 DOI: 10.3390/toxins15120668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Revised: 11/12/2023] [Accepted: 11/22/2023] [Indexed: 12/23/2023] Open
Abstract
The escalating prevalence of antibiotic-resistant bacteria poses an immediate and grave threat to public health. Antimicrobial peptides (AMPs) have gained significant attention as a promising alternative to conventional antibiotics. Animal venom comprises a diverse array of bioactive compounds, which can be a rich source for identifying new functional peptides. In this study, we identified a toxin peptide, Lycotoxin-Pa1a (Lytx-Pa1a), from the transcriptome of the Pardosa astrigera spider venom gland. To enhance its functional properties, we employed an in silico approach to design a novel hybrid peptide, KFH-Pa1a, by predicting antibacterial and cytotoxic functionalities and incorporating the amino-terminal Cu(II)- and Ni(II) (ATCUN)-binding motif. KFH-Pa1a demonstrated markedly superior antimicrobial efficacy against pathogens, including multidrug-resistant (MDR) Pseudomonas aeruginosa, compared to Lytx-Pa1a. Notably, KFH-Pa1a exerted several distinct mechanisms, including the disruption of the bacterial cytoplasmic membrane, the generation of intracellular ROS, and the cleavage and inhibition of bacterial DNA. Additionally, the hybrid peptide showed synergistic activity when combined with conventional antibiotics. Our research not only identified a novel toxin peptide from spider venom but demonstrated in silico-based design of hybrid AMP with strong antimicrobial activity that can contribute to combating MDR pathogens, broadening the utilization of biological resources by incorporating computational approaches.
Collapse
Affiliation(s)
- Min Kyoung Shin
- Department of Life Science, Dongguk University-Seoul, Goyang 10326, Republic of Korea; (M.K.S.); (H.-R.P.); (I.-W.H.); (K.-B.B.); (B.-Y.J.); (S.-H.L.); (J.W.O.)
| | - Hye-Ran Park
- Department of Life Science, Dongguk University-Seoul, Goyang 10326, Republic of Korea; (M.K.S.); (H.-R.P.); (I.-W.H.); (K.-B.B.); (B.-Y.J.); (S.-H.L.); (J.W.O.)
| | - In-Wook Hwang
- Department of Life Science, Dongguk University-Seoul, Goyang 10326, Republic of Korea; (M.K.S.); (H.-R.P.); (I.-W.H.); (K.-B.B.); (B.-Y.J.); (S.-H.L.); (J.W.O.)
| | - Kyung-Bin Bu
- Department of Life Science, Dongguk University-Seoul, Goyang 10326, Republic of Korea; (M.K.S.); (H.-R.P.); (I.-W.H.); (K.-B.B.); (B.-Y.J.); (S.-H.L.); (J.W.O.)
| | - Bo-Young Jang
- Department of Life Science, Dongguk University-Seoul, Goyang 10326, Republic of Korea; (M.K.S.); (H.-R.P.); (I.-W.H.); (K.-B.B.); (B.-Y.J.); (S.-H.L.); (J.W.O.)
| | - Seung-Ho Lee
- Department of Life Science, Dongguk University-Seoul, Goyang 10326, Republic of Korea; (M.K.S.); (H.-R.P.); (I.-W.H.); (K.-B.B.); (B.-Y.J.); (S.-H.L.); (J.W.O.)
| | - Jin Wook Oh
- Department of Life Science, Dongguk University-Seoul, Goyang 10326, Republic of Korea; (M.K.S.); (H.-R.P.); (I.-W.H.); (K.-B.B.); (B.-Y.J.); (S.-H.L.); (J.W.O.)
| | - Jung Sun Yoo
- Species Diversity Research Division, National Institute of Biological Resources, Incheon 22689, Republic of Korea;
| | - Jung-Suk Sung
- Department of Life Science, Dongguk University-Seoul, Goyang 10326, Republic of Korea; (M.K.S.); (H.-R.P.); (I.-W.H.); (K.-B.B.); (B.-Y.J.); (S.-H.L.); (J.W.O.)
| |
Collapse
|
5
|
Oliveira LD, Nachtigall PG, Vialla VL, Campos PF, Costa-Neves AD, Zaher H, Silva NJD, Grazziotin FG, Wilkinson M, Junqueira-de-Azevedo ILM. Comparing morphological and secretory aspects of cephalic glands among the New World coral snakes brings novel insights on their biological roles. Toxicon 2023; 234:107285. [PMID: 37683698 DOI: 10.1016/j.toxicon.2023.107285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Revised: 09/05/2023] [Accepted: 09/06/2023] [Indexed: 09/10/2023]
Abstract
Oral and other cephalic glands have been surveyed by several studies with distinct purposes. Despite the wide diversity and medical relevance of the New World coral snakes, studies focusing on understanding the biological roles of the glands within this group are still scarce. Specifically, the venom glands of some coral snakes were previously investigated but all other cephalic glands remain uncharacterized. In this sense, performing morphological and molecular analysis of these glands may help better understand their biological role. Here, we studied the morphology of the venom, infralabial, rictal, and harderian glands of thirteen species of Micrurus and Micruroides euryxanthus. We also performed a molecular characterization of these glands from selected species of Micrurus using transcriptomic and proteomic approaches. We described substantial morphological variation in the cephalic glands of New World coral snakes and structural evidence for protein-secreting cells in the inferior rictal glands. Our molecular analysis revealed that the venom glands, as expected, are majorly devoted to toxin production, however, the infralabial and inferior rictal glands also expressed some toxin genes at low to medium levels, despite the marked morphological differences. On the other hand, the harderian glands were dominated by the expression of lipocalins, but do not produce toxins. Our integrative analysis, including the prediction of biological processes and pathways, helped decipher some important traits of cephalic glands and better understand their biology.
Collapse
Affiliation(s)
- Leonardo de Oliveira
- Laboratório de Toxinologia Aplicada, Centre of Toxins, Immune-Response and Cell Signaling (CeTICS), Instituto Butantan, São Paulo, 05503-900, Brazil; Herpetology, The Natural History Museum, London, SW7 5BD, United Kingdom.
| | - Pedro Gabriel Nachtigall
- Laboratório de Toxinologia Aplicada, Centre of Toxins, Immune-Response and Cell Signaling (CeTICS), Instituto Butantan, São Paulo, 05503-900, Brazil
| | - Vincent Louis Vialla
- Laboratório de Toxinologia Aplicada, Centre of Toxins, Immune-Response and Cell Signaling (CeTICS), Instituto Butantan, São Paulo, 05503-900, Brazil
| | - Pollyanna F Campos
- Laboratório de Toxinologia Aplicada, Centre of Toxins, Immune-Response and Cell Signaling (CeTICS), Instituto Butantan, São Paulo, 05503-900, Brazil
| | | | - Hussam Zaher
- Museu de Zoologia da Universidade de São Paulo, Avenida Nazaré 481, Ipiranga, 04263-000, São Paulo, Brazil
| | - Nelson Jorge da Silva
- Programa de Pós-Graduação em Ciências Ambientais e Saúde, Pontifícia Universidade Católica de Goiás, Goiânia, Goiás, 74605-140, Brazil
| | - Felipe G Grazziotin
- Laboratório de Coleções Zoológicas, Instituto Butantan, São Paulo, 05503-900, Brazil
| | - Mark Wilkinson
- Herpetology, The Natural History Museum, London, SW7 5BD, United Kingdom
| | - Inácio L M Junqueira-de-Azevedo
- Laboratório de Toxinologia Aplicada, Centre of Toxins, Immune-Response and Cell Signaling (CeTICS), Instituto Butantan, São Paulo, 05503-900, Brazil
| |
Collapse
|
6
|
Rosales-García RA, Rautsaw RM, Hofmann EP, Grünwald CI, Franz-Chavez H, Ahumada-Carrillo IT, Ramirez-Chaparro R, de la Torre-Loranca MA, Strickland JL, Mason AJ, Holding ML, Borja M, Castañeda-Gaytan G, Myers EA, Sasa M, Rokyta DR, Parkinson CL. Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium. J Mol Evol 2023; 91:514-535. [PMID: 37269364 PMCID: PMC10995822 DOI: 10.1007/s00239-023-10115-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Accepted: 05/03/2023] [Indexed: 06/05/2023]
Abstract
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.
Collapse
Affiliation(s)
| | - Rhett M Rautsaw
- Department of Biological Sciences, Clemson University, 190 Collings St., Clemson, SC, 29634, USA
| | - Erich P Hofmann
- Department of Biological Sciences, Clemson University, 190 Collings St., Clemson, SC, 29634, USA
- Science Department, Cape Fear Community College, Wilmington, NC, 28401, USA
| | | | - Hector Franz-Chavez
- Herp.mx A.C., Colima, Mexico
- Biodiversa A. C., Chapala, Jalisco, 45900, Mexico
| | | | | | | | - Jason L Strickland
- Department of Biological Sciences, Clemson University, 190 Collings St., Clemson, SC, 29634, USA
- Department of Biology, University of South Alabama, Mobile, AL, 36688, USA
| | - Andrew J Mason
- Department of Biological Sciences, Clemson University, 190 Collings St., Clemson, SC, 29634, USA
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University, Columbus, OH, 43210, USA
| | - Matthew L Holding
- Department of Biological Sciences, Clemson University, 190 Collings St., Clemson, SC, 29634, USA
- Life Sciences Institute, University of Michigan, Ann Arbor, MI, 48109, USA
| | - Miguel Borja
- Facultad de Ciencias Biológicas, Universdad Juárez del Estado de Durango, Gómez Palacio, Durango, 35010, Mexico
| | - Gamaliel Castañeda-Gaytan
- Facultad de Ciencias Biológicas, Universdad Juárez del Estado de Durango, Gómez Palacio, Durango, 35010, Mexico
| | - Edward A Myers
- Department of Biological Sciences, Clemson University, 190 Collings St., Clemson, SC, 29634, USA
| | - Mahmood Sasa
- Centro Investigaciones en Biodiversidad y Ecología Tropical and Instituto Clodomiro Picado, Universidad de Costa Rica, San José, Costa Rica
| | - Darin R Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Christopher L Parkinson
- Department of Biological Sciences, Clemson University, 190 Collings St., Clemson, SC, 29634, USA.
| |
Collapse
|
7
|
Smith EG, Surm JM, Macrander J, Simhi A, Amir G, Sachkova MY, Lewandowska M, Reitzel AM, Moran Y. Micro and macroevolution of sea anemone venom phenotype. Nat Commun 2023; 14:249. [PMID: 36646703 PMCID: PMC9842752 DOI: 10.1038/s41467-023-35794-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2022] [Accepted: 01/03/2023] [Indexed: 01/18/2023] Open
Abstract
Venom is a complex trait with substantial inter- and intraspecific variability resulting from strong selective pressures acting on the expression of many toxic proteins. However, understanding the processes underlying toxin expression dynamics that determine the venom phenotype remains unresolved. By interspecific comparisons we reveal that toxin expression in sea anemones evolves rapidly and that in each species different toxin family dictates the venom phenotype by massive gene duplication events. In-depth analysis of the sea anemone, Nematostella vectensis, revealed striking variation of the dominant toxin (Nv1) diploid copy number across populations (1-24 copies) resulting from independent expansion/contraction events, which generate distinct haplotypes. Nv1 copy number correlates with expression at both the transcript and protein levels with one population having a near-complete loss of Nv1 production. Finally, we establish the dominant toxin hypothesis which incorporates observations in other venomous lineages that animals have convergently evolved a similar strategy in shaping their venom.
Collapse
Affiliation(s)
- Edward G Smith
- University of North Carolina at Charlotte, Department of Biological Sciences, Charlotte, NC, USA. .,School of Life Sciences, University of Warwick, Coventry, United Kingdom.
| | - Joachim M Surm
- Department of Ecology, Evolution and Behavior, Alexander Silberman Institute of Life Sciences, Faculty of Science, The Hebrew University of Jerusalem, Jerusalem, Israel.
| | - Jason Macrander
- University of North Carolina at Charlotte, Department of Biological Sciences, Charlotte, NC, USA.,Florida Southern College, Biology Department, Lakeland, FL, USA
| | - Adi Simhi
- Department of Ecology, Evolution and Behavior, Alexander Silberman Institute of Life Sciences, Faculty of Science, The Hebrew University of Jerusalem, Jerusalem, Israel.,The Hebrew University of Jerusalem, The School of Computer Science & Engineering, Jerusalem, Israel
| | - Guy Amir
- Department of Ecology, Evolution and Behavior, Alexander Silberman Institute of Life Sciences, Faculty of Science, The Hebrew University of Jerusalem, Jerusalem, Israel.,The Hebrew University of Jerusalem, The School of Computer Science & Engineering, Jerusalem, Israel
| | - Maria Y Sachkova
- Department of Ecology, Evolution and Behavior, Alexander Silberman Institute of Life Sciences, Faculty of Science, The Hebrew University of Jerusalem, Jerusalem, Israel.,Sars International Centre for Marine Molecular Biology, University of Bergen, Bergen, Norway
| | - Magda Lewandowska
- Department of Ecology, Evolution and Behavior, Alexander Silberman Institute of Life Sciences, Faculty of Science, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Adam M Reitzel
- University of North Carolina at Charlotte, Department of Biological Sciences, Charlotte, NC, USA
| | - Yehu Moran
- Department of Ecology, Evolution and Behavior, Alexander Silberman Institute of Life Sciences, Faculty of Science, The Hebrew University of Jerusalem, Jerusalem, Israel.
| |
Collapse
|
8
|
Schield DR, Perry BW, Adams RH, Holding ML, Nikolakis ZL, Gopalan SS, Smith CF, Parker JM, Meik JM, DeGiorgio M, Mackessy SP, Castoe TA. The roles of balancing selection and recombination in the evolution of rattlesnake venom. Nat Ecol Evol 2022; 6:1367-1380. [PMID: 35851850 PMCID: PMC9888523 DOI: 10.1038/s41559-022-01829-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Accepted: 06/15/2022] [Indexed: 02/02/2023]
Abstract
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.
Collapse
Affiliation(s)
- Drew R Schield
- Department of Biology, University of Texas at Arlington, Arlington, TX, USA.
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, USA.
| | - Blair W Perry
- Department of Biology, University of Texas at Arlington, Arlington, TX, USA
- School of Biological Sciences, Washington State University, Pullman, WA, USA
| | - Richard H Adams
- Department of Biological and Environmental Sciences, Georgia College and State University, Milledgeville, GA, USA
| | | | | | | | - Cara F Smith
- School of Biological Sciences, University of Northern Colorado, Greeley, CO, USA
| | - Joshua M Parker
- Life Science Department, Fresno City College, Fresno, CA, USA
| | - Jesse M Meik
- Department of Biological Sciences, Tarleton State University, Stephenville, TX, USA
| | - Michael DeGiorgio
- Department of Electrical Engineering and Computer Science, Florida Atlantic University, Boca Raton, FL, USA
| | - Stephen P Mackessy
- School of Biological Sciences, University of Northern Colorado, Greeley, CO, USA
| | - Todd A Castoe
- Department of Biology, University of Texas at Arlington, Arlington, TX, USA.
| |
Collapse
|
9
|
Myers EA, Strickland JL, Rautsaw RM, Mason AJ, Schramer TD, Nystrom GS, Hogan MP, Yooseph S, Rokyta DR, Parkinson CL. De Novo Genome Assembly Highlights the Role of Lineage-Specific Duplications in the Evolution of Venom in Fea’s Viper. Genome Biol Evol 2022; 14:6603630. [PMID: 35670514 PMCID: PMC9256536 DOI: 10.1093/gbe/evac082] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/25/2022] [Indexed: 11/12/2022] Open
Abstract
Despite the medical significance to humans and important ecological roles filled by vipers, few high-quality genomic resources exist for these snakes outside of a few genera of pitvipers. Here we sequence, assemble, and annotate the genome of Fea’s Viper (Azemiops feae). This taxon is distributed in East Asia and belongs to a monotypic subfamily, sister to the pitvipers. The newly sequenced genome resulted in a 1.56 Gb assembly, a contig N50 of 1.59 Mb, with 97.6% of the genome assembly in contigs >50 Kb, and a BUSCO completeness of 92.4%. We found that A. feae venom is primarily composed of phospholipase A2 (PLA2) proteins expressed by genes that likely arose from lineage-specific PLA2 gene duplications. Additionally, we show that renin, an enzyme associated with blood pressure regulation in mammals and known from the venoms of two viper species including A. feae, is expressed in the venom gland at comparative levels to known toxins and is present in the venom proteome. The cooption of this gene as a toxin may be more widespread in viperids than currently known. To investigate the historical population demographics of A. feae, we performed coalescent-based analyses and determined that the effective population size has remained stable over the last 100 kyr. This suggests Quaternary glacial cycles likely had minimal influence on the demographic history of A. feae. This newly assembled genome will be an important resource for studying the genomic basis of phenotypic evolution and understanding the diversification of venom toxin gene families.
Collapse
Affiliation(s)
- Edward A. Myers
- Department of Biological Sciences, Clemson University , Clemson, SC 29634, USA
| | - Jason L. Strickland
- Department of Biological Sciences, Clemson University , Clemson, SC 29634, USA
- Department of Biology, University of South Alabama , Mobile, AL 36688, USA
| | - Rhett M. Rautsaw
- Department of Biological Sciences, Clemson University , Clemson, SC 29634, USA
| | - Andrew J. Mason
- Department of Biological Sciences, Clemson University , Clemson, SC 29634, USA
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University , Columbus, OH 43210, USA
| | - Tristan D. Schramer
- Department of Biological Sciences, Clemson University , Clemson, SC 29634, USA
| | - Gunnar S. Nystrom
- Department of Biological Science, Florida State University , Tallahassee, FL 32306, USA
| | - Michael P. Hogan
- Department of Biological Science, Florida State University , Tallahassee, FL 32306, USA
| | - Shibu Yooseph
- Department of Computer Science, Genomics and Bioinformatics Cluster, University of Central Florida, 4000 Central Florida Blvd , Orlando, FL, 32816, USA
| | - Darin R. Rokyta
- Department of Biological Science, Florida State University , Tallahassee, FL 32306, USA
| | - Christopher L. Parkinson
- Department of Biological Sciences, Clemson University , Clemson, SC 29634, USA
- Department of Forestry and Environmental Conservation, Clemson University , Clemson, SC 29634, USA
| |
Collapse
|
10
|
VenomMaps: Updated species distribution maps and models for New World pitvipers (Viperidae: Crotalinae). Sci Data 2022; 9:232. [PMID: 35614080 PMCID: PMC9132920 DOI: 10.1038/s41597-022-01323-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Accepted: 04/19/2022] [Indexed: 11/19/2022] Open
Abstract
Beyond providing critical information to biologists, species distributions are useful for naturalists, curious citizens, and applied disciplines including conservation planning and medical intervention. Venomous snakes are one group that highlight the importance of having accurate information given their cosmopolitan distribution and medical significance. Envenomation by snakebite is considered a neglected tropical disease by the World Health Organization and venomous snake distributions are used to assess vulnerability to snakebite based on species occurrence and antivenom/healthcare accessibility. However, recent studies highlighted the need for updated fine-scale distributions of venomous snakes. Pitvipers (Viperidae: Crotalinae) are responsible for >98% of snakebites in the New World. Therefore, to begin to address the need for updated fine-scale distributions, we created VenomMaps, a database and web application containing updated distribution maps and species distribution models for all species of New World pitvipers. With these distributions, biologists can better understand the biogeography and conservation status of this group, researchers can better assess vulnerability to snakebite, and medical professionals can easily discern species found in their area. Measurement(s) | Species Distributions | Technology Type(s) | Geographic Information System • Species Distribution Model (MaxEnt/kuenm) | Factor Type(s) | Occurrence Records • Environmental Data | Sample Characteristic - Organism | Crotalinae | Sample Characteristic - Location | North America • South America |
Collapse
|
11
|
Differences in PLA2 Constitution Distinguish the Venom of Two Endemic Brazilian Mountain Lanceheads, Bothrops cotiara and Bothrops fonsecai. Toxins (Basel) 2022; 14:toxins14040237. [PMID: 35448846 PMCID: PMC9028134 DOI: 10.3390/toxins14040237] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Revised: 02/15/2022] [Accepted: 02/23/2022] [Indexed: 02/01/2023] Open
Abstract
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.
Collapse
|
12
|
Hofmann EP, Rautsaw RM, Mason AJ, Strickland JL, Parkinson CL. Duvernoy's Gland Transcriptomics of the Plains Black-Headed Snake, Tantilla nigriceps (Squamata, Colubridae): Unearthing the Venom of Small Rear-Fanged Snakes. Toxins (Basel) 2021; 13:336. [PMID: 34066626 PMCID: PMC8148590 DOI: 10.3390/toxins13050336] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Revised: 05/03/2021] [Accepted: 05/04/2021] [Indexed: 12/13/2022] Open
Abstract
The venoms of small rear-fanged snakes (RFS) remain largely unexplored, despite increased recognition of their importance in understanding venom evolution more broadly. Sequencing the transcriptome of venom-producing glands has greatly increased the ability of researchers to examine and characterize the toxin repertoire of small taxa with low venom yields. Here, we use RNA-seq to characterize the Duvernoy's gland transcriptome of the Plains Black-headed Snake, Tantilla nigriceps, a small, semi-fossorial colubrid that feeds on a variety of potentially dangerous arthropods including centipedes and spiders. We generated transcriptomes of six individuals from three localities in order to both characterize the toxin expression of this species for the first time, and to look for initial evidence of venom variation in the species. Three toxin families-three-finger neurotoxins (3FTxs), cysteine-rich secretory proteins (CRISPs), and snake venom metalloproteinases (SVMPIIIs)-dominated the transcriptome of T. nigriceps; 3FTx themselves were the dominant toxin family in most individuals, accounting for as much as 86.4% of an individual's toxin expression. Variation in toxin expression between individuals was also noted, with two specimens exhibiting higher relative expression of c-type lectins than any other sample (8.7-11.9% compared to <1%), and another expressed CRISPs higher than any other toxin. This study provides the first Duvernoy's gland transcriptomes of any species of Tantilla, and one of the few transcriptomic studies of RFS not predicated on a single individual. This initial characterization demonstrates the need for further study of toxin expression variation in this species, as well as the need for further exploration of small RFS venoms.
Collapse
Affiliation(s)
- Erich P. Hofmann
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA; (E.P.H.); (R.M.R.); (A.J.M.); (J.L.S.)
| | - Rhett M. Rautsaw
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA; (E.P.H.); (R.M.R.); (A.J.M.); (J.L.S.)
| | - Andrew J. Mason
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA; (E.P.H.); (R.M.R.); (A.J.M.); (J.L.S.)
| | - Jason L. Strickland
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA; (E.P.H.); (R.M.R.); (A.J.M.); (J.L.S.)
| | - Christopher L. Parkinson
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA; (E.P.H.); (R.M.R.); (A.J.M.); (J.L.S.)
- Department of Forestry and Environmental Conservation, Clemson University, Clemson, SC 29634, USA
| |
Collapse
|
13
|
Holding ML, Strickland JL, Rautsaw RM, Hofmann EP, Mason AJ, Hogan MP, Nystrom GS, Ellsworth SA, Colston TJ, Borja M, Castañeda-Gaytán G, Grünwald CI, Jones JM, Freitas-de-Sousa LA, Viala VL, Margres MJ, Hingst-Zaher E, Junqueira-de-Azevedo ILM, Moura-da-Silva AM, Grazziotin FG, Gibbs HL, Rokyta DR, Parkinson CL. Phylogenetically diverse diets favor more complex venoms in North American pitvipers. Proc Natl Acad Sci U S A 2021; 118:e2015579118. [PMID: 33875585 PMCID: PMC8092465 DOI: 10.1073/pnas.2015579118] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Matthew L Holding
- Department of Biological Sciences, Clemson University, Clemson, SC 29634;
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Jason L Strickland
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
| | - Rhett M Rautsaw
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
| | - Erich P Hofmann
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
| | - Andrew J Mason
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University, Columbus, OH 43210
| | - Michael P Hogan
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Gunnar S Nystrom
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Schyler A Ellsworth
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Timothy J Colston
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Miguel Borja
- Facultad de Ciencias Biológicas, Universidad Juárez del Estado de Durango, C.P. 35010 Gómez Palacio, Dgo., Mexico
| | - Gamaliel Castañeda-Gaytán
- Facultad de Ciencias Biológicas, Universidad Juárez del Estado de Durango, C.P. 35010 Gómez Palacio, Dgo., Mexico
| | | | - Jason M Jones
- HERP.MX A.C., Villa del Álvarez, Colima 28973, Mexico
| | | | - Vincent Louis Viala
- Laboratório de Toxinologia Aplicada, Instituto Butantan, São Paulo 05503-900, Brazil
- Center of Toxins, Immune-Response and Cell Signaling, São Paulo 05503-900, Brazil
| | - Mark J Margres
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138
| | | | - Inácio L M Junqueira-de-Azevedo
- Laboratório de Toxinologia Aplicada, Instituto Butantan, São Paulo 05503-900, Brazil
- Center of Toxins, Immune-Response and Cell Signaling, São Paulo 05503-900, Brazil
| | - Ana M Moura-da-Silva
- Laboratório de Imunopatologia, Instituto Butantan, São Paulo 05503-900, Brazil
- Instituto de Pesquisa Clínica Carlos Borborema, Fundação de Medicina Tropical Doutor Heitor Vieira Dourado, Manaus 69040, Brazil
| | - Felipe G Grazziotin
- Laboratório de Coleções Zoológicas, Instituto Butantan, São Paulo 05503-900, Brazil
| | - H Lisle Gibbs
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University, Columbus, OH 43210
| | - Darin R Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Christopher L Parkinson
- Department of Biological Sciences, Clemson University, Clemson, SC 29634;
- Department of Forestry and Environmental Conservation, Clemson University, Clemson, SC 29634
| |
Collapse
|
14
|
Nachtigall PG, Rautsaw RM, Ellsworth SA, Mason AJ, Rokyta DR, Parkinson CL, Junqueira-de-Azevedo ILM. ToxCodAn: a new toxin annotator and guide to venom gland transcriptomics. Brief Bioinform 2021; 22:6235957. [PMID: 33866357 DOI: 10.1093/bib/bbab095] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2021] [Revised: 02/15/2021] [Accepted: 03/03/2021] [Indexed: 01/23/2023] Open
Abstract
MOTIVATION Next-generation sequencing has become exceedingly common and has transformed our ability to explore nonmodel systems. In particular, transcriptomics has facilitated the study of venom and evolution of toxins in venomous lineages; however, many challenges remain. Primarily, annotation of toxins in the transcriptome is a laborious and time-consuming task. Current annotation software often fails to predict the correct coding sequence and overestimates the number of toxins present in the transcriptome. Here, we present ToxCodAn, a python script designed to perform precise annotation of snake venom gland transcriptomes. We test ToxCodAn with a set of previously curated transcriptomes and compare the results to other annotators. In addition, we provide a guide for venom gland transcriptomics to facilitate future research and use Bothrops alternatus as a case study for ToxCodAn and our guide. RESULTS Our analysis reveals that ToxCodAn provides precise annotation of toxins present in the transcriptome of venom glands of snakes. Comparison with other annotators demonstrates that ToxCodAn has better performance with regard to run time ($>20x$ faster), coding sequence prediction ($>3x$ more accurate) and the number of toxins predicted (generating $>4x$ less false positives). In this sense, ToxCodAn is a valuable resource for toxin annotation. The ToxCodAn framework can be expanded in the future to work with other venomous lineages and detect novel toxins.
Collapse
Affiliation(s)
- Pedro G Nachtigall
- Laboratório de Toxinologia Aplicada, CeTICS, Instituto Butantan, São Paulo, SP 05503-900, Brazil
| | - Rhett M Rautsaw
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA
| | - Schyler A Ellsworth
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA
| | - Andrew J Mason
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University, Columbus, OH 43210 USA
| | - Darin R Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA
| | - Christopher L Parkinson
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA
- Department of Forestry and Environmental Conservation, Clemson University, Clemson, SC 29634, USA
| | | |
Collapse
|
15
|
Abstract
Although oral venom systems are ecologically important characters, how they originated is still unclear. In this study, we show that oral venom systems likely originated from a gene regulatory network conserved across amniotes. This network, which we term the “metavenom network,” comprises over 3,000 housekeeping genes coexpressed with venom and play a role in protein folding and modification. Comparative transcriptomics revealed that the network is conserved between venom glands of snakes and salivary glands of mammals. This suggests that while these tissues have evolved different functions, they share a common regulatory core, that persisted since their common ancestor. We propose several evolutionary mechanisms that can utilize this common regulatory core to give rise to venomous animals from their nonvenomous ancestors. Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.
Collapse
|
16
|
Lima EOVD, Tasima LJ, Hatakeyama DM, Serino-Silva C, Rodrigues CFB, Galizio NDC, Chiarelli T, Nishiduka ES, Rocha MMTD, Sant'Anna SS, Grego KF, Tashima AK, Tanaka-Azevedo AM, Morais-Zani KD. Snake venom color and L-amino acid oxidase: An evidence of long-term captive Crotalus durissus terrificus venom plasticity. Toxicon 2021; 193:73-83. [PMID: 33515573 DOI: 10.1016/j.toxicon.2021.01.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Revised: 12/08/2020] [Accepted: 01/20/2021] [Indexed: 12/01/2022]
Abstract
The venom color variation of Crotalus durissus terrificus (Cdt) is attributed to the presence of the toxin L-amino acid oxidase (LAAO). During the venom milking routine of Instituto Butantan, we have noticed that most venoms of captive Cdt specimens show a yellowish color, while most venoms of wild specimens are white. Here we describe a comparative analysis of long-term captive (LTC) and recently wild-caught (RWC) Cdt, focusing on LAAO variation. For the identification of LAAO in individual venoms, four different approaches were employed: evaluation of the enzymatic activity, SDS-PAGE, Western blotting, and ELISA. In addition, mass spectrometry analysis was performed using pooled samples. Although some variation among these methodologies was observed, it was possible to notice that the presence of LAAO was significantly higher in the venom of LTC individuals. LAAO was identified in 60-80% LTC specimens and in only 10-12% of RWC specimens. Furthermore, this enzyme accounts for 5.6% of total venom proteins of LTC Cdt pooled venom, while it corresponds to only 0.7% of RWC Cdt pooled venom. These findings strongly suggest that captive maintenance increases the expression of LAAO in Cdt venom.
Collapse
Affiliation(s)
| | - Lídia Jorge Tasima
- Laboratório de Herpetologia, Instituto Butantan, São Paulo, Brazil; Interunidades em Biotecnologia, Universidade de São Paulo, Instituto de Pesquisas Tecnológicas, Instituto Butantan, São Paulo, Brazil
| | - Daniela Miki Hatakeyama
- Laboratório de Herpetologia, Instituto Butantan, São Paulo, Brazil; Interunidades em Biotecnologia, Universidade de São Paulo, Instituto de Pesquisas Tecnológicas, Instituto Butantan, São Paulo, Brazil
| | - Caroline Serino-Silva
- Laboratório de Herpetologia, Instituto Butantan, São Paulo, Brazil; Interunidades em Biotecnologia, Universidade de São Paulo, Instituto de Pesquisas Tecnológicas, Instituto Butantan, São Paulo, Brazil
| | - Caroline Fabri Bittencourt Rodrigues
- Laboratório de Herpetologia, Instituto Butantan, São Paulo, Brazil; Interunidades em Biotecnologia, Universidade de São Paulo, Instituto de Pesquisas Tecnológicas, Instituto Butantan, São Paulo, Brazil
| | - Nathália da Costa Galizio
- Laboratório de Herpetologia, Instituto Butantan, São Paulo, Brazil; Interunidades em Biotecnologia, Universidade de São Paulo, Instituto de Pesquisas Tecnológicas, Instituto Butantan, São Paulo, Brazil
| | - Tassia Chiarelli
- Departamento de Bioquímica, Universidade Federal de São Paulo, São Paulo, Brazil
| | | | | | | | | | | | - Anita Mitico Tanaka-Azevedo
- Laboratório de Herpetologia, Instituto Butantan, São Paulo, Brazil; Interunidades em Biotecnologia, Universidade de São Paulo, Instituto de Pesquisas Tecnológicas, Instituto Butantan, São Paulo, Brazil
| | - Karen de Morais-Zani
- Laboratório de Herpetologia, Instituto Butantan, São Paulo, Brazil; Interunidades em Biotecnologia, Universidade de São Paulo, Instituto de Pesquisas Tecnológicas, Instituto Butantan, São Paulo, Brazil.
| |
Collapse
|
17
|
Freitas-de-Sousa LA, Nachtigall PG, Portes-Junior JA, Holding ML, Nystrom GS, Ellsworth SA, Guimarães NC, Tioyama E, Ortiz F, Silva BR, Kunz TS, Junqueira-de-Azevedo ILM, Grazziotin FG, Rokyta DR, Moura-da-Silva AM. Size Matters: An Evaluation of the Molecular Basis of Ontogenetic Modifications in the Composition of Bothrops jararacussu Snake Venom. Toxins (Basel) 2020; 12:toxins12120791. [PMID: 33322460 PMCID: PMC7763748 DOI: 10.3390/toxins12120791] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Revised: 12/03/2020] [Accepted: 12/07/2020] [Indexed: 12/13/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Luciana A. Freitas-de-Sousa
- Programa de Pós-Graduação em Ciências-Toxinologia, Laboratório de Imunopatologia, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (N.C.G.); (E.T.)
- Correspondence: (L.A.F.-d.-S.); (A.M.M.-d.-S.); Tel.: +55-11-2627-9779 (A.M.M.-d.-S.)
| | - Pedro G. Nachtigall
- Laboratório Especial de Toxinologia Aplicada, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (P.G.N.); (I.L.M.J.-d.-A.)
| | - José A. Portes-Junior
- Laboratório de Coleções Zoológicas, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (J.A.P.-J.); (F.O.); (B.R.S.); (T.S.K.); (F.G.G.)
| | - Matthew L. Holding
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA; (M.L.H.); (G.S.N.); (S.A.E.); (D.R.R.)
| | - Gunnar S. Nystrom
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA; (M.L.H.); (G.S.N.); (S.A.E.); (D.R.R.)
| | - Schyler A. Ellsworth
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA; (M.L.H.); (G.S.N.); (S.A.E.); (D.R.R.)
| | - Noranathan C. Guimarães
- Programa de Pós-Graduação em Ciências-Toxinologia, Laboratório de Imunopatologia, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (N.C.G.); (E.T.)
| | - Emilly Tioyama
- Programa de Pós-Graduação em Ciências-Toxinologia, Laboratório de Imunopatologia, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (N.C.G.); (E.T.)
| | - Flora Ortiz
- Laboratório de Coleções Zoológicas, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (J.A.P.-J.); (F.O.); (B.R.S.); (T.S.K.); (F.G.G.)
| | - Bruno R. Silva
- Laboratório de Coleções Zoológicas, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (J.A.P.-J.); (F.O.); (B.R.S.); (T.S.K.); (F.G.G.)
| | - Tobias S. Kunz
- Laboratório de Coleções Zoológicas, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (J.A.P.-J.); (F.O.); (B.R.S.); (T.S.K.); (F.G.G.)
| | | | - Felipe G. Grazziotin
- Laboratório de Coleções Zoológicas, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (J.A.P.-J.); (F.O.); (B.R.S.); (T.S.K.); (F.G.G.)
| | - Darin R. Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA; (M.L.H.); (G.S.N.); (S.A.E.); (D.R.R.)
| | - Ana M. Moura-da-Silva
- Programa de Pós-Graduação em Ciências-Toxinologia, Laboratório de Imunopatologia, Instituto Butantan, 05503-900 São Paulo, SP, Brazil; (N.C.G.); (E.T.)
- Instituto de Pesquisa Clínica Carlos Borborema, Fundação de Medicina Tropical Dr. Heitor Vieira Dourado, 69040-000 Manaus, AM, Brazil
- Correspondence: (L.A.F.-d.-S.); (A.M.M.-d.-S.); Tel.: +55-11-2627-9779 (A.M.M.-d.-S.)
| |
Collapse
|
18
|
Lewis JJ, Van Belleghem SM. Mechanisms of Change: A Population-Based Perspective on the Roles of Modularity and Pleiotropy in Diversification. Front Ecol Evol 2020. [DOI: 10.3389/fevo.2020.00261] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
|