Comparative mitogenome analyses uncover mitogenome features and phylogenetic implications of the subfamily Cobitinae
BMC Genomics volume 22, Article number: 50 (2021)
Loaches of Cobitinae, widely distributed in Eurasian continent, have high economic, ornamental and scientific value. However, the phylogeny of Cobitinae fishes within genera or family level remains complex and controversial. Up to now, about 60 Cobitinae mitogenomes had been deposited in GenBank, but their integrated characteristics were not elaborated.
In this study, we sequenced and analyzed the complete mitogenomes of a female Cobits macrostigma. Then we conducted a comparative mitogenome analysis and revealed the conserved and unique characteristics of 58 Cobitinae mitogenomes, including C. macrostigma. Cobitinae mitogenomes display highly conserved tRNA secondary structure, overlaps and non-coding intergenic spacers. In addition, distinct base compositions were observed among different genus and significantly negative linear correlation between AT% and AT-skew were found among Cobitinae, genus Cobitis and Pangio mitogenomes, respectively. A specific 3 bp insertion (GCA) in the atp8-atp6 overlap was identified as a unique feature of loaches, compared to other Cypriniformes fish. Additionally, all protein coding genes underwent a strong purifying selection. Phylogenetic analysis strongly supported the paraphyly of Cobitis and polyphyly of Misgurnus. The strict molecular clock predicted that Cobitinae might have split into northern and southern lineages in the late Eocene (42.11 Ma), furthermore, mtDNA introgression might occur (14.40 Ma) between ancestral species of Cobitis and ancestral species of Misgurnus.
The current study represents the first comparative mitogenomic and phylogenetic analyses within Cobitinae and provides new insights into the mitogenome features and evolution of fishes belonging to the cobitinae family.
Vertebrate mitogenome is a small (16–17 kb) and circular double-stranded molecule . It contains 37 genes including 22 tRNA genes, 13 PCGs and two rRNA genes . It also has two noncoding regions, OL and CR, and the latter contains regulatory elements for controlling the transcription and replication of mtDNA molecule [2, 3]. Due to its unique features, such as high copy numbers in tissues, simple genomic organization, maternal inheritance, almost unambiguous orthology, haploid inheritance and high nucleotide substitution rate [4,5,6], mitogenome has been widely applied in species identification, i.e., DNA barcoding, as well as population genetics, conservation biology, molecular phylogenetics and evolutionary processes [7,8,9,10,11,12,13]. Gene arrangements of fish mitogenomes are generally conserved, only with a few exceptions . However, the genome sequence length, the bias of base composition and start/stop codon, the overlap and IGSs are diverse among different species .
Cobitinae is a subfamily of Cobitidae that was first identified by Hora (1932). To date, it contains 214 species recorded in FishBase, covering 21 genera, such as Cobits, Misgurnus and Paramisgurnus . Loaches of subfamily Cobitinae are bottom-dwelling fishes and widely distributed in Eurasian continent. They usually possess high economic, ornamental and scientific research value. Loach commercial farming, including cobitid loach (M. anguillicaudatus) and large-scale loach (P. dabryanus), occupies a significant position in freshwater aquaculture of Asia, due to their enjoyable taste, high nutritional value, rapid growth and strong adaptation [16,17,18]. In China, loach is used as a diet therapy or folk remedy for patient’s recovery or treatment of many diseases, such as hepatitis, osteomyeitis, carbuncles, and cancers. Many Cobitis populations are mixed diploid-polyploid, even bisexual and unisexual forms co-existing in the same niche [19,20,21]. They are suitable as models to reveal the relationship among hybridization, polyploidization, reproduction, speciation and evolution [21,22,23]. Due to their great diversity, they are also used to trace the biogeographic history of freshwater systems and to reflect geologic events . Cobitinae fishes usually inhabit various benthic habitats in rivers, lakes, streams and ponds . However, dilapidation of the ecological environment has led to a decrease of benthic organisms [26, 27]. Cobitinae fishes are seriously threatened and their wild populations are gradually decreasing . On this account, the diversity of these benthic fishes have been used as a bioindicator to assess the quality of the ecological environment [29, 30]. In addition, many Cobitinae species, such as the “kuhli loaches”, are well-known in Southeast Asia and Europe as ornamental fish for their varied morphological patterns and the ability to ingest bottom organic residues.
Cobitinae fishes are difficult to be classified because of their morphological similarity and high plasticity in morphology . Although the secondary sexual dimorphism is used to define genera, it is not always congruent with the current genera definitions. The molecular phylogeny of Cobitinae fishes has been studied at the genera or family level via one or two mitochondrial and/or nuclear genes [24, 31,32,33,34,35,36], and remains complex and controversial. For example, based on mitochondrial gene cytb and nuclear gene rag-1, Perdices et al. (2016)  reconstructed the phylogenetic relationship of Northern Clade of family Cobitidae that inhabit in Europe, and North and Northwest parts of Asia. The subfamily Cobitinae was divided into Cobitis sensu lato group (Cobitis, Iksookimia, Niwaella and Kichulchoia), Misgurnus sensu lato group (Misgurnus, Paramisgurnus and Koreocobitis), Microcobitis, and Sabanejewia. Although the monophyly of the groups were resolved, the relationships within the groups are discordant with current taxonomic status.
Up to now, about 60 mitogenomes, covering more than 40 species of Cobitinae, have been deposited into GenBank [38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55]. Although a few mitogenomes characteristics were described, the integrated characteristics of Cobitinae mitogenomes are still not well known. In this study, we sequenced the mitogenome of C. macrostigma, the type species of the genus Cobitis , and compared it with other 41 species (57 individuals) to amplify detailed features of the Cobitinae mitogenomes. Additionally, we assembled a large sequence matrix (11,442 bp) of 58 Cobitinae mitogenomes and two outgroups to investigate the phylogenetic status and the origin time of Cobitinae fishes.
General features of C. macrostigma mitogenome
The mitogenome of C. macrostigma was sequenced, annotated and compared with 57 Cobitinae mitogenomes (Table 1). It contains 13 PCGs (nd1–6, nd4l, cox1–3, cytb, atp6 and atp8), 22 tRNA genes, two rRNA genes (12S rRNA and 16S rRNA) and two non-coding regions (OL and CR) (GenBank: MT259034). Gene order and orientation are same to most teleost mitogenomes (Fig. 1, Table 2). PCGs range from 168 bp (atp8) to 1551 bp (cox1) in size, with a total length of 11,427 bp. tRNAs vary from 66 bp (tRNACys(C)) to 76 bp (tRNALys(K)) in size, with a total length of 1557 bp. The length of small encoding subunit 12S rRNA and large subunit 16S rRNA are 952 bp and 1675 bp, respectively. They are flanked by tRNAPhe and tRNALeu(UUR) and interposed by tRNAVal. Among 58 mitogenomes analyzed, the entire mitogenome of C. macrostigma has the highest (99.6%) similarity with C. granoei and lowest (88.2%) with C. sinensis.
Highly conserved tRNAs secondary structure, overlaps and non-coding intergenic spacers among Cobitinae mitogenomes
Cobitinae mitogenomes range from 16,337 bp (L. annandalei) to 16,647 bp (M. anguillicaudatus and C. takatsuensis) in length (Table 1). Their gene composition, gene arrangement and strand bias are highly conserved (Fig.1 and Table 2). Among the 22 tRNAs, due to the absence of DHU arm, tRNAser(AGN) (S1) is the only one that is not folded into the typical clover-leaf secondary structure (Fig. 2a). In the Cobitinae mitogenomes, unmatched base pairs are widespread among tRNAs. Taking C. macrostigma as an example, there are 446 base pairs among the 22 tRNAs, and only one gene (tRNALeu(CUN)) possesses a fully paired stem. In the 425 base pairs of other 21 tRNAs, there are 43 (10.1%) unmatched base pairs that contain 28 noncanonical matches of G-U and 15 other mismatches, including A-C (7), A-A (1), C-C (2), C-U (2), and U-U (3) (Fig. 2a). Most of them are located in the acceptor, DHU and anticodon stems.
We also compared the gene overlaps and IGSs among 58 Cobitinae mitogenomes. Two long overlaps (atp8-atp6 and nd4l-nd4) and two long IGSs (OL and tRNAAsp-cox2) were found in Cobitinae mitogenomes. Highly conserved motifs “ATGCTAA” and “ATGGCAATAA” were found in the overlapped junctions between nd4l and nd4, and between atp8 and atp6, respectively (Fig. 3a). There are also several small overlaps between adjacent tRNA genes, such as tRNAIle - tRNAGln and tRNAThr - tRNAPro. OL is located within the five gene cluster (WANCY) (Table 2, Fig.1) and its secondary structure shows a stable stem-loop hairpin, which is strengthened by six C-G base pairs (Fig. 2b). Among the 31 bp of OL, the C-G base pairs on stems are highly conserved while the loops in the middle are variable (Fig. 3b). Another long IGS, between tRNAAsp and cox2, is also conserved in the 5′ and 3′ end, and highly variable in the middle.
CR, located between tRNAPro and tRNAPhe, is the most variable region in Cobitinae mitogenomes and ranges from 872 bp (Lepidocephalus macrochir) to 990 bp (C. takatsuensis) (Supplementary Table 2) . Three domains are conserved and can be recognized in Cobitinae mitogenomes (Fig. 3c). They are terminal associated sequences (TAS), the central conserved-blocks (CSB-D, CSB-E and CSB-F) and conserved sequence blocks (CSB-1, CSB − 2 and CSB-3).
Usage bias of start and stop codon, codon distributions and relative synonymous codons in Cobitinae mitogenomes
The typical start codon ATG is conservative and is used in 12 PCGs, while GTG is only used in cox1 in 98% (57/58) analyzed Cobitinae mitogenomes except one individual of M. anguillicaudatus (No. 11) (Fig. 4, Supplementary Table 3). Five types of stop codons were found, containing three canonical (TAA, TAG and AGA) and two truncated stop codons (TA- and T--) (Fig. 4). The two truncated termination codons are used in nd2, cox2, atp6, cox3, nd3, nd4 and cytb, the 3′ -ends of which are followed by a tRNA gene encoded with the same strand.
The codon distribution and relative synonymous codon usage (RSCU) of 58 Cobitinae mitogenomes were analyzed. Our results show that codon distribution is largely coincident among these Cobitinae mitogenomes (Supplementary Figure S1). As shown by six representative species of Cobitinae, the codons encoding Leu(CUN), Ala and Thr are the three most frequently present, while those encoding Cys are rare (Fig. 5a). Compared to the other five Cobitinae species, P. anguillaris uses more codons of Leu(CUN) and less codons of Leu(UUR). The patterns of RSCU are also consistent among the analyzed species (Fig. 5b). Degenerated codons are biased to use more A/T than G/C in the 3rd position of PCGs, which results in the content of A + T is higher than G + C in the 3rd position of Cobitinae PCGs. For example, the codons for Arginine CCA and the codes for Tryptophan UGU are prevalent, while their other synonymous codons are relatively less used.
A + T %, AT-skew and their linear correlations of Cobitinae mitogenomes
The A + T content and AT-skew of whole mitogenomes, PCGs, tRNAs, rRNAs and CR were calculated (Fig. 6a-b). The 58 Cobitinae mitogenomes all exhibit AT bias, and the A + T content is the lowest (54.8 ± 0.6%) in tRNAs and the highest (66.3 ± 0.9%) in CR (Fig. 6a, Supplementary Table 2). The AT-skew values are the largest and positive in rRNAs, while they are the smallest in PCGs and most are negative except Canthophrys gongota, Acantopsis choirorhynchos, P. cuneovirgata, P. kuhlii, P. oblonga, and Kottelatlimia pristes (Fig. 6, Supplementary Table 2). These results indicate that PCGs are biased towards using T not A in most Cobitinae mitogenomes. To examine whether the A + T content and AT-skew are different in three codon position of PCGs, we also selected the six Cobitinae species for a more detailed analysis. The A + T content shows 1st < 2nd <3rd in the three position of PCGs in all analyzed fishes. Meanwhile, the AT-skew of 1st and 3rd are positive while 2nd is negative (Table 3). This is due to the bias usage of relative synonymous codons (Fig. 5b). In all analyzed Cobitinae mitogenomes, CRs possess more A and C with all AT-skew values positive (0.002–0.112) and GC-skew negative (− 0.341−− 0.101) (Supplementary Table 2).
The correlations of Cobitinae mitogenomes (yA1 = − 0.0166x – 0.9047, R2 = 0.5991) genus Cobits (yA2 = − 0.012x + 0.5786, R2 = 0.5197) and Pangio (yA3 = = − 0.0466x + 2.5813, R2 = 0.5486) were calculated between A + T % versus AT-skew. All of them showed negative linear correlations, implying that AT-skew becomes more positive with the increasing of A + T content (Fig. 6c). The similar negative linear correlations were also found in G + C % versus GC-skew (Fig. 6d).
Non-synonymous and synonymous substitutions
To better understand the role of selective pressure and evolutionary relations of Cobitinae fishes, the ω or dN/dS value of each PCG was calculated (Fig. 7). All the PCGs evolved under a purifying selection (ω < 0.5). The atp8 gene showed the highest ω value (ω = 0.12) and the cox family genes were lowest (ω = 0.02 ± 0.01). This phenomenon is also found in most Metazoa , but the fold change (> 10 fold) is particularly high in Cobitinae. The lower ω value represents less variations in amino acids. Thus, cox1, cox3 and cytb are potential barcoding markers for Cobitinae species identification.
Phylogenetic analysis of Cobitinae fishes
Molecular phylogenetic analyses were performed using 13 PCGs from 58 Cobitinae mitogenomes, belonging to 41 species from 14 genera. The ML and BI analyses generated similar topology with high bootstrap support / posterior probability values. Each tree was similarly divided into two main clades: Cobitis-Misgurnus-other genera (clade I) and Pangio-Lepidocephalichthys-other genera (clade II) (Fig. 8 and Supplementary Figure S2). Clade I included all analyzed species of Cobitis, Paramisgurnus and Misgurnus, and five species from other genus (I. longicorpa, K. multifasciata, N. delicata, K. naktongensis, and Microcobitis sp.). Four Pangio species, five Lepidocephalichthys species and other five species (K. pristes, A. choirorhynchos, A. gracilentus, L. macrochir, and C. gongota) were clustered into Clade II, among which the analyzed species of genus Pangio and Lepidocephalichthys formed two well-supported (pp = 1.00) monophyletic groups respectively. In addition, Pangio is the sister genus to Lepidocephalichthys.
The BI phylogenetic tree confirmed that Cobitis was a paraphyletic group, since Misgurnus clade A, N. delicate, I. longicorpa, and K. multifasciata shared the common ancestor with the all 15 Cobitis species analyzed in this study, with high posterior probability values (pp = 1.00). The species of Misgurnus were separated into two independent lineages: the majority of M. anguillicaudatus individuals (12/14) and M. bipartitus clustering with the Cobitis species (Misgurnus clade A), and two M. anguillicaudatus individuals, M. mizolepis, M. mohoity, and M. nikolskyi gathering with P. dabryanus and K. naktongensis (Misgurnus clade B).
Divergence time estimation of Cobitinae fishes
The combination of strict clock model and Yule process tree prior provided the best fit to the data sets (Supplementary Table 4). The chronogram with divergence time of Cobitinae lineages was estimated based on the cytB mutation rate (0.68% per million years) (Fig. 9). The first split of Cobitinae lineages was estimated to have occurred in the late Eocene (42.11 Ma, 95% HPD: 36.35–47.86 Ma), then separated into clade I (northern clade) and clade II (southern lineages). Cobitis-Iksookimia-Kichulchoia-Niwaella lineage diverged from the rest of northern clade lineage during the Oligocene (30.07 Ma, 95% HPD: 25.55–34.69 Ma), similar to the previous described , then diversified and further radiated after 4.94 Ma. The mtDNA introgression between ancestral species of Cobitis and ancestral species of Misgurnus seems to have taken place in the Middle Miocene (14.40 Ma, 95% HPD: 12.30–16.54 Ma). C. macrostigma appeared about 0.36 Ma (95% HPD: 0.06–0.55 Ma) in the Pleistocene. Pangio-Lepidocephalichthys-other genera (southern lineages) might originate about 40.45 Ma. In southern lineages, Pangio was estimated to have occurred about 20.14–29.88 Ma, and the divergence times of the four species analyzed in this study are congruent with the previous described dating .
In this study, we conducted a comparative mitogenome analysis and revealed the conserved and unique characteristics of 58 Cobitinae mitogenomes. Cobitinae mitogenomes display highly conserved tRNA secondary structure, overlaps and non-coding intergenic spacers. Among the 22 tRNAs, tRNAser(AGN) (S1) is the only one that is not folded into the typical clover-leaf secondary structure (Fig. 2a). Loss of stem in S1 is common character among Cobitinae and other metazoan mitogenomes [57, 58]. Similarly, the widespread unmatched base pairs among Cobitinae tRNAs is also a conserved feature in the eukaryote mitogenome [59,60,61]. Although their functions are not clear in fish, the unmatched base pairs are considered as the current state of evolutionary and irreversible process, which might be caused by tRNA editing .
Like other cyprinid fishes [14, 63], two long overlaps and two long IGSs were found in Cobitinae mitogenomes. The motif “ATGCTAA” in nd4l-nd4 was conserved in vertebrates, including fish, turtle and human [14, 63,64,65,66]. However, in comparison with the conserved motif (ATGATAA) in other Cypriniformes fishes, there is a specific 3 bp insertion (GCA) in the atp8-atp6 overlap motif of Cobitinae and other loaches [67,68,69], indicating this insertion is a characteristic feature of loaches. IGSs are important for transcription and associated with gene rearrangement in insects [70,71,72]. It is commonly assumed that IGS had a rapid nucleotide substitution rate under relaxed selection . Moreover, Cobitinae mitogenomes share highly conserved sequences in IGSs that are immediately adjacent to tRNAs, such as “CTTTCCCGCC”, “AAGGCGGGA” and “AGC”. Whether these conserved sequences have a function or not and how they act awaits further investigation. As the longest IGSs, CR plays an important role in controlling the transcription and replication of mtDNA molecule by several domains and motifs [74, 75]. Although significant length variation were found in CR of vertebrate , the three domains can also be recognized in Cobitinae mitogenomes. Furthermore, the AT-skew and GC-skew of CR might reflect the strand asymmetry [77,78,79]. In teleost, the skew inversion of CR was only found in the mitogenomes of Albula glossodonta and Bathygadus antrode, showing a reversed strand asymmetry . The normal Cobitinae mitogenomes CR skewness indicates that the strand asymmetry is not reversed.
The phylogenetic analyses show the monophyly of the genus Pangio and Lepidocephalichthys, consistent with the previous study . However, Cobitis, the biggest genus of Cobitinae , is a complex and controversial paraphyletic group. Similar to the trees constructed by cyt b [25, 80, 81], Iksookimia, Kichulchoia and Niwaella species were nested within Cobitis, implying a close relationship among them. Perdices  proposed that these species of Iksookimia, Kichulchoia, and Niwaella might belong to genus Cobitis, as morphologically specialized species derived from a local Cobitis species. However, this assumption awaits more morphological, karyological and molecular investigation. In addition, our phylogenetic analysis confirmed the assumption that M. mizolepis and P. dabryanus are conspecific [33, 80] and the different lineages under the species name C. striata and C. takatsuensis might actually represent different species.
The species of Misgurnus were separated into two independent clade and clustered into Cobitis species and P. dabryanus-K. naktongensis, respectively. The same results were observed in the trees based on the cyt b  and 13 PCGs from 28 cobitidae species . However, all Misgurnus and Koreocobitis species were grouped into a monophyletic clade when their phylogenetic relationships were constructed by nuclear gene rag-1 . This incongruity between mitochondrial and nuclear gene trees was explained by the different evolutionary rate of markers, hybridization or introgression . It is commonly believed that hybridization and subsequent mtDNA introgression might occur between ancestral species of Cobitis and ancestral species of Misgurnus [35, 80]. In this study, we collected 14 mitogenomes from M. anguillicaudatus, which were divided into two genetically divergent clades. The similar phenomenon has been reported by several previous studies, which is explained by hybridization and mtDNA introgression [34, 35, 47, 83, 84]. Considering that M. anguillicaudatus clustered into the clade of Misgurnus and Koreocobitis by nuclear analyses , we supposed that the 12 mitogenomes (No. 1–12) of M. anguillicaudatus in Misgurnus clade A could be considered as the introgressed mtDNA type because of their close relationship with Cobitis species, whereas the other two individuals in Misgurnus clade B retained the original M. anguillicaudatus mitogenomes. M. anguillicaudatus with introgressed mtDNA type spread over most of East Asia, including China, Japan and Korea. M. anguillicaudatus shows extensive ploidy variability in nature. Besides most common diploid individuals (2n = 50), triploid (3n = 75) and tetraploid (4n =100) have been frequently recorded in some localities of China and Japan [21, 47, 85, 86]. Rare pentaploid (5n = 125) and even hexaploid (6n = 150) individuals were found in the Yangtze River basin . All of M. anguillicaudatus polyploids analyzed in this study belonged to the introgressed mtDNA type. Since mtDNA is inherited maternally, these polyploids might have originated from the diploid M. anguillicaudatus with introgressed mtDNA. Further analyses are needed to confirm this hypothesis of inter-genus mtDNA introgression based on a large-scale sampling with quantitative morphological features, definite ploidy, and more genes from both mitochondria and nuclear genomes.
The first split of Cobitinae lineages was estimated to have occurred in the late Eocene (42.11 Ma, 95% HPD: 36.35–47.86 Ma), separating northern clade and southern lineages, consistent with reconstruction dates of the paleo-drainages of East Asia [35, 88]. Cobitinae fishes in Clade I and Clade II, nominated as “northern clade” and “southern lineages” respectively, show a distinct disjunctive distribution with a small area of sympatry in Vietnam . Consistent with their locations, the northern clade spread to most of East Asia, Siberia and Europe, while the southern lineages distribute across the Indian subcontinent and Southeast Asia after their isolation. The nodes within northern clade and southern lineage appear asynchronous, implying that some local dominant factors, rather than large-scale events, might shape the evolution within northern or southern lineage.
This study represents the first comparative mitogenome and phylogenetic analyses within Cobitinae. The conserved and unique characteristics of 58 Cobitinae mitogenomes were revealed. We observed distinct base compositions among different genus and identified a specific 3 bp insertion (GCA) in the atp8-atp6 overlap as a unique feature of loaches. ML and BI analyses both strongly support the paraphyly of Cobitis and polyphyly of Misgurnus. In addtion, Cobitinae might have split into northern and southern lineages in the late Eocene (42.11 Ma), and a mtDNA introgression between Cobitis and Misgurnus might have occured about 14.40 Ma. The current study provides new insights into the mitogenome features and evolution of Cobitinae fishes.
Sampling, sequencing and assembly
The C. macrostigma analyzed in this study was caught from the Yangtze River in Yibin City, Sichuan Province, China (N: 28°46′6.01″, E: 104°38′13.99″) in October 2018 and five individual were transported to the laboratory (National Aquatic Biological Resource Center, NABRC) in oxygen-rich water. It possesses 5–9 large and round spot in the midline of lateral body side  (Fig. 1). Before sampling, they were reared in a square and glass recirculating freshwater tanks with a volume of about 100 L, at 22 °C on a 14 h (hour) light/10 h dark cycle for morphological identification. After deep and overdosed anesthesia with styrylpyridine (a common anaesthetic used in fish, 30-50 mg/L; aladdin, China), one healthy one-year-old female fish, 7 cm in length and 1.8 g in weight, was euthanized by immediately cutting off the spinal cord adjacent to the head. Total DNA was extracted according to the Ezup Column Animal Genomic DNA Kit technical manual (Sangon, Shanghai, China). PCR primers were designed based on the conserved sequences between the mitogenomes of C. granoei (GenBank: NC_023473.1) and C. sinensis (GenBank: NC_007229.1). 742–2495 bp DNA were amplified by using High Fidelity DNA Polymerase (Yeasen, Shanghai) (Supplementary Table 1). To obtain accurate sequences, we chose a cloning strategy. According to manual, PCR amplicon was purified, ligated ESI-Blunt vector (Yeasen, Shanghai) and transfected into 5α Chemically Competent Cell (Tsingke Biological Technology, Beijing). The positive clones were sequenced by Quintara Biosciences (Wuhan, China). The segments, longer than 1500 bp, were sequenced using the primer walking sequencing strategy. The resulting DNA sequences were assembled using DNAStar (DNASTAR Inc., USA) . Other 57 Cobitinae mitogenomes were download from NCBI GenBank database [38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55].
Gene annotation and bioinformatic analyses
tRNA genes and their secondary structures were predicted with MITOS  and tRNAscan-SE 2.0 with default parameters . All 13 PCGs and two rRNA genes were annotated by comparison with the sequences of other Cobitinae fishes in GenBank (https://blast.ncbi.nlm.nih.gov/). The mtDNA maps were drawn using CGView Server V1.0 . The sequence logos of gene overlaps and non-coding IGSs were drawn using WebLogo 3.7.4 . The base composition, codon distributions and relative synonymous codons usage were calculated using DNAStar (DNASTAR Inc., USA) , MEGA 7.0  and Microsoft Excel 2010. Skewness was measured using the formulas: AT-skew = (A% - T%) / (A% + T%) and GC-skew = (G% - C%) / (G% + C%) . The silimlarity of the sequences was calculated in MEGA 7.0  under p-distance and NCBI-BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi).
The phylogenetic analysis was performed based on 13 PCGs of 58 Cobitinae mitogenomes. Sinorhodeus microlepis and Rhodeus shitaiensis were chosen as the outgroups (Table 1). Each of the 13 gene sequences was separately aligned using Muscle v3.8.31  and concatenated into a sequence matrix by PhyloSuite v1.2.2 . Then PartitionFinder2  was used to find the best partitioning strategy and to calculate the best-fit evolutionary models for each subset. For the alignment, a scheme with eight partitions was selected and GTR + G + I was chosen as the best-fit evolutionary model for each partition. Phylogenetic trees were constructed by the maximum likelihood (ML) method and bayesian inference (BI). The ML method was implemented in RAxML v8.2.12 . Each partition scheme was run with the GTRGAMMAI model, and 1000 rapid bootstrapping replications were set to evaluate the bootstrap support values and search for the best-scoring ML tree. The BI phylogeny was performed in MrBayes v3.1.2  with the “unlink” and “prest ratepr = variable” model parameters. 10,000,000 generations were run in two independent runs of four independent Markov Chain Monte Carlo (MCMC) chains, and were sampled every 1000 generations. The convergence of the BI analyses was investigated using Tracer v1.7.1 software. The first 2500 trees were discarded as conservative burn-in, and the rests were used to generate a majority rule consensus tree.
In cobitid fishes, 0.680% (divergence per pairwise comparison per Ma) was calculated and suggested for the mutation rates of cytb gene . In this study, BEAST v1.10.4  was used to estimate the divergence time with the rate (0.68%). GTR + G + I was chosen as the best fit model by PartitionFinder2 . The best-fit clock type and tree prior were selected from two clock models (strict clock and uncorrelated relaxed clock) and four tree priors (Yule process, Exponential growth, Constant size and Bayesian skyline) by comparing the marginal likelihood values estimated by path sampling . The analyses were simultaneously run for 20,000,000 generations, with parameters sampled every 1000, then the first 25% of the trees were discarded as burn-in. Tracer v1.5  and Figtree were used to assess the convergence and view trees, respectively.
Availability of data and materials
C. macrostigma mitochondrial genome has been deposited in GenBank under the accession numbers MT259034. The 59 mitogenomes from Cobitinae species, Sinorhodeus microlepis and Rhodeus shitaiensis were downloaded from GenBank. Their accession numbers and references were listed in Table 1. Other supporting results are included within the article and its additional files.
Million years ago
Protein coding gene
- OL :
Origin of L-strand replication
- cytb :
- rag-1 :
Recombination activating gene 1
- nd1–6 :
NADH dehydrogenase subunit 1–6
- nd4l :
NADH dehydrogenase subunit 4 L
- cox1–3 :
Cytochrome oxidase subunit I-III
- atp6 :
ATPase subunit 6
- atp8 :
ATPase subunit 8
Non-coding intergenic spacer
Relative synonymous codon usage
- ω or dN/dS:
Non-synonymous and synonymous substitutions
Highest posterior density
Satoh TP, Miya M, Mabuchi K, Nishida M. Structure and variation of the mitochondrial genome of fishes. BMC Genomics. 2016;17:719.
Cuppari A, Fernandez-Millan P, Battistini F, Tarres-Sole A, Lyonnais S, Iruela G, et al. DNA specificities modulate the binding of human transcription factor a to mitochondrial DNA control region. Nucleic Acids Res. 2019;47(12):6519–37.
Shadel GS, Clayton DA. Mitochondrial DNA maintenance in vertebrates. Annu Rev Biochem. 1997;66:409–35.
Lin CP, Danforth BN. How do insect nuclear and mitochondrial gene substitution patterns differ? Insights from Bayesian analyses of combined datasets. Mol Phylogenet Evol. 2004;30(3):686–702.
Elson JL, Lightowlers RN. Mitochondrial DNA clonality in the dock: can surveillance swing the case? Trends Genet. 2006;22(11):603–7.
Curole JP, Kocher TD. Mitogenomics: digging deeper with complete mitochondrial genomes. Trends Ecol Evol. 1999;14(10):394–8.
Cole TL, Ksepka DT, Mitchell KJ, Tennyson AJD, Thomas DB, Pan HL, et al. Mitogenomes uncover extinct penguin taxa and Reveal Island formation as a key driver of speciation. Mol Biol Evol. 2019;36(4):784–97.
Baron S, van der Merwe NA, Maritz-Olivier C. The genetic relationship between R. microplus and R. decoloratus ticks in South Africa and their population structure. Mol Phylogenet Evol. 2018;129:60–9.
Min-Shan Ko A, Zhang Y, Yang MA, Hu Y, Cao P, Feng X, et al. Mitochondrial genome of a 22,000-year-old giant panda from southern China reveals a new panda lineage. Curr Biol. 2018;28(12):R693–4.
Wei M, Yu P, Yang Y, Wan Q. The complete mitochondrial genome of Leptobotia taeniaps (Cypriniformes: Cobitidae). Mitochondrial DNA part a, DNA mapping, sequencing, and analysis. Mitochondrial DNA A DNA Mapp Seq Anal. 2016;27(3):1707–8.
Luo H, Kong X, Chen S, Shi W. Mechanisms of gene rearrangement in 13 bothids based on comparison with a newly completed mitogenome of the threespot flounder, Grammatobothus polyophthalmus (Pleuronectiformes: Bothidae). BMC Genomics. 2019;20(1):792.
Dong S, Zhao C, Zhang S, Zhang L, Wu H, Liu H, et al. Mitochondrial genomes of the early land plant lineage liverworts (Marchantiophyta): conserved genome structure, and ongoing low frequency recombination. BMC Genomics. 2019;20(1):953.
Ponts N, Gautier C, Gouzy J, Pinson-Gadais L, Foulongne-Oriol M, Ducos C, et al. Evolution of Fusarium tricinctum and Fusarium avenaceum mitochondrial genomes is driven by mobility of introns and of a new type of palindromic microsatellite repeats. BMC Genomics. 2020;21(1):358.
Yu P, Zhou L, Zhou XY, Yang WT, Zhang J, Zhang XJ, et al. Unusual AT-skew of Sinorhodeus microlepis mitogenome provides new insights into mitogenome features and phylogenetic implications of bitterling fishes. Int J Biol Macromol. 2019;129:339–50.
Froese R, Pauly D. FishBase. World wide web electronic publication. Wwwfishbaseorg; 2019. version 04.
Gao J, Koshio S, Nguyen BT, Wang WM, Cao XJ. Comparative studies on lipid profiles and amino acid composition of wild and cultured dojo loach Misgurnus anguillicaudatus obtained from southern Japan. Fish Sci. 2012;78(6):1331–6.
Chen J, Chen L. Effects of chitosan-supplemented diets on the growth performance, nonspecific immunity and health of loach fish (Misgurnus anguillicadatus). Carbohydr Polym. 2019;225:115227.
Gao Y, He J, He Z, Li Z, Zhao B, Mu Y, et al. Effects of fulvic acid on growth performance and intestinal health of juvenile loach Paramisgurnus dabryanus (Sauvage). Fish Shellfish Immunol. 2017;62:47–56.
Juchno D, Jablonska O, Boron A, Kujawa R, Leska A, Grabowska A, et al. Ploidy-dependent survival of progeny arising from crosses between natural allotriploid Cobitis females and diploid C-taenia males (Pisces, Cobitidae). Genetica. 2014;142(4):351–9.
Morishima K, Yoshikawa H, Arai K. Diploid clone produces unreduced diploid gametes but tetraploid clone generates reduced diploid gametes in the Misgurnus loach. Biol Reprod. 2012;86(2):33.
Zhou L, Gui J. Natural and artificial polyploids in aquaculture. Aquaculture Fish. 2017;2(3):103–11.
Janko K, Culling MA, Rab P, Kotlik P. Ice age cloning--comparison of the quaternary evolutionary histories of sexual and clonal forms of spiny loaches (Cobitis; Teleostei) using the analysis of mitochondrial DNA variation. Mol Ecol. 2005;14(10):2991–3004.
Janko K, Bohlen J, Lamatsch D, Flajshans M, Epplen JT, Rab P, et al. The gynogenetic reproduction of diploid and triploid hybrid spined loaches (Cobitis: Teleostei), and their ability to establish successful clonal lineages--on the evolution of polyploidy in asexual vertebrates. Genetica. 2007;131(2):185–94.
Bohlen J, Šlechtová V, Tan HH, Britz R. Phylogeny of the southeast Asian freshwater fish genus Pangio (Cypriniformes; Cobitidae). Mol Phylogenet Evol. 2011;61(3):854–65.
Chen YX, Chen YF. Three new species of cobitid fish (Teleostei, Cobitidae) from the river Xinjiang and the river Le'anjiang, tributaries of Lake Poyang of China, with remarks on their classification. Folia Zool. 2013;62(2):83–95.
Peng G, Zhu B, Yang D, Su L, Shi H, Li D. Microplastics in sediments of the Changjiang estuary, China. Environ Pollut. 2017;225:283–90.
Wang J, Wang M, Ru S, Liu X. High levels of microplastic pollution in the sediments and benthic organisms of the South Yellow Sea, China. Sci Total Environ. 2019;651(Pt 2):1661–9.
Buj I, Caleta M, Marcic Z, Sanda R, Vukic J, Mrakovcic M. Different histories, different DestiniesImpact of evolutionary history and population genetic structure on extinction risk of the Adriatic Spined loaches (genus Cobitis; Cypriniformes, Actinopterygii). PLoS One. 2015;10(7):e0131580.
Zhao S, Feng C, Quan W, Chen X, Niu J, Shen Z. Role of living environments in the accumulation characteristics of heavy metals in fishes and crabs in the Yangtze River estuary, China. Mar Pollut Bull. 2012;64(6):1163–71.
Allert AL, Fairchild JF, Schmitt CJ, Besser JM, Brumbaugh WG, Olson SJ. Effects of mining-derived metals on riffle-dwelling benthic fishes in Southeast Missouri, USA. Ecotoxicol Environ Saf. 2009;72(6):1642–51.
Perdices A, Bohlen J, Doadrio I. The molecular diversity of adriatic spined loaches (Teleostei, Cobitidae). Mol Phylogenet Evol. 2008;46(1):382–90.
Doadrio I, Perdices A. Phylogenetic relationships among the Ibero-African cobitids (Cobitis, cobitidae) based on cytochrome b sequence data. Mol Phylogenet Evol. 2005;37(2):484–93.
Slechtova V, Bohlen J, Tan HH. Families of Cobitoidea (Teleostei; Cypriniformes) as revealed from nuclear genetic data and the position of the mysterious genera Barbucca, Psilorhynchus, Serpenticobitis and Vaillantella. Mol Phylogenet Evol. 2007;44(3):1358–65.
Morishima K, Nakamura-Shiokawa Y, Bando E, Li YJ, Boron A, Khan MM, et al. Cryptic clonal lineages and genetic diversity in the loach Misgurnus anguillicaudatus (Teleostei: Cobitidae) inferred from nuclear and mitochondrial DNA analyses. Genetica. 2008;132(2):159–71.
Slechtova V, Bohlen J, Perdices A. Molecular phylogeny of the freshwater fish family Cobitidae (Cypriniformes: Teleostei): delimitation of genera, mitochondrial introgression and evolution of sexual dimorphism. Mol Phylogenet Evol. 2008;47(2):812–31.
Tang QY, Shi LX, Liu F, Yu D, Liu HZ. Evolution and phylogenetic application of the MC1R gene in the Cobitoidea (Teleostei: Cypriniformes). Zool Res. 2016;37(5):281–9.
Perdices A, Bohlen J, Slechtova V, Doadrio I. Molecular evidence for multiple origins of the European Spined loaches (Teleostei, Cobitidae). PLoS One. 2016;11(1):e0144628.
Saitoh K, Sado T, Mayden RL, Hanzawa N, Nakamura K, Nishida M, et al. Mitogenomic evolution and interrelationships of the Cypriniformes (Actinopterygii: Ostariophysi): the first evidence toward resolution of higher-level relationships of the world's largest freshwater fish clade based on 59 whole mitogenome sequences. J Mol Evol. 2006;63(6):826–41.
Miya M, Sato Y, Fukunaga T, Sado T, Poulsen JY, Sato K, et al. MiFish, a set of universal PCR primers for metabarcoding environmental DNA from fishes: detection of more than 230 subtropical marine species. R Soc Open Sci. 2015;2(7):150088.
Kim KY, Lee SY, Bang IC, Nam YK. Complete mitogenome sequence of an endangered freshwater fish, Iksookimia choii (Teleostei; Cypriniformes; Cobitidae). Mitochondrial DNA. 2008;19(5):438–45.
Huang S, Tomljanovic T, Tian X, Wang Y, Cao X. The complete mitochondrial genome of natural Cobitis elongatoides (Cypriniformes: Cobitidae). Mitochondrial DNA part a, DNA mapping, sequencing, and analysis. Mitochondrial DNA A DNA Mapp Seq Anal. 2016;27(1):189–90.
Song JR, You P. The complete mitochondrial genome of Cobitis granoei (Cyprinformes: Cobitidae). Mitochondrial DNA part a, DNA mapping, sequencing, and analysis. Mitochondrial DNA A DNA Mapp Seq Anal. 2016;27(1):60–1.
Park HK, Kim KS, Kim KY, Bang IC. The full-length mitochondrial genome of Cobitis nalbanti (Teleostei: Cypriniformes: Cobitidae). Mitochondrial DNA Part B-Resources. 2018;3(2):872–3.
Saitoh K, Chen WJ, Mayden RL. Extensive hybridization and tetrapolyploidy in spined loach fish. Mol Phylogenet Evol. 2010;56(3):1001–10.
Saitoh K, Miya M, Inoue JG, Ishiguro NB, Nishida M. Mitochondrial genomics of ostariophysan fishes: perspectives on phylogeny and biogeography. J Mol Evol. 2003;56(4):464–72.
Yu YY, Li YH, Li RW, Wang WM, Zhou XY. Mitochondrial genome of the natural tetraploid loach Misgurnus anguillicaudatus. Mitochondrial DNA. 2014;25(2):115–6.
Zhou X, Yu Y, Li Y, Wu J, Zhang X, Guo X, et al. Comparative analysis of mitochondrial genomes in distinct nuclear ploidy loach Misgurnus anguillicaudatus and its implications for polyploidy evolution. PLoS One. 2014;9(3):e92033.
Yu YY, Li YH, Li RW, Wang WM, Zhou XY. Complete mitochondrial genome of the natural hexaploid loach, Misgurnus anguillicaudatus (Teleostei: Cypriniformes: Cobitididae). Mitochondrial DNA. 2014;25(2):100–1.
Zhang X, Wang W, Huang S, Chen G, Bai X, Cao X. The complete mitochondrial genomes of natural diploid and tetraploid loaches Misgurnus anguillicaudatus (Cypriniformes: Cobitidae). Mitochondrial DNA. 2014;25(3):196–7.
Zeng L, Wang J, Sheng J, Gu Q, Hong Y. Molecular characteristics of mitochondrial DNA and phylogenetic analysis of the loach (Misgurnus anguillicaudatus) from the Poyang Lake. Mitochondrial DNA. 2012;23(3):187–200.
He S, Gu X, Mayden RL, Chen WJ, Conway KW, Chen Y. Phylogenetic position of the enigmatic genus Psilorhynchus (Ostariophysi: Cypriniformes): evidence from the mitochondrial genome. Mol Phylogenet Evol. 2008;47(1):419–25.
Huang S, Tian X, Wang W, Song W, Zhang X, Bai X, et al. The complete mitochondrial genome of natural Misgurnus bipartitus (Cypriniformes: Cobitidae). Mitochondrial DNA. 2015;26(5):680–1.
Yu YY, Song W, Wang YZ, Wang WM, Zhou XY. Complete mitochondrial genome of the Amur weatherfish, Misgurnus mohoity (Teleostei: Cypriniformes: Cobitididae). Mitochondrial DNA. 2015;26(2):310–2.
Zhu M, Liu F. Mitochondrial genome sequence of Paramisgurnus dabryanus from the Yellow River estuary: implication for Cobitidae phylogeny. Mitochondrial DNA part a, DNA mapping, sequencing, and analysis. Mitochondrial DNA A DNA Mapp Seq Anal. 2016;27(4):3039–40.
Li F, Shao KT, Lin YS, Chang CH. The complete mitochondrial genome of the Rhodeus shitaiensis (Teleostei, Cypriniformes, Acheilognathidae). Mitochondrial DNA. 2015;26(2):301–2.
Castellana S, Vicario S, Saccone C. Evolutionary patterns of the mitochondrial genome in Metazoa: exploring the role of mutation and selection in mitochondrial protein coding genes. Genome Biol Evol. 2011;3:1067–79.
Zhao Q, Wang J, Wang MQ, Cai B, Zhang HF, Wei JF. Complete mitochondrial genome of dinorhynchus dybowskyi (hemiptera: pentatomidae: asopinae) and phylogenetic analysis of pentatomomorpha species. J Insect Sci. 2018;18(2):1.
Lavrov DV, Brown WM. Trichinella spiralis mtDNA: a nematode mitochondrial genome that encodes a putative ATP8 and normally structured tRNAS and has a gene arrangement relatable to those of coelomate metazoans. Genetics. 2001;157(2):621–37.
Smith BA, Jackman JE. Saccharomyces cerevisiae Thg1 uses 5’-pyrophosphate removal to control addition of nucleotides to tRNA (his). Biochemistry-Us. 2014;53(8):1380–91.
Betat H, Long Y, Jackman JE, Morl M. From end to end: tRNA editing at 5’- and 3’-terminal positions. Int J Mol Sci. 2014;15(12):23975–98.
Gray MW. Evolutionary origin of RNA editing. Biochemistry-Us. 2012;51(26):5235–42.
Lavrov DV, Brown WM, Boore JL. A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus. Proc Natl Acad Sci U S A. 2000;97(25):13738–42.
Broughton RE, Milam JE, Roe BA. The complete sequence of the zebrafish (Danio rerio) mitochondrial genome and evolutionary patterns in vertebrate mitochondrial DNA. Genome Res. 2001;11(11):1958–67.
Yu P, Ding S, Yang Q, Bi Z, Chen L, Liu X, et al. Complete sequence and characterization of the paradise fish Macropodus erythropterus (Perciformes: Macropodusinae) mitochondrial genome. Mitochondrial DNA Part B. 2016;1(1):54–5.
Yu P, Yang X, Zhou W, Yang W, Zhou L, Liu X, et al. Comparative mitogenomic and phylogenetic analysis of Apalone spinifera and Apalone ferox (Testudines: Trionychidae). Genetica. 2019;147(2):165–76.
Ermini L, Olivieri C, Rizzi E, Corti G, Bonnal R, Soares P, et al. Complete mitochondrial genome sequence of the Tyrolean iceman. Curr Biol. 2008;18(21):1687–93.
Yu P, Ding S, Yang Q, Li X, Wan Q. The complete mitochondrial genome of Sinibotia robusta (Cypriniformes: Cobitidae). Mitochondrial DNA part a, DNA mapping, sequencing, and analysis. Mitochondrial DNA A DNA Mapp Seq Anal. 2016;27(5):3471–2.
Kanu UC, Zhao G, Xie P, Li Y, Lei D, Niu J, et al. The complete mtDNA genome of Triplophysa strauchii (Cypriniformes, Balitoridae, Cobitoidea): genome charaterization and phylogenetic analysis. Mitochondrial DNA part a, DNA mapping, sequencing, and analysis. Mitochondrial DNA A DNA Mapp Seq Anal. 2016;27(4):2637–8.
Wu J, He Y, Ren H, Zhang Y, Du Z, Xie M, et al. The complete mitochondrial genome sequence of Beaufortia szechuanensis (Cypriniformes, Balitoridae). Mitochondrial DNA part a, DNA mapping, sequencing, and analysis. Mitochondrial DNA A DNA Mapp Seq Anal. 2016;27(4):2535–6.
Taanman JW. The mitochondrial genome: structure, transcription, translation and replication. Biochim Biophys Acta. 1999;1410(2):103–23.
Mao M, Valerio A, Austin AD, Dowton M, Johnson NF. The first mitochondrial genome for the wasp superfamily Platygastroidea: the egg parasitoid Trissolcus basalis. Genom Natl Res Council Can. 2012;55(3):194–204.
Rodovalho Cde M, Lyra ML, Ferro M, Bacci M Jr. The mitochondrial genome of the leaf-cutter ant Atta laevigata: a mitogenome with a large number of intergenic spacers. PLoS One. 2014;9(5):e97117.
Ghikas DV, Kouvelis VN, Typas MA. Phylogenetic and biogeographic implications inferred by mitochondrial intergenic region analyses and ITS1–5.8S-ITS2 of the entomopathogenic fungi Beauveria bassiana and B. brongniartii. BMC Microbiol. 2010;10:174.
Zhang DX, Hewitt GM. Insect mitochondrial control region: a review of its structure, evolution and usefulness in evolutionary studies. Biochem Syst Ecol. 1997;25(2):99–120.
Fonseca MM, Posada D, Harris DJ. Inverted replication of vertebrate mitochondria. Mol Biol Evol. 2008;25(5):805–8.
Rand DM. Endotherms, ectotherms, and mitochondrial genome-size variation. J Mol Evol. 1993;37(3):281–95.
Wei SJ, Shi M, Chen XX, Sharkey MJ, van Achterberg C, Ye GY, et al. New views on strand asymmetry in insect mitochondrial genomes. PLoS One. 2010;5(9):e12708.
Wei SJ, Shi M, Sharkey MJ, van Achterberg C, Chen XX. Comparative mitogenomics of Braconidae (Insecta: hymenoptera) and the phylogenetic utility of mitochondrial genomes with special reference to Holometabolous insects. BMC Genomics. 2010;11:371.
Perna NT, Kocher TD. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995;41(3):353–8.
Perdices A, Ozeren CS, Erkakan F, Freyhof J. Diversity of spined loaches from asia minor in a phylogenetic context (teleostei: cobitidae). PLoS One. 2018;13(10):1.
Tan XC, Li P, Wu TJ, Yang J. Cobitis xui, a new species of spined loach (Teleostei: Cobitidae) from the Pearl River drainage in southern China. Zootaxa. 2019;4604(1):4601.
Saitoh K, Kim IS, Lee EH. Mitochondrial gene introgression between spined loaches via hybridogenesis. Zool Sci. 2004;21(7):795–8.
Koizumi N, Takemura T, Watabe K, Mori A. Genetic variation and diversity of Japanese loach inferred from mitochondrial DNA. Transact Japanese Soc Irrigation. 2010;77:7–16.
Kitagawa T, Fujii Y, Koizumi N. Origin of the two major distinct mtDNA clades of the Japanese population of the oriental weather loach Misgurnus anguillicaudatus (Teleostei: Cobitidae). Folia Zool. 2011;60(4):343–9.
Arai K. Genetics of the loach, Misgurnus anguillicaudatus: recent progress and perspective. Folia Biol (Krakow). 2003;51:107–17.
XJ Y, T Z, YC L, K L, M Z. Chromosomes of Chinese freshwater fishes. Beijing: Science Press; 1989.
Abbas K, Li MY, Wang WM, Zhou XY. First record of the natural occurrence of hexaploid loach Misgurnus anguillicaudatus in Hubei Province, China. J Fish Biol. 2009;75(2):435–41.
Clark MK, Schoenbohm LM, Royden LH, Whipple KX, Burchfiel BC, Zhang X, et al. Surface uplift, tectonics, and erosion of eastern Tibet from large-scale drainage patterns. Tectonics. 2004;23(1):1.
Chen Q, Zheng B. Systematic synopsis of chinese fishes, vol. 1. Beijing, China: Reaktion books; 1987.
Burland TG. DNASTAR’s Lasergene sequence analysis software. Methods Mol Biol. 2000;132:71–91.
Bernt M, Donath A, Juhling F, Externbrink F, Florentz C, Fritzsch G, et al. MITOS: improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 2013;69(2):313–9.
Lowe TM, Chan PP. tRNAscan-SE on-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016;44(W1):W54–7.
Grant JR, Stothard P. The CGView server: a comparative genomics tool for circular genomes. Nucleic Acids Res. 2008;36:W181–4.
Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14(6):1188–90.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Zhang D, Gao F, Jakovlic I, Zou H, Zhang J, Li WX, et al. PhyloSuite: an integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol Ecol Resour. 2020;20(1):348–55.
Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587–9.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.
Pang S, Stones RJ, Ren MM, Liu XG, Wang G, Xia HJ, et al. V3.1: MrBayes on graphics processing units for protein sequence data. Mol Biol Evol. 2015;32(9):2496–7.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.
Baele G, Lemey P, Bedford T, Rambaut A, Suchard MA, Alekseyenko AV. Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. Mol Biol Evol. 2012;29(9):2157–67.
Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior summarization in Bayesian Phylogenetics using tracer 1.7. Syst Biol. 2018;67(5):901–4.
The research was supported by the Wuhan Branch, Supercomputing Centre, Chinese Academy of Sciences, China.
This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB31000000) and the China Agriculture Research System (CARS-45-07).
Ethics approval and consent to participate
The C. macrostigma analyzed in this study was caught from the Yangtze River in Yibin City, Sichuan Province, China and reared in the National Aquatic Biological Resource Center (NABRC). The acquisition of experimental fish complies with the laws of Fishery Administration of the Ministry of Agriculture and Rural Affairs of the People’s Republic of China. We confirm that C. macrostigma is not an endangered or protected species (http://www.iucnredlist.org). The protocol, including the research.
question, key design features, and analysis plan, was provided to the Animal Care and Use Committee of Institute of Hydrobiology, Chinese Academy of Sciences before the study, and all procedures in this research performed with the approval of the Committee. No ethics approval was required for the public sequence data used in this study.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of primers used to amplify the mitogenome of C. macrostigma.
Length, base composition and skewness of Cobitinae fish mitogenomes.
Start and stop codons of 13 PCGs in Cobitinae mitogenomes.
Marginal likelihood values of different combinations of clock model and tree prior.
Codon distribution (A) and relative synonymous codon usage (B) of PCGs in the 58 Cobitinae mitogenomes. CDpT = codons per thousand codons.
Phylogenetic tree constructed by ML methods, based on 13 PCGs of 58 Cobitinae mitogenomes. Sinorhodeus microlepis and Rhodeus shitaiensis were chosen as outgroups. Node numbers represent the bootstrap value.
About this article
Cite this article
Yu, P., Zhou, L., Yang, WT. et al. Comparative mitogenome analyses uncover mitogenome features and phylogenetic implications of the subfamily Cobitinae. BMC Genomics 22, 50 (2021). https://doi.org/10.1186/s12864-020-07360-w
- Mitochondrial genome
- mtDNA introgression
- Divergence time