Isolation, identification and characterisation of Pseudomonas koreensis CM-01 isolated from diseased Malaysian mahseer (Tor tambroides)

Pseudomonas species are one of the most threatening fish pathogens which reside a wide range of environments. In this study, the dominant bacteria were isolated from diseased Malaysian mahseer (Tor tambroides) and tentatively named CM-01. It was identified as Pseudomonas koreensis based on its biochemical, morphological, genetic and physiological information. Its pathogenicity was found to be correlated with twelve virulence genes identified including iron uptake, protease, acylhomoserine lactone synthase gacS/gacA component regulation system, type IV secretion system, hydrogen cyanide production, exolysin, alginate biosynthesis, flagella and pili. The median lethal dose (LD50) for the CM-01 isolate on Malaysian mahseer was documented at 5.01 × 107 CFU/mL. The experimental infection revealed that CM-01 led to significant histological lesions in the fish, ultimately resulting in death. These lesions comprise necrosis, tissue thickening and aggregation. Drug sensitivity tests had shown its susceptibility to beta-lactam combination agents and further suggest its drug of choice. Its growing features had shown its growth at optimal temperature and pH. To the best of our knowledge, this is the first report of P. koreensis linked to diseased T. tambroides. STATEMENT OF RELEVANCE: In this research, a novel strain of Pseudomonas koreensis, CM-01 was isolated from diseased T. tambroides for the first time. The antimicrobial susceptibility, pathogenicity, virulence genes and growth characteristics of CM-01 were studied. These findings established a scientific foundation for the recognition of P. koreensis and the management of fish infections caused by this pathogen.

Toxicogenomics of Five Cytostatics in Fathead Minnow (Pimephales promelas) Larvae

In this study, the toxicogenomic effects of five cytostatics (tamoxifen, methotrexate, capecitabine, cyclophosphamide, and ifosfamide) on fathead minnow (Pimephales promelas) larvae were evaluated. Post-fertilization eggs were exposed to increasing concentrations of the drugs for six days. The expression levels of two genetic biomarkers for toxicity and four thyroid hormone-related gene pathways were measured. Interestingly, the results showed that all concentrations of the five cytostatics affect the transcription levels of both toxicity biomarker genes. Additionally, the thyroid hormone-related genes had different expression levels than the control, with the most significant changes observed in those larvae exposed to cyclophosphamide and ifosfamide. While a previous study found no effects on fish morphology, this study suggests that the five cytostatics modify subtle molecular responses of P. promelas, highlighting the importance of assessing multibiological level endpoints throughout the lifecycle of animals to understand the full portrait of potential effects of cytostatics and other contaminants.

Regulation of sirt1 and foxO1 in glucose metabolism of Megalobrama amblycephala

Both silent information regulator 2 homolog 1 (sirt1) and forkhead box transcription factor 1 (foxO1) are crucial transcription factors involved in glucolipid metabolism and energy regulation. The presnt study aimed to understand their regulatory roles in glucose metabolism. Molecular cloning and sequencing of sirt1 gene of Megalobrama amblycephala (masirt1) was conducted and cellular localization of both the factors were analysed. Their effects and action patterns in the glucose metabolism of Megalobrama amblycephala (M. amblycephala) were investigated through acute and long-term glucose tolerance assays. The results revealed that the full-length masirt1 cDNA sequence was 2350 bp and closely related to Sinocyclocheilus rhinocerous. Sirt1 and foxO1 were found to be mutually dependent and localized in the nucleus. Acute glucose tolerance tests revealed that the expression levels of both factors in the liver of M. amblycephala showed an initial increase followed by a decrease. Plasma glucose levels in M. amblycephala significantly increased at 2 and 12 h (P < 0.05). In a long-term breeding experiment with high-sugar feeding, the expressions of the sirt1 and foxO1 genes in the kidney and intestine of M. amblycephala exhibited synergistic changes. The 51WS groups had significantly higher levels of sirt1 and foxO1 gene expression in the kidney and intestine compared to the 0WS and 17WS groups (P < 0.05). Overall, masirt1 is evolutionarily highly conserved, and the interaction site of sirt1 and foxO1 is located in the nucleus. In long-term hyperglycemic regulation, sirt1 and foxO1 exhibit synergistic regulatory effects in the kidney and intestine of M. amblycephala. This study provides insights into how sirt1 and foxO1 regulate glucose metabolism in M. amblycephala.

Genomic insights into the seawater adaptation in Cyprinidae

BACKGROUND

Cyprinidae, the largest fish family, encompasses approximately 367 genera and 3006 species. While they exhibit remarkable adaptability to diverse aquatic environments, it is exceptionally rare to find them in seawater, with the Far Eastern daces being of few exceptions. Therefore, the Far Eastern daces serve as a valuable model for studying the genetic mechanisms underlying seawater adaptation in Cyprinidae.

RESULTS

Here, we sequenced the chromosome-level genomes of two Far Eastern daces (Pseudaspius brandtii and P. hakonensis), the two known cyprinid fishes found in seawater, and performed comparative genomic analyses to investigate their genetic mechanism of seawater adaptation. Demographic history reconstruction of the two species reveals that their population dynamics are correlated with the glacial-interglacial cycles and sea level changes. Genomic analyses identified Pseudaspius-specific genetic innovations related to seawater adaptation, including positively selected genes, rapidly evolving genes, and conserved non-coding elements (CNEs). Functional assays of Pseudaspius-specific variants of the prolactin (prl) gene showed enhanced cell adaptation to greater osmolarity. Functional assays of Pseudaspius specific CNEs near atg7 and usp45 genes suggest that they exhibit higher promoter activity and significantly induced at high osmolarity.

CONCLUSIONS

Our results reveal the genome-wide evidence for the evolutionary adaptation of cyprinid fishes to seawater, offering valuable insights into the molecular mechanisms supporting the survival of migratory fish in marine environments. These findings are significant as they contribute to our understanding of how cyprinid fishes navigate and thrive in diverse aquatic habitats, providing useful implications for the conservation and management of marine ecosystems.

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

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

Molecular Identification of Acrossocheilus jishouensis (Teleostei: Cyprinidae) and Its Complete Mitochondrial Genome

Molecular identification, such as DNA barcoding, is a useful tool that is widely applied in distinguishing species. To identify the cyprinid Acrossocheilus jishouensis, which was previously known to be restricted to only its type locality, we conducted molecular identification of this species based on 23 samples in five localities. Molecular identification based on the mitochondrial COI gene sequence showed that the morphologically similar samples from the five populations were all A. jishouensis, as the mean genetic distances between populations were very small (0.1-1.6%); thus, the distribution of this species was substantially expanded. The whole mitochondrial genome of one sample was also assembled, which was 16,594 bp in length and consisted of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region. All PCGs began with ATG except the COI gene, which started with GTG; seven PCGs used the complete stop codon TAA, while four terminated in T(AA) and two ended with TAG. The overall base composition reflected a higher proportion of A+T than G+C and a positive AT-skew and negative GC-skew pattern except for the opposite in ND6. Phylogenetic relationships inferred using BI and ML methods revealed that both Acrossocheilus and Onychostoma were nonmonophyletic, which indicated that the traditional diagnoses between these two genera need to be assessed further. The results of this study not only expanded the known distribution ranges of A. jishouensis, but also provided a valuable data resource for future molecular and evolutionary studies of Acrossocheilus and other cyprinids in Barbinae.

Uncoordinated 51-like kinase 1a/b and 2 in fish Megalobrama amblycephala: Molecular cloning, functional characterization, and their potential roles in glucose metabolism

Uncoordinated (Unc) 51-like kinase (ulk1) and ulk2 are closely involved in autophagy activation, but little is known about their roles in regulating glucose homeostasis. In this study, the genes of ulk1a, ulk1b and ulk2 were cloned and characterized in fish Megalobrama amblycephala. All the three genes shared the approximate N-terminal kinase domain and the C-terminal Atg1-like_tMIT domain structure, while the amino acid sequence identity of them are different between M. amblycephala and other vertebrates. Their transcripts were widely observed in various tissues (brain, muscle, gill, heart, spleen, eye, liver, intestine, abdominal adipose and kidney), but differed in tissue expression patterns. During the glucose tolerance test and the insulin tolerance test, the up-regulated transcriptions of ulk1a, ulk1b and ulk2 were all found despite some differences in the temporal patterns. At the same time, the activities of glycolytic enzymes like hexokinase and phosphofructokinase both showed parallel increases. Furthermore, the feeding of a high-carbohydrate diet decreased the transcriptions of ulk1a, ulk1b and ulk2. Collectively, this study demonstrated that ulk1a, ulk1b and ulk2 in M. amblycephala had similar molecular characterizations, but with different conservation and tissue expression patterns. In addition, ulk1/2 might play important roles in maintaining the glucose homeostasis in fish through regulating the glycolytic pathway.

Transcriptome analysis of Schizothorax oconnori (Cypriniformes: Cyprinidae) oocytes: The role of K+ in promoting yolk globule fusion and regulating oocyte maturation

Schizothorax oconnori (S. oconnori) is an economically important fish in Tibet. Oocyte maturation is a physiological process that is of great significance to reproduction and seed production in S. oconnori, yet little is currently known regarding the molecular mechanisms of oocyte development in this species. To identify candidate genes involved in reproduction of female fish, a combination of PacBio and Illumina HiSeq technologies was employed to provide deep coverage of the oocyte transcriptome. Transcriptome analysis revealed several candidate genes that are potentially involved in the regulation of oocyte maturation in S. oconnori, including GIRK1, CHRM3, NPY2R, GABRA3, GnRH3, mGluR1α, GPER1, GDF9, HSP90, and ESR2. Genes that are significantly expressed during oocyte maturation mainly contribute to the GPCR signaling pathway and the estrogen signaling pathway. Neurotransmitter (Ach, NPY, and GABA) and peptide hormone (GnRH3) binding to G protein-coupled receptors (GPCRs) frees G-protein βγ subunits to interact with the G protein-gated inward rectifier K+ channel 1 (GIRK1). This process helps release K+ from granulosa cells to maturing oocytes, allowing yolk globule fusion. This mechanism may play an important role in oocyte maturation in S. oconnori. In conclusion, this study provides a valuable basis for deciphering the reproductive system in S. oconnori during the oocyte maturation process.

Complete mitogenome of Lepidopygopsis typus, an evolutionarily-distinct, endangered cyprinid fish from the Western Ghats Biodiversity Hotspot: Phylogenetic relationships and implications for conservation

The mitogenome of Lepidopygopsis typus, an evolutionarily distinct, endangered, cyprinid fish from the Western Ghats Biodiversity Hotspot, was characterized. Total length of the mitogenome was 16,729 bp, and comprised of 13 protein coding, 22 tRNA and two rRNA genes. Thirteen protein coding genes contributed a total nucleotide length of 11,408 bp, which coded for a total of 3794 amino acids and accounting for 68.2 % of the mitogenome. A maximum likelihood phylogenetic tree based on mitogenomes of members of cyprinid subfamilies Torinae, Schizothoracinae and Schizopygopsinae, recovered L. typus in a well-resolved clade of Torinae with members of the Western Ghats endemic genus Hypselobarbus as sister taxa. Selection analysis suggested that the branch of L. typus and Hypselobarbus spp. was under intensified selection with reference to other members of Torinae. There were 19 codons under diversifying selection in L. typus, which could be the result of positive selection for adapting to high-altitude, upstream tributaries of Periyar River, where the species has a restricted distribution. With respect to Hypselobarbus spp., L. typus was under relaxed selection with about 68 % of codons experiencing neutral evolution. Restricted distribution, low population size and relaxation in selection can likely trigger extinction in L. typus, and therefore urgent conservation and monitoring plans are required to secure the future of this evolutionary distinct, and globally endangered species.

Ultrafast sound production mechanism in one of the smallest vertebrates

Motion is the basis of nearly all animal behavior. Evolution has led to some extraordinary specializations of propulsion mechanisms among invertebrates, including the mandibles of the dracula ant and the claw of the pistol shrimp. In contrast, vertebrate skeletal movement is considered to be limited by the speed of muscle, saturating around 250 Hz. Here, we describe the unique propulsion mechanism by which Danionella cerebrum, a miniature cyprinid fish of only 12 mm length, produces high amplitude sounds exceeding 140 dB (re. 1 µPa, at a distance of one body length). Using a combination of high-speed video, micro-computed tomography (micro-CT), RNA profiling, and finite difference simulations, we found that D. cerebrum employ a unique sound production mechanism that involves a drumming cartilage, a specialized rib, and a dedicated muscle adapted for low fatigue. This apparatus accelerates the drumming cartilage at over 2,000 g, shooting it at the swim bladder to generate a rapid, loud pulse. These pulses are chained together to make calls with either bilaterally alternating or unilateral muscle contractions. D. cerebrum use this remarkable mechanism for acoustic communication with conspecifics.

Evolutionary mechanisms and practical significance of reproductive success and clonal diversity in unisexual vertebrate polyploids

Unisexual reproduction is generally relevant to polyploidy, and unisexual vertebrates are often considered an evolutionary "dead end" due to the accumulation of deleterious mutations and absence of genetic diversity. However, some unisexual polyploids have developed strategies to avoid genomic decay, and thus provide ideal models to unveil unexplored evolutionary mechanisms, from the reproductive success to clonal diversity creation. This article reviews the evolutionary mechanisms for overcoming meiotic barrier and generating genetic diversity in unisexual vertebrates, and summarizes recent research advancements in the polyploid Carassius complex. Gynogenetic gibel carp (Carassius gibelio) is a unique amphitriploid that has undergone a recurrent autotriploidy and has overcome the bottleneck of triploid sterility via gynogenesis. Recently, an efficient strategy in which ploidy changes, including from amphitriploid to amphitetraploid, then from amphitetraploid to novel amphitriploid, drive unisexual-sexual-unisexual reproduction transition and clonal diversity has been revealed. Based on this new discovery, multigenomic reconstruction biotechnology has been used to breed a novel strain with superior growth and stronger disease resistance. Moreover, a unique reproduction mode that combines both abilities of ameiotic oogenesis and sperm-egg fusion, termed as ameio-fusiongensis, has been discovered, and it provides an efficient approach to synthesize sterile allopolyploids. In order to avoid ecological risks upon escape and protect the sustainable property rights of the aquaculture seed industry, a controllable fertility biotechnology approach for precise breeding is being developed by integrating sterile allopolyploid synthesis and gene-editing techniques. This review provides novel insights into the origin and evolution of unisexual vertebrates and into the attempts being made to exploit new breeding biotechnologies in aquaculture.

Development of novel SSR markers and validation by assessing the genetic diversity of endangered Deccan mahseer, Tor khudree

BACKGROUND

The Deccan mahseer, Tor khudree (Sykes, 1839) is a potential game and food fish species belonging to the family cyprinidae and is categorized as endangered. Its distribution is restricted to southern part of India, specifically to Peninsular Rivers. This study is first to assess the genetic diversity and differentiation in Tor khudree by developing novel simple sequence repeat (SSR) markers.

METHODS AND RESULTS

Low depth next generation sequencing followed by sequence analysis in MISA software identified 187,649 SSRs. The novel fourteen validated SSR loci were used for population genetic analysis. All of the SSR loci were highly informative with mean PIC > 0.5. High mean allelic richness (9.29) observed heterozygosity (0.98) and expected heterozygosity (0.79) were observed across the loci. However, genetic differentiation was low but significant (0.052). Negative FIS values were observed in both locus-wise and populations indicating the presence of high heterozygosity. Intrapopulation variation was found to be high (96.29%). The population structure revealed two genetic stocks.

CONCLUSIONS

The results from the present study including the highly polymorphic markers developed would be a useful resource for further research on population genetics and conservation genetics of the species.

Chromosomes of Asian cyprinid fishes: Novel insight into the chromosomal evolution of Labeoninae (Teleostei, Cyprinidae)

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.

A high-quality reference genome for the common creek chub, Semotilus atromaculatus

Creek chub (Semotilus atromaculatus) is a leuciscid minnow species commonly found in anthropogenically disturbed environments, making it an excellent model organism to study human impacts on aquatic systems. Genomic resources for creek chub and other leuciscid species are currently limited. However, advancements in DNA sequencing now allow us to create genomic resources at a historically low cost. Here, we present a high quality, 239 contig reference genome for the common creek chub, created with PacBio HiFi sequencing. We compared the assembly quality of two pipelines: Pacific Biosciences' Improved Phase Assembly (873 contigs) and Hifiasm (239 contigs). Quality and completeness of this genome is comparable to the zebrafish (Danioninae) and fathead minnow (Leuciscidae) genomes. The creek chub genome is highly syntenic to the zebrafish and fathead minnow genomes, and while our assembly does not resolve into the expected 25 chromosomes, synteny with zebrafish suggests that each creek chub chromosome is likely represented by 1-4 large contigs in our assembly. This reference genome is a valuable resource that will enhance genomic biodiversity studies of creek chub and other nonmodel leuciscid species common to disturbed environments.

Population structure and adaptability analysis of Schizothorax o'connori based on whole-genome resequencing

BACKGROUND

Schizothorax o'connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o'connori using the whole-genome resequencing method.

RESULTS

The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o'connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism.

CONCLUSION

Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o'connori. There was no obvious genetic differentiation at the genome level between the S. o'connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o'connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

Chromosome-level genome assembly of Acrossocheilus fasciatus using PacBio sequencing and Hi-C technology

Acrossocheilus fasciatus (Cypriniformes, Cyprinidae) is emerged as a newly commercial stream fish in the south of China with high economic and ornamental value. In this study, a chromosome-level reference genome of A. fasciatus was assembled using PacBio, Illumina and Hi-C sequencing technologies. As a result, a high-quality genome was generated with a size of 879.52 Mb (accession number: JAVLVS000000000), scaffold N50 of 32.7 Mb, and contig N50 of 32.7 Mb. The largest and smallest scafford was 60.57 Mb and 16 kb, respectively. BUSCO analysis showed a completeness score of 98.3%. Meanwhile, the assembled sequences were anchored to 25 pseudo-chromosomes with an integration efficiency of 96.95%. Additionally, we found approximately 390.91 Mb of repetitive sequences that accounting for 44.45% of the assembled genome, and predicted 24,900 protein-coding genes. The available genome reported in the present study provided a crucial resource to further investigate the regulation mechanism of genetic diversity, sexual dimorphism and evolutionary histories.

Transcriptomic and proteomic strategies to reveal the mechanism of Gymnocypris przewalskii scale development

BACKGROUND

Fish scales are typical products of biomineralization and play an important role in the adaptation of fish to their environment. The Gymnocypris przewalskii scales are highly specialized, with scales embedded in only specific parts of the dermis, such as the areas around the anal fin and branchiostegite, making G. przewalskii an ideal material for biomineralization research. In this study, we aimed to unveil genes and pathways controlling scale formation through an integrated analysis of both transcriptome and proteome, of which G. przewalskii tissues of the dorsal skin (no scales) and the rump side skin (with scales) were sequenced. The sequencing results were further combined with cellular experiments to clarify the relationship between genes and signaling pathways.

RESULTS

The results indicated the following: (1) a total of 4,904 differentially expressed genes were screened out, including 3,294 upregulated genes and 1,610 downregulated genes (with a filtering threshold of |log2Fold-Change|> 1 and p-adjust < 0.05). The identified differentially expressed genes contained family members such as FGF, EDAR, Wnt10, and bmp. (2) A total of 535 differentially expressed proteins (DEPs) were filtered out from the proteome, with 204 DEPs downregulated and 331 DEPs upregulated (with a filtering threshold of |Fold-Change|> 1.5 and p < 0.05). (3) Integrated analyses of transcriptome and proteome revealed that emefp1, col1a1, col6a2, col16a1, krt8, and krt18 were important genes contributing to scale development and that PI3K-AKT was the most important signaling pathway involved. (4) With the use of the constructed G. przewalskii fibroblast cell line, emefp1, col1a1, col6a2, col16a1, krt8, and krt18 were confirmed to be positively regulated by the PI3K-AKT signaling pathway.

CONCLUSION

This study provides experimental evidence for PI3K-AKT controlled scale development in G. przewalskii and would benefit further study on stress adaptation, scale biomineralization, and the development of skin appendages.

A new cavefish of Sinocyclocheilus (Teleostei: Cypriniformes: Cyprinidae) from the Nanpanjiang River in Guizhou, China

A new species, Sinocyclocheilus xingyiensis, is described based on specimens collected from a karst cave in Guizhou Province, China. The authors used an integrated taxonomic approach, including morphological and molecular data, to identify the new species as a member of the Sinocyclocheilu angularis group, and it can be distinguished from all other members of this group by a combination of the following features: two pairs of long barbels and long pectoral fins, 42-46 lateral-line scales, 7 (13-14) on outer (inner) side of the first gill arch and 35 (14-15 + 4 + 16 - 17) vertebrae. Phylogenetic analyses based on the cytochrome b (cyt b) gene fragment suggest that S. xingyiensis is a sister lineage to Sinocyclocheilus flexuosdorsalis. The genetic distance (Kimura 2-parameter) between the S. xingyiensis and S. angularis groups of Sinocyclocheilus species based on cyt b gene fragment ranged from 1.2% to 15.4%.

First record of the genus Schizothorax in rivers of the Central Himalayas in China

A total of 40 specimens of the genus Schizothorax (Cyprinidae: Schizothoracinae), including 10 matured males, 19 matured females and 11 juveniles, were captured from the Kirong Tsangpo River in China, which is located in the southern slope of the Central Himalayas. These specimens are identified as Schizothorax richardsonii (Grey, 1832) based on morphological characters and molecular analyses using mitochondrial Cyt b gene sequences. The Kirong population of S. richardsonii is relatively isolated from other populations in the Himalayas and has low genetic diversity. This is the first record of the genus Schizothorax fish in rivers of the Central Himalayas in China. As S. richardsonii is a vulnerable species on the IUCN Red List, a protection plan should be conducted to reduce the impact of anthropogenic disturbance by monitoring the natural population dynamics and assessing the ecological determinants of its distribution.

Genetic erosion in an endangered desert fish during a megadrought despite long-term supportive breeding

Human water use combined with a recent megadrought have reduced river and stream flow through the southwest United States and led to periodic drying of formerly perennial river segments. Reductions in snowmelt runoff and increased extent of drying collectively threaten short-lived, obligate aquatic species, including the endangered Rio Grande silvery minnow (Hybognathus amarus). This species is subject to boom-and-bust population dynamics, under which large fluctuations in abundance are expected to lower estimates of effective population size and erode genetic diversity over time. Rates of diversity loss are also affected by additions of hatchery-origin fish used to supplement the wild population. We used demographic and genetic data from wild and hatchery individuals to examine the relationship of genetic diversity and effective population size to abundance over the last two decades. Genetic diversity was low during the early 2000s, but diversity and demographic metrics stabilized after the hatchery program was initiated and environmental conditions improved. Yet, from 2017 onward, allelic diversity declined (Cohen's d = 1.34) and remained low despite hatchery stocking and brief wild population recovery. Across the time series, single-sample estimates of effective population size based on linkage disequilibrium (LD Ne ) were positively associated (r = 0.53) with wild abundance and total abundance, but as the proportion of hatchery-origin spawners increased, LD Ne declined (r = -0.55). Megadrought limited wild spawner abundance and precluded refreshment of hatchery brood stocks with wild fish; hence, we predict a riverine population increasingly dominated by hatchery-origin individuals and accelerated loss of genetic diversity despite supplementation. We recommend an adaptive and accelerated management plan that integrates river flow management and hatchery operations to slow the pace of genetic diversity loss exacerbated by megadrought.

Transcriptome Analysis Reveals the Molecular Mechanisms Underlying Growth Superiority in a Novel Gymnocypris Hybrid, Gymnocypris przewalskii ♀ × Gymnocypris eckloni ♂

Artificial hybrid breeding can optimize parental traits to cultivate excellent hybrids with enhanced economic value. In this study, we investigated the growth performance and transcriptomes of Gymnocypris przewalskii (♀) and Gymnocypris eckloni (♂) and their F1 hybrid fishes. Hatched individuals of G. przewalskii (GP) and G. eckloni (GE) of the same size and their F1 hybrids (GH) were separately cultured for eight months in three cement tanks (n = 3). The growth indexes were measured, which showed that the growth rate of the groups was GE > GH > GP, while the survival rate was GH > GE > GP. The RNA-Seq data analysis of the muscles from the three Gymnocypris fish strains revealed that gene transcription has a significant impact on F1 hybrid fish and its parents. The differentially expressed genes (DEGs) in GH show less differences with GP, but more with GE. qRT-PCR was used to confirm the expression profiles of the chosen DEGs, and the results showed positive correlations with the RNA-seq data. KEGG enrichment results indicated that the DEGs were related to a variety of molecular functions, such as glycolysis/gluconeogenesis, arachidonic acid formation, citrate cycle, and the MAPK, PI3K-Akt, or mTOR signal pathways. Subsequent analysis indicated that there may be a significant correlation between the differential expression of IGF2 and a difference in the growth of GE and GP.

Four new complete mitochondrial genomes of Gobioninae fishes (Teleostei: Cyprinidae) and their phylogenetic implications

The subfamily Gobioninae is one of the most diverse fish groups within Cyprinidae. Their taxonomy and phylogenetic relationships are not completely resolved. In this study, the complete mitochondrial genomes (mitogenome) of four Gobioninae species (Microphysogobio elongatus, Microphysogobio chinssuensis, Gobio rivuloides and Rhinogobio nasutus) were sequenced and compared. The mitogenomes of four species ranges from 16603 bp to 16609 bp in length, consisting of 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a control region. Most PCGs had significant codon usage bias. Except for the tRNASer (GCT), all the nucleotide substitutions of tRNA loops higher than the stems could fold into a stable secondary structure. The nucleotide compositions of Gobioninae mitogenome were biased toward A/T, and NAD4 was subjected to low purification selection and had a faster evolution rate among 13 PCGs. Bayesian inference and maximum likelihood phylogenetic analyses showed the consistent results. The four sequenced species clustered together with their congener species. However, more samples and mitogenome data are needed to untangle the phylogenetic relationships among genera Microphysogobio, Romanogobio, Hugobio, Biwia and Platysmacheilus.

Genetic approaches reveal a healthy population and an unexpectedly recent origin for an isolated desert spring fish

Foskett Spring in Oregon's desert harbors a historically threatened population of Western Speckled Dace (Rhinichthys klamathensis). Though recently delisted, the dace's recruitment depends upon regular removal of encroaching vegetation. Previous studies assumed that Foskett Dace separated from others in the Warner Valley about 10,000 years ago, thereby framing an enigma about the population's surprising ability to persist for so long in a tiny habitat easily overrun by plants. To investigate that persistence and the effectiveness of interventions to augment population size, we assessed genetic diversity among daces inhabiting Foskett Spring, a refuge at Dace Spring, and three nearby streams. Analysis revealed a robust effective population size (Ne) of nearly 5000 within Foskett Spring, though Ne in the Dace Spring refuge is just 10% of that value. Heterozygosity is slightly lower than expected based on random mating at all five sites, indicating mild inbreeding, but not at a level of concern. These results confirm the genetic health of Foskett Dace. Unexpectedly, genetic differentiation reveals closer similarity between Foskett Dace and a newly discovered population from Nevada's Coleman Creek than between Foskett Dace and dace elsewhere in Oregon. Demographic modeling inferred Coleman Creek as the ancestral source of Foskett Dace fewer than 1000 years ago, much more recently than previously suspected and possibly coincident with the arrival of large herbivores whose grazing may have maintained open water suitable for reproduction. These results solve the enigma of persistence by greatly shortening the duration over which Foskett Dace have inhabited their isolated spring.

Cloning, expression and functional characterization of recombinant tumor necrosis factor α1 (TNFα1) from white crucian carp in gut immune regulation

TNFα is one of important cytokines belonging to TNF superfamily, which can exhibit a pleiotropic effect in immune modulation, homeostasis as well as pathogenesis. However, its immunoregulatory function on mucosal immunity in fish gut are still unclear. In this study, we aimed to investigated the immunoregulatory role of TNFα1 in midgut of white crucian carp (WCC). WCC-TNFα1 sequence and its deduced structure were firstly identified in WCC. Then, tissue-specific analysis revealed that high-level WCC-TNFα1 expression was detected in gill. After Aeromonas hydrophila and lipopolysaccharide (LPS) stimulated, increased trends of WCC-TNFα1 expressions were detected in immune-related tissues and cultured fish cells, respectively. WCC anal-intubated with WCC-TNFα1 fusion protein showed the increased levels of edema and fuzzy appearance in impaired villi, along with atrophy and reduction of goblet cells (GC). Moreover, the expression levels of tight junction (TJ) genes and mucin genes were consistently lower than those of the control (P < 0.05). WCC-TNFα1 treatment could sharply decrease antioxidant status in midgut, while the expression levels of caspase (CASP) genes, unfolded protein response (UPR) genes and redox response genes increased dramatically. Our results suggested that WCC-TNFα1 could exhibit a detrimental effect on antioxidant and mucosal immune regulation in midgut of WCC.

Environmental DNA unveiling the fish community structure and diversity features in the Yangtze River basin

Fish, as top predators in aquatic ecosystems, play an important role in maintaining the structure and functioning of these ecosystems, making their diversity a topic of great interest. This study focused on the Yangtze River Basin to investigate the fish community structure and diversity using environmental DNA (eDNA) technology. The results showed that a total of 71616 fish operational taxonomic units (OTUs) and 90 fish belonging to 23 families were detected, with the Cyprinidae family being the dominant group, followed by the Cobitidae, Amblycipitidae, etc. Compared to historical traditional morphological fish surveys, the quantity of fish detected using eDNA was relatively low, but the overall distribution pattern of fish communities was generally consistent. The highest fish Shannon-Wiener diversity index in the Yangtze River Basin sites reaches 2.60 with an average value of 1.25. The fish diversity index was higher in the downstream compared to the middle and upstream regions, and there were significant differences among different sampling sites. Significant environmental factors influencing α-diversity included chlorophyll-a, chemical oxygen demand, dissolved oxygen, total nitrogen, and elevation. Non-metric multidimensional scaling (NMDS) analysis revealed significant differences in fish community composition between the upstream and middle/lower reaches of the Yangtze River, while the composition of fish communities in the middle and lower reaches was more similar. Redundancy analysis (RDA) indicated that total organic carbon (TOC) was positively correlated with fish community distribution in the upstream, while water temperature and NO3-N were negatively correlated with fish distribution in the upstream. NH3-N and CODMn were negatively correlated with fish distribution in the middle and downstream regions, indicating a relatively severe water pollution in these areas. Additionally, fish communities in the Yangtze River displayed a typical distance decay pattern.

Chromosome-level genome assembly of chub mackerel (Scomber japonicus) from the Indo-Pacific Ocean

Chub mackerels (Scomber japonicus) are a migratory marine fish widely distributed in the Indo-Pacific Ocean. They are globally consumed for their high Omega-3 content, but their population is declining due to global warming. Here, we generated the first chromosome-level genome assembly of chub mackerel (fScoJap1) using the Vertebrate Genomes Project assembly pipeline with PacBio HiFi genomic sequencing and Arima Hi-C chromosome contact data. The final assembly is 828.68 Mb with 24 chromosomes, nearly all containing telomeric repeats at their ends. We annotated 31,656 genes and discovered that approximately 2.19% of the genome contained DNA transposon elements repressed within duplicated genes. Analyzing 5-methylcytosine (5mC) modifications using HiFi reads, we observed open/close chromatin patterns at gene promoters, including the FADS2 gene involved in Omega-3 production. This chromosome-level reference genome provides unprecedented opportunities for advancing our knowledge of chub mackerels in biology, industry, and conservation.

Post-Subfunctionalization Functions of HIF-1αA and HIF-1αB in Cyprinid Fish: Fine-Tuning Mitophagy and Apoptosis Regulation Under Hypoxic Stress

Hypoxia-inducible factor 1 (HIF-1) is a crucial transcriptional factor that can restore oxygen balance in the body by regulating multiple vital activities. Two HIF-1α copies were retained in cyprinid fish after experiencing a teleost-specific genome duplication. How the "divergent collaboration" of HIF-1αA and HIF-1αB proceeds in regulating mitophagy and apoptosis under hypoxic stress in cells of cyprinid fish remains unclear. In this study, zebrafish HIF-1αA/B expression plasmids were constructed and transfected into the epithelioma papulosum cyprini cells and were subjected to hypoxic stress. HIF-1αA induced apoptosis through promoting ROS generation and mitochondrial depolarization when cells were subjected to oxygen deficiency. Conversely, HIF-1αB was primarily responsible for mitophagy induction, prompting ATP production to mitigate apoptosis. HIF-1αA did not induce mitophagy in the mitochondria and lysosomes co-localization assay but it was involved in the regulation of different mitophagy pathways. Over-expression of HIF-1αA increased the expression of bnip3, fundc1, Beclin1, and foxo3, suggesting it has a dual role in mitochondrial autophagy and cell death. Each duplicated copy also experienced functional divergence and target shifting in the regulation of complexes in the mitochondrial electron transport chain (ETC). Our findings shed light on the post-subfunctionalization function of HIF-1αA and HIF-1αB in zebrafish to fine-tune regulation of mitophagy and apoptosis following hypoxia exposure.

Muscle Transcriptome Sequencing Revealed Thermal Stress-Responsive Regulatory Genes in Farmed Rohu, Labeo rohita (Hamilton, 1822)

Rohu, Labeo rohita, is one of the most important aquaculture species in the Indian subcontinent. Understanding the molecular-level physiological responses to thermal stress or climate change is essential. In the present work, transcriptome sequencing was carried out in the muscle tissue of the rohu in response to heat stress (35 °C) in comparison with the control (28 °C). A total of 125 Gb of sequence data was generated, and the raw-reads were filtered and trimmed, which resulted in 484 million quality reads. Reference-based assembly of reads was performed using L. rohita genome, and a total of 90.17% of reads were successfully mapped. A total of 37,462 contigs were assembled with an N50 value of 1854. The differential expression analysis revealed a total of 107 differentially expressed genes (DEGs) (15 up-, 37 down-, and 55 neutrally regulated) as compared to the control group (Log2FC > 2, P < 0.05). Gene enrichment analysis of DEGs indicates that transcripts were associated with molecular, biological, and cellular activities. The randomly selected differentially expressed transcripts were validated by RT-qPCR and found consistent expression patterns in line with the RNA-seq data. Several transcripts such as SERPINE1(HSP47), HSP70, HSP90alpha, Rano class II histocompatibility A beta, PGC-1 and ERR-induced regulator, proto-oncogene c-Fos, myozenin2, alpha-crystallin B chain-like protein, angiopoietin-like protein 8, and acetyl-CoA carboxylases have been identified in muscle tissue of rohu that are associated with stress/immunity. This study identified the key biomarker SERPINE1 (HSP47), which showed significant upregulation (~ 2- to threefold) in muscle tissue of rohu exposed to high temperature. This study can pave a path for the identification of stress-responsive biomarkers linked with thermal adaptations in the farmed carps.

Comparative analyses reveal potential genetic mechanisms for high-altitude adaptation of Schizopygopsis fishes based on chromosome-level genomes

The schizothoracine fishes, widely distributed in the Qinghai-Tibetan Plateau and its adjacent areas, are considered as ideal models for investigation of high-altitude adaptation. Schizophygopsis are one group of the highly specialized schizothoracine fishes, and the genetic basis for their high-altitude adaptation is poorly understood. In this study, we performed comparative genomics analyses to investigate the potential genetic mechanisms for high-altitude adaptation of Schizopygopsis malacanthus and Schizopygopsis pylzovi based on the chromosome-level genomes. Functional enrichment analysis revealed that many expanded gene families in Schizopygopsis were associated with immune response while many contracted gene families were functionally associated with olfaction. Among the 123 positively selected genes (PSGs), angpt2a was detected in HIF-1 signaling pathway and possibly related to the hypoxia adaptation of Schizopygopsis. Furthermore, two PSGs cox15 and ndufb10 were distributed in thermogenesis, and there was a Schizopygopsis-specific missense mutation in cox15 (Gln115Glu), which possibly contributed to the cold temperature adaptation of the Schizopygopsis. Kyoto Encyclopedia of Genes and Genomes enrichment of the PSGs revealed three significant pathways including metabolic pathways, cell cycle, and homologous recombination and Gene Ontology enrichment analysis of the PSGs revealed several categories associated with DNA repair, cellular response to DNA damage stimulus, and metabolic process. Chromosome-scale characterization of olfactory receptor (OR) repertoires indicated that Schizopygopsis had the least number of OR genes, and the OR gene contraction was possibly caused by the limited food variety and the environmental factors such as lower air pressure, lower humidity, and lower temperature. Our study will help expand our understanding of the potential adaptive mechanism of Schizopygopsis to cope with the high-altitude conditions.

Deterministic processes dominate microbial community assembly in artificially bred Schizothorax wangchiachii juveniles after being released into wild

Fish artificial breeding and release is an important method to restore wild populations of endemic fish species around the world. Schizothorax wangchiachii (SW) is an endemic fish in the upper Yangtze River and is one of the most important species for the artificial breeding and release program implemented in the Yalong River drainage system in China. It is unclear how artificially bred SW adapts to the changeable wild environment post-release, after being in a controlled and very different artificial environment. Thus, the gut samples were collected and analyzed for food composition and microbial 16S rRNA in artificially bred SW juveniles at day 0 (before release), 5, 10, 15, 20, 25, and 30 after release to the lower reaches of the Yalong River. The results indicated that SW began to ingest periphytic algae from the natural habitat before day 5, and this feeding habit is gradually stabilized at day 15. Prior to release, Fusobacteria are the dominant bacteria in the gut microbiota of SW, while Proteobacteria and Cyanobacteria generally are the dominant bacteria after release. The results of microbial assembly mechanisms illustrated that deterministic processes played a more prominent role than stochastic processes in the gut microbial community of artificially bred SW juveniles after releasing into the wild. Overall, the present study integrates the macroscopic and microscopic methods to provide an insight into the food and gut microbial reorganization in the released SW. This study will be an important research direction to explore the ecological adaptability of artificially bred fish after releasing into the wild.

Chromosome-level genome assembly and whole-genome resequencing of topmouth culter (Culter alburnus) provide insights into the intraspecific variation of its semi-buoyant and adhesive eggs

Topmouth culter (Culter alburnus) is an ecologically and economically important species belonging to the subfamily Culterinae that is native to and widespread in East Asia. Intraspecific variation of semi-buoyant and adhesive eggs in topmouth culter provides an ideal opportunity to investigate the genetic mechanisms of spawning habits underlying the adaptive radiation of cyprinids in East Asia. In this study, we present a chromosome-level genome assembly of topmouth culter and re-sequenced 158 individuals from six locations in China covering three geographical groups and two egg type variations. The topmouth culter genome size was 1.05 Gb, with a contig N50 length of 17.8 Mb and anchored onto 24 chromosomes. Phylogenetic analysis showed that the divergence time of the Culterinae was coinciding with the time of initiation of the Asian monsoon intensification. Gene family evolutionary analysis indicated that the expanded gene families in topmouth culter were associated with dietary adaptation. Population-level genetic analysis indicated clear differentiation among the six populations, which were clustered into three distinct clusters, consistent with their geographical divergence. The historical effective population size of topmouth culter correlated with the Tibetan Plateau uplifting according to the demographic history reconstruction. A selective sweep analysis between adhesive and semi-buoyant egg populations revealed the genes associated with the hydration and adhesiveness of eggs, indicating divergent selection towards different hydrological environments. This study offers a high-resolution genetic resource for further studies on evolutionary adaptation, genetic breeding and conservation of topmouth culter, providing insights into the molecular mechanisms for egg type variation of East Asian cyprinids.

Unveiling the molecular identity of the diminutive cyprinid, Horadandia brittani (Teleostei: Cyprinidae), a species endemic to Southern India

BACKGROUND

Horadandia brittani is a small cyprinid fish species initially discovered in the coastal floodplains of southern India. For almost 50 years, the genus Horadandia was monotypic with a single species confined to Sri Lanka. In 1992, a new species H. brittani was described from south-western India. Despite being described as a separate species, H. brittani was later considered a synonym of H. atukorali, but in 2013, researchers recognized it as a distinct species based on morphological differences. Despite this clarification, there was still a need to validate the identity of H. brittani and determine its evolutionary relationship with its closely related species using DNA sequences.

METHODS

To address the uncertainties surrounding the identity of H. brittani, the present study utilized molecular techniques to generate DNA sequences. Sample collection involved obtaining specimens of H. brittani from their natural habitats. Subsequently, DNA was extracted from the collected samples, and the mitochondrial cytochrome C oxidase (COI) gene was amplified using appropriate methods.

RESULTS

The analysis of DNA sequences obtained from the COI gene revealed significant genetic distinctions between H. brittani and H. atukorali. The genetic distance values between these two species ranged from 3.21 to 3.63%, clearly indicating that these two species are genetically separate entities. The study successfully established the phylogenetic relationships between H. brittani and H. atukorali based on the COI gene sequences, further confirming the validity of H. brittani as a distinct and separate species.

CONCLUSION

The findings of this study conclusively demonstrate that H. brittani is a valid and separate species, distinct from H. atukorali. The genetic analysis based on mitochondrial COI gene sequences provided strong evidence for the differentiation between these two species. The molecular data generated in this research can be used to identify H. brittani quickly and accurately in the future.

Biogeography and evolutionary history of Puntius sensu lato (Teleostei: Cyprinidae) in Sri Lanka

Sri Lanka's biota is derived largely from Southeast Asian lineages which immigrated via India following its early-Eocene contact with Laurasia. The island is now separated from southeastern India by the 30 km wide Palk Strait which, during sea-level low-stands, was bridged by the 140 km-wide Palk Isthmus. Consequently, biotic ingress and egress were mediated largely by the climate of the isthmus. Because of their dependence on perennial aquatic habitats, freshwater fish are useful models for biogeographic studies. Here we investigate the timing and dynamics of the colonization of-and diversification on-Sri Lanka by a group of four closely-related genera of cyprinid fishes (Puntius sensu lato). We construct a molecular phylogeny based on two mitochondrial and two nuclear gene markers, conduct divergence timing analyses and ancestral-range estimations to infer historical biogeography, and use haplotype networks to discern phylogeographic patterns. The origin of Puntius s.l. is dated to ~ 20 Ma. The source of diversification of Puntius s.l. is Sri Lanka-Peninsular India. Species confined to perhumid rainforests show strong phylogeographic structure, while habitat generalists show little or no such structure. Ancestral range estimations for Plesiopuntius bimaculatus and Puntius dorsalis support an 'Out of Sri Lanka' scenario. Sri Lankan Puntius s.l. derive from multiple migrations across the Palk Isthmus between the early Miocene and the late Pleistocene. Species dependent on an aseasonal climate survived aridification in rainforest refugia in the island's perhumid southwest and went on to recolonize the island and even southern India when pluvial conditions resumed. Our results support an historical extinction of Sri Lanka's montane aquatic fauna, followed by a recent partial recolonization of the highlands, showing also that headwater stream capture facilitated dispersal across basin boundaries.

Characterization and Expression Analysis of Genes from Megalobrama amblycephala Encoding Hemoglobins with Extracellular Microbicidal Activity

Hemoglobin (Hb) usually comprises two α and two β subunits, forming a tetramer responsible for oxygen transportation and storage. Few studies have elucidated fish hemoglobin immune functions. Megalobrama amblycephala is a freshwater-cultured fish prevalent in China. We identified two M. amblycephala hemoglobin subunits and analyzed their expression patterns and antibacterial activities. The respective full-length cDNA sequences of the M. amblycephala Hb α (MaHbα) and β (MaHbβ) subunits were 588 and 603 bp, encoding 143 and 148 amino acids. MaHbα and MaHbβ were highly homologous to hemoglobins from other fish, displaying typical globin-like domains, most heme-binding sites, and tetramer interface regions highly conserved in teleosts. In phylogenetic analyses, the hemoglobin genes from M. amblycephala and other cypriniformes clustered into one branch, and those from other fishes and mammals clustered into other branches, revealing fish hemoglobin conservation. These M. amblycephala Hb subunits exhibit different expression patterns in various tissues and during development. MaHbα is mainly expressed in the blood and brain, while MaHbβ gene expression is highest in the muscle. MaHbα expression was detectable and abundant post-fertilization, with levels fluctuating during the developmental stages. MaHbβ expression began at 3 dph and gradually increased. Expression of both M. amblycephala Hb subunits was down-regulated in most examined tissues and time points post-Aeromonas hydrophila infection, which might be due to red blood cell (RBC) and hematopoietic organ damage. Synthetic MaHbα and MaHbβ peptides showed excellent antimicrobial activities, which could inhibit survival and growth in five aquatic pathogens. Two M. amblycephala hemoglobin subunits were identified, and their expression patterns and antibacterial activities were analyzed, thereby providing a basis for the understanding of evolution and functions of fish hemoglobins.

Vascular Endothelial Growth Factor A (VEGFA) Regulates Hepatic Lipid and Glycogen Metabolism in Schizothorax prenanti

Vascular endothelial growth factor A (VEGFA) plays important roles in angiogenesis, inflammatory response as well as energy metabolism in mammals. However, its effect on glycolipid metabolism in fish has not been reported. In this study, we cloned and characterized the vegfa gene of Schizothorax prenanti (S. prenanti). vegfa expression was significantly higher in liver and muscle than that in other tissues. Then, the VEGFA recombinant protein was expressed in Escherichia coli and obtained after purification. VEGFA i.p. injection significantly increased the serum glucose and TG content compared with the control group. Moreover, VEGFA protein aggravated the glycogen and lipid deposition in the liver of S. prenanti. In addition, we found that VEGFA treatment increased hepatocyte glycogen and lipid droplet content and increased the levels of pAMPKα (T172). Furthermore, AMPKα inhibition attenuated the ability of VEGFA to induce TG and glycogen accumulation. These results demonstrate that VEGFA regulates hepatic lipid and glycogen metabolism through AMPKα in S. prenanti, which may contribute to a better understanding of VEGFA functions in the glycolipid metabolism of fish.

The Mitochondrial Genome of Linichthys laticeps (Cypriniformes: Cyprinidae): Characterization and Phylogeny

Mitochondrial genomes (mitogenomes) have been widely used in phylogenetic analysis and evolutionary biology. The Labeoninae is the largest subfamily of Cypriniformes and has great economic importance and ecological value. In this study, we sequenced, annotated, and characterized the complete mitogenome of Linichthys laticeps and then constructed the phylogenetic tree with previously published Labeoninae mitogenomes. The mitogenome of L. laticeps was 16,593 bp in length, with an A + T content of 57.1%. The mitogenome contained a standard set of 37 genes and a control region with the same order and orientation of genes as most fish mitogenomes. Each protein-coding gene (PCG) was initiated by an initial ATG codon, excluding COI, that began with a GTG codon. Furthermore, most of the PCGs were terminated by a conventional stop codon (TAA/TAG), while an incomplete termination codon (TA/T) was detected in 7 of the 13 PCGs. Most tRNA genes in L. laticeps were predicted to fold into the typical cloverleaf secondary structures. The Ka/Ks (ω) values for all PCGs were below one. The phylogenetic relationships of 96 Labeoninae mitogenomes indicated that Labeoninae was not a monophyletic group and L. laticeps was closely related to the genera Discogobio and Discocheilus. Overall, our study provided the first complete annotated mitogenome of L. laticeps, which filled a knowledge gap in Labeoninae and extended the understanding of the taxonomy and mitogenomic phylogeny of the subfamily Labeoninae.

Effects of hypoxia on the gill morphological structure, apoptosis and hypoxia-related gene expression in blunt snout bream (Megalobrama amblycephala)

The blunt snout bream (Megalobrama amblycephala) is a typical hypoxia-sensitive fish, and hypoxia stress leads to reduced vitality and yield during aquaculture. To explore the specific adaptation mechanism under hypoxia, the blunt snout bream was treated with hypoxia (DO = 2.0 ± 0.1 mg/L) for 24 h, followed by 3 h of recovery. Our results depicted that the gill filament structure of blunt snout bream changed after hypoxia. During hypoxia for 24 h, the gill filament structure was altered, including a more than 80% expansion of the lamellar respiratory surface area and a proportionate apoptosis decrease in interlamellar cell mass (ILCM) volume. Thus, the water-blood diffusion distance was shortened to less than 46%. During hypoxia for 24 h, the activity of ROS in gill tissue increased significantly (p < 0.05), while the mitochondrial membrane potential decreased significantly (p < 0.05). During hypoxia, mRNA expression level of anti-apoptotic gene Bcl-2 in the gills of blunt snout bream decreased significantly (p < 0.05), while the expression of pro-apoptotic gene Bax mRNA increased significantly (p < 0.05). Thus, the ratio of Bax/Bcl-2 mRNA increased in the gills of blunt snout bream to promote the activity of Caspase-3. Together, our results indicated hypoxia-induced apoptosis in the gills of blunt snout bream through the mitochondrial pathway. In addition, a decreased expression of Phd1 and an increased expression of Hif-1α in gills under hypoxia stress indicates that blunt snout bream may cope with hypoxia-induced apoptosis by enhancing the HIF pathway. These results provide new insights into fish's adaptation strategies and mechanisms of hypoxia.

Chromosome-level reference genomes of two imperiled desert fishes: spikedace (Meda fulgida) and loach minnow (Tiaroga cobitis)

North American minnows (Cypriniformes: Leuciscidae) comprise a diverse taxonomic group, but many members, particularly those inhabiting deserts, face elevated extinction risks. Despite conservation concerns, leuciscids remain under sampled for reference assemblies relative to other groups of freshwater fishes. Here, we present 2 chromosome-scale reference genome assemblies spikedace (Meda fulgida) and loach minnow (Tiaroga cobitis) using PacBio, Illumina and Omni-C technologies. The complete assembly for spikedace was 882.1 Mb in total length comprised of 83 scaffolds with N50 = 34.8 Mb, L50 = 11, N75 = 32.3 Mb, and L75 = 18. The complete assembly for loach minnow was 1.3 Gb in total length comprised of 550 scaffolds with N50 = 48.6 Mb, L50 = 13, N75 = 42.3 Mb, and L75 = 20. Completeness assessed via Benchmarking Universal Single-Copy Orthologues (BUSCO) metrics using the Actinopterygii BUSCO database showed ∼97% for spikedace and ∼98% for loach minnow of complete BUSCO proportions. Annotation revealed approximately 32.58 and 29.04% of spikedace and loach minnow total genome lengths to be comprised of protein-coding genes, respectively. Comparative genomic analyses of these endangered and co-distributed fishes revealed widespread structural variants, gene family expansions, and evidence of positive selection in both genomes.

Screening for IBs-relative genes by transcriptome analysis and generation IBs-less mutants in Culter alburnus

Intermuscular bones (IBs), distributed specifically in the myosepta on both sides of lower teleosts, negatively affect palatability and processing. Recent research in zebrafish and several economically important farmed fishes has led to the breakthrough discovery of the mechanism of IBs formation and generation of IBs-loss mutants. This study explored the ossification patterns of IBs in juvenile Culter alburnus. Besides, some key genes and bone-related signaling pathways were identified by transcriptomic data. Furthermore, PCR microarray validation revealed that claudin1 could potentially regulate IBs formation. Additionally, we created several IBs-reduced mutants of C. alburnus by loss of the function of bone morphogenetic proteins 6 (bmp6) gene using CRISPR/Cas9 editing. These results suggested that CRISPR/Cas9-mediated bmp6 knockout was promising approach for breeding IBs-free strain in other cyprinids.

Comprehensive transcripts analysis based on single-molecule real-time sequencing and Illumina sequencing provides insights into the mining of Toll-like receptor family in Schizothorax lissolabiatus

Schizothorax lissolabiatus is an economically important cold-water fish species in southwestern China. Because of water pollution and habitat destruction, the number of wild populations has dramatically decreased. In this study, we used PacBio single-molecule real-time (SMRT) sequencing and Illumina sequencing to generate the first full-length transcriptome and transcriptome, respectively. A total of 19 310 polished consensus reads (PC) were obtained, with an average length of 1379 bp and an N50 length of 1485 bp. Meanwhile, 12 253 transcripts were successfully annotated as known homologous genes. The pathway annotation indicated that the enrichment and expression of most genes were mainly related to membrane, signal transduction and binding, and immune response. Furthermore, we identified 16 Toll-like receptors (TLRs) by mining the data from the transcripts. Phylogeny analysis showed that S. lissolabiatus TLR genes (slTLRs) supported the classification of TLRs into six families as in other vertebrates. Selection pressure analyses showed that 16 slTLRs revealed purification selection at the overall evolutionary selection. Further, positive selection signals were still detected in eight slTLRs, and most of the positive selection sites were located in the leucine-rich repeat region (LRR domain) associated with the recognition of pathogenic microorganisms, indicating that the function of these slTLR genes may be affected. Tissue specific expression analysis showed all slTLRs are present in kidney, spleen and liver but the relative expression varied among tissues. In conclusion, this study not only provided a valuable resource of transcripts for further research on S. lissolabiatus, but also contributed to improve the current understanding of the evolutionary history of immune-related genes and the TLR gene family in S. lissolabiatus.

Changes in the liver of Tinca tinca under successive domestication using an integrated multi-omics approach

Domestication is the process of modifying the phenotype of a population through anthropic selection from human perspectives. Successive generations of domestication have influenced the physiological characteristics of tench Tinca tinca. In current study, we investigated gene and protein expression alterations in the liver of fifth-generation (F5). A total of 420 genes were found to be upregulated and 351 genes were downregulated, while 410 proteins were upregulated and 279 proteins were downregulated in domesticated T. tinca (DT). The integrated analysis of omics data revealed a total of 55 genes/proteins exhibiting consistent upregulation and 12 genes/proteins displaying consistent downregulation in DT. The upregulated genes/proteins in DT, such as SSR1, DERLIN2, OS9, DNAJB11, and HYOU1, exhibit enrichment in the protein processing in the endoplasmic reticulum pathway. Additionally, upregulated genes/proteins such as IL2RB, F13B, and IRF3 are associated with immune response. Conversely, downregulated genes/proteins in DT, including HSD11B1, CYP24A1, and COMT, play roles in hormone metabolism. These findings indicate that domestication can have a substantial impact on the physiological modifications related to protein processing, immune response, and hormone metabolism in DT. These adaptations potentially enhance their ability to thrive in artificial aquaculture environments, leading to improved growth and development. The exploration of genetic changes in DT will not only improve aquaculture practices but also provide significant insights into the broader process of domestication and its effects on physiological functions.

Apical and mechanistic effects of 6PPD-quinone on different life-stages of the fathead minnow (Pimephales promelas)

N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine-quinone (6PPD-quinone) is an emerging contaminant of concern that is generated through the environmental oxidation of the rubber tire anti-degradant 6PPD. Since the initial report of 6PPD-quinone being the cause of urban runoff mortality syndrome of Coho salmon, numerous species have been identified as either sensitive or insensitive to acute lethality caused by 6PPD-quinone. In sensitive species, acute lethality might be caused by uncoupling of mitochondrial respiration in gills. However, little is known about effects of 6PPD-quinone on insensitive species. Here we demonstrate that embryos of fathead minnows (Pimephales promelas) are insensitive to exposure to concentrations as great as 39.97 μg/L for 168 h, and adult fathead minnows are insensitive to exposure to concentrations as great as 9.4 μg/L for 96 h. A multi-omics approach using a targeted transcriptomics array, (EcoToxChips), and proton nuclear magnetic resonance (1H NMR) was used to assess responses of the transcriptomes and metabolomes of gills and livers from adult fathead minnows exposed to 6PPD-quinone for 96 h to begin to identify sublethal effects of 6PPD-quinone. There was little agreement between results of the EcoToxChip and metabolomics analyses, likely because genes present on the EcoToxChip were not representative of pathways suggested to be perturbed by metabolomic analysis. Changes in abundances of transcripts and metabolites in livers and gills suggest that disruption of one‑carbon metabolism and induction of oxidative stress might be occurring in gills and livers, but that tissues differ in their sensitivity or responsiveness to 6PPD-quinone. Overall, several pathways impacted by 6PPD-quinone were identified as candidates for future studies of potential sublethal effects of this chemical.

Red blood cells of Labeo rohita express Toll-like receptors, NOD- like receptors, interleukins, and interferon-I in response to Gram-negative bacterial infections and lipopolysaccharide stimulations

Red blood cells (RBCs) are the most abundant cell types in the circulatory system of vertebrates. In fish, RBCs retain their nuclei throughout their lifetime and remain transcriptionally and translationally active. While their primary function is typically associated with gas exchange, recent reports indicate that nucleated red blood cells can play a significant role in regulating the body's innate immune response. The current article describes the innate immune role of red blood cells in rohu (Labeo rohita), a freshwater fish species that holds significant commercial importance in India and South-East Asian nations. From the whole blood and mucosal surface RBCs have been isolated through density gradient centrifugation with HiSep™LSM 1077 (density 1.007 ± 0.0010) and their purity has been confirmed by the Giemsa staining followed by microscopical observations. Toll-like receptors (TLR2, 3, 4, 5) and nucleotide oligomerization domain (NOD)-like receptors (NOD1 and NOD2) in RBCs of rohu fingerlings were observed to be significantly activated (P < 0.05) on infection with Aeromonas hydrophila and Edwardsiella tarda. This activation resulted in increased expression of interleukins (IL-8, IL-1β) and interferon (IFN)-I genes. The activation of TLR4, NOD1 and NOD2, as well as the expression of interleukins and IFN-I genes have been observed in both in vivo and in vitro stimulation of rohu RBCs with lipopolysaccharides. These findings highlight the importance of fish RBCs in enhancing innate immunity against various pathogenic invasions in rohu.

Comparative transcriptomic analysis of reproductive characteristics of reciprocal hybrid lineages derived from hybridization between Megalobrama amblycephala and Culter alburnus

BACKGROUND

Distant hybridization is an important breeding technique for creating new strains with superior traits by integrating two different genomes. Successful hybridization of Megalobrama amblycephala (Blunt snout bream, BSB, 2n = 48) and Culter alburnus (Topmouth culter, TC, 2n = 48) was achieved to establish hybrid lineages (BT and TB), which provide valuable materials for exploring the mechanisms of distant hybridization fertility. In this study, the gonadal tissue transcriptomes of BSB, TC, BT-F1, and TB-F1 were sequenced using Illumina high-throughput sequencing technology to analyze the reproductive characteristics of BT and TB.

RESULTS

Differential gene expression analysis showed that the differentially expressed genes in BT vs BSB and BT vs TC were mainly enriched in signaling pathways not directly associated with meiosis. While, the differentially expressed genes of TB vs BSB and TB vs TC were mainly enriched in pathways related to meiosis, and most of them were down-regulated, indicating that meiosis is suppressed in TB. Under-dominance (UD) genes were enriched in pathways related to meiosis and DNA repair in TB. Over-dominance (OD) genes were enriched in MAPK signaling pathway, expression level dominance-BSB (ELD-B) genes were enriched in pathways related to steroid hormone synthesis and expression level dominance-TC (ELD-T) genes were not significantly enriched in any pathway in both BT and TB.

CONCLUSIONS

These results suggest that meiotic progression may not be affected in BT, whereas it is clearly inhibited in TB. Offspring of M. amblycephala maternal parent may have better genomic compatibility and fertility. Our study provides important information on the molecular mechanisms of breaking reproductive isolation in distantly hybridized fertile lineages.

[Mitochondrial genome sequence characteristics and phylogenetic analysis of Schizothorax argentatus]

Schizothorax argentatus that only distributes in the Ili River basin in Xinjiang is one of the rare and endangered species of schizothorax in China, thus has high scientific and economic values. In this study, the complete mitochondrial genome sequence of S. argenteus with a length of 16 580 bp was obtained by high-throughput sequencing. The gene compositions and arrangement were similar to those of typical vertebrates. It contained 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a non-coding region (D-loop). The nucleotide compositions were A (30.25%), G (17.28%), C (27.20%), and T (25.27%), respectively, showing obvious AT bias and anti-G bias. Among the tRNA genes, only tRNA-Ser(GCU) could not form a typical cloverleaf structure due to the lack of dihydrouracil arm. The AT-skew and GC-skew values of the ND6 gene were fluctuating the most, suggesting that the gene may experience different selection and mutation pressures from other genes. The mitochondrial control region of S. argenteus contained three different domains, i.e., termination sequence region (ETAS), central conserved region (CSB-F, CSB-E, CSB-D, and CSB-B), and conserved sequence region (CSB1, CSB2, and CSB3). The conserved sequence fragment TT (AT) nGTG, which was ubiquitous in Cypriniformes, was identified at about 50 bp downstream CSB3. Phylogenetic relationships based on the complete mitochondrial genome sequence of 28 Schizothorax species showed that S. argenteus had differentiated earlier and had a distant relationship with other species, which may be closely related to the geographical location and the hydrological environment where it lives.

Metagenomic analysis reveals hidden links between gut microbes and habitat adaptation among cave and surface dwelling Sinocyclocheilus species

Intestinal microbes are closely related to vital host functions such as digestion and nutrient absorption, which play important roles in enhancing host adaptability. As a natural "laboratory", caves provide an outstanding model for understanding the significance of gut microbes and feeding habits in the habitat adaptability of hosts. However, research on the relationship between gut microbes, feeding habits, and the adaptability of troglobites remains insufficient. In this study, we compared the characteristics of the intestinal microbes of Sinocyclocheilus cavefish and surface fish and further established the relationship between intestinal and habitat microbes. Furthermore, we conducted environmental DNA (eDNA) (metabarcoding) analysis of environmental samples to clarify the composition of potential food resources in the habitats of the Sinocyclocheilus cavefish and surface fish. Results showed that the structure of the Sinocyclocheilus gut microbes was more related to ecological type (habitat type) than phylogenetic relationships. While horizontal transfer of habitat microbes was a source of gut microbes, hosts also showed strong selection for inherent microbes as dominant microorganisms. Differences in the composition and structure of gut microbes, especially dominant microbes, may enhance the adaptability of the two Sinocyclocheilus fish types from the perspectives of food intake, nutrient utilization, and harmful substance metabolism, suggesting that food resources, predation patterns, intestinal flora, digestive and absorptive capacity, and feeding habits and preferences are linked to habitat adaptability. These results should facilitate our understanding of the significance of fish gut microbes to habitat adaptation and provide a new perspective for studying the adaptive mechanisms of cavefish.

Population structure and genetic diversity of the endangered fish black shinner Pseudopungtungia nigra (Cyprinidae) in Korea: a wild and restoration population

The black shinner Pseudopungtungia nigra Mori, 1935 is an endangered fish endemic to Korea. It lives in the narrow basin of the Geumgang River, Mangyeonggang River, and Ungcheoncheon Stream, which flow into the West Sea of Korea. One population of P. nigra in Ungcheoncheon Stream has been locally exterminated once; it is now inhabiting the upper reaches of the dam through a restoration program. Efforts to identify and understand the genetic structure of these populations are important for conservation planning. Here, we analyzed genetic diversity using 21 microsatellite markers for 9 populations. The mean number of alleles ranged from 4.4 to 8.1, mean allelic richness ranged from 4.6 to 7.8, mean observed heterozygosity ranged from 0.519 to 0.702, and mean expected heterozygosity ranged from 0.540 to 0.763. All groups had recent and historical bottlenecks (P < 0.05, M-ratio < 0.68). Three groups [YD (2019), OC and UC] had significant inbreeding index values, suggesting that they were engaged in inbreeding. We observed a moderate level of genetic differentiation between MG and the rest of the population (FST = 0.135 to 0.168, P < 0.05). The genetic structure exhibited a fitting constant K = 2, along with separation between MG and the remaining populations. With respect to genetic flow, YD (2019), OC, CG, and ND shifted to the UC population (0.263 to 0.278). The genetic flow of each population was transferred only within the population; there was no gene flow among populations, except for the Ungcheoncheon Stream population. This study shows that the Ungcheoncheon Stream population needs conservation efforts to increase its genetic diversity, and the Geumgang River populations needs a conservation plan that considers the possibility of conservation and evolution through gene exchange among the populations.

An efficient approach to synthesize sterile allopolyploids through the combined reproduction mode of ameiotic oogenesis and sperm-egg fusion in the polyploid Carassius complex

The association between polyploidy and reproduction transition, which is an intriguing issue in evolutionary genetics, can also be exploited as an approach for genetic improvement in agriculture. Recently, we generated novel amphitriploids (NA3n) by integrating the genomes of the gynogenetic Carassius gibelio and sexual C. auratus, and found gynogenesis was recovered in most NA3n females (NA3n♀I). Here, we discovered a unique reproduction mode, termed ameio-fusiongenesis, which combines the abilities of both ameiotic oogenesis and sperm-egg fusion, in a few NA3n females (NA3n♀II). These females inherited ameiotic oogenesis to produce unreduced eggs from gynogenetic C. gibelio and sperm-egg fusion from sexual C. auratus. Subsequently, we utilized this unique reproduction mode to generate a group of synthetic alloheptaploids by crossing NA3n♀II with Megalobrama amblycephala. They contained all chromosomes of maternal NA3n♀II and a chromosomal set of paternal M. amblycephala. Intergenomic chromosome translocations between NA3n♀II and M. amblycephala were also observed in a few somatic cells. Primary oocytes of the alloheptaploid underwent severe apoptosis owing to incomplete double-strand break repair at prophase I. Although spermatocytes displayed similar chromosome behavior at prophase I, they underwent apoptosis due to chromosome separation failure at metaphase I. Therefore, the alloheptaploid females and males were all sterile. Finally, we established a sustainable clone for the large-scale production of NA3n♀II and developed an efficient approach to synthesize diverse allopolyploids containing genomes of different cyprinid species. These findings not only broaden our understanding of reproduction transition but also offer a practical strategy for polyploidy breeding and heterosis fixing.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.