A Recurrent Motif: Diversity and Evolution of ShKT Domain Containing Proteins in the Vampire Snail Cumia reticulata

Functional Characterization of ShK Domain-Containing Protein in the Plant-Parasitic Nematode Bursaphelenchus xylophilus

ShK domain-containing proteins are peptides found in different parasitic and venomous organisms. From a previous transcriptomic dataset from Bursaphelenchus xylophilus, a plant-parasitic nematode that infects forest tree species, we identified 96 transcripts potentially as ShK domain-containing proteins with unknown function in the nematode genome. This study aimed to characterize and explore the functional role of genes encoding ShK domain-containing proteins in B. xylophilus biology. We selected and functionally analyzed nine candidate genes that are putatively specific to B. xylophilus. In situ hybridization revealed expression of one B. xylophilus ShK in the pharyngeal gland cells, suggesting their delivery into host cells. Most of the transcripts are highly expressed during infection and showed a significant upregulation in response to peroxide products compared to the nematode catalase enzymes. We reported, for the first time, the potential involvement of ShK domain genes in oxidative stress, suggesting that these proteins may have an important role in protecting or modulating the reactive oxygen species (ROS) activity of the host plant during parasitism.

Collaborative Expression: Transcriptomics of Conus virgo Suggests Contribution of Multiple Secretory Glands to Venom Production

Venomous marine gastropods of the family Conidae are among the most diversified predators in marine realm-in large due to their complex venoms. Besides being a valuable source of bioactive neuropeptides conotoxins, cone-snails venoms are an excellent model for molecular evolution studies, addressing origin of key innovations. However, these studies are handicapped by scarce current knowledge on the tissues involved in venom production, as it is generally assumed the sole prerogative of the venom gland (VG). The role of other secretory glands that are present in all Conus species (salivary gland, SG) or only in some species (accessory salivary gland, ASG) remains poorly understood. Here, for the first time, we carry out a detailed analysis of the VG, SG, and ASG transcriptomes in the vermivorous Conus virgo. We detect multiple transcripts clusters in both the SG and ASG, whose annotations imply venom-related functions. Despite the subsets of transcripts highly-expressed in the VG, SG, and ASG being very distinct, SG expresses an L-, and ASG-Cerm08-, and MEFRR- superfamily conotoxins, all previously considered specific for VG. We corroborate our results with the analysis of published SG and VG transcriptomes from unrelated fish-hunting C. geographus, and C. striatus, possibly fish-hunting C. rolani, and worm-hunting Conus quercinus. In spite of low expression levels of conotoxins, some other specific clusters of putative venom-related peptides are present and may be highly expressed in the SG of these species. Further functional studies are necessary to determine the role that these peptides play in envenomation. In the meantime, our results show importance of routine multi-tissue sampling both for accurate interpretation of tissue-specific venom composition in cone-snails, and for better understanding origin and evolution of venom peptides genes.

Whole Genome Duplication and Gene Evolution in the Hyperdiverse Venomous Gastropods

The diversity of venomous organisms and the toxins they produce have been increasingly investigated, but taxonomic bias remains important. Neogastropods, a group of marine predators representing almost 22% of the known gastropod diversity, evolved a wide range of feeding strategies, including the production of toxins to subdue their preys. However, whether the diversity of these compounds is at the origin of the hyperdiversification of the group and how genome evolution may correlate with both the compounds and species diversities remain understudied. Among the available gastropods genomes, only eight, with uneven quality assemblies, belong to neogastropods. Here, we generated chromosome-level assemblies of two species belonging to the Tonnoidea and Muricoidea superfamilies (Monoplex corrugatus and Stramonita haemastoma). The two obtained high-quality genomes had 3 and 2.2 Gb, respectively, and 92-89% of the total assembly conformed 35 pseudochromosomes in each species. Through the analysis of syntenic blocks, Hox gene cluster duplication, and synonymous substitutions distribution pattern, we inferred the occurrence of a whole genome duplication event in both genomes. As these species are known to release venom, toxins were annotated in both genomes, but few of them were found in homologous chromosomes. A comparison of the expression of ohnolog genes (using transcriptomes from osphradium and salivary glands in S. haemastoma), where both copies were differentially expressed, showed that most of them had similar expression profiles. The high quality of these genomes makes them valuable reference in their respective taxa, facilitating the identification of genome-level processes at the origin of their evolutionary success.

An informatic workflow for the enhanced annotation of excretory/secretory proteins of Haemonchus contortus

Major advances in genomic and associated technologies have demanded reliable bioinformatic tools and workflows for the annotation of genes and their products via comparative analyses using well-curated reference data sets, accessible in public repositories. However, the accurate in silico annotation of molecules (proteins) encoded in organisms (e.g., multicellular parasites) which are evolutionarily distant from those for which these extensive reference data sets are available, including invertebrate model organisms (e.g., Caenorhabditis elegans - free-living nematode, and Drosophila melanogaster - the vinegar fly) and vertebrate species (e.g., Homo sapiens and Mus musculus), remains a major challenge. Here, we constructed an informatic workflow for the enhanced annotation of biologically-important, excretory/secretory (ES) proteins ("secretome") encoded in the genome of a parasitic roundworm, called Haemonchus contortus (commonly known as the barber's pole worm). We critically evaluated the performance of five distinct methods, refined some of them, and then combined the use of all five methods to comprehensively annotate ES proteins, according to gene ontology, biological pathways and/or metabolic (enzymatic) processes. Then, using optimised parameter settings, we applied this workflow to comprehensively annotate 2591 of all 3353 proteins (77.3%) in the secretome of H. contortus. This result is a substantial improvement (10-25%) over previous annotations using individual, "off-the-shelf" algorithms and default settings, indicating the ready applicability of the present, refined workflow to gene/protein sequence data sets from a wide range of organisms in the Tree-of-Life.

Structural and evolutionary insights into astacin metallopeptidases

The astacins are a family of metallopeptidases (MPs) that has been extensively described from animals. They are multidomain extracellular proteins, which have a conserved core architecture encompassing a signal peptide for secretion, a prodomain or prosegment and a zinc-dependent catalytic domain (CD). This constellation is found in the archetypal name-giving digestive enzyme astacin from the European crayfish Astacus astacus. Astacin catalytic domains span ∼200 residues and consist of two subdomains that flank an extended active-site cleft. They share several structural elements including a long zinc-binding consensus sequence (HEXXHXXGXXH) immediately followed by an EXXRXDRD motif, which features a family-specific glutamate. In addition, a downstream SIMHY-motif encompasses a "Met-turn" methionine and a zinc-binding tyrosine. The overall architecture and some structural features of astacin catalytic domains match those of other more distantly related MPs, which together constitute the metzincin clan of metallopeptidases. We further analysed the structures of PRO-, MAM, TRAF, CUB and EGF-like domains, and described their essential molecular determinants. In addition, we investigated the distribution of astacins across kingdoms and their phylogenetic origin. Through extensive sequence searches we found astacin CDs in > 25,000 sequences down the tree of life from humans beyond Metazoa, including Choanoflagellata, Filasterea and Ichtyosporea. We also found < 400 sequences scattered across non-holozoan eukaryotes including some fungi and one virus, as well as in selected taxa of archaea and bacteria that are pathogens or colonizers of animal hosts, but not in plants. Overall, we propose that astacins originate in the root of Holozoa consistent with Darwinian descent and that the latter genes might be the result of horizontal gene transfer from holozoan donors.

ShK-Domain-Containing Protein from a Parasitic Nematode Modulates Drosophila melanogaster Immunity

A key component to understanding host–parasite interactions is the molecular crosstalk between host and parasite. Excreted/secreted products (ESPs) released by parasitic nematodes play an important role in parasitism. They can directly damage host tissue and modulate host defense. Steinernema carpocapsae, a well-studied parasite of insects releases approximately 500 venom proteins as part of the infection process. Though the identity of these proteins is known, few have been studied in detail. One protein family present in the ESPs released by these nematodes is the ShK family. We studied the most abundant ShK-domain-containing protein in S. carpocapsae ESPs, Sc-ShK-1, to investigate its effects in a fruit fly model. We found that Sc-ShK-1 is toxic under high stress conditions and negatively affects the health of fruit flies. We have shown that Sc-ShK-1 contributes to host immunomodulation in bacterial co-infections resulting in increased mortality and microbial growth. This study provides an insight on ShK-domain-containing proteins from nematodes and suggests these proteins may play an important role in host–parasite interactions.

Vexitoxins: conotoxin-like venom peptides from predatory gastropods of the genus Vexillum

Venoms of predatory marine cone snails are intensely studied because of the biomedical applications of the neuropeptides that they contain, termed conotoxins. Meanwhile some gastropod lineages have independently acquired secretory glands strikingly similar to the venom gland of cone snails, suggesting that they possess similar venoms. Here we focus on the most diversified of these clades, the genus Vexillum. Based on the analysis of a multi-species proteo-transcriptomic dataset, we show that Vexillum species indeed produce complex venoms dominated by highly diversified short cysteine-rich peptides, vexitoxins. Vexitoxins possess the same precursor organization, display overlapping cysteine frameworks and share several common post-translational modifications with conotoxins. Some vexitoxins show sequence similarity to conotoxins and adopt similar domain conformations, including a pharmacologically relevant inhibitory cysteine knot motif. The Vexillum envenomation gland (gL) is a notably more recent evolutionary novelty than the conoidean venom gland. Thus, we hypothesize lower divergence between vexitoxin genes, and their ancestral 'somatic' counterparts compared to that in conotoxins, and we find support for this hypothesis in the evolution of the vexitoxin cluster V027. We use this example to discuss how future studies on vexitoxins can inform the origin of conotoxins, and how they may help to address outstanding questions in venom evolution.

The cnidarian parasite Ceratonova shasta utilizes inherited and recruited venom-like compounds during infection

Background

Cnidarians are the most ancient venomous organisms. They store a cocktail of venom proteins inside unique stinging organelles called nematocysts. When a cnidarian encounters chemical and physical cues from a potential threat or prey animal, the nematocyst is triggered and fires a harpoon-like tubule to penetrate and inject venom into the prey. Nematocysts are present in all Cnidaria, including the morphologically simple Myxozoa, which are a speciose group of microscopic, spore-forming, obligate parasites of fish and invertebrates. Rather than predation or defense, myxozoans use nematocysts for adhesion to hosts, but the involvement of venom in this process is poorly understood. Recent work shows some myxozoans have a reduced repertoire of venom-like compounds (VLCs) relative to free-living cnidarians, however the function of these proteins is not known.

Methods

We searched for VLCs in the nematocyst proteome and a time-series infection transcriptome of Ceratonova shasta, a myxozoan parasite of salmonid fish. We used four parallel approaches to detect VLCs: BLAST and HMMER searches to preexisting cnidarian venom datasets, the machine learning tool ToxClassifier, and structural modeling of nematocyst proteomes. Sequences that scored positive by at least three methods were considered VLCs. We then mapped their time-series expressions in the fish host and analyzed their phylogenetic relatedness to sequences from other venomous animals.

Results

We identified eight VLCs, all of which have closely related sequences in other myxozoan datasets, suggesting a conserved venom profile across Myxozoa, and an overall reduction in venom diversity relative to free-living cnidarians. Expression of the VLCs over the 3-week fish infection varied considerably: three sequences were most expressed at one day post-exposure in the fish’s gills; whereas expression of the other five VLCs peaked at 21 days post-exposure in the intestines, coinciding with the formation of mature parasite spores with nematocysts. Expression of VLC genes early in infection, prior to the development of nematocysts, suggests venoms in C. shasta have been repurposed to facilitate parasite invasion and proliferation within the host. Molecular phylogenetics suggested some VLCs were inherited from a cnidarian ancestor, whereas others were more closely related to sequences from venomous non-Cnidarian organisms and thus may have gained qualities of venom components via convergent evolution. The presence of VLCs and their differential expression during parasite infection enrich the concept of what functions a “venom” can have and represent targets for designing therapeutics against myxozoan infections.

The Origin and Evolution of Antistasin-like Proteins in Leeches (Hirudinida, Clitellata)

Bloodfeeding is employed by many parasitic animals and requires specific innovations for efficient feeding. Some of these innovations are molecular features that are related to the inhibition of hemostasis. For example, bloodfeeding insects, bats, and leeches release proteins with anticoagulatory activity through their salivary secretions. The antistasin-like protein family, composed of serine protease inhibitors with one or more antistasin-like domains, is tightly linked to inhibition of hemostasis in leeches. However, this protein family has been recorded also in non-bloodfeeding invertebrates, such as cnidarians, mollusks, polychaetes, and oligochaetes. The present study aims to 1) root the antistasin-like gene tree and delimit the major orthologous groups, 2) identify potential independent origins of salivary proteins secreted by leeches, and 3) identify major changes in domain and/or motif structure within each orthologous group. Five clades containing leech antistasin-like proteins are distinguishable through rigorous phylogenetic analyses based on nine new transcriptomes and a diverse set of comparative data: the trypsin + leukocyte elastase inhibitors clade, the antistasin clade, the therostasin clade, and two additional, unnamed clades. The antistasin-like gene tree supports multiple origins of leech antistasin-like proteins due to the presence of both leech and non-leech sequences in one of the unnamed clades, but a single origin of factor Xa and trypsin + leukocyte elastase inhibitors. This is further supported by three sequence motifs that are exclusive to antistasins, the trypsin + leukocyte elastase inhibitor clade, and the therostasin clade, respectively. We discuss the implications of our findings for the evolution of this diverse family of leech anticoagulants.