Gene duplication of coagulation factor V and origin of venom prothrombin activator in Pseudonaja textilis snake

ptFVa (Pseudonaja Textilis Venom-Derived Factor Va) Retains Structural Integrity Following Proteolysis by Activated Protein C

Supplemental Digital Content is available in the text.

Objective:

The Australian snake venom ptFV (Pseudonaja textilis venom-derived factor V) variant retains cofactor function despite APC (activated protein C)-dependent proteolysis. Here, we aimed to unravel the mechanistic principles by determining the role of the absent Arg306 cleavage site that is required for the inactivation of FVa (mammalian factor Va).

Approach and Results:

Our findings show that in contrast to human FVa, APC-catalyzed proteolysis of ptFVa at Arg306 and Lys507 does not abrogate ptFVa cofactor function. Remarkably, the structural integrity of APC-proteolyzed ptFVa is maintained indicating that stable noncovalent interactions prevent A2-domain dissociation. Using Molecular Dynamics simulations, we uncovered key regions located in the A1 and A2 domain that may be at the basis of this remarkable characteristic.

Conclusions:

Taken together, we report a completely novel role for uniquely adapted regions in ptFVa that prevent A2 domain dissociation. As such, these results challenge our current understanding by which strict regulatory mechanisms control FVa activity.

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.

Deep venomics of the Pseudonaja genus reveals inter- and intra-specific variation

Australian elapid venom remains an under-investigated resource of novel bioactive peptides. In this study, the venom gland transcriptomes and proteomes of the Australian western brown snakes, Pseudonaja aspidorhyncha and Pseudonaja nuchalis, were compared to Pseudonaja textilis. A deep venomics strategy incorporating high throughput 454 pyrosequencing gave a total of 200,911 raw reads for the three venoms. Subsequent annotation identified 5716 transcripts from 20 different toxin families with inter-specific variation between species observed in eight of the less abundant families. Integration of each venom proteome with the corresponding annotated reads identified 65 isoforms from six toxin families; high sequence coverage highlighted subtle differences between sequences and intra and inter-specific variation between species. High quality MS/MS data identified unusual glycoforms with natriuretic peptides from P. aspidorhyncha and P. nuchaliscontaining O-linked trisaccharides with high homology to the glycosylated region of TNPc. Molecular evolutionary assessments indicated the accelerated evolution of all toxin families with the exception of both natriuretic peptides and P. aspidorhyncha PLA2s that were found to be evolutionarily constrained under purifying selection pressures. This study has revealed a wide range of novel peptide sequences from six bioactive peptide families and highlights the subtle differences between toxins in these closely related species.

BIOLOGICAL SIGNIFICANCE

Mining Australia's vastly untapped source of toxins from its venomous creatures has been significantly advanced by employing deep venomics methodology. Technological advances in transcriptome analysis using next generation sequencing platforms and proteome analysis by highly sensitive tandem mass spectrometry allowed a more comprehensive interrogation of three underinvestigated brown snake (Pseudonaja) venoms uncovering many novel peptide sequences that are unique to these closely related species. This generic strategy will provide invaluable information when applied to other venomous snakes for a deeper understanding of venom composition, envenomation, venom evolution, as well as identifying research tools and drug leads.

Gene duplications are extensive and contribute significantly to the toxic proteome of nematocysts isolated from Acropora digitifera (Cnidaria: Anthozoa: Scleractinia)

Background

Gene duplication followed by adaptive selection is a well-accepted process leading to toxin diversification in venoms. However, emergent genomic, transcriptomic and proteomic evidence now challenges this role to be at best equivocal to other processess . Cnidaria are arguably the most ancient phylum of the extant metazoa that are venomous and such provide a definitive ancestral anchor to examine the evolution of this trait.

Methods

Here we compare predicted toxins from the translated genome of the coral Acropora digitifera to putative toxins revealed by proteomic analysis of soluble proteins discharged from nematocysts, to determine the extent to which gene duplications contribute to venom innovation in this reef-building coral species. A new bioinformatics tool called HHCompare was developed to detect potential gene duplications in the genomic data, which is made freely available (https://github.com/rgacesa/HHCompare).

Results

A total of 55 potential toxin encoding genes could be predicted from the A. digitifera genome, of which 36 (65 %) had likely arisen by gene duplication as evinced using the HHCompare tool and verified using two standard phylogeny methods. Surprisingly, only 22 % (12/55) of the potential toxin repertoire could be detected following rigorous proteomic analysis, for which only half (6/12) of the toxin proteome could be accounted for as peptides encoded by the gene duplicates. Biological activities of these toxins are dominatedby putative phospholipases and toxic peptidases.

Conclusions

Gene expansions in A. digitifera venom are the most extensive yet described in any venomous animal, and gene duplication plays a significant role leading to toxin diversification in this coral species. Since such low numbers of toxins were detected in the proteome, it is unlikely that the venom is evolving rapidly by prey-driven positive natural selection. Rather we contend that the venom has a defensive role deterring predation or harm from interspecific competition and overgrowth by fouling organisms. Factors influencing translation of toxin encoding genes perhaps warrants more profound experimental consideration.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1976-4) contains supplementary material, which is available to authorized users.

Transcriptome and venom proteome of the box jellyfish Chironex fleckeri

BACKGROUND

The box jellyfish, Chironex fleckeri, is the largest and most dangerous cubozoan jellyfish to humans. It produces potent and rapid-acting venom and its sting causes severe localized and systemic effects that are potentially life-threatening. In this study, a combined transcriptomic and proteomic approach was used to identify C. fleckeri proteins that elicit toxic effects in envenoming.

RESULTS

More than 40,000,000 Illumina reads were used to de novo assemble ∼ 34,000 contiguous cDNA sequences and ∼ 20,000 proteins were predicted based on homology searches, protein motifs, gene ontology and biological pathway mapping. More than 170 potential toxin proteins were identified from the transcriptome on the basis of homology to known toxins in publicly available sequence databases. MS/MS analysis of C. fleckeri venom identified over 250 proteins, including a subset of the toxins predicted from analysis of the transcriptome. Potential toxins identified using MS/MS included metalloproteinases, an alpha-macroglobulin domain containing protein, two CRISP proteins and a turripeptide-like protease inhibitor. Nine novel examples of a taxonomically restricted family of potent cnidarian pore-forming toxins were also identified. Members of this toxin family are potently haemolytic and cause pain, inflammation, dermonecrosis, cardiovascular collapse and death in experimental animals, suggesting that these toxins are responsible for many of the symptoms of C. fleckeri envenomation.

CONCLUSIONS

This study provides the first overview of a box jellyfish transcriptome which, coupled with venom proteomics data, enhances our current understanding of box jellyfish venom composition and the molecular structure and function of cnidarian toxins. The generated data represent a useful resource to guide future comparative studies, novel protein/peptide discovery and the development of more effective treatments for jellyfish stings in humans. (Length: 300).

Restriction and recruitment-gene duplication and the origin and evolution of snake venom toxins

Snake venom has been hypothesized to have originated and diversified through a process that involves duplication of genes encoding body proteins with subsequent recruitment of the copy to the venom gland, where natural selection acts to develop or increase toxicity. However, gene duplication is known to be a rare event in vertebrate genomes, and the recruitment of duplicated genes to a novel expression domain (neofunctionalization) is an even rarer process that requires the evolution of novel combinations of transcription factor binding sites in upstream regulatory regions. Therefore, although this hypothesis concerning the evolution of snake venom is very unlikely and should be regarded with caution, it is nonetheless often assumed to be established fact, hindering research into the true origins of snake venom toxins. To critically evaluate this hypothesis, we have generated transcriptomic data for body tissues and salivary and venom glands from five species of venomous and nonvenomous reptiles. Our comparative transcriptomic analysis of these data reveals that snake venom does not evolve through the hypothesized process of duplication and recruitment of genes encoding body proteins. Indeed, our results show that many proposed venom toxins are in fact expressed in a wide variety of body tissues, including the salivary gland of nonvenomous reptiles and that these genes have therefore been restricted to the venom gland following duplication, not recruited. Thus, snake venom evolves through the duplication and subfunctionalization of genes encoding existing salivary proteins. These results highlight the danger of the elegant and intuitive “just-so story” in evolutionary biology.

A limited role for gene duplications in the evolution of platypus venom

Gene duplication followed by adaptive selection is believed to be the primary driver of venom evolution. However, to date, no studies have evaluated the importance of gene duplications for venom evolution using a genomic approach. The availability of a sequenced genome and a venom gland transcriptome for the enigmatic platypus provides a unique opportunity to explore the role that gene duplication plays in venom evolution. Here, we identify gene duplication events and correlate them with expressed transcripts in an in-season venom gland. Gene duplicates (1,508) were identified. These duplicated pairs (421), including genes that have undergone multiple rounds of gene duplications, were expressed in the venom gland. The majority of these genes are involved in metabolism and protein synthesis not toxin functions. Twelve secretory genes including serine proteases, metalloproteinases, and protease inhibitors likely to produce symptoms of envenomation such as vasodilation and pain were detected. Only 16 of 107 platypus genes with high similarity to known toxins evolved through gene duplication. Platypus venom C-type natriuretic peptides and nerve growth factor do not possess lineage-specific gene duplicates. Extensive duplications, believed to increase the potency of toxic content and promote toxin diversification, were not found. This is the first study to take a genome-wide approach in order to examine the impact of gene duplication on venom evolution. Our findings support the idea that adaptive selection acts on gene duplicates to drive the independent evolution and functional diversification of similar venom genes in venomous species. However, gene duplications alone do not explain the "venome" of the platypus. Other mechanisms, such as alternative splicing and mutation, may be important in venom innovation.

Enzymatic toxins from snake venom: structural characterization and mechanism of catalysis

Snake venoms are cocktails of enzymes and non-enzymatic proteins used for both the immobilization and digestion of prey. The most common snake venom enzymes include acetylcholinesterases, l-amino acid oxidases, serine proteinases, metalloproteinases and phospholipases A(2) . Higher catalytic efficiency, thermal stability and resistance to proteolysis make these enzymes attractive models for biochemists, enzymologists and structural biologists. Here, we review the structures of these enzymes and describe their structure-based mechanisms of catalysis and inhibition. Some of the enzymes exist as protein complexes in the venom. Thus we also discuss the functional role of non-enzymatic subunits and the pharmacological effects of such protein complexes. The structures of inhibitor-enzyme complexes provide ideal platforms for the design of potent inhibitors which are useful in the development of prototypes and lead compounds with potential therapeutic applications.

Procoagulant adaptation of a blood coagulation prothrombinase-like enzyme complex in australian elapid venom

The macromolecular enzyme complex prothrombinase serves an indispensable role in blood coagulation as it catalyzes the conversion of prothrombin to thrombin, a key regulatory enzyme in the formation of a blood clot. Interestingly, a virtually identical enzyme complex is found in the venom of some Australian elapid snakes, which is composed of a cofactor factor Va-component and a serine protease factor Xa-like subunit. This review will provide an overview of the identification and characterization of the venom prothrombinase complex and will discuss the rationale for its powerful procoagulant nature responsible for the potent hemostatic toxicity of the elapid venom.

Inhibitory sequences within the B-domain stabilize circulating factor V in an inactive state

Blood coagulation factor V circulates as a procofactor with little or no procoagulant activity. It is activated to factor Va by thrombin following proteolytic removal of a large central B-domain. Although this reaction is well studied, the mechanism by which bond cleavage and B-domain release facilitate the transition to the active cofactor state has not been defined. Here we show that deletion or substitution of specific B-domain sequences drives the expression of procoagulant function without the need for proteolytic processing. Conversion to the constitutively active cofactor state is related, at least in part, to a cluster of amino acids that is highly basic and well conserved across the vertebrate lineage. Our findings demonstrate that discrete sequences in the B-domain serve to stabilize the inactive procofactor state, with proteolysis primarily functioning to remove these inhibitory constraints. These unexpected results provide new insight into the mechanism of factor V activation.

Structure of two genes encoding parallel prothrombin activators in Tropidechis carinatus snake: gene duplication and recruitment of factor X gene to the venom gland

BACKGROUND

Proteins with new function originate through gene duplication followed by divergence. In nature, occurrence of structurally and functionally similar proteins performing highly diverse physiological roles within an organism is rare. Several Australian elapid snakes have two parallel prothrombin activating systems with distinct physiological roles. For example, in Tropidechis carinatus, trocarin D exists in the venom and acts as toxin, whereas coagulation factor (F) X exists in plasma and plays a role in hemostasis.

RESULTS

Here, we show that FX and the trocarin D genes are expressed in a highly tissue-specific manner in T. carinatus. To understand their origin, recruitment and evolution, we determined the complete structure of their genes. Both genes have eight exons with identical exon-intron boundaries. All the introns are 92-99% identical with the exception of intron 1, indicating a recent gene duplication event. The first intron of the trocarin D gene is also nearly identical to that of the FX gene, except for two deletions (255 and 1406 bp) and three insertions (214, 1975, and 2174 bp). The third insertion has a potential scaffold/matrix attached region. The putative promoter of the snake FX gene shares similar cis-elements compared with those of human and murine FX genes. Interestingly, the trocarin D promoter has a 264-bp insertion with core promoter sequences and cis-elements that are known to induce high-level expression. This insertion might be responsible for switching from constitutive expression of the FX gene to inducible expression of trocarin D. Thus, we named this segment as VERSE (Venom Recruitment/Switch Element).

CONCLUSION

To our knowledge, this is the first molecular evidence for the recruitment of a duplicated gene for expression in venom glands by a simple insertion.

Molecular evolution caught in action: gene duplication and evolution of molecular isoforms of prothrombin activators in Pseudonaja textilis (brown snake)

BACKGROUND

The evolution of structurally and functionally similar proteins with highly diverse physiological roles within a single organism is of great interest. Australian elapid snakes offer an excellent opportunity to study the molecular evolution of prothrombin activators. Venom from Pseudonaja textilis contains pseutarin C, a group C prothrombin activator. Its enzymatic subunit is structurally and functionally similar to mammalian factor (F) Xa, whereas its non-enzymatic subunit is similar to FVa. As vertebrates, the snakes also contain a system to activate prothrombin in their own blood during injury. These hemostatic factors are produced in the liver.

RESULTS

Here we describe the presence of two molecular forms of FX expressed in the liver of P. textilis. Both isoforms have molecular signatures and domain architecture of FX. However, one isoform shows approximately 94% sequence identity with the snake FX from Tropidechis carinatus, whereas the other is much closer (90% identity) to the catalytic subunit of pseutarin C (PCCS). Real-time polymerase chain reaction reveals that the latter isoform is expressed approximately 56 000 times lower in the liver of P. textilis. However, the isoforms are not expressed in the venom gland.

CONCLUSION

A detailed analysis of deletions and insertions along with the sequence indicates that the second isoform is an intermediate caught in the evolution of venom prothrombin activator from the blood coagulation FX. Thus, this isoform represents a 'molecular fossil' and reveals the likely evolutionary path of recruitment of FX in the venom gland.