Venom gland transcriptomes of two elapid snakes (Bungarus multicinctus and Naja atra) and evolution of toxin genes

<p>Abstract</p> <p>Background</p> <p>Kraits (genus <it>Bungarus</it>) and cobras (genus <it>Naja</it>) are two representative toxic genera of elapids in the old world. Although they are closely related genera and both of their venoms are very toxic, the compositions of their venoms are very different. To unveil their detailed venoms and their evolutionary patterns, we constructed venom gland cDNA libraries and genomic bacterial artificial chromosome (BAC) libraries for <it>Bungarus multicinctus </it>and <it>Naja atra</it>, respectively. We sequenced about 1500 cDNA clones for each of the venom cDNA libraries and screened BAC libraries of the two snakes by blot analysis using four kinds of toxin probes; <it>i.e</it>., three-finger toxin (3FTx), phospholipase A2 (PLA2), kunitz-type protease inhibitor (Kunitz), and natriuretic peptide (NP).</p> <p>Results</p> <p>In total, 1092 valid expressed sequences tags (ESTs) for <it>B. multicinctus </it>and 1166 ESTs for <it>N. atra </it>were generated. About 70% of these ESTs can be annotated as snake toxin transcripts. 3FTx (64.5%) and <it>β </it>bungarotoxin (25.1%) comprise the main toxin classes in <it>B. multicinctus</it>, while 3FTx (95.8%) is the dominant toxin in <it>N. atra</it>. We also observed several less abundant venom families in <it>B. multicinctus </it>and <it>N. atra</it>, such as PLA2, C-type lectins, and Kunitz. Peculiarly a cluster of NP precursors with tandem NPs was detected in <it>B. multicinctus</it>. A total of 71 positive toxin BAC clones in <it>B. multicinctus </it>and <it>N. atra </it>were identified using four kinds of toxin probes (3FTx, PLA2, Kunitz, and NP), among which 39 3FTx-postive BACs were sequenced to reveal gene structures of 3FTx toxin genes.</p> <p>Conclusions</p> <p>Based on the toxin ESTs and 3FTx gene sequences, the major components of <it>B. multicinctus </it>venom transcriptome are neurotoxins, including long chain alpha neurotoxins (<it>α</it>-ntx) and the recently originated <it>β </it>bungarotoxin, whereas the <it>N. atra </it>venom transcriptome mainly contains 3FTxs with cytotoxicity and neurotoxicity (short chain <it>α</it>-ntx). The data also revealed that tandem duplications contributed the most to the expansion of toxin multigene families. Analysis of nonsynonymous to synonymous nucleotide substitution rate ratios (<it>dN</it>/<it>dS</it>) indicates that not only multigene toxin families but also other less abundant toxins might have been under rapid diversifying evolution.</p

Diversification rates and phenotypic evolution in venomous snakes (Elapidae)

Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original author and source are credited.The relationship between rates of diversification and of body size change (a common proxy for phenotypic evolution) was investigated across Elapidae, the largest radiation of highly venomous snakes. Time-calibrated phylogenetic trees for 175 species of elapids (more than 50% of known taxa) were constructed using seven mitochondrial and nuclear genes. Analyses using these trees revealed no evidence for a link between speciation rates and changes in body size. Two clades (Hydrophis, Micrurus) show anomalously high rates of diversification within Elapidae, yet exhibit rates of body size evolution almost identical to the general elapid ‘background’ rate. Although correlations between speciation rates and rates of body size change exist in certain groups (e.g. ray-finned fishes, passerine birds), the two processes appear to be uncoupled in elapid snakes. There is also no detectable shift in diversification dynamics associated with the colonization of Australasia, which is surprising given that elapids appear to be the first clade of venomous snakes to reach the continent

Diversification of a single ancestral gene into a successful toxin superfamily in highly venomous Australian funnel-web spiders

Background: Spiders have evolved pharmacologically complex venoms that serve to rapidly subdue prey and deter predators. The major toxic factors in most spider venoms are small, disulfide-rich peptides. While there is abundant evidence that snake venoms evolved by recruitment of genes encoding normal body proteins followed by extensive gene duplication accompanied by explosive structural and functional diversification, the evolutionary trajectory of spider-venom peptides is less clear. Results: Here we present evidence of a spider-toxin superfamily encoding a high degree of sequence and functional diversity that has evolved via accelerated duplication and diversification of a single ancestral gene. The peptides within this toxin superfamily are translated as prepropeptides that are posttranslationally processed to yield the mature toxin. The N-terminal signal sequence, as well as the protease recognition site at the junction of the propeptide and mature toxin are conserved, whereas the remainder of the propeptide and mature toxin sequences are variable. All toxin transcripts within this superfamily exhibit a striking cysteine codon bias. We show that different pharmacological classes of toxins within this peptide superfamily evolved under different evolutionary selection pressures. Conclusions: Overall, this study reinforces the hypothesis that spiders use a combinatorial peptide library strategy to evolve a complex cocktail of peptide toxins that target neuronal receptors and ion channels in prey and predators. We show that the ω-hexatoxins that target insect voltage-gated calcium channels evolved under the influence of positive Darwinian selection in an episodic fashion, whereas the κ-hexatoxins that target insect calcium-activated potassium channels appear to be under negative selection. A majority of the diversifying sites in the ω-hexatoxins are concentrated on the molecular surface of the toxins, thereby facilitating neofunctionalisation leading to new toxin pharmacology. © 2014 Pineda et al.; licensee BioMed Central Ltd

Computational Biology in Costa Rica: The Role of a Small Country in the Global Context of Bioinformatics

Introduction: The successful development of high throughput methods for DNA sequencing, transcriptomics, proteomics, and other –omics, has contributed to the emergence of novel possibilities for the examination of complex biological systems through computational analysis. These fields have witnessed unprecedented advances in high income countries. Nevertheless, the role of other nations needs to be examined in order to delineate their contribution within the global context of bioinformatics. Previous articles have focused on the expansion of Computational Biology in Brazil and Mexico [1],[2], two of the largest Latin American countries, and which have shown political commitment to foster their scientific development. Costa Rica is a small Central American country with a population of 4 million, with its territory 164 and 38 times smaller than Brazil and Mexico, respectively. Thus, it is interesting to visualize the possibilities and challenges of this low-income country in the context of the global bioinformatics endeavor.UCR::Vicerrectoría de Investigación::Unidades de Investigación::Ciencias de la Salud::Instituto Clodomiro Picado (ICP

Network adaptation improves temporal representation of naturalistic stimuli in drosophila eye: II Mechanisms

Retinal networks must adapt constantly to best present the ever changing visual world to the brain. Here we test the hypothesis that adaptation is a result of different mechanisms at several synaptic connections within the network. In a companion paper (Part I), we showed that adaptation in the photoreceptors (R1-R6) and large monopolar cells (LMC) of the Drosophila eye improves sensitivity to under-represented signals in seconds by enhancing both the amplitude and frequency distribution of LMCs' voltage responses to repeated naturalistic contrast series. In this paper, we show that such adaptation needs both the light-mediated conductance and feedback-mediated synaptic conductance. A faulty feedforward pathway in histamine receptor mutant flies speeds up the LMC output, mimicking extreme light adaptation. A faulty feedback pathway from L2 LMCs to photoreceptors slows down the LMC output, mimicking dark adaptation. These results underline the importance of network adaptation for efficient coding, and as a mechanism for selectively regulating the size and speed of signals in neurons. We suggest that concert action of many different mechanisms and neural connections are responsible for adaptation to visual stimuli. Further, our results demonstrate the need for detailed circuit reconstructions like that of the Drosophila lamina, to understand how networks process information

Tissue Localization and Extracellular Matrix Degradation by PI, PII and PIII Snake Venom Metalloproteinases: Clues on the Mechanisms of Venom-Induced Hemorrhage

20 páginas, 4 figuras, 3 tablas y 7 tablas en material suplementario.Snake venom hemorrhagic metalloproteinases (SVMPs) of the PI, PII and PIII classes were compared in terms of tissue localization and their ability to hydrolyze basement membrane components in vivo, as well as by a proteomics analysis of exudates collected in tissue injected with these enzymes. Immunohistochemical analyses of co-localization of these SVMPs with type IV collagen revealed that PII and PIII enzymes co-localized with type IV collagen in capillaries, arterioles and post-capillary venules to a higher extent than PI SVMP, which showed a more widespread distribution in the tissue. The patterns of hydrolysis by these three SVMPs of laminin, type VI collagen and nidogen in vivo greatly differ, whereas the three enzymes showed a similar pattern of degradation of type IV collagen, supporting the concept that hydrolysis of this component is critical for the destabilization of microvessel structure leading to hemorrhage. Proteomic analysis of wound exudate revealed similarities and differences between the action of the three SVMPs. Higher extent of proteolysis was observed for the PI enzyme regarding several extracellular matrix components and fibrinogen, whereas exudates from mice injected with PII and PIII SVMPs had higher amounts of some intracellular proteins. Our results provide novel clues for understanding the mechanisms by which SVMPs induce damage to the microvasculature and generate hemorrhage.This work was performed in partial fulfillment of the requirements for the PhD degree for Cristina Herrera at Universidad de Costa Rica.Peer reviewe

Deathly Drool: Evolutionary and Ecological Basis of Septic Bacteria in Komodo Dragon Mouths

Komodo dragons, the world's largest lizard, dispatch their large ungulate prey by biting and tearing flesh. If a prey escapes, oral bacteria inoculated into the wound reputedly induce a sepsis that augments later prey capture by the same or other lizards. However, the ecological and evolutionary basis of sepsis in Komodo prey acquisition is controversial. Two models have been proposed. The “bacteria as venom” model postulates that the oral flora directly benefits the lizard in prey capture irrespective of any benefit to the bacteria. The “passive acquisition” model is that the oral flora of lizards reflects the bacteria found in carrion and sick prey, with no relevance to the ability to induce sepsis in subsequent prey. A third model is proposed and analyzed here, the “lizard-lizard epidemic” model. In this model, bacteria are spread indirectly from one lizard mouth to another. Prey escaping an initial attack act as vectors in infecting new lizards. This model requires specific life history characteristics and ways to refute the model based on these characteristics are proposed and tested. Dragon life histories (some details of which are reported here) prove remarkably consistent with the model, especially that multiple, unrelated lizards feed communally on large carcasses and that escaping, wounded prey are ultimately fed on by other lizards. The identities and evolutionary histories of bacteria in the oral flora may yield the most useful additional insights for further testing the epidemic model and can now be obtained with new technologies

Discovery of a Distinct Superfamily of Kunitz-Type Toxin (KTT) from Tarantulas

BACKGROUND: Kuntiz-type toxins (KTTs) have been found in the venom of animals such as snake, cone snail and sea anemone. The main ancestral function of Kunitz-type proteins was the inhibition of a diverse array of serine proteases, while toxic activities (such as ion-channel blocking) were developed under a variety of Darwinian selection pressures. How new functions were grafted onto an old protein scaffold and what effect Darwinian selection pressures had on KTT evolution remains a puzzle. PRINCIPAL FINDINGS: Here we report the presence of a new superfamily of ktts in spiders (TARANTULAS: Ornithoctonus huwena and Ornithoctonus hainana), which share low sequence similarity to known KTTs and is clustered in a distinct clade in the phylogenetic tree of KTT evolution. The representative molecule of spider KTTs, HWTX-XI, purified from the venom of O. huwena, is a bi-functional protein which is a very potent trypsin inhibitor (about 30-fold more strong than BPTI) as well as a weak Kv1.1 potassium channel blocker. Structural analysis of HWTX-XI in 3-D by NMR together with comparative function analysis of 18 expressed mutants of this toxin revealed two separate sites, corresponding to these two activities, located on the two ends of the cone-shape molecule of HWTX-XI. Comparison of non-synonymous/synonymous mutation ratios (omega) for each site in spider and snake KTTs, as well as PBTI like body Kunitz proteins revealed high Darwinian selection pressure on the binding sites for Kv channels and serine proteases in snake, while only on the proteases in spider and none detected in body proteins, suggesting different rates and patterns of evolution among them. The results also revealed a series of key events in the history of spider KTT evolution, including the formation of a novel KTT family (named sub-Kuntiz-type toxins) derived from the ancestral native KTTs with the loss of the second disulfide bridge accompanied by several dramatic sequence modifications. CONCLUSIONS/SIGNIFICANCE: These finding illustrate that the two activity sites of Kunitz-type toxins are functionally and evolutionally independent and provide new insights into effects of Darwinian selection pressures on KTT evolution, and mechanisms by which new functions can be grafted onto old protein scaffolds

Adaptive Evolution of the Venom-Targeted vWF Protein in Opossums that Eat Pitvipers

The rapid evolution of venom toxin genes is often explained as the result of a biochemical arms race between venomous animals and their prey. However, it is not clear that an arms race analogy is appropriate in this context because there is no published evidence for rapid evolution in genes that might confer toxin resistance among routinely envenomed species. Here we report such evidence from an unusual predator-prey relationship between opossums (Marsupialia: Didelphidae) and pitvipers (Serpentes: Crotalinae). In particular, we found high ratios of replacement to silent substitutions in the gene encoding von Willebrand Factor (vWF), a venom-targeted hemostatic blood protein, in a clade of opossums known to eat pitvipers and to be resistant to their hemorrhagic venom. Observed amino-acid substitutions in venom-resistant opossums include changes in net charge and hydrophobicity that are hypothesized to weaken the bond between vWF and one of its toxic snake-venom ligands, the C-type lectin-like protein botrocetin. Our results provide the first example of rapid adaptive evolution in any venom-targeted molecule, and they support the notion that an evolutionary arms race might be driving the rapid evolution of snake venoms. However, in the arms race implied by our results, venomous snakes are prey, and their venom has a correspondingly defensive function in addition to its usual trophic role