Complete Mitogenome of the : Comparative Analysis and Phylogenetic Relationships among the Members of .

In the last decade, the phylogenetic relationships within the genus have become controversial, due to a lack of molecular data. The mitochondrial genome plays a vital role in the reconstruction of phylogenetic relationships and in revealing the molecular evolution of bony fishes. Herein, we obtained the complete mitogenome of via HiFi reads of the Pacbio Sequel II system and DNBSEQ short-reads. We compared all available mitogenomes of the genus and reconstructed the phylogeny of Nemacheilidae, based on the mitogenomes, using maximum likelihood (ML) methods. The results show that the complete mitogenome sequence of was circular and 16,568 bp in length, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA), 2 ribosomal RNA (rRNA), and a typical control region (D-loop). The most common start codons were ATG, except for , and TAA/TAG were the stop codons for all PCGs. In total, 677 SNPs and 9 INDELs have been found by comparing the sequence divergence between this study and previous reports. Purity selection was found in all PCGs. Phylogeny was inferred by analyzing the 13 PCGs and the concatenated nucleotide sequences of 30 mitogenomes. The phylogenetic analyses based on the nucleotides of the 13 PCGs supported the assumption that the genus can be divided into 4 main clades and demonstrated that and are closely related species for the first time. This study laid the foundation for further study on the mitogenome and phylogeny of Nemacheilidae fishes.