ABSTRACT
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.
ABSTRACT
Whitmania pigra is a unique, fluid-sucking ectoparasite and an anticoagulant medical leech. The codon usage bias (CUB) is the nonuniform usage of synonymous codons in which some codons are more preferred than others. Here, we performed a comprehensive analysis of CUB of genes in W. pigra, analyzing 140 780 transcripts, 59 553 unigenes, and 20 304 qualified coding sequences (CDSs) from the transcriptomic data of W. pigra. The effective number of codons values suggested that the CUB was low in these genes. We recognized profoundly favored codons in W. pigra that have a G/C-ending. Parity rule two-bias plots suggested that both mutation pressure and natural selection might have influenced the CUB. However, neutrality plots revealed that natural selection might have played a major role while mutation pressure might have played a minor role in shaping the CUB. We applied principal component analysis to relative synonymous codon usage values for divided CDSs based on GC content and codon-ending bases. Codon usage in W. pigra had a general inclination toward C-ending codons and natural selection rather than mutation pressure is the dominant force in the genetic evolution of W. pigra. To our knowledge, this is the first study to describe a complete codon usage analysis of W. pigra; this will increase the understanding of CUB and evolution in W. pigra. The analysis of codon usage patterns in W. pigra aids in understanding its evolution and genetic architecture.
Subject(s)
Codon Usage , Evolution, Molecular , Leeches/genetics , Mutation , Selection, Genetic , AnimalsABSTRACT
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.
