Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome.

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Molecular approaches for structural characterization of a new potassium channel blocker from Tityus stigmurus venom: cDNA cloning, homology modeling, dynamic simulations and docking.

Potassium channels are involved in the maintenance of resting membrane potential, control of cardiac and neuronal excitability, neurotransmitters release, muscle contractility and hormone secretion. The Tityus stigmurus scorpion is widely distributed in Northeastern Brazil and known to cause severe human envenomations, inducing pain, hypoesthesia, edema, erythema, paresthesia, headaches and vomiting. Most potassium channel blocking peptides that have been purified from scorpion venoms contain 30-40 amino acids with three or four disulfide bridges. These peptides belong to α-KTx subfamily. On the other hand, the ß-KTx subfamily is poorly characterized, though it is very representative in some scorpion venoms. A transcriptomic approach of T.stigmurus scorpions developed by our group revealed the repertoire of possible molecules present in the venom, including many toxins of the ß-KTx subfamily. One of the ESTs found, named TSTI0003C has a cDNA sequence of 538 bp codifying a mature protein with 47 amino acid residues, corresponding to 5299 Da. This ß-KTx peptide is a new member of the BmTXKß-related toxins, and was here named TstKMK. The three-dimensional structure of this potassium channel toxin of the T. stigmurus scorpion was obtained by computational modeling and refined by molecular dynamic simulations. Furthermore, we have made docking simulations using a Shaker kV-1.2 potassium channel from rats as receptor model and proposed which amino acid residues and interactions could be involved in its blockade.

The complete mitochondrial genome of Bothrops jararaca (Reptilia, Serpentes, Viperidae).

The complete mitochondrial genome, containing 17,526 bp, was determined from the pitviper Bothrops jararaca. It is the first mitogenome for the most medically important genus of snake in Latin America. This mitogenome has common snake mitochondrial features such as a duplicated control region that has nearly identical sequences at two different locations of the mitogenome and a translocation of tRNA-Leu (UUR). Besides, we found a translocation of the tRNA-Pro compared to Colubridae snakes. Finally, an unusual possible duplication containing a tRNA-Phe was observed for the first time and may represent a marker of the genus.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option