Evidence from mitochondrial genomics supports the lower Mesozoic of South Asia as the time and place of basal divergence of cypriniform fishes (Actinopterygii: Ostariophysi)

Naming the other cousin: A new goldie barb (Cyprinidae: Smiliogastrininae) from the northeast escarpment in South Africa, with proposed taxonomic rearrangement of the goldie barb group in southern Africa

A growing body of evidence indicates that the global diversity of freshwater fishes has not been fully documented. Studies of freshwater fishes that were previously thought to be morphologically variable have revealed the existence of deeply divergent lineages, with many distinct species. In southern Africa a number of Enteromius species exhibit either exceedingly wide or divided distribution patterns that should be rare for freshwater fishes with limited dispersal opportunities between river systems. One such species is the sidespot barb, Enteromius neefi. As currently defined, E. neefi has a disjunct distribution that is divided between rivers in the northeast escarpment in South Africa and Eswatini, and tributaries of the Upper Zambezi in Zambia and southern Congo in the Democratic Republic of Congo, with a large geographic gap between these two populations. With the use of molecular and morphological methods, the level of divergence between the two populations was examined, and a new species was described from the Steelpoort River in the Limpopo River system of South Africa. Findings from this study provide further evidence for a number of taxonomic problems within the goldie barbs of southern Africa, and some taxonomic rearrangements are proposed for this group.

Pleistocene Landscape Dynamics Drives Lineage Divergence of a Temperate Freshwater Fish Gobio rivuloides in Coastal Drainages of Northern China

Understanding historical processes underlying lineage distribution patterns is a primary goal of phylogeography. We selected Gobio rivuloides (Cypriniformes: Gobionidae) as a model to improve our knowledge about how intraspecific genetic divergence of freshwater fishes arises in coastal drainages of northern China via statistical analysis using cytochrome b gene. The time-calibrated phylogeny of G. rivuloides showed the divergence of two major lineages (I and II) at ~0.98 Ma (million years ago). Lineage I can be divided into two sub-lineages (I-A and I-B) with a divergence time of ~0.83 Ma. Sub-lineage I-A inhabits the Amur River, and sub-lineage I-B lives in the Luan River and Liao River. Lineage II is distributed in the Yellow River and Hai River, with close genetic relationships between the two drainages, and can be split into two sub-lineages (II-C and II-D) with a divergence time of ~0.60 Ma. Our findings indicate that the splitting of lineages and sub-lineages could be attributed to geographic isolation caused by the formation of the Bohai Sea, river capture, and the episodic hydrologic closing of a paleolake during the late Lower-Middle Pleistocene. It is also the first report we know of displaying a clear phylogeographic break for freshwater fishes across coastal drainages in northern China.

The complete mitochondrial genomes of the Leuciscus baicalensis and Rutilus rutilus: a detailed genomic comparison among closely related species of the Leuciscinae subfamily

Cyprinidae is the largest family in the order of freshwater fish Cypriniformes. Increased subfamily members of Cyprinidae have been suggested to be re-classified for decades. In this study, we sequenced the mitochondrial genomes (mitogenomes) of Leuciscus baicalensis and Rutilus rutilus collected from northwest China and compared with other closely related species to determine their associated family or subfamily. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of Leuciscus baicalensis and Rutilus rutilus and characterized the mitogenomes by the gene structure, gene order, and the secondary structures of the 22 tRNA genes. We compared mitogenome features of Leuciscinae with other subfamilies in Cyprinidae. We used the analytic Bayesian Information and Maximum Likelihood methods to determine phylogenetic trees of 13 PCGs. The mitogenomes of Leuciscus baicalensis and Rutilus rutilus were 16,607 bp and 16,606 bp, respectively. Organization and location of these genes were consistent with already studied Leuciscinae fishes. Synonymous codon usage was conservative in Leuciscinae as compared with other subfamilies in Cyprinidae. Phylogenetic analysis indicated that Leuciscinae was a monophyletic group, and genus Leuciscus was a paraphyletic group. Our approach, for the first time, of studying comparative mitochondrial genomics and phylogenetics together provided a supportive platform to the analysis of population genetics and phylogeny for Leuciscinae. Our results indicated a promising potential of comparative mitochondrial genomics in the manifestation of phylogenetic relationships between fishes, leading us to a suggestion that mitogenomes should be routinely considered in clarifying phylogenetics of family and subfamily members of fish.

Integrative approach for discovering of the new species within the genus Allocreadium Looss, 1900 (Trematoda: Allocreadiidae) and framing of biogeographical hypotheses for the genus

In July 2012 new Allocreadium species was isolated from Carassius gibelio caught in the Arsenyevka River, Primorsky region, Russia. Analyses on the morphometrics and internal organs' topology revealed that these worms are morphologically closest with A. isoporum but both species are independent according to high genetic distances based on the 28S gene fragment (5.434±0.0073%). Unlike A. isoporum found earlier in Europe, the new species named A. pseudoisoporum sp. nov. has a shorter body length and the vitellarium not reaching the posterior end of the body at some distance and its anterior border is on the level of the ventral sucker. Allocreadium pseudoisoporum sp. nov. differs from seven species previously found in the Russian Far East with the following features: smaller size of the body, suckers' ratio, range values, and topology of internal structures. Newly localities in the Pavlovka River and the Artyomovka River were discovered for A. khankaiensis. Morphological variability of the worms from the Pavlovka River was observed in comparison with A. khankaiensis from the Komissarovka River. Using scanning electron microscope, we examined external surfaces of three species (A. pseudoisoporum sp. nov., A. khankaiensis, A. hemibarbi) and observed structures reminiscent sensory receptors. This study was aimed to describe species diversity of allocreadiids inhabiting the south of Primorsky region, and to clarify phylogenetic relationships between the species from the genus Allocreadium Looss, 1900 using molecular genetic methods. The phylogenetic Bayesian tree based on the 28S gene showed a clear separation of ten Allocreadium species and confirmed the validity of A. pseudoisoporum sp. nov. Allocreadium pseudoisoporum sp. nov. is most similar to A. gotoi (genetic distances - 3.578±0.0051% in 28S, and 18.777±0.0149% in cox1), and represented the earliest divergent lineage in Allocreadium clade on the phylogenetic tree based on the 28S gene, thereby indicating its proximity to the ancestral node. Also, dichotomous keys for 25 Palearctic species of Allocreadium were prepared based on the morphology of the adult worms.

Determination and analysis of the complete mitochondrial genome of Barilius barila (Cypriniformes: Danionidae: Chedrinae)

Cyprinid fish Barilius barila found in the Irrawaddy water system is a valuable fishery resource and has been listed as Least Concern by the IUCN. This study determined the complete mitochondrial genome of B. barila from Yunnan, China, for the first time. Circular molecule of B. barila mitogenome was sequenced to be 16,560 bp in length, with the typical gene structure of 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and two noncoding areas (control region and the origin of L-strand replication). Overall nucleotides composition appeared to be 27.5% A, 24.8% T, 19.2% G, and 28.6% C, with a slight AT (52.3%) bias. The topology of the phylogenetic tree showed that B. barila was well grouped with Opsarius caudiocellatus, and clustered together with the genus Opsarius instead of Barilius, revealing that it was more reasonable for Barilius barila to belong to Opsarius rather than Barilius.

The minnow Phoxinus lumaireul (Leuciscidae) shifts the Adriatic-Black Sea basin divide in the north-western Dinaric Karst region

Karst landscapes are characterized by intermittent and sinking streams. The most common method used to study underground hydrological connections in karst is tracing tests. However, a more biologically oriented approach has been suggested: analysis of the genetic structure of aquatic organisms. Biological tracers can be sought among trogloxenes, that is, surface species that occasionally enter caves and groundwater. One such example is the fish genus Phoxinus, which exhibits high genetic diversity and complex phylogeography in the Balkan Peninsula. In the north-western Dinaric Karst, the complex hydrological network was digitalized in 2020. Contemporaneously, Phoxinus lumaireul populations in the Slovenian Dinaric Karst were intensively sampled and analysed for fragments of two mitochondrial genes and one nuclear gene. The derived phylogeographic structure and data on hydrological connections were compared to evaluate support for three alternative scenarios: The genetic structure (1) is a consequence of the ongoing geneflow through underground connections, (2) reflects a previous hydrological network or (3) is an outcome of anthropogenic translocations. The results suggest that the first two scenarios seem to have played a major role, while the third has not had profound effects on the genetic composition. Comparison between the genetic structure of Slovenian Dinaric Karst sampling sites and that of hydrologically isolated reference sampling sites indicated a greater genetic connectivity in the former. Moreover, the range of Adriatic (1a) and Black Sea (1c) haplotypes does not correspond to the Adriatic-Black Sea basin divide but is shifted northwards.

Unauthorized Stocking of an Endangered Bitterling Acheilognathus typus in an Irrigation Pond Detected and Substantiated by Biological and Human Lines of Evidence

Acheilognathus typus, an endangered bitterling, was captured in an irrigation pond in the northern part of Niigata prefecture, Japan, in 2019. This bitterling species had once occupied that region. Its absence for years indicated the possible extinction of the bitterling there. We expected that the recently captured individuals are an unknown remnant stock of that endangered species found through an extensive survey. Mitochondrial genotyping, however, revealed that the recently captured individuals had a common haplotype with those from Kashimadai, Miyagi prefecture, Japan. The uniqueness of the haplotype to these two ponds indicates a stock identity between them. In the pond in Kashimadai, the illegal activity of catching the bitterlings by a person from the northern Niigata prefecture was detected in 2015. We conclude that the bitterling from the northern Niigata prefecture was stocked from a pond in Kashimadai. Our report is the first example of unauthorized bitterling stocking substantiated from both biological and human lines of evidence.

Phylogenetic relationships, origin and historical biogeography of the genus Sprattus (Clupeiformes: Clupeidae)

The genus Sprattus comprises five species of marine pelagic fishes distributed worldwide in antitropical, temperate waters. Their distribution suggests an ancient origin during a cold period of the earth’s history. In this study, we evaluated this hypothesis and corroborated the non-monophyly of the genus Sprattus, using a phylogenetic approach based on DNA sequences of five mitochondrial genome regions. Sprattus sprattus is more closely related to members of the genus Clupea than to other Sprattus species. We also investigated the historical biogeography of the genus, with the phylogenetic tree showing two well-supported clades corresponding to the species distribution in each hemisphere. Time-calibrated phylogenetic analyses showed that an ancient divergence between Northern and Southern Hemispheres occurred at 55.8 MYBP, followed by a diversification in the Oligocene epoch in the Northern Hemisphere clade (33.8 MYBP) and a more recent diversification in the Southern Hemisphere clade (34.2 MYBP). Historical biogeography analyses indicated that the most recent common ancestor (MRCA) likely inhabited the Atlantic Ocean in the Southern Hemisphere. These results suggest that the ancestral population of the MRCA diverged in two populations, one was dispersed to the Northern Hemisphere and the other across the Southern Hemisphere. Given that the Eocene was the warmest epoch since the Paleogene, the ancestral populations would have crossed the tropics through deeper cooler waters, as proposed by the isothermal submergence hypothesis. The non-monophyly confirmed for the genus Sprattus indicates that its systematics should be re-evaluated.

Morphological and molecular characterization of myxobolids (Cnidaria, Myxozoa) infecting cypriniforms (Actinopterygii, Teleostei) endemic to the Iberian Peninsula

The Iberian Peninsula provides a unique freshwater ecosystem for native and endemic cypriniforms to thrive. Despite cypriniforms being hosts to multiple myxobolids worldwide, little research has been performed in this geographic location. In this study, the examination of three Iberian endemic cypriniforms showed that myxosporean richness in the Iberian Peninsula is underestimated, with three new and one known myxobolid species being reported based on morphological and molecular data (SSU). Myxobolus arcasii n. sp. is described from the kidney and gonads of the "bermejuela" Achondrostoma arcasii, M. duriensis n. sp. from the gills of the Northern straight-mouth nase Pseudochondrostoma duriense, and Thelohanellus paludicus n. sp. from the intestine of the Southern Iberian spined-loach Cobitis paludica. Myxobolus pseudodispar Gorbunova, 1936 is further reported from several organs of P. duriense, and from the spleen of A. arcasii. The occurrence of M. pseudodispar in endemic Iberian species reveals that host-shift followed its co-introduction with central European leuciscids into this geographic location. Several other myxobolids originally described from barbels in central Europe have also been reported from the Iberian endemic cypriniform Luciobarbus bocagei. Nonetheless, except for M. musculi, the identification of these myxobolids in L. bocagei is here shown to be dubious and require molecular confirmation. Phylogenetic analyses reveal M. arcasii n. sp. and M. duriensis n. sp. clustering within different lineages of leuciscid-infecting species, showing that myxobolids entered Leuciscidae as hosts multiple times during their evolution. Constituting the first myxobolid reported from the subfamily Cobitinae, Thelohanellus paludicus n. sp. stands alone in the tree topology.

Phylogenetic relationships and classification of the Holarctic family Leuciscidae (Cypriniformes: Cyprinoidei)

The phylogenetic relationships and classification of the freshwater fish order Cypriniformes, like many other species-rich groups of vertebrates, has evolved over time with some consistency and inconsistencies of relationships across various studies. Within Cypriniformes, the Holarctic family Leuciscidae is one of the most widely distributed and highly diverse monophyletic groups of cyprinoids. Despite several studies conducted on this group, alternative hypotheses exist as to the composition and relationships within Leuciscidae. Here we assess the extent, composition, phylogenetic relationships, and taxonomy of this highly diverse group of fishes, using multiple mitochondrial and nuclear loci and a comprehensive and dense taxonomic sampling. Analyses of 418 specimens (410 species) resolve a well-supported Leuciscidae including 362 specimens (358 taxa) in six well-supported subfamilies/major clades: Pseudaspininae/Far East Asian clade (FEA); Laviniinae/North American Western clade (WC); Plagopterinae/North American Creek Chub-Plagopterin clade (CC-P); Leuciscinae/Eurasian Old World clade (OW) (minus Phoxinus) plus North American Notemigonus; Phoxininae/Eurasian Phoxinus clade (PHX); and Pogonichthyinae/North American clade (NA) including all remaining leuciscids. Within Leuciscidae, neither the traditional phoxinins (Phoxinus, FEA, Nearctic genera) nor all Nearctic genera (minus Notemigonus) are resolved as monophyletic; whereas the WC and CC-P form two independent lineages from remaining North American cyprinoids. A close relationship exists between Eurasian Phoxinus, NA, and OW clades, while FEA is the sister group to all remaining Leuciscidae. Major lineages resolved within these six subfamilies are mostly congruent with some previous studies. Our results suggests a complex evolutionary history of this diverse and widespread group of fishes.

Comparative transcriptomics of cyprinid minnows and carp in a common wild setting: a resource for ecological genomics in freshwater communities

Comparative transcriptomics can now be conducted on organisms in natural settings, which has greatly enhanced understanding of genome–environment interactions. Here, we demonstrate the utility and potential pitfalls of comparative transcriptomics of wild organisms, with an example from three cyprinid fish species (Teleostei:Cypriniformes). We present extensively filtered and annotated transcriptome assemblies that provide a valuable resource for studies of genome evolution (e.g. polyploidy), ecological and morphological diversification, speciation, and shared and unique responses to environmental variation in cyprinid fishes. Our results and analyses address the following points: (i) ‘essential developmental genes’ are shown to be ubiquitously expressed in a diverse suite of tissues across later ontogenetic stages (i.e. juveniles and adults), making these genes are useful for assessing the quality of transcriptome assemblies, (ii) the influence of microbiomes and other exogenous DNA, (iii) potentially novel, species-specific genes, and (iv) genomic rearrangements (e.g. whole genome duplication). The data we present provide a resource for future comparative work in cypriniform fishes and other taxa across a variety of sub-disciplines, including stress response, morphological diversification, community ecology, ecotoxicology, and climate change.

Phylogenomic Perspective on the Relationships and Evolutionary History of the Major Otocephalan Lineages

The phylogeny of otocephalan fishes is the subject of broad controversy based on morphological and molecular evidence. The primary unresolved issue pertaining to this lineage relates to the origin of Characiphysi, especially the paraphyly of Characiformes. The considerable uncertainty associated with this lineage has precluded a greater understanding of the origin and evolution of the clade. Herein, a phylogenomic approach was applied to resolve this debate. By analyzing 10 sets of transcriptomic data generated in this study and 12 sets of high-throughput data available in public databases, we obtained 1,110 single-copy orthologous genes (935,265 sites for analysis) from 22 actinopterygians, including 14 otocephalan fishes from six orders: Clupeiformes, Gonorynchiformes, Cypriniformes, Siluriformes, Characiformes, and Gymnotiformes. Based on a selection of 125 nuclear genes screened from single-gene maximum likelihood (ML) analyses and sequence bias testing, well-established relationships among Otocephala were reconstructed. We suggested that Gymnotiformes are more closely related to Characiformes than to Siluriformes and Characiformes are possibly paraphyletic. We also estimated that Otocephala originated in the Early-Late Jurassic, which postdates most previous estimations, and hypothesized scenarios of the early historical biogeographies of major otocephalan lineages.

Evolution of the connection patterns of the cephalic lateral line canal system and its use to diagnose opsariichthyin cyprinid fishes (Teleostei, Cyprinidae)

The cephalic lateral line canal systems were compared among 12 species of the cyprinid tribe Opsariichthyini. All species were characterized by the separation of the supraorbital canal from both the infraorbital and the temporal canals, and the left side of the supratemporal canal from the right side of the canal. In species of Candidia, Opsariichthys, Parazacco, and Zacco, and Nipponocyprissieboldii the temporal canal was separated from the preoperculomandibular canal. In Nipponocypristemminckii and N.koreanus, the temporal canal was connected to the preoperculomandibular canal. Separation of the left and right sides of the supratemporal canal is a possible synapomorphy of the opsariichthyin cyprinids. Opsariichthysuncirostris and O.bidens are unique among the opsariichthyins in that the connection between the infraorbital and temporal canals is retarded. The variation in arrangement of the cephalic lateral line canal system can be used as diagnostic characters for the opsariichthyin species.

Rhodopsin gene copies in Japanese eel originated in a teleost-specific genome duplication

BACKGROUND

Gene duplication is considered important to increasing the genetic diversity in animals. In fish, visual pigment genes are often independently duplicated, and the evolutionary significance of such duplications has long been of interest. Eels have two rhodopsin genes (rho), one of which (freshwater type, fw-rho) functions in freshwater and the other (deep-sea type, ds-rho) in marine environments. Hence, switching of rho expression in retinal cells is tightly linked with eels' unique life cycle, in which they migrate from rivers or lakes to the sea. These rho genes are apparently paralogous, but the timing of their duplication is unclear due to the deep-branching phylogeny. The aim of the present study is to elucidate the evolutionary origin of the two rho copies in eels using comparative genomics methods.

RESULTS

In the present study, we sequenced the genome of Japanese eel Anguilla japonica and reconstructed two regions containing rho by de novo assembly. We found a single corresponding region in a non-teleostean primitive ray-finned fish (spotted gar) and two regions in a primitive teleost (Asian arowana). The order of ds-rho and the neighboring genes was highly conserved among the three species. With respect to fw-rho, which was lost in Asian arowana, the neighboring genes were also syntenic between Japanese eel and Asian arowana. In particular, the pattern of gene losses in ds-rho and fw-rho regions was the same as that in Asian arowana, and no discrepancy was found in any of the teleost genomes examined. Phylogenetic analysis supports mutual monophyly of these two teleostean synteny groups, which correspond to the ds-rho and fw-rho regions.

CONCLUSIONS

Syntenic and phylogenetic analyses suggest that the duplication of rhodopsin gene in Japanese eel predated the divergence of eel (Elopomorpha) and arowana (Osteoglossomorpha). Thus, based on the principle of parsimony, it is most likely that the rhodopsin paralogs were generated through a whole genome duplication in the ancestor of teleosts, and have remained till the present in eels with distinct functional roles. Our result indicates, for the first time, that teleost-specific genome duplication may have contributed to a gene innovation involved in eel-specific migratory life cycle.

The complete mitochondrial genome of Cichla ocellaris

In this study, the complete mitochondrial (mt) genome of Cichlao cellaris was determined. The mt genome has a length of 16,526 bp and encodes 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and one AT-rich non-coding region (D-loop). The gene arrangement was similar to those of typical fishes. The total base composition of the mt genome was 25.1% T, 30.8% C, 29.3% A, and 14.8% G. Of the 13 PCGs, 12 genes start with an ATG codon, except COX1 is with GTG. And among these 13 genes, seven (ND1, COX1, ATP8, ND4L, ND5, ND6, and Cytb) used TAA or TAG as the termination codon, whereas six (ND2, COX2, ATP6, COX3, ND3, and ND4) have incomplete stop codon T. Its control region is atypical in being short at 850 bp. This mt genome sequence data will be useful for phylogenetic and systematic analyses within the family Cichlaidae.

Barcode index number, taxonomic rank and modes of speciation: examples from fish

Species delimitation by DNA sequence data or DNA barcoding is successful, as confirmed by the vast BOLD data base. However, the theory that would explain this fact has not been developed yet. An approach based on Barcoding Index Number (BIN), suggested in the assignment, allows delimiting of taxa of three ranks (species, genera, and families) and statistical validation with a high precision of delimiting (over 80%), as well as shows for majority of Co-1-based single gene trees good correspondence between their topology and conventional taxa content for analyzed fish species (R² ≈ 0.84-0.98). Knowledge of deviations from these data can help to find out new taxa and improve biodiversity description. It is concluded that delimiting is successful for bulk of cases because the geographic mode of speciation prevails in nature. It takes a long time for new taxa to form in isolation, which allows accumulation of random mutations and many different nucleotide substitutions between them that can be detected by molecular markers and give unique DNA barcodes. The use of BIN approach, described here, can aid greatly in making this important question clearer especially under wider examination of other organisms.

Complete mitochondrial genome of the jellyfish, Nemopilema nomurai (Cnidaria: Scyphozoa) and the phylogenetic relationship in the related species

The complete mitochondrial genome of giant jellyfish Nemopilema nomurai collected from East China Sea was determined by next-generation sequencing. The mitogenome is a circular molecule 17,024 bp in length, including 13 protein-coding genes (including Cox 1, Cox2, Atp 8, Atp 6, ND5, ND 6, ND3, ND41, ND1,ND4, Cytb), 6 tRNAs (tRNA-Trp, tRNA-Met, tRNA-Val, tRNA-Arg, tRNA-Glu, tRNA-Asn), 2 rRNA genes (small subunit RNA and large subunit RNA), and 1 putative control region. The phylogenetic tree in the related species showed that giant jellyfish is close to Hydraoligaclis.

Complete mitochondrial genome sequence of E. suratensis revealed by next generation sequencing

The complete mitochondrial genome of Etroplus suratensis, the Green chromide cichlid, was determined for the first time through NGS method. The genome is 16,467 bp (Accession no. KU301747) in length and consisted of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and one control region. Organization of genes and their order are in accordance with other vertebrates. The overall base composition on plus strand was A: 28.3%, G: 15.2%, C: 30.9%, T: 25.6%, and the A + T content 53.9%. The control region contains a putative termination-associated sequence and three conserved sequence blocks. This mitogenome sequence data would play an important role in population genetics and phylogenetics of cichlid fish of India.

Resolving Cypriniformes relationships using an anchored enrichment approach

BACKGROUND

Cypriniformes (minnows, carps, loaches, and suckers) is the largest group of freshwater fishes in the world (~4300 described species). Despite much attention, previous attempts to elucidate relationships using molecular and morphological characters have been incongruent. In this study we present the first phylogenomic analysis using anchored hybrid enrichment for 172 taxa to represent the order (plus three out-group taxa), which is the largest dataset for the order to date (219 loci, 315,288 bp, average locus length of 1011 bp).

RESULTS

Concatenation analysis establishes a robust tree with 97 % of nodes at 100 % bootstrap support. Species tree analysis was highly congruent with the concatenation analysis with only two major differences: monophyly of Cobitoidei and placement of Danionidae.

CONCLUSIONS

Most major clades obtained in prior molecular studies were validated as monophyletic, and we provide robust resolution for the relationships among these clades for the first time. These relationships can be used as a framework for addressing a variety of evolutionary questions (e.g. phylogeography, polyploidization, diversification, trait evolution, comparative genomics) for which Cypriniformes is ideally suited.

Mitogenomic perspectives on the origin of Tibetan loaches and their adaptation to high altitude

Tibetan loaches are the largest group of Tibetan fishes and are well adapted to the Tibetan Plateau. To investigate the origin of Tibetan loaches and their adaptations to the Tibetan Plateau, we determined 32 complete mitochondrial genomes that included 29 Tibetan loach species, two Barbatula species and Schistura longus. By combining these newly determined sequences with other previously published mitochondrial genomes, we assembled a large mitogenomic data set (11,433 bp) of 96 species in the superfamily Cobitoidea, to investigate the phylogenetic status of the genus Triplophysa. The resulting phylogeny strongly supported that the genus Triplophysa forms a monophyletic group within Nemacheilidae. Our molecular dating time suggests that the lineage leading to the Tibetan loaches and other loaches diverged approximately 23.5 Ma, which falls within the period of recent major uplifts of the Tibetan Plateau in the Early Miocene. Selection analyses revealed that the mitochondrial protein-coding genes of Tibetan loaches have larger ratios of nonsynonymous to synonymous substitutions than do those of non-Tibetan loaches, indicating that Tibetan loaches accumulated more nonsynonymous mutations than non-Tibetan loaches and exhibited rapid evolution. Two positively selected sites were identified in the ATP8 and ND1 genes.

Molecular systematics and DNA barcoding of Altai osmans, oreoleuciscus (pisces, cyprinidae, and leuciscinae), and their nearest relatives, inferred from sequences of cytochrome b (Cyt-b), cytochrome oxidase c (Co-1), and complete mitochondrial genome

Mitochondrial DNA (mtDNA) at the protein-coding Cyt-b gene along with data retrieved from GenBank for Co-1 gene fragments and complete mitochondrial genome (mitogenome) of Altai osmans and the nearest relatives of Leuciscinae fish species were compared for the estimation of variability and phylogenetic tree building. Phylogenetic trees were built by four techniques: Bayesian (BA), maximum likelihood (ML), maximum parsimony (MP), and neighbor-joining (NJ). Resolution of Cyt-b trees for species of two genera (Oreoleuciscus and Phoxinus) was quite distinct at all the approaches. For Tribolodon, the single gene trees were not well resolved; however, the mitogenome tree was resolved. Species identification on per individual basis (DNA barcoding) was high for both Cyt-b and Co-1 genes. The trees built using the data for 13 protein mitochondrial genes revealed a complicated phylogenetic pattern within the subfamily Leuciscinae. Scores of the average p-distances at three taxonomic levels were considerably different: (1) 1.16 ± 0.96, (2) 8.21 ± 1.01, and (3) 16.41 ± 0.85 for Cyt-b and (1) 1.04 ± 0.78, (2) 8.30 ± 0.92, and (3) 10.74 ± 0.79 for 13 protein genes of mitogenome, where (1) is intraspecies, (2) is intragenus, and (3) is intrasubfamily levels. Data on mitogenome distances were summarized for the taxonomic hierarchy for the first time. A concordant increase in distance score with growth of the rank of taxa (having the minimum score at the intraspecies level), both for a single gene and the whole mitogenome, substantiates the concept that speciation in the subfamily Leuciscinae in most cases follows the geographic mode. The distinct clustering of Altai osmans, Oreoleuciscus potanini and O. humilis, in the Cyt-b and Co-1 gene trees with small overall genetic distances, obtained for both genes, allows us to consider these taxa as separate but genetically sister species.

Assessment of yellowtail kingfish (Seriola lalandi) as a surrogate host for the production of southern bluefin tuna (Thunnus maccoyii) seed via spermatogonial germ cell transplantation

Germ cell transplantation is an innovative technology for the production of interspecies surrogates, capable of facilitating easier and more economical management of large-bodied broodstock, such as the bluefin tuna. The present study explored the suitability of yellowtail kingfish (Seriola lalandi) as a surrogate host for transplanted southern bluefin tuna (Thunnus maccoyii) spermatogonial cells to produce tuna donor-derived gametes upon sexual maturity. Germ cell populations in testes of donor T. maccoyii males were described using basic histology and the molecular markers vasa and dead-end genes. The peripheral area of the testis was found to contain the highest proportions of dead-end-expressing transplantable Type A spermatogonia. T. maccoyii Type A spermatogonia-enriched preparations were transplanted into the coelomic cavity of 6–10-day-old post-hatch S. lalandi larvae. Fluorescence microscopy and polymerase chain reaction analysis detected the presence of tuna cells in the gonads of the transplanted kingfish fingerlings at 18, 28, 39 and 75 days after transplantation, indicating that the transplanted cells migrated to the genital ridge and had colonised the developing gonad. T. maccoyii germ cell-derived DNA or RNA was not detected at later stages, suggesting that the donor cells were not maintained in the hosts’ gonads.

Karyotype differentiation in 19 species of river loach fishes (Nemacheilidae, Teleostei): extensive variability associated with rDNA and heterochromatin distribution and its phylogenetic and ecological interpretation

BACKGROUND

Loaches of the family Nemacheilidae are one of the most speciose elements of Palearctic freshwater ichthyofauna and have undergone rapid ecological adaptations and colonizations. Their cytotaxonomy is largely unexplored; with the impact of cytogenetical changes on this evolutionary diversification still unknown. An extensive cytogenetical survey was performed in 19 nemacheilid species using both conventional (Giemsa staining, C- banding, Ag- and Chromomycin A3/DAPI stainings) and molecular (fluorescence in situ hybridization with 5S rDNA, 45S rDNA, and telomeric (TTAGGG)n probes) methods. A phylogenetic tree of the analysed specimens was constructed based on one mitochondrial (cytochrome b) and two nuclear (RAG1, IRBP) genes.

RESULTS

Seventeen species showed karyotypes composed of 2n = 50 chromosomes but differentiated by fundamental chromosome number (NF = 68-90). Nemachilichthys ruppelli (2n = 38) and Schistura notostigma (2n = 44-48) displayed reduced 2n with an elevated number of large metacentric chromosomes. Only Schistura fasciolata showed morphologically differentiated sex chromosomes with a multiple system of the XY1Y2 type. Chromomycin A3 (CMA3)- fluorescence revealed interspecific heterogeneity in the distribution of GC-rich heterochromatin including its otherwise very rare association with 5S rDNA sites. The 45S rDNA sites were mostly located on a single chromosome pair contrasting markedly with a pattern of two (Barbatula barbatula, Nemacheilus binotatus, N. ruppelli) to 20 sites (Physoschistura sp.) of 5S rDNA. The cytogenetic changes did not follow the phylogenetic relationships between the samples. A high number of 5S rDNA sites was present in species with small effective population sizes.

CONCLUSION

Despite a prevailing conservatism of 2n, Nemacheilidae exhibited a remarkable cytogenetic variability on microstructural level. We suggest an important role for pericentric inversions, tandem and centric fusions in nemacheilid karyotype differentiation. Short repetitive sequences, genetic drift, founder effect, as well as the involvement of transposable elements in the dispersion of ribosomal DNA sites, might also have played a role in evolutionary processes such as reproductive isolation. These remarkable dynamics of their genomes qualify river loaches as a model for the study of the cytogenetic background of major evolutionary processes such as radiation, endemism and colonization of a wide range of habitats.

An unusual mitochondrial genome structure of the tonguefish, Cynoglossus trigrammus: Control region translocation and a long additional non-coding region inversion

Flatfishes (Pleuronectiformes) exhibit different types of large-scale gene rearrangements. In the present study, the mitochondrial (mt) genome (18,369bp) of a tonguefish, Cynoglossus trigrammus, was determined using de novo mitochondrion genome sequencing. Compared with other flatfishes, the mt genome of C. trigrammus revealed distinct mitogenome architectures that primarily included two striking findings: 1) insertion of an additional long non-coding region (1647bp) making it the second largest genome length among Pleuronectiformes and 2) the translocation of the control region. The reconstructed phylogenetic tree based on 13 mt protein-coding gene sequences recovered the monophyletic suborder Pleuronectoidei and the family Cynoglossidae. These data provide useful information for a better understanding of the mitogenomic diversities and evolution in fish as well as novel genetic markers for studying population genetics and species identification.

'Fish' (Actinopterygii and Elasmobranchii) diversification patterns through deep time

Actinopterygii (ray-finned fishes) and Elasmobranchii (sharks, skates and rays) represent more than half of today's vertebrate taxic diversity (approximately 33000 species) and form the largest component of vertebrate diversity in extant aquatic ecosystems. Yet, patterns of 'fish' evolutionary history remain insufficiently understood and previous studies generally treated each group independently mainly because of their contrasting fossil record composition and corresponding sampling strategies. Because direct reading of palaeodiversity curves is affected by several biases affecting the fossil record, analytical approaches are needed to correct for these biases. In this review, we propose a comprehensive analysis based on comparison of large data sets related to competing phylogenies (including all Recent and fossil taxa) and the fossil record for both groups during the Mesozoic-Cainozoic interval. This approach provides information on the 'fish' fossil record quality and on the corrected 'fish' deep-time phylogenetic palaeodiversity signals, with special emphasis on diversification events. Because taxonomic information is preserved after analytical treatment, identified palaeodiversity events are considered both quantitatively and qualitatively and put within corresponding palaeoenvironmental and biological settings. Results indicate a better fossil record quality for elasmobranchs due to their microfossil-like fossil distribution and their very low diversity in freshwater systems, whereas freshwater actinopterygians are diverse in this realm with lower preservation potential. Several important diversification events are identified at familial and generic levels for elasmobranchs, and marine and freshwater actinopterygians, namely in the Early-Middle Jurassic (elasmobranchs), Late Jurassic (actinopterygians), Early Cretaceous (elasmobranchs, freshwater actinopterygians), Cenomanian (all groups) and the Paleocene-Eocene interval (all groups), the latter two representing the two most exceptional radiations among vertebrates. For each of these events along with the Cretaceous-Paleogene extinction, we provide an in-depth review of the taxa involved and factors that may have influenced the diversity patterns observed. Among these, palaeotemperatures, sea-levels, ocean circulation and productivity as well as continent fragmentation and environment heterogeneity (reef environments) are parameters that largely impacted on 'fish' evolutionary history, along with other biotic constraints.

The complete mitochondrial genome sequence of Beaufortia szechuanensis (Cypriniformes, Balitoridae)

The Beaufortia szechuanensis is endemic to the upper reaches of Yangtze River in China. In this study, we successfully sequenced the mitochondrial genome of the B. szechuanensis, collected from the Qingyi River. The complete mitochondrial genome of B. szechuanensis was a circular molecule of 16,559 bp in length, containing 2 ribosomal RNA (rRNA) genes, 13 protein-coding genes, 22 transfer RNA (tRNA) genes, an origin of light-strand replication (OL) and a control region (D-loop). It is the first complete mitochondrial DNA data in the genus Beaufortia. The cytochrome c oxidase subunit I (COI) sequence of 14 species in Balitoridae was used for phylogenetic analysis. The topology demonstrated that the Balitoridae can be classified into two subfamilies, and the B. szechuanensis belongs to the subfamily Gastromyzoninae.

Complete mitochondrial genome of the Ictiobus cyprinellus (Actinopterygii: Ostariophysi)

The complete mitochondrial genome of Ictiobus cyprinellus was first determined by a polymerase chain reaction-based sequencing method in this study. The mitochondria was 16,611 bp in length, including 13 protein genes, 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and 1 control region. All genes were encoded on the heavy strain except for ND6 and eight tRNA genes. The gene arrangement of I. cyprinellus was similar to those found in the other Catostomidae species. Base composition of the heavy strain was A (29.40%), T (24.86%), C (27.96%), G (17.78%) and with A + T bias of 54.26%.

Complete mitochondrial genome sequence of Labeo fimbriatus (Bloch, 1975)

The complete mitochondrial genome of Labeo fimbriatus, the fringe-lipped carp, was determined for the first time. The genome was 16,614 bp in length and consisted of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and one control region. Organization of genes and their order are in accordance with other vertebrates. The overall base composition on the plus strand was A: 32.35%, G: 15.46%, C: 27.54%, T: 24.647%, and the A + T content 57%. The control region contains a microsatellite, (TA)12, a putative termination-associated sequence and three conserved sequence blocks. This mitogenome sequence data would play an important role in population genetics and phylogenetics of cultivable carps of India.

Wild sex in zebrafish: loss of the natural sex determinant in domesticated strains

Sex determination can be robustly genetic, strongly environmental, or genetic subject to environmental perturbation. The genetic basis of sex determination is unknown for zebrafish (Danio rerio), a model for development and human health. We used RAD-tag population genomics to identify sex-linked polymorphisms. After verifying this "RAD-sex" method on medaka (Oryzias latipes), we studied two domesticated zebrafish strains (AB and TU), two natural laboratory strains (WIK and EKW), and two recent isolates from nature (NA and CB). All four natural strains had a single sex-linked region at the right tip of chromosome 4, enabling sex genotyping by PCR. Genotypes for the single nucleotide polymorphism (SNP) with the strongest statistical association to sex suggested that wild zebrafish have WZ/ZZ sex chromosomes. In natural strains, "male genotypes" became males and some "female genotypes" also became males, suggesting that the environment or genetic background can cause female-to-male sex reversal. Surprisingly, TU and AB lacked detectable sex-linked loci. Phylogenomics rooted on D. nigrofasciatus verified that all strains are monophyletic. Because AB and TU branched as a monophyletic clade, we could not rule out shared loss of the wild sex locus in a common ancestor despite their independent domestication. Mitochondrial DNA sequences showed that investigated strains represent only one of the three identified zebrafish haplogroups. Results suggest that zebrafish in nature possess a WZ/ZZ sex-determination mechanism with a major determinant lying near the right telomere of chromosome 4 that was modified during domestication. Strains providing the zebrafish reference genome lack key components of the natural sex-determination system but may have evolved variant sex-determining mechanisms during two decades in laboratory culture.

Phylogenetic relationships of bitterling fishes (Teleostei: Cypriniformes: Acheilognathinae), inferred from mitochondrial cytochrome B sequences

Bitterling (Teleostei: Acheilognathinae) are small cyprinid fishes with a discrete distribution in East Asia and Europe. We used a complete mitochondrial cytochrome b sequence (1141 bp) from 49 species or subspecies in three genera (Tanakia, Rhodeus, and Acheilognathus), sampled across the major part of their distribution, to elucidate their phylogeny and biogeography, focusing particularly on their origin and dispersal. Based on high support value, the monophyletic Acheilognathinae separated into two major clades, Acheilognathus and Tanakia-Rhodeus. In the latter clade, the monophyly of Rhodeus was poorly supported, though it was topologically nested in Tanakia. On the basis of molecular-clock calibration, both clades diverged in the middle Miocene, with Tanakia-Rhodeus diverging slightly earlier than Acheilognathus. The Tanakia-Rhodeus clade expanded its distribution westward from the Far East, eventually reaching Europe, while Acheilognathus dispersed in the temperate regions of East Asia. A feature common to both clades is that most extant species, including Japanese endemics, appeared by the end of the Pliocene, corresponding with the present delineation of the Japanese archipelago. Autumn-spawning species with an embryonic diapause, unique to bitterling among cyprinid fishes, formed two distinct lineages (barbatulusrhombeus and longipinnis-typus) within Acheilognathus. The estimated time of divergence of the two lineages was approximately from the late Pliocene, a period characterized by glaciations. The timing of divergence suggests that the shift of spawning from spring to autumn, coupled with embryonic diapause, convergently emerged twice in the evolution of bitterling, possibly as an adaptation to the climate of the late Pliocene.

Evolutionary history of two divergent Dmrt1 genes reveals two rounds of polyploidy origins in gibel carp

Polyploidy lineages, despite very rare in vertebrates, have been proposed to play significant role in speciation and evolutionary success, but the occurrence history and consequences are still largely unknown. In this study, we used the conserved Dmrt1 to analyze polyploidy occurrence and evolutionary process in polyploid gibel carp. We identified two divergent Dmrt1 genes and respectively localized the two genes on three homologous chromosomes. Subsequently, the corresponding full-length cDNAs and genomic sequences of Dmrt1 genes were also characterized from the closely related species including Carassius auratus auratus and Cyprinus carpio, and their two Dmrt1 genes were respectively localized on two homologous chromosomes. Significantly, the evolutionary relationship analyses among cDNA and genomic DNA sequences of these Dmrt1 genes revealed two rounds of polyploidy origins in the gibel carp: an early polyploidy might result in an common tetraploid ancestor of Carassius auratus gibelio, Carassius auratus auratus and Cyprinus carpio before 18.49 million years ago (Mya), and an late polyploidy might occur from evolutionary branch of Carassius auratus at around 0.51 Mya, which lead to the occurrence of the hexaploid gibel carp. Therefore, this study provides clear genetic evidence for understanding occurrence time and historical process of polyploidy in polyploid vertebrates.

Clock gene evolution: seasonal timing, phylogenetic signal, or functional constraint?

Genetic determinants of seasonal reproduction are not fully understood but may be important predictors of organism responses to climate change. We used a comparative approach to study the evolution of seasonal timing within a fish community in a natural common garden setting. We tested the hypothesis that allelic length variation in the PolyQ domain of a circadian rhythm gene, Clock1a, corresponded to interspecific differences in seasonal reproductive timing across 5 native and 1 introduced cyprinid fishes (n = 425 individuals) that co-occur in the Rio Grande, NM, USA. Most common allele lengths were longer in native species that initiated reproduction earlier (Spearman's r = -0.70, P = 0.23). Clock1a allele length exhibited strong phylogenetic signal and earlier spawners were evolutionarily derived. Aside from length variation in Clock1a, all other amino acids were identical across native species, suggesting functional constraint over evolutionary time. Interestingly, the endangered Rio Grande silvery minnow (Hybognathus amarus) exhibited less allelic variation in Clock1a and observed heterozygosity was 2- to 6-fold lower than the 5 other (nonimperiled) species. Reduced genetic variation in this functionally important gene may impede this species' capacity to respond to ongoing environmental change.

Stability versus diversity of the dentition during evolutionary radiation in cyprinine fish

Evolutionary radiations, especially adaptive radiations, have been widely studied but mainly for recent events such as in cichlid fish or Anolis lizards. Here, we investigate the radiation of the subfamily Cyprininae, which includes more than 1300 species and is estimated to have originated from Southeast Asia around 55 Ma. In order to decipher a potential adaptive radiation, within a solid phylogenetic framework, we investigated the trophic apparatus, and especially the pharyngeal dentition, as teeth have proved to be important markers of ecological specialization. We compared two tribes within Cyprininae, Poropuntiini and Labeonini, displaying divergent dental patterns, as well as other characters related to their trophic apparatus. Our results suggest that the anatomy of the trophic apparatus and diet are clearly correlated and this explains the difference in dental patterns observed between these two tribes. Our results illustrate the diversity of mechanisms that account for species diversity in this very diverse clade: diversification of dental characters from an ancestral pattern on the one hand, conservation of a basal synapomorphy leading to ecological specialization on the other hand. By integrating morphological, ecological and phylogenetic analyses, it becomes possible to investigate ancient radiation events that have shaped the present diversity of species.

Coelacanth genomes reveal signatures for evolutionary transition from water to land

Coelacanths are known as “living fossils,” as they show remarkable morphological resemblance to the fossil record and belong to the most primitive lineage of living Sarcopterygii (lobe-finned fishes and tetrapods). Coelacanths may be key to elucidating the tempo and mode of evolution from fish to tetrapods. Here, we report the genome sequences of five coelacanths, including four Latimeria chalumnae individuals (three specimens from Tanzania and one from Comoros) and one L. menadoensis individual from Indonesia. These sequences cover two African breeding populations and two known extant coelacanth species. The genome is ∼2.74 Gbp and contains a high proportion (∼60%) of repetitive elements. The genetic diversity among the individuals was extremely low, suggesting a small population size and/or a slow rate of evolution. We found a substantial number of genes that encode olfactory and pheromone receptors with features characteristic of tetrapod receptors for the detection of airborne ligands. We also found that limb enhancers of bmp7 and gli3, both of which are essential for limb formation, are conserved between coelacanth and tetrapods, but not ray-finned fishes. We expect that some tetrapod-like genes may have existed early in the evolution of primitive Sarcopterygii and were later co-opted to adapt to terrestrial environments. These coelacanth genomes will provide a cornerstone for studies to elucidate how ancestral aquatic vertebrates evolved into terrestrial animals.

Morphological and genetic evidence for multiple evolutionary distinct lineages in the endangered and commercially exploited red lined torpedo barbs endemic to the Western Ghats of India

Red lined torpedo barbs (RLTBS) (Cyprinidae: Puntius) endemic to the Western Ghats Hotspot of India, are popular and highly priced freshwater aquarium fishes. Two decades of indiscriminate exploitation for the pet trade, restricted range, fragmented populations and continuing decline in quality of habitats has resulted in their 'Endangered' listing. Here, we tested whether the isolated RLTB populations demonstrated considerable variation qualifying to be considered as distinct conservation targets. Multivariate morphometric analysis using 24 size-adjusted characters delineated all allopatric populations. Similarly, the species-tree highlighted a phylogeny with 12 distinct RLTB lineages corresponding to each of the different riverine populations. However, coalescence-based methods using mitochondrial DNA markers identified only eight evolutionarily distinct lineages. Divergence time analysis points to recent separation of the populations, owing to the geographical isolation, more than 5 million years ago, after the lineages were split into two ancestral stocks in the Paleocene, on north and south of a major geographical gap in the Western Ghats. Our results revealing the existence of eight evolutionarily distinct RLTB lineages calls for the re-determination of conservation targets for these cryptic and endangered taxa.

Evolutionary origin and early biogeography of otophysan fishes (Ostariophysi: Teleostei)

The biogeography of the mega-diverse, freshwater, and globally distributed Otophysi has received considerable attention. This attraction largely stems from assumptions as to their ancient origin, the clade being almost exclusively freshwater, and their suitability as to explanations of trans-oceanic distributions. Despite multiple hypotheses explaining present-day distributions, problems remain, precluding more parsimonious explanations. Underlying previous hypotheses are alternative phylogenies for Otophysi, uncertainties as to temporal diversification and assumptions integral to various explanations. We reexamine the origin and early diversification of this clade based on a comprehensive time-calibrated, molecular-based phylogenetic analysis and event-based approaches for ancestral range inference of lineages. Our results do not corroborate current phylogenetic classifications of otophysans. We demonstrate Siluriformes are never sister to Gymnotiformes and Characiformes are most likely nonmonophyletic. Divergence time estimates specify a split between Cypriniformes and Characiphysi with the fragmentation of Pangea. The early diversification of characiphysans either predated, or was contemporary with, the separation of Africa and South America, and involved a combination of within- and between-continental divergence events for these lineages. The intercontinental diversification of siluroids and characoids postdated major intercontinental tectonic fragmentations (<90 Mya). Post-tectonic drift dispersal events are hypothesized to account for their current distribution patterns.

Complete mitochondrial genome sequence of Cirrhinus mrigala (Hamilton, 1822)

The complete mitochondrial genome of Cirrhinus mrigala was determined using the polymerase chain reaction. The mitogenome (16,594 bp) has the typical vertebrate mitochondrial gene arrangement, including 13 protein-coding genes, 22 tRNA genes, two rRNA genes and one control region. The overall base composition on the heavy strand was as follows: A: 32.0%, G: 15.5%, C: 28.0%, T: 24.55% and the A+T content: 56.5%. The control region contains a dinucleotide repeat motif, (TA)14, a termination-associated sequence and three conserved sequence blocks. These mitogenome sequence data would play an important role in population genetics and the molecular taxonomy of cultivable cyprinids in India.