Comparative mitochondrial genome analysis of Cynoglossidae (Teleost: Pleuronectiformes) and phylogenetic implications

Generally, a teleostean group (e.g., family or genus) owns one type or a set of similar mitochondrial gene arrangement. It is interesting, however, that four different types of gene arrangement have been found in the mitochondrial genome (mitogenome) of Cynoglossidae species. So far, the possible mechanisms of mitogenomic gene rearrangement and its potential implications have aroused widespread attention and caused lots of controversy. Here, a total of 21 Cynoglossidae mitogenomes and a newly sequenced mitogenome of Cynoglossus puncticpes (Pleuronectiformes: Cynoglossidae) were compared. The length ranges from 16 417 bp to 18 369 bp, which is mainly caused by the length heteroplasmy of control region (CR). Further analysis reveals that the difference of tandem repeats acts as a determining factor resulting in the length heterogeneity. Like most gene rearrangements of Cynoglossinae mitogenomes, tRNA-Gln gene encoded by the L-strand has translocated to the H-strand (Q inversion), accompanied by the translocation of CR in C. puncticpes mitogenome. The typical IQM order (tRNA-Ile-Gln-Met) changed to QIM order. Tandem duplication/random loss and mitochondrial recombination were accepted as the most possible models to account for the rearrangements in C. puncticpes mitogenome. Phylogenetic trees showed a strong correlation between the gap spacer in the rearranged QIM area and phylogeny, which provides a fresh idea for phylogenetic studies in future.