Wolbachia depletion blocks transmission of lymphatic filariasis by preventing chitinase-dependent parasite exsheathment.

Lymphatic filariasis is a vector-borne neglected tropical disease prioritized for global elimination. The filarial nematodes that cause the disease host a symbiotic bacterium, Wolbachia, which has been targeted using antibiotics, leading to cessation of parasite embryogenesis, waning of circulating larvae (microfilariae [mf]), and gradual cure of adult infection. One of the benefits of the anti-Wolbachia mode of action is that it avoids the rapid killing of mf, which can drive inflammatory adverse events. However, mf depleted of Wolbachia persist for several months in circulation, and thus patients treated with antibiotics are assumed to remain at risk for transmitting infections. Here, we show that Wolbachia-depleted mf rapidly lose the capacity to develop in the mosquito vector through a defect in exsheathment and inability to migrate through the gut wall. Transcriptomic and Western blotting analyses demonstrate that chitinase, an enzyme essential for mf exsheathment, is down-regulated in Wolbachia-depleted mf and correlates with their inability to exsheath and escape the mosquito midgut. Supplementation of in vitro cultures of Wolbachia-depleted mf with chitinase enzymes restores their ability to exsheath to a similar level to that observed in untreated mf. Our findings elucidate a mechanism of rapid transmission-blocking activity of filariasis after depletion of Wolbachia and adds to the broad range of biological processes of filarial nematodes that are dependent on Wolbachia symbiosis.

Comparative 1D Blue-Native electrophoresis analysis of Plasmodium falciparum and human proteins associated with cytoadherence.

BACKGROUND: To understand more about changes to the molecular components that occur when host endothelium interacts with Plasmodium falciparum-infected erythrocytes, a combined technique of protein separation (1D Blue-Native electrophoresis) and mass spectrometry of infected erythrocytes with endothelial cells (EC) in a co-culture system has been used. METHODS: Native proteins were extracted from co-cultures and identified by mass spectrometry. Proteomic data from different parasite strains, either adhesion proficient (to endothelial cells) or non-adherent, were analysed in parallel to reveal protein associations linked to cytoadherence. Informatic approaches were developed to facilitate this comparison. RESULTS: Blue-Native gel separation and LC/MS/MS identification revealed major differences in samples produced from endothelial cell co-culture with adherent and non-adherent parasite strains. This approach enabled us to identify protein associations seen only with the adhesion proficient parasite strain. CONCLUSIONS: The combination of proteomic and analytical approaches has identified differences between adherent and non-adherent parasite lines in co-culture with EC, providing potential candidates for complexes or associations formed during cytoadherence involved in cell structure, signalling and apoptosis.

What killed Karl Patterson Schmidt? Combined venom gland transcriptomic, venomic and antivenomic analysis of the South African green tree snake (the boomslang), Dispholidus typus.

BACKGROUND: Non-front-fanged colubroid snakes comprise about two-thirds of extant ophidian species. The medical significance of the majority of these snakes is unknown, but at least five species have caused life-threatening or fatal human envenomings. However, the venoms of only a small number of species have been explored. METHODS: A combined venomic and venom gland transcriptomic approach was employed to characterise of venom of Dispholidus typus (boomslang), the snake that caused the tragic death of Professor Karl Patterson Schmidt. The ability of CroFab™ antivenom to immunocapture boomslang venom proteins was investigated using antivenomics. RESULTS: Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families: three-finger toxin (3FTx); phospholipase A2 (PLA2); cysteine-rich secretory proteins (CRISP); snake venom (SV) serine proteinase (SP); C-type lectin-like (CTL); SV metalloproteinases (SVMPs); and disintegrin-like/cysteine-rich (DC) proteolytic fragments. CroFab™ antivenom efficiently immunodepleted some boomslang SVMPs. CONCLUSIONS: The present work is the first to address the overall proteomic profile of D. typus venom. This study allowed us to correlate the toxin composition with the toxic activities of the venom. The antivenomic analysis suggested that the antivenom available at the time of the unfortunate accident could have exhibited at least some immunoreactivity against the boomslang SVMPs responsible for the disseminated intravascular coagulation syndrome that caused K.P. Schmidt's fatal outcome. GENERAL SIGNIFICANCE: This study may stimulate further research on other non-front-fanged colubroid snake venoms capable of causing life-threatening envenomings to humans, which in turn should contribute to prevent fatal human accidents, such as that unfortunately suffered by K.P. Schmidt.

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

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.

VTBuilder: a tool for the assembly of multi isoform transcriptomes.

BACKGROUND: Within many research areas, such as transcriptomics, the millions of short DNA fragments (reads) produced by current sequencing platforms need to be assembled into transcript sequences before they can be utilized. Despite recent advances in assembly software, creating such transcripts from read data harboring isoform variation remains challenging. This is because current approaches fail to identify all variants present or they create chimeric transcripts within which relationships between co-evolving sites and other evolutionary factors are disrupted. We present VTBuilder, a tool for constructing non-chimeric transcripts from read data that has been sequenced from sources containing isoform complexity. RESULTS: We validated VTBuilder using reads simulated from 54 Sanger sequenced transcripts (SSTs) expressed in the venom gland of the saw scaled viper, Echis ocellatus. The SSTs were selected to represent genes from major co-expressed toxin groups known to harbor isoform variants. From the simulated reads, VTBuilder constructed 55 transcripts, 50 of which had a greater than 99% sequence similarity to 48 of the SSTs. In contrast, using the popular assembler tool Trinity (r2013-02-25), only 14 transcripts were constructed with a similar level of sequence identity to just 11 SSTs. Furthermore VTBuilder produced transcripts with a similar length distribution to the SSTs while those produced by Trinity were considerably shorter. To demonstrate that our approach can be scaled to real world data we assembled the venom gland transcriptome of the African puff adder Bitis arietans using paired-end reads sequenced on Illumina's MiSeq platform. VTBuilder constructed 1481 transcripts from 5 million reads and, following annotation, all major toxin genes were recovered demonstrating reconstruction of complex underlying sequence and isoform diversity. CONCLUSION: Unlike other approaches, VTBuilder strives to maintain the relationships between co-evolving sites within the constructed transcripts, and thus increases transcript utility for a wide range of research areas ranging from transcriptomics to phylogenetics and including the monitoring of drug resistant parasite populations. Additionally, improving the quality of transcripts assembled from read data will have an impact on future studies that query these data. VTBuilder has been implemented in java and is available, under the GPL GPU V0.3 license, from http:// http://www.lstmed.ac.uk/vtbuilder .

Dissecting the organ specificity of insecticide resistance candidate genes in Anopheles gambiae: known and novel candidate genes.

BACKGROUND: The elevated expression of enzymes with insecticide metabolism activity can lead to high levels of insecticide resistance in the malaria vector, Anopheles gambiae. In this study, adult female mosquitoes from an insecticide susceptible and resistant strain were dissected into four different body parts. RNA from each of these samples was used in microarray analysis to determine the enrichment patterns of the key detoxification gene families within the mosquito and to identify additional candidate insecticide resistance genes that may have been overlooked in previous experiments on whole organisms. RESULTS: A general enrichment in the transcription of genes from the four major detoxification gene families (carboxylesterases, glutathione transferases, UDP glucornyltransferases and cytochrome P450s) was observed in the midgut and malpighian tubules. Yet the subset of P450 genes that have previously been implicated in insecticide resistance in An gambiae, show a surprisingly varied profile of tissue enrichment, confirmed by qPCR and, for three candidates, by immunostaining. A stringent selection process was used to define a list of 105 genes that are significantly (p ≤0.001) over expressed in body parts from the resistant versus susceptible strain. Over half of these, including all the cytochrome P450s on this list, were identified in previous whole organism comparisons between the strains, but several new candidates were detected, notably from comparisons of the transcriptomes from dissected abdomen integuments. CONCLUSIONS: The use of RNA extracted from the whole organism to identify candidate insecticide resistance genes has a risk of missing candidates if key genes responsible for the phenotype have restricted expression within the body and/or are over expression only in certain tissues. However, as transcription of genes implicated in metabolic resistance to insecticides is not enriched in any one single organ, comparison of the transcriptome of individual dissected body parts cannot be recommended as a preferred means to identify new candidate insecticide resistant genes. Instead the rich data set on in vivo sites of transcription should be consulted when designing follow up qPCR validation steps, or for screening known candidates in field populations.

Gene tree parsimony of multilocus snake venom protein families reveals species tree conflict as a result of multiple parallel gene loss.

The proliferation of gene data from multiple loci of large multigene families has been greatly facilitated by considerable recent advances in sequence generation. The evolution of such gene families, which often undergo complex histories and different rates of change, combined with increases in sequence data, pose complex problems for traditional phylogenetic analyses, and in particular, those that aim to successfully recover species relationships from gene trees. Here, we implement gene tree parsimony analyses on multicopy gene family data sets of snake venom proteins for two separate groups of taxa, incorporating Bayesian posterior distributions as a rigorous strategy to account for the uncertainty present in gene trees. Gene tree parsimony largely failed to infer species trees congruent with each other or with species phylogenies derived from mitochondrial and single-copy nuclear sequences. Analysis of four toxin gene families from a large expressed sequence tag data set from the viper genus Echis failed to produce a consistent topology, and reanalysis of a previously published gene tree parsimony data set, from the family Elapidae, suggested that species tree topologies were predominantly unsupported. We suggest that gene tree parsimony failure in the family Elapidae is likely the result of unequal and/or incomplete sampling of paralogous genes and demonstrate that multiple parallel gene losses are likely responsible for the significant species tree conflict observed in the genus Echis. These results highlight the potential for gene tree parsimony analyses to be undermined by rapidly evolving multilocus gene families under strong natural selection.

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.

Stable high-density and maternally inherited Wolbachia infections in Anopheles moucheti and Anopheles demeilloni mosquitoes.

Wolbachia, a widespread bacterium that can reduce pathogen transmission in mosquitoes, has recently been reported to be present in Anopheles (An.) species. In wild populations of the An. gambiae complex, the primary vectors of Plasmodium malaria in Sub-Saharan Africa, Wolbachia DNA sequences at low density and infection frequencies have been detected. As the majority of studies have used highly sensitive nested PCR as the only method of detection, more robust evidence is required to determine whether Wolbachia strains are established as endosymbionts in Anopheles species. Here, we describe high-density Wolbachia infections in geographically diverse populations of An. moucheti and An. demeilloni. Fluorescent in situ hybridization localized a heavy infection in the ovaries of An. moucheti, and maternal transmission was observed. Genome sequencing of both Wolbachia strains obtained genome depths and coverages comparable to those of other known infections. Notably, homologs of cytoplasmic incompatibility factor (cif) genes were present, indicating that these strains possess the capacity to induce the cytoplasmic incompatibility phenotype, which allows Wolbachia to spread through host populations. These strains should be further investigated as candidates for use in Wolbachia biocontrol strategies in Anopheles aiming to reduce the transmission of malaria.

IR-TEx: An Open Source Data Integration Tool for Big Data Transcriptomics Designed for the Malaria Vector Anopheles gambiae.

IR-TEx is an application written in Shiny (an R package) that allows exploration of the expression of (as well as assigning functions to) transcripts whose expression is associated with insecticide resistance phenotypes in Anopheles gambiae mosquitoes. The application can be used online or downloaded and used locally by anyone. The local application can be modified to add new insecticide resistance datasets generated from multiple -omics platforms. This guide demonstrates how to add new datasets and handle missing data. Furthermore, IR-TEx can be completely and easily recoded to use-omics datasets from any experimental data, making it a valuable resource to many researchers. The protocol illustrates the utility of IR-TEx in identifying new insecticide resistance candidates using the the microsomal glutathione transferase, GSTMS1, as an example. This transcript is upregulated in multiple pyrethroid resistant populations from Côte D'Ivoire and Burkina Faso. The identification of co-correlated transcripts provides further insight into the putative roles of this gene.

The insufficiency of circulating miRNA and DNA as diagnostic tools or as biomarkers of treatment efficacy for Onchocerca volvulus.

Skin snip evaluation for onchocerciasis has insufficient sensitivity when skin microfilarial (mf) densities are low, such as following ivermectin treatment. Mf density is suitable for assessing microfilaricidal efficacy but only serves as an indirect indicator of macrofilaricidal activity. We assessed circulating nucleic acids from Onchocerca volvulus as an alternative to skin snips. We screened a plasma sample set of infected individuals followed up at four, 12 and 21 months after microfilaricidal (ivermectin, n = four), macrofilaricidal (doxycycline, n = nine), or combination treatment (n = five). Two parasite-derived miRNAs, cel-miR-71-5p and bma-lin-4, and O-150 repeat DNA were assessed. Highly abundant DNA repeat families identified in the O. volvulus genome were also evaluated. miRNAs were detected in two of 72 plasma samples (2.8%) and two of 47 samples (4.3%) with microfilaridermia using RT-qPCR. O-150 DNA was detected in eight (44.4%) baseline samples by qPCR and the number of positives declined post-treatment. One doxycycline-treated individual remained O-150 positive. However, only 11 (23.4%) samples with microfilaridermia were qPCR-positive. Analysis by qPCR showed novel DNA repeat families were comparatively less abundant than the O-150 repeat. Circulating parasite-derived nucleic acids are therefore insufficient as diagnostic tools or as biomarkers of treatment efficacy for O. volvulus.

Venomics of the asp viper Vipera aspis aspis from France.

The asp viper Vipera aspis aspis is a venomous snake found in France, and despite its medical importance, the complete toxin repertoire produced is unknown. Here, we used a venomics approach to decipher the composition of its venom. Transcriptomic analysis revealed 80 venom-annotated sequences grouped into 16 gene families. Among the most represented toxins were snake venom metalloproteases (23%), phospholipases A2 (15%), serine proteases (13%), snake venom metalloprotease inhibitors (13%) and C-type lectins (12%). LC-MS of venoms revealed similar profiles regardless of the method of extraction (milking vs defensive bite). Proteomic analysis validated 57 venom-annotated transcriptomic sequences (>70%), including one for each of the 16 families, but also identified 7 sequences not initially annotated as venom proteins, including a serine protease, a disintegrin, a glutaminyl-peptide cyclotransferase, a proactivator polypeptide-like and 3 aminopeptidases. Interestingly, phospholipases A2 were the dominant proteins in the venom, among which included an ammodytoxin B-like sequence, which may explain the reported neurotoxicity following some asp viper envenomations. In total, 87 sequences were retrieved from the Vipera aspis aspis transcriptome and proteome, constituting a valuable resource that will help in understanding the toxinological basis of clinical signs of envenoming and for the mining of useful pharmacological compounds. BIOLOGICAL SIGNIFICANCE: The asp viper (Vipera aspis aspis) causes several hundred envenomations annually in France, including unusual cases with neurological signs, resulting in one death per year on average. Here, we performed a proteotranscriptomic analysis of V. a. aspis venom in order to provide a better understanding of its venom composition. We found that, as in other Vipera species, phospholipase A2 dominates in the venom, and the presence of a sequence related to ammodytoxin B may explain the reported neurotoxicity following some asp viper envenomations. Thus, this study will help in informing the toxinological basis of clinical signs of envenoming.

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

BACKGROUND: Venom 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. RESULTS: Venom gland transcriptomes were created for E. pyramidum leakeyi, E. coloratus and E. carinatus sochureki by sequencing approximately 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. CONCLUSION: The 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.

A Decoy-Receptor Approach Using Nicotinic Acetylcholine Receptor Mimics Reveals Their Potential as Novel Therapeutics Against Neurotoxic Snakebite.

Snakebite is a neglected tropical disease that causes 138,000 deaths each year. Neurotoxic snake venoms contain small neurotoxins, including three-finger toxins (3FTxs), which can cause rapid paralysis in snakebite victims by blocking postsynaptic transmission via nicotinic acetylcholine receptors (nAChRs). These toxins are typically weakly immunogenic and thus are often not effectively targeted by current polyclonal antivenom therapies. We investigated whether nAChR mimics, also known as acetylcholine binding proteins (AChBPs), could effectively capture 3FTxs and therefore be developed as a novel class of snake-generic therapeutics for combatting neurotoxic envenoming. First, we identified the binding specificities of 3FTx from various medically important elapid snake venoms to nAChR using two recombinant nAChR mimics: the AChBP from Lymnaea stagnalis and a humanized neuronal α7 version (α7-AChBP). We next characterized these AChBP-bound and unbound fractions using SDS-PAGE and mass spectrometry. Interestingly, both mimics effectively captured long-chain 3FTxs from multiple snake species but largely failed to capture the highly related short-chain 3FTxs, suggesting a high level of binding specificity. We next investigated whether nAChR mimics could be used as snakebite therapeutics. We showed that while α7-AChBP alone did not protect against Naja haje (Egyptian cobra) venom lethality in vivo, it significantly prolonged survival times when coadministered with a nonprotective dose of antivenom. Thus, nAChR mimics are capable of neutralizing specific venom toxins and may be useful adjunct therapeutics for improving the safety and affordability of existing snakebite treatments by reducing therapeutic doses. Our findings justify exploring the future development of AChBPs as potential snakebite treatments.

Defining the pathogenic threat of envenoming by South African shield-nosed and coral snakes (genus Aspidelaps), and revealing the likely efficacy of available antivenom.

While envenoming by the southern African shield-nosed or coral snakes (genus Aspidelaps) has caused fatalities, bites are uncommon. Consequently, this venom is not used in the mixture of snake venoms used to immunise horses for the manufacture of regional SAIMR (South African Institute for Medical Research) polyvalent antivenom. Aspidelaps species are even excluded from the manufacturer's list of venomous snakes that can be treated by this highly effective product. This leaves clinicians, albeit rarely, in a therapeutic vacuum when treating envenoming by these snakes. This is a significantly understudied small group of nocturnal snakes and little is known about their venom compositions and toxicities. Using a murine preclinical model, this study determined that the paralysing toxicity of venoms from Aspidelaps scutatus intermedius, A. lubricus cowlesi and A. l. lubricus approached that of venoms from highly neurotoxic African cobras and mambas. This finding was consistent with the cross-genus dominance of venom three-finger toxins, including numerous isoforms which showed extensive interspecific variation. Our comprehensive analysis of venom proteomes showed that the three Aspidelaps species possess highly similar venom proteomic compositions. We also revealed that the SAIMR polyvalent antivenom cross-reacted extensively in vitro with venom proteins of the three Aspidelaps. Importantly, this cross-genus venom-IgG binding translated to preclinical (in a murine model) neutralisation of A. s. intermedius venom-induced lethality by the SAIMR polyvalent antivenom, at doses comparable with those that neutralise venom from the cape cobra (Naja nivea), which the antivenom is directed against. Our results suggest a wider than anticipated clinical utility of the SAIMR polyvalent antivenom, and here we seek to inform southern African clinicians that this readily available antivenom is likely to prove effective for victims of Aspidelaps envenoming. BIOLOGICAL SIGNIFICANCE: Coral and shield-nosed snakes (genus Aspidelaps) comprise two species and several subspecies of potentially medically important venomous snakes distributed in Namibia, Botswana, Zimbabwe, Mozambique and South Africa. Documented human fatalities, although rare, have occurred from both A. lubricus and A. scutatus. However, their venom proteomes and the pathological effects of envenomings by this understudied group of nocturnal snakes remain uncharacterised. Furthermore, no commercial antivenom is made using venom from species of the genus Aspidelaps. To fill this gap, we have conducted a transcriptomics-guided comparative proteomics analysis of the venoms of the intermediate shield-nose snake (A. s. intermedius), southern coral snake (A. l. lubricus), and Cowle's shield snake (A. l. cowlesi); investigated the mechanism of action underpinning lethality by A. s. intermedius in the murine model; and assessed the in vitro immunoreactivity of the SAIMR polyvalent antivenom towards the venom toxins of A. l. lubricus and A. l. cowlesi, and the in vivo capability of this antivenom at neutralising the lethal effect of A. s. intermedius venom. Our data revealed a high degree of conservation of the global composition of the three Aspidelaps venom proteomes, all characterised by the overwhelming predominance of neurotoxic 3FTxs, which induced classical signs of systemic neurotoxicity in mice. The SAIMR polyvalent antivenom extensively binds to Aspidelaps venom toxins and neutralised, with a potency of 0.235 mg venom/mL antivenom, the lethal effect of A. s. intermedius venom. Our data suggest that the SAIMR antivenom could be a useful therapeutic tool for treating human envenomings by Aspidelaps species.

Intact protein mass spectrometry reveals intraspecies variations in venom composition of a local population of Vipera kaznakovi in Northeastern Turkey.

We report on the variable venom composition of a population of the Caucasus viper (Vipera kaznakovi) in Northeastern Turkey. We applied a combination of venom gland transcriptomics, de-complexing bottom-up and top-down venomics. In contrast to sole bottom-up venomics approaches and gel or chromatography based venom comparison, our combined approach enables a faster and more detailed comparison of venom proteomes from multiple individuals. In total, we identified peptides and proteins from 15 toxin families, including snake venom metalloproteinases (svMP; 37.8%), phospholipases A2 (PLA2; 19.0%), snake venom serine proteinases (svSP; 11.5%), C-type lectins (CTL; 6.9%) and cysteine-rich secretory proteins (CRISP; 5.0%), in addition to several low abundant toxin families. Furthermore, we identified intraspecies variations of the venom composition of V. kaznakovi, and find these were mainly driven by the age of the animals, with lower svSP abundance detected in juveniles. On the proteoform level, several small molecular weight toxins between 5 and 8 kDa in size, as well as PLA2s, drove the differences observed between juvenile and adult individuals. This study provides novel insights into the venom variability of V. kaznakovi and highlights the utility of intact mass profiling for fast and detailed comparison of snake venom. BIOLOGICAL SIGNIFICANCE: Population level and ontogenetic venom variation (e.g. diet, habitat, sex or age) can result in a loss of antivenom efficacy against snakebites from wide ranging snake populations. The current state of the art for the analysis of snake venoms are de-complexing bottom-up proteomics approaches. While useful, these have the significant drawback of being time-consuming and following costly protocols, and consequently are often applied to pooled venom samples. To overcome these shortcomings and to enable rapid and detailed profiling of large numbers of individual venom samples, we integrated an intact protein analysis workflow into a transcriptomics-guided bottom-up approach. The application of this workflow to snake individuals of a local population of V. kaznakovi revealed intraspecies variations in venom composition, which are primarily explained by the age of the animals, and highlighted svSP abundance to be one of the molecular drivers for the compositional differences observed.

Molecular characterisation of endogenous snake venom metalloproteinase inhibitors.

Viper venoms contain one of the most potent mixtures of proteases in natural existence and yet the venom gland and proteins in this mixture are refractory to degradation. Here we demonstrate that the sub-10-kDa components of venom from two African viper species (Echis ocellatus and Cerastes cerastes cerastes) are predominantly composed of the tri-peptide pyroglutamate-lysine-tryptophan (pEKW). This tripeptide is encoded by tandemly repeating elements and, in E. ocellatus, on the same transcript as a novel, and highly unusual, poly-histidine-poly-glycine peptide (pHpG) also detected in E. ocellatus venom. The pEKW and pHpG peptides inhibit the proteolytic activity of the haemorrhagic snake venom metalloproteinase (SVMP), EoVMP-2, and the haemorrhagic activity of E. ocellatus venom. These results demonstrate that these vipers express abundant transcripts encoding tandemly repeated protease inhibitor cassettes and accumulate significant quantities of peptide inhibitors in venoms to provide a basis for attenuating the proteolytic activity of SVMPs.

The medical threat of mamba envenoming in sub-Saharan Africa revealed by genus-wide analysis of venom composition, toxicity and antivenomics profiling of available antivenoms.

Mambas (genus Dendroaspis) are among the most feared of the medically important elapid snakes found in sub-Saharan Africa, but many facets of their biology, including the diversity of venom composition, remain relatively understudied. Here, we present a reconstruction of mamba phylogeny, alongside genus-wide venom gland transcriptomic and high-resolution top-down venomic analyses. Whereas the green mambas, D. viridis, D. angusticeps, D. j. jamesoni and D. j. kaimosae, express 3FTx-predominant venoms, black mamba (D. polylepis) venom is dominated by dendrotoxins I and K. The divergent terrestrial ecology of D. polylepis compared to the arboreal niche occupied by all other mambas makes it plausible that this major difference in venom composition is due to dietary variation. The pattern of intrageneric venom variability across Dendroaspis represented a valuable opportunity to investigate, in a genus-wide context, the variant toxicity of the venom, and the degree of paraspecific cross-reactivity between antivenoms and mamba venoms. To this end, the immunological profiles of the five mamba venoms were assessed against a panel of commercial antivenoms generated for the sub-Saharan Africa market. This study provides a genus-wide overview of which available antivenoms may be more efficacious in neutralising human envenomings caused by mambas, irrespective of the species responsible. The information gathered in this study lays the foundations for rationalising the notably different potency and pharmacological profiles of Dendroaspis venoms at locus resolution. This understanding will allow selection and design of toxin immunogens with a view to generating a safer and more efficacious pan-specific antivenom against any mamba envenomation. BIOLOGICAL SIGNIFICANCE: The mambas (genus Dendroaspis) comprise five especially notorious medically important venomous snakes endemic to sub-Saharan Africa. Their highly potent venoms comprise a high diversity of pharmacologically active peptides, including extremely rapid-acting neurotoxins. Previous studies on mamba venoms have focused on the biochemical and pharmacological characterisation of their most relevant toxins to rationalize the common neurological and neuromuscular symptoms of envenomings caused by these species, but there has been little work on overall venom composition or comparisons between them. Only very recently an overview of the composition of the venom of two Dendroaspis species, D. angusticeps and D. polylepis, has been unveiled through venomics approaches. Here we present the first genus-wide transcriptomic-proteomic analysis of mamba venom composition. The transcriptomic analyses described in this paper have contributed 29 (D. polylepis), 23 (D. angusticeps), 40 (D. viridis), 25 (D. j. jamesoni) and 21 (D. j. kaimosae), novel full-length toxin sequences to the non-redundant Dendroaspis sequence database. The mamba genus-wide venomic analysis demonstrated that major D. polylepis venom components are Kunitz-fold family toxins. This feature is unique in relation to the relatively conserved three-finger toxin (3FTx)-dominated venom compositions of the green mambas. Venom variation was interpreted in the context of dietary variation due to the divergent terrestrial ecology of D. polylepis compared to the arboreal niche occupied by all other mambas. Additionally, the degree of cross-reactivity conservation of mamba venoms was assessed by antivenomics against a panel of commercial antivenoms generated for the sub-Saharan Africa market. This study provides a genus-wide overview to infer which available antivenoms may be capable of neutralising human envenomings caused by mambas, irrespective of the species responsible. The information gathered in this study lays the foundations for rationalising the pharmacological profiles of mamba venoms at locus resolution. This understanding will contribute to the generation of a safer and more efficacious pan-Dendroaspis therapeutic antivenom against any mamba envenomation.

Identification of cDNAs encoding viper venom hyaluronidases: cross-generic sequence conservation of full-length and unusually short variant transcripts.

The immobilisation of prey by snakes is most efficiently achieved by the rapid dissemination of venom from its site of injection into the blood stream. Hyaluronidase is a common component of snake venoms and has been termed the "venom spreading factor". In the absence of nucleotide or protein sequence data to confirm the functional identity of this venom component, we interrogated a venom gland EST database for the saw-scaled viper, Echis ocellatus (Nigeria), using the gene ontology (GO) term "carbohydrate metabolism". A single hyalurononglucosaminadase-activity matching sequence (EOC00242) was found and used to design PCR primers to acquire the full-length cDNA sequence. Although very different from the bee venom and mammalian hyaluronidase sequences, the E. ocellatus sequence retained all the catalytic, positional and structural residues that characterise this class of carbohydrate metabolising hydrolases. An extraordinarily high level of sequence identity (>95%) was observed in analogous venom gland cDNA sequences isolated (by PCR) from another saw-scaled viper species, E. pyramidum leakeyi (Kenya), and from the sahara horned viper, Cerastes cerastes cerastes (Egypt) and the puff adder, Bitis arietans (Nigeria). Smaller amplicons, lacking hyaluronidase catalytic residues because of 768 bp or 855 bp central deletions, appear to encode either truncated peptides without hyaluronidase activity, or are non-translated transcripts because they lack consensus translation initiating motifs.

The Evolution of Fangs, Venom, and Mimicry Systems in Blenny Fishes.

Venom systems have evolved on multiple occasions across the animal kingdom, and they can act as key adaptations to protect animals from predators [1]. Consequently, venomous animals serve as models for a rich source of mimicry types, as non-venomous species benefit from reductions in predation risk by mimicking the coloration, body shape, and/or movement of toxic counterparts [2-5]. The frequent evolution of such deceitful imitations provides notable examples of phenotypic convergence and are often invoked as classic exemplars of evolution by natural selection. Here, we investigate the evolution of fangs, venom, and mimetic relationships in reef fishes from the tribe Nemophini (fangblennies). Comparative morphological analyses reveal that enlarged canine teeth (fangs) originated at the base of the Nemophini radiation and have enabled a micropredatory feeding strategy in non-venomous Plagiotremus spp. Subsequently, the evolution of deep anterior grooves and their coupling to venom secretory tissue provide Meiacanthus spp. with toxic venom that they effectively employ for defense. We find that fangblenny venom contains a number of toxic components that have been independently recruited into other animal venoms, some of which cause toxicity via interactions with opioid receptors, and result in a multifunctional biochemical phenotype that exerts potent hypotensive effects. The evolution of fangblenny venom has seemingly led to phenotypic convergence via the formation of a diverse array of mimetic relationships that provide protective (Batesian mimicry) and predatory (aggressive mimicry) benefits to other fishes [2, 6]. Our results further our understanding of how novel morphological and biochemical adaptations stimulate ecological interactions in the natural world.