Structural characteristics and evolution of the Protobothrops elegans pancreatic phospholipase A2 gene in contrast with those of Protobothrops genus venom phospholipase A2 genes

Toxinology provides multidirectional and multidimensional opportunities: A personal perspective

In nature, toxins have evolved as weapons to capture and subdue the prey or to counter predators or competitors. When they are inadvertently injected into humans, they cause symptoms ranging from mild discomfort to debilitation and death. Toxinology is the science of studying venoms and toxins that are produced by a wide variety of organisms. In the past, the structure, function and mechanisms of most abundant and/or most toxic components were characterized to understand and to develop strategies to neutralize their toxicity. With recent technical advances, we are able to evaluate and determine the toxin profiles using transcriptomes of venom glands and proteomes of tiny amounts of venom. Enormous amounts of data from these studies have opened tremendous opportunities in many directions of basic and applied research. The lower costs for profiling venoms will further fuel the expansion of toxin database, which in turn will provide greater exciting and bright opportunities in toxin research.

Discovery of the Gene Encoding a Novel Small Serum Protein (SSP) of Protobothrops flavoviridis and the Evolution of SSPs

Small serum proteins (SSPs) are low-molecular-weight proteins in snake serum with affinities for various venom proteins. Five SSPs, PfSSP-1 through PfSSP-5, have been reported in Protobothrops flavoviridis (“habu”, Pf) serum so far. Recently, we reported that the five genes encoding these PfSSPs are arranged in tandem on a single chromosome. However, the physiological functions and evolutionary origins of the five SSPs remain poorly understood. In a detailed analysis of the habu draft genome, we found a gene encoding a novel SSP, SSP-6. Structural analysis of the genes encoding SSPs and their genomic arrangement revealed the following: (1) SSP-6 forms a third SSP subgroup; (2) SSP-5 and SSP-6 were present in all snake genomes before the divergence of non-venomous and venomous snakes, while SSP-4 was acquired only by venomous snakes; (3) the composition of paralogous SSP genes in snake genomes seems to reflect snake habitat differences; and (4) the evolutionary emergence of SSP genes is probably related to the physiological functions of SSPs, with an initial snake repertoire of SSP-6 and SSP-5. SSP-4 and its derivative, SSP-3, as well as SSP-1 and SSP-2, appear to be venom-related and were acquired later.

Insights into the Evolution of a Snake Venom Multi-Gene Family from the Genomic Organization of Echis ocellatus SVMP Genes

The molecular events underlying the evolution of the Snake Venom Metalloproteinase (SVMP) family from an A Disintegrin And Metalloproteinase (ADAM) ancestor remain poorly understood. Comparative genomics may provide decisive information to reconstruct the evolutionary history of this multi-locus toxin family. Here, we report the genomic organization of Echis ocellatus genes encoding SVMPs from the PII and PI classes. Comparisons between them and between these genes and the genomic structures of Anolis carolinensis ADAM28 and E. ocellatus PIII-SVMP EOC00089 suggest that insertions and deletions of intronic regions played key roles along the evolutionary pathway that shaped the current diversity within the multi-locus SVMP gene family. In particular, our data suggest that emergence of EOC00028-like PI-SVMP from an ancestral PII(e/d)-type SVMP involved splicing site mutations that abolished both the 3′ splice AG acceptor site of intron 12* and the 5′ splice GT donor site of intron 13*, and resulted in the intronization of exon 13* and the consequent destruction of the structural integrity of the PII-SVMP characteristic disintegrin domain.

The taxonomic position and the unexpected divergence of the Habu viper, Protobothrops among Japanese subtropical islands

There are four Habu species currently recognized in Japan: Protobothrops flavoviridis from the Amami Islands and the Okinawa Islands, P. tokarensis from the Tokara Islands, P. elegans from the Yaeyama Islands and Ovophis okinabvensis from the Amami Islands and the Okinawa Islands. To clarify their taxonomic positions, we determined the complete mitochondria genome sequence (approx. 17kb) from two specimens from two different islands each for P. flavoviridis, P. tokarensis and P. elegans as well as one specimen of O. okinavensis and reconstructed the molecular phylogeny of Protobothrops using the published sequences of related species. The maximum likelihood tree showed four major species groups within Protbothrops: Group I consisting of P. cornutus, P. dabieshanensis, P. jerdonii and P. xiangchengensis; Group II consisting of P. flavoviridis and P. tokarensis; Group III consisting of P. maolensis, P. mucrosquamatus and P. elegans; Group IV consisting of P. himalayanus and P. kaubacki. Since we observed an unexpected divergence and the paraphyly of the two samples of P. flavoviridis collected from different islands, Amami-Oshima and Okinawajima within the Group II, we expanded the analysis by increasing the number of P. flavoviridis and P. tokarensis collected from 10 islands: Amami-Oshima (5 specimens), Kakeromajima (4) and Tokunoshima (4) from the Amami Islands, Okinawajima (4), Iheyajima (4), Iejima (4), Tokashikijima (4) and Kumejima (4) from the Okinawa Islands, Kodakarajima (P. tokarensis) (4) and Takarajima (P. tokarensis) (4) from the Tokara Islands. The maximum likelihood tree of the 44 samples replicated the significant divergence of P. flavoviridis between the Amami Clade including Amami-Oshima, Kakeromajima and Tokunoshima and the Okinawa Clade including Okinawajima, Iheyajima, Iejima, Tokashikijima and Kumejima. The Amami Clade also include all specimens from the Tokara Islands currently known as an independent species, P. tokarensis, suggesting the paraphyly of the taxon, P. flavoviridis. In contrast, we observed a distinct lineage of the two specimens from the Yaeyama Islands, supporting the validity of the taxon, P. elegans as an independent species. By MCMC method, we estimated the divergence time between the Amami Clade and the Okinawa Clade to be 6.51MYA, suggesting that the vicariance of the two clades preceded the geological separation of the Amami Islands and the Okinawa Islands (∼1.5MYA). As expected from the limited mobility of terrestrial reptiles including snakes, we observed high genetic divergence in Habu mtDNA among Japanese subtropical island populations.

The finding of a group IIE phospholipase A2 gene in a specified segment of Protobothrops flavoviridis genome and its possible evolutionary relationship to group IIA phospholipase A2 genes

The genes encoding group IIE phospholipase A₂, abbreviated as IIE PLA₂, and its 5' and 3' flanking regions of Crotalinae snakes such as Protobothrops flavoviridis, P. tokarensis, P. elegans, and Ovophis okinavensis, were found and sequenced. The genes consisted of four exons and three introns and coded for 22 or 24 amino acid residues of the signal peptides and 134 amino acid residues of the mature proteins. These IIE PLA₂s show high similarity to those from mammals and Colubridae snakes. The high expression level of IIE PLA₂s in Crotalinae venom glands suggests that they should work as venomous proteins. The blast analysis indicated that the gene encoding OTUD3, which is ovarian tumor domain-containing protein 3, is located in the 3' downstream of IIE PLA₂ gene. Moreover, a group IIA PLA₂ gene was found in the 5' upstream of IIE PLA₂ gene linked to the OTUD3 gene (OTUD3) in the P. flavoviridis genome. It became evident that the specified arrangement of IIA PLA₂ gene, IIE PLA₂ gene, and OTUD3 in this order is common in the genomes of humans to snakes. The present finding that the genes encoding various secretory PLA₂s form a cluster in the genomes of humans to birds is closely related to the previous finding that six venom PLA₂ isozyme genes are densely clustered in the so-called NIS-1 fragment of the P. flavoviridis genome. It is also suggested that venom IIA PLA₂ genes may be evolutionarily derived from the IIE PLA₂ gene.