Transcriptomic profiling and the first spatial expression analysis of candidate genes in the salivary gland of the East Asian medicinal leech, Hirudo nipponia

Hirudo nipponia, a blood-sucking leech native to East Asia, possesses a rich repertoire of active ingredients in its saliva, showcasing significant medical potential due to its anticoagulant, anti-inflammatory, and antibacterial effects against human diseases. Despite previous studies on the transcriptomic and proteomic characteristics of leech saliva, which have identified medicinal compounds, our knowledge of tissue-specific transcriptomes and their spatial expression patterns remains incomplete. In this study, we conducted an extensive transcriptomic profiling of the salivary gland tissue in H. nipponia based on de novo assemblies of tissue-specific transcriptomes from the salivary gland, teeth, and general head region. Through gene ontology (GO) analysis and hierarchical clustering, we discovered a novel set of anti-coagulant factors-i.e., Hni-Antistasin, Hni-Ghilanten, Hni-Bdellin, Hni-Hirudin-as well as a previously unrecognized immune-related gene, Hni-GLIPR1 and uncharacterized salivary gland specific transcripts. By employing in situ hybridization, we provided the first visualization of gene expression sites within the salivary gland of H. nipponia. Our findings expand on our understanding of transcripts specifically expressed in the salivary gland of blood-sucking leeches, offering valuable resources for the exploration of previously unidentified substances with medicinal applications.

Cloning, characterization, and heterologous expression of a candidate Hirudin gene from the salivary gland transcriptome of Hirudo nipponia

Hirudin is a pharmacologically active substance in leeches with potent blood anticoagulation properties. Although recombinant hirudin production isolated from Hirudo medicinalis Linnaeus and Hirudinaria manillensis Lesson is known, to our knowledge, this study is the first to report recombinant hirudin expression and production from Hirudo nipponia Whitman. Thus, the present study aimed to clone and characterize the full-length cDNA of a candidate hirudin gene (c16237_g1), which is localized on the salivary gland transcriptome of H. nipponia, and further evaluate its recombinant production using a eukaryotic expression system. The 489-bp cDNA possessed several properties of the hirudin “core” motifs associated with binding to the thrombin catalytic pocket. A fusion expression vector (pPIC9K-hirudin) was constructed and successfully transformed into Pichia pastoris strain GS115 via electroporation. Sodium dodecyl sulphate–polyacrylamide gel electrophoresis and western blot analysis confirmed hirudin expression. The recombinant protein was expressed with a yield of 6.68 mg/L culture. Mass spectrometry analysis further confirmed target protein expression. The concentration and antithrombin activity of purified hirudin were 1.67 mg/mL and 14,000 ATU/mL, respectively. These findings provide a basis for further elucidating the molecular anticoagulation mechanism of hirudin, and address China’s growing market demand for engineered H. nipponia-derived hirudin and hirudin-based drugs.

Molecular cloning and functional analysis of HnSaratin from Hirudo nipponia

In order to finish a bloodmeal successfully, hematophagous organisms often stored a variety of anticoagulant proteins in their salivary glands, such as proteins that inhibit platelet aggregation. When they ingest a bloodmeal, these proteins are injected into the host to prevent the blood from clotting. As one of the origins of leeches used in traditional Chinese medicine, H. nipponia was proved to be clinically effective in treatment of cardiovascular and cerebrovascular diseases. This study cloned the sequence of HnSaratin cDNA derived from salivary glands of H. nipponia. The sequence contains an open reading frame of 387 bp, encoding a protein of 128 amino acids containing a signal peptide of 21 amino acids. After removal of the signal peptide, the molecular mass of mature HnSaratin was 12.37 kDa, with a theoretical isoelectric point (pI) of 3.89. The N-terminal of mature HnSaratin was folded into a globular structure, in which 3 disulfide bonds, a ββαβββ topology and 2 Glu residues that binds collagenous Lys2 were located, and the C-terminal formed a flexible region. The fusion HnSaratin protein was obtained by a prokaryotic expression system. The protein showed anti-platelet aggregation activity, and was observed to prevent blood clotting in rats. The significant high expression of HnSaratin mRNA in salivary glands was induced by bloodmeal ingestion of H. nipponia. Briefly, our work provides theoretical basis for further development and utilization of H. nipponia.

Venom peptides - A comprehensive translational perspective in pain management

Venom peptides have been evolving complex therapeutic interventions that potently and selectively modulate a range of targets such as ion channels, receptors, and signaling pathways of physiological processes making it potential therapeutic. Several venom peptides were deduced in vivo for clinical development targeting pain management, diabetes, cardiovascular diseases, antimicrobial activity. Several contributions have been detailed for a clear perspective for a better understanding of venomous animals, their venom, and their pharmacological effects. Here we unravel and summarize the recent advances in wide venom peptides across varieties of species for their therapeutics prospects.

Functional morphology of the suckers and teeth of the medicinal leech Hirudo verbana Carena, 1820 (Annelida; Clitellata; Hirudinida): A scanning electron microscope study

In this study, the triple jaws and suckers of the leeches belonging to the Hirudo verbana Carena, 1820 (Annelida; Clitellata; Hirudinida) were examined using the stereomicroscopy and scanning electron microscopy. In H. verbana, suckers are seen on the first annulus and last annulus of the body. The mouth opens in the center of the front suction cup, and behind this opening is a movable triple jaw apparatus with many teeth. The posterior sucker disc consists of the last seven body segments and lacks an opening. The shape of the jaw is trignatous. The pharynx is equally located around of the three muscular jaws. The jaws are muscular covered with cuticle and carry a row of teeth arranged at the tip. In this study, it was determined that secretory canal holes were identified between the teeth. The results show that the size of teeth determines long-term bleeding so revealing the structure and working mechanism of the teeth has importance for medicinal leeches. At the same time, the difference of teeth and jaw structures of leeches may be a criterion in the classification of medicinal leeches.

Macrobdella decora: Old World Leech Gut Microbial Community Structure Conserved in a New World Leech

Leeches are found in terrestrial, aquatic, and marine habitats on all continents. Sanguivorous leeches have been used in medicine for millennia. Modern scientific uses include studies of neurons, anticoagulants, and gut microbial symbioses. Hirudo verbana, the European medicinal leech, maintains a gut community dominated by two bacterial symbionts, Aeromonas veronii and Mucinivorans hirudinis, which sometimes account for as much as 97% of the total crop microbiota. The highly simplified gut anatomy and microbiome of H. verbana make it an excellent model organism for studying gut microbial dynamics. The North American medicinal leech, Macrobdella decora, is a hirudinid leech native to Canada and the northern United States. In this study, we show that M. decora symbiont communities are very similar to those in H. verbana. We performed an extensive study using field-caught M. decora and purchased H. verbana from two suppliers. Deep sequencing of the V4 region of the 16S rRNA gene allowed us to determine that the core microbiome of M. decora consists of Bacteroides, Aeromonas, Proteocatella, and Butyricicoccus. The analysis revealed that the compositions of the gut microbiomes of the two leech species were significantly different at all taxonomic levels. The R ² value was highest at the genus and amplicon sequence variant (ASV) levels and much lower at the phylum, class, and order levels. The gut and bladder microbial communities were distinct. We propose that M. decora is an alternative to H. verbana for studies of wild-caught animals and provide evidence for the conservation of digestive-tract and bladder symbionts in annelid models.IMPORTANCE Building evidence implicates the gut microbiome in critical animal functions such as regulating digestion, nutrition, immune regulation, and development. Simplified, phylogenetically diverse models for hypothesis testing are necessary because of the difficulty of assigning causative relationships in complex gut microbiomes. Previous research used Hirudo verbana as a tractable animal model of digestive-tract symbioses. Our data show that Macrobdella decora may work just as well without the drawback of being an endangered organism and with the added advantage of easy access to field-caught specimens. The similarity of the microbial community structures of species from two different continents reveals the highly conserved nature of the microbial symbionts in sanguivorous leeches.

Diet Diversity in Carnivorous Terebrid Snails Is Tied to the Presence and Absence of a Venom Gland

Predator-prey interactions are thought to play a driving role in animal evolution, especially for groups that have developed venom as their predatory strategy. However, how the diet of venomous animals influences the composition of venom arsenals remains uncertain. Two prevailing hypotheses to explain the relationship between diet and venom composition focus on prey preference and the types of compounds in venom, and a positive correlation between dietary breadth and the number of compounds in venom. Here, we examined venom complexity, phylogenetic relationship, collection depth, and biogeography of the Terebridae (auger snails) to determine if repeated innovations in terebrid foregut anatomy and venom composition correspond to diet variation. We performed the first molecular study of the diet of terebrid marine snails by metabarcoding the gut content of 71 terebrid specimens from 17 species. Our results suggest that the presence or absence of a venom gland is strongly correlated with dietary breadth. Specifically, terebrid species without a venom gland displayed greater diversity in their diet. Additionally, we propose a revision of the definition of venom complexity in conoidean snails to more accurately capture the breadth of ecological influences. These findings suggest that prey diet is an important factor in terebrid venom evolution and diversification and further investigations of other understudied organisms, like terebrids, are needed to develop robust hypotheses in this area.

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.

Genome-enabled insights into the biology of thrips as crop pests

BACKGROUND

The western flower thrips, Frankliniella occidentalis (Pergande), is a globally invasive pest and plant virus vector on a wide array of food, fiber, and ornamental crops. The underlying genetic mechanisms of the processes governing thrips pest and vector biology, feeding behaviors, ecology, and insecticide resistance are largely unknown. To address this gap, we present the F. occidentalis draft genome assembly and official gene set.

RESULTS

We report on the first genome sequence for any member of the insect order Thysanoptera. Benchmarking Universal Single-Copy Ortholog (BUSCO) assessments of the genome assembly (size = 415.8 Mb, scaffold N50 = 948.9 kb) revealed a relatively complete and well-annotated assembly in comparison to other insect genomes. The genome is unusually GC-rich (50%) compared to other insect genomes to date. The official gene set (OGS v1.0) contains 16,859 genes, of which ~ 10% were manually verified and corrected by our consortium. We focused on manual annotation, phylogenetic, and expression evidence analyses for gene sets centered on primary themes in the life histories and activities of plant-colonizing insects. Highlights include the following: (1) divergent clades and large expansions in genes associated with environmental sensing (chemosensory receptors) and detoxification (CYP4, CYP6, and CCE enzymes) of substances encountered in agricultural environments; (2) a comprehensive set of salivary gland genes supported by enriched expression; (3) apparent absence of members of the IMD innate immune defense pathway; and (4) developmental- and sex-specific expression analyses of genes associated with progression from larvae to adulthood through neometaboly, a distinct form of maturation differing from either incomplete or complete metamorphosis in the Insecta.

CONCLUSIONS

Analysis of the F. occidentalis genome offers insights into the polyphagous behavior of this insect pest that finds, colonizes, and survives on a widely diverse array of plants. The genomic resources presented here enable a more complete analysis of insect evolution and biology, providing a missing taxon for contemporary insect genomics-based analyses. Our study also offers a genomic benchmark for molecular and evolutionary investigations of other Thysanoptera species.

Medical leech therapy in Ayurveda and biomedicine - A review

Medicinal leech therapy or Hirudino therapy have roots back in ancient civilization. It was a prevalent form of therapy in various ailments. The novel Ayurveda text Sushruta Samhita devoted a complete chapter on hirudino therapy. In the early 20th century this therapy had a major setback due to origin and evolution of antibiotics. There was a discontinuity in the flow of knowledge about this therapy. Then, resumed and revived after few recent decades, due to its contribution in reconstructive surgeries. During this period, the research work on various aspects have been conducted. The present paper summarizes the various aspects of medicinal leech therapy both from Ayurveda text and the present knowledge and to enable the fraternity to use the both source for benefits of humankind.

Leeches in the extreme: Morphological, physiological, and behavioral adaptations to inhospitable habitats

With more than 700 described species, leeches include morphological, physiological, and behavioral diversity and occur in terrestrial and aquatic habitats, including freshwater, estuarine, and marine ecosystems. Leeches inhabit a number of extreme environments, including extremes in temperature, moisture, salinity, pressure, light, and pollution. In some cases, leeches in extreme environments have specialized morphological, physiological, or behavioral adaptations to survive these conditions, yet unique adaptations are not apparent in some species. Leeches that inhabit inhospitable habitats occur in more than one branch or family of leech phylogeny suggesting that there have been independent invasions of environments with extreme conditions. Herein, we review examples of leeches that live in extreme conditions and the exceptional biology that has contributed to leeches being the most extreme annelids.

Draft genome of the European medicinal leech Hirudo medicinalis (Annelida, Clitellata, Hirudiniformes) with emphasis on anticoagulants

The European medicinal leech has been used for medicinal purposes for millennia, and continues to be used today in modern hospital settings. Its utility is granted by the extremely potent anticoagulation factors that the leech secretes into the incision wound during feeding and, although a handful of studies have targeted certain anticoagulants, the full range of anticoagulation factors expressed by this species remains unknown. Here, we present the first draft genome of the European medicinal leech, Hirudo medicinalis, and estimate that we have sequenced between 79–94% of the full genome. Leveraging these data, we searched for anticoagulation factors across the genome of H. medicinalis. Following orthology determination through a series of BLAST searches, as well as phylogenetic analyses, we estimate that fully 15 different known anticoagulation factors are utilized by the species, and that 17 other proteins that have been linked to antihemostasis are also present in the genome. We underscore the utility of the draft genome for comparative studies of leeches and discuss our results in an evolutionary context.

Draft genome sequences of Hirudo medicinalis and salivary transcriptome of three closely related medicinal leeches

BACKGROUND

Salivary cell secretion (SCS) plays a critical role in blood feeding by medicinal leeches, making them of use for certain medical purposes even today.

RESULTS

We annotated the Hirudo medicinalis genome and performed RNA-seq on salivary cells isolated from three closely related leech species, H. medicinalis, Hirudo orientalis, and Hirudo verbana. Differential expression analysis verified by proteomics identified salivary cell-specific gene expression, many of which encode previously unknown salivary components. However, the genes encoding known anticoagulants have been found to be expressed not only in salivary cells. The function-related analysis of the unique salivary cell genes enabled an update of the concept of interactions between salivary proteins and components of haemostasis.

CONCLUSIONS

Here we report a genome draft of Hirudo medicinalis and describe identification of novel salivary proteins and new homologs of genes encoding known anticoagulants in transcriptomes of three medicinal leech species. Our data provide new insights in genetics of blood-feeding lifestyle in leeches.

European Medicinal Leeches-New Roles in Modern Medicine

Before the advent of modern medicine, natural resources were widely used by indigenous populations for the prevention and treatment of diseases. The associated knowledge, collectively described as folk medicine or traditional medicine, was largely based on trial-and-error testing of plant extracts (herbal remedies) and the use of invertebrates, particularly medicinal maggots of the blowfly Lucilia sericata and blood-sucking leeches. The widespread use of traditional medicine in the West declined as scientific advances allowed reproducible testing under controlled conditions and gave rise to the modern fields of biomedical research and pharmacology. However, many drugs are still derived from natural resources, and interest in traditional medicine has been renewed by the ability of researchers to investigate the medical potential of diverse species by high-throughput screening. Likewise, researchers are starting to look again at the benefits of maggot and leech therapy, based on the hypothesis that the use of such animals in traditional medicine is likely to reflect the presence of specific bioactive molecules that can be developed as drug leads. In this review, we consider the modern medical benefits of European medicinal leeches based on the systematic screening of their salivary proteins.

Identification and characterization of hirudin-HN, a new thrombin inhibitor, from the salivary glands of Hirudo nipponia

Transcriptome sequencing data (6.5 Gb) of the salivary glands of the haematophagous leech Hirudo nipponia was obtained by using the BGIseq-500 platform. After identification and analysis, one transcript (Unigene5370) was annotated to hirudin HV3 from Hirudo medicinalis with an e-value of 1e-29 and was named hirudin-HN. This transcript was a new thrombin inhibitor gene belonging to the proteinase inhibitor I14 (hirudin) family. Hirudin-HN, with a 270-bp cDNA, encodes an 89-aa protein containing a 20-aa signal peptide. The mature hirudin-HN protein contains the typical structural characteristics of hirudin, e.g., three conserved disulfide bonds and the PKP and DFxxIP motifs. Proteins (Hir and M-Hir) were obtained via prokaryotic expression, and the mature hirudin-HN protein was shown to have anticoagulant activity and thrombin affinity by using the chromogenic substrate S2238 and surface plasmon resonance (SPR) interaction analysis, respectively. The N-terminal structure of the mature hirudin-HN protein was shown to be important for anticoagulant activity by comparing the activity and thrombin affinity of Hir and M-Hir. The abundances of Hirudin-HN mRNA and protein were higher in the salivary glands of starving animals than in those of feeding or fed leeches. These results provided a foundation for further study on the structure-function relationship of hirudin-HN with thrombin.

Spatiotemporal Expression of Anticoagulation Factor Antistasin in Freshwater Leeches

Antistasin, which was originally discovered in the salivary glands of the Mexican leech Haementeria officinalis, was newly isolated from Helobdella austinensis. To confirm the temporal expression of antistasin during embryogenesis, we carried out semi-quantitative RT-PCR. Hau-antistasin1 was uniquely expressed at stage 4 of the cleavage and was strongly expressed in the late stages of organogenesis, as were other antistasin members. In order to confirm the spatial expression of antistasin, we performed fluorescence in situ hybridization in the late stages of organogenesis. The expression of each antistasin in the proboscis showed a similar pattern and varied in expression in the body. In addition, the spatial expression of antistasin orthologs in different leeches showed the possibility of different function across leech species. Hau-antistasin1 was expressed in the same region as hedgehog, which is a known mediator of signal transduction pathway. Hau-antistasin1 is probably a downstream target of Hedgehog signaling, involved in segment polarity signal pathway.

Phylogenetic Position and Description of a New Species of Medicinal Leech from the Eastern United States

A new species of medicinal leech, Macrobdella mimicus n. sp., is described from specimens collected in Maryland; this is the first description of a North American macrobdellid since 1975. Superficially, the new species resembles the well-known Macrobdella decora, as both species possess 4 accessory pores arranged symmetrically on the ventral surface, yet the new species is distinguished from M. decora in possessing 4-4½ annuli (rather than 3½) between the gonopores and 4 annuli (rather than 5 annuli) between the female gonopore and the first pair of accessory pores. Phylogenetic analyses, based on 2 mitochondrial and 2 nuclear loci for a set of closely related taxa, confirms the placement of the new species within the family Macrobdellidae and places it as the sister taxon to M. decora and M. diplotertia.

The salivary transcriptome of Limnobdella mexicana (Annelida: Clitellata: Praobdellidae) and orthology determination of major leech anticoagulants

Bloodfeeding requires several adaptations that allow the parasite to feed efficiently. Leeches and other hematophagous animals have developed different mechanisms to inhibit hemostasis, one of the main barriers imposed by their hosts. Limnobdella mexicana is a member of the leech family Praobdellidae, a family of host generalists known for their preference to attach on mucosal membranes of mammals, such as those in nasopharyngeal cavities, bladders and ocular orbits. Previous studies have hypothesized a positive relationship between diversity of anticoagulants and diversity of hosts in bloodfeeding leeches. However, orthology determination of putative anticoagulants and the lack of standardization of sequencing effort and method hinder comparisons between publicly available transcriptomes generated in different laboratories. In the present study, we examine the first transcriptome of a praobdellid leech and identify 15 putative anticoagulants using a phylogeny-based inference approach, amino-acid conservation, Pfam domains and BLAST searches. Our phylogenetic analyses suggest that the ancestral leech was able to inhibit factor Xa and that some hirudins that have been reported in previous studies on leech anticoagulants may not be orthologous with the archetypal hirudin.

Transcriptomics and differential gene expression in Whitmania pigra (Annelida: Clitellata: Hirudinida: Hirudinidae): Contrasting feeding and fasting modes

The medicinal utility of leeches has been demonstrated through decades of use in modern hospital settings, mainly as relievers of venous congestion following flap or digit replantation surgery. In the present study, we sequence and annotate (through BLAST- and Gene Ontology-based approaches) the salivary transcriptome of the nonblood feeding hirudinid Whitmania pigra and assess the differential gene expression of anticoagulation factors (through both quantitative real-time PCR [qRT-PCR] and in silico-based methods) during feeding and fasting conditions. This was done in order to evince the diversity of putative anticoagulation factors, as well as estimate the levels of upregulation of genes immediately after feeding. In total, we found sequences with demonstrated orthology (via both phylogenetic analyses and BLAST-based approaches) to seven different proteins that have previously been linked to anticoagulatory capabilities-eglin C, bdellin, granulin, guamerin, hyaluronidase, destabilase I, and lipocalin. All of these were recovered from leeches both in the fasting and in the feeding conditions, but all show signs of upregulation in the feeding leeches. Interestingly, our RNA-seq effort, coupled with a hypergeometric test, indicated that the differentially expressed genes were disproportionately involved in three main immunological pathways (endocytosis, peroxisome regulation, and lysosome regulation). The results and implications of the finding of anticoagulants in this nonblood feeding leech and the putative upregulation of anticoagulation factors after feeding are briefly discussed in an evolutionary context.

Transcriptomic analysis of the salivary gland of medicinal leech Hirudo nipponia

Hirudo nipponia (known as Shui Zhi in Chinese) is a well-known Chinese medicine with numerous active ingredients in its body, especially in its saliva. This native Chinese blood-sucking leech has been used for therapeutic purposes since before 100 AD. Modern Chinese physicians use it for a wide range of diseases. Genomic data and molecular information about the pharmacologically active substances produced by this medicinal leech are presently unavailable despite this organism’s medicinal importance. In this study, we performed transcriptome profiling of the salivary glands of medicinal leech H. nipponia using the Illumina platform. In total, 84,657,362 clean reads were assembled into 50,535 unigenes. The obtained unigenes were compared to public databases. Furthermore, a unigene sequence similarity search and comparisons with the whole transcriptome of medical leech were performed to identify potential proteins. Finally, more than 21 genes were predicted to be involved in anticoagulatory, antithrombotic, antibacterial, anti-inflammatory and antitumor processes, which might play important roles in the treatment of various diseases. This study is the first analysis of a sialotranscriptome in H. nipponia. The transcriptome profile will shed light on its genetic background and provide a useful tool to deepen our understanding of the medical value of H. nipponia.

Analysis of codon usage patterns in Hirudinaria manillensis reveals a preference for GC-ending codons caused by dominant selection constraints

Background

Hirudinaria manillensis is an ephemeral, blood-sucking ectoparasite, possessing anticoagulant capacities with potential medical applications. Analysis of codon usage patterns would contribute to our understanding of the evolutionary mechanisms and genetic architecture of H. manillensis, which in turn would provide insight into the characteristics of other leeches. We analysed codon usage and related indices using 18,000 coding sequences (CDSs) retrieved from H. manillensis RNA-Seq data.

Results

We identified four highly preferred codons in H. manillensis that have G/C-endings. Points generated in an effective number of codons (ENC) plot distributed below the standard curve and the slope of a neutrality plot was less than 1. Highly expressed CDSs had lower ENC content and higher GC content than weakly expressed CDSs. Principal component analysis conducted on relative synonymous codon usage (RSCU) values divided CDSs according to GC content and divided codons according to ending bases. Moreover, by determining codon usage, we found that the majority of blood-diet related genes have undergone less adaptive evolution in H. manillensis, except for those with homologous sequences in the host species.

Conclusions

Codon usage in H. manillensis had an overall preference toward C-endings and indicated that codon usage patterns are mediated by differential expression, GC content, and biological function. Although mutation pressure effects were also notable, the majority of genetic evolution in H. manillensis was driven by natural selection.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4937-x) contains supplementary material, which is available to authorized users.

Marine Leech Anticoagulant Diversity and Evolution

Leeches (Annelida: Hirudinea) possess powerful salivary anticoagulants and, accordingly, are frequently employed in modern, authoritative medicine. Members of the almost exclusively marine family Piscicolidae account for 20% of leech species diversity, and they feed on host groups (e.g., sharks) not encountered by their freshwater and terrestrial counterparts. Moreover, some species of Ozobranchidae feed on endangered marine turtles and have been implicated as potential vectors for the tumor-associated turtle herpesvirus. In spite of their ecological importance and unique host associations, there is a distinct paucity of data regarding the salivary transcriptomes of either of these families. Using next-generation sequencing, we profiled transcribed, putative anticoagulants and other salivary bioactive compounds that have previously been linked to blood feeding from 7 piscicolid species (3 elasmobranch feeders; 4 non-cartilaginous fish feeders) and 1 ozobranchid species (2 samples). In total, 149 putative anticoagulants and bioactive loci were discovered in varying constellations throughout the different samples. The putative anticoagulants showed a broad spectrum of described antagonistic pathways, such as inhibition of factor Xa and platelet aggregation, which likely have similar bioactive roles in marine fish and turtles. A transcript with homology to ohanin, originally isolated from king cobras, was found in Cystobranchus vividus but is otherwise unknown from leeches. Estimation of selection pressures for the putative anticoagulants recovered evidence for both positive and purifying selection along several isolated branches in the gene trees, and positive selection was also estimated for a few select codons in a variety of marine species. Similarly, phylogenetic analyses of the amino acid sequences for several anticoagulants indicated divergent evolution.

Are Fireworms Venomous? Evidence for the Convergent Evolution of Toxin Homologs in Three Species of Fireworms (Annelida, Amphinomidae)

Amphinomids, more commonly known as fireworms, are a basal lineage of marine annelids characterized by the presence of defensive dorsal calcareous chaetae, which break off upon contact. It has long been hypothesized that amphinomids are venomous and use the chaetae to inject a toxic substance. However, studies investigating fireworm venom from a morphological or molecular perspective are scarce and no venom gland has been identified to date, nor any toxin characterized at the molecular level. To investigate this question, we analyzed the transcriptomes of three species of fireworms-Eurythoe complanata, Hermodice carunculata, and Paramphinome jeffreysii-following a venomics approach to identify putative venom compounds. Our venomics pipeline involved de novo transcriptome assembly, open reading frame, and signal sequence prediction, followed by three different homology search strategies: BLAST, HMMER sequence, and HMMER domain. Following this pipeline, we identified 34 clusters of orthologous genes, representing 13 known toxin classes that have been repeatedly recruited into animal venoms. Specifically, the three species share a similar toxin profile with C-type lectins, peptidases, metalloproteinases, spider toxins, and CAP proteins found among the most highly expressed toxin homologs. Despite their great diversity, the putative toxins identified are predominantly involved in three major biological processes: hemostasis, inflammatory response, and allergic reactions, all of which are commonly disrupted after fireworm stings. Although the putative fireworm toxins identified here need to be further validated, our results strongly suggest that fireworms are venomous animals that use a complex mixture of toxins for defense against predators.

Hirudins and hirudin-like factors in Hirudinidae: implications for function and phylogenetic relationships

Haematophagous leeches express a broad variety of bioactive factors that are released from the salivary gland cells into the wound of a host during feeding. Among these, hirudin is probably the best studied factor and, moreover, the only one that has successfully made the transition from nature to clinical use. Many components of the leech saliva still remain either poorly characterized or yet completely unknown. Only recently, a new class of leech-derived factors has been discovered in Hirudo medicinalis, the hirudin-like factors (HLFs). HLFs comprise typical structural features of hirudin but lack others. We were able to verify the expression of HLFs not only in two additional species of the genus Hirudo, but also in Hirudinaria manillensis. Various phylogenetic analyses based on gene and protein sequences support a sister group relationship between hirudins and HLFs. Although potential molecular targets of HLFs remain unknown, the presence of multiple isoforms in individual leeches of different genera points to key functions in the regulation of several processes associated with the blood-sucking life style of leeches.

Host Matters: Medicinal Leech Digestive-Tract Symbionts and Their Pathogenic Potential

Digestive-tract microbiota exert tremendous influence over host health. Host-symbiont model systems are studied to investigate how symbioses are initiated and maintained, as well as to identify host processes affected by resident microbiota. The medicinal leech, Hirudo verbana, is an excellent model to address such questions owing to a microbiome that is consistently dominated by two species, Aeromonas veronii and Mucinivorans hirudinis, both of which are cultivable and have sequenced genomes. This review outlines current knowledge about the dynamics of the H. verbana microbiome. We discuss in depth the factors required for A. veronii colonization and proliferation in the leech crop and summarize the current understanding of interactions between A. veronii and its annelid host. Lastly, we discuss leech usage in modern medicine and highlight how leech-therapy associated infections, often attributable to Aeromonas spp., are of growing clinical concern due in part to an increased prevalence of fluoroquinolone resistant strains.

Comparative Transcriptomic Analyses of Three Species of Placobdella (Rhynchobdellida: Glossiphoniidae) Confirms a Single Origin of Blood Feeding in Leeches

One of the recalcitrant questions regarding the evolutionary history of clitellate annelids involves the feeding preference of the common ancestor of extant rhynchobdellid (proboscis bearing) and arhynchobdellid (jaw bearing) leeches. Whereas early evidence, based on morphological data, pointed towards independent acquisitions of blood feeding in the 2 orders, molecular-based phylogenetic data suggest that the ancestor of modern leeches was a sanguivore. Here, we use a comparative transcriptomic approach in order to increase our understanding of the diversity of anticoagulation factors for 3 species of the genus Placobdella, for which comparative data have been lacking, and inspect these in light of archetypal anticoagulant data for both arhynchobdellid and other rhynchobdellid species. Notwithstanding the varying levels of host specificity displayed by the 3 different species of Placobdella, transcriptomic profiles with respect to anticoagulation factors were largely similar -this despite the fact that Placobdella kwetlumye only retains a single pair of salivary glands, as opposed to the 2 pairs more common in the genus. Results show that 9 different anticoagulant proteins and an additional 5 putative antihemostasis proteins are expressed in salivary secretions of the 3 species. In particular, an ortholog of the archetypal, single-copy, anticoagulant hirudin (not previously available as comparative data for rhynchobdellids) is present in at least 2 of 3 species examined, corroborating the notion of a single origin of blood feeding in the ancestral leech.

Transcripts involved in hemostasis: Exploring salivary complexes from Haementeria vizottoi leeches through transcriptomics, phylogenetic studies and structural features

Throughout evolution, parasites have adapted in order to successfully intervene in the host defense, producing specific peptides and proteins. Interestingly, these peptides and proteins have been exploited as potential drug candidates against several diseases. Furthermore, biotechnology studies and cDNA libraries have remarkably contributed to identify potentially bioactive molecules. In this regard, herein, a cDNA library of salivary complexes from Haementeria vizottoi leeches was constructed, the transcriptome was characterized and a phylogenetic analysis was performed considering antistasin-like and antiplatelet-like proteins. Hundred twenty three transcripts were identified coding for putative proteins involved in animal feeding (representing about 10% of the expression level). These sequences showed similarities with myohemerythrins, carbonic anhydrases, anticoagulants, antimicrobials, proteases and protease inhibitors. The phylogenetic analysis, regarding antistasin-like and antiplatetlet-like proteins, revealed two main clades in the Rhynchobdellida leeches. As expected, the sequences from H. vizottoi have presented high similarities with those types of proteins. Thus, our findings could be helpful not only to identify new coagulation inhibitors, but also to better understand the biological composition of the salivary complexes.

More than just one: multiplicity of Hirudins and Hirudin-like Factors in the Medicinal Leech, Hirudo medicinalis

Blood-sucking leeches like the medicinal leech, Hirudo medicinalis, have been used for medical purposes since ancient times. During feeding, medicinal leeches transfer a broad range of bioactive substances into the host's wound to prevent premature hemostasis and blood coagulation. Hirudin is probably the best known of these substances. Despite its long history of investigation, recombinant production and clinical use, there still exist conflicting data regarding the primary structure of hirudin. Entirely unclear is the potential biological significance of three different subtypes and many isoforms of hirudins that have been characterized so far. Furthermore, there is only incomplete information on their cDNA sequences and no information at all on gene structures and DNA sequences are available in the databases. Our efforts to fill these gaps revealed the presence of multiple hirudin-encoding genes in the genome of Hirudo medicinalis. We have strong evidence for the expression of all three subtypes of hirudin within individual leeches and for the expression of additional hirudins or hirudin-like factors that may have different biological functions and may be promising candidates for new drugs.

Why do we study animal toxins?

Venom (toxins) is an important trait evolved along the evolutionary tree of animals. Our knowledges on venoms, such as their origins and loss, the biological relevance and the coevolutionary patterns with other organisms are greatly helpful in understanding many fundamental biological questions, i.e., the environmental adaptation and survival competition, the evolution shaped development and balance of venoms, and the sophisticated correlations among venom, immunity, body power, intelligence, their genetic basis, inherent association, as well as the cost-benefit and trade-offs of biological economy. Lethal animal envenomation can be found worldwide. However, from foe to friend, toxin studies have led lots of important discoveries and exciting avenues in deciphering and fighting human diseases, including the works awarded the Nobel Prize and lots of key clinic therapeutics. According to our survey, so far, only less than 0.1% of the toxins of the venomous animals in China have been explored. We emphasize on the similarities shared by venom and immune systems, as well as the studies of toxin knowledge-based physiological toxin-like proteins/peptides (TLPs). We propose the natural pairing hypothesis. Evolution links toxins with humans. Our mission is to find out the right natural pairings and interactions of our body elements with toxins, and with endogenous toxin-like molecules. Although, in nature, toxins may endanger human lives, but from a philosophical point of view, knowing them well is an effective way to better understand ourselves. So, this is why we study toxins.

Quo vadis venomics? A roadmap to neglected venomous invertebrates

Venomics research is being revolutionized by the increased use of sensitive -omics techniques to identify venom toxins and their transcripts in both well studied and neglected venomous taxa. The study of neglected venomous taxa is necessary both for understanding the full diversity of venom systems that have evolved in the animal kingdom, and to robustly answer fundamental questions about the biology and evolution of venoms without the distorting effect that can result from the current bias introduced by some heavily studied taxa. In this review we draw the outlines of a roadmap into the diversity of poorly studied and understood venomous and putatively venomous invertebrates, which together represent tens of thousands of unique venoms. The main groups we discuss are crustaceans, flies, centipedes, non-spider and non-scorpion arachnids, annelids, molluscs, platyhelminths, nemerteans, and echinoderms. We review what is known about the morphology of the venom systems in these groups, the composition of their venoms, and the bioactivities of the venoms to provide researchers with an entry into a large and scattered literature. We conclude with a short discussion of some important methodological aspects that have come to light with the recent use of new -omics techniques in the study of venoms.

A Polychaete's powerful punch: venom gland transcriptomics of Glycera reveals a complex cocktail of toxin homologs

Glycerids are marine annelids commonly known as bloodworms. Bloodworms have an eversible proboscis adorned with jaws connected to venom glands. Bloodworms prey on invertebrates, and it is known that the venom glands produce compounds that can induce toxic effects in animals. Yet, none of these putative toxins has been characterized on a molecular basis. Here we present the transcriptomic profiles of the venom glands of three species of bloodworm, Glycera dibranchiata, Glycera fallax and Glycera tridactyla, as well as the body tissue of G. tridactyla. The venom glands express a complex mixture of transcripts coding for putative toxin precursors. These transcripts represent 20 known toxin classes that have been convergently recruited into animal venoms, as well as transcripts potentially coding for Glycera-specific toxins. The toxins represent five functional categories: Pore-forming and membrane-disrupting toxins, neurotoxins, protease inhibitors, other enzymes, and CAP domain toxins. Many of the transcripts coding for putative Glycera toxins belong to classes that have been widely recruited into venoms, but some are homologs of toxins previously only known from the venoms of scorpaeniform fish and monotremes (stonustoxin-like toxin), turrid gastropods (turripeptide-like peptides), and sea anemones (gigantoxin I-like neurotoxin). This complex mixture of toxin homologs suggests that bloodworms employ venom while predating on macroscopic prey, casting doubt on the previously widespread opinion that G. dibranchiata is a detritivore. Our results further show that researchers should be aware that different assembly methods, as well as different methods of homology prediction, can influence the transcriptomic profiling of venom glands.

May salivary gland secretory proteins from hematophagous leeches (Hirudo verbana) reach pharmacologically relevant concentrations in the vertebrate host?

Saliva of hematophagous leeches (Hirudo sp.) contains bioactive proteins which allow the leech proper feeding and storage of ingested blood, but may also exert effects in the host. Leech therapy is used to treat many different ailments in humans, although only a small fraction of salivary proteins are characterized yet. Moreover, we do not know whether complete transfer of salivary proteins stored in the unicellular salivary glands in a leech to the host during feeding may generate concentrations that are sufficiently high to affect physiological processes in the host. Our 3D reconstruction of a portion of internal leech tissue from histological sections revealed that one leech contains approx. 37,000 salivary gland cells. Using tissue slices from pig liver and mouse skeletal muscle for reference, we obtained data for protein densities in leech salivary gland cells. As individual salivary cells are voluminous (67,000 µm³) and the stored proteins are densely packed (approx. 500 µg/mm³), we extrapolated that a single leech may contain up to 1.2 mg of salivary proteins. Analyzing protein extracts of unfed or fed leeches by 2D electrophoresis, we calculated the relative molar amounts of individual salivary proteins in the mass range of 17–60 kDa which may be released from a single leech during feeding. Distribution of these salivary proteins in the host (assumed plasma volume of 5 l) may result in concentrations of individual compounds between 3 and 236 pmol/l. Such concentrations seem sufficiently high to exert biochemical interactions with target molecules in the host.

Diversity and selective pressures of anticoagulants in three medicinal leeches (Hirudinida: Hirudinidae, Macrobdellidae)

Although medicinal leeches have long been used as treatment for various ailments because of their potent anticoagulation factors, neither the full diversity of salivary components that inhibit coagulation, nor the evolutionary selection acting on them has been thoroughly investigated. Here, we constructed expressed sequence tag libraries from salivary glands of two species of medicinal hirudinoid leeches, Hirudo verbana and Aliolimnatis fenestrata, and identified anticoagulant-orthologs through BLASTx searches. The data set then was augmented by the addition of a previously constructed EST library from the macrobdelloid leech Macrobdella decora. The identified orthologs then were compared and contrasted with well-characterized anticoagulants from a variety of leeches with different feeding habits, including non-sanguivorous species. Moreover, four different statistical methods for predicting signatures of positive and negative evolutionary pressures were used for 10 rounds each to assess the level and type of selection acting on the molecules as a whole and on specific sites. In total, sequences showing putative BLASTx-orthology with five and three anticoagulant-families were recovered in the A. fenestrata and H. verbana EST libraries respectively. Selection pressure analyses predicted high levels of purifying selection across the anticoagulant diversity, although a few isolated sites showed signatures of positive selection. This study represents a first attempt at mapping the anticoagulant repertoires in a comparative fashion across several leech families.

Disintegrins from hematophagous sources

Bloodsucking arthropods are a rich source of salivary molecules (sialogenins) which inhibit platelet aggregation, neutrophil function and angiogenesis. Here we review the literature on salivary disintegrins and their targets. Disintegrins were first discovered in snake venoms, and were instrumental in our understanding of integrin function and also for the development of anti-thrombotic drugs. In hematophagous animals, most disintegrins described so far have been discovered in the salivary gland of ticks and leeches. A limited number have also been found in hookworms and horseflies, and none identified in mosquitoes or sand flies. The vast majority of salivary disintegrins reported display a RGD motif and were described as platelet aggregation inhibitors, and few others as negative modulator of neutrophil or endothelial cell functions. This notably low number of reported disintegrins is certainly an underestimation of the actual complexity of this family of proteins in hematophagous secretions. Therefore an algorithm was created in order to identify the tripeptide motifs RGD, KGD, VGD, MLD, KTS, RTS, WGD, or RED (flanked by cysteines) in sialogenins deposited in GenBank database. The search included sequences from various blood-sucking animals such as ticks (e.g., Ixodes sp., Argas sp., Rhipicephalus sp., Amblyommasp.), tabanids (e.g., Tabanus sp.), bugs (e.g., Triatoma sp., Rhodnius prolixus), mosquitoes (e.g., Anopheles sp., Aedes sp., Culex sp.), sand flies (e.g., Lutzomyia sp., Phlebotomus sp.), leeches (e.g., Macrobdella sp., Placobdella sp.) and worms (e.g., Ancylostoma sp.). This approach allowed the identification of a remarkably high number of novel putative sialogenins with tripeptide motifs typical of disintegrins (>450 sequences) whose biological activity remains to be verified. This database is accessible online as a hyperlinked worksheet and displays biochemical, taxonomic, and gene ontology aspects for each putative disintegrin. It is also freely available for download (right click with the mouse) at links http://exon.niaid.nih.gov/transcriptome/RGD/RGD-Peps-WEB.xlsx (web version) and http://exon.niaid.nih.gov/transcriptome/RGD/RGD-sialogenins.zip (stand alone version).

Bacterial symbiont and salivary peptide evolution in the context of leech phylogeny

The evolutionary history of leeches is employed as a general framework for understanding more than merely the systematics of this charismatic group of annelid worms, and serves as a basis for understanding blood-feeding related correlates ranging from the specifics of gut-associated bacterial symbionts to salivary anticoagulant peptides. A variety of medicinal leech families were examined for intraluminal crop bacterial symbionts. Species of Aeromonas and Bacteroidetes were characterized with DNA gyrase B and 16S rDNA. Bacteroidetes isolates were found to be much more phylogenetically diverse and suggested stronger evidence of phylogenetic correlation than the gammaproteobacteria. Patterns that look like co-speciation with limited taxon sampling do not in the full context of phylogeny. Bioactive compounds that are expressed as gene products, like those in leech salivary glands, have 'passed the test' of evolutionary selection. We produced and bioinformatically mined salivary gland EST libraries across medicinal leech lineages to experimentally and statistically evaluate whether evolutionary selection on peptides can identify structure-function activities of known therapeutically relevant bioactive compounds like antithrombin, hirudin and antistasin. The combined information content of a well corroborated leech phylogeny and broad taxonomic coverage of expressed proteins leads to a rich understanding of evolution and function in leech history.