Skip to main content

Advertisement

Springer Nature is making Coronavirus research free. View research | View latest news | Sign up for updates

Transcriptomic response to three osmotic stresses in gills of hybrid tilapia (Oreochromis mossambicus female × O. urolepis hornorum male)

Abstract

Background

Osmotic stress is a widespread phenomenon in aquatic animal. The ability to cope with salinity stress and alkaline stress is quite important for the survival of aquatic species under natural conditions. Tilapia is an important commercial euryhaline fish species. What’s more tilapia is a good experimental material for osmotic stress regulation research, but the molecular regulation mechanism underlying different osmotic pressure of tilapia is still unexplored.

Results

To elucidate the osmoregulation strategy behind its hyper salinity, alkalinity and salinity-alkalinity stress of tilapia, the transcriptomes of gills in hybrid tilapia (Oreochromis mossambicus ♀ × O. urolepis hornorum ♂) under salinity stress (S: 25‰), alkalinity stress(A: 4‰) and salinity-alkalinity stress (SA: S: 15‰, A: 4‰) were sequenced using deep-sequencing platform Illumina/HiSeq-2000 and differential expression genes (DEGs) were identified. A total of 1958, 1472 and 1315 upregulated and 1824, 1940 and 1735 downregulated genes (P-value < 0.05) were identified in the salt stress, alkali stress and saline-alkali stress groups, respectively, compared with those in the control group. Furthermore, Kyoto Encyclopedia of Genes and Genomes pathway analyses were conducted in the significant different expression genes. In all significant DEGs, some of the typical genes involved in osmoregulation, including carbonic anhydrase (CA), calcium/calmodulin-dependent protein kinase (CaM kinase) II (CAMK2), aquaporin-1(AQP1), sodium bicarbonate cotransporter (SLC4A4/NBC1), chloride channel 2(CLCN2), sodium/potassium/chloride transporter (SLC12A2 / NKCC1) and other osmoregulation genes were also identified. RNA-seq results were validated with quantitative real-time PCR (qPCR), the 17 random selected genes showed a consistent direction in both RNA-Seq and qPCR analysis, demonstrated that the results of RNA-seq were reliable.

Conclusions

The present results would be helpful to elucidate the osmoregulation mechanism of aquatic animals adapting to saline-alkali challenge. This study provides a global overview of gene expression patterns and pathways that related to osmoregulation in hybrid tilapia, and could contribute to a better understanding of the molecular regulation mechanism in different osmotic stresses.

Background

Saline-alkali water accounts for a considerable proportion of water resources in the world. Salinity has been long recognized as one of the fundamental factors affecting aquatic species distribution and influencing physiological processes of marine and estuarine organisms, such as survival, hemolymph osmolarity, and tissue water content [1,2,3]. Furthermore, salinity has significant effects on physiology of aquatic organisms, and salinity adaptation is a complicated process that involves a series of physiological responses to the environment with different osmotic regulation requirements [4, 5]. Carbonate alkalinity stress was considered as a major risk factor for fishes surviving in saline-alkaline water [6, 7]. Although many studies have examined osmoregulatory physiology in freshwater and marine populations of fish [8,9,10,11,12], it still remains unknown about the molecular mechanism of osmotic stress tolerance [13]. Therefore, the identification and characterization of the genes and the regulatory factors involved in hyper-osmoregulation are now essential for increasing the production and the efficiency of selective breeding programs for some important fish species.

Tilapia is an important commercial fish in China, which is euryhaline fish species and it is the second most important fish after carps in aquaculture [14]. Owing to ease of aquaculture, marketability and stable market prices, the production of tilapia has quadrupled over the past decade [15], and which is promoted and farmed in more than 100 countries and regions [16]. Owing to its importance in tropical and subtropical aquaculture and its extreme euryhaline ability, tilapia has received a high level of scientific interest and is one of the most popular model species for research on fish osmoregulation [17,18,19]. With the increasing scarcity of freshwater available for aquaculture, tilapia that can live in brackish or hyper-osmotic water would enable expanding the range of culture and increasing the production of global tilapia [20]. In China one new hybrid tilapia strain, named Mo-Ho tilapia “Guangfu No.1”, is a hybrid from Mozambique tilapia Oreochromis mossambicus female × Wami tilapia O. urolepis hornorum male and is tolerant of hyper-osmotic water and grows well in saline pools [21]. It has been widely farmed on the southern coast of China.

With regard to ion-regulation, a set of tissues and organs such as the gills, kidney and gut plays a vital role in the teleosts [22,23,24]. Gill contains complex transport epithelia functions as aquatic gas exchange, acid-base regulation and excretion of nitrogenous wastes [25]. Several studies have shown that the epithelium of gill is reorganized, the distribution, number and size of chloride cells both have changed, considerably to meet the requirements of ion transport and permeability [26, 27]. When subjected to various abiotic stresses from ambient water, such as salt and alkali, gills are the first site for sensing stress signals and initiating signaling cascades at the molecular level in response to adverse environments. Through these responses, fish can regulate gene expression of regulatory, functional proteins, and morphology of gill epithelium [28,29,30] to enhance stress tolerance. It is important for current fish biology researches to study fish response to osmotic stress.

Recent years, with the prosperous development of the next-generation sequencing, high-throughput sequencing has been widely used to identify conserved and novel functional genes and signaling pathways in aquatic animals, such as Litopenaeus vannamei [13, 31], Eriocheir sinensis [32], Palaemon caridean shrimps [33] and Mytilus trossulus [34]. Guo et al identified three possible osmoregulation-related signaling pathways as lipid metabolism related pathways, tight junction pathway and thyroid hormone signaling pathway in Siberian sturgeon Acipenser baeri in response to salinity stress [35]. Lv et al showed that 552 differentially expressed genes including amino acid transport proteins and ion transport enzymes were obtained in Portunus trituberculatus under low salinity stress [36].

In this study, the transcriptome of the hybrid tilapia (O. mossambicus♀ × O. urolepis hornorum♂) gills was sequenced by using Illumina/Hiseq-2000 RNA-seq technology for the first time. The transcriptomic data between control group and salinity challenge group, alkalinity challenge group and salinity-alkalinity challenge group were compared and analyzed to identify osmoregulation-related genes and pathways. The discoveries of this study can be conductive to illustrate the mechanism of osmotic regulation for aquatic animals adapting to salinity challenges, alkalinity challenge and salinity-alkalinity challenge.

Results

Sequencing and the analysis of reads

To identify mRNAs in hybrid tilapia in response to osmotic stress challenges, mRNA libraries derived from 12 groups were constructed and sequenced using the Illumina deep-sequencing platform. Using high-throughput sequencing, 54,418,529, 61,482,829, 48,353,150 and 44,867,651 average raw reads were obtained from gill tissue of control group(C), saline challenge group (S), alkali challenge group (A) and saline-alkali challenge group (SA) in hybrid tilapia, respectively (Table 1). After quality assessment, the control and challenge groups yielded 53,886,140, 60,565,527, 47,545,806 and 44,505,885 average clean reads with a Q30 average percentage 90.99, 88.63, 88.76 and 91.73%, respectively. After mapping of these clean reads and the reference genome of Oreochromis niloticus (https://www.ncbi.nlm.nih.gov/assembly/GCA_001858045.3), average of three sample size (n = 3) of four groups were 38,413,422 (71.63%), 41,942,362 (69.23%), 33,426,853 (70.45%) and 32,895,304 (73.95%) were validated as the unique matches.

Table 1 Number of high-throughput raw reads, clean reads and mapped clean reads generated from tilapia gills mRNA library

Differential gene expression (DGE) analysis

Based on the expression patterns of DEGs in this study, the values of log10FPKM (FPKM: fragments per kilo bases per million reads) were used to cluster the 12 groups, respectively, to measure and identify the similarities among the groups (Fig. 1). Based on the criteria that |logFC| > 1, P-value < 0.05 and FDR ≤ 0.05, the edegR and DESeq2 method were used to conducted differential expression analysis, the difference between these two methods were not significant. According to Sahraeian et al study [37], DESeq2 provided the most accurate differential analysis, so the analysis results of DESeq2 method were used to conduct the following analysis (Fig. 2). 3772 significant DEGs were detected in group S compared with control group, of which there were 1958 up-regulated genes and 1814 down-regulated genes in group S (Additional file 1: Fig. S1A). We also identified 3412 significant DEGs in group A compared with control group including 1472 up-regulated genes and 1940 down-regulated genes in group A (Additional file 1: Fig. S1B). Simultaneously, 3050 significant DEGs were detected in group SA compared with group C, of which there were 1315 up-regulated genes and 1735 down-regulated genes in group SA (Additional file 1: Fig. S1C). In addition, 2504 significant DEGs were detected in group SA compared with group S, of which there were 940 up-regulated genes and 1564 down-regulated genes in group SA (Additional file 1: Fig. S1D). Meanwhile, 609 differentially expressed genes in group SA compared with group A, which include 312 up-regulated genes and 297 down-regulated genes in group SA (Additional file 1: Fig. S1E). Further analysis indicated that a total of 1250 unigenes showed significantly different expression levels in salinity stress (S), alkalinity stress (A) and salinity-alkalinity stress (SA) (Fig. 3a). And a total of 46 unigenes showed significantly different expression levels in the C vs. S, C vs. A, C vs. SA, S vs. SA and A vs. SA groups, respectively (Fig. 3b).

Fig. 1
figure1

Hierarchical cluster analysis of differentially expressed genes based on the data of log ratio fold change. Heat map showing differentially expressed mRNAs upon osmotic stress in hybrid tilapia. The expression heat map was generated with the value of log10FPKM. The color scale shows the levels of differentially expressed genes: red color indicates enhanced expression of mRNA and blue color indicates decrease in expression levels of mRNA

Fig. 2
figure2

The number of differentially expressed DEGs in different groups

Fig. 3
figure3

DEGs number and venn diagram of overlap of the different groups. a: DEGs number and venn diagram of overlap of the C vs. S, C vs. A and C vs. SA groups. b: DEGs number and venn diagram of overlap of the C vs. S, C vs. A, C vs. SA, S vs. SA and A vs. SA groups

Gene ontology (GO) analysis of significant DEGs

In this study, GO functional analysis indicated 2300 (C vs. S), 2110 (C vs. A) and 1960 (C vs. SA) significant genes in total were classified into 40, 35 and 37 sub-categories of three major categories: biological processes, cellular components and molecular function, respectively (Fig. 4). Through the comparative analysis of GO enrichment results between different groups, it was found that some significant GO terms, such as protein binding transcription factor activity (GO:0000988), behavior (GO:0007610), signaling (GO:0023052) and molecular function regulator (GO:0098772) were enriched uniquely in salt tolerant. While in salt-alkalinity group, reproductive process (GO:0022414) and cell aggregation (GO:0098743) were uniquely enriched. We also found some significant GO terms are most likely to play an essential role in regulating osmotic stress tolerance in hybrids, such as immune system process (GO:0002376), response to stimulus (GO:0050896), transporter activity (GO:0005215) and electron carrier activity (GO:0009055). And these results of GO enrichments will provide the research clues for our following research.

Fig. 4
figure4

Gene Ontology (GO) classification of assembled unigenes. a: C vs. S; b: C vs. A; c: C vs. SA; d: S vs. SA; e: A vs. SA

Kyoto encyclopedia of genes and genomes analysis

The top 30 KEGG pathways with the most number of annotated sequences were listed (Additional file 1: Table S1, S2, S3). These pathways consisted of cytokine-cytokine receptor interaction, protein processing in endoplasmic reticulum, autoimmune thyroid disease, pyrimidine metabolism, amino sugar and nucleotide sugar metabolism, p53 signaling pathway, biosynthesis of amino acids, glycolysis / gluconeogenesis, DNA replication, proximal tubule bicarbonate reclamation, proteasome, plycine, serine and threonine metabolism, alanine, aspartate and glutamate metabolism, TNF signaling pathway, IL-17 signaling pathway, and so on.

Gene set enrichment analysis (GSEA) analysis of genes

To date, GSEA is being used to identify specific biological processes involved in disease outcomes and bioinformatics analysis in the medical literature and the biological research. GSEA enrichment tool was used to analysis genes between different groups, and then GO analysis was performed to find which have biological regulate function gene set. After GSEA analysis in this study, multiple GO terms were enriched in salinity group, alkalinity group and salt-alkalinity group (Additional file 1: Table S4, S5, S6) and a vast number genes were annotated to take participated in osmotic stress tolerance.

Identification of significant DEGs and KEGG pathways related to osmoregulation

In addition, in all significant DEGs, some functional genes that related to osmoregulation, including ATP1A (sodium / potassium - transporting ATPase subunit alpha), PRLR (prolactin receptor), NKCC1 (sodium / potassium / chloride transporter), SLC9A3, NHE3 (sodium / hydrogen exchanger), TRPV4 (transient receptor potential cation channel subfamily V member 4), SLC4A2, AE2 (anion exchanger) and CLCN2 (chloride channel 2), were also identified separately. Some genes that may play important roles in osmoregulation in salinity stress, alkalinity stress and salinity-alkalinity stress were listed in Table 2.

Table 2 Osmoregulation-related differentially expressed genes (DEGs) regulated after stressed in group C vs. S, C vs. A and C vs. SA

A comparative analysis of the KEGG pathways enriched in the three groups revealed that 25 of the pathways were significantly enriched in all the three groups, including steroid biosynthesis (ko00100), IL-17 signaling pathway (ko04657), glycine, serine and threonine metabolism (ko00260), Cell cycle (ko04110), and so on. However, only two pathways, proteasome (ko03050) and the p53 signaling pathway (ko04115), were enriched in the salinity group and the alkalinity group. Meanwhile, three pathways, fructose and mannose metabolism (ko00051), glycolysis / gluconeogenesis (ko00010) and alanine, aspartate and glutamate metabolism (ko00250), were enriched in the salinity and saline-alkali group. While in alkali stress group and salinity-alkalinity group, there were 18 pathways were significantly enriched, include natural killer cell mediated cytotoxicity (ko04650), synthesis and degradation of ketone bodies (ko00072), glycosaminoglycan degradation (ko00531), and so on (Table 3).

Table 3 Comparative analysis of the KEGG pathways enriched in the three groups

In the present study, 17 osmoregulation-related signaling pathways as steroid biosynthesis (ko00100), mucin type O-glycan biosynthesis (ko00512), cell cycle (ko04110), proximal tubule bicarbonate reclamation (ko04964), glycine, serine and threonine metabolism (ko00260), alanine, aspartate and glutamate metabolism (ko00250), p53 signaling pathway (ko04115), glycolysis / gluconeogenesis (ko00010), glutathione metabolism (ko00480), ECM-receptor interaction (ko04512), protein processing in endoplasmic reticulum (ko04141), arachidonic acid metabolism (ko00590), TNF signaling pathway (ko04668), mineral absorption (ko04978), cell adhesion molecules (CAMs) (ko04514), valine, leucine and isoleucine biosynthesis (ko00290), pantothenate and CoA biosynthesis (ko00770) were showed in Table 3.

Significant DEGs analysis in different groups

Hierarchical cluster, GO function analysis and KEGG pathway analysis were carried out for significant DEGs in each interval in Venn chart (Fig. 3a). Hierarchical clusters were analyzed of significant DEGs that were expression uniquely in salinity stress (S), alkalinity stress (A), salinity-alkalinity stress (SA), salinity stress and alkalinity stress (S & A), salinity stress and salinity-alkalinity stress (S & SA), alkalinity stress and salinity-alkalinity stress (A & SA), the heat map of significant DEGs were showed in Additional file 1: Fig. S2.

Furthermore, GO analyses were performed to identify similarities and differences among differentially expressed genes only in salinity stress (S), alkalinity stress (A), salinity-alkalinity stress (SA), alkalinity and salinity-alkalinity stress (A, SA), salinity and salinity-alkalinity stress (S, SA) and three osmotic stress shared (S, A, SA). By comparing the analysis results of GO enrichments, we found that transcription factor activity, protein binding (GO:0000988) and behavior (GO:0007610) were uniquely enriched in salinity stress (A). We also found that cell aggregation (GO:0098743), multi-organism process (GO:0051704) and membrane-enclosed lumen (GO:0031974) were uniquely enriched in salinity-alkalinity stress (SA), salinity-alkalinity stress (S, SA) and three osmotic stress (S, A, SA), respectively. All these analysis results will provide research data and clues for the future research directions, and are essential for the research of the molecular mechanism of osmotic regulation in teleost. The results of GO enrichment were listed as Fig. 5.

Fig. 5
figure5

Gene Ontology (GO) classification of differential expression genes in group S, A, SA, S & A, S & SA, A & SA S, and A & SA. a: S; b: A; c: SA; d: S & A; e: S & SA; f: A & SA; g: S, A & SA

Meanwhile, KEGG analyses were also performed to identify similarities and differences among differentially expressed genes uniquely in salinity stress (S), alkalinity stress (A), salinity-alkalinity stress (SA), alkalinity and salinity-alkalinity stress (A & SA), salinity and salinity-alkalinity stress (S & SA) and three osmotic stress shared (S & A & SA) groups. A total of 133, 226 and 58 significant DEGs of only expressed in S, A and SA were annotated in KEGG and grouped into 12, 17 and 13 pathway terms, respectively (Tables 4, 5, and 6). Furthermore, 15, 12, 27 and 29 pathway terms were annotated in S & A, S & SA, A & SA and S & A & SA groups, respectively (Tables 7, 8, 9, and 10). We found that cytokine-cytokine receptor interaction, ECM-receptor interaction, glycolysis / gluconeogenesis, ether lipid metabolism, linoleic acid metabolism, carbon fixation in photosynthetic organisms, beta-alanine metabolism, alpha-linolenic acid metabolism, quorum sensing and styrene degradation pathway were uniquely annotated in salinity stress. However, the pathways enriched only in alkalinity stress are mostly related to immunity and disease. Under the salinity-alkalinity stress, mostly pathways were enriched related to metabolism and immunity. The results suggested that compared with salinity stress, alkalinity stress makes tilapia more susceptible to pathogen infection, which more likely leads to diseases.

Table 4 KEGG pathway analysis of DEGs expressed only in group S, not in other groups
Table 5 KEGG pathway analysis of DEGs expressed only in group A, not in other groups
Table 6 KEGG pathway analysis of DEGs expressed only in group SA, not in other groups
Table 7 KEGG pathway analysis of DEGs expressed only in group S & A, not in other groups
Table 8 KEGG pathway analysis of DEGs expressed only in group S & SA, not in other groups
Table 9 KEGG pathway analysis of DEGs expressed only in group A & SA, not in other groups
Table 10 KEGG pathway analysis of DEGs expressed only in group S & A & SA, not in other groups

The validation of differently expressional genes by qRT-PCR

To validate the veracity and reliability of differentially expressed genes identified by mRNA-Seq, we randomly selected 17 genes for qRT-PCR validation from those with different expression patterns based on functional enrichment and pathway results. Melting-curve analysis revealed a single product for all tested genes. Log2FCs from qRT-PCR were compared with the mRNA-Seq expression analysis results (Fig. 6). The results of qRT-PCR were significantly correlated with the mRNA-Seq results with a correlation coefficient of 0.88 (P < 0.01) demonstrating the credibility of the mRNA-Seq results.

Fig. 6
figure6

Validation of RNA-seq data by qPCR. To validate the RNA-seq data, the relative mRNA levels of 17 randomly selected DEGs in the gills of hybrid tilapia were examined by qPCR. The mRNA levels by qPCR are presented as the fold change compared with the mock-treated control after normalization against β -actin. The relative expression levels from the RNA-seq analysis were calculated as log2FC values. The gene nkai1 was calculated in two different groups. (Note: nkain1–1 was in a group C vs. A, nkain1–2 was in group C vs. SA; trpv4–1 was in group S vs. SA, trpv4–2 was in group C vs. S)

Discussion

Osmoregulation in tilapia is a complex process because of the diverse range of salinities and alkalinities that they are exposed to in their natural habitat. To be euryhaline, tilapia must be able to cope with salt gain and osmotic water loss when it is in hyper-osmotic water stress. Edwards et al discussed the principles of ion and water transport in euryhaline fishes, and discussed the role of including gills, kidney, urinary bladder, gastrointestinal tract, rectal gland of elasmobranchs, skin, opercular membrane, and yolk sac organs in regulated the osmotic stress [38]. N’Golo et al reported the changes of ionocyte morphology and function after transfer from fresh to hypersaline waters in the tilapia Sarotherodon melanotheron [39].

High-throughput RNA-seq is a good method in illuminating the underlying molecular mechanisms of osmoregulation and immunization [40, 41], and the transcriptomic analysis with osmotic stress has been realized in some aquatic vertebrate species and invertebrate species, such as Oreochromis mossambicus [42], Acipenser baeri [35], Oryzias latipes [43], Litopenaeus vannamei [13, 31], Eriocheir sinensis [32], Crassostrea gigas [44] and Mytilus trossulus [34]. However, most of these studies were directed to only one osmotic stress treatment, but the natural water environment is complex and changing, so it is important to perform various osmotic stress treatments on fish. In this study, we conducted high-throughput sequencing of three osmotic stresses, including saline stress, alkaline stress and mixed of saline-alkaline stress on hybrid tilapia, in order to understand the osmoregulation mechanism of tilapia under different osmotic stress.

In the present study, 3772, 3412 and 3050 significantly different expression genes were identified in salinity stress, alkalinity stress and saline-alkaline stress, respectively. GO and KEGG analysis were conducted to identified genes and pathways that related to osmoregulation. GO analysis results indicated that cellular process, cell part and binding were the dominant group in the category of biological processes, cellular components and molecular functioning, respectively. These results were consistent with the researches in aquatic animals such as Pacific White Shrimp Litopenaeus vannamei [31], Siberian sturgeon Acipenser baeri [35], Eriocheir sinensis [32] and Ruditapes philippinarum [39]. Numerous ion transport enzymes and ion transporters associated with osmoregulation were identified in the three osmotic stress groups. Osmoregulation related genes, such as NKA, PRL, PRLR, TRPV4, AQP1, NHE3, and solute carrier family have been reported in many studies. Zhu et al found that the NKAα1 mRNA expression levels in the gills of tilapia reached peak level at 24 h after transfer from freshwater to seawater [45]. Yoko et al research showed that PRLR2 facilitates acclimation responses to increased extracellular osmolality [46]. Yuan et al observed the synchronous expression trend of the renal PRLR with pituitary PRL (5d) and the asynchronous expression peaks between branchial (8d) and renal PRLR (5d) in olive flounder Paralichthys olivaceus, and the mRNA levels of PRL receptor (PRLR) in both gill and kidney displayed a similar trend to the pituitary PRL [47]. Andre et al reported that TRPV4 expression increased by 48 h after transfer from freshwater to seawater and declined as early as 6 h after transfer from seawater to freshwater, indicating TRPV4 expression maybe regulated through the same second messenger pathways involved in hyposmotically-induced PRL release [48]. Muhammad et al suggested that cytoplasmic carbonic anhydrase (ChqCAc) gene is involved in response to changes in pH and in systemic acid–base balance in freshwater crayfish [49]. Furthermore, gene NHE3 was detected and analyzed in Sciaenops ocellatus [50] and Macrobrachium species [51]. Some other genes slc12a1, NKCC2 and AQP1 were also detected and analyzed in Oryzias latipes after transfer from fresh water to seawater [52]. Tatsuya et al also found Na+ / K+ − ATPase activity increased in specific gills of Octopus ocellatus after transfer from 30-ppt normal seawater to 20-ppt salinity [53].

By matching to KEGG pathway database, possible functions of the significant DEGs were analyzed in order to learn more about the physiological responses of hybrid tilapia to simultaneous salinity stress, alkalinity stress and salinity-alkalinity stress. 39, 51 and 58 KEGG pathways were found in salinity stress, alkalinity stress and salinity-alkalinity stress, respectively. 25 pathways were enriched in the three osmotic stress group. Among these pathways, seven pathways were related to metabolism, eight pathways were related to disease and immunity, and the remaining was related to biosynthesis, cell cycle and proximal tubule bicarbonate reclamation. In deal with hyper-osmotic stress and the energetically demanding salt extrusion, signaling related to DNA replication, mis-match repair, protein biosynthesis and degradation, metabolism and cell cycle regulation were increased [42]. Jiang et al research indicated that amino acid metabolism played crucial roles in osmoregulation [54].

In the proximal tubule bicarbonate reclamation pathway of salinity stress, the gene expression of carbonic anhydrase 4 (CA4), aquaporin-1 (AQP1), sodium-coupled neutral amino acid transporter (SNAT3, solute carrier family 38 member 3), malate dehydrogenase (MDH1) and phosphoenolpyruvate carboxykinase (GTP) (PEPCK) were significantly upregulated. However, the expression of sodium bicarbonate cotransporter (NBC1, solute carrier family 4 member 4) was significantly decreased. While in the alkalinity stress, the expression of NBC1 was significantly up-regulated, but the expression of NHE3 was similar with the control group. However, in salinity-alkalinity stress, the expression of NHE3 was significantly down-regulated, NBC1 gene was up-regulated, but the expression of AQP1 was similar to the control group. Those genes plays pivotal role in maintaining the balance of proximal tubular cell and the regulation of blood pH. The renal proximal tubule is responsible for most of the renal sodium, chloride, and bicarbonate reabsorption [55]. Alexander et al showed that NHE3 participates in the reabsorption of Na+, bicarbonate, and water from the proximal tubule, contributing to the maintenance of intravascular volume, blood pressure, and acid-base homeostasis [56]. Masashi et al found that a majority of renal proximal bicarbonate absorption is accomplished by coordinated operation of the apical NHE3 and the basolateral electrogenic Na+/ HCO3 cotransporter NBC1 [57]. Defects in NHE3 or NBC1 can theoretically reduce the reabsorption of filtered HCO3 and cause proximal renal tubular acidosis (pRTA) [58,59,60,61]. Gordon et al reported that AQP1 enhances the reabsorption of HCO3 by the renal proximal tubule by increasing the CO2 permeability of the apical membrane and the macroscopic cotransport of Na+ plus two HCO3 occurs as NBC transports Na+ plus CO32− and AQP1 transports CO2 and H2O [62]. Our results indicated that proximal tubule bicarbonate reclamation pathway was evidently crucial for maintaining systemic acid-base homeostasis in osmotic stress of tilapia.

Nine pathways were uniquely enriched in salinity stress, such as ECM-receptor interaction, protein processing in endoplasmic reticulum, nitrogen metabolism and arachidonic acid metabolism. Among these pathways, arachidonic acid metabolism playing important role in osmoregulation has been reported in Eriocheir sinensis [32], Litopenaeus vannamei [31], Ruditapes philippinarum [39] and Litopenaeus vannamei [63]. Six KEGG pathways were characteristic enriched in alkalinity stress group. Meanwhile, 12 KEGG pathways, including pantothenate and CoA biosynthesis, mineral absorption, valine, leucine and isoleucine biosynthesis, cell adhesion molecules (CAMs), glyoxylate, dicarboxylate metabolism and propanoate metabolism, etc., were uniquely annotated in salinity-alkalinity stress. In the mineral absorption pathway, TRPV6, chloride channel 2 (clcn2), zinc transporter (ZnT1), Na+ / K+-exchanging ATPase and femitin (FTH1) were upregulated significantly. Meanwhile, NHE3 and heme oxygenase 1 (HMOX1) were downregulated significantly. Through the network formed of all genes, this pathway played vital role in ion regulation and maintain the balance of osmotic stress in internal and external of the tilapia.

From KEGG pathway analysis of significant DEGs uniquely in salinity stress group, three enriched pathways related to lipid metabolism including alpha-linolenic acid metabolism, linoleic acid metabolism and ether lipid metabolism were examined. The long-chain polyunsaturated fatty acid (LC-PUFA) level raised in euryhaline and migratory teleosts during salinity acclimation has been reported by many researches, such as Oncorhynchus mykiss [64], Dicentrarchus labrax [65], and Oncorhynchus masou [66]. LC-PUFAs, as structural components of plasma membrane, enhanced the fluidity of cell membranes [67], results in the changes of other physical membrane properties and affects the membrane-protein binding [68]. These results could be used to explain the increase of LC-PUFAs level in euryhaline and migratory teleosts during salinity acclimation. While in alkalinity stress, the enriched pathways were mostly related to disease, pathogen defense and innate immune response, suggesting that tilapia was more vulnerable to alkaline stress than saline stress.

While in salinity-alkalinity stress group, most pathways were related to biosynthesis and metabolism, such as pantothenate and CoA biosynthesis, neomycin, kanamycin and gentamicin biosynthesis, nicotinate and nicotinamide metabolism and 2-Oxocarboxylic acid metabolism. In saline & alkaline, the enriched pathways of significant DEGs mostly related to antigen processing and antibody production. In saline & saline-alkaline, most pathways were related to secretion and metabolism, such as gastric acid secretion, pancreatic secretion, nitrogen metabolism and propanoate metabolism. While in alkaline & saline-alkaline, most pathways were related to immunity, and mismatch repair, nucleotide excision repair, DNA replication, glycine, serine and threonine metabolism, cell cycle, sphingolipid metabolism, glycosaminoglycan degradation, and so on, which were related to osmoregulation. Inokuchi et al showed that the apoptosis of freshwater-type mitochondrial rich cells would be enhanced and stimulated the recruitment of mitochondrial rich-cell complexes to develop hypo-osmoregulatory ability for seawater adaptation after transfer from hypo-osmotic to hyper-osmotic water in Oreochromis mossambicus [69]. Moreover, osmotic stress has been reported to cause oxidative stress, protein damage, DNA double-strand breaks and oxidative base modification [70,71,72]. These results are in consistent with the enrichment of canonical pathways such as mismatch repair, nucleotide excision repair, DNA replication, glycine, serine and threonine metabolism and cell cycle in response to osmotic tolerance of the hybrid tilapia. In saline & alkaline & salinity-alkalinity, 28 pathways were found, of which 10 pathways were related to metabolism, such as glutathione metabolism, glycine, serine and threonine metabolism, purine metabolism and pyrimidine metabolism.

Conclusions

Here, the transcriptome of the hybrid tilapia response to saline stress, alkaline stress and saline-alkaline stress was examined following the application of RNA-seq. Through more than 40 million sequence reads and the expression data of the four libraries, the study provided some useful insights into signal transduction pathways in hybrid tilapia and offered a number of candidate genes as potential markers of tolerance to osmotic stress for tilapia. The present research described for the first time the effects of salinity, alkalinity and salinity-alkalinity on gilled classic osmoregulation-related genes and signaling pathways during acclimation to saline and alkaline environment in tilapia. Our results provide a deep insight into the gill osmotic stress response of the hybrid tilapia to saline stress, alkaline stress and the mix of saline-alkaline stress that will be helpful in identifying reliable osmoregulation target genes and pathways for fish in osmotic stress.

Methods

Sample collection

The gill tissue samples were obtained from the hybrid tilapia (O. mossambicus♀ × O. urolepis hornorum♂) specimens at the Pearl River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou, China. All fishes were acclimated in the laboratory environment (28 °C) for two weeks before the experiment. The average weight of the experimental fishes was 51.27 ± 3.78 g. Three treatments with three replicates each were set up. The tilapias were divided into three groups (90 tilapias per group) and acclimated in saline stress (S, 25 ‰), alkali stress (A, 4 ‰) and saline-alkali stress (SA, S: 15 ‰, A: 4 ‰) at 28 °C, respectively. The control group was maintained in fresh water (saline, 0 ‰; alkali, 0 ‰). After 24 h (h) of treatment, 18 tilapias from each group were randomly selected and anesthetized by immersion in 100 mg/L MS-222 (Sigma, MO, USA) diluted in aquarium water. Following anesthesia, the fishes were decapitated and the gills were dissected immediately and immersed into liquid nitrogen for high-throughput sequencing. After 24 h of osmotic stress, the remaining fishes were released.

RNA isolation, library construction and illumina sequencing

The TRIzol® reagent (Invitrogen, USA) was used to extract the total RNA from four samples according to the manufacturer’s protocol and then treated by DNase I to implement the DNA digestion and obtained pure RNA products. Finally, the integrity of the total RNA was confirmed using 1% agarose gel electrophoresis and an Agilent 2100 Bio analyzer (Agilent Technologies, USA). Equal amounts of RNA from different individuals in the same group were pooled for library construction. Next generation sequencing library preparations were constructed according to the manufacturer’s protocol (NEBNext® Ultra™ RNA Library Prep Kit for Illumina®). Subsequently, the deep sequencing of the four mRNA libraries was performed using Illumina HiSeq 2000 according to the manufacturer’s instructions at the GENEWIZ (Genewiz, Suzhou, China).

Basic analysis of sequencing data

In order to remove technical sequences, including adapters, polymerase chain reaction (PCR) primers, or fragments thereof, and the quality of bases lower than 20, the pass filter data in the fastq format were processed by Trimmomatic (v0.30) to provide high-quality, clean data. Firstly, the reference genome sequences and the gene model annotation files of relative species were downloaded from genome website, such as NCBI. Secondly, Hisat2 (v2.0.1) was used to indexed with reference genome sequence. Finally, we mapped the clean reads to reference sequences of Oreochromis niloticus (National Center for Biotechnology Information, NCBI; accession numbers: GCA_001858045.3, https://www.ncbi.nlm.nih.gov/assembly/GCA_001858045.3) via software Hisat2 (v2.0.1). In the beginning transcripts in fasta format are converted from known gff annotation file and indexed properly. Then, with the file as a reference gene file, HTSeq (v0.6.1) estimated gene and isoform expression levels from the pair-end clean data.

Differential expression analysis and bioinformatics analysis of genes

DESeq2 and edegR method were used to conducted differential expression analysis. FPKM (Fragments Per Kilobases per Millionreads) method can eliminate the influence of gene length and sequencing amount on the calculated gene expression, and the calculated gene expression level can be directly used to compare the expression differences between different genes, so it was used to normalized the data. After adjusting by BH (fdr correction with Benjamini/Hochberg) for controlling the false discovery rate, the P-values < 0.05 and |log2FC| > = 1 were set to detect significant differentially expressed ones [73, 74]. The gene ontology (GO) enrichment analysis indicated that the significant differentially expressed genes with a significant P-value less than 0.05 were enriched for biological processes, molecular functioning, and cellular components through the software GO-TermFinder. KEGG (Kyoto Encyclopedia of Genes and Genomes) is a collection of databases dealing with genomes, biological pathways, diseases, drugs, and chemical substances (https://www.genome.jp/kegg/kegg3.html) [75]. We used scripts in house to enrich significant differential expression gene in KEGG pathways. The Venny online software (http://bioinfogp.cnb.csic.es/tools/venny/) was used to combine the analyses of differentially expressed genes. Gene Set Enrichment Analysis (GSEA) was also used to detect the expression changes of the whole gene set. This enrichment tool can fully detect genes that are not significantly different in expression levels but are plays important biological role at the overall level.

RNA-Seq data validation by qPCR and statistical analysis

To validate and measure the differential expression of mRNAs identified by high-throughput sequencing, 16 differentially expressed mRNAs were randomly selected for real-time quantitative PCR analyses. Real-time quantitative PCR was performed using the SYBR Green PCR Master Mix (Life Technologies, USA) according to the manufacture’s protocol in an ABI 7300 Real-Time PCR System (Applied Biosystems, Foster City, CA) using the following conditions: 50 °C for 2 min, 95 °C for 2 min, 40 cycles of 95 °C for 15 s, 60 °C for 15 s and 72 °C 30s. At the end of qPCR assay, the melting curve was analyzed to evaluate the amplification specificity. All primer sequences are listed in Table 11. The relative expression was calculated using the 2-Ct method and β-actin was used as the internal control to normalize mRNAs expression. The results are presented as the mean ± standard error. The analysis of all data was implemented using the statistical software SPSS 22.0 (SPSS, Chicago, IL, USA) and the expression levels of mRNAs were analyzed using one-way ANOVA method.

Table 11 Primers of different expression genes used for the qRT-PCR analysis

Availability of data and materials

All the data supporting our findings are contained within the manuscript. The clean sequence reads dataset supporting the conclusions of this article are available on NCBI BioProject with accession number PRJNA558948 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA558948) and will be made publicly available upon acceptance of this manuscript for publishing.

Abbreviations

A:

Alkaline stress group

bp:

Base pairs (measuring unit)

C:

Control group

DEGs:

Differentially expressed genes

FPKM:

Fragments Per Kilobases per Millionreads

GO:

Gene ontology

GSEA:

Gene Set Enrichment Analysis

KEGG:

Kyoto encyclopedia of genes and genomes

KO:

KEGG ortholog database

qRT-PCR:

Quantitative Real Time Polymerase Chain Reaction

RNA-Seq:

RNA-sequencing

S:

Saline stress group

SA:

Saline-Alkaline group

References

  1. 1.

    De Lisle PF, Roberts MH Jr. The effect of salinity on cadmium toxicity to the estuarine mysid Mysidopsis bahia: role of chemical speciation. Aquat Toxicol. 1998;12(4):357–70. https://doi.org/10.1016/0166-445X(88)90062-8.

  2. 2.

    Nakamura K, Iwaizumi K, Yamada S. Hemolymph patterns of free amino acids in the brine shrimp Artemia franciscana after three days starvation at different salinities. Comp Biochem Physiol A. 2007;147(1):254–9. https://doi.org/10.1016/j.cbpa.2007.01.008.

  3. 3.

    McFarland K, Ludovic D, Aswani KV. Effect of acute salinity changes on hemolymph osmolality and clearance rate of the nonnative mussel, Perna viridis, and the native oyster, Crassostrea virginica, in Southwest Florida. Aquat Invasions. 2013;8(3):299–310. https://doi.org/10.3391/ai.2013.8.3.06.

  4. 4.

    Chong-Robles J, Charmantier G, Boulo V, Lizárraga-Valdéz J, Enríquez-Paredes L, Giffard-Mena I. Osmoregulation pattern and salinity tolerance of the white shrimp Litopenaeus vannamei (Boone, 1931) during post-embryonic development. Aquaculture. 2014;422(423):261–7. https://doi.org/10.1016/j.aquaculture.2013.11.034.

  5. 5.

    Yang WK, Chung CH, Hui CC, Tang CH, Lee TH. Different expression patterns of renal Na(+)/K(+)-ATPaseα-isoform-like proteins between tilapia and milkfish following salinity challenges. Comp Biochem Physiol B. 2016;202:23–30. https://doi.org/10.1016/j.cbpb.2016.07.008.

  6. 6.

    Yao Z, Lai Q, Zhou K, Rizalita RE, Wang H. Developmental biology of medaka fish (Oryzias latipes) exposed to alkalinity stress. J Appl Ichthyol. 2010;26(3):397–402. https://doi.org/10.1111/j.1439-0426.2009.01360.x.

  7. 7.

    Saha N, Kharbuli ZY, Bhattacharjee A, Goswami C, Häussingerb D. Effect of alkalinity (pH 10) on ureogenesis in the air-breathing walking catfish, Clarias batrachus. Comp Biochem Physiol Part A. 2002;132(2):353–64. https://doi.org/10.1016/S1095-6433(02)00044-2.

  8. 8.

    Brennan RS, Galvez F, Whitehead A. Reciprocal osmotic challenges reveal mechanisms of divergence in phenotypic plasticity in the killifish Fundulus heteroclitus. Exp Biol. 2015;218(8):1212–22. https://doi.org/10.1242/jeb.110445.

  9. 9.

    Kozak GM, Brennan R, Berdan EL, Fuller RC, Whitehead A. Functional and population genomic divergence within and between two species of killifish adapted to different osmotic niches. Evolution. 2014;68(1):63–80. https://doi.org/10.1111/evo.12265.

  10. 10.

    Velotta JP, McCormick SD, O'Neill RJ, Schultz ET. Relaxed selection causes microevolution of seawater osmoregulation and gene expression in landlocked alewives. Oecologia. 2014;175(4):1081–92. https://doi.org/10.1007/s00442-014-2961-3.

  11. 11.

    Velotta JP, Mccormick SD, Schultz ET. Trade-offs in osmoregulation and parallel shifts in molecular function follow ecological transitions to freshwater in the alewife. Evolution. 2015;69(10):2676–88. https://doi.org/10.1111/evo.12774.

  12. 12.

    Whitehead A, Roach JL, Zhang S, Galvez F. Genomic mechanisms of evolved physiological plasticity in killifish distributed along an environmental salinity gradient. Proc Natl Acad Sci U S A. 2011;108(15):6193–8. https://doi.org/10.1073/pnas.1017542108.

  13. 13.

    Hu D, Pan L, Zhao Q, Ren Q. Transcriptomic response to low salinity stress in gills of the Pacific white shrimp Litopenaeus vannamei. Mar Genom. 2015;24:297–304. https://doi.org/10.1016/j.margen.2015.07.003.

  14. 14.

    Xia JH, Bai Z, Meng Z, Zhang Y, Wang L, Liu F, et al. Signatures of selection in tilapia revealed by whole genome re-sequencing. Sci Rep. 2015;5:14168. https://doi.org/10.1038/srep14168.

  15. 15.

    Wang M, Lu MX. Tilapia polyculture: a global review. Aquac Res. 2016;47(8):2363–74. https://doi.org/10.1111/are.12708.

  16. 16.

    Yuan Y, Yuan Y, Dai Y, Gong Y. Economic profitability of tilapia farming in China. Aquacult Int. 2017;25(3):1253–64. https://doi.org/10.1007/s10499-017-0111-8.

  17. 17.

    Hwang PP, Lee TH, Lin LY. Ion regulation in fish gills: recent progress in the cellular and molecular mechanisms. Am J Phys Regul Integr Comp Phys. 2011;301(4):28–47. https://doi.org/10.1152/ajpregu.00047.2011.

  18. 18.

    Takei Y, Hiroi J, Takahashi H, Sakamoto T. Diverse mechanisms for body fluid regulation in teleost fishes. Am J Phys Regul Integr Comp Phys. 2014;307(7):778–92. https://doi.org/10.1152/ajpregu.00104.2014.

  19. 19.

    Perry SF. The chloride cells: structure and function in the gills of freshwater fishes. Annu Rev Physiol. 1997;59:325–47. https://doi.org/10.1146/annurev.physiol.59.1.325.

  20. 20.

    Cnaani A, Hulata G. Improving salinity tolerance in tilapias: past experience and future prospects. Isr Aquacult Bamidgeh IIC J. 2011;63(1):1–21. https://doi.org/10.1071/MF11030.

  21. 21.

    Zhu HP, Liu ZG, Lu MX, Gao FY, Ke XL, Ma DM, et al. Screening and identification of a microsatellite marker associated with sex in Wami tilapia Oreochromis urolepis hornorum. J Genet. 2016;95(2):283–9. https://doi.org/10.1007/s12041-016-0653-y.

  22. 22.

    Claiborne J, Evans DH. Osmotic and ionic regulation in fishes. Osmotic and Ionic Regulation. 2008:295–366. https://doi.org/10.1201/9780849380525.ch8.

  23. 23.

    Whittamore JM. Osmoregulation and epithelial water transport: lessons from the intestine of marine teleost fish. J Comp Physiol B. 2012;182(1):1–39. https://doi.org/10.1007/s00360-011-0601-3.

  24. 24.

    Chourasia TK, D'Cotta H, Baroiller JF, Slosman T, Cnaani A. Effects of the acclimation to high salinity on intestinal ion and peptide transporters in two tilapia species that differ in their salinity tolerance. C Comp Biochem Physiol A. 2018;218:16–23. https://doi.org/10.1016/j.cbpa.2018.01.004.

  25. 25.

    Evans DH, Piermarini PM, Choe KP. The multifunctional fish gill: dominant site of gas exchange, osmoregulation, acid-base regulation, and excretion of nitrogenous waste. Physiol Rev. 2005;85(1):97–177. https://doi.org/10.1152/physrev.00050.2003.

  26. 26.

    Carmona R, García-Gallego M, Sanz A, Domezaín A, Ostos-Garrido MV. Chloride cells and pavement cells in gill epithelia of Acipenser naccarii: ultrastructural modifications in seawater-acclimated specimens. J Fish Biol. 2010;64(2):553–66. https://doi.org/10.1111/j.0022-1112.2004.00321.x.

  27. 27.

    Hirai N, Tagawa M, Kaneko T, Seikai T, Tanaka M. Distributional changes in Branchial chloride cells during freshwater adaptation in Japanese Sea bass Lateolabrax japonicus. Zool Sci. 1999;16(1):43–9. https://doi.org/10.2108/zsj.16.43.

  28. 28.

    Brown P. Gill chloride cell surface-area is greater in freshwater-adapted adult sea trout (Salmo trutta, L.) than those adapted to sea water. J Fish Biol. 1992;40:481–4. https://doi.org/10.1111/j.1095-8649.1992.tb02596.x.

  29. 29.

    Franklin CE. Surface ultrastructural changes in the gills of sockeye salmon (teleostei: Oncorhynchus nerka) during seawater transfer: comparison of successful and unsuccessful seawater adaptation. J Morphol. 2010;206(1):13–23. https://doi.org/10.1002/jmor.1052060103.

  30. 30.

    Kültz D, Jürss K, Jonas L. Cellular and epithelial adjustments to altered salinity in the gill and opercular epithelium of a cichlid fish (Oreochromis mossambicus). Cell Tissue Res. 1995;279(1):65–73. https://doi.org/10.1007/BF00300692.

  31. 31.

    Ke C, Erchao L, Tongyu L, Chang X, Xiaodan W, Heizhao L, et al. Transcriptome and molecular pathway analysis of the Hepatopancreas in the Pacific white shrimp Litopenaeus vannamei under chronic low-salinity stress. PLoS One. 2015;10(7):e0131503. https://doi.org/10.1371/journal.pone.0131503.

  32. 32.

    Li EC, Wang SL, Li C, Wang XD, Chen K, Chen LQ. Transcriptome sequencing revealed the genes and pathways involved in salinity stress of Chinese mitten crab Eriocheir sinensis. Physiol genomics. 2014;46(5):177–90. https://doi.org/10.1152/physiolgenomics.00191.2013.

  33. 33.

    Marie AD, Smith S, Green AJ, Rico C, Lejeusne C. Transcriptomic response to thermal and salinity stress in introduced and native sympatric Palaemon caridean shrimps. Sci Rep. 2017;7(1):13980. https://doi.org/10.1038/s41598-017-13631-6.

  34. 34.

    Lockwood BL, Somero GN. Transcriptomic responses to salinity stress in invasive and native blue mussels (genus Mytilus). Mol Ecol. 2011;20(3):517–29. https://doi.org/10.1111/j.1365-294X.2010.04973.x.

  35. 35.

    Guo B, Tang Z, Wu C, Xu K, Qi P. Transcriptomic analysis reveal an efficient osmoregulatory system in siberian sturgeon Acipenser baeri in response to salinity stress. Sci Rep. 2018;8(1):14353. https://doi.org/10.1038/s41598-018-32771-x.

  36. 36.

    Lv JJ, Liu P, Wang Y, Gao BQ, Chen P, Li J. Transcriptome analysis of Portunus trituberculatus in response to salinity stress provides insights into the molecular basis of osmoregulation. PLoS One. 2013;8(12):e82155. https://doi.org/10.1371/journal.pone.0082155.

  37. 37.

    Sahraeian SME, Mohiyuddin M, Sebra, R., Tilgner, H., Afshar PT, Au KF, et al. Gaining comprehensive biological insight into the transcriptome by performing a broad-spectrum RNA-seq analysis. Nat Commun,2017, 8(1). https://doi.org/10.1038/s41467-017-00050-4.

  38. 38.

    Edwards SL, Marshall WS. Principles and patterns of osmoregulation and Euryhalinities in fishes. Fish Physiology. 2012;32:1–44. https://doi.org/10.1016/B978-0-12-396951-4.00001-3.

  39. 39.

    Ouattara N'G, Bodinier C, Nègre-Sadargues G, D’Cotta H, Messad S, Charmantier G, et al. Changes in gill ionocyte morphology and function following transfer from fresh to hypersaline waters in the tilapia Sarotherodon melanotheron. Aquaculture. 2009;290(1–2):155–64. https://doi.org/10.1016/j.aquaculture.2009.01.025.

  40. 40.

    Nie HT, Jiang LW, Chen P, Huo ZM, Yang F, Yan XW. High throughput sequencing of RNA transcriptomes in Ruditapes philippinarum identifies genes involved in osmotic stress response. Sci Rep. 2017;7(1):4953. https://doi.org/10.1038/s41598-017-05397-8.

  41. 41.

    Dong WQ, Chen YX, Lu WX, Wu B, Qi PZ. Transcriptome analysis of Mytilus coruscus hemocytes in response to Vibrio alginnolyficus infection. Fish and Shellfish Immunol. 2017;70:560–7. https://doi.org/10.1016/j.fsi.2017.08.034.

  42. 42.

    Lam SH, Lui EY, Li Z, Cai S, Sung WK, Mathavan S, et al. Differential transcriptomic analyses revealed genes and signaling pathways involved in iono-osmoregulation and cellular remodeling in the gills of euryhaline Mozambique tilapia Oreochromis mossambicus. BMC Genomics. 2014;15:921. https://doi.org/10.1186/1471-2164-15-921.

  43. 43.

    Wong KS, Ozaki H, Suzuki Y, Iwasaki W, Takei Y. Discovery of osmotic sensitive transcription factors in fish intestine via a transcriptomic approach. BMC Genomics. 2014;15(1):1134. https://doi.org/10.1186/1471-2164-15-1134.

  44. 44.

    Zhao X, Yu H, Kong L, Li Q. Transcriptomic responses to salinity stress in the Pacific oyster Crassostrea gigas. PLoS One. 2012;7(9):e46244. https://doi.org/10.1371/journal.pone.0046244.

  45. 45.

    Zhu HP, Liu ZG, Gao FY, Lu MX, Liu YJ, Su HH, et al. Characterization and expression of Na+ /K+ −ATPase in gills and kidneys of the teleost fish Oreochromis mossambicus, Oreochromis urolepis hornorum and their hybrids in response to salinity challenge. Comp Biochem Physiol A. 2018;224:1–10. https://doi.org/10.1016/j.cbpa.2018.05.017.

  46. 46.

    Yamaguchi Y, Breves JP, Haws MC, Lerner DT, Grau EG, Seale AP. Acute salinity tolerance and the control of two prolactins and their receptors in the Nile tilapia (Oreochromis niloticus) and Mozambique tilapia (O. mossambicus): a comparative study. Gen Comp Endocr. 2018;257:168–76. https://doi.org/10.1016/j.ygcen.2017.06.018.

  47. 47.

    Yuan MZ, Jia QQ, Wang T, Lu Q, Tang LL, Wang YJ, et al. Dynamic responses of prolactin, growth hormone and their receptors to hyposmotic acclimation in the olive flounder Paralichthys olivaceus. Gen Comp Endocr. 2017;254:8–13. https://doi.org/10.1016/j.ygcen.2017.09.005.

  48. 48.

    Seale AP, Watanabe S, Breves JP, Lerner DT, Kaneko T, Grau EG. Differential regulation of trpv4 mRNA levels by acclimation salinity and extracellular osmolality in euryhaline tilapia. Gen Comp Endocr. 2012;178(1):123–30. https://doi.org/10.1016/j.ygcen.2012.04.020.

  49. 49.

    Ali MY, Pavasovic A, Mather PB, Prentis PJ. Analysis, characterisation and expression of gill-expressed carbonic anhydrase genes in the freshwater crayfish Cherax quadricarinatus. Gene. 2015;564(2):176–87. https://doi.org/10.1016/j.gene.2015.03.074.

  50. 50.

    Watson CJ, Nordi WM, Esbaugh AJ. Osmoregulation and branchial plasticity after acute freshwater transfer in red drum, sciaenops ocellatus. Comp Biochem Physiol A. 2014;178:82–9. https://doi.org/10.1016/j.cbpa.2014.08.008.

  51. 51.

    Nesrine BB, Viviane B, Mireille CD, Klaus A, Guy C, Catherine LN. Osmoregulation in larvae and juveniles of two recently separated Macrobrachium species: expression patterns of ion transporter genes. Comp Biochem Physiol A. 2016;195:39–45. https://doi.org/10.1016/j.cbpa.2016.02.005.

  52. 52.

    Hao HH, Li YL, Yung CT, Jiun LH, Pung-Pung H. A new model for fish ion regulation: identification of ionocytes in freshwater- and seawater-acclimated medaka (Oryzias latipes). Cell Tissue Res. 2014;357(1):225–43. https://doi.org/10.1007/s00441-014-1883-z.

  53. 53.

    Tatsuya S, Satoshi O, Yudai N, Chiaki A, Hideya T, Taro W, et al. Osmotic/ionic status of body fluids in the euryhaline cephalopod suggest possible parallel evolution of osmoregulation. Sci Rep. 2015;5:14469. https://doi.org/10.1038/srep14469.

  54. 54.

    Jiang W, Tian X, Fang Z, Li L, Dong S, Li H, et al. Metabolic responses in the gills of tongue sole (Cynoglossus semilaevis) exposed to salinity stress using NMR-based metabolomics. Sci Total Environ. 2019;653:465–74. https://doi.org/10.1016/j.scitotenv.2018.10.404.

  55. 55.

    Imai M, Seldin DW, Kokko JP. Effect of perfusion rate on the fluxes of water, sodium, chloride and urea across the proximal convoluted tubule. Kidney Int. 1977;11(1):18–27. https://doi.org/10.1038/ki.1977.3.

  56. 56.

    Alexander RT, Dimke H, Cordat E. Proximal tubular NHEs: sodium, protons and calcium? Renal Physiology. 2013;305(3):F229–F36. https://doi.org/10.1152/ajprenal.00065.2013.

  57. 57.

    Masashi S, Maria HV, Hideomi Y, Shoko H, Yuehong L, Takashi S, et al. Functional analysis of a novel missense NBC1 mutation and of other mutations causing proximal renal tubular acidosis. Pflug Arch Eur J Phy. 2008;455(4):583–93. https://doi.org/10.1007/s00424-007-0319-y.

  58. 58.

    Lo YF, Yang SS, Seki G, Yamada H, Horita S, Yamazaki O, et al. Severe metabolic acidosis causes early lethality in NBC1 W516X knock-in mice as a model of human isolated proximal renal tubular acidosis. Kidney Int. 2011;79(7):730–41. https://doi.org/10.1038/ki.2010.523.

  59. 59.

    Demirci FY, Chang MH, Mah TS, Romero MF, Gorin MB. Proximal renal tubular acidosis and ocular pathology: a novel missense mutation in the gene (SLC4A4) for sodium bicarbonate cotransporter protein (NBCe1). Mol Vis. 2006;12(35–36):324–30.

  60. 60.

    Dinour D, Chang MH, Satoh J, Smith BL, Angle N, Knecht A, et al. A novel missense mutation in the sodium bicarbonate cotransporter (NBCe1/SLC4A4) causes proximal tubular acidosis and glaucoma through ion transport defects. J Biol Chem. 2004;279(50):52238–46. https://doi.org/10.1074/jbc.%20M406591200.

  61. 61.

    Inatomi J, Horita S, Braverman N, Sekine T, Yamada H, Suzuki Y, et al. Mutational and functional analysis of SLC4A4 in a patient with proximal renal tubular acidosis. Pflug Arch Eur J Phy. 2004;448(4):438–44. https://doi.org/10.1007/s00424-004-1278-1.

  62. 62.

    Cooper GJ, Zhou Y, Bouyer P, Grichtchenko II, Boron WF. Transport of volatile solutes through AQP1. J Physiol. 2010;542(1):17–29. https://doi.org/10.1113/jphysiol.2002.023218.

  63. 63.

    Lu X, Kong J, Luan S, Dai P, Meng X, Cao B, et al. Transcriptome analysis of the hepatopancreas in the pacific white shrimp (Litopenaeus vannamei) under acute ammonia stress. PLoS One. 2016;11(10):e0164396. https://doi.org/10.1371/journal.pone.0164396.

  64. 64.

    Haliloglu HI, Bayir A, Sirkecioglu AN, Aras NM, Atamanalp M. Comparison of fatty acid composition in some tissues of rainbow trout (Oncorhynchus mykiss) living in seawater and freshwater. Food Chem. 2004;86(1):55–9. https://doi.org/10.1016/j.foodchem.2003.08.028.

  65. 65.

    Cordier M, Brichon G, Weber JM, Zwingelstein G. Changes in the fatty acid composition of phospholipids in tissues of farmed sea bass (Dicentrarchus labrax) during an annual cycle. Roles of environmental temperature and salinity. Comp Biochem Physiol B. 2002;133(3):281–8. https://doi.org/10.1016/S1096-4959(02)00149-5.

  66. 66.

    Li HO, Yamada J. Changes of the fatty acid composition in smolts of masu salmon (Oncorhynchus masou), associated with desmoltification and sea-water transfer. Comp Biochem Physiol A. 1992;103:221–6. https://doi.org/10.1016/0300-9629(92)90266-S.

  67. 67.

    Bell MV, Henderson RJ, Sargent JR. The role of polyunsaturated fatty acids in fish. Comp Biochem Physiol B. 1986;83(4):711–9. https://doi.org/10.1016/0305-0491(86)90135-5.

  68. 68.

    Hulbert AJ, Else PL. Membranes as possible pacemakers of metabolism. J Theor Biol. 1999;199(3):257–74. https://doi.org/10.1006/jtbi.1999.0955.

  69. 69.

    Inokuchi M, Kaneko T. Recruitment and degeneration of mitochondrion-rich cells in the gills of Mozambique tilapia Oreochromis mossambicus during adaptation to a hyperosmotic environment. Comp Biochem Physiol A. 2012;162(3):245–51. https://doi.org/10.1016/j.cbpa.2012.03.018.

  70. 70.

    Kültz D, Chakravarty D. Hyperosmolality in the form of elevated NaCl but not urea causes DNA damage in murine kidney cells. Proc Natl Acad Sci U S A. 2001;98(4):1999–2004. https://doi.org/10.1073/pnas.98.4.1999.

  71. 71.

    Zhang Z, Dmitrieva NI, Park JH, Levine RL, Burg MB. High urea and NaCl carbonylate proteins in renal cells in culture and in vivo, and high urea causes 8-oxoguanine lesions in their DNA. Proc Natl Acad Sci U S A. 2004;101(25):9491–6. https://doi.org/10.1073/pnas.0402961101.

  72. 72.

    Kültz D. Hyperosmolality triggers oxidative damage in kidney cells. Proc Natl Acad Sci U S A. 2004;101(25):9177–8. https://doi.org/10.1073/pnas.0403241101.

  73. 73.

    Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8. https://doi.org/10.1038/nmeth.1226.

  74. 74.

    Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7(10):986–95. https://doi.org/10.1101/gr.7.10.986.

  75. 75.

    Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(40):D480–4. https://doi.org/10.1093/nar/gkm882.

Download references

Acknowledgements

We would like to thank the reviewers for their helpful comments on the original manuscript.

Funding

This work was supported by grants from the Science and Technology Program of Guangzhou, China (Grant No. 201803020034), the Central Public-interest Scientific Institution Basal Research Fund, CAFS (Grant No. 2019KQ-003), the Guangdong Provincial Special Fund For Modern Agriculture Industry Technology Innovation Teams (Grant No. 2019KJ150), the National Natural Science Foundation of China (Grant No. 31001108) and the National Freshwater Genetic Resource Center (Grant No. 2019DKA30470). There was no additional external funding received for this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

HS and HZ designed the experiment and drafted the manuscript, HS performed the experiments and analyzed the data. DM, ZL and FG participated in experimental design and coordination and contributed to the manuscript preparation. All authors reviewed and approved the manuscript.

Correspondence to Huaping Zhu.

Ethics declarations

Ethics approval and consent to participate

The hybrid tilapia used in this study was produced by at the Pearl River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou, China. All experiments were performed according to the Guide for the Care and Use of Laboratory Animals of China. This study was also approved by the Animal Care and Use committee of the Centre for Applied Aquatic Genomics at the Chinese Academy of Fishery Sciences.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Su, H., Ma, D., Zhu, H. et al. Transcriptomic response to three osmotic stresses in gills of hybrid tilapia (Oreochromis mossambicus female × O. urolepis hornorum male). BMC Genomics 21, 110 (2020). https://doi.org/10.1186/s12864-020-6512-5

Download citation

Keywords

  • Tilapia
  • Transcriptome
  • Osmoregulation
  • Osmotic stress