1
|
Carmelet-Rescan D, Morgan-Richards M, Trewick SA. Metabolic differentiation of brushtail possum populations resistant and susceptible to plant toxins revealed via differential gene expression. J Comp Physiol B 2024:10.1007/s00360-024-01591-z. [PMID: 39495241 DOI: 10.1007/s00360-024-01591-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2024] [Revised: 09/02/2024] [Accepted: 10/14/2024] [Indexed: 11/05/2024]
Abstract
The Australian brushtail possum (Trichosurus vulpecula) is adapted to a wide range of food plants across its range and is exposed to numerous physiological challenges. Populations that are resistant to the plant toxin sodium fluoroacetate are of particular interest as this compound has been used since the 1940s for vertebrate pest management around the world. Candidate gene identification is an important first step in understanding how spatial populations have responded to local selection resulting in local physiological divergence. We employ differential gene expression of liver samples from wild-caught brushtail possums from toxin-resistant and toxin-susceptible populations to identify candidate genes that might be involved in metabolic pathways associated with toxin-resistance. This allowed us to identify genetic pathways involved in resistance to the plant toxin sodium fluoroacetate in Western Australian possums but not those originally from south eastern Australia. We identified differentially expressed genes in the liver that are associated with cell signalling, encapsulating structure, cell mobility, and tricarboxylic acid cycle. The gene expression differences detected indicate which metabolic pathways are most likely to be associated with sodium fluoroacetate resistance in these marsupials and we provide a comprehensive list of candidate genes and pathways to focus on for future studies.
Collapse
Affiliation(s)
- David Carmelet-Rescan
- Wildlife and Ecology, School of Natural Sciences, Massey University, Private Bag 11-222, Palmerston North, New Zealand.
| | - Mary Morgan-Richards
- Wildlife and Ecology, School of Natural Sciences, Massey University, Private Bag 11-222, Palmerston North, New Zealand
| | - Steven A Trewick
- Wildlife and Ecology, School of Natural Sciences, Massey University, Private Bag 11-222, Palmerston North, New Zealand
| |
Collapse
|
2
|
Hirst SR, Rautsaw RM, VanHorn CM, Beer MA, McDonald PJ, Rosales García RA, Rodriguez Lopez B, Rubio Rincón A, Franz Chávez H, Vásquez-Cruz V, Kelly Hernández A, Storfer A, Borja M, Castañeda-Gaytán G, Frandsen PB, Parkinson CL, Strickland JL, Margres MJ. Where the "ruber" Meets the Road: Using the Genome of the Red Diamond Rattlesnake to Unravel the Evolutionary Processes Driving Venom Evolution. Genome Biol Evol 2024; 16:evae198. [PMID: 39255072 PMCID: PMC11440179 DOI: 10.1093/gbe/evae198] [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: 06/07/2024] [Revised: 08/15/2024] [Accepted: 09/02/2024] [Indexed: 09/12/2024] Open
Abstract
Understanding the proximate and ultimate causes of phenotypic variation is fundamental in evolutionary research, as such variation provides the substrate for selection to act upon. Although trait variation can arise due to selection, the importance of neutral processes is sometimes understudied. We presented the first reference-quality genome of the Red Diamond Rattlesnake (Crotalus ruber) and used range-wide 'omic data to estimate the degree to which neutral and adaptive evolutionary processes shaped venom evolution. We characterized population structure and found substantial genetic differentiation across two populations, each with distinct demographic histories. We identified significant differentiation in venom expression across age classes with substantially reduced but discernible differentiation across populations. We then used conditional redundancy analysis to test whether venom expression variation was best predicted by neutral divergence patterns or geographically variable (a)biotic factors. Snake size was the most significant predictor of venom variation, with environment, prey availability, and neutral sequence variation also identified as significant factors, though to a lesser degree. By directly including neutrality in the model, our results confidently highlight the predominant, yet not singular, role of life history in shaping venom evolution.
Collapse
Affiliation(s)
- Samuel R Hirst
- Department of Integrative Biology, University of South Florida, Tampa, FL, USA
| | - Rhett M Rautsaw
- Department of Integrative Biology, University of South Florida, Tampa, FL, USA
- School of Biological Sciences, Washington State University, Pullman, WA, USA
| | - Cameron M VanHorn
- Department of Integrative Biology, University of South Florida, Tampa, FL, USA
| | - Marc A Beer
- School of Biological Sciences, Washington State University, Pullman, WA, USA
| | - Preston J McDonald
- Department of Integrative Biology, University of South Florida, Tampa, FL, USA
| | | | - Bruno Rodriguez Lopez
- Facultad de Ciencias Biológicas, Universidad Juárez del Estado de Durango, Durango, Mexico
| | - Alexandra Rubio Rincón
- Facultad de Ciencias Biológicas, Universidad Juárez del Estado de Durango, Durango, Mexico
| | | | - Víctor Vásquez-Cruz
- Facultad de Ciencias Biológicas y Agropecuarias, Universidad Veracruzana, Veracruz, Mexico
- PIMVS Herpetario Palancoatl, Veracruz, Mexico
| | | | - Andrew Storfer
- School of Biological Sciences, Washington State University, Pullman, WA, USA
| | - Miguel Borja
- Facultad de Ciencias Biológicas, Universidad Juárez del Estado de Durango, Durango, Mexico
| | | | - Paul B Frandsen
- Department of Plant and Wildlife Sciences, Brigham Young University, Provo, UT, USA
| | | | | | - Mark J Margres
- Department of Integrative Biology, University of South Florida, Tampa, FL, USA
| |
Collapse
|
3
|
Ellsworth SA, Rautsaw RM, Ward MJ, Holding ML, Rokyta DR. Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes. J Mol Evol 2024:10.1007/s00239-024-10191-y. [PMID: 39026042 DOI: 10.1007/s00239-024-10191-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2024] [Accepted: 07/09/2024] [Indexed: 07/20/2024]
Abstract
Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.
Collapse
Affiliation(s)
- Schyler A Ellsworth
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Rhett M Rautsaw
- Department of Integrative Biology, University of South Florida, Tampa, FL, 33620, USA
- School of Biological Sciences, Washington State University, Pullman, WA, 99164, USA
| | - Micaiah J Ward
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Matthew L Holding
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
- Life Sciences Institute, University of Michigan, Ann Arbor, MI, 48109, USA
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, 48109, USA
| | - Darin R Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA.
| |
Collapse
|
4
|
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
|
5
|
Salabi F, Vazirianzadeh B, Baradaran M. Identification, classification, and characterization of alpha and beta subunits of LVP1 protein from the venom gland of four Iranian scorpion species. Sci Rep 2023; 13:22277. [PMID: 38097679 PMCID: PMC10721645 DOI: 10.1038/s41598-023-49556-6] [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: 06/11/2023] [Accepted: 12/09/2023] [Indexed: 12/17/2023] Open
Abstract
Scorpion venoms contain bioactive peptides and proteins. Some, can be used for pharmaceutical purposes. So, identification of venom proteins matters because, in addition to determining the function of the toxins can also be an excellent guide to developing new drugs. Here, we got transcriptome of venom glands from four Iranian scorpion species, including Hemsicorpius lepturus, Mesobuthus eupeus, Andructunus crassicuada, and Hottentotta saulcyi using cDNA library synthesis and high-throughput transcriptomic analysis of the venom glands. In a comparative way, we identified the cDNA encoding isoforms of subunits (alpha and beta) of BotLVP1/BmLVP1-like protein in the venom gland of three species except for H. lepturus. Characterization and structure determination of the LVP1_like proteins combined with gene map analysis provided evidence of the existence of some isoforms of LVP1_like proteins, encoded by genes with two exons and one intron, which can be classified in CSαβ superfamily in the venom gland of three Iranian scorpion species. According to the high similarity with BotLVP1 and BmLVP1, these proteins could also be potent to mediate cholesterol homeostasis. However, further research is needed to prove it, and this study just may lay the foundation lead to light up this way.
Collapse
Affiliation(s)
- Fatemeh Salabi
- Razi Vaccine and Serum Research Institute, Agricultural Research, Education and Extension Organization (AREEO), Ahvaz, Iran
| | - Babak Vazirianzadeh
- Social Determinant of Health Research Center, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran
| | - Masoumeh Baradaran
- Toxicology Research Center, Medical Basic Sciences Research Institute, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran.
| |
Collapse
|
6
|
Koludarov I, Velasque M, Senoner T, Timm T, Greve C, Hamadou AB, Gupta DK, Lochnit G, Heinzinger M, Vilcinskas A, Gloag R, Harpur BA, Podsiadlowski L, Rost B, Jackson TNW, Dutertre S, Stolle E, von Reumont BM. Prevalent bee venom genes evolved before the aculeate stinger and eusociality. BMC Biol 2023; 21:229. [PMID: 37867198 PMCID: PMC10591384 DOI: 10.1186/s12915-023-01656-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2023] [Accepted: 06/29/2023] [Indexed: 10/24/2023] Open
Abstract
BACKGROUND Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected. RESULTS Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification. CONCLUSIONS Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.
Collapse
Affiliation(s)
- Ivan Koludarov
- Justus Liebig University of Gießen, Institute for Insect Biotechnology, Heinrich-Buff-Ring 58, 35392, Giessen, Germany.
- Department of Informatics, Bioinformatics and Computational Biology, i12, Technical University of Munich, Boltzmannstr. 3, Garching, 85748, Munich, Germany.
| | - Mariana Velasque
- Genomics & Regulatory Systems Unit, Okinawa Institute of Science & Technology, Tancha, Okinawa, 1919, Japan
| | - Tobias Senoner
- Department of Informatics, Bioinformatics and Computational Biology, i12, Technical University of Munich, Boltzmannstr. 3, Garching, 85748, Munich, Germany
| | - Thomas Timm
- Protein Analytics, Institute of Biochemistry, Justus Liebig University, Friedrichstrasse 24, 35392, Giessen, Germany
| | - Carola Greve
- LOEWE Centre for Translational Biodiversity Genomics (TBG), Senckenberganlage 25, 60325, Frankfurt, Germany
| | - Alexander Ben Hamadou
- LOEWE Centre for Translational Biodiversity Genomics (TBG), Senckenberganlage 25, 60325, Frankfurt, Germany
| | - Deepak Kumar Gupta
- LOEWE Centre for Translational Biodiversity Genomics (TBG), Senckenberganlage 25, 60325, Frankfurt, Germany
| | - Günter Lochnit
- Protein Analytics, Institute of Biochemistry, Justus Liebig University, Friedrichstrasse 24, 35392, Giessen, Germany
| | - Michael Heinzinger
- Department of Informatics, Bioinformatics and Computational Biology, i12, Technical University of Munich, Boltzmannstr. 3, Garching, 85748, Munich, Germany
| | - Andreas Vilcinskas
- Justus Liebig University of Gießen, Institute for Insect Biotechnology, Heinrich-Buff-Ring 58, 35392, Giessen, Germany
- Fraunhofer Institute for Molecular Biology and Applied Ecology, Department of Bioresources, Ohlebergsweg 12, 35392, Giessen, Germany
| | - Rosalyn Gloag
- Rosalyn Gloag - School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, 2006, Australia
| | - Brock A Harpur
- Brock A. Harpur - Department of Entomology, Purdue University, 901 W. State Street, West Lafayette, IN, 47907, USA
| | - Lars Podsiadlowski
- Leibniz Institute for the Analysis of Biodiversity Change, Zoological Research Museum Alexander Koenig, Centre of Molecular Biodiversity Research, Adenauerallee 160, 53113, Bonn, Germany
| | - Burkhard Rost
- Department of Informatics, Bioinformatics and Computational Biology, i12, Technical University of Munich, Boltzmannstr. 3, Garching, 85748, Munich, Germany
| | - Timothy N W Jackson
- Australian Venom Research Unit, Department of Biochemistry and Pharmacology, University of Melbourne, Grattan Street, Parkville, Viktoria, 3010, Australia
| | | | - Eckart Stolle
- Leibniz Institute for the Analysis of Biodiversity Change, Zoological Research Museum Alexander Koenig, Centre of Molecular Biodiversity Research, Adenauerallee 160, 53113, Bonn, Germany
| | - Björn M von Reumont
- LOEWE Centre for Translational Biodiversity Genomics (TBG), Senckenberganlage 25, 60325, Frankfurt, Germany.
- Faculty of Biological Sciences, Group of Applied Bioinformatics, Goethe University Frankfurt, Max-Von-Laue Str. 13, 60438, Frankfurt, Germany.
| |
Collapse
|
7
|
Ye X, He C, Yang Y, Sun YH, Xiong S, Chan KC, Si Y, Xiao S, Zhao X, Lin H, Mei Y, Yao Y, Ye G, Wu F, Fang Q. Comprehensive isoform-level analysis reveals the contribution of alternative isoforms to venom evolution and repertoire diversity. Genome Res 2023; 33:1554-1567. [PMID: 37798117 PMCID: PMC10620052 DOI: 10.1101/gr.277707.123] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Accepted: 08/08/2023] [Indexed: 10/07/2023]
Abstract
Animal venom systems have emerged as valuable models for investigating how novel polygenic phenotypes may arise from gene evolution by varying molecular mechanisms. However, a significant portion of venom genes produce alternative mRNA isoforms that have not been extensively characterized, hindering a comprehensive understanding of venom biology. In this study, we present a full-length isoform-level profiling workflow integrating multiple RNA sequencing technologies, allowing us to reconstruct a high-resolution transcriptome landscape of venom genes in the parasitoid wasp Pteromalus puparum Our findings demonstrate that more than half of the venom genes generate multiple isoforms within the venom gland. Through mass spectrometry analysis, we confirm that alternative splicing contributes to the diversity of venom proteins, acting as a mechanism for expanding the venom repertoire. Notably, we identified seven venom genes that exhibit distinct isoform usages between the venom gland and other tissues. Furthermore, evolutionary analyses of venom serpin3 and orcokinin further reveal that the co-option of an ancient isoform and a newly evolved isoform, respectively, contributes to venom recruitment, providing valuable insights into the genetic mechanisms driving venom evolution in parasitoid wasps. Overall, our study presents a comprehensive investigation of venom genes at the isoform level, significantly advancing our understanding of alternative isoforms in venom diversity and evolution and setting the stage for further in-depth research on venoms.
Collapse
Affiliation(s)
- Xinhai Ye
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China
- Shanghai Institute for Advanced Study, Zhejiang University, Shanghai 201203, China
- College of Computer Science and Technology, Zhejiang University, Hangzhou 310027, China
| | - Chun He
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China;
| | - Yi Yang
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China
| | - Yu H Sun
- Department of Biology, University of Rochester, Rochester, New York 14627, USA
| | - Shijiao Xiong
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China
| | - Kevin C Chan
- Shanghai Institute for Advanced Study, Zhejiang University, Shanghai 201203, China
| | - Yuxuan Si
- College of Computer Science and Technology, Zhejiang University, Hangzhou 310027, China
| | - Shan Xiao
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China
| | - Xianxin Zhao
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China
| | - Haiwei Lin
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China
| | - Yang Mei
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China
| | - Yufeng Yao
- Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Hangzhou 310020, China
| | - Gongyin Ye
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China
| | - Fei Wu
- Shanghai Institute for Advanced Study, Zhejiang University, Shanghai 201203, China;
- College of Computer Science and Technology, Zhejiang University, Hangzhou 310027, China
| | - Qi Fang
- State Key Laboratory of Rice Biology and Breeding & Ministry of Agricultural and Rural Affairs Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou 310058, China;
| |
Collapse
|
8
|
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
|
9
|
Gable SM, Mendez JM, Bushroe NA, Wilson A, Byars MI, Tollis M. The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order. Genes (Basel) 2023; 14:1387. [PMID: 37510292 PMCID: PMC10379679 DOI: 10.3390/genes14071387] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Revised: 06/28/2023] [Accepted: 06/29/2023] [Indexed: 07/30/2023] Open
Abstract
Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.
Collapse
Affiliation(s)
- Simone M Gable
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ 86011, USA
| | - Jasmine M Mendez
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ 86011, USA
| | - Nicholas A Bushroe
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ 86011, USA
| | - Adam Wilson
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ 86011, USA
| | - Michael I Byars
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ 86011, USA
| | - Marc Tollis
- School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ 86011, USA
| |
Collapse
|
10
|
Heptinstall TC, Strickland JL, Rosales-Garcia RA, Rautsaw RM, Simpson CL, Nystrom GS, Ellsworth SA, Hogan MP, Borja M, Fernandes Campos P, Grazziotin FG, Rokyta DR, Junqueira-de-Azevedo ILM, Parkinson CL. Venom phenotype conservation suggests integrated specialization in a lizard-eating snake. Toxicon 2023; 229:107135. [PMID: 37146732 PMCID: PMC11000244 DOI: 10.1016/j.toxicon.2023.107135] [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: 03/02/2023] [Revised: 04/18/2023] [Accepted: 04/24/2023] [Indexed: 05/07/2023]
Abstract
Biological specialization reduces the size of niche space while increasing efficiency in the use of available resources. Specialization often leads to phenotypic changes via natural selection aligning with niche space constraints. Commonly observed changes are in size, shape, behavior, and traits associated with feeding. One often selected trait for dietary specialization is venom, which, in snakes, often shows variation dependent on diet across and within species. The Neotropical Blunt-headed Treesnake (Imantodes cenchoa) is a highly specialized, rear-fanged, arboreal, lizard hunter that displays a long thin body, enlarged eyes, and a large Duvernoy's gland. However, toxin characterization of I. cenchoa has never been completed. Here, we use RNA-seq and mass spectrometry to assemble, annotate, and analyze the venom gland transcriptomes of four I. cenchoa from across their range. We find a lack of significant venom variation at the sequence and expression levels, suggesting venom conservation across the species. We propose this conservation provides evidence of a specialized venom repertoire, adapted to maximize efficiency of capturing and processing lizards. Importantly, this study provides the most complete venom gland transcriptomes of I. cenchoa and evidence of venom specialization in a rear-fanged snake, giving insight into selective pressures of venom across all snake species.
Collapse
Affiliation(s)
| | - 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; School of Biological Sciences, Washington State University, Pullman, WA, 99164, USA; Department of Integrative Biology, University of South Florida, Tampa, FL, 33620, USA
| | - Cassandra L Simpson
- Department of Biological Sciences, Clemson University, Clemson, SC, 29634, USA
| | - Gunnar S Nystrom
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Schyler A Ellsworth
- 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
| | - Miguel Borja
- Facultad de Ciencias Biológicas, Universdad Juárez del Estado de Durango, Av. Universidad s/n. Fracc. Filadelfia, C.P. 35070, Gómez Palacio, Dgo., Mexico
| | | | - Felipe G Grazziotin
- Laboratório Especial de Colecões Zoológicas, Instituto Butantan, São Paulo, São Paulo, Brazil
| | - 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
|
11
|
Westfall AK, Gopalan SS, Perry BW, Adams RH, Saviola AJ, Mackessy SP, Castoe TA. Single-Cell Heterogeneity in Snake Venom Expression Is Hardwired by Co-Option of Regulators from Progressively Activated Pathways. Genome Biol Evol 2023; 15:evad109. [PMID: 37311204 PMCID: PMC10289209 DOI: 10.1093/gbe/evad109] [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: 03/06/2023] [Revised: 05/31/2023] [Accepted: 06/07/2023] [Indexed: 06/15/2023] Open
Abstract
The ubiquitous cellular heterogeneity underlying many organism-level phenotypes raises questions about what factors drive this heterogeneity and how these complex heterogeneous systems evolve. Here, we use single-cell expression data from a Prairie rattlesnake (Crotalus viridis) venom gland to evaluate hypotheses for signaling networks underlying snake venom regulation and the degree to which different venom gene families have evolutionarily recruited distinct regulatory architectures. Our findings suggest that snake venom regulatory systems have evolutionarily co-opted trans-regulatory factors from extracellular signal-regulated kinase and unfolded protein response pathways that specifically coordinate expression of distinct venom toxins in a phased sequence across a single population of secretory cells. This pattern of co-option results in extensive cell-to-cell variation in venom gene expression, even between tandemly duplicated paralogs, suggesting this regulatory architecture has evolved to circumvent cellular constraints. While the exact nature of such constraints remains an open question, we propose that such regulatory heterogeneity may circumvent steric constraints on chromatin, cellular physiological constraints (e.g., endoplasmic reticulum stress or negative protein-protein interactions), or a combination of these. Regardless of the precise nature of these constraints, this example suggests that, in some cases, dynamic cellular constraints may impose previously unappreciated secondary constraints on the evolution of gene regulatory networks that favors heterogeneous expression.
Collapse
Affiliation(s)
| | | | - Blair W Perry
- Department of Biology, The University of Texas Arlington, Texas, USA
- School of Biological Sciences, Washington State University, Pullman, Washington, USA
| | - Richard H Adams
- Department of Entomology and Plant Pathology, University of Arkansas, Fayetteville, USA
| | - Anthony J Saviola
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver, Aurora, USA
| | - Stephen P Mackessy
- School of Biological Sciences, University of Northern Colorado, Greeley, USA
| | - Todd A Castoe
- Department of Biology, The University of Texas Arlington, Texas, USA
| |
Collapse
|
12
|
Pfenninger M, Foucault Q, Waldvogel AM, Feldmeyer B. Selective effects of a short transient environmental fluctuation on a natural population. Mol Ecol 2023; 32:335-349. [PMID: 36282585 DOI: 10.1111/mec.16748] [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: 02/15/2022] [Revised: 09/21/2022] [Accepted: 10/21/2022] [Indexed: 01/11/2023]
Abstract
Natural populations experience continuous and often transient changes of environmental conditions. These in turn may result in fluctuating selection pressures leading to variable demographic and evolutionary population responses. Rapid adaptation as short-term response to a sudden environmental change has in several cases been attributed to polygenic traits, but the underlying genomic dynamics and architecture are poorly understood. In this study, we took advantage of a natural experiment in an insect population of the non-biting midge Chironomus riparius by monitoring genome-wide allele frequencies before and after a cold snap event. Whole genome pooled sequencing of time series samples revealed 10 selected haplotypes carrying ancient polymorphisms, partially with signatures of balancing selection. By constantly cold exposing genetically variable individuals in the laboratory, we could demonstrate with whole genome resequencing (i) that among the survivors, the same alleles rose in frequency as in the wild, and (ii) that the identified variants additively predicted fitness (survival time) of its bearers. Finally, by simultaneously sequencing the genome and the transcriptome of cold exposed individuals we could tentatively link some of the selected SNPs to the cis- and trans-regulation of genes and pathways known to be involved in cold response of insects, such as cytochrome P450 and fatty acid metabolism. Altogether, our results shed light on the strength and speed of selection in natural populations and the genomic architecture of its underlying polygenic trait. Population genomic time series data thus appear as promising tool for measuring the selective tracking of fluctuating selection in natural populations.
Collapse
Affiliation(s)
- Markus Pfenninger
- Department of Molecular Ecology, Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, Germany.,LOEWE Centre for Translational Biodiversity Genomics, Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, Germany.,Institute for Molecular and Organismic Evolution, Johannes Gutenberg University, Mainz, Germany
| | - Quentin Foucault
- Department of Molecular Ecology, Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, Germany
| | - Ann-Marie Waldvogel
- Department of Ecological Genomics, Institute of Zoology, University of Cologne, Köln, Germany
| | - Barbara Feldmeyer
- Department of Molecular Ecology, Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, Germany
| |
Collapse
|
13
|
Pfenninger M, Foucault Q. Population Genomic Time Series Data of a Natural Population Suggests Adaptive Tracking of Fluctuating Environmental Changes. Integr Comp Biol 2022; 62:1812-1826. [PMID: 35762661 DOI: 10.1093/icb/icac098] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Revised: 06/07/2022] [Accepted: 06/16/2022] [Indexed: 01/05/2023] Open
Abstract
Natural populations are constantly exposed to fluctuating environmental changes that negatively affect their fitness in unpredictable ways. While theoretical models show the possibility of counteracting these environmental changes through rapid evolutionary adaptations, there have been few empirical studies demonstrating such adaptive tracking in natural populations. Here, we analyzed environmental data, fitness-related phenotyping and genomic time-series data sampled over 3 years from a natural Chironomus riparius (Diptera, Insecta) population to address this question. We show that the population's environment varied significantly on the time scale of the sampling in many selectively relevant dimensions, independently of each other. Similarly, phenotypic fitness components evolved significantly on the same temporal scale (mean 0.32 Haldanes), likewise independent from each other. The allele frequencies of 367,446 SNPs across the genome showed evidence of positive selection. Using temporal correlation of spatially coherent allele frequency changes revealed 35,574 haplotypes with more than one selected SNP. The mean selection coefficient for these haplotypes was 0.30 (s.d. = 0.68). The frequency changes of these haplotypes clustered in 46 different temporal patterns, indicating concerted, independent evolution of many polygenic traits. Nine of these patterns were strongly correlated with measured environmental variables. Enrichment analysis of affected genes suggested the implication of a wide variety of biological processes. Thus, our results suggest overall that the natural population of C. riparius tracks environmental change through rapid polygenic adaptation in many independent dimensions. This is further evidence that natural selection is pervasive at the genomic level and that evolutionary and ecological time scales may not differ at all, at least in some organisms.
Collapse
Affiliation(s)
- Markus Pfenninger
- Department Molecular Ecology, Senckenberg Biodiversity and Climate Research Centre, Senckenberganlage 25, 60325 Frankfurt am Main, Germany.,Institute for Molecular and Organismic Evolution, Johannes Gutenberg University, Johann-Joachim-Becher-Weg 7, 55128 Mainz, Germany.,LOEWE Centre for Translational Biodiversity Genomics, Senckenberg Biodiversity and Climate Research Centre, Senckenberganlage 25, 60325 Frankfurt am Main, Germany
| | - Quentin Foucault
- Department Molecular Ecology, Senckenberg Biodiversity and Climate Research Centre, Senckenberganlage 25, 60325 Frankfurt am Main, Germany.,Institute for Molecular and Organismic Evolution, Johannes Gutenberg University, Johann-Joachim-Becher-Weg 7, 55128 Mainz, Germany
| |
Collapse
|
14
|
Holding ML, Trevine VC, Zinenko O, Strickland JL, Rautsaw RM, Mason AJ, Hogan MP, Parkinson CL, Grazziotin FG, Santana SE, Davis MA, Rokyta DR. Evolutionary allometry and ecological correlates of fang length evolution in vipers. Proc Biol Sci 2022; 289:20221132. [PMID: 36300520 PMCID: PMC9449470 DOI: 10.1098/rspb.2022.1132] [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: 06/10/2022] [Accepted: 08/11/2022] [Indexed: 12/14/2022] Open
Abstract
Traits for prey acquisition form the phenotypic interface of predator-prey interactions. In venomous predators, morphological variation in venom delivery apparatus like fangs and stingers may be optimized for dispatching prey. Here, we determine how a single dimension of venom injection systems evolves in response to variation in the size, climatic conditions and dietary ecology of viperid snakes. We measured fang length in more than 1900 museum specimens representing 199 viper species (55% of recognized species). We find both phylogenetic signal and within-clade variation in relative fang length across vipers suggesting both general taxonomic trends and potential adaptive divergence in fang length. We recover positive evolutionary allometry and little static allometry in fang length. Proportionally longer fangs have evolved in larger species, which may facilitate venom injection in more voluminous prey. Finally, we leverage climatic and diet data to assess the global correlates of fang length. We find that models of fang length evolution are improved through the inclusion of both temperature and diet, particularly the extent to which diets are mammal-heavy diets. These findings demonstrate how adaptive variation can emerge among components of complex prey capture systems.
Collapse
Affiliation(s)
- Matthew L. Holding
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA
| | - Vivian C. Trevine
- Laboratório de Coleções Zoológicas, Instituto Butantan, Av. Vital Brazil, 1500, 05503-900 São Paulo, SP, Brazil
| | - Oleksandr Zinenko
- V. N. Karazin Kharkiv National University, 4 Svobody square, Kharkiv 61022, Ukraine
| | - Jason L. Strickland
- Department of Biological Sciences, Clemson University, 190 Collings street, Clemson, SC 29634, USA
| | - Rhett M. Rautsaw
- Department of Biological Sciences, Clemson University, 190 Collings street, Clemson, SC 29634, USA
| | - Andrew J. Mason
- Department of Biological Sciences, Clemson University, 190 Collings street, Clemson, SC 29634, USA
| | - Michael P. Hogan
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA
| | - Christopher L. Parkinson
- Department of Biological Sciences, Clemson University, 190 Collings street, Clemson, SC 29634, USA
| | - Felipe G. Grazziotin
- Laboratório de Coleções Zoológicas, Instituto Butantan, Av. Vital Brazil, 1500, 05503-900 São Paulo, SP, Brazil
| | - Sharlene E. Santana
- Department of Biology and Burke Museum of Natural History and Culture, University of Washington, Seattle, WA 98105, USA
| | - Mark A. Davis
- Illinois Natural History Survey, University of Illinois, Champaign, IL 61820, USA
| | - Darin R. Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA
| |
Collapse
|
15
|
Harrison CM, Colbert J, Richter CJ, McDonald PJ, Trumbull LM, Ellsworth SA, Hogan MP, Rokyta DR, Margres MJ. Using Morphological, Genetic, and Venom Analyses to Present Current and Historic Evidence of Crotalus horridus x adamanteus Hybridization on Jekyll Island, Georgia. SOUTHEAST NAT 2022. [DOI: 10.1656/058.021.0209] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Affiliation(s)
| | - Joseph Colbert
- Jekyll Island Authority Conservation Department, Jekyll Island, GA 31527
| | - Collin J. Richter
- Warnell School of Forestry and Natural Resources, University of Georgia, Athens, GA 30602
| | - Preston J. McDonald
- Department of Integrative Biology, University of South Florida, Tampa, FL 33620
| | - Lauren M. Trumbull
- Department of Integrative Biology, University of South Florida, Tampa, FL 33620
| | - Schyler A. Ellsworth
- Department of Biological Sciences, Florida State University, Tallahassee, FL 33306
| | - Michael P. Hogan
- Department of Biological Sciences, Florida State University, Tallahassee, FL 33306
| | - Darin R. Rokyta
- Department of Biological Sciences, Florida State University, Tallahassee, FL 33306
| | - Mark J. Margres
- Department of Integrative Biology, University of South Florida, Tampa, FL 33620
| |
Collapse
|
16
|
van Thiel J, Khan MA, Wouters RM, Harris RJ, Casewell NR, Fry BG, Kini RM, Mackessy SP, Vonk FJ, Wüster W, Richardson MK. Convergent evolution of toxin resistance in animals. Biol Rev Camb Philos Soc 2022; 97:1823-1843. [PMID: 35580905 PMCID: PMC9543476 DOI: 10.1111/brv.12865] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Revised: 04/22/2022] [Accepted: 04/26/2022] [Indexed: 12/17/2022]
Abstract
Convergence is the phenomenon whereby similar phenotypes evolve independently in different lineages. One example is resistance to toxins in animals. Toxins have evolved many times throughout the tree of life. They disrupt molecular and physiological pathways in target species, thereby incapacitating prey or deterring a predator. In response, molecular resistance has evolved in many species exposed to toxins to counteract their harmful effects. Here, we review current knowledge on the convergence of toxin resistance using examples from a wide range of toxin families. We explore the evolutionary processes and molecular adaptations driving toxin resistance. However, resistance adaptations may carry a fitness cost if they disrupt the normal physiology of the resistant animal. Therefore, there is a trade‐off between maintaining a functional molecular target and reducing toxin susceptibility. There are relatively few solutions that satisfy this trade‐off. As a result, we see a small set of molecular adaptations appearing repeatedly in diverse animal lineages, a phenomenon that is consistent with models of deterministic evolution. Convergence may also explain what has been called ‘autoresistance’. This is often thought to have evolved for self‐protection, but we argue instead that it may be a consequence of poisonous animals feeding on toxic prey. Toxin resistance provides a unique and compelling model system for studying the interplay between trophic interactions, selection pressures and the molecular mechanisms underlying evolutionary novelties.
Collapse
Affiliation(s)
- Jory van Thiel
- Institute of Biology Leiden, Leiden University, Sylviusweg 72, 2333 BE Leiden, The Netherlands
| | - Muzaffar A Khan
- Institute of Biology Leiden, Leiden University, Sylviusweg 72, 2333 BE Leiden, The Netherlands
| | - Roel M Wouters
- Institute of Biology Leiden, Leiden University, Sylviusweg 72, 2333 BE Leiden, The Netherlands
| | - Richard J Harris
- Venom Evolution Lab, School of Biological Sciences, University of Queensland, St Lucia, 4072, Australia
| | - Nicholas R Casewell
- Centre for Snakebite Research & Interventions, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool L3 5QA, U.K
| | - Bryan G Fry
- Venom Evolution Lab, School of Biological Sciences, University of Queensland, St Lucia, 4072, Australia
| | - R Manjunatha Kini
- Department of Biological Sciences, National University of Singapore, Singapore, 117558, Singapore.,Department of Pharmacology, Yong Loo Lin School of Medicine, National University of Singapore, Singapore, 117600, Singapore.,Department of Biochemistry, Medical College of Virginia, Virginia Commonwealth University, Richmond, VA, 23298, U.S.A
| | - Stephen P Mackessy
- School of Biological Sciences, University of Northern Colorado, Greeley, CO, 80639-0017, U.S.A
| | - Freek J Vonk
- Naturalis Biodiversity Center, Darwinweg 2, 2333 CR Leiden, The Netherlands.,Amsterdam Institute of Molecular and Life Sciences, Division of BioAnalytical Chemistry, Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1085, 1081HV Amsterdam, The Netherlands
| | - Wolfgang Wüster
- Molecular Ecology and Fisheries Genetics Laboratory, School of Natural Sciences, Bangor University, Bangor, LL57 2UW, U.K
| | - Michael K Richardson
- Institute of Biology Leiden, Leiden University, Sylviusweg 72, 2333 BE Leiden, The Netherlands
| |
Collapse
|
17
|
Roldán-Padrón O, Cruz-Pérez MS, Castro-Guillén JL, García-Arredondo JA, Mendiola-Olaya E, Saldaña-Gutiérrez C, Herrera-Paniagua P, Blanco-Labra A, García-Gasca T. Hybridization between Crotalus aquilus and Crotalus polystictus Species: A Comparison of Their Venom Toxicity and Enzymatic Activities. BIOLOGY 2022; 11:661. [PMID: 35625389 PMCID: PMC9138290 DOI: 10.3390/biology11050661] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/26/2022] [Revised: 04/15/2022] [Accepted: 04/17/2022] [Indexed: 06/15/2023]
Abstract
Hybridization is defined as the interbreeding of individuals from two populations distinguishable by one or more heritable characteristics. Snake hybridization represents an interesting opportunity to analyze variability and how genetics affect the venom components between parents and hybrids. Snake venoms exhibit a high degree of variability related to biological and biogeographical factors. The aim of this work is to analyze the protein patterns and enzymatic activity of some of the main hemotoxic enzymes in snake venoms, such as serine proteases (trypsin-like, chymotrypsin-like, and elastase-like), metalloproteases, hyaluronidases, and phospholipase A2. The lethal dose of 50 (LD50) of venom from the Crotalus aquilus (Cabf) and Crotalus polystictus (Cpbm) parents and their hybrids in captivity was determined, and phenetic analysis is also conducted, which showed a high similarity between the hybrids and C. polystictus. The protein banding patterns and enzymatic activity analyze by zymography resulted in a combination of proteins from the parental venoms in the hybrids, with variability among them. In some cases, the enzymatic activity is higher in the hybrids with a lower LD50 than in the parents, indicating higher toxicity. These data show the variability among snake venoms and suggest that hybridization is an important factor in changes in protein concentration, peptide variability, and enzymatic activity that affect toxicity and lethality.
Collapse
Affiliation(s)
- Octavio Roldán-Padrón
- Laboratorio de Biología Celular y Molecuar, Facultad de Ciencias Naturales, Universidad Autónoma de Querétaro, Av. de las Ciencias S/N, Juriquilla, Queretaro 76230, Qro, Mexico; (O.R.-P.); (M.S.C.-P.); (C.S.-G.); (P.H.-P.)
| | - Martha Sandra Cruz-Pérez
- Laboratorio de Biología Celular y Molecuar, Facultad de Ciencias Naturales, Universidad Autónoma de Querétaro, Av. de las Ciencias S/N, Juriquilla, Queretaro 76230, Qro, Mexico; (O.R.-P.); (M.S.C.-P.); (C.S.-G.); (P.H.-P.)
| | - José Luis Castro-Guillén
- Departamento de Biotecnología y Bioquímica, Centro de Investigación y Estudios Avanzados del IPN, Km. 9.6 Libramiento Norte Carr. Irapuato-León, Irapuato 36824, Gto, Mexico; (J.L.C.-G.); (E.M.-O.)
| | - José Alejandro García-Arredondo
- Laboratorio de Química Medicinal, Facultad de Química, Universidad Autónoma de Querétaro, Cerro de las Campanas S/N, Centro Universitario, Queretaro 76010, Qro, Mexico;
| | - Elizabeth Mendiola-Olaya
- Departamento de Biotecnología y Bioquímica, Centro de Investigación y Estudios Avanzados del IPN, Km. 9.6 Libramiento Norte Carr. Irapuato-León, Irapuato 36824, Gto, Mexico; (J.L.C.-G.); (E.M.-O.)
| | - Carlos Saldaña-Gutiérrez
- Laboratorio de Biología Celular y Molecuar, Facultad de Ciencias Naturales, Universidad Autónoma de Querétaro, Av. de las Ciencias S/N, Juriquilla, Queretaro 76230, Qro, Mexico; (O.R.-P.); (M.S.C.-P.); (C.S.-G.); (P.H.-P.)
| | - Patricia Herrera-Paniagua
- Laboratorio de Biología Celular y Molecuar, Facultad de Ciencias Naturales, Universidad Autónoma de Querétaro, Av. de las Ciencias S/N, Juriquilla, Queretaro 76230, Qro, Mexico; (O.R.-P.); (M.S.C.-P.); (C.S.-G.); (P.H.-P.)
| | - Alejandro Blanco-Labra
- Departamento de Biotecnología y Bioquímica, Centro de Investigación y Estudios Avanzados del IPN, Km. 9.6 Libramiento Norte Carr. Irapuato-León, Irapuato 36824, Gto, Mexico; (J.L.C.-G.); (E.M.-O.)
| | - Teresa García-Gasca
- Laboratorio de Biología Celular y Molecuar, Facultad de Ciencias Naturales, Universidad Autónoma de Querétaro, Av. de las Ciencias S/N, Juriquilla, Queretaro 76230, Qro, Mexico; (O.R.-P.); (M.S.C.-P.); (C.S.-G.); (P.H.-P.)
| |
Collapse
|
18
|
Mason AJ, Holding ML, Rautsaw RM, Rokyta DR, Parkinson CL, Gibbs HL. Venom gene sequence diversity and expression jointly shape diet adaptation in pitvipers. Mol Biol Evol 2022; 39:6567549. [PMID: 35413123 PMCID: PMC9040050 DOI: 10.1093/molbev/msac082] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA2s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms—such as higher absolute levels of expression—are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey.
Collapse
Affiliation(s)
- Andrew J Mason
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University, Columbus, OH, USA
| | | | - Rhett M Rautsaw
- Department of Biological Sciences, Clemson University, Clemson, SC, USA
| | - Darin R Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL, USA
| | - Christopher L Parkinson
- Department of Biological Sciences, Clemson University, Clemson, SC, USA.,Department of Forestry and Environmental Conservation, Clemson University, Clemson, SC, USA
| | - H Lisle Gibbs
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University, Columbus, OH, USA
| |
Collapse
|
19
|
Brown KE, Kelly JK. Genome-wide association mapping of transcriptome variation in Mimulus guttatus indicates differing patterns of selection on cis- versus trans-acting mutations. Genetics 2022; 220:iyab189. [PMID: 34791192 PMCID: PMC8733635 DOI: 10.1093/genetics/iyab189] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Accepted: 09/28/2021] [Indexed: 11/14/2022] Open
Abstract
We measured the floral bud transcriptome of 151 fully sequenced lines of Mimulus guttatus from one natural population. Thousands of single nucleotide polymorphisms (SNPs) are implicated as transcription regulators, but there is a striking difference in the allele frequency spectrum of cis-acting and trans-acting mutations. Cis-SNPs have intermediate frequencies (consistent with balancing selection) while trans-SNPs exhibit a rare-alleles model (consistent with purifying selection). This pattern only becomes clear when transcript variation is normalized on a gene-to-gene basis. If a global normalization is applied, as is typically in RNAseq experiments, asymmetric transcript distributions combined with "rarity disequilibrium" produce a superabundance of false positives for trans-acting SNPs. To explore the cause of purifying selection on trans-acting mutations, we identified gene expression modules as sets of coexpressed genes. The extent to which trans-acting mutations influence modules is a strong predictor of allele frequency. Mutations altering expression of genes with high "connectedness" (those that are highly predictive of the representative module expression value) have the lowest allele frequency. The expression modules can also predict whole-plant traits such as flower size. We find that a substantial portion of the genetic (co)variance among traits can be described as an emergent property of genetic effects on expression modules.
Collapse
Affiliation(s)
- Keely E Brown
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA
| | - John K Kelly
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA
| |
Collapse
|
20
|
Barua A, Koludarov I, Mikheyev AS. Co-option of the same ancestral gene family gave rise to mammalian and reptilian toxins. BMC Biol 2021; 19:268. [PMID: 34949191 PMCID: PMC8705180 DOI: 10.1186/s12915-021-01191-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2021] [Accepted: 11/11/2021] [Indexed: 12/03/2022] Open
Abstract
Background Evolution can occur with surprising predictability when organisms face similar ecological challenges. For most traits, it is difficult to ascertain whether this occurs due to constraints imposed by the number of possible phenotypic solutions or because of parallel responses by shared genetic and regulatory architecture. Exceptionally, oral venoms are a tractable model of trait evolution, being largely composed of proteinaceous toxins that have evolved in many tetrapods, ranging from reptiles to mammals. Given the diversity of venomous lineages, they are believed to have evolved convergently, even though biochemically similar toxins occur in all taxa. Results Here, we investigate whether ancestral genes harbouring similar biochemical activity may have primed venom evolution, focusing on the origins of kallikrein-like serine proteases that form the core of most vertebrate oral venoms. Using syntenic relationships between genes flanking known toxins, we traced the origin of kallikreins to a single locus containing one or more nearby paralogous kallikrein-like clusters. Additionally, phylogenetic analysis of vertebrate serine proteases revealed that kallikrein-like toxins in mammals and reptiles are genetically distinct from non-toxin ones. Conclusions Given the shared regulatory and genetic machinery, these findings suggest that tetrapod venoms evolved by co-option of proteins that were likely already present in saliva. We term such genes ‘toxipotent’—in the case of salivary kallikreins they already had potent vasodilatory activity that was weaponized by venomous lineages. Furthermore, the ubiquitous distribution of kallikreins across vertebrates suggests that the evolution of envenomation may be more common than previously recognized, blurring the line between venomous and non-venomous animals. Supplementary Information The online version contains supplementary material available at 10.1186/s12915-021-01191-1.
Collapse
Affiliation(s)
- Agneesh Barua
- Ecology and Evolution Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa, Japan.
| | - Ivan Koludarov
- Animal Venomics Group, Justus Leibig University, Giessen, Germany
| | - Alexander S Mikheyev
- Research School of Biology, Australian National University, Canberra, ACT, Australia.
| |
Collapse
|
21
|
Varying Intensities of Introgression Obscure Incipient Venom-Associated Speciation in the Timber Rattlesnake ( Crotalus horridus). Toxins (Basel) 2021; 13:toxins13110782. [PMID: 34822565 PMCID: PMC8625053 DOI: 10.3390/toxins13110782] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2021] [Revised: 10/14/2021] [Accepted: 10/25/2021] [Indexed: 11/16/2022] Open
Abstract
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.
Collapse
|
22
|
Robinson KE, Holding ML, Whitford MD, Saviola AJ, Yates JR, Clark RW. Phenotypic and functional variation in venom and venom resistance of two sympatric rattlesnakes and their prey. J Evol Biol 2021; 34:1447-1465. [PMID: 34322920 DOI: 10.1111/jeb.13907] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Revised: 06/27/2021] [Accepted: 07/13/2021] [Indexed: 12/15/2022]
Abstract
Predator-prey interactions often lead to the coevolution of adaptations associated with avoiding predation and, for predators, overcoming those defences. Antagonistic coevolutionary relationships are often not simple interactions between a single predator and prey but rather a complex web of interactions between multiple coexisting species. Coevolution between venomous rattlesnakes and small mammals has led to physiological venom resistance in several mammalian taxa. In general, viperid venoms contain large quantities of snake venom metalloproteinase toxins (SVMPs), which are inactivated by SVMP inhibitors expressed in resistant mammals. We explored variation in venom chemistry, SVMP expression, and SVMP resistance across four co-distributed species (California Ground Squirrels, Bryant's Woodrats, Southern Pacific Rattlesnakes, and Red Diamond Rattlesnakes) collected from four different populations in Southern California. Our aim was to understand phenotypic and functional variation in venom and venom resistance in order to compare coevolutionary dynamics of a system involving two sympatric predator-prey pairs to past studies that have focused on single pairs. Proteomic analysis of venoms indicated that these rattlesnakes express different phenotypes when in sympatry, with Red Diamonds expressing more typical viperid venom (with a diversity of SVMPs) and Southern Pacifics expressing a more atypical venom with a broader range of non-enzymatic toxins. We also found that although blood sera from both mammals were generally able to inhibit SVMPs from both rattlesnake species, inhibition depended strongly on the snake population, with snakes from one geographic site expressing SVMPs to which few mammals were resistant. Additionally, we found that Red Diamond venom, rather than woodrat resistance, was locally adapted. Our findings highlight the complexity of coevolutionary relationships between multiple predators and prey that exhibit similar offensive and defensive strategies in sympatry.
Collapse
Affiliation(s)
- Kelly E Robinson
- Department of Biology, San Diego State University, San Diego, CA, USA.,Department of Biology, University of Nevada, Reno, NV, USA.,Program in Ecology, Evolution and Conservation Biology, University of Nevada, Reno, NV, USA
| | - Matthew L Holding
- Department of Natural Resources and Environmental Science, University of Nevada, Reno, NV, USA.,Department of Biological Science, Florida State University, Tallahassee, FL, USA
| | - Malachi D Whitford
- Department of Biology, San Diego State University, San Diego, CA, USA.,Ecology Graduate Group, University of California, Davis, CA, USA
| | - Anthony J Saviola
- Department of Molecular Medicine and Neurobiology, The Scripps Research Institute, La Jolla, CA, USA.,Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
| | - John R Yates
- Department of Molecular Medicine and Neurobiology, The Scripps Research Institute, La Jolla, CA, USA
| | - Rulon W Clark
- Department of Biology, San Diego State University, San Diego, CA, USA
| |
Collapse
|
23
|
Calvete JJ, Lomonte B, Saviola AJ, Bonilla F, Sasa M, Williams DJ, Undheim EA, Sunagar K, Jackson TN. Mutual enlightenment: A toolbox of concepts and methods for integrating evolutionary and clinical toxinology via snake venomics and the contextual stance. Toxicon X 2021; 9-10:100070. [PMID: 34195606 PMCID: PMC8234350 DOI: 10.1016/j.toxcx.2021.100070] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Revised: 06/01/2021] [Accepted: 06/07/2021] [Indexed: 12/21/2022] Open
Abstract
Snakebite envenoming is a neglected tropical disease that may claim over 100,000 human lives annually worldwide. Snakebite occurs as the result of an interaction between a human and a snake that elicits either a defensive response from the snake or, more rarely, a feeding response as the result of mistaken identity. Snakebite envenoming is therefore a biological and, more specifically, an ecological problem. Snake venom itself is often described as a "cocktail", as it is a heterogenous mixture of molecules including the toxins (which are typically proteinaceous) responsible for the pathophysiological consequences of envenoming. The primary function of venom in snake ecology is pre-subjugation, with defensive deployment of the secretion typically considered a secondary function. The particular composition of any given venom cocktail is shaped by evolutionary forces that include phylogenetic constraints associated with the snake's lineage and adaptive responses to the snake's ecological context, including the taxa it preys upon and by which it is predated upon. In the present article, we describe how conceptual frameworks from ecology and evolutionary biology can enter into a mutually enlightening relationship with clinical toxinology by enabling the consideration of snakebite envenoming from an "ecological stance". We detail the insights that may emerge from such a perspective and highlight the ways in which the high-fidelity descriptive knowledge emerging from applications of -omics era technologies - "venomics" and "antivenomics" - can combine with evolutionary explanations to deliver a detailed understanding of this multifactorial health crisis.
Collapse
Affiliation(s)
- Juan J. Calvete
- Evolutionary and Translational Venomics Laboratory, Instituto de Biomedicina de Valencia, CSIC, Valencia, Spain
| | - Bruno Lomonte
- Unidad de Proteómica, Instituto Clodomiro Picado, Facultad de Microbiología, Universidad de Costa Rica, San José, Costa Rica
| | - Anthony J. Saviola
- Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
| | - Fabián Bonilla
- Laboratorio de Investigación en Animales Peligrosos (LIAP), Instituto Clodomiro Picado, Facultad de Microbiología, Universidad de Costa Rica, San José, Costa Rica
| | - Mahmood Sasa
- Laboratorio de Investigación en Animales Peligrosos (LIAP), Instituto Clodomiro Picado, Facultad de Microbiología, Universidad de Costa Rica, San José, Costa Rica
- Museo de Zoología, Centro de Investigaciones en Biodiversidad y Ecología Tropical, Universidad de Costa Rica, Costa Rica
| | | | - Eivind A.B. Undheim
- Centre for Biodiversity Dynamics, Department of Biology, NTNU, Trondheim, Norway
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Kartik Sunagar
- Evolutionary Venomics Lab, Centre for Ecological Sciences, Indian Institute of Science, Bangalore, Karnataka, India
| | - Timothy N.W. Jackson
- Australian Venom Research Unit, Department of Pharmacology and Therapeutics, University of Melbourne, Melbourne, Australia
| |
Collapse
|
24
|
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
|
25
|
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: 44] [Impact Index Per Article: 14.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
|
26
|
The Chemosensory Repertoire of the Eastern Diamondback Rattlesnake (Crotalus adamanteus) Reveals Complementary Genetics of Olfactory and Vomeronasal-Type Receptors. J Mol Evol 2021; 89:313-328. [PMID: 33881604 DOI: 10.1007/s00239-021-10007-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Accepted: 03/22/2021] [Indexed: 01/14/2023]
Abstract
Pitviper sensory perception incorporates diverse stimuli through the integration of trichromatic color vision, bifocal heat-sensing, and dual-system chemoperception. Chemoperception, or olfaction, is mediated by chemoreceptors in the olfactory bulb and the vomeronasal organ, but the true genomic complexity of the gene families and their relative contributions is unknown. A full genomic accounting of pitviper chemoperception directly complements our current understanding of their venoms by generating a more complete polyphenic representation of their predatory arsenal. To characterize the genetic repertoire of pitviper chemoperception, we analyzed a full-genome assembly for Crotalus adamanteus, the eastern diamondback rattlesnake. We identified hundreds of genes encoding both olfactory receptors (ORs; 362 full-length genes) and type-2 vomeronasal receptors (V2Rs; 430 full-length genes). Many chemoreceptor genes are organized into large tandem repeat arrays. Comparative analysis of V2R orthologs across squamates demonstrates how gene array expansion and contraction underlies the evolution of the chemoreceptor repertoire, which likely reflects shifts in life history traits. Chromosomal assignments of chemosensory genes identified sex chromosome specific chemoreceptor genes, providing gene candidates underlying observed sex-specific chemosensory-based behaviors. We detected widespread episodic evolution in the extracellular, ligand-binding domains of both ORs and V2Rs, suggesting the diversification of chemoreceptors is driven by transient periods of positive selection. We provide a robust genetic framework for studying pitviper chemosensory ecology and evolution.
Collapse
|
27
|
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
|
28
|
Margres MJ, Rautsaw RM, Strickland JL, Mason AJ, Schramer TD, Hofmann EP, Stiers E, Ellsworth SA, Nystrom GS, Hogan MP, Bartlett DA, Colston TJ, Gilbert DM, Rokyta DR, Parkinson CL. The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype. Proc Natl Acad Sci U S A 2021; 118:e2014634118. [PMID: 33468678 PMCID: PMC7848695 DOI: 10.1073/pnas.2014634118] [Citation(s) in RCA: 41] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Mark J Margres
- Department of Biological Sciences, Clemson University, Clemson, SC 29634;
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138
- Department of Integrative Biology, University of South Florida, Tampa, FL 33620
| | - Rhett M Rautsaw
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
| | - Jason L Strickland
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
- Department of Biology, University of South Alabama, Mobile, AL 36688
| | - Andrew J Mason
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
| | - Tristan D Schramer
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
| | - Erich P Hofmann
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
| | - Erin Stiers
- Department of Biological Sciences, Clemson University, Clemson, SC 29634
| | - Schyler A Ellsworth
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Gunnar S Nystrom
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Michael P Hogan
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Daniel A Bartlett
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - Timothy J Colston
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - David M Gilbert
- Department of Biological Science, Florida State University, Tallahassee, FL 32306
| | - 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
|
29
|
Banker SE, Lemmon AR, Hassinger AB, Dye M, Holland SD, Kortyna ML, Ospina OE, Ralicki H, Lemmon EM. Hierarchical Hybrid Enrichment: Multitiered Genomic Data Collection Across Evolutionary Scales, With Application to Chorus Frogs (Pseudacris). Syst Biol 2021; 69:756-773. [PMID: 31886503 PMCID: PMC7302053 DOI: 10.1093/sysbio/syz074] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2018] [Revised: 10/26/2019] [Accepted: 11/04/2019] [Indexed: 01/25/2023] Open
Abstract
Determining the optimal targets of genomic subsampling for phylogenomics, phylogeography, and population genomics remains a challenge for evolutionary biologists. Of the available methods for subsampling the genome, hybrid enrichment (sequence capture) has become one of the primary means of data collection for systematics, due to the flexibility and cost efficiency of this approach. Despite the utility of this method, information is lacking as to what genomic targets are most appropriate for addressing questions at different evolutionary scales. In this study, first, we compare the benefits of target loci developed for deep- and shallow scales by comparing these loci at each of three taxonomic levels: within a genus (phylogenetics), within a species (phylogeography), and within a hybrid zone (population genomics). Specifically, we target evolutionarily conserved loci that are appropriate for deeper phylogenetic scales and more rapidly evolving loci that are informative for phylogeographic and population genomic scales. Second, we assess the efficacy of targeting multiple-locus sets for different taxonomic levels in the same hybrid enrichment reaction, an approach we term hierarchical hybrid enrichment. Third, we apply this approach to the North American chorus frog genus Pseudacris to answer key evolutionary questions across taxonomic and temporal scales. We demonstrate that in this system the type of genomic target that produces the most resolved gene trees differs depending on the taxonomic level, although the potential for error is substantially lower for the deep-scale loci at all levels. We successfully recover data for the two different locus sets with high efficiency. Using hierarchical data targeting deep and shallow levels: we 1) resolve the phylogeny of the genus Pseudacris and introduce a novel visual and hypothesis testing method that uses nodal heat maps to examine the robustness of branch support values to the removal of sites and loci; 2) estimate the phylogeographic history of Pseudacris feriarum, which reveals up to five independent invasions leading to sympatry with congener Pseudacris nigrita to form replicated reinforcement contact zones with ongoing gene flow into sympatry; and 3) quantify with high confidence the frequency of hybridization in one of these zones between P. feriarum and P. nigrita, which is lower than microsatellite-based estimates. We find that the hierarchical hybrid enrichment approach offers an efficient, multitiered data collection method for simultaneously addressing questions spanning multiple evolutionary scales. [Anchored hybrid enrichment; heat map; hybridization; phylogenetics; phylogeography; population genomics; reinforcement; reproductive character displacement.].
Collapse
Affiliation(s)
- Sarah E Banker
- Department of Biological Science, Florida State University, 319 Stadium Drive, Tallahassee, FL 32306, USA.,Department of Integrative Biology, University of California, Berkeley, #3160 Berkeley, CA 94720-3160, USA
| | - Alan R Lemmon
- Department of Scientific Computing, Florida State University 400 Dirac Science Library, Tallahassee, FL 32306, USA
| | - Alyssa Bigelow Hassinger
- Department of Biological Science, Florida State University, 319 Stadium Drive, Tallahassee, FL 32306, USA.,Department of Evolution, Ecology, and Organismal Biology, Ohio State University, 318 W. 12th Avenue, 300 Aronoff Laboratory, Columbus, OH 43210, USA
| | - Mysia Dye
- Department of Biological Science, Florida State University, 319 Stadium Drive, Tallahassee, FL 32306, USA
| | - Sean D Holland
- Department of Biological Science, Florida State University, 319 Stadium Drive, Tallahassee, FL 32306, USA
| | - Michelle L Kortyna
- Department of Biological Science, Florida State University, 319 Stadium Drive, Tallahassee, FL 32306, USA
| | - Oscar E Ospina
- Department of Biological Science, Florida State University, 319 Stadium Drive, Tallahassee, FL 32306, USA
| | - Hannah Ralicki
- Department of Biological Science, Florida State University, 319 Stadium Drive, Tallahassee, FL 32306, USA.,Department of Biological Science, University of Connecticut, 91 North Eagleville Road, Storrs, CT 06268, USA
| | - Emily Moriarty Lemmon
- Department of Biological Science, Florida State University, 319 Stadium Drive, Tallahassee, FL 32306, USA
| |
Collapse
|
30
|
Hamanaka K, Mori A. Toxicity of venom from the mamushi, Gloydius blomhoffii, (Squamata, Crotalinae) to centipedes. Toxicon 2020; 188:11-15. [PMID: 33039367 DOI: 10.1016/j.toxicon.2020.10.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2020] [Revised: 09/28/2020] [Accepted: 10/05/2020] [Indexed: 02/06/2023]
Abstract
Vipers include approximately 300 species and usually feed on vertebrates, but over 30 species of them occasionally eat centipedes. Centipedes have been also known to occur in stomach contents of a Japanese pit viper, mamushi, Gloydius blomhoffii. Toxicity of the venom of mamushi to small mammals has been well studied, but there is no information concerning its toxicity to arthropods. Here, we studied the toxicity of the raw venom to the red-headed centipede, Scolopendra subspinipes mutilans, by comparing with the toxicity to two other common prey animals, a house mouse, Mus musclus, and a pond frog, Pelophylax nigromaculatus. The lethal doses for mice weighing around 21.5 g and frogs weighing around 3.78 g were less than 5 μl (equivalent to ca. 0.23 and 1.32 μl/g, respectively), which presumably corresponds to an approximate dose of mamushi's one envenomation. On the other hand, centipedes weighing around 1.86 g needed 10-36 μl of venom to die (16.0 μl/g on average). This result suggests that the centipedes are much more resistant to the venom than other prey animals, and it is difficult for mamushi to kill or incapacitate centipedes by the venom of a single envenomation.
Collapse
Affiliation(s)
- Kyosuke Hamanaka
- Department of Zoology, Graduate School of Science, Kyoto University, Sakyo, Kyoto, 606-8502, Japan.
| | - Akira Mori
- Department of Zoology, Graduate School of Science, Kyoto University, Sakyo, Kyoto, 606-8502, Japan
| |
Collapse
|
31
|
Temperature-dependent life history and transcriptomic responses in heat-tolerant versus heat-sensitive Brachionus rotifers. Sci Rep 2020; 10:13281. [PMID: 32764662 PMCID: PMC7411042 DOI: 10.1038/s41598-020-70173-0] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2020] [Accepted: 07/22/2020] [Indexed: 01/12/2023] Open
Abstract
Thermal stress response is an essential physiological trait that determines occurrence and temporal succession in nature, including response to climate change. We compared temperature-related demography in closely related heat-tolerant and heat-sensitive Brachionus rotifer species. We found significant differences in heat response, with the heat-sensitive species adopting a strategy of long survival and low population growth, while the heat-tolerant followed the opposite strategy. In both species, we examined the genetic basis of physiological variation by comparing gene expression across increasing temperatures. Comparative transcriptomic analyses identified shared and opposing responses to heat. Interestingly, expression of heat shock proteins (hsps) was strikingly different in the two species and mirrored differences in population growth rates, showing that hsp genes are likely a key component of a species' adaptation to different temperatures. Temperature induction caused opposing patterns of expression in further functional categories including energy, carbohydrate and lipid metabolism, and in genes related to ribosomal proteins. In the heat-sensitive species, elevated temperatures caused up-regulation of genes related to meiosis induction and post-translational histone modifications. This work demonstrates the sweeping reorganizations of biological functions that accompany temperature adaptation in these two species and reveals potential molecular mechanisms that might be activated for adaptation to global warming.
Collapse
|
32
|
Ochoa A, Broe M, Moriarty Lemmon E, Lemmon AR, Rokyta DR, Gibbs HL. Drift, selection and adaptive variation in small populations of a threatened rattlesnake. Mol Ecol 2020; 29:2612-2625. [PMID: 32557885 DOI: 10.1111/mec.15517] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2019] [Revised: 05/09/2020] [Accepted: 05/21/2020] [Indexed: 01/22/2023]
Abstract
An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long- and short-term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre-date the effects of recent drift, and that functional variation in these loci persists despite small short-term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of "drift debt," a nonequilibrium state where present-day levels of variation overestimate the amount of functional genetic diversity present in future populations.
Collapse
Affiliation(s)
- Alexander Ochoa
- Ohio Biodiversity Conservation Partnership and Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, USA
| | - Michael Broe
- Ohio Biodiversity Conservation Partnership and Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, USA
| | | | - Alan R Lemmon
- Department of Scientific Computing, Florida State University, Tallahassee, FL, USA
| | - Darin R Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL, USA
| | - H Lisle Gibbs
- Ohio Biodiversity Conservation Partnership and Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, USA
| |
Collapse
|
33
|
Spontaneous Tumor Regression in Tasmanian Devils Associated with RASL11A Activation. Genetics 2020; 215:1143-1152. [PMID: 32554701 DOI: 10.1534/genetics.120.303428] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2020] [Accepted: 06/12/2020] [Indexed: 12/30/2022] Open
Abstract
Spontaneous tumor regression has been documented in a small proportion of human cancer patients, but the specific mechanisms underlying tumor regression without treatment are not well understood. Tasmanian devils are threatened with extinction from a transmissible cancer due to universal susceptibility and a near 100% case fatality rate. In over 10,000 cases, <20 instances of natural tumor regression have been detected. Previous work in this system has focused on Tasmanian devil genetic variation associated with the regression phenotype. Here, we used comparative and functional genomics to identify tumor genetic variation associated with tumor regression. We show that a single point mutation in the 5' untranslated region of the putative tumor suppressor RASL11A significantly contributes to tumor regression. RASL11A was expressed in regressed tumors but silenced in wild-type, nonregressed tumors, consistent with RASL11A downregulation in human cancers. Induced RASL11A expression significantly reduced tumor cell proliferation in vitro The RAS pathway is frequently altered in human cancers, and RASL11A activation may provide a therapeutic treatment option for Tasmanian devils as well as a general mechanism for tumor inhibition.
Collapse
|
34
|
Barua A, Mikheyev AS. Toxin expression in snake venom evolves rapidly with constant shifts in evolutionary rates. Proc Biol Sci 2020; 287:20200613. [PMID: 32345154 PMCID: PMC7282918 DOI: 10.1098/rspb.2020.0613] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2020] [Accepted: 03/30/2020] [Indexed: 12/21/2022] Open
Abstract
Key innovations provide ecological opportunity by enabling access to new resources, colonization of new environments, and are associated with adaptive radiation. The most well-known pattern associated with adaptive radiation is an early burst of phenotypic diversification. Venoms facilitate prey capture and are widely believed to be key innovations leading to adaptive radiation. However, few studies have estimated their evolutionary rate dynamics. Here, we test for patterns of adaptive evolution in venom gene expression data from 52 venomous snake species. By identifying shifts in tempo and mode of evolution along with models of phenotypic evolution, we show that snake venom exhibits the macroevolutionary dynamics expected of key innovations. Namely, all toxin families undergo shifts in their rates of evolution, likely in response to changes in adaptive optima. Furthermore, we show that rapid-pulsed evolution modelled as a Lévy process better fits snake venom evolution than conventional early burst or Ornstein-Uhlenbeck models. While our results support the idea of snake venom being a key innovation, the innovation of venom chemistry lacks clear mechanisms that would lead to reproductive isolation and thus adaptive radiation. Therefore, the extent to which venom directly influences the diversification process is still a matter of contention.
Collapse
Affiliation(s)
- Agneesh Barua
- Ecology and Evolution Unit, Okinawa Institute of Science and Technology Graduate University, 1919-1 Tancha, Onna-son, Kunigami-gun, Okinawa-ken 904-0495, Japan
| | - Alexander S. Mikheyev
- Ecology and Evolution Unit, Okinawa Institute of Science and Technology Graduate University, 1919-1 Tancha, Onna-son, Kunigami-gun, Okinawa-ken 904-0495, Japan
- Evolutionary genomics group, Australian National University, Canberra ACT 0200, Australia
| |
Collapse
|
35
|
Ashwood LM, Norton RS, Undheim EAB, Hurwood DA, Prentis PJ. Characterising Functional Venom Profiles of Anthozoans and Medusozoans within Their Ecological Context. Mar Drugs 2020; 18:E202. [PMID: 32283847 PMCID: PMC7230708 DOI: 10.3390/md18040202] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2020] [Revised: 03/31/2020] [Accepted: 04/06/2020] [Indexed: 12/13/2022] Open
Abstract
This review examines the current state of knowledge regarding toxins from anthozoans (sea anemones, coral, zoanthids, corallimorphs, sea pens and tube anemones). We provide an overview of venom from phylum Cnidaria and review the diversity of venom composition between the two major clades (Medusozoa and Anthozoa). We highlight that the functional and ecological context of venom has implications for the temporal and spatial expression of protein and peptide toxins within class Anthozoa. Understanding the nuances in the regulation of venom arsenals has been made possible by recent advances in analytical technologies that allow characterisation of the spatial distributions of toxins. Furthermore, anthozoans are unique in that ecological roles can be assigned using tissue expression data, thereby circumventing some of the challenges related to pharmacological screening.
Collapse
Affiliation(s)
- Lauren M. Ashwood
- School of Biology and Environmental Science, Science and Engineering Faculty, Queensland University of Technology, Brisbane, QLD 4000, Australia
| | - Raymond S. Norton
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, 381 Royal Parade, Parkville, Victoria 3052, Australia
- ARC Centre for Fragment-Based Design, Monash University, Parkville, Victoria 3052, Australia
| | - Eivind A. B. Undheim
- Centre for Biodiversity Dynamics, Department of Biology, Norwegian University of Science and Technology, 7491 Trondheim, Norway
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, PO Box 1066 Blindern, 0316 Oslo, Norway
- Centre for Advanced Imaging, University of Queensland, St Lucia, QLD 4072, Australia
| | - David A. Hurwood
- School of Biology and Environmental Science, Science and Engineering Faculty, Queensland University of Technology, Brisbane, QLD 4000, Australia
- Institute of Future Environments, Queensland University of Technology, Brisbane, QLD 4000, Australia
| | - Peter J. Prentis
- School of Biology and Environmental Science, Science and Engineering Faculty, Queensland University of Technology, Brisbane, QLD 4000, Australia
- Institute of Future Environments, Queensland University of Technology, Brisbane, QLD 4000, Australia
| |
Collapse
|
36
|
Zancolli G, Calvete JJ, Cardwell MD, Greene HW, Hayes WK, Hegarty MJ, Herrmann HW, Holycross AT, Lannutti DI, Mulley JF, Sanz L, Travis ZD, Whorley JR, Wüster CE, Wüster W. When one phenotype is not enough: divergent evolutionary trajectories govern venom variation in a widespread rattlesnake species. Proc Biol Sci 2020; 286:20182735. [PMID: 30862287 DOI: 10.1098/rspb.2018.2735] [Citation(s) in RCA: 47] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.
Collapse
Affiliation(s)
- Giulia Zancolli
- 1 Molecular Ecology and Fisheries Genetics Laboratory, School of Natural Sciences, Bangor University , Bangor LL57 2UW , UK
| | - Juan J Calvete
- 2 Evolutionary and Translational Venomics Laboratory, CSIC , Jaume Roig 11, Valencia 46010 , Spain
| | - Michael D Cardwell
- 3 Department of Biology, San Diego State University , San Diego, CA 92182 , USA
| | - Harry W Greene
- 4 Department of Ecology and Evolutionary Biology, Cornell University , Corson Hall, Ithaca, NY 14853 , USA
| | - William K Hayes
- 5 Department of Earth and Biological Sciences, School of Medicine, Loma Linda University , Loma Linda, CA 92350 , USA
| | - Matthew J Hegarty
- 6 Institute of Biological, Environmental and Rural Sciences, Aberystwyth University , Aberystwyth SY23 3EE , UK
| | - Hans-Werner Herrmann
- 7 Wildlife Conservation and Management, School of Natural Resources and the Environment, University of Arizona , 1064 East Lowell Street (ENR2), Tucson, AZ 85721 , USA
| | - Andrew T Holycross
- 8 Natural History Collections, Arizona State University , 734 W. Alameda Drive, Tempe, AZ 85282 , USA
| | - Dominic I Lannutti
- 9 Department of Biological Sciences, University of Texas at El Paso , 500 W. University, El Paso, TX 79968 , USA
| | - John F Mulley
- 1 Molecular Ecology and Fisheries Genetics Laboratory, School of Natural Sciences, Bangor University , Bangor LL57 2UW , UK
| | - Libia Sanz
- 2 Evolutionary and Translational Venomics Laboratory, CSIC , Jaume Roig 11, Valencia 46010 , Spain
| | - Zachary D Travis
- 5 Department of Earth and Biological Sciences, School of Medicine, Loma Linda University , Loma Linda, CA 92350 , USA
| | - Joshua R Whorley
- 10 Seattle Central College, Science, Technology, Engineering & Mathematics Division , 1701 Broadway Ave. E., Seattle, WA 98122 , USA
| | - Catharine E Wüster
- 1 Molecular Ecology and Fisheries Genetics Laboratory, School of Natural Sciences, Bangor University , Bangor LL57 2UW , UK
| | - Wolfgang Wüster
- 1 Molecular Ecology and Fisheries Genetics Laboratory, School of Natural Sciences, Bangor University , Bangor LL57 2UW , UK
| |
Collapse
|
37
|
Patton AH, Margres MJ, Stahlke AR, Hendricks S, Lewallen K, Hamede RK, Ruiz-Aravena M, Ryder O, McCallum HI, Jones ME, Hohenlohe PA, Storfer A. Contemporary Demographic Reconstruction Methods Are Robust to Genome Assembly Quality: A Case Study in Tasmanian Devils. Mol Biol Evol 2020; 36:2906-2921. [PMID: 31424552 PMCID: PMC6878949 DOI: 10.1093/molbev/msz191] [Citation(s) in RCA: 55] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Reconstructing species’ demographic histories is a central focus of molecular ecology and evolution. Recently, an expanding suite of methods leveraging either the sequentially Markovian coalescent (SMC) or the site-frequency spectrum has been developed to reconstruct population size histories from genomic sequence data. However, few studies have investigated the robustness of these methods to genome assemblies of varying quality. In this study, we first present an improved genome assembly for the Tasmanian devil using the Chicago library method. Compared with the original reference genome, our new assembly reduces the number of scaffolds (from 35,975 to 10,010) and increases the scaffold N90 (from 0.101 to 2.164 Mb). Second, we assess the performance of four contemporary genomic methods for inferring population size history (PSMC, MSMC, SMC++, Stairway Plot), using the two devil genome assemblies as well as simulated, artificially fragmented genomes that approximate the hypothesized demographic history of Tasmanian devils. We demonstrate that each method is robust to assembly quality, producing similar estimates of Ne when simulated genomes were fragmented into up to 5,000 scaffolds. Overall, methods reliant on the SMC are most reliable between ∼300 generations before present (gbp) and 100 kgbp, whereas methods exclusively reliant on the site-frequency spectrum are most reliable between the present and 30 gbp. Our results suggest that when used in concert, genomic methods for reconstructing species’ effective population size histories 1) can be applied to nonmodel organisms without highly contiguous reference genomes, and 2) are capable of detecting independently documented effects of historical geological events.
Collapse
Affiliation(s)
- Austin H Patton
- School of Biological Sciences, Washington State University, Pullman, WA
| | - Mark J Margres
- School of Biological Sciences, Washington State University, Pullman, WA.,Department of Organismic and Evolutionary Biology, Harvard University, MA
| | - Amanda R Stahlke
- Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID
| | - Sarah Hendricks
- Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID
| | - Kevin Lewallen
- Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID
| | - Rodrigo K Hamede
- School of Natural Sciences, University of Tasmania, Hobart, Australia
| | | | - Oliver Ryder
- Institute for Conservation Research, San Diego, CA
| | | | - Menna E Jones
- School of Natural Sciences, University of Tasmania, Hobart, Australia
| | - Paul A Hohenlohe
- Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID
| | - Andrew Storfer
- School of Biological Sciences, Washington State University, Pullman, WA
| |
Collapse
|
38
|
Barua A, Mikheyev AS. Many Options, Few Solutions: Over 60 My Snakes Converged on a Few Optimal Venom Formulations. Mol Biol Evol 2020; 36:1964-1974. [PMID: 31220860 PMCID: PMC6736290 DOI: 10.1093/molbev/msz125] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Gene expression changes contribute to complex trait variations in both individuals and populations. However, the evolution of gene expression underlying complex traits over macroevolutionary timescales remains poorly understood. Snake venoms are proteinaceous cocktails where the expression of each toxin can be quantified and mapped to a distinct genomic locus and traced for millions of years. Using a phylogenetic generalized linear mixed model, we analyzed expression data of toxin genes from 52 snake species spanning the 3 venomous snake families and estimated phylogenetic covariance, which acts as a measure of evolutionary constraint. We find that evolution of toxin combinations is not constrained. However, although all combinations are in principle possible, the actual dimensionality of phylomorphic space is low, with envenomation strategies focused around only four major toxin families: metalloproteases, three-finger toxins, serine proteases, and phospholipases A2. Although most extant snakes prioritize either a single or a combination of major toxin families, they are repeatedly recruited and lost. We find that over macroevolutionary timescales, the venom phenotypes were not shaped by phylogenetic constraints, which include important microevolutionary constraints such as epistasis and pleiotropy, but more likely by ecological filtering that permits a small number of optimal solutions. As a result, phenotypic optima were repeatedly attained by distantly related species. These results indicate that venoms evolve by selection on biochemistry of prey envenomation, which permit diversity through parallelism, and impose strong limits, since only a few of the theoretically possible strategies seem to work well and are observed in extant snakes.
Collapse
Affiliation(s)
- Agneesh Barua
- Okinawa Institute of Science and Technology Graduate University, Onna, Japan
| | - Alexander S Mikheyev
- Okinawa Institute of Science and Technology Graduate University, Onna, Japan.,Evolutionary Genomics Research Group, Ecology and Evolution Unit, Australian National University, Canberra, Australia
| |
Collapse
|
39
|
Modahl CM, Mackessy SP. Venoms of Rear-Fanged Snakes: New Proteins and Novel Activities. Front Ecol Evol 2019. [DOI: 10.3389/fevo.2019.00279] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
|
40
|
Rautsaw RM, Hofmann EP, Margres MJ, Holding ML, Strickland JL, Mason AJ, Rokyta DR, Parkinson CL. Intraspecific sequence and gene expression variation contribute little to venom diversity in sidewinder rattlesnakes ( Crotalus cerastes). Proc Biol Sci 2019; 286:20190810. [PMID: 31266424 DOI: 10.1098/rspb.2019.0810] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Traits can evolve rapidly through changes in gene expression or protein-coding sequences. However, these forms of genetic variation can be correlated and changes to one can influence the other. As a result, we might expect traits lacking differential expression to preferentially evolve through changes in protein sequences or morphological adaptation. Given the lack of differential expression across the distribution of sidewinder rattlesnakes ( Crotalus cerastes), we tested this hypothesis by comparing the coding regions of genes expressed in the venom gland transcriptomes and fang morphology. We calculated Tajima's D and FST across four populations comparing toxin and nontoxin loci. Overall, we found little evidence of directional selection or differentiation between populations, suggesting that changes to protein sequences do not underlie the evolution of sidewinder venom or that toxins are under extremely variant selection pressures. Although low-expression toxins do not have higher sequence divergence between populations, they do have more standing variation on which selection can act. Additionally, we found significant differences in fang length among populations. The lack of differential expression and sequence divergence suggests sidewinders-given their generalist diet, moderate gene flow and environmental variation-are under stabilizing selection which functions to maintain a generalist phenotype. Overall, we demonstrate the importance of examining the relationship between gene expression and protein-coding changes to understand the evolution of complex traits.
Collapse
Affiliation(s)
- Rhett M Rautsaw
- 1 Department of Biological Sciences, Clemson University , Clemson, SC 29634 , USA
| | - Erich P Hofmann
- 1 Department of Biological Sciences, Clemson University , Clemson, SC 29634 , USA
| | - Mark J Margres
- 1 Department of Biological Sciences, Clemson University , Clemson, SC 29634 , USA
| | - Matthew L Holding
- 1 Department of Biological Sciences, Clemson University , Clemson, SC 29634 , USA.,3 Department of Biological Science, Florida State University , Tallahassee, FL 32306 , USA
| | - Jason L Strickland
- 1 Department of Biological Sciences, Clemson University , Clemson, SC 29634 , USA
| | - Andrew J Mason
- 1 Department of Biological Sciences, Clemson University , Clemson, SC 29634 , USA
| | - Darin R Rokyta
- 3 Department of Biological Science, Florida State University , Tallahassee, FL 32306 , USA
| | - Christopher L Parkinson
- 1 Department of Biological Sciences, Clemson University , Clemson, SC 29634 , USA.,2 Department of Forestry and Environmental Conservation, Clemson University , Clemson, SC 29634 , USA
| |
Collapse
|
41
|
Margres MJ, Patton A, Wray KP, Hassinger ATB, Ward MJ, Lemmon EM, Lemmon AR, Rokyta DR. Tipping the Scales: The Migration-Selection Balance Leans toward Selection in Snake Venoms. Mol Biol Evol 2019; 36:271-282. [PMID: 30395254 DOI: 10.1093/molbev/msy207] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The migration-selection interaction is the strongest determinant of whether a beneficial allele increases in frequency within a population. Results of empirical studies examining the role of gene flow in an adaptive context, however, have largely been equivocal, with examples of opposing outcomes being repeatedly documented (e.g., local adaptation with high levels of gene flow vs. gene swamping). We compared neutral genomic and venom expression divergence for three sympatric pit vipers with differing ecologies to determine if and how migration-selection disequilibria result in local adaptation. We specifically tested whether neutral differentiation predicted phenotypic differentiation within an isolation-by-distance framework. The decoupling of neutral and phenotypic differentiation would indicate selection led to adaptive divergence irrespective of migration, whereas a significant relationship between neutral and venom expression differentiation would provide evidence in favor of the constraining force of gene flow. Neutral differentiation and geographic distance predicted phenotypic differentiation only in the generalist species, indicating that selection was the predominant mechanism in the migration-selection balance underlying adaptive venom evolution in both specialists. Dispersal is thought to be a stronger influence on genetic differentiation than specialization, but our results suggest the opposite. Greater specialization may lead to greater diversification rates, and extreme spatial and temporal variation in selective pressures can favor generalist phenotypes evolving under strong stabilizing selection. Our results are consistent with these expectations, suggesting that the equivocal findings of studies examining the role of gene flow in an adaptive context may be explained by ecological specialization theory.
Collapse
Affiliation(s)
- Mark J Margres
- Department of Biological Science, Florida State University, Tallahassee, FL.,School of Biological Sciences, Washington State University, Pullman, WA.,Department of Biological Sciences, Clemson University, Clemson, SC
| | - Austin Patton
- School of Biological Sciences, Washington State University, Pullman, WA
| | - Kenneth P Wray
- Department of Biological Science, Florida State University, Tallahassee, FL
| | - Alyssa T B Hassinger
- Department of Biological Science, Florida State University, Tallahassee, FL.,Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH
| | - Micaiah J Ward
- Department of Biological Science, Florida State University, Tallahassee, FL
| | | | - Alan R Lemmon
- Department of Scientific Computing, Florida State University, Tallahassee, FL
| | - Darin R Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL
| |
Collapse
|
42
|
Evidence for Snake Venom Plasticity in a Long-Term Study with Individual Captive Bothrops atrox. Toxins (Basel) 2019; 11:toxins11050294. [PMID: 31137619 PMCID: PMC6563259 DOI: 10.3390/toxins11050294] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2019] [Revised: 05/15/2019] [Accepted: 05/21/2019] [Indexed: 12/14/2022] Open
Abstract
Variability in snake venom composition has been frequently reported and correlated to the adaptability of snakes to environmental conditions. Previous studies report plasticity for the venom phenotype. However, these observations are not conclusive, as the results were based on pooled venoms, which present high individual variability. Here we tested the hypothesis of plasticity by influence of confinement and single diet type in the venom composition of 13 adult specimens of Bothrops atrox snakes, maintained under captivity for more than three years. Individual variability in venom composition was observed in samples extracted just after the capture of the snakes. However, composition was conserved in venoms periodically extracted from nine specimens, which presented low variability restricted to the less abundant components. In a second group, composed of four snakes, drastic changes were observed in the venom samples extracted at different periods, mostly related to snake venom metalloproteinases (SVMPs), the core function toxins of B. atrox venom, which occurred approximately between 400 and 500 days in captivity. These data show plasticity in the venom phenotype during the lifetime of adult snakes maintained under captive conditions. Causes or functional consequences involved in the phenotype modification require further investigations.
Collapse
|
43
|
Harris RJ, Jenner RA. Evolutionary Ecology of Fish Venom: Adaptations and Consequences of Evolving a Venom System. Toxins (Basel) 2019; 11:E60. [PMID: 30678265 PMCID: PMC6409815 DOI: 10.3390/toxins11020060] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2018] [Revised: 01/14/2019] [Accepted: 01/18/2019] [Indexed: 01/21/2023] Open
Abstract
Research on venomous animals has mainly focused on the molecular, biochemical, and pharmacological aspects of venom toxins. However, it is the relatively neglected broader study of evolutionary ecology that is crucial for understanding the biological relevance of venom systems. As fish have convergently evolved venom systems multiple times, it makes them ideal organisms to investigate the evolutionary ecology of venom on a broader scale. This review outlines what is known about how fish venom systems evolved as a result of natural enemy interactions and about the ecological consequences of evolving a venom system. This review will show how research on the evolutionary ecology of venom in fish can aid in understanding the evolutionary ecology of animal venoms more generally. Further, understanding these broad ecological questions can shed more light on the other areas of toxinology, with applications across multiple disciplinary fields.
Collapse
Affiliation(s)
- Richard J Harris
- Venom Evolution Lab, School of Biological Sciences, The University of Queensland, St Lucia, Brisbane, Queensland 4072, Australia.
| | - Ronald A Jenner
- Department of Life Sciences, the Natural History Museum, Cromwell Road, SW7 5BD London, UK.
| |
Collapse
|
44
|
Healy K, Carbone C, Jackson AL. Snake venom potency and yield are associated with prey-evolution, predator metabolism and habitat structure. Ecol Lett 2019; 22:527-537. [PMID: 30616302 DOI: 10.1111/ele.13216] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2018] [Revised: 08/15/2018] [Accepted: 12/01/2018] [Indexed: 01/17/2023]
Abstract
Snake venom is well known for its ability to incapacitate and kill prey. Yet, potency and the amount of venom available varies greatly across species, ranging from the seemingly harmless to those capable of killing vast numbers of potential prey. This variation is poorly understood, with comparative approaches confounded by the use of atypical prey species as models to measure venom potency. Here, we account for such confounding issues by incorporating the phylogenetic similarity between a snake's diet and the species used to measure its potency. In a comparative analysis of 102 species we show that snake venom potency is generally prey-specific. We also show that venom yields are lower in species occupying three dimensional environments and increases with body size corresponding to metabolic rate, but faster than predicted from increases in prey size. These results underline the importance of physiological and environmental factors in the evolution of predator traits.
Collapse
Affiliation(s)
- Kevin Healy
- Department of Zoology, School of Natural Sciences, Trinity College Dublin, Dublin 2, Ireland.,School of Biology, University of St Andrews, St Andrews, KY16 9TH, UK.,School of Natural Sciences, National University of Ireland Galway, Galway, Ireland
| | - Chris Carbone
- Institute of Zoology, Zoological Society of London, London, UK
| | - Andrew L Jackson
- Department of Zoology, School of Natural Sciences, Trinity College Dublin, Dublin 2, Ireland
| |
Collapse
|
45
|
Kraaijeveld K, Oostra V, Liefting M, Wertheim B, de Meijer E, Ellers J. Regulatory and sequence evolution in response to selection for improved associative learning ability in Nasonia vitripennis. BMC Genomics 2018; 19:892. [PMID: 30526508 PMCID: PMC6288879 DOI: 10.1186/s12864-018-5310-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2018] [Accepted: 11/26/2018] [Indexed: 12/14/2022] Open
Abstract
Background Selection acts on the phenotype, yet only the genotype is inherited. While both the phenotypic and genotypic response to short-term selection can be measured, the link between these is a major unsolved problem in evolutionary biology, in particular for complex behavioural phenotypes. Results Here we characterize the genomic and the transcriptomic basis of associative learning ability in the parasitic wasp Nasonia vitripennis and use gene network analysis to link the two. We artificially selected for improved associative learning ability in four independent pairs of lines and identified signatures of selection across the genome. Allele frequency diverged consistently between the selected and control lines in 118 single nucleotide polymorphisms (SNPs), clustering in 51 distinct genomic regions containing 128 genes. The majority of SNPs were found in regulatory regions, suggesting a potential role for gene expression evolution. We therefore sequenced the transcriptomes of selected and control lines and identified 36 consistently differentially expressed transcripts with large changes in expression. None of the differentially expressed genes also showed sequence divergence as a result of selection. Instead, gene network analysis showed many of the genes with consistent allele frequency differences and all of the differentially expressed genes to cluster in a single co-expression network. At a functional level, both genomic and transcriptomic analyses implicated members of gene networks known to be involved in neural plasticity and cognitive processes. Conclusions Taken together, our results reveal how specific cognitive abilities can readily respond to selection via a complex interplay between regulatory and sequence evolution. Electronic supplementary material The online version of this article (10.1186/s12864-018-5310-9) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Ken Kraaijeveld
- Department of Ecological Science, Faculty of Earth and Life Sciences, Vrije Universiteit, De Boelelaan 1085, 1081, HV, Amsterdam, The Netherlands.
| | - Vicencio Oostra
- Department of Genetics, Evolution and Environment, University College London, Darwin Building, Gower Street, WC1E 6BT, London, UK
| | - Maartje Liefting
- Department of Ecological Science, Faculty of Earth and Life Sciences, Vrije Universiteit, De Boelelaan 1085, 1081, HV, Amsterdam, The Netherlands
| | - Bregje Wertheim
- Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, The Netherlands
| | - Emile de Meijer
- Leiden Genome Technology Center, Department of Human Genetics, Leiden University Medical Center, Leiden, The Netherlands
| | - Jacintha Ellers
- Department of Ecological Science, Faculty of Earth and Life Sciences, Vrije Universiteit, De Boelelaan 1085, 1081, HV, Amsterdam, The Netherlands
| |
Collapse
|
46
|
Margres MJ, Ruiz-Aravena M, Hamede R, Jones ME, Lawrance MF, Hendricks SA, Patton A, Davis BW, Ostrander EA, McCallum H, Hohenlohe PA, Storfer A. The Genomic Basis of Tumor Regression in Tasmanian Devils (Sarcophilus harrisii). Genome Biol Evol 2018; 10:3012-3025. [PMID: 30321343 PMCID: PMC6251476 DOI: 10.1093/gbe/evy229] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/10/2018] [Indexed: 02/06/2023] Open
Abstract
Understanding the genetic basis of disease-related phenotypes, such as cancer susceptibility, is crucial for the advancement of personalized medicine. Although most cancers are somatic in origin, a small number of transmissible cancers have been documented. Two such cancers have emerged in the Tasmanian devil (Sarcophilus harrisii) and now threaten the species with extinction. Recently, cases of natural tumor regression in Tasmanian devils infected with the clonally contagious cancer have been detected. We used whole-genome sequencing and FST-based approaches to identify the genetic basis of tumor regression by comparing the genomes of seven individuals that underwent tumor regression with those of three infected individuals that did not. We found three highly differentiated candidate genomic regions containing several genes related to immune response and/or cancer risk, indicating that the genomic basis of tumor regression was polygenic. Within these genomic regions, we identified putative regulatory variation in candidate genes but no nonsynonymous variation, suggesting that natural tumor regression may be driven, at least in part, by differential host expression of key loci. Comparative oncology can provide insight into the genetic basis of cancer risk, tumor development, and the pathogenicity of cancer, particularly due to our limited ability to monitor natural, untreated tumor progression in human patients. Our results support the hypothesis that host immune response is necessary for triggering tumor regression, providing candidate genes that may translate to novel treatments in human and nonhuman cancers.
Collapse
Affiliation(s)
- Mark J Margres
- School of Biological Sciences, Washington State University
| | - Manuel Ruiz-Aravena
- School of Natural Sciences, University of Tasmania, Hobart, Tasmania, Australia
| | - Rodrigo Hamede
- School of Natural Sciences, University of Tasmania, Hobart, Tasmania, Australia.,Centre for Integrative Ecology, Deakin University, Waurn Ponds, Victoria, Australia
| | - Menna E Jones
- School of Natural Sciences, University of Tasmania, Hobart, Tasmania, Australia
| | | | - Sarah A Hendricks
- Department of Biological Sciences, Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow
| | - Austin Patton
- School of Biological Sciences, Washington State University
| | - Brian W Davis
- Department of Veterinary Integrative Biosciences, Texas A&M University, College Station.,Cancer Genetics and Comparative Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland
| | - Elaine A Ostrander
- Cancer Genetics and Comparative Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland
| | - Hamish McCallum
- School of Environment, Griffith University, Nathan, Queensland, Australia
| | - Paul A Hohenlohe
- Department of Biological Sciences, Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow
| | - Andrew Storfer
- School of Biological Sciences, Washington State University
| |
Collapse
|
47
|
Hofmann EP, Rautsaw RM, Strickland JL, Holding ML, Hogan MP, Mason AJ, Rokyta DR, Parkinson CL. Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (Crotalus cerastes) lineages reveal little differential expression despite individual variation. Sci Rep 2018; 8:15534. [PMID: 30341342 PMCID: PMC6195556 DOI: 10.1038/s41598-018-33943-5] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Accepted: 10/03/2018] [Indexed: 11/24/2022] Open
Abstract
Changes in gene expression can rapidly influence adaptive traits in the early stages of lineage diversification. Venom is an adaptive trait comprised of numerous toxins used for prey capture and defense. Snake venoms can vary widely between conspecific populations, but the influence of lineage diversification on such compositional differences are unknown. To explore venom differentiation in the early stages of lineage diversification, we used RNA-seq and mass spectrometry to characterize Sidewinder Rattlesnake (Crotalus cerastes) venom. We generated the first venom-gland transcriptomes and complementary venom proteomes for eight individuals collected across the United States and tested for expression differences across life history traits and between subspecific, mitochondrial, and phylotranscriptomic hypotheses. Sidewinder venom was comprised primarily of hemorrhagic toxins, with few cases of differential expression attributable to life history or lineage hypotheses. However, phylotranscriptomic lineage comparisons more than doubled instances of significant expression differences compared to all other factors. Nevertheless, only 6.4% of toxins were differentially expressed overall, suggesting that shallow divergence has not led to major changes in Sidewinder venom composition. Our results demonstrate the need for consensus venom-gland transcriptomes based on multiple individuals and highlight the potential for discrepancies in differential expression between different phylogenetic hypotheses.
Collapse
Affiliation(s)
- Erich P Hofmann
- Clemson University, Department of Biological Sciences, Clemson, SC, 29634, USA
| | - Rhett M Rautsaw
- Clemson University, Department of Biological Sciences, Clemson, SC, 29634, USA
| | - Jason L Strickland
- Clemson University, Department of Biological Sciences, Clemson, SC, 29634, USA
- University of Central Florida, Department of Biology, Orlando, FL, 32816, USA
| | - Matthew L Holding
- Clemson University, Department of Biological Sciences, Clemson, SC, 29634, USA
- Florida State University, Department of Biological Science, Tallahassee, FL, 32306, USA
| | - Michael P Hogan
- Florida State University, Department of Biological Science, Tallahassee, FL, 32306, USA
| | - Andrew J Mason
- Clemson University, Department of Biological Sciences, Clemson, SC, 29634, USA
| | - Darin R Rokyta
- Florida State University, Department of Biological Science, Tallahassee, FL, 32306, USA
| | - Christopher L Parkinson
- Clemson University, Department of Biological Sciences, Clemson, SC, 29634, USA.
- Clemson University, Department of Forestry and Environmental Conservation, Clemson, SC, 29634, USA.
| |
Collapse
|
48
|
Durban J, Sasa M, Calvete JJ. Venom gland transcriptomics and microRNA profiling of juvenile and adult yellow-bellied sea snake, Hydrophis platurus, from Playa del Coco (Guanacaste, Costa Rica). Toxicon 2018; 153:96-105. [DOI: 10.1016/j.toxicon.2018.08.016] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Revised: 08/13/2018] [Accepted: 08/29/2018] [Indexed: 02/08/2023]
|
49
|
Ward MJ, Ellsworth SA, Hogan MP, Nystrom GS, Martinez P, Budhdeo A, Zelaya R, Perez A, Powell B, He H, Rokyta DR. Female-biased population divergence in the venom of the Hentz striped scorpion (Centruroides hentzi). Toxicon 2018; 152:137-149. [PMID: 30096334 DOI: 10.1016/j.toxicon.2018.07.026] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2018] [Revised: 06/27/2018] [Accepted: 07/31/2018] [Indexed: 11/16/2022]
Abstract
Sex-biased genes are expressed at higher levels in one sex and contribute to phenotypic differences between males and females, as well as overall phenotypic variation within and among populations. Venom has evolved primarily for predation and defense, making venom expression a highly variable phenotype as a result of local adaptation. Several scorpion species have shown both intraspecific and intersexual venom variation, and males have been observed using venom in courtship and mating, suggesting the existence of venom-specific, sex-biased genes that may contribute to population divergence. We used reversed-phase high-performance liquid chromatography (RP-HPLC), Agilent protein bioanalyzer chips, nano-liquid chromatography mass spectrometry (nLC/MS/MS), and median lethal dose (LD50) assays in fruit flies (Drosophila melanogaster) and banded crickets (Gryllodes sigillatus) to investigate proteomic and functional venom variation within and among three Florida populations of the Hentz striped scorpion (Centruroides hentzi). We found significant venom variation among populations, with females, not males, being responsible for this divergence. We also found significant variation in venom expression within populations, with males contributing more to within population variation than females. Our results provide evidence that male and female scorpions experience different natural and sexual selective pressures that have led to the expression of sex-biased venom genes and that these genes may be consequential in population divergence.
Collapse
Affiliation(s)
- Micaiah J Ward
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Schyler A Ellsworth
- 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
| | - Gunnar S Nystrom
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Paul Martinez
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Amisha Budhdeo
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Roxana Zelaya
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Alexander Perez
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Barclay Powell
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA
| | - Huan He
- Institute of Molecular Biophysics and College of Medicine, Florida State University, Tallahassee, FL, 32306, USA
| | - Darin R Rokyta
- Department of Biological Science, Florida State University, Tallahassee, FL, 32306, USA.
| |
Collapse
|
50
|
Modahl CM, Frietze S, Mackessy SP. Transcriptome-facilitated proteomic characterization of rear-fanged snake venoms reveal abundant metalloproteinases with enhanced activity. J Proteomics 2018; 187:223-234. [PMID: 30092380 DOI: 10.1016/j.jprot.2018.08.004] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2018] [Revised: 07/24/2018] [Accepted: 08/03/2018] [Indexed: 12/18/2022]
Abstract
High-throughput technologies were used to identify venom gland toxin expression and to characterize the venom proteomes of two rear-fanged snakes, Ahaetulla prasina (Asian Green Vine Snake) and Borikenophis portoricensis (Puerto Rican Racer). Sixty-nine complete toxin-coding transcripts from 12 venom protein superfamilies (A. prasina) and 50 complete coding transcripts from 11 venom protein superfamilies (B. portoricensis) were identified in the venom glands. However, only 18% (A. prasina) and 32% (B. portoricensis) of the translated protein isoforms were detected in the proteome of these venoms. Both venom gland transcriptomes and venom proteomes were dominated by P-III metalloproteinases. Three-finger toxins, cysteine-rich secretory proteins, and C-type lectins were present in moderate amounts, but other protein superfamilies showed very low abundances. Venoms contained metalloproteinase activity comparable to viperid snake venom levels, but other common venom enzymes were absent or present at negligible levels. Western blot analysis showed metalloproteinase and cysteine-rich secretory protein epitopes shared with the highly venomous Boomslang (Dispholidus typus). The abundance of metalloproteinases emphasizes the important trophic role of these toxins. Comprehensive, transcriptome-informed definition of proteomes and functional characterization of venom proteins in rear-fanged snake families help to elucidate toxin evolution and provide models for protein structure-function analyses.
Collapse
Affiliation(s)
- Cassandra M Modahl
- School of Biological Sciences, University of Northern Colorado, 501 20th St., Greeley, CO 80639-0017, USA; Department of Biological Sciences, Faculty of Science, National University of Singapore, Singapore 117543, Singapore
| | - Seth Frietze
- Department of Biomedical and Health Sciences, University of Vermont, 302 Rowell, Burlington, VT 05405, USA
| | - Stephen P Mackessy
- School of Biological Sciences, University of Northern Colorado, 501 20th St., Greeley, CO 80639-0017, USA.
| |
Collapse
|