Gill Transcriptome Sequencing and De Novo Annotation of Acanthogobius ommaturus in Response to Salinity Stress

Unraveling the molecular mechanisms of fish physiological response to freshwater salinization: A comparative multi-tissue transcriptomic study in a river polluted by potash mining

Freshwater salinization is an escalating global environmental issue that threatens freshwater biodiversity, especially fish populations. This study aims to uncover the molecular basis of salinity physiological responses in a non-native minnow species (Phoxinus septimaniae x P. dragarum) exposed to saline effluents from potash mines in the Llobregat River, Barcelona, Spain. Employing high-throughput mRNA sequencing and differential gene expression analyses, brain, gills, and liver tissues collected from fish at two stations (upstream and downstream of saline effluent discharge) were examined. Salinization markedly influenced global gene expression profiles, with the brain exhibiting the most differentially expressed genes, emphasizing its unique sensitivity to salinity fluctuations. Pathway analyses revealed the expected enrichment of ion transport and osmoregulation pathways across all tissues. Furthermore, tissue-specific pathways associated with stress, reproduction, growth, immune responses, methylation, and neurological development were identified in the context of salinization. Rigorous validation of RNA-seq data through quantitative PCR (qPCR) underscored the robustness and consistency of our findings across platforms. This investigation unveils intricate molecular mechanisms steering salinity physiological response in non-native minnows confronting diverse environmental stressors. This comprehensive analysis sheds light on the underlying genetic and physiological mechanisms governing fish physiological response in salinity-stressed environments, offering essential knowledge for the conservation and management of freshwater ecosystems facing salinization.

Rapid adaptive and acute stimulatory responses to low salinity stress in Pacific white shrimp (Litopenaeus vannamei): Insights from integrative transcriptomic and proteomic analysis

The Pacific white shrimp (Litopenaeus vannamei) is a euryhaline crustacean capable of tolerating a wide range of ambient salinity, but the strategies of hepatopancreas to rapid adaptive or acute stimulatory responses to extremely low salinity fluctuations remains unclear. In this study, we integrated transcriptomic and proteomic analyses on the hepatopancreas derived from rapid adaptative (RA) and acute stimulatory (AS) responses to extremely low salinity stress (0.3 ppt) to unveil specific regulatory mechanisms. The RA group displayed normal epithelial cells and tubule structures, while the AS group showed histological changes and lesions. A total of 754 and 649 differentially expressed genes (DEGs) were identified in RA and AS treatments, respectively. For proteome, a total of 206 and 66 differentially expressed proteins (DEPs) were obtained in the RA/CT and AS/CT comparison groups, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were conducted among the DEGs and DEPs, revealing that metabolic related pathways were significantly enriched pathways in both comparison groups. In addition, correlation analysis of transcriptomic and proteomic results showed that 20 and 3 pairs of DEGs/DEPs were identified in RA vs. CT and AS vs. CT comparison groups, respectively. This study is the first report on the rapid adaptive and acute stimulatory transcriptomic and proteomic responses of L. vannamei to extremely low salinity, shedding light on the mechanisms underlying osmoregulation in euryhaline crustaceans.

Effects of low salinity stress on osmoregulation and gill transcriptome in different populations of mud crab Scylla paramamosain

Animals living in estuaries suffer from rapid and continuous salinity fluctuations, while the global warming and extreme precipitation aggravate this situation. Osmoregulation is important for estuarine animals adapt to salinity fluctuations. The present study investigated the effects of low salinity stress on osmoregulation and gill transcriptome in two populations of mud crab from Hangzhou Bay and Zhangzhou Bay of China, respectively. Crabs were transferred from salinity 25 ppt to 5 ppt for 96 h. Edematous swelling in gill filaments was caused by low salinity stress and was more serious in Zhangzhou Bay population. Gill Na⁺/K⁺-ATPase activity increased (p < 0.01) in both populations under the low salinity stress and was significantly higher (p < 0.01) in Hangzhou Bay population than in Zhangzhou Bay population. According to transcriptome analysis, there were 191 genes differentially expressed under the low salinity stress in gill tissue of both populations. Several ion transport and energy metabolism related pathways, as well as the arginine and proline metabolism pathway, were enriched by these genes. On the other hand, 272 genes were identified to differentially express between two populations under the low salinity stress, but not under the control salinity. The enrichment analysis showed that these genes were mainly related to ion transport, energy metabolism, osmolytes metabolism and methyltransferase activity. In conclusion, the present study suggested that mud crab exploited a combination of extracellular anisosmotic regulation and intracellular isosmotic regulation for osmoregulation under the low salinity stress. Hangzhou Bay population showed a greater osmoregulatory capacity, which is probably due to the enhanced ion transport, energy supply, and osmolytes regulation. Meanwhile, epigenetic modification might also contribute to an inherent osmoregulation ability for Hangzhou Bay population to response to salinity fluctuation rapidly.

Transcriptome analysis reveals the high temperature induced damage is a significant factor affecting the osmotic function of gill tissue in Siberian sturgeon (Acipenser baerii)

BACKGROUND

Maintaining osmotic equilibrium plays an important role in the survival of cold-water fishes. Heat stress has been proven to reduce the activity of Na⁺/K⁺-ATPase in the gill tissue, leading to destruction of the osmotic equilibrium. However, the mechanism of megatemperature affecting gill osmoregulation has not been fully elucidated.

RESULTS

In this study, Siberian sturgeon (Acipenser baerii) was used to analyze histopathological change, plasma ion level, and transcriptome of gill tissue subjected to 20℃, 24℃and 28℃. The results showed that ROS level and damage were increased in gill tissue with the increasing of heat stress temperature. Plasma Cl^- level at 28℃ was distinctly lower than that at 20℃ and 24℃, while no significant difference was found in Na⁺ and K⁺ ion levels among different groups. Transcriptome analysis displayed that osmoregulation-, DNA-repair- and apoptosis-related terms or pathways were enriched in GO and KEGG analysis. Moreover, 194 osmoregulation-related genes were identified. Amongst, the expression of genes limiting ion outflow, occluding (OCLN), and ion absorption, solute carrier family 4, member 2 (AE2) solute carrier family 9, member 3 (NHE3) chloride channel 2 (CLC-2) were increased, while Na⁺/K⁺-ATPase alpha (NKA-a) expression was decreased after heat stress.

CONCLUSIONS

This study reveals for the first time that the effect of heat stress on damage and osmotic regulation in gill tissue of cold-water fishes. Heat stress increases the permeability of fish's gill tissue, and induces the gill tissue to keep ion balance through active ion absorption and passive ion outflow. Our study will contribute to research of global-warming-caused effects on cold-water fishes.

Transcriptome and molecular regulatory mechanisms analysis of gills in the black tiger shrimp Penaeus monodon under chronic low-salinity stress

Background: Salinity is one of the main influencing factors in the culture environment and is extremely important for the survival, growth, development and reproduction of aquatic animals. Methods: In this study, a comparative transcriptome analysis (maintained for 45 days in three different salinities, 30 psu (HC group), 18 psu (MC group) and 3 psu (LC group)) was performed by high-throughput sequencing of economically cultured Penaeus monodon. P. monodon gill tissues from each treatment were collected for RNA-seq analysis to identify potential genes and pathways in response to low salinity stress. Results: A total of 64,475 unigenes were annotated in this study. There were 1,140 upregulated genes and 1,531 downregulated genes observed in the LC vs. HC group and 1,000 upregulated genes and 1,062 downregulated genes observed in the MC vs. HC group. In the LC vs. HC group, 583 DEGs significantly mapped to 37 signaling pathways, such as the NOD-like receptor signaling pathway, Toll-like receptor signaling pathway, and PI3K-Akt signaling pathway; in the MC vs. HC group, 444 DEGs significantly mapped to 28 signaling pathways, such as the MAPK signaling pathway, Hippo signaling pathway and calcium signaling pathway. These pathways were significantly associated mainly with signal transduction, immunity and metabolism. Conclusions: These results suggest that low salinity stress may affect regulatory mechanisms such as metabolism, immunity, and signal transduction in addition to osmolarity in P. monodon. The greater the difference in salinity, the more significant the difference in genes. This study provides some guidance for understanding the low-salt domestication culture of P. monodon.

Regulating Strategies of Transcription and Alternative Splicing for Cold Tolerance Harpadon nehereus Fish

In recent years, Harpadon nehereus gradually become a dominant species with great potential for exploitation in the East China Sea, and it is worth investigating whether H. nehereus would tolerate cold stress to continue to expand into the colder northern waters. The molecular regulation level is favorable evidence to explore the cold tolerance of H. nehereus, a total of 6,650, 1,936, and 2,772 differentially expressed genes (DEGs) in transcription regulation, and 4,409, 1,250, and 2,303 differential alternative splicing genes (DASGs) in alternative splicing regulation were identified in H. nehereus at 13, 15, and 17°C, respectively, importantly, 47 genes were identified as the key candidate genes for cold tolerance in H. nehereus. In transcription regulation, up-regulated DEGs were enriched in metabolic process terms and ribosome, spliceosome pathway, etc., while down-regulated DEGs were enriched in signal transduction terms, focal adhesion, proteoglycans in cancer pathway, etc., at 13, 15, and 17°C, respectively. In alternative splicing regulation, spliceosome, mRNA surveillance pathway, etc., were significantly enriched in DASGs. In a word, H. nehereus adapts to cold environments mainly through transcription and translation, transmembrane transport, protein modification, etc., while cold stress may also induce some diseases in H. nehereus.

Transcriptome profiles revealed high- and low-salinity water altered gill homeostasis in half-smooth tongue sole (Cynoglossus semilaevis)

Salinity is an important environmental factor that affects fish growth, development, and reproduction. As euryhaline fish, half-smooth tongue sole (Cynoglossus semilaevis) are a suitable species for deciphering the salinity adaptation mechanism of fish; however, the molecular mechanisms underlying low- and high-salinity responses remain unclear. In this study, RNA-seq was applied to characterize the genes and regulatory pathways involved in C. semilaevis gill responses to high- (32 ppt), low- (8 ppt), and control-salinity (24 ppt) water. Gills were rich in mitochondria-rich cells (MRCs) in high salinity. Compared with control, 2137 and 218 differentially expressed genes (DEGs) were identified in low and high salinity, respectively. The enriched functions of most DEGs were metabolism, ion transport, regulation of cell cycle, and immune response. The DEGs involved in oxidative phosphorylation, citrate cycle, and fatty acid metabolism were down-regulated in low salinity. For ion transport, high and low salinity significantly altered the expressions of prlr, ca12, and cftr. In cell cycle arrest and cellular repair, gadd45b, igfbp5, and igfbp2 were significantly upregulated in high and low salinity. For immune response, il10, il34, il12b, and crp increased in high and low salinity. Our findings suggested that alterations in material and energy metabolism, ions transport, cell cycle arrest, cellular repair, and immune response, are required to maintain C. semilaevis gill homeostasis under high and low salinity. This study provides insight into the divergence of C. semilaevis osmoregulation mechanisms acclimating to high and low salinity, which will serve as reference for the healthy culture of C. semilaevis.

RNA Sequencing Analysis Reveals Divergent Adaptive Response to Hypo- and Hyper-Salinity in Greater Amberjack (Seriola dumerili) Juveniles

Simple Summary

The gill tanscriptomes of greater amberjack (Seriola dumerili) reared under different salinity stress were analyzed. The regulatory networks of salinity-related pathways were explored through Kyoto Encyclopedia of Gene and Genome (KEGG) pathway enrichment and bioinformatics analyses. This will be of great value in understanding the molecular basis of salinity adaptation in greater amberjack.

Abstract

Salinity significantly affects physiological and metabolic activities, breeding, development, survival, and growth of marine fish. The greater amberjack (Seriola dumerili) is a fast-growing species that has immensely contributed to global aquaculture diversification. However, the tolerance, adaptation, and molecular responses of greater amberjack to salinity are unclear. This study reared greater amberjack juveniles under different salinity stresses (40, 30, 20, and 10 ppt) for 30 days to assess their tolerance, adaptation, and molecular responses to salinity. RNA sequencing analysis of gill tissue was used to identify genes and biological processes involved in greater amberjack response to salinity stress at 40, 30, and 20 ppt. Eighteen differentially expressed genes (DEGs) (nine upregulated and nine downregulated) were identified in the 40 vs. 30 ppt group. Moreover, 417 DEGs (205 up-regulated and 212 down-regulated) were identified in the 20 vs. 30 ppt group. qPCR and transcriptomic analysis indicated that salinity stress affected the expression of genes involved in steroid biosynthesis (ebp, sqle, lss, dhcr7, dhcr24, and cyp51a1), lipid metabolism (msmo1, nsdhl, ogdh, and edar), ion transporters (slc25a48, slc37a4, slc44a4, and apq4), and immune response (wnt4 and tlr5). Furthermore, KEGG pathway enrichment analysis showed that the DEGs were enriched in steroid biosynthesis, lipids metabolism, cytokine–cytokine receptor interaction, tryptophan metabolism, and insulin signaling pathway. Therefore, this study provides insights into the molecular mechanisms of marine fish adaptation to salinity.

Characterization and analysis of transcriptome complexity using SMRT-Seq combined with RNA-Seq for a better understanding of Acanthogobius ommaturus in response to temperature stress

Acanthogobius ommaturus is a eurythermic fish, which is widely distributed in coastal, estuarine and bay waters of China, Japan and Korea. Due to the lack of whole genomic information, full-length transcriptome of A. ommaturus was firstly generated by single molecule real-time sequencing (SMRT-seq) in this study. A total of 49,833 full-length non-redundant transcripts (FLNRTs), 2255 alternative splices, 46,856 simple sequence repeats, 5094 long non-coding RNAs and 2708 transcription factors were obtained. Additionally, FLNRTs were used as reference sequences for the following transcriptome analysis of the temperature stress (7 °C, 14 °C, 21 °C (control), 28 °C and 35 °C). GO and KEGG enrichment analysis using GSEA were performed on all genes in 10 response modules which were screened out by WGCNA. Enrichment analysis showed that protein degradation, immune response and energy metabolism play an active role in the temperature stress of A. ommaturus. The differentially expressed hub genes (DEHGs) in response modules were closely related to adhesion, vascular remodeling and disease. The results of this study provided the first systematical full-length transcriptome profile of A. ommaturus and characterized its temperature stress responses, which will serve as the foundation for further exploring the molecular mechanism of the temperature stress in fish.

Gills full-length transcriptomic analysis of osmoregulatory adaptive responses to salinity stress in Coilia nasus

Salinity changes will threaten the survival of aquatic animals. However, osmoregulatory mechanism of Coilia nasus has not been explored. Oxford Nanopore Technologies (ONT) sequencing was performed in C. nasus gills during hypotonic and hyperosmotic stress. 23.8 G clean reads and 27,659 full-length non-redundant sequences were generated via ONT sequencing. Alternative splicing, alternative polyadenylation, transcript factors, and long noncoding RNA were identified. During hypotonic stress, 58 up-regulated differentially expressed genes (DEGs) and 36 down-regulated DEGs were identified. During hypertonic stress, 429 up-regulated DEGs and 480 down-regulated DEGs were identified. These DEGs were associated with metabolism, cell cycle, and transport. The analysis of these DEGs indicated that carbohydrate and fatty acid metabolism were activated to provide energy for cell cycle and transport during hypotonic and hypertonic stress. Cell cycle was also promoted during hypotonic and hypertonic stress. To resist hypotonic stress, polyamines metabolism, ion absorption and water transport from extra-cellular to intra-cellular were promoted, while ion secretion was inhibited. During hypotonic stress, glutamine, alanine, proline, and inositol metabolism were activated. Ion absorption and water transport from intra-cellular to extra-cellular were inhibited. Moreover, different transcript isoforms generated from the same gene performed different expression patterns during hypotonic and hypertonic stress. These findings will be beneficial to understand osmoregulatory mechanism of Coilia nasus.

RNA-seq analyses of Marine Medaka (Oryzias melastigma) reveals salinity responsive transcriptomes in the gills and livers

Increasing salinity levels in marine and estuarine ecosystems greatly influence developmental, physiological and molecular activities of inhabiting fauna. Marine medaka (Oryzias melastigma), a euryhaline research model, has extraordinary abilities to survive in a wide range of aquatic salinity. To elucidate how marine medaka copes with salinity differences, the responses of Oryzias melastigma after being transferred to different salt concentrations [0 practical salinity units (psu), 15 psu, 30 psu (control), 45 psu] were studied at developmental, histochemical and transcriptome levels in the gill and liver tissues. A greater number of gills differentially expressed genes (DEG) under 0 psu (609) than 15 psu (157) and 45 psu (312), indicating transcriptomic adjustments in gills were more sensitive to the extreme hypotonic environment. A greater number of livers DEGs were observed in 45 psu (1,664) than 0 psu (87) and L15 psu (512), suggesting that liver was more susceptible to hypertonic environment. Further functional analyses of DEGs showed that gills have a more immediate response, mainly in adjusting ion balance, immune and signal transduction. In contrast, DEGs in livers were involved in protein synthesis and processing. We also identified common DEGs in both gill and liver and found they were mostly involved in osmotic regulation of amino sugar and nucleotide sugar metabolism and steroid biosynthesis. Additionally, salinity stresses showed no significant effects on most developmental and histochemical parameters except increased heartbeat with increasing salinity and decreased glycogen after transferred from stable conditions (30 psu) to other salinity environments. These findings suggested that salinity-stress induced changes in gene expressions could reduce the effects on developmental and histochemical parameters. Overall, this study provides a useful resource for understanding the molecular mechanisms of fish responses to salinity stresses.

Transcriptome Profiling Reveals a Divergent Adaptive Response to Hyper- and Hypo-Salinity in the Yellow Drum, Nibea albiflora

Simple Summary

Global warming and certain climate disasters (typhoon, tsunami, etc.) can lead to fluctuation in seawater salinity that causes salinity stress in fish. The aim of this study was to investigate the functional genes and relevant pathways in response to salinity stress in the yellow drum. Genes and pathways related to signal transduction, osmoregulation, and metabolism may be involved in the adaptive regulation to salinity in the yellow drum. Additionally, the genes under salinity stress were mainly divided into three expression trends. Our results provided novel insights into further study of the salinity adaptability of euryhaline fishes.

Abstract

The yellow drum (Nibea albiflora) is an important marine economic fish that is widely distributed in the coastal waters of the Northwest Pacific. In order to understand the molecular regulatory mechanism of the yellow drum under salinity stress, in the present study, transcriptome analysis was performed under gradients with six salinities (10, 15, 20, 25, 30, and 35 psu). Compared to 25 psu, 907, 1109, 1309, 18, and 243 differentially expressed genes (DEGs) were obtained under 10, 15, 20, 30, and 35 psu salinities, respectively. The differential gene expression was further validated by quantitative real-time PCR (qPCR). The results of the tendency analysis showed that all DEGs of the yellow drum under salinity fluctuation were mainly divided into three expression trends. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that the PI3K-Akt signaling pathway, Jak-STAT signaling pathway as well as the glutathione metabolism and steroid biosynthesis pathways may be the key pathways for the salinity adaptive regulation mechanism of the yellow drum. G protein-coupled receptors (GPCRs), the solute carrier family (SLC), the transient receptor potential cation channel subfamily V member 6 (TRPV6), isocitrate dehydrogenase (IDH1), and fructose-bisphosphate aldolase C-B (ALDOCB) may be the key genes in the response of the yellow drum to salinity stress. This study explored the transcriptional patterns of the yellow drum under salinity stress and provided fundamental information for the study of salinity adaptability in this species.