Draft Genome of the Asian Buffalo Leech Hirudinaria manillensis

Comparative genomics of two Asian medicinal leeches Hirudo nipponia and Hirudo tianjinensis: With emphasis on antithrombotic genes and their corresponding proteins

Leeches secrete various biologically active substances which have important medical and pharmaceutical values in antithrombotic treatments. Here, we provide a high quality genome of two Asian medicinal leeches Hirudo nipponia and Hirudo tianjinensis, based on which, we identified 22 antithrombotic gene families, including fourteen coagulation inhibitors, four platelet aggregation inhibitors, three fibrinolysis enhancers, and one tissue penetration enhancer. The total numbers of antithrombotic genes were similar between H. nipponia (N = 86) and H. tianjinensis (N = 83). Molecular evolution analysis showed that no significant differences were detected between the two species in any of the three selection indices (dN, dS, and dN/dS), nor in the number of sites under positive/purifying selection. RNA-Seq based gene expression analysis showed that the overall expression patterns of the antithrombotic gene families were not significantly deviated between the two species. Our results indicated that there were rather close similarities between the two leeches on genomic characteristics, especially for the molecular evolution and expression of antithrombotic genes. Our study provides the most comprehensive collection of antithrombotic biomacromolecules from the two Asian medicinal leeches to date. These results will greatly facilitate the research and application of leech derivatives for medical and pharmaceutical purposes of thrombosis.

A Chromosome-Level Genome Assembly of the Non-Hematophagous Leech Whitmania pigra (Whitman 1884): Identification and Expression Analysis of Antithrombotic Genes

Despite being a non-hematophagous leech, Whitmania pigra is widely used in traditional Chinese medicine for the treatment of antithrombotic diseases. In this study, we provide a high quality genome of W. pigra and based on which, we performed a systematic identification of the potential antithrombotic genes and their corresponding proteins. We identified twenty antithrombotic gene families including thirteen coagulation inhibitors, three platelet aggregation inhibitors, three fibrinolysis enhancers, and one tissue penetration enhancer. Unexpectedly, a total of 79 antithrombotic genes were identified, more than a typical blood-feeding Hirudinaria manillensis, which had only 72 antithrombotic genes. In addition, combining with the RNA-seq data of W. pigra and H. manillensis, we calculated the expression levels of antithrombotic genes of the two species. Five and four gene families had significantly higher and lower expression levels in W. pigra than in H. manillensis, respectively. These results showed that the number and expression level of antithrombotic genes of a non-hematophagous leech are not always less than those of a hematophagous leech. Our study provides the most comprehensive collection of antithrombotic biomacromolecules from a non-hematophagous leech to date and will significantly enhance the investigation and utilization of leech derivatives in thrombosis therapy research and pharmaceutical applications.

Revisiting the Asian Buffalo Leech (Hirudinaria manillensis) Genome: Focus on Antithrombotic Genes and Their Corresponding Proteins

Leeches are well-known annelids due to their obligate blood-feeding habits. Some leech species secrete various biologically active substances which have important medical and pharmaceutical value in antithrombotic treatments. In this study, we provided a high-quality genome of the Asian buffalo leech (Hirudinaria manillensis), based on which we performed a systematic identification of potential antithrombotic genes and their corresponding proteins. Combining automatic and manual prediction, we identified 21 antithrombotic gene families including fourteen coagulation inhibitors, three platelet aggregation inhibitors, three fibrinolysis enhancers, and one tissue penetration enhancer. A total of 72 antithrombotic genes, including two pseudogenes, were identified, including most of their corresponding proteins forming three or more disulfide bonds. Three protein families (LDTI, antistasin, and granulin) had internal tandem repeats containing 6, 10, and 12 conserved cysteines, respectively. We also measured the anticoagulant activities of the five identified hirudins (hirudin_Hman1 ~ hirudin_Hman5). The results showed that three (hirudin_Hman1, hirudin_Hman2, and hirudin_Hman5), but not the remaining two, exhibited anticoagulant activities. Our study provides the most comprehensive collection of antithrombotic biomacromolecules from a leech to date. These results will greatly facilitate the research and application of leech derivatives for medical and pharmaceutical purposes in the treatment of thrombotic diseases.

Cloning and functional identification of pmKPI cDNA in Poecilobdella manillensis

BACKGROUND

Kazal-type serine protease inhibitors play a role in physiological processes such as blood coagulation and fibrinolysis. The amino acid residues at the P1 site are different, and they inhibit different types of proteases. The inhibitory mechanism of the protease in the salivary glands of Poecilobdella manillensis is still unclear.

METHODS AND RESULTS

Based on cloning, prokaryotic expression and bioinformatics analysis, we studied the role of Kazal-type serine protease inhibitors in P. manillensis and analyzed their expression by quantitative real-time PCR. The results suggested that the recombinant protein was successfully expressed in the supernatant when a prokaryotic expression vector was constructed and induced with 0.2 mmol/L IPTG at 37 °C for 4 h, and the enzymatic activity was determined. The mature protein encodes 91 amino acids and has a relative molecular weight of 9929.32 Da, and after removing the signal peptide, the theoretical isoelectric point was 8.79. It is an unstable protein without a transmembrane domain. The mature protein contains two Kazal-type domains, in which all P1 residues are Lys, consisting of an α helix and three antiparallel β sheets. The upregulated expression of the mRNA was induced after a meal was provided, and the results showed an increasing and then decreasing trend.

CONCLUSIONS

Taken together, the results indicate that mature proteins from P. manillensis inhibit thrombin activity, laying the foundation for the subsequent in-depth study of the function of genes encoding Kazal-type serine protease inhibitors.

Molecular mechanisms underlying hematophagia revealed by comparative analyses of leech genomes

BACKGROUND

Leeches have been used in traditional Chinese medicine since prehistoric times to treat a spectrum of ailments, but very little is known about their physiological, genetic, and evolutionary characteristics.

FINDINGS

We sequenced and assembled chromosome-level genomes of 3 leech species (bloodsucking Hirudo nipponia and Hirudinaria manillensis and nonbloodsucking Whitmania pigra). The dynamic population histories and genome-wide expression patterns of the 2 bloodsucking leech species were found to be similar. A combined analysis of the genomic and transcriptional data revealed that the bloodsucking leeches have a presumably enhanced auditory sense for prey location in relatively deep fresh water. The copy number of genes related to anticoagulation, analgesia, and anti-inflammation increased in the bloodsucking leeches, and their gene expressions responded dynamically to the bloodsucking process. Furthermore, the expanded FBN1 gene family may help in rapid body swelling of leeches after bloodsucking, and the expanded GLB3 gene family may be associated with long-term storage of prey blood in a leech's body.

CONCLUSIONS

The high-quality reference genomes and comprehensive datasets obtained in this study may facilitate innovations in the artificial culture and strain optimization of leeches.

Make it double: identification and characterization of a Tandem-Hirudin from the Asian medicinal leech Hirudinaria manillensis

Haematophagous leeches express a broad variety of secretory proteins in their salivary glands, among them are hirudins and hirudin-like factors. Here, we describe the identification, molecular and initial functional characterization of Tandem-Hirudin (TH), a novel salivary gland derived factor identified in the Asian medicinal leech, Hirudinaria manillensis. In contrast to the typical structure of hirudins, TH comprises two globular domains arranged in a tandem-like orientation and lacks the elongated C-terminal tail. Similar structures of thrombin inhibitors have so far been identified only in kissing bugs and ticks. Expression of TH was performed in both cell-based and cell-free bacterial systems. A subsequent functional characterization revealed no evidence for a thrombin-inhibitory potency of TH.

Comprehensive Transcriptome Sequencing Analysis of Hirudinaria manillensis in Different Growth Periods

Medical leeches are widely been used in biochemical and clinical medical studies, helping to restore blood circulation to grafted or severely injured tissue. Mostly, adult leeches are being used in the traditional pharmacopeia, but the gene expression profiling of leeches in different growth periods is not well-reported. So, in this study, we used transcriptome analysis to analyze the comparative gene expression patterns of Hirudinaria manillensis (H. manillensis) in different growth periods, including larval, young, and adult stages. We constructed 24 cDNA libraries from H. manillensis larval, young, and adult stages, and about 54,639,118 sequences were generated, 18,106 mRNA transcripts of which 958 novel mRNAs and 491 lncRNAs were also assembled as well. Furthermore, the results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the differentially upregulated genes from the larval to adult stages were enriched in pathways such as cilium, myofibril, contractile fiber, cytoskeleton proteins, dilated cardiomyopathy, adrenergic signaling in cardiomyocytes, etc. Moreover, in the adult stages, a significant increase in the expression of the Hirudin-HM (HIRM2) genes was detected. In addition, our comparative transcriptome profiling data from different growth stages of H. manillensis also identified a large number of DEGs and DElncRNAs which were tentatively found to be associated with the growth of H. manillensis; as it grew, the muscle-related gene expression increased, while the lipid metabolism and need for stimulation and nutrition-related genes decreased. Similarly, the higher expression of HIRM2 might attribute to the high expression of protein disulfide isomerase gene family (PDI) family genes in adulthood, which provides an important clue that why adult leeches rather than young leeches are widely used in clinical therapeutics and traditional Chinese medicine.

Affinity ultrafiltration and UPLC-HR-Orbitrap-MS based screening of thrombin-targeted small molecules with anticoagulation activity from Poecilobdella manillensis

This study aims to screen potential anticoagulant components from leeches, a representative animal-sourced traditional Chinese medicine using thrombin (THR)-targeted ultrafiltration combined with ultrahigh performance liquid chromatography and high-resolution Orbitrap mass spectrometry (UPLC-HR-Orbitrap-MS). As a result, five small molecules in leech extract were discovered to interact with THR for the first time. Among them, two new compounds were isolated and their structures were identified by IR, HR-MS and NMR data. Furthermore, their THR inhibitory activity was confirmed with IC₅₀ values of 4.74 and 8.31 μM, respectively. In addition, molecular docking analysis showed that the active (catalytic) site of THR might be the possible binding site of the two hits. Finally, reverse screening analysis indicated that LTA4-H, ACE and ALOX5AP were potential anticoagulant targets of the two new compounds. This study will broaden our understanding of the medicinal substance basis in leeches and further contribute to the discovery and development of clinical anticoagulant drugs from leeches.

High-quality genome assembly of Metaphire vulgaris

Earthworms enrich the soil and protect the health of their ecological environment. Previous studies on these invertebrates determined their protein content, hormone secretions, medicinal value, and ecological habits, but their whole genomic sequence remains incomplete. We performed whole genome sequencing of Metaphire vulgaris (Chen, 1930), which belongs to the genus Metaphire of the family Megascolecidae. The genome assembly was 729 Mb, with a N50 contig size of 4.2 Mb. In total, 559 contigs were anchored to 41 chromosomes according to the results of Hi-C (High-throughput Chromosome Conformation Capture) technology, which was confirmed by karyological analysis. A comparison of the genomic sequences and genes indicated that there was a whole-genome duplication in M. vulgaris followed by several chromosome fusion events. Hox genes and lumbrokinase genes were identified as partial clusters surrounding the genome. Our high-quality genome assembly of M. vulgaris will provide valuable information for gene function and evolutionary studies in earthworms.

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.