Robert A. Harrison

The king cobra genome reveals dynamic gene evolution and adaptation in the snake venom system

Significance Snake venoms are toxic protein cocktails used for prey capture. To investigate the evolution of these complex biological weapon systems, we sequenced the genome of a venomous snake, the king cobra, and assessed the composition of venom gland expressed genes, small RNAs, and secreted venom proteins. We show that regulatory components of the venom secretory system may have evolved from a pancreatic origin and that venom toxin genes were co-opted by distinct genomic mechanisms. After co-option, toxin genes important for prey capture have massively expanded by gene duplication and evolved under positive selection, resulting in protein neofunctionalization. This diverse and dramatic venom-related genomic response seemingly occurs in response to a coevolutionary arms race between venomous snakes and their prey. Snakes are limbless predators, and many species use venom to help overpower relatively large, agile prey. Snake venoms are complex protein mixtures encoded by several multilocus gene families that function synergistically to cause incapacitation. To examine venom evolution, we sequenced and interrogated the genome of a venomous snake, the king cobra (Ophiophagus hannah), and compared it, together with our unique transcriptome, microRNA, and proteome datasets from this species, with data from other vertebrates. In contrast to the platypus, the only other venomous vertebrate with a sequenced genome, we find that snake toxin genes evolve through several distinct co-option mechanisms and exhibit surprisingly variable levels of gene duplication and directional selection that correlate with their functional importance in prey capture. The enigmatic accessory venom gland shows a very different pattern of toxin gene expression from the main venom gland and seems to have recruited toxin-like lectin genes repeatedly for new nontoxic functions. In addition, tissue-specific microRNA analyses suggested the co-option of core genetic regulatory components of the venom secretory system from a pancreatic origin. Although the king cobra is limbless, we recovered coding sequences for all Hox genes involved in amniote limb development, with the exception of Hoxd12. Our results provide a unique view of the origin and evolution of snake venom and reveal multiple genome-level adaptive responses to natural selection in this complex biological weapon system. More generally, they provide insight into mechanisms of protein evolution under strong selection.

The Global Snake Bite Initiative: an antidote for snake bite

Clinicians have for a long time witnessed the tragedy of injury, disability, and death from snake bite that is a daily occurrence in many parts of Africa, Asia, and Latin America. To many people living in these regions, including some of the world’s poorest communities, snake bite is an ever present occupational risk and environmental hazard, an additional penalty of poverty. Like malaria, dengue, tuberculosis, and parasitic diseases, the risk of snake bite is always present. Unlike many of these other public health risks, however, the burden of human suff ering caused by snake bite remains un-recognised, invisible, and unheard by the global public health community, forgotten by development agencies and governments alike. The problem is so underrated that it was only added to WHO’s list of neglected tropical diseases in April, 2009.Yet an estimated 5·4–5·5 million people are bitten by snakes each year,

Snake envenoming: a disease of poverty.

Background Most epidemiological and clinical reports on snake envenoming focus on a single country and describe rural communities as being at greatest risk. Reports linking snakebite vulnerability to socioeconomic status are usually limited to anecdotal statements. The few reports with a global perspective have identified the tropical regions of Asia and Africa as suffering the highest levels of snakebite-induced mortality. Our analysis examined the association between globally available data on snakebite-induced mortality and socioeconomic indicators of poverty. Methodology/Principal Findings We acquired data on (i) the Human Development Index, (ii) the Per Capita Government Expenditure on Health, (iii) the Percentage Labour Force in Agriculture and (iv) Gross Domestic Product Per Capita from publicly available databases on the 138 countries for which snakebite-induced mortality rates have recently been estimated. The socioeconomic datasets were then plotted against the snakebite-induced mortality estimates (where both datasets were available) and the relationship determined. Each analysis illustrated a strong association between snakebite-induced mortality and poverty. Conclusions/Significance This study, the first of its kind, unequivocally demonstrates that snake envenoming is a disease of the poor. The negative association between snakebite deaths and government expenditure on health confirms that the burden of mortality is highest in those countries least able to deal with the considerable financial cost of snakebite.

Coevolution of diet and prey-specific venom activity supports the role of selection in snake venom evolution

The processes that drive the evolution of snake venom variability, particularly the role of diet, have been a topic of intense recent research interest. Here, we test whether extensive variation in venom composition in the medically important viper genus Echis is associated with shifts in diet. Examination of stomach and hindgut contents revealed extreme variation between the major clades of Echis in the proportion of arthropod prey consumed. The toxicity (median lethal dose, LD50) of representative Echis venoms to a natural scorpion prey species was found to be strongly associated with the degree of arthropod feeding. Mapping the results onto a novel Echis phylogeny generated from nuclear and mitochondrial sequence data revealed two independent instances of coevolution of venom toxicity and diet. Unlike venom LD50, the speed with which venoms incapacitated and killed scorpions was not associated with the degree of arthropod feeding. The prey-specific venom toxicity of arthropod-feeding Echis may thus be adaptive primarily by reducing venom expenditure. Overall, our results provide strong evidence that variation in snake venom composition results from adaptive evolution driven by natural selection for different diets, and underscores the need for a multi-faceted, integrative approach to the study of the causes of venom evolution.

Ending the drought: New strategies for improving the flow of affordable, effective antivenoms in Asia and Africa

The development of snake antivenoms more than a century ago should have heralded effective treatment of the scourge of snakebite envenoming in impoverished, mostly rural populations around the world. That snakebite still exists today, as a widely untreated illness that maims, kills and terrifies men, women and children in vulnerable communities, is a cruel anachronism. Antivenom can be an effective, safe and affordable treatment for snakebites, but apathy, inaction and the politicisation of public health have marginalised both the problem (making snakebite arguably the most neglected of all neglected tropical diseases) and its solution. For lack of any coordinated approach, provision of antivenoms has been pushed off the public health agenda, leading to an incongruous decline in demand for these crucial antidotes, excused and fed by new priorities, an absence of epidemiological data, and a poor regulatory framework. These factors facilitated the infiltration of poor quality products that degrade user confidence and undermine legitimate producers. The result is that tens of thousands are denied an essential life-saving medicine, allowing a toll of human suffering that is a summation of many individual catastrophes. No strategy has been developed to address this problem and to overcome the intransigence and inaction responsible for the global tragedy of snakebite. Attempts to engage with the broader public health community through the World Health Organisation (WHO), GAVI, and other agencies have failed. Consequently, the toxinology community has taken on a leadership role in a new approach, the Global Snakebite Initiative, which seeks to mobilise the resources, skills and experience of scientists and clinicians for whom venoms, toxins, antivenoms, snakes and snakebites are already fields of interest. Proteomics is one such discipline, which has embraced the potential of using venoms in bio-discovery and systems biology. The fields of venomics and antivenomics have recently evolved from this discipline, offering fresh hope for the victims of snakebites by providing an exciting insight into the complexities, nature, fundamental properties and significance of venom constituents. Such a rational approach brings with it the potential to design new immunising mixtures from which to raise potent antivenoms with wider therapeutic ranges. This addresses a major practical limitation in antivenom use recognised since the beginning of the 20th century: the restriction of therapeutic effectiveness to the specific venom immunogen used in production. Antivenomic techniques enable the interactions between venoms and antivenoms to be examined in detail, and if combined with functional assays of specific activity and followed up by clinical trials of effectiveness and safety, can be powerful tools with which to evaluate the suitability of current and new antivenoms for meeting urgent regional needs. We propose two mechanisms through which the Global Snakebite Initiative might seek to end the antivenom drought in Africa and Asia: first by establishing a multidisciplinary, multicentre, international collaboration to evaluate currently available antivenoms against the venoms of medically important snakes from specific nations in Africa and Asia using a combination of proteomic, antivenomic and WHO-endorsed preclinical assessment protocols, to provide a validated evidence base for either recommending or rejecting individual products; and secondly by bringing the power of proteomics to bear on the design of new immunising mixtures to raise Pan-African and Pan-Asian polyvalent antivenoms of improved potency and quality. These products will be subject to rigorous clinical assessment. We propose radically to change the basis upon which antivenoms are produced and supplied for the developing world. Donor funding and strategic public health alliances will be sought to make it possible not only to sustain the financial viability of antivenom production partnerships, but also to ensure that patients are relieved of the costs of antivenom so that poverty is no longer a barrier to the treatment of this important, but grossly neglected public health emergency.

Medically important differences in snake venom composition are dictated by distinct postgenomic mechanisms

Significance The toxic composition of snake venom varies between species. Such variation can have major medical implications for the treatment of human snakebite victims. Venom variation is largely attributed to differences in toxin-encoding genes present in the genome or venom gland of snakes. Here, we demonstrate that mechanisms affecting the transcription, translation, and posttranslational modification of toxins also significantly contribute to the diversity of venom protein composition. Venom variation observed between related snake species is therefore the result of a complex interaction between a variety of genetic and postgenomic factors acting on toxin genes. Ultimately, this variation results in significant differences in venom-induced pathology and lethality and can undermine the efficacy of antivenom therapies used to treat human snakebite victims. Variation in venom composition is a ubiquitous phenomenon in snakes and occurs both interspecifically and intraspecifically. Venom variation can have severe outcomes for snakebite victims by rendering the specific antibodies found in antivenoms ineffective against heterologous toxins found in different venoms. The rapid evolutionary expansion of different toxin-encoding gene families in different snake lineages is widely perceived as the main cause of venom variation. However, this view is simplistic and disregards the understudied influence that processes acting on gene transcription and translation may have on the production of the venom proteome. Here, we assess the venom composition of six related viperid snakes and compare interspecific changes in the number of toxin genes, their transcription in the venom gland, and their translation into proteins secreted in venom. Our results reveal that multiple levels of regulation are responsible for generating variation in venom composition between related snake species. We demonstrate that differential levels of toxin transcription, translation, and their posttranslational modification have a substantial impact upon the resulting venom protein mixture. Notably, these processes act to varying extents on different toxin paralogs found in different snakes and are therefore likely to be as important as ancestral gene duplication events for generating compositionally distinct venom proteomes. Our results suggest that these processes may also contribute to altering the toxicity of snake venoms, and we demonstrate how this variability can undermine the treatment of a neglected tropical disease, snakebite.

Comparative venom gland transcriptome surveys of the saw-scaled vipers (Viperidae: Echis) reveal substantial intra-family gene diversity and novel venom transcripts

BackgroundVenom variation occurs at all taxonomical levels and can impact significantly upon the clinical manifestations and efficacy of antivenom therapy following snakebite. Variation in snake venom composition is thought to be subject to strong natural selection as a result of adaptation towards specific diets. Members of the medically important genus Echis exhibit considerable variation in venom composition, which has been demonstrated to co-evolve with evolutionary shifts in diet. We adopt a venom gland transcriptome approach in order to investigate the diversity of toxins in the genus and elucidate the mechanisms which result in prey-specific adaptations of venom composition.ResultsVenom gland transcriptomes were created for E. pyramidum leakeyi, E. coloratus and E. carinatus sochureki by sequencing ~1000 expressed sequence tags from venom gland cDNA libraries. A standardised methodology allowed a comprehensive intra-genus comparison of the venom gland profiles to be undertaken, including the previously described E. ocellatus transcriptome. Blast annotation revealed the presence of snake venom metalloproteinases, C-type lectins, group II phopholipases A2, serine proteases, L-amino oxidases and growth factors in all transcriptomes throughout the genus. Transcripts encoding disintegrins, cysteine-rich secretory proteins and hyaluronidases were obtained from at least one, but not all, species. A representative group of novel venom transcripts exhibiting similarity to lysosomal acid lipase were identified from the E. coloratus transcriptome, whilst novel metallopeptidases exhibiting similarity to neprilysin and dipeptidyl peptidase III were identified from E. p. leakeyi and E. coloratus respectively.ConclusionThe comparison of Echis venom gland transcriptomes revealed substantial intrageneric venom variation in representations and cluster numbers of the most abundant venom toxin families. The expression profiles of established toxin groups exhibit little obvious association with venom-related adaptations to diet described from this genus. We suggest therefore that alterations in isoform diversity or transcript expression levels within the major venom protein families are likely to be responsible for prey specificity, rather than differences in the representation of entire toxin families or the recruitment of novel toxin families, although the recruitment of lysosomal acid lipase as a response to vertebrate feeding cannot be excluded. Evidence of marked intrageneric venom variation within the medically important genus Echis strongly advocates further investigations into the medical significance of venom variation in this genus and its impact upon antivenom therapy.

Snake Venomics of African Spitting Cobras: Toxin Composition and Assessment of Congeneric Cross-Reactivity of the Pan-African EchiTAb-Plus-ICP Antivenom by Antivenomics and Neutralization Approaches

Venomic analysis of the venoms of Naja nigricollis, N. katiensis, N. nubiae, N. mossambica, and N. pallida revealed similar compositional trends. The high content of cytotoxins and PLA(2)s may account for the extensive tissue necrosis characteristic of the envenomings by these species. The high abundance of a type I α-neurotoxin in N. nubiae may be responsible for the high lethal toxicity of this venom (in rodents). The ability of EchiTAb-Plus-ICP antivenom to immunodeplete and neutralize the venoms of African spitting cobras was assessed by antivenomics and neutralization tests. It partially immunodepleted 3FTx and PLA(2)s and completely immunodepleted SVMPs and CRISPs in all venoms. The antivenom neutralized the dermonecrotic and PLA(2) activities of all African Naja venoms, whereas lethality was eliminated in the venoms of N. nigricollis, N. mossambica, and N. pallida but not in those of N. nubiae and N. katiensis. The lack of neutralization of lethality of N. nubiae venom may be of medical relevance only in relatively populous areas of the Saharan region. The impaired activity of EchiTAb-Plus-ICP against N. katiensis may not represent a major concern. This species is sympatric with N. nigricollis in many regions of Africa, although very few bites have been attributed to it.

Domain Loss Facilitates Accelerated Evolution and Neofunctionalization of Duplicate Snake Venom Metalloproteinase Toxin Genes

Gene duplication is a key mechanism for the adaptive evolution and neofunctionalization of gene families. Large multigene families often exhibit complex evolutionary histories as a result of frequent gene duplication acting in concordance with positive selection pressures. Alterations in the domain structure of genes, causing changes in the molecular scaffold of proteins, can also result in a complex evolutionary history and has been observed in functionally diverse multigene toxin families. Here, we investigate the role alterations in domain structure have on the tempo of evolution and neofunctionalization of multigene families using the snake venom metalloproteinases (SVMPs) as a model system. Our results reveal that the evolutionary history of viperid (Serpentes: Viperidae) SVMPs is repeatedly punctuated by domain loss, with the single loss of the cysteine-rich domain, facilitating the formation of P-II class SVMPs, occurring prior to the convergent loss of the disintegrin domain to form multiple P-I SVMP structures. Notably, the majority of phylogenetic branches where domain loss was inferred to have occurred exhibited highly significant evidence of positive selection in surface-exposed amino acid residues, resulting in the neofunctionalization of P-II and P-I SVMP classes. These results provide a valuable insight into the mechanisms by which complex gene families evolve and detail how the loss of domain structures can catalyze the accelerated evolution of novel gene paralogues. The ensuing generation of differing molecular scaffolds encoded by the same multigene family facilitates gene neofunctionalization while presenting an evolutionary advantage through the retention of multiple genes capable of encoding functionally distinct proteins.

Combined snake venomics and venom gland transcriptomic analysis of the ocellated carpet viper, Echis ocellatus

Snakebite in Africa causes thousands of deaths annually and considerable permanent physical disability. The saw-scaled viper, Echis ocellatus, represents the single most medically important snake species in West Africa. To provide a detailed compositional analysis of the venom of E. ocellatus for designing novel toxin-specific immunotherapy and to delineate sequence structure-function relationships of individual toxins, we characterised the venom proteome and the venom gland transcriptome. Whole E. ocellatus venom was fractionated by reverse-phase HPLC, followed by analysis of each chromatographic fraction using a combination of SDS-PAGE, N-terminal sequencing, MALDI-TOF mass fingerprinting, and CID-MS/MS of tryptic peptides. This analysis identified around 35 distinct proteins of molecular masses in the range of 5.5-110 kDa belonging to 8 different toxin families (disintegrin, DC-fragment, phospholipase A(2), cysteine-rich secretory protein, serine proteinase, C-type lectin, l-amino acid oxidase, and Zn(2+)-dependent metalloprotease). Comparison of the toxin composition of E. ocellatus venom determined using a proteomic approach, with the predicted proteome derived from assembly of 1000 EST sequences from a E. ocellatus venom gland cDNA library, shows some differences. Most notably, peptides derived from 26% of the venom proteins could not be ascribed an exact match in the transcriptome. Similarly, 64 (67%) out of the 95 putative toxin clusters reported in the transcriptome did not match to peptides detected in the venom proteome. These data suggest that the final composition of venom is influenced by transcriptional and post-translational mechanisms that may be more complex than previously appreciated. This, in turn, emphasises the value of combining proteomic and transcriptomic approaches to acquire a more complete understanding of the precise composition of snake venom, than would be gleaned from using one analysis alone. From a clinical perspective, the large amount of SVMPs (66.5% of the total venom proteins) is consistent with the haemorrhagic pathology associated with E. ocellatus envenoming. More significantly, whilst the proteomic analysis confirms the majority of these metalloproteinases (58%) belong to the SVMP PIII class, MS/MS derived peptide sequencing also demonstrates a major constituent (32%) of E. ocellatus venom is a PIV-SVMP with a quaternary structure comprising a 48 kDa (Q2UXQ4 or Q2UXQ5) PIII-SVMP subunit, and two 14-16 kDa C-type lectin-like domains [EOC_00087 and EOC_00124] which display similarity to echicetin alpha [P81017] and beta [P81996] subunits.