Molecular Phylogeny and Conservation Priorities of the Subfamily Acheilognathinae (Teleostei: Cyprinidae)

Novel insights into adaptive evolution based on the unusual AT-skew in Acheilognathus gracilis mitogenome and phylogenetic relationships of bitterling

Acheilognathus gracilis, a bitterling species, distribute in lower reaches of Yangtze River. They are identified as the top-priority bitterling species for conservation as having high evolutionary distinctiveness and are at risk of extinction. In present study, we first sequenced the complete mitogenome of A. gracilis and analyzed its phylogenetic position using 13 PCGs. The A. gracilis mitogenome is 16,774 bp in length, including 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, a control region and the origin of the light strand replication. The overall base composition of A. gracilis in descending order is T 27.9 %, A 27.7 %, C 26.1 % and G 18.3 %, shows a unusual AT-skew with slightly negative. Further investigation revealed A. gracilis uses excess T over A in NADH dehydrogenase 5 (nd5), whereas the most of other bitterlings are biased toward to use A not T, implying there is likely to be unique strategy of adaptive evolution in A. gracilis. We also compared 13 PCGs of 30 bitterling mitogenomes and the results exhibit highly conservative. Phylogenetic trees constructed by 13 PCGs strongly support the monophyly of Acheilognathus and the paraphyly of Rhodeus and Tanakia. Current results will provide valuable information for follow-up research on conservation of species facing with serious population decline and can provide novel insights into the phylogenetic analysis and evolutionary biology research.

Mitochondrial genome of Acheilognathusbarbatulus (Cypriniformes, Cyprinidae, Acheilognathinae): characterisation and phylogenetic analysis

Acheilognathusbarbatulus is distributed in Yangtze River, Yellow River and Pearl River systems in China. Genome data can help to understand the phylogenetic relationships of A.barbatulus, but its complete mitochondrial genome has not been published. We determined the complete mitochondrial genome structure and characteristics of this species and constructed a comprehensive phylogenetic tree, based on mitochondrial genome data of several species of Acheilognathus, Rhodeus and Pseudorasboraparva. The complete length of the mitochondrial genome of A.barbatulus is 16726 bp. The genome is a covalently closed double-stranded circular molecule containing 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, a D-loop and a light strand replication initiation region. The base composition of the complete mitochondrial genome is A (29.33%) > T (27.6%) > C (26.12%) > G (16.95%), showing a strong AT preference and anti-G bias. All 13 PCGs have different degrees of codon preference, except for cytochrome c oxidase 1, which uses GTG as the start codon. All the PCGs use ATG as the start codon and the stop codon is dominated by TAG. The encoded amino acids Leu and Ser exist in two types, whereas the rest are all present as one type, except for tRNASer (GCT), which lacks the D-arm and has an incomplete secondary structure, all other tRNAs can be folded to form a typical cloverleaf secondary structure. Based on the 13 PCG tandems, the Maximum Likelihood and Bayesian trees were constructed, based on the concatenated sequence of 13 PCGs for the genera Acheilognathus and Rhodeus, with Pseudorasboraparva as the outgroup. Acheilognathusbarbatulus, Acheilognathustonkinensis and Acheilognathuscf.macropterus were clustered together and the most closely related. The results of this study enrich the mitochondrial genomic data of Acheilognathus and provide molecular and genetic base information for species conservation, molecular identification and species evolution of Acheilognathinae.

Complete mitochondrial genome of Rhodeus cyanorostris (Teleostei, Cyprinidae): characterization and phylogenetic analysis

Rhodeuscyanorostris Li, Liao & Arai, 2020 is a freshwater fish that is endemic to China and restricted to Chengdu City in Sichuan Province. This study is the first to sequence and characterize the complete mitochondrial genome of R.cyanorostris. The mitogenome of R.cyanorostris is 16580 bp in length, including 13 protein-coding genes, two rRNA genes, 22 tRNA genes, and a control region (D-loop). The base composition of the sequence is 28.5% A, 27.6% C, 26.4% T, and 17.5% G, with a bias toward A+T. The genome structure, nucleotide composition, and codon usage of the mitogenome of R.cyanorostris are consistent with those of other species of Rhodeus. To verify the molecular phylogeny of the genus Rhodeus, we provide new insights to better understand the taxonomic status of R.cyanorostris. The phylogenetic trees present four major clades based on 19 mitogenomic sequences from 16 Rhodeus species. Rhodeuscyanorostris exhibits the closest phylogenetic relationship with R.pseudosericeus, R.amarus, and R.sericeus. This study discloses the complete mitochondrial genome sequence of R.cyanorostris for the first time and provides the most comprehensive phylogenetic reconstruction of the genus Rhodeus based on whole mitochondrial genome sequences. The information obtained in this study will provide new insights for conservation, phylogenetic analysis, and evolutionary biology research.

Individual copy number variation and extensive diversity between major MHC-DAB1 allelic lineages in the European bitterling

Polymorphism of the major histocompatibility complex (MHC), DAB1 gene was characterized for the first time in the European bitterling (Rhodeus amarus), a freshwater fish employed in studies of host-parasite coevolution and mate choice, taking advantage of newly designed primers coupled with high-throughput amplicon sequencing. Across 221 genotyped individuals, we detected 1–4 variants per fish, with 28% individuals possessing 3–4 variants. We identified 36 DAB1 variants, and they showed high sequence diversity mostly located within predicted antigen-binding sites, and both global and codon-specific excess of non-synonymous mutations. Despite deep divergence between two major allelic lineages, functional diversity was surprisingly low (3 supertypes). Overall, these findings suggest the role of positive and balancing selection in promotion and long-time maintenance of DAB1 polymorphism. Further investigations will clarify the role of pathogen-mediated selection to drive the evolution of DAB1 variation.

Comparative phylogenetic analysis of Schizothorax plagiostomus and Schizothorax esocinus with other members of subfamilies of cyprinidae on the basis of complete mitochondrial genome and 12S, 16S ribosomal RNA from Northren areas of Pakistan

We assessed the relationship of Schizothoracinae species with other subfamilies Alburninae, Xenocyprinae, Cultrinae and Squaliobarbinae of family Cyprinidae by creating the phylogenetic trees using complete mitogenome and 12S and 16S RNA. Our representative species show the great affiliation with other but separated from a group composed of Metzia mesembrinum, Metzia longinasus, Metzia lineata and Metzia formosae of subfamily Alburninae while other subfamilies formed distinct group. The members of subfamily Schizothoracinae shows separate line of evolution from subfamily Barbinae.

Hybridization between an endangered freshwater fish and an introduced congeneric species and consequent genetic introgression

Artificial transplantation of organisms and consequent invasive hybridization can lead to the extinction of native species. In Matsuyama, Japan, a native bitterling fish, Tanakia lanceolata, is known to form hybrids with another bitterling species, T. limbata, which was recently introduced from western Kyushu, Japan. These bitterlings spawn in the gills of two freshwater unionid species, Pronodularia japanensis and Nodularia douglasiae nipponensis, which have rapidly declined on the Matsuyama Plain in the past 30 years. To gauge the effect of invasive hybridization, we determined the genetic introgression between T. lanceolata and T. limbata and analyzed the morphology of these species and their hybrids to infer their niche overlap. We collected adult individuals of Tanakia spp. and genotyped them based on six microsatellite loci and mitochondrial cytochrome b sequences. We analyzed their meristic characters and body shapes by geometric morphometrics. We found that 10.9% of all individuals collected were hybrids. Whereas T. lanceolata were more densely distributed downstream and T. limbata were distributed upstream, their hybrids were widely distributed, covering the entire range of native T. lanceolata. The body height and anal fin length of T. limbata were greater than those of T. lanceolata, but their hybrids were highly morphologically variable, covering both parental morphs, and were widely distributed in the habitats of both parental species. Hybridization has occurred in both directions, but introduced T. limbata females and native T. lanceolata males are more likely to have crossed. This study shows that invasive hybridization with the introduced T. limbata is a potential threat to the native population of T. lanceolata via genetic introgression and replacement of its niche in streams.

Unusual AT-skew of Sinorhodeus microlepis mitogenome provides new insights into mitogenome features and phylogenetic implications of bitterling fishes

Sinorhodeus microlepis (S. microlepis) is recently described as a new species and represents a new genus Sinorhodeu of the subfamily Acheilognathinae. In this study, we first sequenced the complete mitogenome of S. microlepis and compared with the other 29 bitterling mitogenomes. The S. microlepis mitogenome is 16,591 bp in length and contains 37 genes. Gene distribution pattern is identical among 30 bitterling mitogenomes. A significant linear correlation between A+T% and AT-skew were found among 29 bitterling mitogenomes, except S. microlepis shows unusual AT-skew with slightly negative in tRNAs and PCGs. Bitterling mitogenomes exhibit highly conserved usage bias of start codon, relative synonymous codons and amino acids, overlaps and non-coding intergenic spacers. Phylogenetic trees constructed by 13 PCGs strongly support the polyphyly of the genus Acheilognathus and the paraphyly of Rhodeus and Tanakia. Together with the unusual characters of S. microlepis mitogenomes and phylogenetic trees, S. microlepis should be a sister species to the genus Rhodeu that might diverge about 13.69 Ma (95% HPD: 12.96-14.48 Ma).

Hybridization between two bitterling fish species in their sympatric range and a river where one species is native and the other is introduced

The distributions of two bitterling fish (subfamily: Acheilognathinae), Tanakia lanceolata and T. limbata, overlap in western Japan. Acheilognathinae fish lay their eggs in the gills of freshwater bivalves, and the early juvenile stage develops in the gills. Populations of freshwater bivalves are declining worldwide, which has limited the number of spawning substrate for bitterlings. T. limbata has been artificially introduced to some rivers in Ehime, Japan, where it coexists with native T. lanceolata, and some hybrids have been observed. We collected both species from several sites in western Japan, and from the Kunichi River system in Ehime, and analyzed genetic population structure based on six microsatellite loci and sequences of the mitochondrial cytochrome b gene. Structure analysis identified three genetically distinct populations: T. lanceolata, T. limbata “West Kyushu”, and T. limbata “Setouchi”. Two clades of T. limbata were also supported by molecular phylogenetic analyses based on cytochrome b. Hybrids in Ehime originated mostly from interbreeding between male T. lanceolata and female T. limbata “West Kyushu”, and made up 10.2% of all collected fish, suggesting that hybrids occurred frequently between females of colonizing species and males of native species. On the other hand, interspecific hybrids were detected at rates of 40.0%, 20.0%, and 17.6% in the Ima River (Fukuoka), Midori River (Kumamoto), and Kase River (Saga), respectively, which are naturally sympatric regions. We found a few T. limbata “Setouchi” in the Midori and Kase Rivers, which were supposed to be introduced from other regions, coexisting with native T. limbata “West Kyushu”, and this cryptic invasion may have triggered the interspecific hybridization. These results suggest that artificial introduction of a fish species, a decline in the unionid population, and degradation of habitat have caused broad hybridization of bitterlings in western Japan.

Mitochondrial genome of Paracheilognathus imberbis (Cypriniformes: Cyprinidae: Acheilognathinae)

In this study, the complete mitochondrial genome (mitogenome) sequence of Paracheilognathus imberbis (Cypriniformes: Cyprinidae: Acheilognathinae) was determined by long PCR and primer-walking methods. The complete mitochondrial genome is 16,819 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes as well as a displacement loop (D-loop). The overall base composition of the genome is A(29.73%), T(27.25%), C(26.63%) and G(17.40%), respectively. The mitogenome of P. imberbis displayed novel gene order arrangement compared with published Rhodeus sinensis to date. The mitogenome would contribute to resolving phylogenetic position and interrelationships of Acheilognathinae.

The complete mitogenome of Acheilognathus barbatus (Cypriniformes; Cyprinidae)

The complete mitochondrial genome was sequenced from the freshwater cyprinid fish, Acheilognathus barbatus. The mitochondrial genome sequence was 16,770 bp in size, and the gene order and contents were identical with the congeneric species A. macropterus, A. typus and A. yamatsutae. Six genes (COII, ND2, COIII, ND3, ND4 and Cytb) had an incomplete stop codon as shown in other cyprinid fishes. Base composition of the genome is A (29.1%), T (26.9%), C (26.6%) and G (17.4%) with an A + T rich feature (56%) as that of other vertebrate mitochondrial genomes.

Phylogenetic relationships of Acheilognathidae (Cypriniformes: Cyprinoidea) as revealed from evidence of both nuclear and mitochondrial gene sequence variation: evidence for necessary taxonomic revision in the family and the identification of cryptic species

Bitterlings are relatively small cypriniform species and extremely interesting evolutionarily due to their unusual reproductive behaviors and their coevolutionary relationships with freshwater mussels. As a group, they have attracted a great deal of attention in biological studies. Understanding the origin and evolution of their mating system demands a well-corroborated hypothesis of their evolutionary relationships. In this study, we provide the most comprehensive phylogenetic reconstruction of species relationships of the group based on partitioned maximum likelihood and Bayesian methods using DNA sequence variation of nuclear and mitochondrial genes on 41 species, several subspecies and three undescribed species. Our findings support the monophyly of the Acheilognathidae. Two of the three currently recognized genera are not monophyletic and the family can be subdivided into six clades. These clades are further regarded as genera based on both their phylogenetic relationships and a reappraisal of morphological characters. We present a revised classification for the Acheilognathidae with five genera/lineages: Rhodeus, Acheilognathus (new constitution), Tanakia (new constitution), Paratanakia gen. nov., and Pseudorhodeus gen. nov. and an unnamed clade containing five species currently referred to as "Acheilognathus". Gene trees of several bitterling species indicate that the taxa are not monophyletic. This result highlights a potentially dramatic underestimation of species diversity in this family. Using our new phylogenetic framework, we discuss the evolution of the Acheilognathidae relative to classification, taxonomy and biogeography.