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.

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.

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.

Leeches are well known for blood-feeding habits and are widely used for medicinal purposes as they secrete various antithrombotic substances. However, some leeches exhibit non-hematophagous habits and their significance for medicinal use is controversial. Here we provide the chromosome-level genomes of two non-hematophagous leeches, Whitmania acranulata and Whitmania laevis, and, in combination with previous results from Whitmania pigra, we compared these genomes with an emphasis on antithrombotic biomolecules. All three species had the same chromosome number of 11. The genome size, repeat site percentage, and number of protein-coding genes of W. laevis (173.87 Mb, 28.28%, 23,818) were similar to those of W. pigra (169.37 Mb, 27.02%, 24,156), whereas these values of W. acranulata (181.72 Mb, 29.55%, 27,069) were higher than those of the other two leeches. W. laevis was a monophyletic clade of W. pigra, whereas W. acranulata had a paraphyletic relationship with W. pigra. The number of antithrombotic genes in W. laevis (N = 76) was similar to that of W. pigra (N = 79), whereas W. acranulata (N = 102) had apparently more such genes. Of the 21 gene families, 9 and 11 were differentially expressed in W. acranulata and W. laevis compared to W. pigra, respectively. The expression profiles of the antithrombotic gene families were more similar between W. acranulata and W. laevis. Although there were several cases of gene loss or pseudogenization, most antithrombotic genes of the three Whitmania species were intact and transcribable. These results provide valuable insights into the evolution of non-hematophagous leeches and development of antithrombotic drugs.

Vitellogenin (Vg), a storage protein, has been significantly studied for its egg yolk precursor role in oviparous animals. Recent studies found that vitellogenin and its Vg‐like homologs were fundamentally involved in many other biological processes in social insects such as female caste differences and oxidative stress resilience. In this study, we conducted the first large‐scale molecular evolutionary analyses of vitellogenin coding genes (Vg) and Vg‐like genes of bumble bees, a primitively eusocial insect belonging to the genus Bombus. We obtained sequences for each of the four genes (Vg, Vg‐like‐A, Vg‐like‐B, and Vg‐like‐C) from 27 bumble bee genomes (nine were newly sequenced in this study), and sequences from the two closest clades of Bombus, including five Apis species and five Tetragonula species. Our molecular evolutionary analyses show that in bumble bee, the conventional Vg experienced strong positive selection, while the Vg‐like genes showed overall relaxation of purifying selection. In Apis and Tetragonula; however, all four genes were found under purifying selection. Furthermore, the conventional Vg showed signs of strong positive selection in most subgenera in Bombus, apart from the obligate parasitic subgenus Psithyrus which has no caste differentiation. Together, these results indicate that the conventional Vg, a key pleiotropic gene in social insects, is the most rapidly evolving copy, potentially due to its multiple known social functions for both worker and queen castes. This study shows that concerted evolution and purifying selection shaped the evolution of the Vg gene family following their ancient gene duplication and may be the leading forces behind the evolution of new potential protein function enabling functional social pleiotropy.