Transcriptomic analysis reveals putative osmoregulation mechanisms in the kidney of euryhaline turbot Scophthalmus maximus responded to hypo-saline seawater*

CUI Wenxiao MA Aijun HUANG Zhihui WANG Xin’an SUN Zhibin LIU Zhifeng ZHANG Wei YANG Jingkun ZHANG Jinsheng QU Jiangbo

1 College of Fisheries and Life Science, Shanghai Ocean University, Ministry of Education, Shanghai 201306, China

2 Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, China

3 Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences; Shandong Key Laboratory of Marine Fisheries Biotechnology and Genetic Breeding; Qingdao Key Laboratory for Marine Fish Breeding and Biotechnology, Qingdao 266071,China

4 Yantai Tianyuan Aquatic Limited Corporation, Yantai 264003, China

Abstract Turbot harbor a relatively remarkable ability to adapt to opposing osmotic challenges and are an excellent model species to study the physiological adaptations of f lounder associated with osmoregulatory plasticity. The kidney transcriptome of turbot treated 24 h in water of hypo-salinity (salinity 5) and seawater(salinity 30) was sequenced and characterized. In silico analysis indicated that all unigenes had signif icant hits in seven databases. The functional annotation analysis of the transcriptome showed that the immune system and biological processes associated with digestion, absorption, and metabolism played an important role in the osmoregulation of turbot in response to hypo-salinity. Analysis of biological processes associated with inorganic channels and transporters indicated that mineral absorption and bile secretion contributed to iono-osmoregulation resulting in cell volume regulation and cell phenotypic plasticity. Moreover, we analyzed and predicted the mechanisms of canonical signaling transduction. Biological processes involved in renin secretion, ECM-receptor interaction, adherens junction, and focal adhesion played an important role in the plasticity phenotype in hypo-stress, while the signal transduction network composed of the MAPK signaling pathway and PI3K-Akt signaling pathway with GABAergic synapse, worked in hypoosmoregulation signal transduction in the turbot. In addition, analysis of the tissue specif icity of targeted gene expression using qPCR during salinity stress was carried out. The results showed that the kidney, gill,and spleen were vital regulating organs of osmotic pressure, and the osmoregulation pattern of euryhaline f ish diff ered among species.

K eyword: osmoregulation mechanisms; Scophthalmus maximus; hypo-saline; signal transduction; cell volume regulation; gene expression

1 INTRODUCTION

As an inherent physicochemical property of water,salinity is one of the main factors limiting species distribution (Brennan et al., 2015). Euryhaline f ish possess particular osmotic regulation mechanisms that enable them to sense and compensate for salinity f luctuation in the marine environment (Edwards and Marshall, 2012). They can perceive osmosensory signaling that activate control system, then coordinate eff ectors to counteract osmotic stress (Fiol and Kültz,2007). The adaptability of euryhaline f ish to the salinity in a dramatically changing aquatic environment is essential for survival, and provides a classical model for investigating osmotic regulation(Lam et al., 2014). The osmoregulatory mechanisms that euryhaline f ish harbor include complex functional networks encoded by widely distributed functional proteins (Kültz, 2015). Using transcriptomics approaches, it was shown that the functional distributions of diff erentially expressed genes (DEGs)found in the kidney is wide than in the gill of the European eel (Anguillaanguilla) (Kalujnaia et al.,2007).

Turbot (Scophthalmusmaximus) is a euryhaline f lounder with great commercial value that is native to the North Atlantic and its saltwater-freshwater border(Pereiro et al., 2012). Turbot have a relatively remarkable ability to adapt to opposing osmotic challenges and have been proven an excellent model species to study physiological adaptations of f lounder associated with osmoregulatory plasticity. There have been many studies on the adaptability of turbot to salinity mainly in the aspects of physiology (Imsland et al., 2001; Tang et al., 2006; Dietz et al., 2013), but the osmoregulatory mechanism at the transcriptional level is still unknown. Osmoregulation in teleosts is mainly mediated by the kidneys, gills, and intestine(Marshall and Grosell, 2005). Turbot kidney is an important organ of osmotic regulation, which plays an irreplaceable role in osmotic regulation (Prodocimo et al., 2008). Kidney filtration and reabsorption is important for osmoregulation as any excess uptake of divalent ions is excreted by the kidney (Beyenbach,1986, 2004). Moreover, some osmoregulatory hormones in teleost fishes, such as renin and cortisol,are secreted by the kidney and are critical for body fluid regulation through its eff ect on tubular water reabsorption in the kidney (Babey et al., 2011). In addition, there are a large number of special osmoregulatory eff ectors, osmotic stress signaling and ion channels acting on osmoregulation in the kidney (Kasprowicz, 2011).

In this study, we employed next-generation sequencing (NGS) technology to sequence and assemble the de novo kidney transcriptome of euryhaline turbot in response to a hypo-osmotic environment. First, the transcripts were assembled to unigenes the Trinity transcriptome assembler.Subsequently, GO and KEGG enrichment analyses revealed unigenes that are involved in diff erent biological process and signaling pathways. The aims of this study were to identify transduction processes of osmotic pressure signaling and eff ector mechanism response to low salinity stress. Although variations in transcript abundance does not mean the variations of protein contents, nevertheless, this study identif ied many crucial biological processes and signaling pathways involved in osmoregulation from the perspective of transcriptome in kidneys under low salinity stress. To the best of our knowledge, this was the f irst NGS study to provide a transcriptome-wide perspective on the mechanisms of osmoregulation in turbot, especially the process of signal transduction in kidneys of euryhaline f ish in response to hyposalinity.

2 MATERIAL AND METHOD

2.1 Material

2.1.1 Ethics statement

The handing of turbot was performed in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The study procedures were approved by the Experimental Animal Ethics Committee in Chinese Academy of Fishery Sciences, China.

2.1.2 Fish

The healthy juvenile turbot (12.52±0.34 cm) were from Yantai Tianyuan Aquaculture Co., Ltd. in Shandong Province, China. Healthy f ish were maintained for 2 weeks in 150-L tanks. Water was aerated to keep normoxic oxygen saturation at all times and the water was not exchanged. Diluted seawater (salinity of 5, 10) was obtained by mixing freshwater with seawater (salinity of 30) taken from a deep well, whereas hypersaline water (salinity of 50)was obtained by mixing seawater crystals with seawater. We sampled intestine, kidney, gill, liver,skin, spleen, and brain tissues from turbot that had been treated in their respective diluted seawater or seawater, and hypersaline environments for 24 h.Before sampling, f ish individuals were anesthetized with MS-222. Three biological replicates were set in separate aquarium with same aquatic condition in each experiment. The fresh tissues were immersed in RNAstor immediately at 4°C overnight and then frozen at -80°C until the next work.

2.2 Method

2.2.1 Transcriptome sequencing

Total RNA was extracted from kidney tissue(salinity of 5 and 30) samples (n=3) using Trizol reagent (Invitrogen, Carlsbad, CA, USA) and RNA quality and quantity were evaluated in a Qubit and a Nanodrop, respectively. After checking RNA purity and concentration, the integrity of RNA samples was assessed using an Agilent 2100 instrument (Agilent,Santa Clara, CA, USA). RNAs from individual f ish in each group were used to construct cDNA libraries.The six cDNA libraries were constructed and subjected to pair-end sequencing on an Illumina Hiseq 2000 platform at Novogene Bioinformatics Technology Co. Ltd. (Beijing, China).

2.2.2 Quality control, assembly, and gene annotation

Raw reads were produced through base calling and stored in fastq format. The raw data were cleaned after f iltration by removing the adaptors, poly-N of more than 10%, and low-quality sequences (base quality score ＜Q20). Transcriptome assembly was performed using Trinity software (Grabherr et al.,2011) as described for de novo transcriptome assembly without the reference genome. The k-mer value was set to 25 and min_kmer_cov set to 2, by default. The longest assembled sequences of each transcription set were def ined as unigenes. Genes encoding protein domains were identif ied by searching against non-redundant (NR), NT, Swiss-Prot, PFAM databases through the BLASTp. GO functions for all unigenes were classif ied using Gene2go v2.5 (Götz et al., 2008) with an E-value threshold of 1e-6. The pathway annotations and classif ications of unigenes were performed based on KEGG and COG databases using the BLAST, with an E-value threshold of 1e-10 and 1e-3, respectively.

Trinity was used to assemble clean reads into the sequences by RSEM (Li and Dewey, 2011). The expression of genes was normalized and calculated by the FPKM (Fragments Per kb per Million fragments) method (Trapnell et al., 2010). Diff erential expression of unigenes was analyzed using the MAplot-based method with random sampling model(MARS) in DEGseq (Anders and Huber, 2010). ThePvalue was adjusted by means of theqvalue. A padj＜0.05 was set as the threshold for signif icantly diff erential expression, according to biological repetition. Functional categories analysis on Go database of the DEGs was performed using GOseq(Young et al., 2010) based on a hypergeometric distribution test. GO terms with false discovery rates(FDR) ≤0.05 were def ined as signif icantly enriched in DEGs. Using KOBAS (2.0) software, pathway enrichment analysis of the DEGs based on KEGG was performed (Mao et al., 2005).

2.2.3 Fluorescence quantitative PCR and statistical analysis

Total RNA was extracted from tissues treated with four gradients of salinity using an RNA prep pure Tissue Kit (Tiangen, Beijing, China). First strand cDNA was synthesized by reverse transcription using Quant Reverse Transcriptase. All primers for qPCR were designed using Primer5.0 software. Primers are listed in Supplementary Table S1. Real-time PCR was performed on a 185-5096 CFX96 Real-Time PCR Detection System (Beijing yue da Biotechnology Co.Ltd., Beijing, China) in a 15-μL reaction using 1 μL cDNA, 10 μL 2×TransStart Top Green qPCR SuperMix, 0.4 μL forward and reverse primers, 7.8 μL RNase-Free ddH2O and 0.4 μL 50×Passive Reference Dye 1. Cycling conditions were 94°C for 30 s followed by 40 cycles of 95°C for 5 s, and 60°C for 30 s. The blastN and blastX tools were used for a homogeneity search available online at http://blast.ncbi.nlm.nih.gov/Blast. All quantitative data shown are the means± S.E.M. of multiple samples (n=3). The expression of genes was analyzed by the 2-ΔΔCtmethod(Schmittgen and Zakrajsek, 2000). Statistics software SPSS19.0 (SPSS, Chicago, IL, USA) was used to test signif icance. The signif icance threshold was set atP＜0.05 for all tests.

3 RESULT AND DISCUSSION

3.1 Sequencing and de novo assembly of kidney transcriptomes

The results showed that in total an average of 42 652 674 clean reads were obtained from the six libraries of turbot samples under salinity stress.193 336 unigenes were obtained from 226 176 transcripts assembled with clean reads (Fig.1). As shown, the length distribution of all unigenes ranged from 201 to 25 839 bp, with an average length of 786 bp.

3.2 Functional annotation and classif ication of the transcripts

In order to determine their function, all unigenes were blast analyzed against publicly available databases. All unigenes had signif icant hits in the NCBI NR, NT, PFAM, Swiss Prot, KOG, GO, and KEGG databases (Table 1). Overall, 59 007 unigenes(30.52%) were annotated in at least one database, and 47 375 (24.50%), and 19 098 (9.88%) unigenes were annotated in the NT and KOG databases, respectively.According to the results annotated in the GO database,34 997 unigenes were divided into 54 subcategories(Fig.2). A total of 19 098 unigenes were annotated into 26 diff erent processes on the KOG database, and the largest category was represented by “signal transduction mechanisms”, which was followed by“general function prediction only”, and“posttranslational modif ication, protein turnover,chaperones” (Fig.3). The total annotated unigenes were classif ied into f ive major functional categories:cellular processes, environmental information processing, genetic information processing,metabolism, and organismal systems, which were further divided into 4, 3, 4, 12, and 9 subcategories,respectively (Fig.4), and the largest category was represented by “signal transduction”, which was followed by “endocrine system”.

Fig.1 Summary statistics of sequencing output and de novo assembly of kidney transcriptome

Fig.3 Categories of KOG

Table 1 Summary of unigene annotations

3.3 Assessment of the reliability of transcriptional data

To assess the DEGs obtained by RNA sequencing for the reliability of transcriptional data, qPCR was carried out with the total RNA used for RNA sequencing, to examine the expression of 10 DEGs in kidneys under salinities of 5 and 30, respectively,which included six upregulated unigenes and four downregulated unigenes (Fig.5a). Among the 10 selected genes, the relative expression of seven genes between brackish water- and seawater-treated kidneys matched with transcriptomes data, except genesAPOM,AATP6V1E1andH4. Our results showed that expression levels of the ten genes were obviously changed after salinity stress. This suggested that the real-time PCR data for hypo-salinity treated kidneys were in good concordance with the transcriptome data (Fig.5b).

3.4 Functional annotation analysis of transcriptomes

Analysis of the biological functions of DEGs showed that 23 signif icant enrichment pathways were upregulated and 16 of these were in the immunological category (Supplementary Table S2). There were 41 signif icant enrichment pathways in which genes were downregulated, other than legionellosis; these belonged to energetic metabolism related to digestion and absorption (Supplementary Table S3). We presumed that the immune system functions in hypoosmotic stress to defend turbot against changes of diversif ied external environmental factors involving in environmental salinity, temperature and autochthonous microorganisms (Narnaware et al.,2000). However, the mRNA expressions of genes in energetic metabolism pathways were reduced, which suggest SW acclimation was associated with decreased respiratory activity and energy production in kidney (Kalujnaia et al., 2007).

Fig.4 Categories of KEGG

Fig.5 Expression of 10 genes in kidney

In hypo-osmotic brackish water, turbot need to reduce salt loss and water gain with kidney (Farrell,2011). To maintain normal physiological functions,the signaling pathways related to cell morphology and cell volume regulation are upregulated (Lam et al.,2014). Moreover, the signaling pathways associated with maintaining cell morphology, by increasing ionocytes with the capacity to reabsorb dissolved salt lost to the f iltrate, were upregulated in hypo-osmotic stress. Consequently, signaling pathways and biological processes that were enriched in the kidney involved apoptosis, the cell cycle, oocyte meiosis, the p53 signaling pathway, regulation of actin cytoskeleton,adherens junction, focal adhesion, gap junction,signaling pathways regulating pluripotency of stem cells, and tight junctions (Supplementary Table S4).These signaling pathways and biological processes could play important roles in cell growth and death,movement, and cellular interactions in response to osmotic stress, which has been documented in other cell types (Kalujnaia et al., 2007; Hoff mann et al.,2009; Lam et al., 2014). In agreement with previous studies (Kültz and Burg, 1998; Kültz and Avila, 2001),our analysis indicated that transcripts associated with MAPK signaling, an important pathway in the osmotic stress response, were abundant in the hypo-salinity treated kidney of turbot.

Fig.6 Interaction of nine signaling pathways

3.5 Analysis of signal transduction in hypo-osmotic regulation

A total of 45 signal pathways were detected in 246 physiological processes of turbot hypoosmoregulation, of which 25 signal pathways belonged to the category of signal transduction(Supplementary Table S5). Diff erential expression analysis showed that the expression prof iles of the genes encoding core protein, as well as its upstream ligands and downstream acting protein genes, of nine out of 45 signaling pathways, changed. The results revealed that the ErbB, FoxO, Jak-STAT, MAPK,Notch, PI3K-Akt, phosphatidylinositol, Ras, and cGMP-PKG signaling pathway, which belong to the category of signal transduction, may play an important role in signal transduction of ionic and osmotic regulation, and the control of extracellular f luid volume in the kidney of turbot. Among these signal pathways, the ErbB, FoxO, Jak-STAT, MAPK, Notch,Ras, and cGMP-PKG signaling pathway belong to signal transduction mediated by enzyme-coupled receptors, while the PI3K-Akt signaling pathway and phosphatidylinositol signaling pathways are mediated by G-protein coupled receptors (Bonni et al., 1997;Vojtek and Der, 1998; Grant et al., 2002; Arden, 2004;Ehebauer et al., 2006). However, signal transduction using an enzyme- coupled receptor and G-protein coupled receptor require that the extracellular ligand(signal molecule) binds to the cell receptor to form a complex that acts on the target protein (Bakir and Sezerman, 2006).

Interplay between signaling pathways (Fig.6)indicated that the MAPK signaling pathway participated in the biological function of seven signaling pathways other than the phosphatidylinositol signaling pathway, and the PI3K-Akt signaling pathway interacted with f ive signaling pathways except for the phosphatidylinositol signaling pathway,Notch signaling pathway, and MAPK signaling pathway. Therefore, the importance of the MAPK signaling pathway (Fouchs et al., 2010) and PI3K-Akt signaling pathway (Xu et al., 2015) in the osmotic adjustment of turbot in a hypo-osmotic environment was also demonstrated, and also showed that two signaling pathways may integrate and precisely control multiple signals in cells to form a signaling network system and f inally control the appropriate response.

Analysis of the ligands of biological processes and signaling pathways enriched in the endocrine system suggested that the physiological productions of biological processes were involved in insulin secretion,melanogenesis, renin secretion, renin-angiotensin system, thyroid hormone synthesis, and signaling molecules of signaling pathways involving glucagon signaling pathway, GnRH signaling pathway, insulin signaling pathway, oxytocin signaling pathway, PPAR signaling pathway, prolactin signaling pathway, and the thyroid hormone signaling pathway, and all were through proteins. In addition, the results of DEGs showed that only the mRNA of genes encoding core proteins of renin secretion had changed. The analysis showed that renin was likely to be generated in hypoosmoregulation. Our f indings corroborated those of a previous report on the renin-angiotensin-aldosterone system, which occupy a crucial position in maintaining of homeostasis electrolyte and f luid, and in enhancing reabsorption of Na+from kidney (Kurtz A, Wagner,1999). Signal transduction in the organism requires interactions between the signal molecule and the receptor (Bhattacharya, 2010). In the signal transduction process of hypo-osmoregulation in turbot, biological processes that were enriched in signaling molecules and interactions involved cell adhesion molecules (CAMs), cytokine-cytokine receptor interactions, extracellular matrix (ECM)-receptor interactions, and neuroactive ligand-receptor interactions. The analysis of the amount and abundance of genes expressed in physiological processes suggested the importance of ECM-receptor interaction,cytokine-cytokine receptor interaction, neuroactive ligand-receptor interaction, and CAMs in signal transduction decreased in hypo-osmoregulation.

Cell signal transduction is a key process of cellto-cell communication essential for the coordination of cell activities, the control of cell growth and division, generation of tissues, and formation of morphology in multicellular organisms (Rappolee and Armant, 2009). Cells can communicate with each other not only by secreting chemical signals but also by cell contact and gap junction signaling(Kurt et al., 2011). The analysis of the enrichment of biological processes indicated that tight junction,adherens junction, focal adhesion, gap junction,signaling pathways regulating pluripotency of stem cells, were enriched in the cellular community. In addition, according to the mode of action, all played a role in cell-cell junctions, while focal adhesion also acted on the connection of cells to the ECM;however, all of the linkages acted on transmembrane transport and signal transduction. Analyses of the DEGs and enrichment of biological processes revealed that the mRNA of genes encoding the core proteins of adherens junctions and focal adhesion,important biological processes in the hypo-osmotic stress response, had changed. Therefore, the gap junctions that exist widely in animal tissues (Shibata et al., 2001) may play an important role in hypoosmoregulation signal transduction in turbot.

The conduction of nerve impulses is completed through the synapses in neurons-neurons or neuronseff ector cells (such as muscle cells). It was found that signaling pathways involving cholinergic synapse,dopaminergic synapse, GABAergic synapse,glutamatergic synapse, and serotonergic synapse were enriched in the nervous system. More importantly, all synapses are chemical. However,chemical synapses are biological junctions through which the signals of neurons can be exchanged with each other and with non-neuronal cells such as those in muscles or glands, whereas electrical synapses directly transfer electrical signals from one cell to another through gap junctions. Therefore, gap junctions associated with electrical synaptic signaling may not play a signif icant role in hypo-osmoregulation signal transduction in turbot. Analysis of the genes of chemical post-synaptic cell receptors and receptorrelated proteins showed diff erential expressions of genes associated with neurotransmitter receptors and receptor-related proteins of the GABAergic synapse,among the f ive chemical synapses. The diff erential expression underscored the importance of GABAergic synapses in hypo-osmoregulation signal transduction in turbot.

In summary, on the transcriptional level, biological processes involved in renin secretion, ECM-receptor interaction, adherens junction, and focal adhesion play an important role in the plasticity phenotype in hypo-stress, while the signal transduction network composed of the MAPK signaling pathway and PI3KAkt signaling pathway with the GABAergic synapse work in hypo-osmoregulation signal transduction in turbot.

3.6 Analysis of biological processes associated with inorganic channels and transporters

There are some studies on osmoregulation mechanisms of euryhaline f ish kidney, and basic schemas involving in ion channels and transporters for salinity stress responses have been proposed(Kültz, 2015). By analyzing the roles of signal pathways and biological processes that were signif icantly enriched in genes, biological processes involving mineral absorption (Fig.7) and bile secretion (Fig.8) associated with iono-osmoregulation were identif ied. Several homologs of key ion channels and transporters that aid in absorbing salt, such asSLC5A1(SGLT1),SLC9A3(NHE-3), andSLC26A6(PAT1), were identif ied as downregulated in brackishtreated kidney on mineral absorption, which conf irmed previous reports (Inokuchi et al., 2008; Hummel et al., 2011; Karaica et al., 2015). Genes that were related to water molecule secretion,AQP1(aquaporin),AQP4,AQP8, and the genes related to Na+uptake,ASBT(SLC10A2),SLC5A1(SGLT1), andSLC9A3(NHE-3), were down-regulated among 12 DEGs during the physiological process of bile secretion.The potential roles for aquaporin isoforms in water or solute transport in the gill of other f ish have been discussed (Cutler and Cramb, 2002; Nishimoto et al.,2007; Cutler et al., 2009). In agreement with previous studies observed in higher vertebrates, SLC10A2 underscores the role of the ileal Na+/bile acid cotransporter in the intestinal reclamation of bile acids (Oelkers et al., 1997). However, the mRNA of genesNPT2b,B0AT1(SLC6A19),NHE1(SLC9A1),andNBC(SLC4A5) related to Na+uptake were independent of the two physiological processes,presumably due to the diversity of gene function.

3.7 Analysis of the tissue specif icity of targeted gene expression subjected to salinity stress

Fig.7 Schematic view of inorganic ion channels and transporters involved in mineral absorption encoded by genes identif ied as diff erentially expressed in the kidney of turbot

We performed expression analyses using qPCR on crucial genes associated with biological processes or pathways functioning in hypo-osmotic acclimation.The relative expression levels as determined by realtime quantitative PCR are summarized in Fig.9. Our results exhibited slightly similar expression patterns in response to salt stress. The lowest expression level of six genes was found in intestinal tissue in which the transcripts of the six genes were accumulated at the highest quantity at salinity 30 (control group).This indicates that the intestinal tissue was vulnerable to osmotic pressure and may play a role in osmoregulation in response to salt stress. The expression trend of geneNKAα1f irst increased and then decreased, and was up-regulated to the highest level under salinity-10 treatment in all tissues except for intestinal tissue. The expression patterns of genes were diff erent between tissues and salinity and suggested the specif icity of the gene expression with organization and salinity. The highest expression level of genePRLwas found in spleen, followed by kidney and liver, while the highest expression level of genePAT1was found in the gill. As a freshwateradapting hormone, PRL mediates hydromineral balance of euryhaline f ish acclimated to freshwater(Pickford and Phillips, 1959), which is diff erent from our results. As is well known, thePAT1gene works in the absorption of Cl-in gills (Lee et al., 2006), which is consistent with transcriptome sequencing data and the results of real-time quantitative PCR.Apolipoprotein (APO) is the major protein component of high-density lipoprotein particles in serum and participates in the reverse transport of cholesterol from tissues to the liver for excretion (Fielding et al.,1972). Considering that f ish exploit lipids as the major energy source (after proteins), lipid metabolism is important for osmoregulation and homeostasis in f ish.The amplif ication of geneAPOMmRNA was found to be highest in gill, liver, and kidney, moderate in spleen, and weak in other tissues. For geneHMDH,the expression levels in spleen and kidney were signif icantly higher than in other tissues, which proved the obvious tissue specif icity of geneHMDHexpression and the importance of kidney and spleen in osmotic pressure regulation. The functions of HMDH in vertebrate endocrine systems (Goward and Nicholls, 1994) is consistent with that of turbot transcriptome analyses, which showed that the gene encoding core proteins of renin secretion changed. In adult zebraf ish, geneFABP6was expressed in intestine, ovary, liver, heart, and kidney (Alves-Costa et al., 2008). Our results showed that the expression of geneFABP6mRNA was higher in the gills, liver,spleen, and kidney, and was low in other tissues,presumably because of the f ish species and tissue specif icity of gene expression. After our analyses, we hypothesized that the kidney, gill, and spleen were vital regulating organs of osmotic pressure, and the species specif icity of genes related with osmoregulation expression showed that the osmoregulation pattern of euryhaline f ishes diff ered between species.

Fig.8 Schematic view of inorganic ion channels and transporters involved in bile secretion encoded by genes identif ied as diff erentially expressed in the kidney of turbot

4 CONCLUSION

The kidney transcriptome of turbot responded to brackish water was sequenced and characterized. In silico analysis indicated that the sequence of transcripts of various abundance and length were represented in the assembled contigs. The functional annotation analysis of transcriptomes showed that the immune system and the biological processes associated with digestion, absorption, and metabolism played an important role in osmoregulation in turbot responded to hypo-salinity. The biological processes of mineral absorption and bile secretion contributed to iono-osmoregulation and resulted in cell volume regulation and cell phenotypic plasticity. Of importance, we have analyzed and predicted signal transduction networks composed of the MAPK signaling pathway and PI3K-Akt signaling pathway with the GABAergic synapse worked in hypoosmoregulation signal transduction in turbot. To the best of our knowledge, this was the f irst NGS study that provided a transcriptome-wide perspective about the mechanisms of osmoregulation in kidneys of turbot responded to hypo-salinity. In addition, an analysis of the tissue specif ic expression prof iles of targeted genes subjected to salinity stress showed that the kidney, gill, and spleen were vital as regulating organs under osmotic pressure, and the osmoregulation patterns of euryhaline f ishes diff ered between species.Consistent with other studies, our studies indicated that many of the genes and signal pathways related with iono-osmoregulation were identif ied in the view of transcriptome. Finally, this study has provided an unprecedented global transcriptomic view of kidney iono-osmoregulation.

Fig.9 Expression of target genes in seven tissues of S. maximus after osmotic stress

5 DATA AVAILABILITY STATEMENT

The raw reads for the next generation sequencing data have been deposited in NCBI Sequence Read Archive (SRA) with the GenBank number SRP153594.