Limited evidence for extensive genetic differentiation between X and Y chromosomes in Hybognathus amarus (Cypriniformes: Leuciscidae)

Sex determination systems and genetic sex differentiation across fishes are highly diverse but are unknown for most Cypriniformes, including Rio Grande silvery minnow (Hybognathus amarus). In this study, we aimed to detect and validate sex-linked markers to infer sex determination system and to demonstrate the utility of combining several methods for sex-linked marker detection in nonmodel organisms. To identify potential sex-linked markers, Nextera-tagmented reductively amplified DNA (nextRAD) libraries were generated from 66 females, 64 males, and 60 larvae of unknown sex. These data were combined with female and male de novo genomes from Nanopore long-read sequences. We identified five potential unique male nextRAD-tags and one potential unique male contig, suggesting an XY sex determination system. We also identified two single-nucleotide polymorphisms (SNPs) in the same contig with values of FST, allele frequencies, and heterozygosity conforming with expectations of an XY system. Through PCR we validated the marker containing the sex-linked SNPs and a single nextRAD-tag sex-associated marker but it was not male specific. Instead, more copies of this locus in the male genome were suggested by enhanced amplification in males. Results are consistent with an XY system with low differentiation between sex-determining regions. Further research is needed to confirm the level of differentiation between the sex chromosomes. Nonetheless, this study highlighted the power of combining reduced representation and whole-genome sequencing for identifying sex-linked markers, especially when reduced representation sequencing does not include extensive variation between sexes, either because such variation is not present or not captured.

The transcriptome sequencing analysis reveals immune mechanisms of soybean fermented powder on the loach (Misgurnus anguillicaudatus) in response to Lipopolysaccharide (LPS) infection

The loach (Misgurnus anguillicaudatus), a small commercial fish that is widely cultivated for its high-quality protein, vitamins, minerals, and essential amino acid, is a member of the genus Misgurnus and the family Cyprinidae. In this study, we gave the LPS-injected loach fermented soybean meal and used transcriptome sequencing to investigate the impact of the fermented soybean powder on the loach's immune system. 3384 up-regulated genes and 12116 down-regulated genes were found among the 15500 differentially expressed genes, according to the results. The differentially expressed genes were shown to be involved in cellular processes, metabolic processes, cellular anatomical entities, and binding, according to the Go functional annotation. Meanwhile, the KEGG enrichment analysis indicated that the soybean fermented powder treated groups showed significant differences in DNA replication, Nucleotide excision repair, Fanconi anemia pathway, and Base excision repair pathways, suggesting that these pathways are closely related to the enhancement of the immune function of loach by soybean fermented powder. The particular conclusions not exclusively can provide a new conception for the rational utilization of soybean fermented powder but also can provide theoretical guidance for the subsequent healthy breeding of loach.

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.

Metabolomic and transcriptomic profiling reveals the effect of dietary protein and lipid levels on growth performance in loach (Paramisgurnus dabryanus)

The subject of this study was to explore the optimum requirements of loach (Paramisgurnus dabryanus) regarding dietary proteins and lipids and discuss the underlying mechanism. We designed nine diets to determine the effects of different levels of dietary crude protein (CP: 30%, 35%, and 40%) and ether extract (EE: 6%, 10%, and 14%) on the growth performance and metabolism of P. dabryanus. In total, 2160 healthy P. dabryanus (5.19 ± 0.01 g) were divided into nine groups with four replications at 60 fish per barrel stocking density. The trial lasted for eight weeks. Serum and liver samples were gathered for metabolomic and transcriptomic analyses. The results showed that the specific growth rate of P. dabryanus in the CP40EE10 group was the fastest and notably higher than that in other groups (P< 0.05). Analysis of the metabolome results found that the mTOR signaling pathway, glycerophospholipid metabolism, D-arginine and D-ornithine metabolism were significantly enriched pathways in the CP40EE10 group compared with the other groups (P< 0.05). Moreover, the transcriptomic analysis of differentially expressed genes (DEGs) showed that the expression of ARG (arginase) involved in protein synthesis was significantly upregulated in the CP40EE10 group compared to the slowest growing group (P< 0.05). Additionally, the expression of SPLA2 (secretory phospholipase A2) involved in lipid metabolism and FBP (fructose-1,6-bisphosphatase) involved in glucose metabolism were all significantly downregulated in the CP30EE6 group compared with the CP40EE10 group (P< 0.05). Furthermore, the analysis of differentially expressed metabolites (DEMs) and DEGs co-enriched in the KEGG pathway revealed that the significantly enriched pathways were arginine and proline metabolism, glycerophospholipid metabolism, and glycolysis/gluconeogenesis in CP40EE10 compared with other groups (P< 0.05). We conclude that including 40% CP and 10% EE in the P. dabryanus diet could result in a better growth rate. We hypothesized from metabolomic and transcriptomic analyses that the CP40EE10 diet might promote the growth of P. dabryanus by promoting protein synthesis, lipid metabolism, and energy production.

The first cavefish in the Dinaric Karst? Cave colonization made possible by phenotypic plasticity in Telestes karsticus

Cave animals are an excellent model system for studying adaptive evolution. At present, however, little is known about the mechanisms that enable surface colonizers to survive in the challenging environment of caves. One possibility is that these species have the necessary genetic background to respond with plastic changes to the pressures of underground habitats. To gain insight into this process, we conducted a comparative study with the fish species Telestes karsticus, which occurs in a hydrological system consisting of an interconnected stream and a cave. Results showed that T. karsticus resided year-round and spawned in Sušik cave, making it the first known cavefish in the Dinaric Karst. Cave and surface populations differed in morphological and physiological characteristics, as well as in patterns of gene expression without any evidence of genetic divergence. To test whether observed trait differences were plastic or genetic, we placed adult fish from both populations under light/dark or constant dark conditions. Common laboratory conditions erased all morphometric differences between the two morphs, suggesting phenotypic plasticity is driving the divergence of shape and size in wild fish. Lighter pigmentation and increased fat deposition exhibited by cave individuals were also observed in surface fish kept in the dark in the laboratory. Our study also revealed that specialized cave traits were not solely attributed to developmental plasticity, but also arose from adult responses, including acclimatization. Thus, we conclude that T. karsticus can adapt to cave conditions, with phenotypic plasticity playing an important role in the process of cave colonization.

Loss of behavioral stress response in blind cavefish reduces energy expenditure

The stress response is essential for animal self-defense and survival. However, species may exhibit stress response variation depending on their specific environmental and selection pressures. Blind cavefish dwell in cave environments, which differ markedly in stressors and resource availability compared to surface aquatic environments. However, whether blind cavefish exhibit differences in stress response as an adaptation to their cave environments remains unclear. Here, we investigated differences in stress response in six closely related Triplophysa species, including three blind cavefish (T. longibarbata, T. jiarongensis, and T. rosa) and three normal-sighted river fish (T. nasobarbatula, T. dongsaiensis, and T. bleekeri). Results showed that blind cavefish exhibited a range of distinct behavioral responses compared to sighted river fish, including greater levels of activity, shorter duration of freezing, absence of erratic movements or thrashing behavior, and opposite behavioral trends over time. Furthermore, the cavefish species demonstrated attenuated increases in metabolic rate in response to stressors related to novel environments. Cave-dwelling T. rosa also exhibited lower basal hypothalamic-pituitary-inter-renal (HPI) axis-related gene expression levels and stress hormone concentrations compared to river-dwelling T. bleekeri. These results suggest that blind cavefish may have lost their behavioral stress response, potentially mediated by a reduction in basal activity of the HPI axis, thus enabling the conservation of energy by reducing unnecessary expenditure in energy-limited caves.

Toll-interacting protein negatively regulated innate immune response via NF-κB signal pathway in blunt snout bream, Megalobrama amblycephala

Toll-interacting protein (Tollip) is an important negative regulator of Toll-like receptor-mediated innate immunity by preventing excessive proinflammatory responses. The structure and function of Tollip have been well identified in mammals, but the piscine Tollip remains poorly understood. In the present study, a homologue of Tollip was identified and characterized from blunt snout bream (named MaTollip), which was composed of an 831 bp open reading frame encoding a protein of 276 amino acids. Phylogenetic analysis indicated that MaTollip is a novel member of Tollip family and possessed the highest similarity to that of grass carp (99.28%). Multiple alignment of amino acid sequence showed that MaTOLLIP shared a high degree of structural conservation, including a TBD domain, a C2 domain and a CUE domain, with its counterparts from other vertebrates. With regard to tissue-specific expression without immune challenge, MaTollip was constitutively expressed in a wide range of normal tissues, with the highest in the head-kidney and the lowest in the intestine. MaTollip expression in the head-kidney was strongly upregulated upon LPS stimulation and A. hydrophila infection. Fluorescence microscopic analysis revealed that the green fluorescent protein-TOLLIP was localized predominantly in the cytoplasm of EPC cells in a dot-like state. When MaTollip was overexpressed in HEK-293T and EPC cells, it could significantly inhibit the activity of nuclear factor-κB (NF-κB) promoter in a dose dependent manner. MaTollip overexpression in MAF cells lowered drastically the transcriptional expression level of lipopolysaccharide-induced proinflammatory cytokines (IL-1β, IL-6 and IL-8), whereas they were dramatically promoted by MaTollip knock down with siRNA. Taken together, this study demonstrated that MaTollip played a pivotal role in mediating host innate immune response to pathogen invasion, and unveiled the involvement of MaTollip in NF-κB-mediated transcription of inflammation genes, which paved the way for further studies of immune negative regulation mechanisms mediated by Tollip in fish.

Effects of 17α-Methyltestosterone on the Transcriptome and Sex Hormones in the Brain of Gobiocypris rarus

17α-Methyltestosterone (MT), a synthetic environmental endocrine disruptor with androgenic effects, has been shown to disrupt the reproductive system and inhibit germ cell maturation in Gobiocypris rarus. To further investigate the regulation of gonadal development by MT through the hypothalamic-pituitary-gonadal (HPG) axis, G. rarus were exposed to 0, 25, 50, and 100 ng/L of MT for 7, 14, and 21 days. We analyzed its biological indicators, gonadotropin-releasing hormone (GnRH), gonadotropins, reproduction-related gene expression, and brain tissue transcriptome profiles. We found a significant decrease in the gonadosomatic index (GSI) in G. rarus males exposed to MT for 21 days compared to the control group. GnRH, follicle-stimulating hormone (FSH), and luteinizing hormone (LH) levels, as well as the expressions of the gnrh3, gnrhr1, gnrhr3, fshβ, and cyp19a1b genes, were significantly reduced in the brains of both male and female fish when exposed to 100 ng/L MT for 14 days compared to the controls. Therefore, we further constructed four RNA-seq libraries from 100 ng/L MT-treated groups of male and female fish, obtaining 2412 and 2509 DEGs in male and female brain tissue, respectively. Three common pathways were observed to be affected in both sexes after exposure to MT, namely, nicotinate and nicotinamide metabolism, focal adhesion, and cell adhesion molecules. Furthermore, we found that MT affected the PI3K/Akt/FoxO3a signaling pathway through the upregulation of foxo3 and ccnd2, and the downregulation of pik3c3 and ccnd1. Therefore, we hypothesize that MT interferes with the levels of gonadotropin-releasing hormone (GnRH, FSH, and LH) in G. rarus brains through the PI3K/Akt/FoxO3a signaling pathway, and affects the expression of key genes in the hormone production pathway (gnrh3, gnrhr1 and cyp19a1b) to interfere with the stability of the HPG axis, thus leading to abnormal gonadal development. This study provides a multidimensional perspective on the damaging effects of MT on fish and confirms that G. rarus is a suitable model animal for aquatic toxicology.

Complete Mitogenome of the Triplophysa bombifrons: Comparative Analysis and Phylogenetic Relationships among the Members of Triplophysa

In the last decade, the phylogenetic relationships within the genus Triplophysa have become controversial, due to a lack of molecular data. The mitochondrial genome plays a vital role in the reconstruction of phylogenetic relationships and in revealing the molecular evolution of bony fishes. Herein, we obtained the complete mitogenome of Triplophysa bombifrons via HiFi reads of the Pacbio Sequel II system and DNBSEQ short-reads. We compared all available mitogenomes of the Triplophysa genus and reconstructed the phylogeny of Nemacheilidae, based on the mitogenomes, using maximum likelihood (ML) methods. The results show that the complete mitogenome sequence of T. bombifrons was circular and 16,568 bp in length, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA), 2 ribosomal RNA (rRNA), and a typical control region (D-loop). The most common start codons were ATG, except for cox1, and TAA/TAG were the stop codons for all PCGs. In total, 677 SNPs and 9 INDELs have been found by comparing the sequence divergence between this study and previous reports. Purity selection was found in all PCGs. Phylogeny was inferred by analyzing the 13 PCGs and the concatenated nucleotide sequences of 30 mitogenomes. The phylogenetic analyses based on the nucleotides of the 13 PCGs supported the assumption that the Triplophysa genus can be divided into 4 main clades and demonstrated that T. bombifrons and T. tenuis are closely related species for the first time. This study laid the foundation for further study on the mitogenome and phylogeny of Nemacheilidae fishes.

Transcriptomic analysis of intermuscular bone development in barbel steed (Hemibarbus labeo)

Intermuscular bones (IBs), which are little, bony spicules in muscle, are embedded in lower teleosts' myosepta. Despite the importance of studying IB development in freshwater aquaculture species, the genes associated with IB development need to be further explored. In the present study, we identified four stages of IB development in barbel steed (Hemibarbus labeo), namely stage 1: IBs have not emerged, stage 2: a few small IBs have emerged in the tail, stage 3: longer IBs gradually emerged in the tail and stage 4: all of the IBs in the tail are mature and long, via Alizarin red staining. Subsequently, we used the HiseqXTen platform to sequence and de novo assemble the transcriptome of epaxial muscle (between 35th and 40th myomere) of barbel steed at 29 days (stage 1) and 42 days (stage 3) after hatching. A total of 190,814 unigenes were obtained with an average length and N50 of 648 bp and 1027 bp, respectively. We found 2174 differentially expressed genes (DEGs) between stages 1 and 3, of which 378 and 1796 were up- and down-regulated, respectively. Functional enrichment analysis showed that several DEGs functioned in ossification, positive regulation of osteoblast differentiation, osteoblast differentiation, and BMP signaling pathway, and were further enriched in signal pathway, including osteoclast differentiation, TGF-β signaling pathway, cytokine-cytokine receptor interaction, Jak-STAT signaling pathway, and other KEEG pathways. In conclusion, we identified genes that may be related to IB development, such as kazal type serine peptidase inhibitor domain 1 (KAZALD1), extracellular matrix protein 1 (ECM1), tetranectin, bone morphogenetic protein 1 (bmp1), acid phosphatase 5 (ACP5), collagen type XI alpha 1 chain (COL11A1), matrix metallopeptidase 9 (MMP9), pannexin-3 (PANX3), sp7 transcription factor (Sp7), and c-x-c motif chemokine ligand 8 (CXCL8), by comparing the transcriptomes of epaxial muscle before and after IB ossification. This study provided a theoretical basis for identifying the molecular mechanisms underlying IB development in fish.

Three new species of Homatula (Teleostei: Nemacheilidae) from Yunnan, China, with comments on habitat conservation

Based on the morphological comparisons and molecular results, three new species of Homatula, i.e., H. geminusclathratus sp. nov., H. microcephala sp. nov., and H. longibarbatus sp. nov., have been described and named from the Lancang-jiang (the upper Mekong River) and the Chuan-he (the upper Black River, a tributary of the Red River) basins. The loaches of Homatula from the Lancang-jiang and the Chuan-he can be distinguished via morphology, genetics, and geographic distribution. All of the 10 recorded species distributed in the Nu-jiang (the upper Salween River), the Lancang-jiang, and the upper Black River share the following combination of character states: whole body, except head, densely scaled; lateral line complete; and a short adipose crest along the dorsal midline of the caudal peduncle, anteriorly not reaching vertically through the anal-fin origin. Species with these characters are called the densely-scaled group of Homatula. The three newly described species belong to the densely-scaled group of Homatula. Based on molecular phylogenetics, these Homatula species form a monophyletic group that can be divided into two clades, the densely-scaled group and the non-densely-scaled group. The densely-scaled group of Homatula includes 13 species occurring between the Nu-jiang and the upper Black River. The non-densely-scaled group is non-monophyletic and includes 14 species that are distributed in the Red, Pearl, Yangtze, and Yellow River basins. Species of the non-densely-scaled group are clustered into four sub-clades that are constrained to the four river basins. Homatula exclusively inhabits mountain streams with rapid or gentle currents, vauclusian springs, underground rivers connected to streams, and ditches near villages and farmland. No specimens of Homatula were collected from the main streams of Lixian-jiang, Lancang-jiang, and Nu-jiang as well as their large tributaries. Small environmental changes in the habitat of Homatula, such as water pollution or extensive human use, can lead to species/population extinction. Effective conservation of rare and endemic fishes, like loaches of Homatula, entails systematic observations and more targeted protection.

A new set of microsatellite markers for Phoxinus lumaireul senso lato, Phoxinus marsilii and Phoxinus krkae for population and molecular taxonomic studies

Minnows of the genus Phoxinus are common and an often highly abundant fish species in Palearctic freshwater habitats. Phoxinus species have a complex evolutionary history, phylogenetic relationships are not well understood and there are a number of unresolved taxonomic problems. There are currently 23 different mitochondrial genetic lineages identified in the genus Phoxinus, 13 of which are recognized as valid species. The taxonomic status of these lineages requires resolution, including the degree to which they can interbreed. Suitable nuclear molecular markers for studies of population divergence and interbreeding between morphotypes and mitochondrial lineages are lacking for Phoxinus species. Therefore, the authors developed a set of microsatellite markers using genomic information from Phoxinus lumaireul and tested their suitability for this and two related species, Phoxinus krkae and Phoxinus marsilii. Out of 16 microsatellite candidate loci isolated, 12 were found to be in Hardy-Weinberg equilibrium when tested on two P. lumaireul senso lato populations. Seven loci amplified across the three species, enabling the study of intraspecific genetic diversity and population structure within P. marsilii and P. krkae. The markers were able to clearly resolve differences among the three tested species, including the recently described P. krkae, and are therefore suitable for the detection of introgression and hybridization among populations consisting of mixtures of two or more of P. lumaireul s. l., P. marsilii and P. krkae.

Towards the phylogenetic placement of the enigmatic African genus Prolabeops Schultz, 1941

The small cyprinid genus Prolabeops Schultz, 1941 is restricted to the Nyong and Sanaga River systems in Cameroon. In the past, the genus had been suggested to be either a member of the Labeoninae, Torinae or the Smiliogastrinae mainly on the basis of morphological similarities, and it is nowadays considered as incertae sedis within the Cypriniformes. This study provides the first attempt to reveal the phylogenetic position of Prolabeops using molecular data. For this purpose, the authors sequenced a large fraction of the mitochondrial genome (c. 13,600 bp), including all mitochondrial protein coding genes, of two Prolabeops melanhypopterus specimens and an additional four Enteromius specimens. The large-scale phylogenetic analysis was based on an alignment including all mitochondrial protein coding genes of 902 specimens representing c. 899 cypriniform species. Prolabeops was clearly recovered within the African Smiliogastrinae, forming a weakly supported clade together with Enteromius jae, Enteromius hulstaerti and Barboides gracilis. The study data underline the urgent need of a thorough taxonomic revision of the small African barbs collectively placed in the genus Enteromius.

Cryptic genetic divergence of the red dwarf rasbora, Microrasbora rubescens, in and around Inle Lake: implications for the origin of endemicity in the ancient lake in Myanmar

Inle Lake, an ancient lake located in the Shan Plateau of Myanmar, is a biogeographically attractive region with high fish endemism. Some endemic species inhabit the lake as well as the surrounding areas. The genetic and ecological relationships between populations in the lake and surrounding areas provide important insights into the process underlying ichthyofaunal formation in Inle Lake. In this study, the authors focused on red dwarf rasbora Microrasbora rubescens, an endemic genus and species in this region, and estimated its population structure and evolutionary scenario based on genome-wide polymorphism, mtDNA and geometric morphometric analyses using samples from Inle Lake and three areas surrounding the lake. The results showed that M. rubescens comprises at least three genetically divergent lineages (Inle, Heho and Hopong) with distinct geographic structures consistent with nuclear and mtDNA data. In contrast, there was no clear regional differentiation in morphology. The divergence time estimation based on mtDNA suggests that the Hopong lineage diverged at 2.7 Ma and the Inle and Heho lineages diverged at 1.9 Ma - consistent with the nuclear DNA results. The deep divergence observed in the endemic species supports the ancient history of ichthyofaunal development in this region. The distinct regional differentiation and morphological conservatism of this species might have been shaped by niche conservatism in stagnant water environments that limit dispersal and morphological diversification. Future comprehensive genetic and morphological analyses and comparisons for other native species should reveal the geographic and ecological processes that shaped the ichthyofauna in this region.

Multilocus Phylogeography and Population Genetic Analyses of Opsariichthys hainanensis Reveal Pleistocene Isolation Followed by High Gene Flow around the Gulf of Tonkin

The ichthyofauna of continental islands is characterized by immigration through a land bridge due to fluctuating sea levels. Hainan Island is adjacent to the southern margin of mainland China and provides opportunities for understanding the origin and diversification of freshwater fishes. The aim of our study was to evaluate the level of genetic variation and phylogeographic structure of Opsariichthys hainanensis on Hainan Island and mainland China, using mtDNA cyt b gene (1140 bp) and D-loop (926 bp), nuclear RAG1 gene (1506 bp), and 12 microsatellite loci. Mitochondrial phylogenetic analysis identified five major lineages according to the geographical distribution from different populations. We suggested that two dispersal events occurred: the population in the Changhua River migrated to the Red River (Lineage B), and the populations in the South Hainan region moved northwards to the North Hainan region. However, populations in Northwest Hainan Island dispersed to the populations around the Gulf of Tonkin (Lineage A1) and populations in Northeast Hainan Island dispersed to the populations in mainland China (Lineage A2). Our results indicated that the populations of O. hainanensis suffered a bottleneck event followed by a recent population expansion supported by the ABC analysis. We suggest that O. hainanensis populations were found mostly in the lowlands and a lack of suitable freshwater habitat in southern mainland China and Hainan during the Last Interglacial period, and then expansion occurred during the Last Glacial Maximum.

Regulation of Rhesus glycoprotein-related genes in large-scale loach Paramisgurnus dabryanus during ammonia loading

Waterborne ammonia is one of the crucial issues that limited production and animal health in aquaculture. Ammonia-tolerant varieties are highly desired in intensive fish farming. Screening for the key regulatory genes of ammonia tolerance is essential for variety breeding. According to the previous hypothesis, Rh glycoproteins play an important role in ammonia excretion in teleosts. However, the ammonia defensive mechanisms are not well described at present for large-scale loach (Paramisgurnus dabryanus), a typical air-breathing and commercially important fish in East Asia. Here we show that the transcription of Rh glycoprotein-related genes was significantly affected by ammonia exposure in this species. Probit analysis showed that 96 h-LC₅₀ of NH₄Cl at 23 ℃ and pH 7.2 was 92.64 mmol/L. A significant increase of Rhcg expression in gills was observed after 48 h of 60 mmol/L and 36 h of 80 mmol/L NH₄Cl exposure, suggesting that Rhcg present on the apical side of the branchial epithelium facilitates NH₃ excretion out of gills. A high concentration of acute ammonia exposure induced elevated Rhbg transcript in the gills of large-scale loaches, while a slight change in Rhbg expression was observed in response to lower ammonia, suggesting that transcriptions of Rhbg genes are activated by a considerably high level of ambient ammonia to eliminate excessive endogenous nitrogen. The Rhag mRNA level in gills of large-scale loaches increased markedly with the prolonging of exposure time from 0 to 36 h of ammonia loading, suggesting Rhag localized in gills may be primarily associated with ammonia handling. During 7-21 days of ammonia exposure, the expression of most Rh glycoproteins-related genes in the gills decreased, indicating that the functional role of Rh glycoproteins is not primarily associated with ammonia defense over a long period (more than 7 days). Although a significant transcript of Rhbg was found in the skin of a large-scale loach, the lack of Rhcg and down-regulation of Rhag may indicate that the skin is not an essential location of ammonia excretion, at least when submerged to high levels of ammonia in the environment. In conclusion, Rh glycoproteins localized in gills as ammonia transporters play a momentous role in ammonia detoxification in this species during acute ammonia loading. However, it does not show a positive function during long-term ammonia exposure. Furthermore, the physiological function of Rh glycoproteins localized in the skin is still unclear and deserves further study.

Characterization of aquaporin-1ab (Aqp1ab) mRNA in mud loach (Misgurnus mizolepis) exposed to heavy metal and immunostimulant stimuli

Aquaporins (AQPs) are key proteins that regulate fluid homeostasis in cells via modulating osmotic water transport. In the present study, we identified three variants of Aqp1ab transcript (mmAQP1ab x1, mmAQP1ab x2, and mmAQP1ab x3) in mud loaches (Misgurnus mizolepis), and their expression patterns were examined in response to heavy metal and immunostimulant exposure. Mud loach Aqp1ab gene has a somewhat different organizational structure (i.e. five exons interrupted by four introns) compared to most other teleostean Aqp1ab orthologues, which have four exons. The 5'-flanking regulatory region of Aqp gene showed diverse transcription factor binding motifs, particularly those associated with stress/immune responses. Developmental expression patterns indicated that Aqp1ab mRNA was maternally inherited, presumably important for fine-tuning gene expression during embryonic and early larval developments. Expression of mud loach Aqp1ab mRNA was significantly and differentially modulated in several tissues (intestine, kidneys, spleen, and liver) in response to various heavy metal treatments. In addition, Aqp1ab gene expression was highly induced in response to immune challenge (LPS and polyI:C injections). Collectively, our results suggested that AQPs are multifunctional effectors playing diverse roles in cellular pathways relevant to immune and/or stress adaptation responses, in addition to their involvement in osmoregulation.

Exposure to Contaminated River Water is Associated with Early Hatching and Dysregulation of Gene Expression in Early Life Stages of the Endangered Copper Redhorse (Moxostoma hubbsi)

The copper redhorse (Moxostoma hubbsi) is an endangered fish that spawns exclusively in the Richelieu River (Quebec, Canada). Tributaries of the Richelieu are contaminated with high levels of current-use pesticides, which may impact early-life stage (ELS) copper redhorse and other native fishes. We assessed the effects of exposure to contaminated river water on ELS copper redhorse and river redhorse (Moxostoma carinatum), a related fish that shares the copper redhorse's spawning grounds and nursery habitat. A riverside flow-through system was used to expose copper and river redhorse embryos (1000 each) to Richelieu River water or laboratory water as a control. Fish were maintained until 14 days posthatch, and water samples were taken daily for chemical analysis. Following a heavy rain event, concentrations of two neonicotinoid pesticides, clothianidin and thiamethoxam, exceeded water quality guidelines for aquatic life (20 ng/L). Using nontargeted screening, we tentatively identified an additional 24 pharmaceutical and personal care products and 23 pesticides in river water. Effects of river water on ELS fish were observed in both species, but the copper redhorse appeared to be more sensitive. Fish exposed to river water hatched 10.7 (copper redhorse) and 2.4 (river redhorse) cumulative degree days earlier than controls. Copper redhorse survival was significantly lower in river water (73 ± 16%) compared to laboratory water (93 ± 3%), whereas river redhorse survival was similar between treatments (84 ± 6% and 89 ± 4%, respectively). Sequencing of copper redhorse larvae RNA revealed 18 differentially expressed genes (DEGs) following 14 days of exposure to river water. Eight up-regulated DEGs were linked to immune function and injury response, and seven down-regulated DEGs were involved with digestion and nutrient absorption. The present study provided valuable data on the effects of ELS exposure to a real-world mixture of contaminants in two fish species of concern. Environ Toxicol Chem 2022;41:1950-1966. © 2022 SETAC.

Re-description of the loach species Homatula variegata (Dabry de Thiersant, 1874) (Pisces: Nemacheilidae) from the middle Huang-He basin in Shaanxi Province of Central China

Despite the wide recognition of Homatula variegata (Dabry de Thiersant, 1874) as a valid Chinese loach, its type locality, identity and distribution still remain contentious. A molecular phylogenetic analysis based on the mitochondrial cytochrome b (cyt b) gene for samples from all known range of the species showed that three distinct species are involved. Morphological comparisons, coupled with examination on the types, confirmed that H. variegata s. str. is characterized by having a sparsely scaled predorsal body, an adipose crest along the dorsal midline of the caudal peduncle that extends forward to the vertical through the posterior end of the anal-fin base and a broadly rounded caudal fin. It is found merely in the Wei-He of the Huang-He basin. The Jinsha-Jiang population, previously misidentified as H. variegata, represents a distinct species, for which Homatula oxygnathra is the available name. Homatula laxiclathra, endemic to the Wei-He, is also a valid species distinct from H. variegata.

Identification, expression patterns, evolutionary characteristics and recombinant protein activities analysis of CD209 gene from Megalobrama amblycephala

CD209 is a type II transmembrane protein in the C-type lectin family, which is involved in the regulation of innate and adaptive immune system. Although it has been widely studied in mammals, but little has been reported about fish CD209 genes. In the present study, Megalobrama amblycephala CD209 (MaCD209) gene was cloned and characterized, its expression patterns, evolutionary characteristics, agglutinative and bacteriostatic activities were analyzed. These results showed that the open reading frame (ORF) of MaCD209 gene was 795 bp, encoding 264 aa, and the calculated molecular weight of the encoded protein was 29.7 kDa. MaCD209 was predicted to contain 2 N-glycosylation sites, 1 functional domain (C-LECT-DC-SIGN-like) and 1 transmembrane domain. Multiple sequence alignment showed that the amino acid sequence of MaCD209 was highly homologous with that of partial fishes, especially the highly conserved C-LECT-DC-SIGN-like domain and functional sites of CD209. Phylogenetic analysis showed that the CD209 genes from M. amblycephala and other cypriniformes fishes were clustered into one group, which was reliable and could be used for evolutionary analysis. Then, adaptive evolutionary analysis of teleost CD209 was conducted, and several positively selected sites were identified using site and branch-site models. Quantitative real-time PCR analysis showed that MaCD209 gene was highly expressed in the liver and heart. Moreover, the expression of MaCD209 was significantly induced upon Aeromonas hydrophila infection, with the peak levels at 4 h or 12 h post infection. The immunohistochemical analysis also revealed increased distribution of MaCD209 protein post bacterial infection. In addition, recombinant MaCD209 (rMaCD209) protein was prepared using a pET32a expression system, which showed excellent bacterial binding and agglutinative activities in a Ca²⁺-independent manner. However, rMaCD209 could only inhibit the proliferation of Escherichia coli rather than A. hydrophila. In conclusion, this study identified the MaCD209 gene, detected its expression and evolutionary characteristics, and evaluated the biological activities of rMaCD209 protein, which would provide a theoretical basis for understanding the evolution and functions of fish CD209 genes.

Enrichment of polystyrene microplastics induces histological damage, oxidative stress, Keap1-Nrf2 signaling pathway-related gene expression in loach juveniles (Paramisgurnus dabryanus)

Polystyrene microplastics (PS-MPs, particle size<5 mm) cause great harm to aquatic organisms. However, their precise effects are not completely understood. In China, placing plastic film at the pond bottom has become an important loach aquaculture mode. In this mode, MPs will affect loach health. This study investigated the enrichment of PS-MPs and its effects on the growth, liver histomorphology, antioxidant enzymes, and Keap1-Nrf2 signaling pathway-related gene expression in loach juveniles (Paramisgurnus dabryanus). The loach juveniles were raised at the concentration of 1000 μg/L fluorescent polystyrene microplastics (PS-MPs) with particle size of 0.5 µm or 5 µm for seven days, the results showed that fluorescent PS-MPs were found to be enriched in liver, intestine, and gill, and the enrichment amount was higher in liver than in gill and intestine (P < 0.05). Furthermore, the enrichment amount of different-sized PS-MPs was different in liver, gill, and intestine. The loach juveniles were cultured for 21 days in the water of the concentration of 100 or 1000 μg/L PS-MPs with particle size of 0.5 µm or 5 µm, the results showed that the survival rate, weight gain rate, and specific growth rate of loach juveniles were significantly reduced. The histological analysis revealed that PS-MPs caused liver damage. The activities of superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GSH-PX), and acetylcholinesterase (AChE) were decreased with the extended exposure to PS-MPs. Generally, the expressions of Nrf2 and Keap1 showed the similar change trend. From 7-14 day, the expression trend of oxidative stressed-related genes was not completely consistent with that of Nrf2 gene, but on day 21, the gene expression trend of oxidative stress-related SOD, CAT, and GSH-PX in the downstream of Keap1-Nrf2 signaling pathway was roughly consistent with that of Nrf2 gene. Basically, the change trends of these three gene expression were similar to those of their corresponding enzyme activities. This study provides theoretical basis for the toxicological effects of PS-MPs on freshwater fish.

De novo gonad transcriptome analysis of elongate loach (Leptobotia elongata) provides novel insights into sex-related genes

Elongate loach (Leptobotia elongata) is endemic in middle and upper reaches of the Yangtze River in China. Because of many anthropogenic factors such as overfishing and dam construction, the loach has become a highly endangered species. So far, the genomic resources which benefit for species conservation and utilization are still lacking in elongate loach. Therefore, the first gonad transcriptome analysis of the loach was conducted in this study, providing novel insights into sex-related genes. A total of 286,800,660 clean reads with a total length of 42.02 Gb were obtained. 18,975 differentially expressed genes (DEGs) were identified, where 12,976 DEGs, especially Sox9a, Sox9b, Igf2 and Fgfr2, were upregulated in the testis, and 5999 DEGs, especially Zp3, Eg2, Plk1, Ccnb1, Cdc20 and Mos, were upregulated in the ovary. Meanwhile, some testis-specific genes (i.e. Cald1, Atp1a, Muc2 and Ca2) and ovary-specific genes (i.e. Ca4, Tuba, Acp5, Ccna, Larp6 and Nop4) were identified and verified. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs, we found a series of enrichment pathways related to reproduction in elongate loach, such as the MAPK signaling pathway, oxytocin signaling pathway and oocyte meiosis pathway. Twelve DEGs were randomly selected to verify RNA-seq results by qPCR. In conclusion, this study provides a data source to study the molecular characteristics and regulatory mechanisms of sex-related genes in elongate loach, which has a potential to improve the resource protection and aquaculture production of the loach.

Molecular cloning, characterization and expression analysis of P53 from high latitude fish Phoxinus lagowskii and its response to hypoxia

As an intermediate link between multiple cellular stresses and cellular responses, p53, together with its upstream and downstream regulators and related genes, constitutes a complex network that regulates cellular stresses and cellular responses. However, no studies have investigated p53 in Phoxinus lagowskii, particularly the expression of p53 under different hypoxic conditions. In the present study, the cDNA of p53 from the Phoxinus lagowskii was cloned by the combination of homology cloning and rapid amplification of cDNA ends (RACE) approaches. The full-length cDNA of Pl-p53 was 1878 bp, including an open reading frame (ORF) of 1116 bp encoding a polypeptide of 371 amino acids with a predicted molecular weight of 41.22 kDa and a theoretical isoelectric point of 7.38. Quantitative real-time (qRT) PCR assays revealed that Pl-p53 was commonly expressed in all tissues examined, with highest expression in the heart. In addition, we investigated the expression of Pl-p53 in different tissues under different hypoxic conditions. In the short-term hypoxia group, Pl-p53 expression was down-regulated in both the brain and heart. The Pl-p53 expression was significantly elevated at 6 h in the muscle and liver, and was significantly up-regulated at 24 h in spleen. These results suggest that Pl-p53 plays different regulatory roles and provide a theoretical basis for the changes of p53 in fish facing hypoxic environments.