Effects of fasting and environmental factors on appetite regulators in pond loach Misgurnus anguillicaudatus

The pond loach (Misgurnus anguillicaudatus) is an important aquaculture freshwater species, used as an ornamental fish, food source for humans and angling bait. Pond loaches are resistant to fasting and extreme environmental conditions, including temperature and low oxygen levels. Little is known about how these factors affect the feeding physiology and the endocrine regulation of feeding of loaches. In this study, we examined the effects of fasting, as well as increased temperature and decreased oxygen levels on food intake and transcript levels of appetite regulators. Fasted fish had lower blood glucose levels, and lower expression levels of intestine CCK and PYY, and brain CART1, but had higher levels of brain orexin and ghrelin than fed fish. Fish held at 30 °C had higher food intake, glucose levels, and mRNA levels of intestine CCK and PYY, and brain CART2, but lower brain orexin levels than fish at 20 °C. Fish held at low oxygen levels had a lower food intake, higher intestine CCKa and ghrelin, and brain orexin, CART2 and ghrelin mRNA expression levels than fish held at high O2 levels. Our results suggest that fasting and high temperatures increase the expression of orexigenic and anorexigenic factors respectively, whereas the increase in expression of both orexigenic and anorexigenic factors in low O2 environments might not be related to their role in feeding, but possibly to protection from tissue damage. The results of our study might shed new light on how pond loaches are able to cope with extreme environmental conditions such as low food availability, extreme temperatures and hypoxia.

Genetic characterization, structural analysis, and detection of positive selection in small heat shock proteins of Cypriniformes and Clupeiformes

In the current investigation, a total of 42 full-length, non-redundant small heat shock proteins (sHsp) were detected in Cyprinus carpio, Labeo rohita, Danio rerio, Salmo salar, Oncorhynchus mykiss, and Clupea harengus. The sHsp genes were classified into three groups based on phylogenetic analysis. All the sHsps were shown to have higher aliphatic index values, which is an indication that these proteins are more thermally stable. The hydrophilic nature of sHsps was deduced from the fact that all fish species had negative GRAVY scores. In all of the representative fish species, sHsp genes were assigned to distinct chromosomes in an inconsistent and unequal manner. Segmental duplications are the main events that have contributed to the expansion of the sHsp genes in all species. We were also able to determine the selective pressure that was placed on particular codons and discovered several significant coding sites within the coding region of sHsps. Eventually, diversifying positive selection was found to be connected with evolutionary changes in sHsp proteins, which showed that gene evolution controlled the fish adaption event in response to environmental conditions. Clarification of the links between sHsps and environmental stress in fish will be achieved through rigorous genomic comparison, which will also yield substantial new insights.

Environmental DNA unveils deep phylogeographic structure of a freshwater fish

Phylogeography bears an important part in ecology and evolution. However, current phylogeographic studies are largely constrained by limited numbers of individual samples. Using an environmental DNA (eDNA) assay for phylogeographic analyses, this study provides detailed information regarding the history of Siberian stone loach Barbatula toni, a primary freshwater fish across the whole range of Hokkaido, Japan. Based on an eDNA metabarcoding on 293 river water samples, we detected eDNA from B. toni in 189 rivers. A total of 51 samples, representing the entire island, were then selected from the B. toni eDNA-positive sample set for the subsequent analyses. To elucidate the phylogeographic structure of B. toni, newly developed eDNA metabarcoding primers (Barba-cytb-F/R) were applied to these samples, specifically targeting their haplotypic variation in cytochrome b. After a bioinformatic processing to mitigate haplotypic false positives, a total of 50 eDNA haplotypes were identified. Two regionally restricted, genetically distinct lineages of the species were revealed as a result of phylogeographic analyses on the haplotypes and tissue-derived DNA from B. toni. According to a molecular clock analysis, they have been genetically isolated for at least 1.5 million years, suggesting their ancient origin and colonisation of Hokkaido, presumably in the glacial periods. These results demonstrate how freshwater fishes can alter their distributions over evolutionary timescales and how eDNA assay can deepen our understanding of phylogeography.

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.

A cyclical switch of gametogenic pathways in hybrids depends on the ploidy level

The cellular and molecular mechanisms governing sexual reproduction are conserved across eukaryotes. Nevertheless, hybridization can disrupt these mechanisms, leading to asexual reproduction, often accompanied by polyploidy. In this study, we investigate how ploidy level and ratio of parental genomes in hybrids affect their reproductive mode. We analyze the gametogenesis of sexual species and their diploid and triploid hybrids from the freshwater fish family Cobitidae, using newly developed cytogenetic markers. We find that diploid hybrid females possess oogonia and oocytes with original (diploid) and duplicated (tetraploid) ploidy. Diploid oocytes cannot progress beyond pachytene due to aberrant pairing. However, tetraploid oocytes, which emerge after premeiotic genome endoreplication, exhibit normal pairing and result in diploid gametes. Triploid hybrid females possess diploid, triploid, and haploid oogonia and oocytes. Triploid and haploid oocytes cannot progress beyond pachytene checkpoint due to aberrant chromosome pairing, while diploid oocytes have normal pairing in meiosis, resulting in haploid gametes. Diploid oocytes emerge after premeiotic elimination of a single-copied genome. Triploid hybrid males are sterile due to aberrant pairing and the failure of chromosomal segregation during meiotic divisions. Thus, changes in ploidy and genome dosage may lead to cyclical alteration of gametogenic pathways in hybrids.

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.

lncRNA-miRNA-mRNA network in kidney transcriptome of Labeo rohita under hypersaline environment

The present study describes the kidney transcriptome of Labeo rohita, a freshwater fish, exposed to gradually increased salinity concentrations (2, 4, 6 and 8ppt). A total of 10.25 Gbps data was generated, and a suite of bioinformatics tools, including FEELnc, CPC2 and BLASTn were employed for identification of long non-coding RNAs (lncRNAs) and micro RNAs (miRNAs). Our analysis revealed a total of 170, 118, 99, and 269 differentially expressed lncRNA and 120, 118, 99, and 124 differentially expressed miRNAs in 2, 4, 6 and 8 ppt treatment groups respectively. Two competing endogenous RNA (ceRNA) networks were constructed i.e. A* ceRNA network with up-regulated lncRNAs and mRNAs, down-regulated miRNAs; and B* ceRNA network vice versa. 2ppt group had 131 and 83 lncRNA-miRNA-mRNA pairs in A* and B* networks, respectively. 4ppt group featured 163 pairs in A* network and 191 in B* network, while the 6ppt had 103 and 105 pairs. 8ppt group included 192 and 174 pairs. These networks illuminate the intricate RNA interactions in freshwater fish to varying salinity conditions.

Genome-wide identification and structural analysis of the BMP gene family in Triplophysa dalaica

BACKGROUND

Bone morphogenetic proteins (BMPs) are part of the transforming growth factor beta (TGF-β) superfamily and play crucial roles in bone development, as well as in the formation and maintenance of various organs. Triplophysa dalaica, a small loach fish that primarily inhabits relatively high elevations and cooler water bodies, was the focus of this study. Understanding the function of BMP genes during the morphogenesis of T. dalaica helps to clarify the mechanisms of its evolution and serves as a reference for the study of BMP genes in other bony fishes. The data for the T. dalaica transcriptome and genome used in this investigation were derived from the outcomes of our laboratory sequencing.

RESULTS

This study identified a total of 26 BMP genes, all of which, except for BMP1, possess similar TGF-β structural domains. We conducted an analysis of these 26 BMP genes, examining their physicochemical properties, subcellular localization, phylogenetic relationships, covariance within and among species, chromosomal localization, gene structure, conserved motifs, conserved structural domains, and expression patterns. Our findings indicated that three BMP genes were associated with unstable proteins, while 11 BMP genes were located within the extracellular matrix. Furthermore, some BMP genes were duplicated, with the majority being enriched in the GO:0008083 pathway, which is related to growth factor activity. It was hypothesized that genes within the BMP1/3/11/15 subgroup (Group I) play a significant role in the growth and development of T. dalaica. By analyzing the expression patterns of proteins in nine tissues (gonad, kidney, gill, spleen, brain, liver, fin, heart, and muscle), we found that BMP genes play diverse regulatory roles during different stages of growth and development and exhibit characteristics of division of labor.

CONCLUSIONS

This study contributes to a deeper understanding of BMP gene family member expression patterns in high-altitude, high-salinity environments and provides valuable insights for future research on the BMP gene family in bony fishes.

Across the highest mountain on earth: discordant phylogeographic patterns and recent dispersal of Tibetan stone loaches (Triplophysa) in the Himalayas

Phylogeographic congruence among co-distributed taxa is regarded as an inherent inference to vicariance events. Nonetheless, incongruent patterns of contemporary lineage divergence among taxa indicated that species differ in their response to common past events. To investigate the role of past events, ecological traits and lineage diversification time in shaping the contemporary phylogeographic patterns, comparative analyses were conducted for Tibetan stone loaches in the Himalayas using three gene markers and two ecological traits (depth of caudal peduncle in their length and presence/absence of posterior chamber of the air bladder). By a thorough sampling in two flanks of the Himalayas, the authors detected that phylogenetic breaks were spatially discordant and divergences of populations were also temporally asynchronous in co-distributed loaches. Estimated divergence time using fossil-calibrated node dating indicated that the Tibetan stone loaches colonised into the south flank of the Himalayas until the Pleistocene. The demographic expansions were also disconcerted between populations in north and south flanks, or east and west Himalayas. Ongoing gene flows between populations in north and south sides implied that the Himalayas do not strictly impede dispersal of cold-adapted species. The results highlight that the quaternary climatic oscillation, in conjunction with ecological traits and lineage diversification time, shaped contemporary phylogenetic patterns of stone loaches in the Himalayas and provide new insights into the biodiversity and composition of species in the Himalayas and surrounding region.

Comparative transcriptome profile reveals insight into the antibacterial immunity mechanism of the loach (Misgurnus anguillicaudatus) fed with soybean fermented broth during lipopolysaccharide (LPS) exposure

Loach (Misgurnus anguillicaudatus) is a common freshwater commercial fish species in China. The meat of this fish is a good source of protein and other nutrients that are needed for human health. Aquaculture challenges such as diseases and pest susceptibility, excessive density, and nutritional deficiency result in low production of loach rather than increased demand. Due to a lack of knowledge about the immune system of loaches, we carried out this study to better understand its antibacterial molecular mechanism. Here, we performed RNA sequencing from liver tissue obtained from soya bean-fermented fed loach after subjecting it to the LPS challenge. The results revealed a total of 18,399 differentially expressed genes (DEGs) in the LPS-treated and control groups. There were 7482 DEGs that were upregulated and 10,917 DEGs were downregulated. The enrichment analysis of DEGs revealed that the majority of DEGs were found to be abundant in the pathways of DNA replication, spliceosome, nucleotide exception repair, cell cycle, and Herpes simplex virus 1 infection. Furthermore, qRT-PCR analysis of 21 selected DEGs demonstrated that the transcriptomic data is extremely reliable. Overall, this study provides insight into the molecular features and control mechanisms of genes that affect loach growth. The availability of this information will also contribute to the enhancement of the breeding and protection of loach resources.

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%.

Adaptive evolution of mitochondrial genomes in Triplophysa cavefishes

Extreme conditions in caves pose survival challenges for cave dwellers, who gradually develop adaptive survival features. Cavefishes are one of the most successful animals among cave dwellers. Triplophysa cavefishes are an important group of cavefishes, and they show remarkable adaptability to the extreme environments of caves. However, there is a limited understanding of their adaptation mechanisms. In this study, eight complete mitochondrial genomes of Triplophysa cavefishes were newly obtained, and their genomic characteristics, including the base composition, base bias, and codon usage, were analyzed. Phylogenetic analysis was carried out based on 13 mitochondrial protein-coding genes from 44 Nemacheilidae species. This showed that Triplophysa cavefishes and non-cavefishes separate into two reciprocally monophyletic clades, suggesting a single origin of the cave phenotype. Positive selection analysis strongly suggested that the selection pressure in cavefishes is higher than that in non-cavefishes. Furthermore, the ND5 gene in cavefishes showed evidence of positive selection, which suggests that the gene may play an important role in the adaptation of cavefishes to the cave environment. Protein structure analysis of the ND5 subunit implied that the sites of positive selection in cavefishes might allow them to acquire lower ND5 protein stability, compared to that in non-cavefishes, which might help the accumulation of nonsynonymous (mildly deleterious) mutations. Together, our study revealed the genetic signatures of cave adaptation in Triplophysa cavefishes from the perspective of energy metabolism.

Cloning, tissue distribution, and effects of different circadian rhythms on the mRNA expression levels of circadian clock genes Per1a and Per1b in Phoxinus lagowskii

This study describes the cloning and characterization of Period 1a and Period 1b genes and the analysis of their mRNA and protein expression in Amur minnow (Phoxinus lagowskii) after exposure to different light cycles. The full-length P. lagowskii Per1a and Per1b genes encode proteins consisting of 1393 and 1409 amino acids, and share high homology with the per1 genes of other freshwater fish species. The Per1a and Per1b genes were widely expressed within the brain, eye, and peripheral tissues. The acrophase of the Per1a gene in the pituitary gland occurred during the dark phase at ZT15 (zeitgeber time 15, 12 L: 12 D) and ZT18 (8 L, 16 D), whereas the acrophase of the Per1b gene in the pituitary gland was observed during the light phase. Our study suggests that the expression of Per1a and Per1b in P. lagowskii varied depending on differences in circadian rhythm patterns. The results of our dual-luciferase reporter assays demonstrated that the P. lagowskii Per1b gene enhances the activation of NF-κB. This study is the first to examine the circadian clock gene Per1a and Per1b in the high-latitude fish P. lagowskii, offering valuable insights into the effects of different light periods on this fish species.

Contrasting genetic and morphological differentiation among geographical lineages of a stenotopic miniature rasborine, Boraras maculatus, in Peninsular Malaysia

The variability in the stenotopic miniature rasborine Boraras maculatus (Cypriniformes: Danionidae: Rasborinae) across acidic-water habitats of Peninsular Malaysia (PM) was investigated using two molecular markers (the mitochondrial cytochrome c oxidase subunit I [COI] gene and the nuclear rhodopsin gene), as well as morphological evidence. Molecular phylogenetic analyses revealed differentiation among populations of B. maculatus in PM with the distinction of four allopatric lineages. Each of them was recognized as a putative species by automatic species delimitation methods. These lineages diverged from each other between 7.4 and 1.9 million years ago. A principal component analysis (PCA) was conducted to examine the multivariate variation in 11 morphometric measurements among three of these lineages. PCA results showed a significant overlap in morphological characteristics among these lineages. Additionally, a photograph-based machine learning approach failed to fully differentiate these lineages, suggesting limited morphological differentiation. B. maculatus represents a case of morphological stasis in a stenotopic miniature species. Strong habitat preference, coupled with long-term habitat fragmentation, may explain why each lineage of B. maculatus has a restricted distribution and did not disperse to other regions within and outside of PM, despite ample possibilities when the Sunda shelf was emerged and drained by large paleodrainages for most of the past 7 million years. The conservation status of B. maculatus and its peat swamp habitats are discussed, and it is concluded that peat swamps comprise several evolutionary units. Each of these units is considered a conservation unit and deserves appropriate protection.

Complete mitochondrial genome and its phylogenetic implications of Rhinogobio nasutus, an endemic species from the Yellow River

BACKGROUND

Rhinogobio nasutus, an endemic species from the Yellow River, is listed under the second class of the National Key Protected Wildlife List in China due to its dramatically decreased population. Despite its important status, the mitochondrial genes and phylogenetic relationships of R. nasutus are unknown.

METHODS AND RESULTS

The complete mitochondrial genome of R. nasutus was sequenced, assembled, and annotated for the first time. The mitochondrial genome was 16,609 bp in length, consisting of 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and 1 non-coding control region. The gene order in the mitochondrial genome of R. nasutus was identical to that of other Rhinogobios species. Analysis of synonymous and non-synonymous nucleotide substitutions showed that the Ka/Ks ratio in all tested protein-coding genes was less than 1, indicating that these genes were evolving under purifying selection. Further phylogenetic analysis showed that R. nasutus was first clustered with R. typus, then grouped with the other two Rhinogobio species, indicating the phylogenetically close relationship between R. nasutus and R. typus.

CONCLUSIONS

This was the first genomic resource developed for R. nasutus, which could not only improve our understanding of its phylogenetic status, but also serve as a genomic tool for the development of genetic markers that will be used in conservation and evolutionary genetics studies.

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

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

The New Mitochondrial Genome of Hemiculterella wui (Cypriniformes, Xenocyprididae): Sequence, Structure, and Phylogenetic Analyses

Hemiculterella wui is an endemic small freshwater fish, distributed in the Pearl River system and Qiantang River, China. In this study, we identified and annotated the complete mitochondrial genome sequence of H. wui. The mitochondrial genome was 16,619 bp in length and contained 13 protein coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region. The nucleotide composition of the mitochondrial genome was 29.9% A, 25.3% T, 27.4% C, and 17.5% G, respectively. Most PCGs used the ATG start codon, except COI and ATPase 8 started with the GTG start codon. Five PCGs used the TAA termination codon and ATPase 8 ended with the TAG stop codon, and the remaining seven genes used two incomplete stop codons (T and TA). Most of the tRNA genes showed classical cloverleaf secondary structures, except that tRNASer(AGY) lacked the dihydrouracil loop. The average Ka/Ks value of the ATPase 8 gene was the highest, while the average Ka/Ks value of the COI gene was the lowest. Phylogenetic analyses showed that H. wui has a very close relationship with Pseudohemiculter dispar and H. sauvagei. This study will provide a valuable basis for further studies of taxonomy and phylogenetic analyses in H. wui and Xenocyprididae.

Integrated metabolomic and transcriptomic responses to heat stress in a high-altitude fish, Triplophysa siluroides

Species in Triplophysa display strong adaptability to the extreme environment of the plateau, thus offering an ideal model to study the molecular mechanism of fish adaptation to environmental stress. In the present study, we conducted integrated analysis of the transcriptome and metabolism of liver tissue in Triplophysa siluroides under heat stress (28 °C) and control (10 °C) conditions to identify heat stress-induced genes, metabolites and pathways. RNA-Seq identified 2373 differentially expressed genes, which consisted of 1360 upregulated genes and 1013 downregulated genes, in the heat stress group vs. the control group. Genes in the heat shock protein (Hsp) family, including Hsp40, Hsp70, Hsp90 and other Hsps, were strongly upregulated by heat stress. Pathway enrichment analysis revealed that the PI3K/AKT/mTOR and protein processing in the endoplasmic reticulum (ER) pathways were significantly affected by heat stress. Metabolism sequencing identified a total of 155 differentially abundant metabolites, including 118 significantly upregulated metabolites and 37 downregulated metabolites. Combined analysis of the transcriptome and metabolism results showed that ubiquitin-dependent proteolysis and purine metabolism pathways were enhanced in response to acute heat stress to protect cells from damage under stress conditions. The results of this study may contribute to our understanding of the underlying molecular mechanism of the heat stress response in cold-water fish.

Oxynoemacheilus marmaraensis, a new species from the Susurluk River, Türkiye (Teleostei: Nemacheilidae)

Oxynoemacheilus marmaraensis, new species, is restricted to the Susurluk River. It is distinguished from all the named species of Oxynoemacheilus in the northwestern Anatolian by the flank with a vermiculate pattern and the presence of a suborbital groove in males, and no axillary lobe at the base of the pelvic fin. It also differs from the closest species, Oxynoemacheilus kentritensis, by having 58 nucleotide substitution sites. The genetic distance is 10.49% between O. marmaraensis and O. kentritensis. Phylogenetic analyses and species delimitation tests (Poisson tree processes and assemble species by automatic partitioning) support the validity of O. marmaraensis as a distinct species.

Environmental DNA metabarcoding of water samples as a tool for monitoring Iberian freshwater fish composition

Environmental DNA (eDNA) metabarcoding has been increasingly used to monitor the community assemblages of a wide variety of organisms. Here, we test the efficacy of eDNA metabarcoding to assess the composition of Iberian freshwater fishes, one of the most endangered groups of vertebrates in Spain. For this purpose, we sampled 12 sampling sites throughout one of Spain's largest basins, the Duero, which is home to approximately 70% of the genera and 30% of the primary freshwater fish in Spain. We sampled these sampling sites in the summer by using electrofishing, a traditional sampling method, and eDNA metabarcoding of river water samples using the mitochondrial 12S rRNA gene (12S) as a marker. We also resampled four of these sampling sites in autumn by eDNA. We compared the results obtained through eDNA metabarcoding with those of electrofishing surveys (ones conducted for the present study and past ones) and assessed the suitability of 12S as an eDNA metabarcoding marker for this group of freshwater fishes. We found that the 12S fragment, analysed for 25 Iberian species, showed sufficient taxonomic resolution to be useful for eDNA approaches, and even showed population-level differences in the studied populations across the tissue samples for Achondrostoma arcasii. In most cases, a greater number of species was detected through eDNA metabarcoding than through electrofishing. Based on our results, eDNA metabarcoding is a powerful tool to study the freshwater fish composition in the Iberian Peninsula and to unmask cryptic diversity. However, we highlight the need to generate a local genetic database for 12S gene for such studies and to interpret the results with caution when studying only mitochondrial DNA. Finally, our survey shows that the high detection sensitivity of eDNA metabarcoding and the non-invasiveness of this method allows it to act as a detection system for species of low abundance, such as early invasive species or species in population decline, two key aspects of conservation management of Spanish freshwater fishes.

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.

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.