Transcriptome profiling of the low-salinity stress responses in the gills of the juvenile Pseudopleuronectes yokohamae

Interplay of gene expression and regulators under salinity stress in gill of Labeo rohita

BACKGROUND

Labeo rohita is the most preferred freshwater carp species in India. The concern of increasing salinity concentration in freshwater bodies due to climate change may greatly impact the aquatic environment. Gills are one of the important osmoregulatory organs and have direct contact with external environment. Hence, the current study is conducted to understand the gill transcriptomic response of L. rohita under hypersalinity environment.

RESULTS

Comprehensive analysis of differentially expressed long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and mRNAs was performed in gills of L. rohita treated with 2, 4, 6 and 8ppt salinity concentrations. Networks of lncRNA-miRNA-mRNA revealed involvement of 20, 33, 52 and 61 differentially expressed lncRNAs, 11, 13, 26 and 21 differentially expressed miRNAs in 2, 4, 6 and 8ppt groups between control and treatment respectively. These lncRNA-miRNA pairs were regulating 87, 214, 499 and 435 differentially expressed mRNAs (DE mRNAs) in 2, 4, 6 and 8ppt treatments respectively. Functional analysis of these genes showed enrichment in pathways related to ion transportation and osmolyte production to cope with induced osmotic pressure due to high salt concentration. Pathways related to signal transduction (MAPK, FOXO and phosphatidylinositol signaling), and environmental information processing were also upregulated under hypersalinity. Energy metabolism and innate immune response pathways also appear to be regulated. Protein turnover was high at 8ppt as evidenced by enrichment of the proteasome and aminoacyl tRNA synthesis pathways, along with other enriched KEGG terms such as apoptosis, cellular senescence and cell cycle.

CONCLUSION

Altogether, the RNA-seq analysis provided valuable insights into competitive endogenous (lncRNA-miRNA-mRNA) regulatory network of L. rohita under salinity stress. L. rohita is adapting to the salinity stress by means of upregulating protein turnover, osmolyte production and removing the damaged cells using apoptotic pathway and regulating the cell growth and hence diverting the essential energy for coping with salinity stress.

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 Revealed Osmoregulation Related Regulatory Networks and Hub Genes in the Gills of Hilsa shad, Tenualosa ilisha, during the Migratory Osmotic Stress

Tenualosa ilisha (Hilsa shad), an anadromous fish, usually inhabits coastal and estuarine waters, and migrates to freshwater for spawning. In this study, large-scale gill transcriptome analyses from three salinity regions, i.e., fresh, brackish and marine water, revealed 3277 differentially expressed genes (DEGs), out of which 232 were found to be common between marine vs freshwater and brackish vs freshwater. These genes were mapped into 54 KEGG Pathways, and the most significant of these were focal adhesion, adherens junction, tight junction, and PI3K-Akt signaling pathways. A total of 24 osmoregulatory genes were found to be differentially expressed in different habitats. The gene members of slc16 and slc2 families showed a dissimilar pattern of expressions, while two claudin genes (cldn11 & cldn10), transmembrane tm56b, and voltage-gated potassium channel gene kcna10 were downregulated in freshwater samples, as compared to that of brackish and marine environment. Protein-protein interaction (PPI) network analysis of 232 DEGs showed 101 genes to be involved in PPI, while fn1 gene was found to be interacting with the highest number of genes (36). Twenty-five hub genes belonged to 12 functional groups, with muscle structure development with seven genes, forming the major group. These results provided valuable information about the genes, potentially involved in the molecular mechanisms regulating water homeostasis in gills, during migration for spawning and low-salinity adaptation in Hilsa shad. These genes may form the basis for the bio-marker development for adaptation to the stress levied by major environmental changes, due to hatchery/culture conditions.

The Effect of Salinity Stress on Enzyme Activities, Histology, and Transcriptome of Silver Carp (Hypophthalmichthys molitrix)

A 56-day study was performed to examine the effect of freshwater (FW) and brackish water (BW 6‱ salinity) on the antioxidant ability, Na+/K+-ATPase (NKA) activities, histology, and transcriptome of the gill and kidney tissue in juvenile silver carp (Hypophthalmichthys molitrix). The results show that when juvenile silver carp were exposed to 6‱ salinity, the activities of superoxide dismutase (SOD) and catalase (CAT) were shown to be substantially increased (p < 0.05), while glutathione peroxidase (GSH-PX) activities in gill were not significantly affected (p < 0.05). In kidney tissue, SOD, CAT, and GSH-PX, enzyme activities peaked at 24, 8, and 4 h, respectively, but were not significantly different compared with the control group (p < 0.05). In addition, significant effects of salinity were observed for the NKA level in both the gills and kidney tissues (p < 0.05). The gill filaments of juvenile silver carp under the BW group all underwent adverse changes within 72 h, such as cracks and ruptures in the main part of the gill filaments, bending of the gill lamellae and enlargement of the gaps, and an increase in the number of mucus and chloride-secreting cells. Transcriptome sequencing showed 171 and 261 genes in the gill and kidney tissues of juvenile silver carp compared to the BW group, respectively. Based on their gene ontology annotations, transcripts were sorted into four functional gene groups, each of which may play a role in salt tolerance. Systems involved in these processes include metabolism, signal transduction, immunoinflammatory response, and ion transport. The above findings indicate that the regulation processes in juvenile silver carp under brackish water conditions are complex and multifaceted. These processes and mechanisms shed light on the regulatory mechanism of silver carp osmolarity and provide a theoretical foundation for future research into silver carp growth in brackish water aquaculture area.

Integrated application of multi-omics approach and biochemical assays provides insights into physiological responses to saline-alkaline stress in the gills of crucian carp (Carassius auratus)

Given the decline of freshwater resources in recent years, the accessible space for freshwater aquaculture is rapidly shrinking, and aquaculture in saline-alkaline water has become a critical approach to meet the rising demand. However, the molecular mechanism behind the adverse effects of saline-alkaline water on fish and the regulatory mechanism in fish tolerance remains unclear. Here, adult crucian carp (Carassius auratus) were exposed to 60 mmol/L NaHCO₃ for 30 days. It was observed that long-term carbonate alkalinity (CA) exposure not only caused gill oxidative stress but also changed the levels of several physiological parameters associated with ammonia transport, including blood ammonia, urea nitrogen (BUN), glutamine (Gln), and glutamine synthetase (GS). According to the metabolomics study, differential metabolites (DMs) engaged in various metabolic pathways, such as glycerophospholipid metabolism, sphingolipid metabolism, and arachidonic acid metabolism. In addition, transcriptomics data showed that differentially expressed genes (DEGs) were closely related to ammonia transport, apoptosis, and immunological response. In general, comprehensive multi-omics and biochemical analysis revealed that crucian carp might adopt Rh glycoprotein as a carrier to mediate ammonia transport and increase glutamine and urea synthesis under long-term high saline-alkaline stress to mitigate the adverse effects of blocked ammonia excretion. Simultaneously, saline-alkaline stress caused the destruction of the antioxidant system and the disorder of lipid metabolism in the crucian carp gills, which induced apoptosis and immunological response. To our knowledge, this is the first study to investigate fish's molecular and metabolic mechanisms under saline-alkaline stress using integrated metabolomics, transcriptomics, and biochemical assays. Overall, the results of this study provided new insights into the molecular mechanism behind the adverse effects of saline-alkaline water on fish and the regulatory mechanism in fish tolerance.

Transcriptome Analysis Reveals the Molecular Response to Salinity Challenge in Larvae of the Giant Freshwater Prawn Macrobrachium rosenbergii

Salinity is a crucial factor influencing the growth, development, immunity, and reproduction of aquatic organisms; however, little is known about the molecular mechanism of the response to salinity challenge in larvae of the giant freshwater prawn Macrobrachium rosenbergii. Herein, larvae cultured in three treatment groups with salinities of 10, 13, and 16‰ (S10, S13, and S16) were collected, and then transcriptome analysis was conducted by RNA-seq. A total of 6,473, 3,830 and 3,584 differentially expressed genes (DEGs) were identified in the S10 vs. S13 comparison, S10 vs. S16 comparison and S13 vs. S16 comparison, respectively. These genes are involved in osmoregulation, energy metabolism, molting, and the immune response. qPCR analysis was used to detect the expression patterns of 16 DEGs to verify the accuracy of the transcriptome data. Protein–protein interaction (PPI) analysis for DEGs and microsatellite marker screening were also conducted to reveal the molecular mechanism of salinity regulation. Together, our results will provide insight into the molecular genetic basis of adaptation to salinity challenge for larvae of M. rosenbergii.

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

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.

Gill physiological and transcriptomic response of the threatened freshwater mussel Solenaia oleivora to salinity shift

Solenaia oleivora, a freshwater shellfish endemic to China, is becoming one of the most threatened freshwater mussels owing to water pollution, habitat fragmentation, and overfishing. Hence, exploring its response to different environmental factors is important for its conservation. In this work, we investigated the physiological and transcriptomic response of S. oleivora to increased salinity. We found that increased salinity caused the death of S. oleivora. High salinity caused shrinking and deformation of gill filaments, reduced gill cilia, and induced cell apoptosis in gills. The activities of superoxide dismutase (SOD), catalase (CAT), acid phosphatase (ACP), alkaline phosphatase (AKP), as well as glutathione (GSH) content were increased at the beginning of salinity stress (3-12 h), while SOD and ACP activities decreased at 48 h. Transcriptome data revealed that high salinity stress (48 h) induced 766 differentially expressed genes (DEGs). Among these DEGs, the majority of the stress response and ion transport-related genes were up-regulated, while most of the immune-related genes were down-regulated. In conclusion, these findings suggest that the antioxidant and immune functions of S. oleivora can be inhibited by high salinity, which may be one of the main reasons for its low survival rate under conditions of increasing salinity.

Differential Expression of the Apolipoprotein AI Gene in Spotnape Ponyfish (Nuchequula nuchalis) Inhabiting Different Salinity Ranges at the Top of the Estuary and in the Deep-Bay Area of Gwangyang Bay, South Korea

Spotnape ponyfish (Nuchequula nuchalis) is a dominant species that is broadly distributed from estuarine to deep-bay areas, reflecting a euryhaline habitat. Apolipoprotein AI (ApoAI) is a main component of plasma lipoproteins and has crucial roles in lipid metabolism and the defense immune system. In this study, we characterized the N. nuchalis ApoAI gene and analyzed the expression of the ApoAI transcript in N. nuchalis collected at various sites in the estuary and the deep-bay area which have different salinities. Owing to the fish’s mobility, we conducted stable isotope analyses to confirm the habitat characteristics of N. nuchalis. Carbon and nitrogen isotope ratios (δ¹³C and δ¹⁵N) from N. nuchalis indicated different feeding sources and trophic levels in the estuarine and deep-bay habitats. The characterized N. nuchalis ApoAI displayed residual repeats that formed a pair of alpha helices, indicating that the protein belongs to the apolipoprotein family. In the phylogenetic analysis, there was no sister group of N. nuchalis ApoAI among the large clades of fish species. The transcriptional expression level of ApoAI was higher in N. nuchalis inhabiting the deep-bay area with a high salinity (over 31 psu) than in N. nuchalis inhabiting the top of the estuary with a low salinity (6~15 psu). In addition, the expression patterns of N. nuchalis ApoAI were positively correlated with environmental factors (transparency, pH, TC, and TIC) in the high salinity area. These results suggest that ApoAI gene expression can reflect habitat characteristics of N. nuchalis which traverse broad salinity ranges and is associated with functional roles of osmoregulation and lipid metabolism for fish growth and development.

Role of Fucoidan on the Growth Behavior and Blood Metabolites and Toxic Effects of Atrazine in Nile Tilapia Oreochromis niloticus (Linnaeus, 1758)

Simple Summary

Toxic derivatives reach the ponds and cages where fish are grown, and the continuous exposure to these contaminants proved to impair the healthy status of several finfish species. In some countries famous for cultivating rice and corn, atrazine (ATZ) is massively applied to protect plants from invaders. Many functional additives are permitted for application in the aquaculture sector as natural alternatives for chemotherapies. In this study, the toxicity impacts of ATZ and the protective role of fucoidan were investigated on the health performance of Nile tilapia. Long-term exposure to ATZ resulted in low growth rate, impaired hepato-renal function, intestinal inflammation, and oxidative stress in Nile tilapia. However, the obtained results soundly support fucoidan’s potential role to cope with the impacts of ATZ on Nile tilapia.

Abstract

Waterborne herbicides are stressful agents that threaten the productivity and safety of finfish species. In this study, the toxicity impacts of atrazine (ATZ) and the protective role of fucoidan were investigated on the health performance of Nile tilapia. For 40 days, the total number of 180 Nile tilapia was assigned in four groups (triplicates, 15 fish/replicate), where the first (control) and third groups were offered the control diet, while the second and fourth groups were offered a fucoidan (FCN). Further, in the third and fourth groups, the water was mixed with atrazine (ATZ) at 1.39 mg/L daily. The growth rate, FCR, and survival rate were markedly enhanced by fucoidan but severely declined by ATZ exposure (p < 0.05). The morphological structure of the intestine in the control fish revealed normal structure, while fucoidan-treated groups showed eminent enhancement and branching of the intestinal villi. The intestine of ATZ-treated fish revealed deterioration and the intestinal mucosa, inflammatory cell infiltration, and separation of lining epithelium. The highest Hb, PCV, RBCs, WBCs, total protein, and albumin were observed in Nile tilapia fed fucoidan, but the worst values were seen in ATZ-intoxicated fish (p < 0.05). The liver-related enzymes (ALT and AST) and kidney function (urea and creatinine) showed impaired values by ATZ toxicity and were regulated by dietary fucoidan. Meanwhile, fish fed fucoidan and exposed to ATZ had lower total cholesterol and triglyceride values than fish exposed to ATZ without fucoidan feeding (p < 0.05). The SOD, CAT, GPx, cortisol, and glucose levels were increased in ATZ-exposed fish and reduced by fucoidan (p < 0.05). However, the level of malondialdehyde (MDA) was reduced in fucoidan-fed fish and increased in ATZ-exposed fish (p < 0.05). Altogether, dietary fucoidan is required in fish diets to alleviate the impacts of ATZ-induced toxicity.