Search

The family Cyprinidae is the largest freshwater fish group with 377 genera and over 3,000 described species. However, this group of fish has very limited cytogenetics and advanced molecular cytogenetics information. Therefore, in this study the karyotypes and other chromosomal characteristics of 15 species in the tribe Systomini (Cyprininae) were examined using Ag-NOR staining along with fluorescence in situ hybridization (5S and 18S rDNA). All species share a similar karyotype (2n = 50; NF = 88-100) in both sexes and no differentiated sex chromosome was observed. Chromosomes bearing NOR sites ranged from one to four pairs among the species, mostly mapped adjacent to telomeres in the short arms of distinct pairs in all analyzed species. This difference indicates an extensive rearrangement of chromosomes including genomic differences. The use of the 5S and 18S rDNA probe confirmed the Ag-NOR sites interstitially located in the telomeric regions of distinct chromosomes, characterizing an interspecies variation of these sites. In most of its analyzed species, the signals of 18S rDNA probe corresponded to the Ag-NOR regions, except in Barbonymus altus, B. gonionotus, B. schwanenfeldii and Puntius brevis having these signals on the same as Ag-NOR regions and other sites., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Buasriyot et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The lake minnow Eupallasella percnurus is a small leuciscid fish. In Poland, this species has been in a continuous decline since the mid-20th century and is presently considered as a extremely endangered. According to Polish law, E. percnurus is a strictly protected species that requires active conservation measures. In Poland, one the most common and effective measure of active protection E. percnurus is initiation of new populations. For this purpose, in 2004-2012, juvenile individuals originating from aquaculture conditions were translocated to group of isolated water bodies not inhabited by this species. The juveniles were offspring of parental fish belonging to the same local population, which is extinct at present. Five of those attempts were successful. The aim of the present study was to assess the genetic variation in a group new populations and compare genetic variation indicators with 13 old populations that had existed for decades. The polymorphism of 13 microsatellite markers was investigated, significance of differences in the genetic variation indicators between the groups were tested using a one-way analysis of variance (ANOVA). The mean values of all summary statistics under study, i.e. observed heterozygosity, expected heterozygosity and the total number of alleles, were higher in the group of new populations compared to almost all old ones. A similar dependence was found for Garza-Williamson M values, where the mean for the group of new populations was higher than in almost all old populations. Our results indicate that all recently established E. percnurus populations have not yet experienced any extensive founder effects or bottlenecks. They have preserved a large part of the genetic variability typical of their maternal population, which might also have been relatively high. This feature of new populations, may give them a relatively high ability to adapt to changing environments in the future., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Wolnicki et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The Labeoninae subfamily is a highly diversified but demonstrably monophyletic lineage of cyprinid fishes comprising five tribes and six incertae sedis genera. This widely distributed assemblage contains some 48 genera and around 480 recognized species distributed in freshwaters of Africa and Asia. In this study, the karyotypes and other chromosomal properties of five Labeoninae species found in Thailand Labeo chrysophekadion (Labeonini) and Epalzeorhynchos bicolor, Epalzeorhynchos munense, Henicorhynchus siamensis, Thynnichthys thynnoides (´Osteochilini´) were examined using conventional and molecular cytogenetic protocols. Our results confirmed a diploid chromosome number (2n) invariably 2n = 50, but the ratio of uni- and bi-armed chromosomes was highly variable among their karyotypes, indicating extensive structural chromosomal rearrangements. Karyotype of L. chrysophekadion contained 10m+6sm+20st+14a, 32m+10sm+8st for H. siamensis, 20m+12sm+10st+8a in E. bicolor, 20m+8sm+8st+14a in E. munense, and 18m+24sm+8st in T. thynnoides. Except for H. siamensis, which had four sites of 5S rDNA sites, other species under study had only one chromosome pair with those sites. In contrast, only one pair containing 18S rDNA sites were found in the karyotypes of three species, whereas two sites were found in that of E. bicolor. These cytogenetic patterns indicated that the cytogenomic divergence patterns of these labeonine species largely corresponded to the inferred phylogenetic tree. In spite of the 2n stability, diverse patterns of rDNA and microsatellite distribution as well as their various karyotype structures demonstrated significant evolutionary differentiation of Labeoninae genomes as exemplified in examined species. Labeoninae offers a traditional point of view on the evolutionary forces fostering biological diversity, and the recent findings add new pieces to comprehend the function of structural chromosomal rearrangements in adaption and speciation., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Khensuwan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The mitochondrial genome plays a crucial role in the balance of energy and heat production in organisms and, thus, may be under natural selection due to its potential role in adaptive divergence and speciation. Here, we studied natural selection on the mitogenome of Acrossocheilus (Cypriniformes: Cyprinidae), a genus of fish that inhabits a broad latitudinal distribution ranging from the tropics and subtropics through temperate regions. Specifically, we used 25 published mitogenome sequences of Acrossocheilus species to investigate phylogenetic relationships in this genus and detected signals of positive selection on 13 protein-coding, mitochondrial genes. We found that relaxed purifying selection and genetic drift were the predominant evolutionary forces acting on the analyzed mitogenomes. However, we also found evidence of diversifying selection on some codons, indicating episodes of positive selection. Additionally, we analyzed the mitogenomic data within an environmental modeling framework and found that the Ka/Ks ratio of ATP6 may correlated with a mean diurnal temperature range (p = 0.0449), while the Ka/Ks ratio of COX2 may correlated with precipitation during the driest month (p = 0.00761). These results suggest that the mitogenomes of Acrossocheilus species may be involved in evolutionary adaptations to different habitats. Based on this, we believe that our study provides a new insight into the role of the mitochondrial genome of Acrossocheilus species in adaptation to different environments. During our study, we also discovered several cases of paraphyly and polyphyly among accessions of species and their putative synonyms. Thus, our study suggests that a careful reassessment of the taxonomy of Acrossocheilus is using high-quality molecular data merited., Competing Interests: The authors have declared that no competing interests exist.

Vertebrate pigmentation is a fundamentally important, multifaceted phenotype. Zebrafish, Danio rerio, has been a valuable model for understanding genetics and development of pigment pattern formation due to its genetic and experimental tractability, advantages that are shared across several Danio species having a striking array of pigment patterns. Here, we use the sister species D. quagga and D. kyathit, with stripes and spots, respectively, to understand how natural genetic variation impacts phenotypes at cellular and organismal levels. We first show that D. quagga and D. kyathit phenotypes resemble those of wild-type D. rerio and several single locus mutants of D. rerio, respectively, in a morphospace defined by pattern variation along dorsoventral and anteroposterior axes. We then identify differences in patterning at the cellular level between D. quagga and D. kyathit by repeated daily imaging during pattern development and quantitative comparisons of adult phenotypes, revealing that patterns are similar initially but diverge ontogenetically. To assess the genetic architecture of these differences, we employ reduced-representation sequencing of second-generation hybrids. Despite the similarity of D. quagga to D. rerio, and D. kyathit to some D. rerio mutants, our analyses reveal a complex genetic basis for differences between D. quagga and D. kyathit, with several quantitative trait loci contributing to variation in overall pattern and cellular phenotypes, epistatic interactions between loci, and abundant segregating variation within species. Our findings provide a window into the evolutionary genetics of pattern-forming mechanisms in Danio and highlight the complexity of differences that can arise even between sister species. Further studies of natural genetic diversity underlying pattern variation in D. quagga and D. kyathit should provide insights complementary to those from zebrafish mutant phenotypes and more distant species comparisons., Competing Interests: The authors have declared that no competing interests exist.

Loss in intraspecific diversity can alter ecosystem functions, but the underlying mechanisms are still elusive, and intraspecific biodiversity-ecosystem function (iBEF) relationships have been restrained to primary producers. Here, we manipulated genetic and functional richness of a fish consumer (Phoxinus phoxinus) to test whether iBEF relationships exist in consumer species and whether they are more likely sustained by genetic or functional richness. We found that both genotypic and functional richness affected ecosystem functioning, either independently or interactively. Loss in genotypic richness reduced benthic invertebrate diversity consistently across functional richness treatments, whereas it reduced zooplankton diversity only when functional richness was high. Finally, losses in genotypic and functional richness altered functions (decomposition) through trophic cascades. We concluded that iBEF relationships lead to substantial top-down effects on entire food chains. The loss of genotypic richness impacted ecological properties as much as the loss of functional richness, probably because it sustains "cryptic" functional diversity., Competing Interests: The authors have declared that no competing interest exist.

We have gained considerable insight into the mechanisms which recognize and repair DNA damage, but how they adapt to extreme environmental challenges remains poorly understood. Cavefish have proven to be fascinating models for exploring the evolution of DNA repair in the complete absence of UV-induced DNA damage and light. We have previously revealed that the Somalian cavefish Phreatichthys andruzzii, lacks photoreactivation repair via the loss of light, UV and ROS-induced photolyase gene transcription mediated by D-box enhancer elements. Here, we explore whether other systems repairing UV-induced DNA damage have been similarly affected in this cavefish model. By performing a comparative study using P. andruzzii and the surface-dwelling zebrafish, we provide evidence for a conservation of sunlight-regulated Nucleotide Excision Repair (NER). Specifically, the expression of the ddb2 gene which encodes a key NER recognition factor is robustly induced following exposure to light, UV and oxidative stress in both species. As in the case of the photolyase genes, D-boxes in the ddb2 promoter are sufficient to induce transcription in zebrafish. Interestingly, despite the loss of D-box-regulated photolyase gene expression in P. andruzzii, the D-box is required for ddb2 induction by visible light and oxidative stress in cavefish. However, in the cavefish ddb2 gene this D-box-mediated induction requires cooperation with an adjacent, highly conserved E2F element. Furthermore, while in zebrafish UV-induced ddb2 expression results from transcriptional activation accompanied by stabilization of the ddb2 mRNA, in P. andruzzii UV induces ddb2 expression exclusively via an increase in mRNA stability. Thus, we reveal plasticity in the transcriptional and post transcriptional mechanisms regulating the repair of sunlight-induced DNA damage under long-term environmental challenges., Competing Interests: The authors have declared that no competing interests exist.

The phenomenon of a massive vertebral deformity was recorded in the radiating Labeobarbus assemblage from the middle reaches of the Genale River (south-eastern Ethiopia, East Africa). Within this sympatric assemblage, five trophic morphs-generalized, lipped, piscivorous and two scraping feeders-were reported between 1993 and 2019. In 2009, a new morph with prevalence of ~10% was discovered. The new morph, termed 'short', had an abnormally shortened vertebral column and a significantly deeper body. This type of deformity is common in farmed Atlantic salmon and other artificially reared fish, but is rare in nature. In the Genale Labeobarbus assemblage, the deformity was present exclusively within the generalized and lipped morphs. The short morph had between seven and 36 deformed (compressed and/or fused) vertebrae. Their body depth was positively correlated with number of deformed vertebrae. In another collection in 2019, the short morph was still present at a frequency of 11%. Various environmental and genetic factors could contribute to the development of this deformity in the Genale Labeobarbus, but based on the available data, it is impossible to confidently identify the key factor(s). Whether the result of genetics, the environment, or both, this deep-bodied phenotype is assumed to be an anti-predator adaptation, as there is evidence of its selective advantage in the generalized morph. The Genale monstrosity is the first reported case of a massive deformity of the vertebral column in a natural population of African fishes., Competing Interests: The authors have declared that no competing interests exist.

In many aquatic species, alteration of habitats and human-induced barriers shape the population's genetic structure in rivers with longitudinal connectivity. The golden mahseer, Tor putitora (GM) is an endangered and sensitive cyprinid species. It is considered an indicator of a healthy freshwater ecosystem and is found in cold-water habitats. Therefore, it is crucial to understand how longitudinal connectivity and anthropogenic factors affect the diversity and population genetic structure of GM. The population genetic structure, gene flow and demography of the GM in four Himalayan rivers were investigated by mitochondrial cytochrome b gene (cyt b) as well as microsatellite genotyping. The results showed overall high mtDNA diversity (hd: 0.795) couple with low nucleotide diversity (π: 0.0012) in all GM populations. We also found significant levels of observed heterozygosity (ranging from 0.618 to 0.676), with three genetic clusters. The mtDNA and microsatellite analysis suggested that there are close genetic relationships between the Bhagirathi and Ganga populations; whereas, significant level of genetic differentiation was observed with that of Alaknanda and Yamuna populations. Haplotype distribution, unimodal distribution graph and results of the neutrality test indicated a sign of recent population growth in the GM population. Analysis of molecular variance (AMOVA) and spatial molecular variance (SAMOVA) revealed existence of genetic structures in GM populations. In addition, spatial genetic analysis detected a significant correlation between the pairwise genetic and geographical distances for the entire study area (Mantel test, rM = 0.126; P = 0.010). Considering the significant level of heterozygosity, high rate of unidirectional migration and the intra-population structuring in Alaknanda and Yamuna, it is crucial to propose an effective conservation plan for the GM populations. In general, dams obstruct continuous water flow and create isolated microhabitats. Therefore, we recommend the establishment of microscale protected areas near GM breeding sites and construction of fish pass to maintain the genetic connectivity of fish species that enhance viable populations., Competing Interests: The authors have declared that no competing interests exist.

Vision is very important to fish as it is required for foraging food, fighting competitors, fleeing from predators, and finding potential mates. Vertebrates express opsin genes in photoreceptor cells to receive visual signals, and the variety of light levels in aquatic habits has driven fish to evolve multiple opsin genes with expression profiles that are highly plastic. In this study, red shiners (Cyprinella lutrensis) were exposed to four water turbidity treatments and their opsin genes were cloned to elucidate how opsin gene expression could be modulated by ambient light conditions. Opsin gene cloning revealed that these fish have single RH1, SWS1, SWS2 and LWS genes and two RH2 genes. Phylogenetic analysis also indicated that these two RH2 opsin genes-RH2A and RH2B -are in-paralogous. Using quantitative PCR, we found evidence that opsin expression is plastic in adults. Elevated proportional expression of LWS in the cone under ambient light and turbid treatment indicated that the red shiner's visual spectrum displays a red shift in response to increased turbidity., Competing Interests: The authors have declared that no competing interests exist.

Starvation is a common stress in fish. The underlying molecular mechanisms associated with growth depression caused by feeding restriction and compensatory growth are not well understood. We investigated the effect of fasting and refeeding on the transcriptome profiles of brain in juvenile S. hollandi using RNA-seq. A total of 4.73 × 108 raw reads were obtained from nine brain samples. De novo transcriptome assembly identified 387,085 unigenes with 2.1×109 nucleotides. A total of 936 annotated unigenes showed significantly differential expression among the control, fasting, and fasting-refeeding groups. The down-regulated differentially expressed genes (DEGs) during fasting were mainly associated with cell cycle, DNA replication, and mitosis. The up-regulated DEGs were mainly related to glucose and lipid metabolism, material transportation, and transcription factors. Most decreased DEGs during fasting were restored to normal levels after refeeding. Comparing with the control group, genes associated with protein synthesis, stimulus response, and carbohydrate metabolism were significantly over-expressed and pro-opio melanocortin (POMC) was down-regulated during the refeeding period. In conclusion, fish mobilized stored energetic materials and reduced energy consumption to prolong survival during fasting. After refeeding, the down-regulation of DEGs, e.g., POMC may be associated with compensatory growth. Up-regulation of DEGs related to ribosomal protein, stimulus response, and carbohydrate metabolism may contribute to eliminate negative effect of starvation on brain. This study provided the first transcriptome data related with impact of short-time starvation and refeeding in S. hollandi brains., Competing Interests: The authors have declared that no competing interests exist.

Tench (Tinca tinca L.) has great economic potential due to its high rate of fecundity and long-life span. Population genetic studies based on allozymes, microsatellites, PCR-RFLP and sequence analysis of genes and DNA fragments have revealed the presence of Eastern and Western phylogroups. However, the lack of genomic resources for this species has complicated the development of genetic markers. In this study, the tench transcriptome and genome were sequenced by high-throughput sequencing. A total of 60,414 putative SNPs were identified in the tench transcriptome using a computational pipeline. A set of 96 SNPs was selected for validation and a total of 92 SNPs was validated, resulting in the highest conversion and validation rate for a non-model species obtained to date (95.83%). The validated SNPs were used to genotype 140 individuals belonging to two tench breeds (Tabor and Hungarian), showing low (FST = 0.0450) but significant (<0.0001) genetic differentiation between the two tench breeds. This implies that set of validated SNPs array can be used to distinguish the tench breeds and that it might be useful for studying a range of associations between DNA sequence and traits of importance. These genomic resources created for the tench will provide insight into population genetics, conservation fish stock management, and aquaculture., Competing Interests: The authors have declared that no competing interests exist.

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., Competing Interests: The authors have declared that no competing interests exist.

The Sturgeon chub (Macrhybopsis gelida) is a cyprinid fish native to the Missouri and Mississippi River basins of the U.S. Suspected long-term declines in the size of its distribution have prompted a review of its conservation status by the U.S. Fish and Wildlife Service, a process which depends on reliable methods to delineate the distribution and status of extant populations. To facilitate monitoring of Sturgeon chub populations, we developed a quantitative PCR assay to detect Sturgeon chub DNA in environmental samples. The assay consistently detected Sturgeon chub DNA in concentrations as low as 2 copies per reaction, and did not amplify DNA from non-target fish species that are sympatric in the upper Missouri River basin. Field tests of this assay with environmental samples successfully detected Sturgeon chub from sites known to be occupied. This assay offers an extremely sensitive methodology that can be applied to determine the range of Sturgeon chub, regardless of variation in habitat characteristics., Competing Interests: The authors have declared that no competing interests exist.

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., Competing Interests: The authors have declared that no competing interests exist.

Biological invasion by non-native subspecies or populations is one of the most serious threats to ecosystems, because these species might be easily established in the introduced area and can negatively affect native populations through competition and hybridization. Pale chub Opsariichthys platypus, one of the most common fish in East Asia, exhibits clear genetic differentiation among regional populations; however, introgression and subsequent loss of genetic integrity have been occurring throughout Japan due to the artificial introduction of non-native conspecifics. In this study, we developed a simple SNP genotyping method to discriminate between native and non-native mitochondrial DNA (mtDNA) haplotypes in pale chub using real-time PCR assay. We then investigated the distribution patterns of non-native pale chub in Tokai region, located in the center of Honshu Island, Japan and developed a predictive model of the occurrence of non-natives to reveal the factors influencing their invasion. The specificity and accuracy of the genotyping method were confirmed by using samples whose haplotypes were determined previously. Extensive occurrence of non-native haplotypes in Tokai region was detected by this method. In addition, our models suggested that the presence of non-natives varied greatly depending on the river system, and was positively influenced by the impounded water areas. Our method could accurately distinguish between native and non-native haplotypes of pale chub in Japan and suggested key environmental factors associated with the presence of non-natives. This approach can greatly reduce experimental costs be a great contribution for quantitative investigation.

Four species of Devario are recorded from Bangladesh: D. aequipinnatus, D. anomalus, D. coxi, new species, and D. devario. Devario aequipinnatus has a wide distribution in northern India and Bangladesh. Devario coxi, from southeastern Bangladesh near Cox's Bazar, differs from D. aequipinnatus in mtDNA (COI, p-distance 1.8%), colouration, proportional measurements, and meristics. The minor morphological differences and low frequency of overlapping meristics suggest relatively recent separation of D. coxi from other D. aequipinnatus. Devario anomalus occurs only in southeastern Bangladesh and is here reported from localities in addition to the type locality. It differs from the similar D. xyrops in adjacent Myanmar by slender body shape and by 2.3% p-distance in the COI gene. Specimens of D. anomalus from the Sangu River were found to have the mitochondrial genome of D. aequipinnatus from Bangladesh, but agree with other D. anomalus in the nuclear RAG1 gene. Devario devario has a wide distribution on the Indian Peninsula and border regions; in Bangladesh it is restricted in distribution to the Ganga, Brahmaputra, and Meghna drainages. Reports of D. assamensis and D. malabaricus from Bangladesh are misidentifications. Perilampus ostreographus M'Clelland, 1839, is tentatively synonymized with D. aequipinnatus. Phylogenetic analysis of 14 species of striped devarios based on the COI gene results in a polytomy with four unresolved clades. Devario deruptotalea from the Chindwin basin is the sister group of D. aequipinnatus+D. coxi. Devario devario is the sistergroup of D. xyrops+D. anomalus.

Current knowledge on the biological responses of freshwater fish under projected scenarios of climate change remains limited. Here, we examine differences in the protein configuration of two endemic Iberian freshwater fish species, Squalius carolitertii and the critically endangered S. torgalensis that inhabit in the Atlantic-type northern and in the Mediterranean-type southwestern regions, respectively. We performed protein structure modeling of fourteen genes linked to protein folding, energy metabolism, circadian rhythms and immune responses. Structural differences in proteins between the two species were found for HSC70, FKBP52, HIF1α and GPB1. For S. torgalensis, besides structural differences, we found higher thermostability for two proteins (HSP90 and GBP1), which can be advantageous in a warmer environment. Additionally, we investigated how these species might respond to projected scenarios of 3° climate change warming, acidification (ΔpH = -0.4), and their combined effects. Significant changes in gene expression were observed in response to all treatments, particularly under the combined warming and acidification. While S. carolitertii presented changes in gene expression for multiple proteins related to folding (hsp90aa1, hsc70, fkbp4 and stip1), only one such gene was altered in S. torgalensis (stip1). However, S. torgalensis showed a greater capacity for energy production under both the acidification and combined scenarios by increasing cs gene expression and maintaining ldha gene expression in muscle. Overall, these findings suggest that S. torgalensis is better prepared to cope with projected climate change. Worryingly, under the simulated scenarios, disturbances to circadian rhythm and immune system genes (cry1aa, per1a and gbp1) raise concerns for the persistence of both species, highlighting the need to consider multi-stressor effects when evaluating climate change impacts upon fish. This work also highlights that assessments of the potential of endangered freshwater species to cope with environmental change are crucial to help decision-makers adopt future conservation strategies.

The Chinese lake gudgeon Sarcocheilichthys sinensis is a small cyprinid fish with great aquaculture potential both for its edible and ornamental values. Nevertheless, available genomic and transcriptomic information for this fish is extremely deficient. In this study, a normalized cDNA library was constructed using 13 mixed tissues of an adult male S. sinensis, and was sequenced by the Illumina HiSeq2500 platform. De novo assembly was performed using 38,911,511 obtained clean reads, and a total of 147,282 unigenes with an average length of 900 bp were finally achieved. 96.2% of these unigenes were annotated in 9 public databases, and 16 segments of growth-related genes were identified for future studies. In addition, 28,493 unigenes were assigned to 61 subcategories of Gene Ontology (GO), and 10,483 unigenes were assigned to 25 categories of Cluster of Orthologous Group (COG). Moreover, 14,943 unigenes were classified into 225 pathways of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 30,666 microsatellites were detected from 17,627 unigenes with an average distribution density of 1:2405 bp. This transcriptome data set will be valuable for researches on discovery, expression and evolution on genes of interest. Meanwhile, the identified microsatellites would be useful tools for genetic and genomic studies in S. sinensis., Competing Interests: The authors have declared that no competing interests exist.

The lake minnow (Eupallasella percnurus) is critically endangered. In this paper we characterize the genetic properties of this fish over its range of occurrence in Poland and propose the use of this knowledge in its active protection. Twelve populations of lake minnow from across its range in Poland were investigated. 13 microsatellite loci were investigated to evaluate genetic variation and distance among populations. The magnitude of the genetic bottleneck or founder effects was investigated. In the studied populations, the allelic diversity and heterozygosity showed that genetic variation in this species is low. At most loci, only 2-3 alleles per population were detected. The average number of alleles detected across all loci was 35, and ranged from 24 to 53. The average observed heterozygosity (Ho) across all investigated loci was 0.38 (range 0.21-0.59); the average expected heterozygosity (He) was 0.36 (range 0.18-0.55). The populations remained in Hardy-Weinberg equilibrium. The average Garza-Williamson M index value for all populations was low (0.47), suggesting a reduction in genetic variation due to a founder effect or a genetic bottleneck. Genetic distance among populations was high or very high (FST range: 0.20-0.64; δμ2 range: 1.32-16.98); this was likely a consequence of low gene flow among isolated populations, a founder effect or other genetic bottleneck, and strong genetic drift. The large genetic differences among the investigated lake minnow populations are likely to also exist among other populations of this species, and knowledge of these differences should inform active protection programs based on translocation of wild or cultivated fish of this species. The method presented here can potentially be applied to any population of lake minnows or closely related species., Competing Interests: The authors have declared that no competing interests exist.

Percocypris pingi is a near threatened cyprinid species, which has suffered a dramatic decline due to anthropogenic factors. As one response to this decline, hatchery release for P. pingi has been conducted in the lower reaches of the Yalong River since 2012. To understand the conservation status of this species and the potential impact of the release of hatchery-reared fish, we studied the genetic diversity and population structure of wild and hatchery populations of P. pingi. Two hatchery populations (Jinping [JPH] and Ya'an [YAH]) and two wild populations (Muli [MLW] and Woluo [WLW]) of P. pingi were analyzed based on microsatellite markers and the mitochondrial DNA control region. The results showed that P. pingi possesses moderate levels of genetic diversity, with observed heterozygosities ranging from 0.657 to 0.770 and nucleotide diversities ranging from 0.00212 to 0.00491. Our results also suggested WLW harbors considerable proportions of genetic diversity in this species and serves as a refuge for P. pingi during anthropogenic disturbance, thus playing an important role for the conservation of P. pingi populations. Microsatellite and mitochondrial markers both indicated close genetic relationships between YAH and MLW, JPH and WLW, respectively. The results to some extent reflected the geographical provenances for original broodstocks of the two hatchery populations, which provide some practical guidance for hatchery release of P. pingi. The existence of remarkable genetic divergence distributed along limited geographical range (approximately 10 kilometers) suggests the two wild populations should be regarded at least as two distinct evolutionary significant units (ESUs) and management units (MUs). Considering reduced intra-population genetic variation in hatchery population for release and significant genetic compositions of the two hatchery populations, some appropriate breeding strategies were proposed to benefit conservation of P. pingi., Competing Interests: This study was funded in part by Yalong River Hydropower Development Company, Ltd (http://www.ehdc.com.cn/faces/en/overview2.html). This company is an enterprise engaged in development of hydropower projects and generation of hydroelectric energy. At the same time, it has been doing environment protection work associated with the hydropower development. Weixiong Gan, Rukui Zeng and Longjun Deng are employed by Yalong River Hydropower Development Company, Ltd. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.

Hybridization and incomplete lineage sorting are common confounding factors in phylogeny and speciation resulting in mitonuclear disparity. Mitochondrial introgression, a particular case of hybridization, may, in extreme cases, lead to replacement of the mitochondrial genome of one species with that of another (mitochondrial capture). We investigated mitochondrial introgression involving two species of the cyprinid genus Squalius in the western Peloponnese region of Greece using molecular and morphological data. We found evidence of complete mitochondrial introgression of Squalius keadicus into two populations recognized as Squalius peloponensis from the Miras and Pamissos River basins and a divergence of mitochondrial genomes of S. keadicus from the Evrotas basin from that of the introgressed populations dating from the Pleistocene. Secondary contact among basins is a possible factor in connection of the species and the introgression event. Morphological analyses support the hypothesis of mitochondrial introgression, as S. keadicus was different from the other three populations recognized as S. peloponensis, although significant differences were found among the four populations. Isolation by geographical barriers arose during Pleistocene in the western Peloponnese were the source of the evolution of the two reciprocally monophyletic subclades found in the S. keadicus mitochondrial clade, and the morphological differences found among the four populations. Along with the lack of structure in the nuclear genome in the three populations ascribed to S. peloponensis, this suggests an incipient speciation process occurring in these Squalius species in the western Peloponnese., Competing Interests: The authors have declared that no competing interests exist.

Loach minnow (Rhinichthys cobitis) and spikedace (Meda fulgida) are legally protected with the status of Endangered under the U.S. Endangered Species Act and are endemic to the Gila River basin of Arizona and New Mexico. Efficient and sensitive methods for monitoring these species' distributions are critical for prioritizing conservation efforts. We developed quantitative PCR assays for detecting loach minnow and spikedace DNA in environmental samples. Each assay reliably detected low concentrations of target DNA without detection of non-target species, including other cyprinid fishes with which they co-occur., Competing Interests: The authors have declared that no competing interests exist.

Capoeta damascina was earlier considered by many authors as one of the most common freshwater fish species found throughout the Levant, Mesopotamia, Turkey, and Iran. However, owing to a high variation in morphological characters among and within its various populations, 17 nominal species were described, several of which were regarded as valid by subsequent revising authors. Capoeta damascina proved to be a complex of closely related species, which had been poorly studied. The current study aims at defining C. damascina and the C. damascina species complex. It investigates phylogenetic relationships among the various members of the C. damascina complex, based on mitochondrial and nuclear DNA sequences. Phylogenetic relationships were projected against paleogeographical events to interpret the geographic distribution of the taxa under consideration in relation to the area's geological history. Samples were obtained from throughout the geographic range and were subjected to genetic analyses, using two molecular markers targeting the mitochondrial cytochrome oxidase I (n = 103) and the two adjacent divergence regions (D1-D2) of the nuclear 28S rRNA genes (n = 65). Six closely related species were recognized within the C. damascina complex, constituting two main lineages: A western lineage represented by C. caelestis, C. damascina, and C. umbla and an eastern lineage represented by C. buhsei, C. coadi, and C. saadii. The results indicate that speciation of these taxa is rather a recent event. Dispersal occurred during the Pleistocene, resulting in present-day distribution patterns. A coherent picture of the phylogenetic relationships and evolutionary history of the C. damascina species complex is drawn, explaining the current patterns of distribution as a result of paleogeographic events and ecological adaptations.

Schizothorax prenanti (S. prenanti) is mainly distributed in the upstream regions of the Yangtze River and its tributaries in China. This species is indigenous and commercially important. However, in recent years, wild populations and aquacultures have faced the serious challenges of germplasm variation loss and an increased susceptibility to a range of pathogens. Currently, the genetics and immune mechanisms of S. prenanti are unknown, partly due to a lack of genome and transcriptome information. Here, we sought to identify genes related to immune functions and to identify molecular markers to study the function of these genes and for trait mapping. To this end, the transcriptome from spleen tissues of S. prenanti was analyzed and sequenced. Using paired-end reads from the Illumina Hiseq2500 platform, 48,517 transcripts were isolated from the spleen transcriptome. These transcripts could be clustered into 37,785 unigenes with an N50 length of 2,539 bp. The majority of the unigenes (35,653, 94.4%) were successfully annotated using non-redundant nucleotide sequence analysis (nt), and the non-redundant protein (nr), Swiss-Prot, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. KEGG pathway assignment identified more than 500 immune-related genes. Furthermore, 7,545 putative simple sequence repeats (SSRs), 857,535 single nucleotide polymorphisms (SNPs), and 53,481 insertion/deletion (InDels) were detected from the transcriptome. This is the first reported high-throughput transcriptome analysis of S. prenanti, and it provides valuable genetic resources for the investigation of immune mechanisms, conservation of germplasm, and molecular marker-assisted breeding of S. prenanti.

The blunt snout bream (Megalobrama amblycephala) is an important freshwater aquaculture species, but it is sensitive to hypoxia. No transcriptome data related to growth and hypoxia response are available for this species. In this study, we performed de novo transcriptome sequencing for the liver and gills of the fast-growth family and slow-growth family derived from 'Pujiang No.1' F10 blunt snout bream that were under hypoxic stress and normoxia, respectively. The fish were divided into the following 4 groups: fast-growth family under hypoxic stress, FH; slow-growth family under hypoxic stress, SH; fast-growth family under normoxia, FN; and slow-growth family under normoxia, SN. A total of 185 million high-quality reads were obtained from the normalized cDNA of the pooled samples, which were assembled into 465,582 contigs and 237,172 transcripts. A total of 31,338 transcripts from the same locus (unigenes) were annotated and assigned to 104 functional groups, and 23,103 unigenes were classified into seven main categories, including 45 secondary KEGG pathways. A total of 22,255 (71%) known putative unigenes were found to be shared across the genomes of five model fish species and mammals, and a substantial number (9.4%) of potentially novel genes were identified. When 6,639 unigenes were used in the analysis of differential expression (DE) genes, the number of putative DE genes related to growth pathways in FH, SH, SN and FN was 159, 118, 92 and 65 in both the liver and gills, respectively, and the number of DE genes related to hypoxic response was 57, 33, 23 and 21 in FH, FN, SH and SN, respectively. Our results suggest that growth performance of the fast-growth family should be due to complex mutual gene regulatory mechanisms of these putative DE genes between growth and hypoxia.

Ten microsatellite loci were characterized for 34 locations from roundtail chub (Gila robusta complex) to better resolve patterns of genetic variation among local populations in the lower Colorado River basin. This group has had a complex taxonomic history and previous molecular analyses failed to identify species diagnostic molecular markers. Our results supported previous molecular studies based on allozymes and DNA sequences, which found that most genetic variance was explained by differences among local populations. Samples from most localities were so divergent species-level diagnostic markers were not found. Some geographic samples were discordant with current taxonomy due to admixture or misidentification; therefore, additional morphological studies are necessary. Differences in spatial genetic structure were consistent with differences in connectivity of stream habitats, with the typically mainstem species, G. robusta, exhibiting greater genetic connectedness within the Gila River drainage. No species exhibited strong isolation by distance over the entire stream network, but the two species typically found in headwaters, G. nigra and G. intermedia, exhibited greater than expected genetic similarity between geographically proximate populations, and usually clustered with individuals from the same geographic location and/or sub-basin. These results highlight the significance of microevolutionary processes and importance of maintaining local populations to maximize evolutionary potential for this complex. Augmentation stocking as a conservation management strategy should only occur under extreme circumstances, and potential source populations should be geographically proximate stocks of the same species, especially for the headwater forms.

Changing drainage patterns have played a significant role in the evolution of western North American aquatic taxa. Relict dace, Relictus solitarius, is a Great Basin endemic cyprinid with a native range that is restricted to four valleys in eastern Nevada. Relictus solitarius now occupies spring systems that are the remnants of Pleistocene-era pluvial lakes, although it may have occurred in the area for much longer. Here we use mitochondrial DNA sequence data to assess range-wide genetic diversity of R. solitarius, and to estimate divergence times to determine whether pluvial drainages played an important role in shaping intraspecific genetic diversity. Genetic diversification within R. solitarius began during the early to mid-Pleistocene, separating populations within two sets of valleys (Butte/Ruby and Goshute/Steptoe). Additional diversification in each of the two sets of valleys occurred more recently, in the mid- to late-Pleistocene. Holocene desiccation has further isolated populations, and each population sampled contains unique mtDNA haplotypes. Pluvial drainage patterns did contribute to the genetic structure observed within R. solitarius, but most of the intraspecific diversification does not appear to be associated with the Last Glacial Maximum. Holocene desiccation has also contributed to the observed genetic structure. The relict dace populations we sampled are all unique, and we recommend that future management efforts should strive to preserve as much of the genetic diversity as possible.

Previous studies have shown that assortative mating acts as a driver of speciation by countering hybridization between two populations of the same species (pre-zygotic isolation) or through mate choice among the hybrids (hybrid speciation). In both speciation types, assortative mating promotes speciation over a transient hybridization stage. We studied mate choice in a hybrid vertebrate complex, the allopolyploid fish Squalius alburnoides. This complex is composed by several genomotypes connected by an intricate reproductive dynamics. We developed a model that predicts the hybrid complex can persist when females exhibit particular mate choice patterns. Our model is able to reproduce the diversity of population dynamic outcomes found in nature, namely the dominance of the triploids and the dominance of the tetraploids, depending on female mate choice patterns and frequency of the parental species. Experimental mate choice trials showed that females exhibit the preferences predicted by the model. Thus, despite the known role of assortative mating in driving speciation, our findings suggest that certain mate choice patterns can instead hinder speciation and support the persistence of hybrids over time without speciation or extinction.

How allopolyploids are able not only to cope but profit from their condition is a question that remains elusive, but is of great importance within the context of successful allopolyploid evolution. One outstanding example of successful allopolyploidy is the endemic Iberian cyprinid Squalius alburnoides. Previously, based on the evaluation of a few genes, it was reported that the transcription levels between diploid and triploid S. alburnoides were similar. If this phenomenon occurs on a full genomic scale, a wide functional ''diploidization'' could be related to the success of these polyploids. We generated RNA-seq data from whole juvenile fish and from adult livers, to perform the first comparative quantitative transcriptomic analysis between diploid and triploid individuals of a vertebrate allopolyploid. Together with an assay to estimate relative expression per cell, it was possible to infer the relative sizes of transcriptomes. This showed that diploid and triploid S. alburnoides hybrids have similar liver transcriptome sizes. This in turn made it valid to directly compare the S. alburnoides RNA-seq transcript data sets and obtain a profile of dosage responses across the S. alburnoides transcriptome. We found that 64% of transcripts in juveniles' samples and 44% in liver samples differed less than twofold between diploid and triploid hybrids (similar expression). Yet, respectively 29% and 15% of transcripts presented accurate dosage compensation (PAA/PA expression ratio of 1 instead of 1.5). Therefore, an exact functional diploidization of the triploid genome does not occur, but a significant down regulation of gene expression in triploids was observed. However, for those genes with similar expression levels between diploids and triploids, expression is not globally strictly proportional to gene dosage nor is it set to a perfect diploid level. This quantitative expression flexibility may be a strong contributor to overcome the genomic shock, and be an immediate evolutionary advantage of allopolyploids.

Fathead minnow and zebrafish are among the most intensively studied fish species in environmental toxicogenomics. To aid the assessment and interpretation of subtle transcriptomic effects from treatment conditions of interest, better characterization and understanding are needed for natural variation in gene expression among fish individuals from lab cultures. Leveraging the transcriptomics data from a number of our toxicogenomics studies conducted over the years, we conducted a meta-analysis of nearly 600 microarrays generated from the ovary tissue of untreated, reproductively mature fathead minnow and zebrafish samples. As expected, there was considerable batch-to-batch transcriptomic variation; this "batch-effect" appeared to differentially impact subsets of fish transcriptomes in a nonsystematic way. Temporally more closely spaced batches tended to share a greater transcriptomic similarity among one another. The overall level of within-batch variation was quite low in fish ovary tissue, making it a suitable system for studying chemical stressors with subtle biological effects. The observed differences in the within-batch variability of gene expression, at the levels of both individual genes and pathways, were probably both technical and biological. This suggests that biological interpretation and prioritization of genes and pathways targeted by experimental conditions should take into account both their intrinsic variability and the size of induced transcriptional changes. There was significant conservation of both the genomes and transcriptomes between fathead minnow and zebrafish. The high degree of conservation offers promising opportunities in not only studying fish molecular responses to environmental stressors by a comparative biology approach, but also effective sharing of a large amount of existing public transcriptomics data for developing toxicogenomics applications.

The blunt snout bream (Megalobrama amblycephala) is an important freshwater aquaculture fish throughout China. Because of widespread introductions of this species to many regions, the genetic diversity of wild and natural populations is now threatened. In the present study, SRAP (sequence-related amplified polymorphism) markers were used to assess genetic diversity of blunt snout bream. Three natural populations (Liangzi Lake, Poyang Lake and Yuni Lake, one cultured population (Nanxian) and one genetic strain ('Pujiang No. 1') of blunt snout bream were screened with 88 SRAP primer combinations, of which 13 primer pairs produced stable and reproducible amplification patterns. In total, 172 bands were produced, of which 132 bands were polymorphic. Nei's gene diversity (h) and Shannon's information index (I) values provided evidence of differences in genetic diversity among the five populations (Poyang Lake>Liangzi Lake>Nanxian>'Pujiang No. 1'>Yuni Lake). Based on cluster analysis conducted on genetic distance values, the five blunt snout bream populations were divided into three groups, Poyang Lake and Liangzi Lake (natural populations), Nanxian and 'Pujiang No. 1' (cultured population and genetically selected strain), and Yuni Lake (natural population). Significant genetic differentiation was found among the five populations using analysis of molecular variance (AMOVA), with more genetic divergence existing among populations (55.49%), than within populations (44.51%). This molecular marker technique is a simple and efficient method to quantify genetic diversity within and among fish populations, and is employed here to help manage and conserve germplasm variability of blunt snout bream and to support the ongoing selective breeding programme for this fish.

Although stream-dwelling gudgeons (Cyprinidae, genus: Gobio) are widespread in Central Europe, the taxonomy of this group and the distribution of its species are still unexplored in detail. The aims of our study are to ascertain taxonomic composition and distribution of the former Gobio gobio superspecies in the inner area of the Carpathian Basin. Since the presence of cryptic species is suspected in this area, we examined the taxonomic and phylogenetic relationships of Central European Gobio taxa by sequencing the mitochondrial DNA control region (mtCR). Additionally, we characterized the genetic structure of 27 stream-dwelling gudgeon populations of this area by Amplified Fragment Length Polymorphism (AFLP). Results of mtCR analysis proved the presence of three species already known as G. obtusirostris (dominant in NW-Hungary), G. gobio (sporadic) and G. carpathicus (sporadic). Additionally, the analysis revealed the existence of one doubtful taxon, G. sp1 (dominant in NE-Hungary), and a new isolated haplogroup (dominant in SW-Hungary). Although Network analysis showed significant detachment among haplogroups, their genetic distances were quite small. Therefore Bayesian phylogenetic analysis showed weak nodal support for the branching pattern both for newly described haplotypes, and for the already accepted species. AFLP data showed distinct population structure and a clear pattern of isolation was revealed by distance of stocks. At the same time, level of separation was not affected by the altitudinal position of sites. Moreover we found three major clusters of populations which were separated according to hydrographic regions, and corresponded to the findings of mtCR analysis. Our results suggest the on-going speciation of gudgeons in the Carpathian Basin, however the separation of haplogroups seems to only be an intermediate phase. The discovered natural pattern seems to be only slightly influenced by anthropogenic impacts. Additionally our results put into question the suitability of the recently accepted within Gobio genus taxonomy.

Blunt snout bream (Megalobrama amblycephala) are susceptible to hepatic steatosis when maintained in modern intensive culture systems. The aim of this study was to investigate the potential roles of microRNAs (miRNAs) in diet-induced hepatic steatosis in this species. MiRNAs, small non-coding RNAs that regulate gene expression at the posttranscriptional level, are involved in diverse biological processes, including lipid metabolism. Deep sequencing of hepatic small RNA libraries from blunt snout bream fed normal-fat and high-fat diets identified 202 (193 known and 9 novel) miRNAs, of which 12 were differentially expressed between the normal-fat and high-fat diet groups. Quantitative stem-loop reverse transcriptase-polymerase chain reaction analyses confirmed the upregulation of miR-30c and miR-30e-3p and the downregulation of miR-145 and miR-15a-5p in high-fat diet-fed fish. Bioinformatics tools were used to predict the targets of these verified miRNAs and to explore potential downstream gene ontology biological process categories and Kyoto Encyclopedia of Genes and Genomes pathways. Six putative lipid metabolism-related target genes (fetuin-B, Cyp7a1, NADH dehydrogenase (ubiquinone) 1 beta subcomplex subunit 2, 3-oxoacid CoA transferase 1b, stearoyl-CoA desaturase, and fatty-acid synthase) were identified as having potential important roles in the development of diet-induced hepatic steatosis in blunt snout bream. The results presented here are a foundation for future studies of miRNA-controlled lipid metabolism regulatory networks in blunt snout bream.

As a result of specific adaptations and habitat preferences strongly rheophilic fish species may show high levels of endemism. Many temperate rheophilic fish species were subjected to a series of range contractions during the Pleistocene, and then successfully expanded during the Holocene, colonising previously abandoned areas. The Carpathian barbel (Barbus carpathicus Kotlík, Tsigenopoulos, Ráb et Berrebi 2002) occurs in the montane streams in three basins of the main Central European rivers in the northern part of the Carpathian range. We used genetic variation within 3 mitochondrial and 9 microsatellite loci to determine a pattern of postglacial expansion in B. carpathicus. We found that overall genetic variation within the species is relatively low. Estimate of time to the most recent common ancestor (tMRCA) of mitochondrial sequences falls within the Holocene. The highest levels of genetic variation found in upper reaches of the Tisa river in the Danube basin suggest that glacial refugia were located in the south-eastern part of the species range. Our data suggest that the species crossed different watersheds at least six times as three genetically distinct groups (probably established in different expansion episodes) were found in northern part of the species range. Clines of genetic variation were observed in both the Danube and Vistula basins, which probably resulted from subsequent bottlenecks while colonizing successive habitats (south eastern populations) or due to the admixture of genetically diverse individuals to a previously uniform population (Vistula basin). Therefore, B. carpathicus underwent both demographic breakdowns and expansions during the Holocene, showing its distribution and demography are sensitive to environmental change. Our findings are important in the light of the current human-induced habitats alterations.

Background: The evolution of the Yunnan Plateau's drainages network during the Pleistocene was dominated by the intense uplifts of the Qinghai-Tibetan Plateau. In the present study, we investigated the association between the evolutionary histories of three main drainage systems and the geographic patterns of genetic differentiation of Poropuntius huangchuchieni., Methodology/principal Findings: We sequenced the complete sequences of mitochondrial control region for 304 specimens and the sequences of Cytochrome b gene for 15 specimens of the species P. huangchuchieni and 5 specimens of Poropuntius opisthoptera. Phylogenetic analysis identified five major lineages, of which lineages MK-A and MK-B constrained to the Mekong River System, lineages RL and LX to the Red River System, and lineage SW to the Salween River System. The genetic distance and network analysis detected significant divergences among these lineages. Mismatch distribution analysis implied that the population of P. huangchuchieni underwent demographic stability and the lineage MK-B, sublineages MK-A1 and LX-1 underwent a recent population expansion. The divergence of the 5 major lineages was dated back to 0.73-1.57 MYA., Conclusions/significance: Our results suggest that P. opisthoptera was a paraphyletic group of P. huangchuchieni. The phylogenetic pattern of P. huangchuchieni was mostly associated with the drainage's structures and the geomorphological history of the Southwest Yunnan Plateau. Also the differentiation of the major lineages among the three drainages systems coincides with the Kunlun-Yellow River Movement (1.10-0.60 MYA). The genetic differentiation within river basins and recent demographical expansions that occurred in some lineages and sublineages are consistent with the palaeoclimatic oscillations during the Pleistocene. Additionally, our results also suggest that the populations of P. huangchuchieni had keep long term large effective population sizes and demographic stability in the recent evolutionary history, which may be responsible for the high genetic diversity and incomplete lineages sorting of Poropuntius huangchuchieni.

Local adaptation towards divergent ecological conditions often results in genetic differentiation and adaptive phenotypic divergence. To illuminate the ecological distinctiveness of the schizothoracine fish, we studied a Gymnocypris species complex consisting of three morphs distributed across four bodies of water (the Yellow River, Lake Qinghai, the Ganzi River and Lake Keluke) in the Northeast Tibetan Plateau. We used a combination of mitochondrial (16S rRNA and Cyt b) and nuclear (RAG-2) genetic sequences to investigate the phylogeography of these morphs based on a sample of 277 specimens. Analysis of gill rakers allowed for mapping of phenotypic trajectories along the phylogeny. The phylogenetic and morphological analyses showed that the three sparsely rakered morphs were present at two extremes of the phylogenetic tree: the Yellow River morphs were located at the basal phylogenetic split, and the Lake Keluke and Ganzi River morphs at the peak, with the densely rakered Lake Qinghai morphs located between these two extremes. Age estimation further indicated that the sparsely rakered morphs constituted the oldest and youngest lineages, whereas the densely rakered morph was assigned to an intermediate-age lineage. These results are most compatible with the process of evolutionary convergence or reversal. Disruptive natural selection due to divergent habitats and dietary preferences is likely the driving force behind the formation of new morphs, and the similarities between their phenotypes may be attributable to the similarities between their forms of niche tracking associated with food acquisition. This study provides the first genetic evidence for the occurrence of convergence or reversal in the schizothoracine fish of the Tibetan Plateau at small temporal scales.

This study used phylogenetic analyses of mitochondrial cytochrome b sequences to investigate genetic diversity within three broadly co-distributed freshwater fish genera (Galaxias, Pseudobarbus and Sandelia) to shed some light on the processes that promoted lineage diversification and shaped geographical distribution patterns. A total of 205 sequences of Galaxias, 177 sequences of Pseudobarbus and 98 sequences of Sandelia from 146 localities across nine river systems in the south-western Cape Floristic Region (South Africa) were used. The data were analysed using phylogenetic and haplotype network methods and divergence times for the clades retrieved were estimated using *BEAST. Nine extremely divergent (3.5-25.3%) lineages were found within Galaxias. Similarly, deep phylogeographic divergence was evident within Pseudobarbus, with four markedly distinct (3.8-10.0%) phylogroups identified. Sandelia had two deeply divergent (5.5-5.9%) lineages, but seven minor lineages with strong geographical congruence were also identified. The Miocene-Pliocene major sea-level transgression and the resultant isolation of populations in upland refugia appear to have driven widespread allopatric divergence within the three genera. Subsequent coalescence of rivers during the Pleistocene major sea-level regression as well as intermittent drainage connections during wet periods are proposed to have facilitated range expansion of lineages that currently occur across isolated river systems. The high degree of genetic differentiation recovered from the present and previous studies suggest that freshwater fish diversity within the south-western CFR may be vastly underestimated, and taxonomic revisions are required.

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.

Two cyprinid species, Parachondrostoma toxostoma, an endemic threatened species, and Chondrostoma nasus, an invasive species, live in sympatry in southern France and form two sympatric zones where the presence of intergeneric hybrids is reported. To estimate the potential threat to endemic species linked to the introduction of invasive species, we focused on the DAB genes (functional MHC IIB genes) because of their adaptive significance and role in parasite resistance. More specifically, we investigated (1) the variability of MHC IIB genes, (2) the selection pattern shaping MHC polymorphism, and (3) the extent to which trans-species evolution and intergeneric hybridization affect MHC polymorphism. In sympatric areas, the native species has more diversified MHC IIB genes when compared to the invasive species, probably resulting from the different origins and dispersal of both species. A similar level of MHC polymorphism was found at population level in both species, suggesting similar mechanisms generating MHC diversity. In contrast, a higher number of DAB-like alleles per specimen were found in invasive species. Invasive species tended to express the alleles of two DAB lineages, whilst native species tended to express the alleles of only the DAB3 lineage. Hybrids have a pattern of MHC expression intermediate between both species. Whilst positive selection acting on peptide binding sites (PBS) was demonstrated in both species, a slightly higher number of positively selected sites were identified in C. nasus, which could result from parasite-mediated selection. Bayesian clustering analysis revealed a similar pattern of structuring for the genetic variation when using microsatellites or the MHC approach. We confirmed the importance of trans-species evolution for MHC polymorphism. In addition, we demonstrated bidirectional gene flow for MHC IIB genes in sympatric areas. The positive significant correlation between MHC and microsatellites suggests that demographic factors may contribute to MHC variation on a short time scale.

Fierce predatory freshwater fishes, the species of Percocypris (Cyprinidae, Teleostei) inhabit large rivers or lakes, and have a specific distribution pattern. Only a single species or subspecies occurs in each large-scale drainage basin of the Southeastern Tibetan Plateau. In this study, the molecular phylogenetic relationships for all but one of the described subspecies/species of Percocypris were investigated based on three mitochondrial genes (16S; COI; Cyt b) and one nuclear marker (Rag2). The results of Maximum Likelihood and Bayesian Inference analyses show that Percocypris is a strongly supported monophyletic group and that it is the sister group of Schizothorax. Combined with analyses of morphological characters, our results suggest that Percocypris needs to be reclassified, and we propose that six species be recognized, with corresponding distributions in five main drainages (including one lake). In addition, based on the results of the estimation of divergence times and ancestral drainages, we hypothesize that Percocypris likely originated in the early Miocene from a paleo-connected drainage system containing the contemporary main drainages of the Southeastern Tibetan Plateau. This study suggests that vicariance (due to the uplift of the Tibetan Plateau modifying the large-scale morphologies of drainage basins in the Southeastern Tibetan Plateau) has played an important role in the speciation of the genus. Furthermore, external morphological characters (such as the length of the fins) and an internal trait (the position of pterygiophore) appear to be correlated with different habitats in rivers and the lake.

Leuciscine fishes represent an important component of freshwater ichthyofauna endemic to northern Mediterranean areas. This lineage shows high intra-specific morphological variability and exhibits high levels of hybridization, two characteristics that contribute to systematic uncertainties, misclassification of taxa and, potentially, the mismanagement of biodiversity. This study focused on brook chub, Squalius lucumonis, an endemic taxon of Central Italy. The taxonomic status of this species has long been questioned, and a hybrid origin from sympatric leusciscines (S. squalus x Rutilus rubilio, or S. squalus x Telestes muticellus) has been hypothesised. A phenotypic (evaluating shape and meristic counts) and genetic (using mitochondrial and nuclear markers) investigation of these four taxa was conducted to test species delimitation in sympatric areas and to evaluate the taxonomic status of S. lucumonis. One hundred and forty-five individuals of all four taxa were collected within streams of the lowest portion of the Tiber River basin and analysed; this region encompasses a large portion of the S. lucumonis distribution. The different morphological and genetic approaches were individually examined, compared, and then combined in a quantitative model to both investigate the limits of each approach and to identify cases of misclassification. The results obtained confirm the cladogenetic non-hybrid origin of S. lucumonis, highlight the need for immediate conservation actions and emphasise the value of an integrated approach in the study of leuciscines evolution.

Background: Amur ide (Leuciscus waleckii) is an economically and ecologically important species in Northern Asia. The Dali Nor population inhabiting Dali Nor Lake, a typical saline-alkaline lake in Inner Mongolia, is well-known for its adaptation to extremely high alkalinity. Genome information is needed for conservation and aquaculture purposes, as well as to gain further understanding into the genetics of stress tolerance. The objective of the study is to sequence the transcriptome and obtain a well-assembled transcriptome of Amur ide., Results: The transcriptome of Amur ide was sequenced using the Illumina platform and assembled into 53,632 cDNA contigs, with an average length of 647 bp and a N50 length of 1,094 bp. A total of 19,338 unique proteins were identified, and gene ontology and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses classified all contigs into functional categories. Open Reading Frames (ORFs) were detected from 34,888 (65.1%) of contigs with an average length of 577 bp, while 9,638 full-length cDNAs were identified. Comparative analyses revealed that 31,790 (59.3%) contigs have a significant similarity to zebrafish proteins, and 27,096 (50.5%), 27,524 (51.3%) and 27,996 (52.2%) to teraodon, medaka and three-spined stickleback proteins, respectively. A total of 10,395 microsatellites and 34,299 SNPs were identified and classified. A dN/dS analysis on unigenes was performed, which identified that 61 of the genes were under strong positive selection. Most of the genes are associated with stress adaptation and immunity, suggesting that the extreme alkaline-saline environment resulted in fast evolution of certain genes., Conclusions: The transcriptome of Amur ide had been deeply sequenced, assembled and characterized, providing a valuable resource for a better understanding of the Amur ide genome. The transcriptome data will facilitate future functional studies on the Amur ide genome, as well as provide insight into potential mechanisms for adaptation to an extreme alkaline-saline environment.

Background: The taxonomic validity of Northeast Indian endemic Mahseer species, Tor progeneius and Neolissochilus hexastichus, has been argued repeatedly. This is mainly due to disagreements in recognizing the species based on morphological characters. Consequently, both the species have been concealed for many decades. DNA barcoding has become a promising and an independent technique for accurate species level identification. Therefore, utilization of such technique in association with the traditional morphotaxonomic description can resolve the species dilemma of this important group of sport fishes., Methodology/principal Findings: Altogether, 28 mahseer specimens including paratypes were studied from different locations in Northeast India, and 24 morphometric characters were measured invariably. The Principal Component Analysis with morphometric data revealed five distinct groups of sample that were taxonomically categorized into 4 species, viz., Tor putitora, T. progeneius, Neolissochilus hexagonolepis and N. hexastichus. Analysis with a dataset of 76 DNA barcode sequences of different mahseer species exhibited that the queries of T. putitora and N. hexagonolepis clustered cohesively with the respective conspecific database sequences maintaining 0.8% maximum K2P divergence. The closest congeneric divergence was 3 times higher than the mean conspecific divergence and was considered as barcode gap. The maximum divergence among the samples of T. progeneius and T. putitora was 0.8% that was much below the barcode gap, indicating them being synonymous. The query sequences of N. hexastichus invariably formed a discrete and a congeneric clade with the database sequences and maintained the interspecific divergence that supported its distinct species status. Notably, N. hexastichus was encountered in a single site and seemed to be under threat., Conclusion: This study substantiated the identification of N. hexastichus to be a true species, and tentatively regarded T. progeneius to be a synonym of T. putitora. It would guide the conservationists to initiate priority conservation of N. hexastichus and T. putitora.

The Squalius alburnoides complex (Steindachner) is one of the most intricate hybrid polyploid systems known in vertebrates. In this complex, the constant switch of the genome composition in consecutive generations, very frequently involving a change on the ploidy level, promotes repetitive situations of potential genomic shock. Previously in this complex, it was showed that in response to the increase in genome dosage, triploids hybrids could regulate gene expression to a diploid state. In this work we compared the small RNA profiles in the different genomic compositions interacting in the complex in order to explore the miRNA involvement in gene expression regulation of triploids. Using high-throughput arrays and sequencing technologies we were able to verify that diploid and triploid hybrids shared most of their sequences and their miRNA expression profiles were high correlated. However, an overall view indicates an up-regulation of several miRNAs in triploids and a global miRNA expression in triploids higher than the predicted from an additive model. Those results point to a participation of miRNAs in the cellular functional stability needed when the ploidy change.

Idiopathic autism, caused by genetic susceptibility interacting with unknown environmental triggers, has increased dramatically in the past 25 years. Identifying environmental triggers has been difficult due to poorly understood pathophysiology and subjective definitions of autism. The use of antidepressants by pregnant women has been associated with autism. These and other unmetabolized psychoactive pharmaceuticals (UPPs) have also been found in drinking water from surface sources, providing another possible exposure route and raising questions about human health consequences. Here, we examined gene expression patterns of fathead minnows treated with a mixture of three psychoactive pharmaceuticals (fluoxetine, venlafaxine & carbamazepine) in dosages intended to be similar to the highest observed conservative estimates of environmental concentrations. We conducted microarray experiments examining brain tissue of fish exposed to individual pharmaceuticals and a mixture of all three. We used gene-class analysis to test for enrichment of gene sets involved with ten human neurological disorders. Only sets associated with idiopathic autism were unambiguously enriched. We found that UPPs induce autism-like gene expression patterns in fish. Our findings suggest a new potential trigger for idiopathic autism in genetically susceptible individuals involving an overlooked source of environmental contamination.

Background: Paedocypris is a newly established genus of fish in Southeast Asia. Paedocypris is characterized by several unique features, including a tiny adult size (thus named miniature fish or minifish), fragmentary habitats of acidic peat blackwater swamps, an unusual reproduction mode and truncated development. These peculiarities lend themselves excellent for studying chromosome evolution and rapid speciation in vertebrates but also make them highly controversial for the phylogenetic position., Methodology and Principal Findings: We have established an organ procedure to prepare chromosome spreads from tiny organs of minifish and performed a cytogenetic study on two species of the genus Paedocypris, namely P. carbunculus (Pc) and P. sp. "Singkep" (Ps). We found 30 and 34 chromosomes in diploid cells of Pc and Ps, respectively, which are unusual in teleost fishes. The diploid metaphase has 5 pairs of metacentrics and 7 pairs of subtelocentrics in Pc compared to 3 pairs of metacentrics and 11 pairs of subtelocentrics in Ps, whereas the haploid metaphase contains 5 metacentrics and 7 subtelocentrics in Pc compared to 3 metacentrics and 11 subtelocentrics Ps. Chromosome behavior in first meiosis revealed the presence of a chromosomal ring consisting of 2 metacentrics in Pc, suggesting that centric fusion rather than fission was responsible for the karyotypic evolution from Ps to Pc. Flow cytometry revealed that Pc had a 45% nuclear staining intensity relative to medaka whose genome is 700 Mb in size and contains 0.81 pg DNA. The Pc genome should have 315 Mb in length and 0.36 pg of DNA, which represent one of the smallest values in vertebrates, suggesting genome miniaturization in this organism., Conclusions: Our data demonstrate that gross chromosome rearrangements and genome miniaturization have accompanied the evolution of Paedocypris fishes. Our data also place Paedocypris outside currently described taxa of the Cypriniformes.

Investigating processes and timing of differentiation of organisms is critical in the understanding of the evolutionary mechanisms involved in microevolution, speciation, and macroevolution that generated the extant biodiversity. From this perspective, the Telestes genus is of special interest: the Telestes species have a wide distribution range across Europe (from the Danubian district to Mediterranean districts) and have not been prone to translocation. Molecular data (mtDNA: 1,232 bp including the entire Cyt b gene; nuclear genome: 11 microsatellites) were gathered from 34 populations of the Telestes genus, almost encompassing the entire geographic range. Using several phylogenetic and molecular dating methods interpreted in conjunction with paleoclimatic and geomorphologic evidence, we investigated the processes and timing of differentiation of the Telestes lineages. The observed genetic structure and diversity were largely congruent between mtDNA and microsatellites. The Messinian Salinity Crisis (Late Miocene) seems to have played a major role in the speciation processes of the genus. Focusing on T. souffia, a species occurring in the Danube and Rhone drainages, we were able to point out several specific events from the Pleistocene to the Holocene that have likely driven the differentiation and the historical demography of this taxon. This study provides support for an evolutionary history of dispersal and vicariance with unprecedented resolution for any freshwater fish in this region.

Divergent natural selection rooted in differential resource use can generate and maintain intraspecific eco-morphological divergence (i.e., resource polymorphism), ultimately leading to population splitting and speciation. Differing bottom environments create lake habitats with different benthos communities, which may cause selection in benthivorous fishes. Here, we document the nature of eco-morphological and genetic divergence among local populations of the Japanese gudgeon Sarcocheilichthys (Cyprinidae), which inhabits contrasting habitats in the littoral zones (rocky vs. pebbly habitats) in Lake Biwa, a representative ancient lake in East Asia. Eco-morphological analyses revealed that Sarcocheilichthys variegatus microoculus from rocky and pebbly zones differed in morphology and diet, and that populations from rocky environments had longer heads and deeper bodies, which are expected to be advantageous for capturing cryptic and/or attached prey in structurally complex, rocky habitats. Sarcocheilichthys biwaensis, a rock-dwelling specialist, exhibited similar morphologies to the sympatric congener, S. v. microoculus, except for body/fin coloration. Genetic analyses based on mitochondrial and nuclear microsatellite DNA data revealed no clear genetic differentiation among local populations within/between the gudgeon species. Although the morphogenetic factors that contribute to morphological divergence remain unclear, our results suggest that the gudgeon populations in Lake Biwa show a state of resource polymorphism associated with differences in the bottom environment. This is a novel example of resource polymorphism in fish within an Asian ancient lake, emphasizing the importance and generality of feeding adaptation as an evolutionary mechanism that generates morphological diversification.

The Asian cyprinid fish, the topmouth gudgeon (Pseudorasbora parva), was introduced into Europe in the 1960s. A highly invasive freshwater fish, it is currently found in at least 32 countries outside its native range. Here we analyse a 700 base pair fragment of the mitochondrial cytochrome b gene to examine different models of colonisation and spread within the invasive range, and to investigate the factors that may have contributed to their invasion success. Haplotype and nucleotide diversity of the introduced populations from continental Europe was higher than that of the native populations, although two recently introduced populations from the British Isles showed low levels of variability. Based on coalescent theory, all introduced and some native populations showed a relative excess of nucleotide diversity compared to haplotype diversity. This suggests that these populations are not in mutation-drift equilibrium, but rather that the relative inflated level of nucleotide diversity is consistent with recent admixture. This study elucidates the colonisation patterns of P. parva in Europe and provides an evolutionary framework of their invasion. It supports the hypothesis that their European colonisation was initiated by their introduction to a single location or small geographic area with subsequent complex pattern of spread including both long distance and stepping-stone dispersal. Furthermore, it was preceded by, or associated with, the admixture of genetically diverse source populations that may have augmented its invasive-potential.