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.

Repurposing the orphan drug nitisinone to control the transmission of African trypanosomiasis.

Tsetse transmit African trypanosomiasis, which is a disease fatal to both humans and animals. A vaccine to protect against this disease does not exist so transmission control relies on eliminating tsetse populations. Although neurotoxic insecticides are the gold standard for insect control, they negatively impact the environment and reduce populations of insect pollinator species. Here we present a promising, environment-friendly alternative to current insecticides that targets the insect tyrosine metabolism pathway. A bloodmeal contains high levels of tyrosine, which is toxic to haematophagous insects if it is not degraded and eliminated. RNA interference (RNAi) of either the first two enzymes in the tyrosine degradation pathway (tyrosine aminotransferase (TAT) and 4-hydroxyphenylpyruvate dioxygenase (HPPD)) was lethal to tsetse. Furthermore, nitisinone (NTBC), an FDA-approved tyrosine catabolism inhibitor, killed tsetse regardless if the drug was orally or topically applied. However, oral administration of NTBC to bumblebees did not affect their survival. Using a novel mathematical model, we show that NTBC could reduce the transmission of African trypanosomiasis in sub-Saharan Africa, thus accelerating current disease elimination programmes.

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.

The transcription factor Maf-S regulates metabolic resistance to insecticides in the malaria vector Anopheles gambiae.

BACKGROUND: Malaria control in Africa is dependent upon the use insecticides but intensive use of a limited number of chemicals has led to resistance in mosquito populations. Increased production of enzymes that detoxify insecticides is one of the most potent resistance mechanisms. Several metabolic enzymes have been implicated in insecticide resistance but the processes controlling their expression have remained largely elusive. RESULTS: Here, we show that the transcription factor Maf-S regulates expression of multiple detoxification genes, including the key insecticide metabolisers CYP6M2 and GSTD1 in the African malaria vector Anopheles gambiae. Attenuation of this transcription factor through RNAi induced knockdown reduced transcript levels of these effectors and significantly increased mortality after exposure to the pyrethroid insecticides and DDT (permethrin: 9.2% to 19.2% (p = 0.015), deltamethrin: 3.9% to 21.6% (p = 0.036) and DDT: 1% to 11.7% (p = <0.01), whilst dramatically decreasing mortality induced by the organophosphate malathion (79.6% to 8.0% (p = <0.01)). Additional genes regulated by Maf-S were also identified providing new insight into the role of this transcription factor in insects. CONCLUSION: Maf-S is a key regulator of detoxification genes in Anopheles mosquitoes. Disrupting this transcription factor has opposing effects on the mosquito's response to different insecticide classes providing a mechanistic explanation to the negative cross resistance that has been reported between pyrethroids and organophosphates.

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.

Intraspecific venom variation in the medically important puff adder (Bitis arietans): Comparative venom gland transcriptomics, in vitro venom activity and immunological recognition by antivenom.

BACKGROUND: Variation in snake venoms is well documented, both between and within species, with intraspecific venom variation often correlated with geographically distinct populations. The puff adder, Bitis arietans, is widely distributed across sub-Saharan Africa and into the Arabian Peninsula where it is considered a leading cause of the ~310,000 annual snakebites across the region, with its venom capable of causing substantial morbidity and mortality. Despite its medical importance and wide geographic distribution, there is little known about venom variation between different B. arietans populations and the potential implications of this variation on antivenom efficacy. METHODOLOGY: We applied a range of analyses, including venom gland transcriptomics, in vitro enzymatic assays and reverse phase chromatography to comparatively analyse B. arietans venoms originating from Nigeria, Tanzania, and South Africa. Immunological assays and in vitro enzymatic neutralisation assays were then applied to investigate the impact of venom variation on the potential efficacy of three antivenom products; SAIMR Polyvalent, EchiTAb-Plus and Fav-Afrique. FINDINGS: Through the first comparison of venom gland transcriptomes of B. arietans from three geographically distinct regions (Nigeria, Tanzania, and South Africa), we identified substantial variation in toxin expression. Findings of venom variation were further supported by chromatographic venom profiling, and the application of enzymatic assays to quantify the activity of three pathologically relevant toxin families. However, the use of western blotting, ELISA, and in vitro enzymatic inhibition assays revealed that variation within B. arietans venom does not appear to substantially impact upon the efficacy of three African polyvalent antivenoms. CONCLUSIONS: The large distribution and medical importance of B. arietans makes this species ideal for understanding venom variation and the impact this has on therapeutic efficacy. The findings in this study highlight the likelihood for considerable venom toxin variation across the range of B. arietans, but that this may not dramatically impact upon the utility of treatment available in the region.

Diverse RNA viruses of parasitic nematodes can elicit antibody responses in vertebrate hosts.

Parasitic nematodes have an intimate, chronic and lifelong exposure to vertebrate tissues. Here we mined 41 published parasitic nematode transcriptomes from vertebrate hosts and identified 91 RNA viruses across 13 virus orders from 24 families in ~70% (28 out of 41) of parasitic nematode species, which include only 5 previously reported viruses. We observe widespread distribution of virus-nematode associations across multiple continents, suggesting an ancestral acquisition event and host-virus co-evolution. Characterization of viruses of Brugia malayi (BMRV1) and Onchocerca volvulus (OVRV1) shows that these viruses are abundant in reproductive tissues of adult parasites. Importantly, the presence of BMRV1 RNA in B. malayi parasites mounts an RNA interference response against BMRV1 suggesting active viral replication. Finally, BMRV1 and OVRV1 were found to elicit antibody responses in serum samples from infected jirds and infected or exposed humans, indicating direct exposure to the immune system.

Different PfEMP1-expressing Plasmodium falciparum variants induce divergent endothelial transcriptional responses during co-culture.

The human malaria parasite Plasmodium falciparum is responsible for the majority of mortality and morbidity caused by malaria infection and differs from other human malaria species in the degree of accumulation of parasite-infected red blood cells in the microvasculature, known as cytoadherence or sequestration. In P. falciparum, cytoadherence is mediated by a protein called PfEMP1 which, due to its exposure to the host immune system, undergoes antigenic variation resulting in the expression of different PfEMP1 variants on the infected erythrocyte membrane. These PfEMP1s contain various combinations of adhesive domains, which allow for the differential engagement of a repertoire of endothelial receptors on the host microvasculature, with specific receptor usage associated with severe disease. We used a co-culture model of cytoadherence incubating human brain microvascular endothelial cells with erythrocytes infected with two parasite lines expressing different PfEMP1s that demonstrate different binding profiles to vascular endothelium. We determined the transcriptional profile of human brain microvascular endothelial cells (HBMEC) following different incubation periods with infected erythrocytes, identifying different transcriptional profiles of pathways previously found to be involved in the pathology of severe malaria, such as inflammation, apoptosis and barrier integrity, induced by the two PfEMP1 variants.

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.

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.

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.