Torix tukubana is a poorly understood proboscidate leech species, generally an ectoparasite on amphibian species. In this study, the complete mitochondrial genome (mitogenome) of T. tukubana was sequenced using next-generation sequencing (NGS), and the essential characteristics, gene arrangement, and phylogenetic relationship were analyzed. The results showed that the T. tukubana mitogenome was 14,814 bp in length, consisting of 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, and 1 control region (CR). The mitogenome composition presented a strong A + T bias (73.6%). All tRNAs had the typical clover structure except the trnS1 (TCT), whose dihydrouridine (DHU) arm was short, having only one complementary base pair. Additionally, 8 gene order patterns were identified among 25 known Hirudinea species, and T. tukubana was identical to the Hirudinea ground pattern. A phylogenetic analysis based on 13 PCGs indicated that all the studied species clustered into three main clades. The relationships among Hirudinea species were basically consistent with their gene arrangement results, but different from their morphological taxonomy. T. tukubana was in the monophyletic group of Glossiphoniidae, a finding consistent with previous research. Our results provided the essential characteristics of the T. tukubana mitogenome. As the first complete mitogenome of Torix, it could offer valuable information for a systematic understanding of the Hirudinea species.
Keywords: Torix tukubana; comparative analyses; leech; mitochondrial genome; phylogenetics.