The complete mitochondrial genome sequence of Whitmania pigra (Annelida, Hirudinea): the first representative from the class Hirudinea

The mitochondrial genome is a significant tool for investigating the evolutionary history of metazoan animals. The currently available mitochondrial genome data in GenBank is limited to understand the detail evolutionary relationship among the metazoan animals, especially in the phylum Annelida. Here we present the mitochondrial gene organization, gene order and codon usage of the leech Whitmania pigra (Annelida), which is the first representative from the class Hirudinea. It is a circular molecule of 14,426bp, and encodes 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. All 37 genes of W. pigra mitochondrial genome are transcribed from the same strand, which is identical to studied annelids, two echiurans, two sipunculans and many other lophotrochozoans. Five conserved gene clusters can be found in mitochondrial genomes of nine studied annelids, including (1) cox1-N-cox2; (2) cox3-Q-nad6-cob-W-atp6; (3) H-nad5-F-E-P-T-nad4L-nad4; (4) srRNA-V-lrRNA; and (5) nad3-S(1)-nad2. Compared with that of other studied annelids, translocations of transfer RNAs were found in the gene arrangement of W. pigra mitochondrial genome. Phylogenetic analysis strongly support that the species from Hirudinina and Oligochaeta form a monophyletic group Clitellata (BPM=100, BPP=100), which is consistent with previous research based on morphological and other molecular data. Both gene order data and amino acid sequences reveal that echiurans are derived annelids and sipunculans should be clustered with annelids and echiurans.