The loach haplotype-resolved genome and the identification of Mex3a involved in fish air breathing.

Fish air breathing is crucial for the transition of vertebrates from water to land. So far, the genes involved in fish air breathing have not been well identified. Here, we performed gene enrichment analysis of positively selected genes (PSGs) in loach (Misgurnus anguillicaudatus, an air-breathing fish) in comparison to Triplophysa tibetana (a non-air-breathing fish), haplotype-resolved genome assembly of the loach, and gene evolutionary analysis of air-breathing and non-air-breathing fishes and found that the PSG mex3a originated from ancient air-breathing fish species. Deletion of Mex3a impaired loach air-breathing capacity by inhibiting angiogenesis through its interaction with T-box transcription factor 20. Mex3a overexpression significantly promoted angiogenesis. Structural analysis and point mutation revealed the critical role of the 201st amino acid in loach Mex3a for angiogenesis. Our findings innovatively indicate that the ancient mex3a is a fish air-breathing gene, which holds significance for understanding fish air breathing and provides a valuable resource for cultivating hypoxia-tolerant fish varieties.

Chromosome-scale and haplotype-resolved genome assembly of the autotetraploid Misgurnus anguillicaudatus.

In nature, diploids and tetraploids are two common types of polyploid evolution. Misgurnus anguillicaudatus (mud loach) is a remarkable fish species that exhibits both diploid and tetraploid forms. However, reconstructing the four haplotypes of its autotetraploid genome remains unresolved. Here, we generated the first haplotype-resolved, chromosome-level genome of autotetraploid M. anguillicaudatus with a size of 4.76 Gb, contig N50 of 6.78 Mb, and scaffold N50 of 44.11 Mb. We identified approximately 2.9 Gb (61.03% of genome) of repetitive sequences and predicted 91,485 protein-coding genes. Moreover, allelic gene expression levels indicated the absence of significant dominant haplotypes within the autotetraploid loach genome. This genome will provide a valuable biological model for unraveling the mechanisms of polyploid formation and evolution, adaptation to environmental changes, and benefit for aquaculture applications and biodiversity conservation.

Chromosome-level genome assemblies of vulnerable male and female elongate loach (Leptobotia elongata).

Endemic to the upper and middle reaches of the Yangtze River in China, elongate loach (Leptobotia elongata) has become a vulnerable species mainly due to overfishing and habitat destruction. Thus far, no genome data of this species are reported. As a result, lacking of such genomic information has restricted practical conservation and utilization of this economic fish. Here, we constructed chromosome-level genome assemblies for both male and female elongate loach by integration of MGI, PacBio HiFi and Hi-C sequencing technologies. Two primary genome assemblies (586-Mb and 589-Mb) were obtained for female and male fishes, respectively. Indeed, 98.22% and 98.61% of the contig sequences were anchored onto 25 chromosomes, with identification of 26.22% and 25.92% repeat contents in both assembled genomes. Meanwhile, a total of 25,215 and 25,253 protein-coding genes were annotated, of which 97.41% and 98.8% could be predicted with functions. Taken together, our genome data presented here provide a valuable genomic resource for in-depth evolutionary and functional research, as well as molecular breeding and conservation of this economic fish species.

Ecological genetics of isolated loach populations indicate compromised adaptive potential.

Many endangered species live in fragmented and isolated populations with low genetic variability, signs of inbreeding, and small effective population sizes - all features elevating their extinction risk. The flat-headed loach (Oreonectes platycephalus), a small noemacheilid fish, is widely across southern China, but only in the headwaters of hillstreams; as a result, they are spatially isolated from conspecific populations. We surveyed single nucleotide polymorphisms in 16 Hong Kong populations of O. platycephalus to determine whether loach populations from different streams were genetically isolated from each other, showed low levels of genetic diversity, signs of inbreeding, and had small contemporary effective population sizes. Estimates of average observed heterozygosity (HO = 0.0473), average weighted nucleotide diversity (πw = 0.0546) and contemporary effective population sizes (Ne = 10.2 ~ 129.8) were very low, and several populations showed clear signs of inbreeding as judged from relatedness estimates. The degree of genetic differentiation among populations was very high (average FST = 0.668), even over short geographic distances (<1.5 km), with clear patterns of isolation by distance. These results suggest that Hong Kong populations of O. platycephalus have experienced strong genetic drift and loss of genetic variability because sea-level rise after the last glaciation reduced connectedness among paleodrainages, isolating populations in headwaters. All this, together with the fact that the levels of genetic diversity and contemporary effective population sizes within O. platycephalus populations are lower than most other freshwater fishes, suggests that they face high local extinction risk and have limited capacity for future adaptation.

Schistura sonarengaensis, a new species of cave-dwelling loach (Teleostei: Nemacheilidae) from Meghalaya, northeast India.

A new species of nemacheilid loach, Schistura sonarengaensis sp. nov., is described from three cave-dwelling populations (Barak-Surma-Meghna drainage) in the South Garo Hills district of Meghalaya, India. The new species possesses prominent eyes but is easily distinguished from all the congeners of the genus Schistura from Barak-Surma-Meghna and adjacent rivers drainages of northeast Indian (except S. syngkai) in having 13-26 vertically elongated to circular mid-lateral black blotches (brownish in life) overlayered on a grayish-black mid-lateral stripe on a dull white or pale-beige (golden brown in life) body. However, it can be easily distinguished from S. syngkai in having a complete (vs. incomplete) lateral line with more 72-89 (vs. 19-42) lateral-line pored scales, greater pre-dorsal length (48.5-53.1 vs. 41.9-44.1 %SL), a wider body at dorsal-fin origin (11.3-16.7 vs. 9.4-10.3 %SL), greater dorsal (18.1-21.1 vs. 15.1-17.0 %SL) and lateral (20.9-24.1 vs. 17.4-18.9 %SL) head length, a wider head (14.5-18.5 vs. 11.6-13.3 %SL), and moderately forked (vs. emarginated) caudal fin. Further, molecular analysis confirms the distinctiveness of S. sonarengaensis sp. nov. from its congeners found in northeast India by significant divergences with uncorrected genetic distance ranging from 3.7% to 17.3% in the mitochondrial cytochrome c oxidase subunit I (COI) gene dataset. The phylogenetic position of the new species with its sister species was evaluated using maximum likelihood and Bayesian analysis. The species delimitation approaches assemble species by automatic partitioning (ASAP) and Poisson tree processes (PTP) utilized for testing species assignments consistently identified our test group as a distinct species from its sister species. Although the new species lacks typical morphological adaptations usually associated with a subterranean life, such as complete absence (or vestigial presence) of eyes and pigmentation, it exhibits a reduction of pigmentation when compared to the epigean congeners.

Molecular characterization, spatiotemporal expression, and background adaptation regulation of tyrosinase in loach (Misgurnus anguillicaudatus).

The Poyang Lake region is home to large-blackspot loaches (LBL), small-blackspot loaches (SBL), and non-blackspot loaches (NBL), Misgurnus anguillicaudatus. To investigate the impact of tyrosinase on spot development, the complementary DNAs (cDNA) of tyrosinase in M. anguillicaudatus (designated as Matyr) were cloned using the rapid amplification of cDNA ends (RACE)-PCR method. The full-length cDNA for Matyr was 2020 bp, and the open-reading frame comprised 1617 bp, encoding a predicted protein with 538 amino acids. Phylogenetic studies revealed that MaTyr was first grouped with Tyr of Triplophysa tibetana and Leptobotia taeniops, and then Tyr of other cyprinid fish. The quantitative reverse-transcription-PCR results show that Matyr was highly expressed in the muscle, caudal fin, and dorsal skin. The Matyr gene's messenger RNA expression pattern steadily increased from the fertilized ovum period to the somitogenesis period, and from the muscle effect stage to 6 days after fertilization, it considerably increased (p < 0.01). The Matyr hybridization signals with similar location could be found in all developmental stages of three kinds of loaches using whole-mount in situ hybridization (WISH) technology and were the strongest during the organ development period and melanin formation period. Dot hybridization signals in LBLs rapidly spread to the back of the body beginning at the period when the eyes first formed melanin, and their dimensions were larger than those of NBLs during the same time period. The body color of loaches could change reversibly with black/white background adaptation. The α-msh, mitfa, and tyr are mainly expressed in loaches adapted with a black background. Tyr gene could be involved in the development of blackspots and body color polymorphism, and contribute to organ development in the loach.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.