The evolution and structure of snake venom phosphodiesterase (svPDE) highlight its importance in venom actions

  1. Cheng-Tsung Pan
  2. Chien-Chu Lin
  3. I-Jin Lin
  4. Kun-Yi Chien
  5. Yeong-Shin Lin  Is a corresponding author
  6. Hsiao-Han Chang  Is a corresponding author
  7. Wen-Guey Wu  Is a corresponding author
  1. National Tsing Hua University, Taiwan
  2. Chang Gung University, Taiwan
  3. National Yang Ming Chiao Tung University, Taiwan

Abstract

For decades, studies of snake venoms focused on the venom-ome-specific toxins (VSTs). VSTs are dominant soluble proteins believed to contribute to the main venomous effects and emerged into gene clusters for fast adaptation and diversification of snake venoms. However, the conserved minor venom components, such as snake venom phosphodiesterase (svPDE), remain largely unexplored. Here, we focus on svPDE by genomic and transcriptomic analysis across snake clades and demonstrate that soluble svPDE is co-opted from the ancestral membrane-attached ENPP3 (ectonucleotide pyrophosphatase/phosphodiesterase 3) gene by replacing the original 5' exon with the exon encoding a signal peptide. Notably, the exons, promoters and transcription/translation starts have been replaced multiple times during snake evolution, suggesting the evolutionary necessity of svPDE. The structural and biochemical analyses also show that svPDE shares the similar functions with ENPP family, suggesting its perturbation to the purinergic signaling and insulin transduction in venomous effects.

Data availability

The CDS and peptides of N. atra svPDE revealed in this study have been presented in Figure 1-figure supplement 1. The signal peptides of svPDE in different snake species identified with the public genome assemblies have been shown in Figure 1-figure supplement 3. The newly identified first exon for each species could be retrieved by searching its corresponding genome assembly (NCBI accessions in Supplementary File 1).

Article and author information

Author details

  1. Cheng-Tsung Pan

    Institute of Bioinformatics and Structural Biology, National Tsing Hua University, Hsinchu, Taiwan
    Competing interests
    The authors declare that no competing interests exist.
  2. Chien-Chu Lin

    Institute of Bioinformatics and Structural Biology, National Tsing Hua University, Hsinchu, Taiwan
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0647-427X
  3. I-Jin Lin

    Institute of Bioinformatics and Structural Biology, National Tsing Hua University, Hsinchu, Taiwan
    Competing interests
    The authors declare that no competing interests exist.
  4. Kun-Yi Chien

    Department of Biochemistry and Molecular Biology, Chang Gung University, Taoyuan, Taiwan
    Competing interests
    The authors declare that no competing interests exist.
  5. Yeong-Shin Lin

    Institute of Bioinformatics and Systems Biology, National Yang Ming Chiao Tung University, Hsinchu, Taiwan
    For correspondence
    yslin@nycu.edu.tw
    Competing interests
    The authors declare that no competing interests exist.
  6. Hsiao-Han Chang

    Institute of Bioinformatics and Structural Biology, National Tsing Hua University, Hsinchu, Taiwan
    For correspondence
    hhchang@life.nthu.edu.tw
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8016-1530
  7. Wen-Guey Wu

    Institute of Bioinformatics and Structural Biology, National Tsing Hua University, Hsinchu, Taiwan
    For correspondence
    wgwu@life.nthu.edu.tw
    Competing interests
    The authors declare that no competing interests exist.

Funding

Ministry of Education (Higher Education SPROUT Project)

  • I-Jin Lin
  • Hsiao-Han Chang
  • Wen-Guey Wu

Ministry of Education (Yushan Scholar Program)

  • Cheng-Tsung Pan
  • Hsiao-Han Chang

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Ethics

Animal experimentation: One male Taiwan cobra, raised by Tainan World Snake King Educational Farm, was donated by Tainan County Government and the permit was granted by Hsinchu City Animal Protection and Health Inspection office for this study. The animal use protocol was approved by the Institutional Animal Care and Use Committee of National Tsing Hua University (IACUC Protocol No. 11003H010 and Approval No. 110010).

Copyright

© 2023, Pan et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,291
    views
  • 182
    downloads
  • 2
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Cheng-Tsung Pan
  2. Chien-Chu Lin
  3. I-Jin Lin
  4. Kun-Yi Chien
  5. Yeong-Shin Lin
  6. Hsiao-Han Chang
  7. Wen-Guey Wu
(2023)
The evolution and structure of snake venom phosphodiesterase (svPDE) highlight its importance in venom actions
eLife 12:e83966.
https://doi.org/10.7554/eLife.83966

Share this article

https://doi.org/10.7554/eLife.83966

Further reading

    1. Evolutionary Biology
    Mattias Siljestam, Claus Rueffler
    Research Article

    The majority of highly polymorphic genes are related to immune functions and with over 100 alleles within a population, genes of the major histocompatibility complex (MHC) are the most polymorphic loci in vertebrates. How such extraordinary polymorphism arose and is maintained is controversial. One possibility is heterozygote advantage (HA), which can in principle maintain any number of alleles, but biologically explicit models based on this mechanism have so far failed to reliably predict the coexistence of significantly more than ten alleles. We here present an eco-evolutionary model showing that evolution can result in the emergence and maintenance of more than 100 alleles under HA if the following two assumptions are fulfilled: first, pathogens are lethal in the absence of an appropriate immune defence; second, the effect of pathogens depends on host condition, with hosts in poorer condition being affected more strongly. Thus, our results show that HA can be a more potent force in explaining the extraordinary polymorphism found at MHC loci than currently recognized.

    1. Evolutionary Biology
    Matthew Osmond, Graham Coop
    Research Article

    Spatial patterns in genetic diversity are shaped by individuals dispersing from their parents and larger-scale population movements. It has long been appreciated that these patterns of movement shape the underlying genealogies along the genome leading to geographic patterns of isolation by distance in contemporary population genetic data. However, extracting the enormous amount of information contained in genealogies along recombining sequences has, until recently, not been computationally feasible. Here we capitalize on important recent advances in genome-wide gene-genealogy reconstruction and develop methods to use thousands of trees to estimate per-generation dispersal rates and to locate the genetic ancestors of a sample back through time. We take a likelihood approach in continuous space using a simple approximate model (branching Brownian motion) as our prior distribution of spatial genealogies. After testing our method with simulations we apply it to Arabidopsis thaliana. We estimate a dispersal rate of roughly 60km2 per generation, slightly higher across latitude than across longitude, potentially reflecting a northward post-glacial expansion. Locating ancestors allows us to visualize major geographic movements, alternative geographic histories, and admixture. Our method highlights the huge amount of information about past dispersal events and population movements contained in genome-wide genealogies.