Although snake venoms have been the subject of intense research, primarily because of their tremendous potential as a bioresource for design and development of therapeutic compounds, some specific groups of snakes, such as the genus Atractaspis, have been completely neglected. To date only limited number of toxins, such as sarafotoxins have been well characterized from this lineage. In order to investigate the molecular diversity of venom from Atractaspis aterrima—the slender burrowing asp, we utilized a high-throughput transcriptomic approach completed with an original bioinformatics analysis pipeline. Surprisingly, we found that Sarafotoxins do not constitute the major ingredient of the transcriptomic cocktail; rather a large number of previously well-characterized snake venom-components were identified. Notably, we recovered a large diversity of three-finger toxins (3FTxs), which were found to have evolved under the significant influence of positive selection. From the normalized and non-normalized transcriptome libraries, we were able to evaluate the relative abundance of the different toxin groups, uncover rare transcripts, and gain new insight into the transcriptomic machinery. In addition to previously characterized toxin families, we were able to detect numerous highly-transcribed compounds that possess all the key features of venom-components and may constitute new classes of toxins.
References
[1]
Fry, B.G.; Casewell, N.R.; Wüster, W.; Vidal, N.; Young, B.; Jackson, T.N. The structural and functional diversification of the Toxicofera reptile venom system. Toxicon 2012, 60, 434–448, doi:10.1016/j.toxicon.2012.02.013.
[2]
Ducancel, F. Endothelin-like peptides. Cell. Mol. Life Sci. 2005, 62, 2828–2839, doi:10.1007/s00018-005-5286-x.
[3]
Kochva, E. Atractaspis (serpentes, Atractaspididae) the burrowing asp. A mulstidisciplinary mini review. Bull. Nat. Hist. Mus. Zool. 2002, 68, 91–99.
[4]
Golani, I.; Kochva, E. Biting behaviour of Atractaspis. Copeia 1988, 792–797, doi:10.2307/1445406.
[5]
Weiser, E.; Wollberg, Z.; Kochva, E.; Lee, S.Y. Cardiotoxic effects of the venom of the burrowing asp, Atractaspis engaddensis (Atractaspididae, Ophidia). Toxicon 1984, 22, 767–774, doi:10.1016/0041-0101(84)90159-4.
[6]
Kurnik, D.; Haviv, Y.; Kochva, E. A snake bite by the Burrowing Asp, Atractaspis engaddensis. Toxicon 1999, 37, 223–227, doi:10.1016/S0041-0101(98)00166-4.
[7]
Kloog, Y.; Ambar, I.; Sokolovsky, M.; Kochva, E.; Wollberg, Z.; Bdolah, A. Sarafotoxin, a novel vasoconstrictor peptide: Phosphoinositide hydrolysis in rat heart and brain. Science 1988, 242, 268–270.
[8]
Ducancel, F.; Matre, V.; Dupont, C.; Lajeunesse, E.; Wollberg, Z.; Bdolah, A.; Kochva, E.; Boulain, J.C.; Ménez, A. Cloning and sequence analysis of cDNAs encoding precursors of sarafotoxins. Evidence for an unusual “rosary-type” organization. J. Biol. Chem. 1993, 268, 3052–3055.
[9]
sHayashi, M.A.F.; Ligny-Lemaire, C.; Wollberg, Z.; Wery, M.; Galat, A.; Ogawa, T.; Muller, B.H.; Lamthanh, H.; Doljansky, Y.; Bdolah, A.; et al. Long-sarafotoxins: Characterization of a new family of endothelin-like peptides. Peptides 2004, 25, 1243–1251, doi:10.1016/j.peptides.2004.05.010.
[10]
Quinton, L.; Le Caer, J.-P.; Phan, G.; Ligny-Lemaire, C.; Bourdais-Jomaron, J.; Ducancel, F.; Chamot-Rooke, J. Characterization of toxins within crude venoms by combined use of Fourier transform mass spectrometry and cloning. Anal. Chem. 2005, 77, 6630–6639, doi:10.1021/ac050575k.
[11]
Kochva, E.; Bdolah, A.; Wollberg, Z. Sarafotoxins and endothelins: Evolution, structure and function. Toxicon 1993, 31, 541–568, doi:10.1016/0041-0101(93)90111-U.
[12]
Yanagisawa, M.; Kurihara, H.; Kimura, S.; Tomobe, Y.; Kobayashi, M.; Mitsui, Y.; Yazaki, Y.; Goto, K.; Masaki, T. A novel potent vasoconstrictor peptide produced by vascular endothelial cells. Nature 1988, 332, 411–415, doi:10.1038/332411a0.
Fry, B.G. From genome to evenome”enMolecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins. Genome Res. 2005, 15, 403–420, doi:10.1101/gr.3228405.
[16]
Fry, B.G.; Undheim, E.A.B.; Ali, S.A.; Jackson, T.N.W.; Debono, J.; Scheib, H.; Ruder, T.; Morgenstern, D.; Cadwallader, L.; Whitehead, D.; et al. Squeezers and leaf-cutters: Differential diversification and degeneration of the venom system in toxicoferan reptiles. Mol. Cell. Proteomics 2013, 12, 1881–1899, doi:10.1074/mcp.M112.023143.
[17]
Mourier, G.; Hajj, M.; Cordier, F.; Zorba, A.; Gao, X.; Coskun, T.; Herbet, A.; Marcon, E.; Beau, F.; Delepierre, M.; et al. Pharmacological and structural characterization of long-sarafotoxins, a new family of endothelin-like peptides: Role of the C-terminus extension. Biochimie 2012, 94, 461–470, doi:10.1016/j.biochi.2011.08.014.
[18]
Fry, B.G.; Scheib, H.; van der Weerd, L.; Young, B.; McNaughtan, J.; Ramjan, S.F.R.; Vidal, N.; Poelmann, R.E.; Norman, J.A. Evolution of an arsenal: Structural and functional diversification of the venom system in the advanced snakes (Caenophidia). Mol. Cell. Proteomics MCP 2008, 7, 215–246.
[19]
Fry, B.G.; Scheib, H.; de L M Junqueira de Azevedo, I.; Silva, D.A.; Casewell, N.R. Novel transcripts in the maxillary venom glands of advanced snakes. Toxicon 2012, 59, 696–708, doi:10.1016/j.toxicon.2012.03.005.
[20]
Derrien, T.; Guigi, R.; Johnson, R. The long non-coding RNAs: A new (P)layer in the ydark matterh. Front. Genet. 2012, 1, doi:10.3389/fgene.2011.00107.
[21]
Terrat, Y.; Biass, D.; Dutertre, S.; Favreau, P.; Remm, M.; St？cklin, R.; Piquemal, D.; Ducancel, F. High-resolution picture of a venom gland transcriptome: Case study with the marine snail Conus consors. Toxicon 2012, 59, 34–46, doi:10.1016/j.toxicon.2011.10.001.
[22]
Dutertre, S.; Jin, A.; Kaas, Q.; Jones, A.; Alewood, P.F.; Lewis, R.J. Deep venomics reveals the mechanism for expanded peptide diversity in cone snail venom. Mol. Cell. Proteomics 2013, 12, 312–329, doi:10.1074/mcp.M112.021469.
[23]
Becker, A.; Dowdle, E.B.; Hechler, U.; Kauser, K.; Donner, P.; Schleuning, W.D. Bibrotoxin, a novel member of the endothelin/sarafotoxin peptide family, from the venom of the burrowing asp Atractaspis bibroni. FEBS Lett. 1993, 315, 100–103, doi:10.1016/0014-5793(93)81142-M.
[24]
Fry, B.G.; Wüster, W.; Kini, R.M.; Brusic, V.; Khan, A.; Venkataraman, D.; Rooney, A.P. Molecular evolution and phylogeny of elapid snake venom three-finger toxins. J. Mol. Evol. 2003, 57, 110–129, doi:10.1007/s00239-003-2461-2.
[25]
Fry, B.G.; Lumsden, N.G.; Wüster, W.; Wickramaratna, J.C.; Hodgson, W.C.; Kini, R.M. Isolation of a neurotoxin (alpha-colubritoxin) from a nonvenomous colubrid: Evidence for early origin of venom in snakes. J. Mol. Evol. 2003, 57, 446–452, doi:10.1007/s00239-003-2497-3.
[26]
Borodovsky, M.; Lomsadze, A. Eukaryotic gene prediction using GeneMark.hmm-E and GeneMark-ES. Curr. Protoc. Bioinformatics 2011, doi:10.1002/0471250953.bi0406s35.
[27]
Durban, J.; Juárez, P.; Angulo, Y.; Lomonte, B.; Flores-Diaz, M.; Alape-Giraz, A.; Sasa, M.; Sanz, L.; Gutiérrez, J.M.; Dopazo, J.; et al. Profiling the venom gland transcriptomes of Costa Rican snakes by 454 pyrosequencing. BMC Genomics 2011, 12, doi:10.1186/1471-2164-12-259.
[28]
Petersen, T.N.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. Nat. Methods 2011, 8, 785–786, doi:10.1038/nmeth.1701.
[29]
Jungo, F.; Bougueleret, L.; Xenarios, I.; Poux, S. The UniProtKB/Swiss-Prot Tox-Prot program: A central hub of integrated venom protein data. Toxicon 2012, 60, 551–557.
[30]
Li, W.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659, doi:10.1093/bioinformatics/btl158.
[31]
Halary, S.; Leigh, J.W.; Cheaib, B.; Lopez, P.; Bapteste, E. Network analyses structure genetic diversity in independent genetic worlds. Proc. Natl. Acad. Sci. USA 2010, 107, 127–132.
[32]
Smoot, M.E.; Ono, K.; Ruscheinski, J.; Wang, P.-L.; Ideker, T. Cytoscape 2.8: New features for data integration and network visualization. Bioinformatics 2011, 27, 431–432.
[33]
Goldman, N.; Yang, Z. A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol. Biol. Evol. 1994, 11, 725–736.
[34]
Yang, Z. Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol. Biol. Evol. 1998, 15, 568–573, doi:10.1093/oxfordjournals.molbev.a025957.
[35]
Yang, Z. PAML 4: Phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 2007, 24, 1586–1591, doi:10.1093/molbev/msm088.
[36]
Nielsen, R.; Yang, Z. Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics 1998, 148, 929–936.
[37]
Yang, Z.; Wong, W.S.W.; Nielsen, R. Bayes empirical bayes inference of amino acid sites under positive selection. Mol. Biol. Evol. 2005, 22, 1107–1118, doi:10.1093/molbev/msi097.
Kosakovsky Pond, S.L.; Murrell, B.; Fourment, M.; Frost, S.D.W.; Delport, W.; Scheffler, K. A random effects branch-site model for detecting episodic diversifying selection. Mol. Biol. Evol. 2011, 28, 3033–3043.
[42]
Kelley, L.A.; Sternberg, M.J.E. Protein structure prediction on the Web: A case study using the Phyre server. Nat. Protoc. 2009, 4, 363–371, doi:10.1038/nprot.2009.2.
[43]
DeLano, WL. The PyMOL Molecular Graphics System; DeLano Scientific: San Carlos, CA, USA, 2002.
[44]
Armon, A.; Graur, D.; Ben-Tal, N. ConSurf: An algorithmic tool for the identification of functional regions in proteins by surface mapping of phylogenetic information. J. Mol. Biol. 2001, 307, 447–463, doi:10.1006/jmbi.2000.4474.
[45]
Rokyta, D.R.; Lemmon, A.R.; Margres, M.J.; Aronow, K. The venom-gland transcriptome of the eastern diamondback rattlesnake (Crotalus adamanteus). BMC Genomics 2012, 13, doi:10.1186/1471-2164-13-312.
[46]
Trevisan-Silva, D.; Gremski, L.H.; Chaim, O.M.; da Silveira, R.B.; Meissner, G.O.; Mangili, O.C.; Barbaro, K.C.; Gremski, W.; Veiga, S.S.; Senff-Ribeiro, A. Astacin-like metalloproteases are a gene family of toxins present in the venom of different species of the brown spider (genus Loxosceles). Biochimie 2010, 92, 21–32, doi:10.1016/j.biochi.2009.10.003.
[47]
Kawaguchi, M.; Yasumasu, S.; Shimizu, A.; Hiroi, J.; Yoshizaki, N.; Nagata, K.; Tanokura, M.; Iuchi, I. Purification and gene cloning of Fundulus heteroclitus hatching enzyme. A hatching enzyme system composed of high choriolytic enzyme and low choriolytic enzyme is conserved between two different teleosts, Fundulus heteroclitus and medaka Oryzias latipes. FEBS J. 2005, 272, 4315–4326.
