1. Introduction
Freshwater ecosystems provide a wide range of natural resources for fisheries, aquaculture, production, and recreational activities, which are essential for human survival [1]. However, human activities such as urbanization, commercial pollution, and hydrological disturbances severely affect freshwater ecosystems such as rivers, lakes, and wetlands [2,3]. Salinization also exerts pressure on freshwater ecological resources. Global warming has led to rising sea and estuarine levels, and the continued intrusion of saline water into the upper reaches of rivers is one of the main causes of increased salinity in freshwater ecosystems [4,5,6]. Furthermore, human industrial activities, road de-icing salts, and agricultural irrigation have accelerated freshwater resource salinization [7,8]. Groundwater salt rises to the surface via capillary action in arid and semiarid areas, and salts deposited on the soil surface after water evaporation also cause soil salinization [9]. Freshwater ecosystem salinization is a global problem affecting over 100 countries and poses a significant threat to agriculture, negatively impacting biodiversity, species ecosystems, food security, and human economic welfare [10]. Therefore, studying the response of freshwater species to saline stress and the effective exploitation of saline soils are urgent and essential tasks for scientists.
Salinity is an important environmental factor affecting metabolism, growth, reproduction, and physiological processes [11,12,13]. Studies have shown that salinity in aquaculture waters above the species acclimation threshold will lead to histopathological changes in organisms and can also cause immune disorders, osmoregulatory imbalances, disturbances in the gut microbiota, and increased mortality, among other harmful effects [14,15,16]. Aquatic animals respond to the osmotic stress caused by increased salinity through heightened energy metabolism [17]. However, excessive energy metabolism alters oxidative levels, producing reactive oxygen species (ROS) accumulation, possibly triggering oxidative stress injury [18]. Furthermore, salinity stress affects the ion transport processes and transport-related enzyme activities, leading to disturbances in the osmotic regulation of salt ions in organisms [19]. Among all aquatic animals, crustaceans are the most sensitive to salinity fluctuations due to their osmoregulatory system’s inefficiency, especially in narrowly saline species [20]. Although the effects of salinity fluctuations on crustaceans have become a popular topic, most of them are only detected using macroscopic indicators (e.g., enzyme activities and changes in gene expression) [21,22,23], and the microscopic internal regulatory mechanisms are still unclear. Therefore, elucidating these regulatory mechanisms will help us understand whether crustaceans can cope with further salinity fluctuations, is important for screening salt-tolerant species of aquatic animals, and provides a scientific basis for exploring the aquaculture of aquatic species in saline soils.
High-throughput RNA sequencing technology has been widely used in various organisms and has been successfully used to analyze many aquatic species in response to salinity stress. Jiang et al. established a transcriptome database of pufferfish under long-term low-salinity stress and screened three solute carrier family (SLC) genes related to ion transport [24]. Liu et al. analyzed the liver transcript levels of halibut under low-salinity stress and showed that salinity stress affects the liver lipid metabolism of halibut [25]. Through transcriptome sequencing, Chen et al. identified various genes involved in osmoregulation under different salinity stresses in Oratosquilla oratoria [26]. Although considerable progress has been made in the studying of transcript levels in salinity-stressed aquatic animals, relatively few studies have been conducted on freshwater, narrowly saline crustaceans, especially at the transcript level.
Procambarus clarkii is a freshwater crayfish native to the United States and Mexico and was introduced to China from Japan in the 1920s [27]. Its delicious meat and high protein content has made it widely welcomed by consumers. Crayfish production was ~2.89 million tons in 2022 in China, making it the largest crustacean aquaculture species in China [28]. Therefore, in the present study, crayfish were subjected to salinity stress and investigated using transcriptome sequencing technology. The response of the crayfish to different salinity stresses, especially immunity, metabolism, ion transport, and osmoregulation, was analyzed to illustrate the mechanism of crayfish resistance to salinity. The results of this study are intended to deepen our understanding of the mechanisms by which freshwater organisms respond to salinity stress and provide useful references for the healthy culture of crayfish and the utilization of saline soils.
2. Materials and Methods
2.1. Experimental Animal
Healthy P. clarkii, with a body length of 6.0 ± 0.5 cm and a weight of 15.0 ± 2.0 g, were purchased from Jian Li, China. They were temporarily reared in the recirculating water system of the Shenzhen Test Base of the South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences (Shenzhen, China), and fed 3% of their body weight with commercial feed (China Hubei Tianbang Feed, Tianmen, China). The water temperature was maintained at 25 ± 1 °C, with a pH of 7.5 and dissolved oxygen levels of 6.5 ± 0.5 mg/L.
2.2. Exposure and Sample Collection
According to the safe salinity (6 ppt) of P. clarkii (body weight about 15 g) and our previous results of the 96 h semi-lethal salinity experiment (21 ppt) [29,30], the crayfish were marked with different labels as follows: control group (NC, 0 ppt), low-salt group (LS, 6 ppt), and high-salt group (HS, 18 ppt). A total of 270 healthy crayfish were selected and divided into three groups of 90 crayfish each. The experimental salt water was prepared in proportions of seawater and fresh water, and the water was not replaced during the test. The hepatopancreas tissues of three crayfish of each group were dissociated from the carapace at 6 h, 24 h, 72 h, respectively, and frozen in liquid nitrogen for transcriptome sequencing.
2.3. RNA Extraction, Library Preparation, and Sequencing
The Trizol reagent kit (Invitrogen, Carlsbad, CA, USA) was used to extract total RNA following the manufacturer’s protocol. RNA quality was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and RNase-free agarose gel electrophoresis. The NEBNext® UltraTM RNA Library Prep Kit (Ipswich, MA, USA) was used to construct mRNA-seq libraries for all the samples, which were then sequenced on an Illumina HiSeqTM 4000 by Gene Denovo Biotechnology Co. (Guangzhou, China). High-quality reads were produced by fastp software (version 0.18.0) to remove reads that contained adapters [31], poly nucleotides (N), and low-quality bases. Thus, contigs and unigenes were obtained by a short reads assembling program—Trinity software (version: r20140413p1) [32].
2.4. Data Analysis
2.4.1. Gene Differential Expression and Functional Annotation
The unigene expression was calculated and normalized using RPKM (reads per kb per million reads) [33]. Differential expression analysis of RNAs was performed between the two groups using the DESeq2_1.20.0 software (and between the two samples using edgeR) [34,35]. Genes were considered differentially expressed if they had a false discovery rate (FDR) parameter of less than 0.05 and an absolute fold change of ≥2. All unigenes were used to analyze gene function with Blastx alignment based on the following databases: NCBI non-redundant protein (Nr) database, Swiss-Prot protein database, Kyoto Encyclopedia of Genes, Genomes (KEGG) database, COG/KOG database, and Gene Ontology (GO) database.
2.4.2. GO, KEGG Enrichment, and Gene Expression Pattern Analysis
To explore how P. clarkii responds to salinity fluctuations, we performed GO classification and KEGG pathway annotation of all the differentially expressed genes obtained. Firstly, we used the Blast2GO program (v2.3.5) to identify the characteristics and biochemical metabolic pathways of the DEGs and their products, and the threshold parameter for significant differences was Q ≤ 0.05. Subsequently, we mapped the unigenes to the KEGG database to find the pathways to which they belonged and further determined the unigenes of the ion transport-related pathway.
The expression levels were then normalized to different values to illustrate the expression pattern of the DEGs. In brief, the expression level of each sample was normalized to 0, log2 (v1/v0), log2 (v2/v0), and the values were clustered by Short Time-series Expression Miner software 1.3.8 (STEM 1.3.8). The cluster profile was constructed by three rules as follows: At first, the maximum unit change between the different time points would be 1. Then, the maximum output profiles should be set to 20 (similar expression pattern profiles should be merged). At last, the minimum ratio of the fold change in the DEGs should be more than 2.0. The profiles with p-value ≤ 0.05 were considered as significant differences.
2.4.3. Ion Transport-Related Gene Selection and Function Prediction
To understand the differential expression of ion transport-related genes at seven salinity time points, based on the detailed information of the Nr-annotated DEGs, we selected unigenes whose annotation information was functionally annotated as “ion channel” and “transporter”. Next, the functional characteristics of all ion transport genes were described using the GeneCard database (
2.5. Verification by Quantitative Real-Time PCR (qRT-PCR)
To ensure the reliability of the RNA sequencing results, we randomly selected 9 DEGs for verification by quantitative RT-PCR at two different salinity time points. Using 18S rRNA as a reference gene, Primer Premier 6.0 software was used to design the gene-specific primer pairs (Supplementary Materials Table S1). The experiments were performed on a Roche Light Cycler 480 thermal cycler (Roche Applied Science, Penzberg, Germany). The total volume of the reaction solution was 12.5 μL containing 1 μL cDNA, 0.5 μL upstream and downstream primers, 6.25 μL 2×SYBR Premix ExTaq, and 4.25 μL sterile water. qRT-PCR was performed as follows: Five qPCR replicates for each gene, with each replicate using independently prepared RNA samples, and each sample from three technical replicates. The experimental procedure was 95 °C for 2 min. Running conditions were 95 °C for 2 min, 1 cycle, then 95 °C for 15 s, and 60 °C for 1 min, 40 cycles. Amplicons were verified by melting curve analysis and relative expression was determined using the comparative CT method (2−ΔΔCT) [36].
2.6. Statistical Analysis
All biological experiments were repeated three times independently, and the data analyzed were expressed as mean ± standard deviation (SD). Significant differences were analyzed using Duncan’s multiple comparison test followed by GraphPad Prism software 9. p < 0.05 was defined as statistically significant.
3. Results
3.1. RNA-Seq Data
A total of 884,590,210 raw reads were obtained from seven groups of crayfish hepatopancreas treated at different salinities and time points (Table 1). After data filtering, we obtained 882,552,854 clean reads, accounting for an average of 99.76% of the raw reads. The average Q20% and Q30% were 96.62% and 91.18%, respectively, and the average GC content was 51.84% (Table 1). Clean reads were assembled using Trinity software, and 52,533 unigenes were generated with an average length of 1000 bp, the maximum length of 32,384 bp, and a minimum length of 201 bp. The number and length of N50 were 8596 and 1688 bp, respectively, with an average GC content of 43.18% (Table 2). Based on the protein databases, the gene annotation results showed that 16,957 unigenes were annotated. Of these, 16,893, 8908, 12,210, 8234, and 4816 unigenes showed significant matches with sequences in the Nr, Swiss-Prot, KEGG, COG/KOG, and GO databases, respectively (Supplementary Materials Table S2).
3.2. Differential Expression Responses of P. clarkii Exposed to Different Salinity Gradients
To investigate the genes that respond to changes in salinity, we compared DEG counts between the six salinity treatment groups and the control group. Compared with the control group, 292 (184 up and 108 down), 622 (280 up and 342 down), 614 (425 up and 189 down), 324 (214 up and 110 down), 2545 (1368 up and 1177 down), and 613 (299 up and 314 down) unigenes were obtained at LS-6h, LS-24, LS-72, HS-6h, HS-24h, and HS-72h, respectively (Figure 1). We then performed GO term enrichment analysis on the DEGs from the six treatment groups to explore how salinity fluctuations affected physiological processes in crayfish (Figure 2). The results indicated that differentially expressed genes were mainly enriched in biological processes (mainly single organisms and metabolic processes), cellular components (mainly cells and cell parts), and molecular functions (mainly catalytic activity and binding).
KEGG enrichment analysis revealed the top 20 pathways in which the DEGs were enriched (Figure 3). Of the three low-salt treatment groups, DEGs were mainly enriched in pyruvate metabolism (Q = 0.010) and glycolysis/gluconeogenesis (Q = 0.035) pathways within the LS-6h group. The DEGs in the LS-24h group were primarily enriched in the pentose and glucuronate interconversion (Q = 0.000), ascorbate and aldarate metabolism (Q = 0.000), drug metabolism—cytochrome P450 (Q = 0.000), metabolism of xenobiotics by cytochrome P450 (Q = 0.000), and retinol metabolism (Q = 0.000) pathways. However, fewer DEGs were significantly enriched in the LS-72h treatment group, mainly in drug metabolism—other enzymes (Q = 0.033). In the high-salinity treatment groups, the DEGs were significantly enriched in 19 pathways at HS-24h, including the pentose and glucuronate interconversions (Q = 0.000) and metabolic pathways (Q = 0.000). Furthermore, DEGs in the ECM–receptor interaction (Q = 0.002) and ascorbate and aldarate metabolism (Q = 0.028) pathways were significantly enriched in the HS-6h and HS-72h treatment groups.
3.3. Trend Analysis of DEGs under Salinity Treatment
To visually observe the expression patterns of transcripts in relation to salinity, we analyzed the expression trends of genes in the different salinity treatment groups and performed hierarchical clustering of genes enriched in different pathways (Figure 4 and Figure 5).
Based on the normalization of gene expression trends, 20 expression patterns of 1231 genes were identified under low salt stress (Figure 4A). Among them, seven categories of gene expression patterns showed an overall post-change upregulation trend with increasing salinity stress duration, including profiles of 17 (316 genes), 15 (107 genes), 18 (61 genes), 19 (44 genes), 12 (31 genes), 10 (15 genes), and 6 (10 genes). Similarly, the following seven categories of gene expression patterns showed an overall post-change downregulation trend: profile 1 (165 genes), profile 7 (79 genes), profile 2 (61 genes), profile 13 (28 genes), profile 9 (23 genes), profile 4 (8 genes), and profile 0 (7 genes). In addition, six categories of genes, profile 16 (101 genes), profile 8 (44 genes), profile 3 (36 genes), profile 5 (27 genes), profile 14 (22 genes), and profile 11 (12 genes), showed an overall unchanged expression trend. Furthermore, 36 genes that were significantly differentially expressed in different pathways and signaling pathways were identified using cluster analysis. These included metabolic pathways (glucose, lipid, protein, vitamin, ammonia, and nitrogen metabolism) (Figure 4B), immune and antioxidant pathways, and ABC transporter signaling pathways (Figure 4C).
Similarly, 3133 genes with 20 different expression trends were characterized under high salt stress (Figure 5A). Of these, seven categories of gene expression patterns showed a general up-regulation trend, including profiles 18 (579 genes), 17 (212 genes), 12 (64 genes), 19 (33 genes), 15 (33 genes), 10 (19 genes), and 6 (7 genes). Profiles 1 (525 genes), 7 (248 genes), 2 (193 genes), 0 (151 genes), 9 (62 genes), 4 (57 genes), and 13 (6 genes) showed an overall post-change downregulation trend. Moreover, six gene categories, Profile 11 (379 genes), Profile 16 (265 genes), Profile 8 (173 genes), Profile 14 (58 genes), Profile 3 (40 genes), and Profile 5 (29 genes) showed overall unchanged expression trends. A total of 96 significant differentially expressed genes were clustered into 14 pathway classes. In addition to the metabolic, immune, and antioxidant pathways (Figure 5B), several important signaling pathways were included (Foxo, Hippo, MAPK, Notch, Wnt Toll, and Imd signaling pathways, phosphatidylinositol signaling system, and longevity regulating pathway) (Figure 5C).
3.4. Analysis of Genes Related to Ion Exchange in P. clarkii
Based on the screening of transcript data and gene annotation, we identified 36 solute carrier family (SLC) genes related to ion transport that were differentially expressed in the six salinity treatment groups (Table 3 and Figure 6). These genes include SLC41A1 (Magnesium ion transporter), SLC40A1 (Iron ion transporter), SLC39A11/9/2 (Zinc transporter), SLC30A9/7/6/5/2/1 (Zinc transporter), SLC26A11/6 (chloride–bicarbonate exchanger), SLC22A8/6/5/3/15B (Organic cation transporter), SLC13A5 (Na+-dependent carboxylate and sulfate transporter), SLC12A9/6/2 (Cation-coupled Cl-Cotransporter), SLC10A/6/3/2 (Ileal sodium/bile acid cotransporter), SLC9A9/8/3R1/2 (Sodium/hydrogen exchanger), SLC8A2/1 (Sodium–calcium exchanger), SLC5A3/1 (Sodium-dependent vitamin transporter), and SLC4A10/1 (Sodium bicarbonate cotransporter).
3.5. Quantitative Real-Time PCR
To confirm the RNA-Seq results, nine DEGs were randomly selected for RT-qPCR validation in two salinity groups (Figure 7). Among these, five genes were upregulated (SLC5A1, SOD-4, DUOX, CHER, and TPI1B) (Figure 7A,C,E,F,H), and four genes were downregulated (AO2, HSP70, RGN, and SORD) (Figure 7B,D,G,I). The gene expression trend was similar to that of the transcriptome sequencing results, indicating that the sequencing results were reliable.
4. Discussion
Salinity is an important environmental factor, and changes in salinity affect aquatic animal growth and development, metabolism, immunity, and numerous physiological processes [11,12]. In their natural habitats, crustaceans are constantly exposed to a wide range of salinities, making their survival extremely difficult because they lack effective osmoregulators [37]. Understanding the mechanisms by which crustaceans respond to environmental salinity fluctuations is an increasingly important scientific issue. Therefore, in the present study, a freshwater economic crayfish species was used as a research model, and its response pathways to different salinities and ion osmoregulation mechanisms were investigated by transcriptome sequencing.
In this study, we set a low salt concentration of 6 ppt and a high salt concentration of 18 ppt stress on crayfish. The results showed that the differentially expressed genes showed a tendency to increase and then decrease with the prolongation of time under both high and low salt stress, with the most remarkable response to the low salt stress of 24 h. Although crayfish temporarily adapted to these stresses, they still posed significant stress challenges, which explains why crayfish only survive for a short period in many lake inlets [38]. To further determine which pathways were enriched for functional genes under salinity stress, we performed GO functional enrichment analyses on six comparator groups. These showed that the upregulated and downregulated genes were significantly enriched in metabolic processes, cellular components, catalytic activity, and binding processes. Crustaceans respond to environmental stress by compensating for energy loss through metabolic pathways [39,40,41]. A study of the response of mud crabs to salinity fluctuations also showed that enzyme catalysis regulates ionic changes and osmotic pressure [42]. Numerous ion transporters bind and transport salt ions to alleviate the difference in salt ion concentration inside and outside the cell membrane [19]. Similarly, KEGG enrichment analyses showed that metabolic pathways responded most strongly to salinity stress, which was consistent with the results of the GO analyses, revealing that metabolic processes are important in the response to salinity fluctuations in crayfish. Notably, differentially expressed genes were also significantly enriched in ECM–receptor interactions under high salinity stress (Q = 0.002). Therefore, we hypothesized that high salt levels induce a cellular stress response that activates functional genes to regulate cellular responses and ultimately prevent cell death [43,44].
Under varying salinity pressures, crayfish require the synergistic expression of numerous genes to respond to the hypertonic external environment. In our study, 1132 genes were categorized into 20 different expression patterns under low salt (6 ppt) stress, while there were 3133 genes under high salt (18 ppt) stress. This demonstrates that high salt levels induce a more dramatic response and that crayfish require more genes to be differentially expressed to cope with osmotic stress [15]. Gene clustering analyses indicated that metabolic pathways were important in the response to salinity stress in crayfish. This was consistent with our GO and KEGG enrichment analyses, suggesting that metabolic pathways may be necessary for crayfish to compensate for energy loss due to hypertonic environments. Furthermore, changes in immune- and antioxidant-related genes were also found in the salinity results, possibly suggesting that salinity stress leads to immune disorders in crayfish; similar findings have been reported in various aquatic animals [15,45,46]. MAPKs are key signal transducers from the cell’s surface to the inside of the cell and are activated to control many key processes in cell physiology, including cell growth, differentiation, and adaptation to environmental pressures [47]. Our results showed that the differential expression of genes related to the MAPK signaling pathway suggests that the MAPK signaling pathway may be essential for organisms to cope with environmental stress. Two ATP-binding protein genes were differentially expressed under low-salinity stress, implying that crayfish must rapidly translocate ATP to counteract hyperosmotic stress [48]. Notably, crayfish require additional signaling pathways to respond to high-salinity stress. The FoxO signaling pathway is an important biological regulatory network, and the activation of the FoxO signaling pathway can have important roles in metabolic regulation, disease genesis, and cell proliferation [49]. We discovered that FoxO signaling pathway genes were differentially expressed at different time points during high salinity stress, suggesting that high salinity stress may lead to reduced metabolism and increased fatty acid oxidation and glycolysis in crayfish to counteract metabolic diseases; a similar report has been published on Oreochromis mossambicus [50]. The Hippo, Notch, and phosphatidylinositol signaling pathways play important roles in the cell differentiation, proliferation, and apoptosis pathways [51,52,53]. The related pathway genes in our results were also differentially expressed at different time points under high salinity stress, suggesting that high salinity stress influences the survival of cells, and that the organism adjusts the cell cycle and promotes cell apoptosis to adapt to the hyperosmotic environment outside the cell membrane. Additionally, the Wnt signaling pathway and Toll and IMD signaling pathway genes were also differentially expressed, which are related to intercellular communication and pathogen recognition, and are similarly expressed in many aquatic animals during increased environmental stress [54,55,56]. Notably, longevity regulatory pathway genes were differentially expressed under high salt stress, which is rare in crustaceans. The longevity regulatory pathway is a complex regulatory network involving numerous signaling pathways, such as the mitochondrial, insulin, and reproductive pathways [57]. It has been reported in higher animals [58] but is less frequently reported in aquatic animals in response to environmental stress. Further, our results revealed that most genes under high salt stress had a significant degree of differential change (sharply upregulated or downregulated) at HS-24h and a tendency to regress at HS-72h. This suggests that high salinity treatment of the crayfish for a period of time (at 24 h) resulted in peak stress, and that the crayfish required more rapid changes in gene expression to regulate body functions in response to the hypertonic environment and produced transient acclimatization and survival, followed by gene expression callbacks. However, failure to acclimatize to high salinity may have led to mortality in the crayfish [59]. Similar findings were reported for the transcript of Oreochromis niloticus in response to salinity stress [60].
In aquatic animals, osmotic homeostasis under salinity stress is regulated by multiple inorganic ion exchanges both inside and outside the cell membrane [61,62]. This study revealed the differential expression of genes related to ion exchange in crayfish under low and high salinity gradients in seven groups. To visualize the ion regulatory roles of these genes, we simplified all the genes related to ion exchange into a single diagram (Figure 8). Sodium (Na+) and chloride ions (Cl−) account for more constituents of the hemolymph of crustaceans [11,42,62]. Thus, when crayfish are exposed to salinity stress, osmotic pressure is balanced by the preferential uptake or secretion of Na+ and Cl−. Our results revealed that SLC12A2/6/9, which are mainly responsible for the transport of sodium chloride ions [63], were differentially expressed under salinity stress; in particular, SLC12A2 was significantly upregulated at HS-24h, suggesting that crayfish under acute salinity stress may first reduce salt loss through the active uptake of extramembrane sodium chloride and potassium ions. Notably, the upregulated expression of SLC9A9/8/2/3R1 may generate a steep electric potential difference by enhancing H+ release (hydration), which, in turn, promotes Na+ uptake [64]. Under low-salinity conditions, we found that the expression of SLC8A2/1 increased. New evidence suggests that SLC8A2/1 is involved in the removal (Ca2+) from cells, followed by the release of only 2Na+ [65]. We also found that SLC4A10/1 and SLC26A11/6 were initially upregulated and subsequently downregulated under salinity stress. Earlier studies have suggested that SLC4A10/1 may increase the co-transport of Na+ and bicarbonate (HCO3−), allowing more Na+ to be transported into the plasma [66]. Simultaneously, SLC26A11/6 can promote Cl− uptake and HCO3− secretion [67]. Specifically, a high positive transepithelial potential would be generated in the plasma owing to the rapid influx of large amounts of Na+ [68]. To drive Na+ out of the plasma, SLC4A10/1 was downregulated at HS-72h to slow this trend [62,66]. Additionally, seven sodium-dependent biomolecule co-transporter genes (SLC13/10/5 family) were also identified, among which SLC5A1 (sodium-dependent glucose transporter) expression was significantly elevated under salinity stress, which once again corroborated that rapid transmembrane transport and constant catabolism of glucose to provide ATP quickly are essential for crayfish to cope with osmotic stress [68]. Finally, several genes involved in the transport of Fe2+, Zn2+, Mg2+, and other cations were identified. Among them, Mg2+, Fe2+, and Zn2+ may activate other ion channels by altering catalytic processes or cellular structures [69,70]. Previous studies have suggested that some inorganic cations may play a role in the regulation of cell morphology during osmotic control, indirectly influencing ion secretion and uptake [71]. Overall, we hypothesized that changes in the salt stress levels dominate the ion transport patterns of P. clarkii. The synergistic action of various ion transport channels triggers a complex regulatory strategy that ultimately improves tolerance to salt stress.
5. Conclusions
Freshwater salinization is a pressing global environmental problem. Excessive salinity severely affects the metabolism, osmoregulation, behavior, and reproduction of organisms, ultimately leading to a serious loss of biodiversity in freshwater ecosystems. In this study, RNA-Seq analyses were performed on P. clarkii under acute high and low salt stress, and more differentially expressed genes were identified under high-salt stress than under low salt stress. This evidence suggests that high salt concentrations may induce more intense osmotic stress than low salt concentrations. Our results also suggest that salinity stress alters the rates of glucose, lipid, and amino acid metabolism and may ultimately enhance osmoregulation in crayfish by regulating the rates of energy metabolism and ion transmembrane transport. Moreover, immune- and antioxidant-related pathway genes were differentially expressed under salinity stress, suggesting that salinity stress induces immune disorders in crayfish. Importantly, cell proliferation, differentiation, apoptosis, and related signaling pathways (Foxo, Wnt, Hippo, Notch signaling pathway, etc.) responded to high salinity stress in crayfish compared to low-salinity stress, revealing that regulating the cell replication cycle and accelerating apoptosis may be necessary for crayfish to respond to high salinity stress. We identified 36 solute carrier family (SLC) genes related to ion transport, indicating possible ion exchange mechanisms in crayfish exposed to salinity stress. The results of the study elaborate on the response mechanism of crayfish to salinity, provide a molecular theoretical basis for the domestication of crayfish for low-salt adaptation, and provide data support for the development and utilization of saline soils.
L.L., J.-H.H., S.-G.J. and L.-S.Y. conceived and designed the experiments; L.L., L.-S.Y., Q.-B.Y., F.-L.Z. and J.-H.H. performed the experiments; S.J. and Y.-D.L. contributed to sample collection and data analysis; L.L. analyzed the data and wrote the paper; L.-S.Y. assisted with writing and proofreading. The entire study has not been published, accepted for publication, or under consideration for publication elsewhere. All authors have read and agreed to the published version of the manuscript.
All experiments in this study were approved by the Animal Care and Use Committee of the South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences (no. SCSFRI96-253) and performed according to the regulations and guidelines established by this committee.
Not applicable.
Data are contained within the article or
We thank all the staff in Laboratory of Aquaculture and Genetic breeding, South China Sea Fisheries Research Institute, for their help during the experiments. We greatly appreciate the technical assistance provided by Wen-Ya Wei for his help with qRT-PCR assays. Special thanks to Dong-Liang Liu for his help in preparing the materials for this experiment.
The authors declare no conflicts of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Figure 2. GO functional annotations of DEGs for P. clarkii exposed to salinity fluctuation; (A): NCvsHS-6h, (B): NCvsHS-24h, (C): NCvsHS-72h; (D): NCvsLS-6h, (E): NCvsLS-24h, (F): NCvsLS-72h.
Figure 3. KEGG enrichment analysis results of DEGs of P. clarkii under salt stress; (A): NCvsLS-6h, (B): NCvsLS-24h, (C): NCvsLS-72h; (D): NCvsHS-6h, (E): NCvsHS-24h, (F): NCvsHS-72h.
Figure 4. Analysis of gene expression patterns under low salt stress based on normalization of gene expression trends. (A): gene expression trend analysis; (B): gene expression clustering pattern analysis of metabolic pathways; (C): gene expression clustering pattern analysis of other signaling pathways. (The color scale on the right ranges from the lowest (blue) to the highest (red) expression level). The detailed gene information in the right range is shown in Supplementary Materials Table S3.
Figure 5. Analysis of gene expression patterns under high salt stress based on normalization of gene expression trends. (A): gene expression trend analysis; (B): gene expression clustering pattern analysis of metabolic pathways; (C): gene expression clustering pattern analysis of other signaling pathways. (The color scale on the right ranges from the lowest (blue) to the highest (red) expression level). The detailed gene information in the right range is shown in Supplementary Materials Table S4.
Figure 6. Chordal graph of 36 solute carrier family genes under different salinity treatments. N1–N36 represent the genes of Number 1–36 in Table 3, the genes correspond to the treatment groups one by one, the wider the chord width, the higher the gene expression.
Figure 7. Validation of DEGs with qRT-PCR. Relative mRNA expression level < 1 indicates downregulation, and relative mRNA expression level > 1 indicates upregulation. (A): SLC5A1, (B): AO2, (C): SOD-4, (D): HSP70, (E): DUOX, (F): CHER, (G): RGN, (H): TPI1B, (I): SORD. The data of this study are presented as mean ± SD of three parallel measurements.
Figure 8. Schematic of ion exchange-related genes in P. clarkii exposed to salinity environment.
Sequencing sample data and quality control.
Sample | Raw Data | Clean Data (%) | Q20 (%) | Q30 (%) | GC (%) |
---|---|---|---|---|---|
NC-1 | 43,079,546 | 42,982,162 (99.77%) | 96.51% | 90.75% | 51.68% |
NC-2 | 38,593,886 | 38,496,476 (99.75%) | 96.17% | 90.06% | 51.31% |
NC-3 | 43,843,540 | 43,725,484 (99.73%) | 96.80% | 91.58% | 51.91% |
LS-6h-1 | 44,961,252 | 44,810,234 (99.66%) | 96.41% | 90.82% | 51.87% |
LS-6h-2 | 40,666,280 | 40,580,188 (99.79%) | 96.61% | 91.15% | 52.15% |
LS-6h-3 | 45,471,462 | 45,335,054 (99.70%) | 96.56% | 91.11% | 50.69% |
LS-24h-1 | 36,343,018 | 36,273,260 (99.81%) | 96.75% | 91.39% | 52.11% |
LS-24h-2 | 37,333,288 | 37,262,908 (99.81%) | 96.61% | 91.05% | 51.67% |
LS-24h-3 | 43,514,756 | 43,385,520 (99.70%) | 96.48% | 91.04% | 50.90% |
LS-72h-1 | 41,937,250 | 41,860,300 (99.82%) | 96.72% | 91.26% | 50.97% |
LS-72h-2 | 39,298,830 | 39,218,856 (99.80%) | 96.96% | 91.88% | 52.03% |
LS-72h-3 | 39,176,456 | 39,097,798 (99.80%) | 96.30% | 90.59% | 52.02% |
HS-6h-1 | 45,496,514 | 45,380,462 (99.74%) | 96.26% | 90.52% | 50.76% |
HS-6h-2 | 37,971,054 | 37,894,032 (99.80%) | 96.69% | 91.36% | 53.46% |
HS-6h-3 | 44,680,448 | 44,583,988 (99.78%) | 96.63% | 90.88% | 51.51% |
HS-24h-1 | 40,967,824 | 40,872,802 (99.77%) | 96.24% | 90.51% | 52.03% |
HS-24h-2 | 41,390,576 | 41,287,814 (99.75%) | 96.22% | 90.49% | 52.74% |
HS-24h-3 | 46,509,192 | 46,401,214 (99.77%) | 98.00% | 94.48% | 52.18% |
HS-72h-1 | 41,216,070 | 41,140,894 (99.82%) | 96.42% | 90.79% | 52.28% |
HS-72h-2 | 48,341,970 | 48,254,082 (99.82%) | 96.89% | 91.75% | 52.41% |
HS-72h-3 | 43,796,998 | 43,709,326 (99.80%) | 96.70% | 91.35% | 52.01% |
Quality statistics of unigenes assembly.
Genes | GC (%) | N50 | N50 | Max Length | Min Length | Average Length | Total |
---|---|---|---|---|---|---|---|
52,533 | 43.18 | 8596 | 1688 | 32,394 | 201 | 1000 | 52,553,090 |
The 36 SLC family (solute carrier family) genes related to ion transport were differentially expressed in six salinity groups.
Number | Symbol | Description | Rpkm (Mean) | ||||||
---|---|---|---|---|---|---|---|---|---|
NC | LS- | LS- | LS- | HS- | HS- | HS- | |||
1 | SLC41A1 | Solute carrier family 41 member 1 (Magnesium ion transporter) | 1.21 | 1.11 | 1.96 | 0.96 | 1.12 | 3.69 | 0.83 |
2 | SLC40A1 | Solute carrier family 40 member 1 (Iron ion transporter) | 1.50 | 0.76 | 0.39 | 0.86 | 0.89 | 0.61 | 0.68 |
3 | SLC39A9 | Solute carrier family 39 member 9 (Zinc transporter ZIP9) | 9.99 | 11.04 | 11.65 | 11.11 | 12.65 | 12.88 | 13.12 |
4 | SLC39A11 | Solute carrier family 39 member 11 (Zinc transporter ZIP11) | 6.87 | 7.49 | 8.40 | 8.30 | 8.70 | 8.62 | 6.25 |
5 | SLC39A1 | Solute carrier family 39 member 1 (Zinc transporter ZIP1) | 169.84 | 96.51 | 63.39 | 68.11 | 90.07 | 62.21 | 83.81 |
6 | SLC30A9 | Solute carrier family 30 member 9 (Zinc transporter 9) | 1.37 | 1.12 | 1.24 | 1.58 | 1.24 | 1.10 | 1.10 |
7 | SLC30A7 | Solute carrier family 30 member 7 (Zinc transporter 7) | 1.96 | 1.11 | 1.14 | 1.67 | 1.46 | 2.03 | 1.44 |
8 | SLC30A6 | Solute carrier family 30 member 6 (Zinc transporter 6) | 0.89 | 1.05 | 1.06 | 0.94 | 1.11 | 1.38 | 0.93 |
9 | SLC30A5 | Solute carrier family 30 member 5 (Zinc transporter 5) | 0.85 | 1.20 | 1.20 | 0.94 | 1.31 | 1.30 | 1.36 |
10 | SLC30A2 | Solute carrier family 30 member 2 (Zinc transporter 2) | 14.44 | 12.71 | 10.37 | 12.75 | 14.66 | 13.42 | 9.99 |
11 | SLC30A1 | Solute carrier family 30 member 1 (Zinc transporter 1) | 2.88 | 2.39 | 1.52 | 1.82 | 2.47 | 1.63 | 1.70 |
12 | SLC26A6 | Solute carrier family 26 member 6 (chloride-bicarbonate exchanger) | 2.50 | 2.82 | 2.77 | 2.97 | 4.08 | 4.83 | 2.83 |
13 | SLC26A11 | Solute carrier family 26 member 11 (chloride-bicarbonate exchanger) | 1.89 | 2.31 | 2.34 | 2.35 | 3.34 | 2.39 | 1.55 |
14 | SLC22A8 | Solute carrier family 22 member 8 (Organic cation transporter) | 9.04 | 10.59 | 12.87 | 8.80 | 11.79 | 4.89 | 11.21 |
15 | SLC22A6 | Solute carrier family 22 member 6 (Organic cation transporter) | 22.11 | 24.18 | 14.19 | 15.98 | 17.75 | 19.86 | 14.23 |
16 | SLC22A5 | Solute carrier family 22 member 5 (Organic cation transporter) | 0.79 | 0.86 | 0.76 | 0.63 | 0.68 | 0.54 | 0.47 |
17 | SLC22A3 | Solute carrier family 22 member 3 (Organic cation transporter) | 8.86 | 11.13 | 8.70 | 6.63 | 9.70 | 12.58 | 10.05 |
18 | SLC22A15B | Solute carrier family 22 member 15B (Organic ion transporter) | 11.15 | 11.08 | 9.40 | 11.31 | 11.34 | 10.87 | 10.04 |
19 | SLC13A5 | Solute carrier family 13 member 5 (Na+-dependent carboxylate and sulfate transporter) | 15.67 | 12.67 | 9.84 | 13.95 | 11.28 | 34.18 | 15.09 |
20 | SLC12A9 | Solute carrier family 12 member 9 (Cation-coupled Cl− cotransporter) | 10.11 | 9.59 | 10.29 | 9.13 | 11.05 | 8.53 | 6.02 |
21 | SLC12A6 | Solute carrier family 12 member 6 (Cation-coupled Cl− cotransporter) | 5.00 | 4.57 | 8.13 | 5.15 | 7.42 | 5.91 | 4.65 |
22 | SLC12A2 | Solute carrier family 12 member 2 (Na+/K+/2Cl− cotransporter) | 13.86 | 25.66 | 17.36 | 11.79 | 17.32 | 74.29 | 26.06 |
23 | SLC10A6 | Solute carrier family 10 member 6 (Ileal sodium/bile acid cotransporter) | 56.69 | 37.14 | 40.56 | 46.31 | 43.40 | 16.17 | 27.24 |
24 | SLC10A3 | Solute carrier family 10 member 3 (Ileal sodium/bile acid cotransporter) | 3.08 | 2.52 | 3.42 | 1.81 | 2.65 | 2.11 | 1.50 |
25 | SLC10A2 | Solute carrier family 10 member 2 (Sodium/bile acid cotransporter) | 51.12 | 89.37 | 36.31 | 50.32 | 54.28 | 34.86 | 36.86 |
26 | SLC9A8 | Solute carrier family 9 member 8 (Sodium/hydrogen exchanger 8) | 1.12 | 1.71 | 1.74 | 1.75 | 1.84 | 3.16 | 1.95 |
27 | SLC9A9 | Solute carrier family 9 member 9 (Sodium/hydrogen exchanger 9) | 2.64 | 2.88 | 3.02 | 2.51 | 3.28 | 4.20 | 2.28 |
28 | SLC9A3R1 | Solute carrier family 9 member 3 regulator 1 (Na+/H+exchange regulatory cofactor NHE-RF1) | 29.13 | 40.50 | 32.40 | 28.00 | 37.29 | 58.53 | 32.90 |
29 | SLC9A2 | Solute carrier family 9 member 2 (Na+/H+ exchanger) | 0.13 | 0.34 | 0.27 | 0.16 | 0.54 | 5.04 | 0.56 |
30 | SLC8A2 | Solute carrier family 8 member 2 (Sodium-calcium exchanger) | 0.68 | 1.60 | 1.72 | 1.22 | 2.14 | 4.33 | 1.48 |
31 | SLC8A1 | Solute carrier family 8 member 1 (Sodium-calcium exchanger) | 0.31 | 0.63 | 0.38 | 0.62 | 0.46 | 0.50 | 0.46 |
32 | SLC5A8 | Solute carrier family 5 member 8 (Sodium-dependent vitamin transporter) | 2.32 | 2.59 | 2.65 | 2.35 | 2.78 | 3.90 | 1.77 |
33 | SLC5A3 | Solute carrier family 5 member 3 (Sodium/myo-inositol cotransporter) | 5.42 | 10.56 | 8.26 | 7.93 | 11.63 | 12.35 | 8.77 |
34 | SLC5A1 | Solute carrier family 5 member 1 (Sodium/glucose cotransporter) | 23.56 | 227.19 | 181.73 | 82.43 | 160.33 | 283.00 | 110.62 |
35 | SLC4A1 | Solute carrier family 4 member 1 (Sodium bicarbonate cotransporter) | 1.68 | 1.79 | 1.72 | 1.50 | 1.84 | 1.56 | 1.00 |
36 | SLC4A10 | Solute carrier family 4 member 10 (Sodium bicarbonate cotransporter) | 2.04 | 2.40 | 2.23 | 1.81 | 2.60 | 1.29 | 1.37 |
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Dudgeon, D.; Arthington, A.H.; Gessner, M.O.; Kawabata, Z.I.; Knowler, D.J.; Lévêque, C.; Naiman, R.J.; Prieur-Richard, A.H.; Soto, D.; Stiassny, M.L.J. et al. Freshwater Biodiversity: Importance, Threats, Status and Conservation Challenges. Biol. Rev. Camb. Philos. Soc.; 2006; 81, pp. 163-182. [DOI: https://dx.doi.org/10.1017/S1464793105006950] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16336747]
2. Dudgeon, D. Multiple Threats Imperil Freshwater Biodiversity in the Anthropocene. Curr. Biol.; 2019; 29, pp. R960-R967. [DOI: https://dx.doi.org/10.1016/j.cub.2019.08.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31593677]
3. Reid, A.J.; Carlson, A.K.; Creed, I.F.; Eliason, E.J.; Gell, P.A.; Johnson, P.T.J.; Kidd, K.A.; MacCormack, T.J.; Olden, J.D.; Ormerod, S.J. et al. Emerging Threats and Persistent Conservation Challenges for Freshwater Biodiversity. Biol. Rev.; 2019; 94, pp. 849-873. [DOI: https://dx.doi.org/10.1111/brv.12480] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30467930]
4. Jeppesen, E.; Søndergaard, M.; Pedersen, A.R.; Jürgens, K.; Strzelczak, A.; Lauridsen, T.L.; Johansson, L.S. Salinity Induced Regime Shift in Shallow Brackish Lagoons. Ecosystems; 2007; 10, pp. 47-57. [DOI: https://dx.doi.org/10.1007/s10021-006-9007-6]
5. Gardner, A.S.; Moholdt, G.; Cogley, J.G.; Wouters, B.; Arendt, A.A.; Wahr, J.; Berthier, E.; Hock, R.; Pfeffer, W.T.; Kaser, G. et al. A Reconciled Estimate of Glacier Contributions to Sea Level Rise: 2003 to 2009. Science; 2013; 340, pp. 852-857. [DOI: https://dx.doi.org/10.1126/science.1234532] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23687045]
6. Müller, B.; Gächter, R. Increasing Chloride Concentrations in Lake Constance: Characterization of Sources and Estimation of Loads. Aquat. Sci.; 2012; 74, pp. 101-112. [DOI: https://dx.doi.org/10.1007/s00027-011-0200-0]
7. Cañedo-Argüelles, M.; Kefford, B.J.; Piscart, C.; Prat, N.; Schäfer, R.B.; Schulz, C.J. Salinisation of Rivers: An Urgent Ecological Issue. Environ. Pollut.; 2013; 173, pp. 157-167. [DOI: https://dx.doi.org/10.1016/j.envpol.2012.10.011]
8. Kaushal, S.S.; Likens, G.E.; Mayer, P.M.; Shatkay, R.R.; Shelton, S.A.; Grant, S.B.; Utz, R.M.; Yaculak, A.M.; Maas, C.M.; Reimer, J.E. et al. The Anthropogenic Salt Cycle. Nat. Rev. Earth Environ.; 2023; 4, pp. 770-784. [DOI: https://dx.doi.org/10.1038/s43017-023-00485-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/38515734]
9. Cañedo-Argüelles, M.; Kefford, B.; Schäfer, R. Salt in Freshwaters: Causes, Effects and Prospects—Introduction to the Theme Issue. Philos. Trans. R. Soc. B Biol. Sci.; 2019; 374, 20180002. [DOI: https://dx.doi.org/10.1098/rstb.2018.0002]
10. Arora, S.; Singh, A.K.; Singh, Y.P. Bioremediation of Salt Affected Soils: An Indian Perspective; Springer: Berlin/Heidelberg, Germany, 2017; ISBN 9783319482576
11. Rind, K.; Beyrend, D.; Blondeau-Bidet, E.; Charmantier, G.; Cucchi, P.; Lignot, J.H. Effects of Different Salinities on the Osmoregulatory Capacity of Mediterranean Sticklebacks Living in Freshwater. J. Zool.; 2017; 303, pp. 270-280. [DOI: https://dx.doi.org/10.1111/jzo.12491]
12. Thomas, D.; Rekha, M.U.; Jani, J.R.; Sreekanth, G.B.; Thiagarajan, G.; Subburaj, R.; Kailasam, M.; Vijayan, K.K. Effects of Salinity Amendments on the Embryonic and Larval Development of a Tropical Brackishwater Ornamental Silver Moony Fish, Monodactylus argenteus (Linnaeus, 1758). Aquaculture; 2021; 544, 737073. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2021.737073]
13. Qin, Y.; Jiang, S.; Huang, J.; Zhou, F.; Yang, Q.; Jiang, S.; Yang, L. C-Type Lectin Response to Bacterial Infection and Ammonia Nitrogen Stress in Tiger Shrimp (Penaeus Monodon). Fish Shellfish Immunol.; 2019; 90, pp. 188-198. [DOI: https://dx.doi.org/10.1016/j.fsi.2019.04.034] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31028898]
14. Chen, J.; Zhou, F.; Huang, J.; Jiang, S. Ammonia and salinity tolerance of Penaeus monodon across eight breeding families. SpringerPlus; 2016; 5, 171. [DOI: https://dx.doi.org/10.1186/s40064-016-1878-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27026868]
15. Xiao, Y.; Zhang, Y.; Xu, W.; Chen, D.; Li, B.; Cheng, Y.; Guo, X. The Effects of Salinities Stress on Histopathological Changes, Serum Biochemical Index, Non-Specific Immune and Transcriptome Analysis in Red Swamp Crayfish Procambarus clarkii. Sci. Total Environ.; 2022; 840, 156502. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2022.156502] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35691360]
16. Zhao, Y.; Luo, L.; Liu, D.L.; Huang, J.H.; Jiang, S.G.; Zhou, F.L.; Yang, Q.B.; Jiang, S.; Li, Y.D.; Tan, L.Q. et al. Effect of Acute Salinity Stress on Metabolism, Antioxidant Status, and Histological Structure of Procambarus clarkii. Aquac. Res.; 2023; 2023, 2748257. [DOI: https://dx.doi.org/10.1155/2023/2748257]
17. Jacquin, L.; Petitjean, Q.; Côte, J.; Laffaille, P.; Jean, S. Effects of Pollution on Fish Behavior, Personality, and Cognition: Some Research Perspectives. Front. Ecol. Evol.; 2020; 8, 86. [DOI: https://dx.doi.org/10.3389/fevo.2020.00086]
18. Bal, A.; Panda, F.; Pati, S.G.; Das, K.; Agrawal, P.K.; Paital, B. Modulation of Physiological Oxidative Stress and Antioxidant Status by Abiotic Factors Especially Salinity in Aquatic Organisms: Redox Regulation under Salinity Stress. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol.; 2021; 241, 108971. [DOI: https://dx.doi.org/10.1016/j.cbpc.2020.108971] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33421636]
19. Lou, F.; Gao, T.; Han, Z. International Journal of Biological Macromolecules Effect of Salinity Fl Uctuation on the Transcriptome of the Japanese Mantis Shrimp Oratosquilla oratoria. Int. J. Biol. Macromol.; 2019; 140, pp. 1202-1213. [DOI: https://dx.doi.org/10.1016/j.ijbiomac.2019.08.223] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31470058]
20. Castaño-Sánchez, A.; Hose, G.C.; Reboleira, A.S.P.S. Salinity and Temperature Increase Impact Groundwater Crustaceans. Sci. Rep.; 2020; 10, 12328. [DOI: https://dx.doi.org/10.1038/s41598-020-69050-7]
21. Pan, L.Q.; Zhang, L.J.; Liu, H.Y. Effects of Salinity and PH on Ion-Transport Enzyme Activities, Survival and Growth of Litopenaeus vannamei Postlarvae. Aquaculture; 2007; 273, pp. 711-720. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2007.07.218]
22. Fabri, L.M.; Moraes, C.M.; Costa, M.I.C.; Garçon, D.P.; Fontes, C.F.L.; Pinto, M.R.; McNamara, J.C.; Leone, F.A. Salinity-Dependent Modulation by Protein Kinases and the FXYD2 Peptide of Gill (Na+, K+)-ATPase Activity in the Freshwater Shrimp Macrobrachium amazonicum (Decapoda, Palaemonidae). Biochim. Biophys. Acta-Biomembr.; 2022; 1864, 183982. [DOI: https://dx.doi.org/10.1016/j.bbamem.2022.183982] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35671812]
23. Koyama, H.; Kamiya, K.; Sasaki, Y.; Yamakawa, R.; Kuniyoshi, H.; Piyapattanakorn, S.; Watabe, S. Cloning of Glutamine Synthetase Gene from Abdominal Muscle of Kuruma Shrimp Marsupenaeus japonicus and Its Expression Profile. Fish. Sci.; 2023; 89, pp. 215-222. [DOI: https://dx.doi.org/10.1007/s12562-022-01658-2]
24. Jiang, J.L.; Xu, J.; Ye, L.; Sun, M.L.; Jiang, Z.Q.; Mao, M.G. Identification of Differentially Expressed Genes in Gills of Tiger Puffer (Takifugu rubripes) in Response to Low-Salinity Stress. Comp. Biochem. Physiol. Part B Biochem. Mol. Biol.; 2020; 243, 110437. [DOI: https://dx.doi.org/10.1016/j.cbpb.2020.110437] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32247057]
25. Liu, Z.; Ma, A.; Yuan, C.; Zhao, T.; Chang, H.; Zhang, J. Transcriptome Analysis of Liver Lipid Metabolism Disorders of the Turbot Scophthalmus maximus in Response to Low Salinity Stress. Aquaculture; 2021; 534, 736273. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2020.736273]
26. Lou, F.; Gao, T.; Han, Z. Transcriptome Analyses Reveal Alterations in Muscle Metabolism, Immune Responses and Reproductive Behavior of Japanese Mantis Shrimp (Oratosquilla oratoria) at Different Cold Temperature. Comp. Biochem. Physiol.-Part D Genom. Proteom.; 2019; 32, 100615. [DOI: https://dx.doi.org/10.1016/j.cbd.2019.100615] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31419604]
27. Huang, J.; Tang, S.; Cai, F.; Lin, Y.; Wu, Z. Microsatellite Evidence of Dispersal Mechanism of Red Swamp Crayfish (Procambarus clarkii) in the Pearl River Basin and Implications for Its Management. Sci. Rep.; 2017; 7, 8272. [DOI: https://dx.doi.org/10.1038/s41598-017-08552-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28811671]
28. Ministry of Agriculture and Rural Affairs of China. China Fishery Statistical Yearbook; Ministry of Agriculture and Rural Affairs of China: Beijing, China, 2023; 34.(In Chinese)
29. Li, H.T.; Zhou, W.Z.; Gao, H.L.; Zhang, G. Combined toxicity test of salinity and alkalinity on Procambarus clarkii. Aquaculture; 2006; 27, 1.(In Chinese)
30. Liu, D.L.; Huang, J.H.; Yang, L.S.; Tan, L.Q. Effects of multiple environmental factors on molting death of Procambarus clarkii and countermeasures. S. China Fish. Sci.; 2020; 16, pp. 29-35. (In Chinese) [DOI: https://dx.doi.org/10.12131/20190128]
31. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. Fastp: An Ultra-Fast All-in-One FASTQ Preprocessor. Bioinformatics; 2018; 34, pp. i884-i890. [DOI: https://dx.doi.org/10.1093/bioinformatics/bty560]
32. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q. et al. Full-Length Transcriptome Assembly from RNA-Seq Data without a Reference Genome. Nat. Biotechnol.; 2011; 29, pp. 644-652. [DOI: https://dx.doi.org/10.1038/nbt.1883]
33. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and Quantifying Mammalian Transcriptomes by RNA-Seq. Nat. Methods; 2008; 5, pp. 621-628. [DOI: https://dx.doi.org/10.1038/nmeth.1226] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18516045]
34. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol.; 2014; 15, 550. [DOI: https://dx.doi.org/10.1186/s13059-014-0550-8]
35. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics; 2009; 26, pp. 139-140. [DOI: https://dx.doi.org/10.1093/bioinformatics/btp616]
36. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods; 2001; 25, pp. 402-408. [DOI: https://dx.doi.org/10.1006/meth.2001.1262]
37. Chen, K.; Li, E.; Li, T.; Xu, C.; Wang, X.; Lin, H.; Qin, J.G.; Chen, L. Transcriptome and Molecular Pathway Analysis of the Hepatopancreas in the Pacific White Shrimp Litopenaeus vannamei under Chronic Low-Salinity Stress. PLoS ONE; 2015; 10, 0131503. [DOI: https://dx.doi.org/10.1371/journal.pone.0131503] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26147449]
38. Bissattini, A.M.; Traversetti, L.; Bellavia, G.; Scalici, M. Tolerance of Increasing Water Salinity in the Red Swamp Crayfish Procambarus clarkii (Girard, 1852). J. Crustac. Biol.; 2015; 35, pp. 682-685. [DOI: https://dx.doi.org/10.1163/1937240X-00002366]
39. Luo, L.; Huang, J.H.; Liu, D.L.; Jiang, S.G.; Zhou, F.L.; Jiang, S.; Yang, Q.B.; Li, Y.D.; Li, T.; Tan, L.Q. et al. Transcriptome Reveals the Important Role of Metabolic Imbalances, Immune Disorders and Apoptosis in the Treatment of Procambarus clarkii at Super High Temperature. Comp. Biochem. Physiol.-Part D Genom. Proteom.; 2021; 37, 100781. [DOI: https://dx.doi.org/10.1016/j.cbd.2020.100781]
40. Ye, L.; Jiang, S.; Zhu, X.; Yang, Q.; Wen, W.; Wu, K. Effects of salinity on growth and energy budget of juvenile Penaeus monodon. Aquaculture; 2009; 290, pp. 140-144. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2009.01.028]
41. Luo, L.; Huang, J.H.; Liu, D.L.; Jiang, S.G.; Zhou, F.L.; Jiang, S.; Yang, Q.B.; Li, Y.D.; Li, T.; Tan, L.Q. et al. Comparative Transcriptome Analysis of Differentially Expressed Genes and Pathways in Procambarus clarkii (Louisiana Crawfish) at Different Acute Temperature Stress. Genomics; 2022; 114, 110415. [DOI: https://dx.doi.org/10.1016/j.ygeno.2022.110415]
42. Wang, H.; Tang, L.; Wei, H.; Lu, J.; Mu, C.; Wang, C. Transcriptomic Analysis of Adaptive Mechanisms in Response to Sudden Salinity Drop in the Mud Crab, Scylla paramamosain. BMC Genom.; 2018; 19, 421. [DOI: https://dx.doi.org/10.1186/s12864-018-4981-6]
43. Lee, M.H.; Atkinson, S.; Murphy, G. Identification of the Extracellular Matrix (ECM) Binding Motifs of Tissue Inhibitor of Metalloproteinases (TIMP)-3 and Effective Transfer to TIMP-1. J. Biol. Chem.; 2007; 282, pp. 6887-6898. [DOI: https://dx.doi.org/10.1074/jbc.M610490200] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17202148]
44. Freedman, B.R.; Bade, N.D.; Riggin, C.N.; Zhang, S.; Haines, P.G.; Ong, K.L.; Janmey, P.A. The (Dys)Functional Extracellular Matrix. Biochim. Biophys. Acta-Mol. Cell Res.; 2015; 1853, pp. 3153-3164. [DOI: https://dx.doi.org/10.1016/j.bbamcr.2015.04.015] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25930943]
45. Cao, Q.; Wang, H.; Fan, C.; Sun, Y.; Li, J.; Cheng, J.; Chu, P.; Yin, S. Environmental Salinity Influences the Branchial Expression of TCR Pathway Related Genes Based on Transcriptome of a Catadromous Fish. Comp. Biochem. Physiol.-Part D Genom. Proteom.; 2021; 38, 100815. [DOI: https://dx.doi.org/10.1016/j.cbd.2021.100815] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33610026]
46. He, Z.; Shou, C.; Han, Z. Transcriptome Analysis of Marbled Rockfish Sebastiscus marmoratus under Salinity Stress. Animals; 2023; 13, 400. [DOI: https://dx.doi.org/10.3390/ani13030400] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36766289]
47. Fang, J.Y.; Richardson, B.C. Biological Effects of the MAPK Pathways. Lancet Oncol.; 2005; 6, pp. 322-327. [DOI: https://dx.doi.org/10.1016/S1470-2045(05)70168-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15863380]
48. Wang, L.; Wang, W.N.; Liu, Y.; Cai, D.X.; Li, J.Z.; Wang, A.L. Two Types of ATPases from the Pacific White Shrimp, Litopenaeus vannamei in Response to Environmental Stress. Mol. Biol. Rep.; 2012; 39, pp. 6427-6438. [DOI: https://dx.doi.org/10.1007/s11033-012-1461-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22311011]
49. Gutekunst, J.; Andriantsoa, R.; Falckenhayn, C.; Hanna, K.; Stein, W.; Rasamy, J.; Lyko, F. Clonal Genome Evolution and Rapid Invasive Spread of the Marbled Crayfish. Nat. Ecol. Evol.; 2018; 2, pp. 567-573. [DOI: https://dx.doi.org/10.1038/s41559-018-0467-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29403072]
50. Su, H.; Ma, D.; Fan, J.; Zhong, Z.; Li, Y.; Zhu, H. Metabolism Response Mechanism in the Gill of Oreochromis Mossambicus under Salinity, Alkalinity and Saline-Alkalinity Stresses. Ecotoxicol. Environ. Saf.; 2023; 251, 114523. [DOI: https://dx.doi.org/10.1016/j.ecoenv.2023.114523] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36638565]
51. Fu, M.; Hu, Y.; Lan, T.; Guan, K.L.; Luo, T.; Luo, M. The Hippo Signalling Pathway and Its Implications in Human Health and Diseases. Signal Transduct. Target. Ther.; 2022; 7, 376. [DOI: https://dx.doi.org/10.1038/s41392-022-01191-9]
52. Bray, S.J. Notch Signalling: A Simple Pathway Becomes Complex. Nat. Rev. Mol. Cell Biol.; 2006; 7, pp. 678-689. [DOI: https://dx.doi.org/10.1038/nrm2009]
53. Powis, G.; Phil, D. Inhibitors of Phosphatidylinositol Signalling as Antiproliferative Agents. Cancer Metastasis Rev.; 1994; 13, pp. 91-103. [DOI: https://dx.doi.org/10.1007/BF00690420] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8143348]
54. Jeffries, K.M.; Connon, R.E.; Verhille, C.E.; Dabruzzi, T.F.; Britton, M.T.; Durbin-Johnson, B.P.; Fangue, N.A. Divergent Transcriptomic Signatures in Response to Salinity Exposure in Two Populations of an Estuarine Fish. Evol. Appl.; 2019; 12, pp. 1212-1226. [DOI: https://dx.doi.org/10.1111/eva.12799] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31293632]
55. Zhou, Q.L.; Wang, L.Y.; Zhao, X.L.; Yang, Y.S.; Ma, Q.; Chen, G. Effects of Salinity Acclimation on Histological Characteristics and MiRNA Expression Profiles of Scales in Juvenile Rainbow Trout (Oncorhynchus mykiss). BMC Genom.; 2022; 23, 300. [DOI: https://dx.doi.org/10.1186/s12864-022-08531-7] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35413812]
56. Tang, D.; Shi, X.; Guo, H.; Bai, Y.; Shen, C.; Zhang, Y.; Wang, Z. Comparative Transcriptome Analysis of the Gills of Procambarus clarkii Provides Novel Insights into the Immune-Related Mechanism of Copper Stress Tolerance. Fish Shellfish Immunol.; 2020; 96, pp. 32-40. [DOI: https://dx.doi.org/10.1016/j.fsi.2019.11.060] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31786343]
57. Lu, J.Y.; Simon, M.; Zhao, Y.; Ablaeva, J.; Corson, N.; Choi, Y.; Yamada, K.L.Y.H.; Schork, N.J.; Hood, W.R.; Hill, G.E. et al. Comparative Transcriptomics Reveals Circadian and Pluripotency Networks as Two Pillars of Longevity Regulation. Cell Metab.; 2022; 34, pp. 836-856.e5. [DOI: https://dx.doi.org/10.1016/j.cmet.2022.04.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35580607]
58. Liu, L.l.; Cao, Z.h.; He, C.l.; Zhong, Y.c.; Liu, W.y.; Zhang, P.; Yang, F.; Xu, Y.j. Ferric Ion Induction of Triggering Receptor Expressed in Myeloid Cells-2 Expression and PI3K/Akt Signaling Pathway in Preosteoclast Cells to Promote Osteoclast Differentiation. Orthop. Surg.; 2020; 12, pp. 1304-1312. [DOI: https://dx.doi.org/10.1111/os.12750] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32729185]
59. Prymaczok, N.C.; Pasqualino, V.M.; Viau, V.E.; Rodríguez, E.M.; Medesani, D.A. Involvement of the Crustacean Hyperglycemic Hormone (CHH) in the Physiological Compensation of the Freshwater Crayfish Cherax quadricarinatus to Low Temperature and High Salinity Stress. J. Comp. Physiol. B Biochem. Syst. Environ. Physiol.; 2016; 186, pp. 181-191. [DOI: https://dx.doi.org/10.1007/s00360-015-0954-0]
60. Xu, Z.; Gan, L.; Li, T.; Xu, C.; Chen, K.; Wang, X.; Qin, J.G.; Chen, L.; Li, E. Transcriptome Profiling and Molecular Pathway Analysis of Genes in Association with Salinity Adaptation in Nile Tilapia Oreochromis niloticus. PLoS ONE; 2015; 10, e0136506. [DOI: https://dx.doi.org/10.1371/journal.pone.0136506] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26305564]
61. Evans, D.H.; Piermarini, P.M.; Choe, K.P. The Multifunctional Fish Gill: Dominant Site of Gas Exchange, Osmoregulation, Acid-Base Regulation, and Excretion of Nitrogenous Waste. Physiol. Rev.; 2005; 85, pp. 97-177. [DOI: https://dx.doi.org/10.1152/physrev.00050.2003]
62. Velotta, J.P.; McCormick, S.D.; O’Neill, R.J.; Schultz, E.T. Relaxed Selection Causes Microevolution of Seawater Osmoregulation and Gene Expression in Landlocked Alewives. Oecologia; 2014; 175, pp. 1081-1092. [DOI: https://dx.doi.org/10.1007/s00442-014-2961-3]
63. Arroyo, J.P.; Kahle, K.T.; Gamba, G. The SLC12 Family of Electroneutral Cation-Coupled Chloride Cotransporters. Mol. Asp. Med.; 2013; 34, pp. 288-298. [DOI: https://dx.doi.org/10.1016/j.mam.2012.05.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23506871]
64. Xu, H.; Ghishan, F.K.; Kiela, P.R. SLC9 Gene Family: Function, Expression, and Regulation. Compr. Physiol.; 2018; 8, pp. 555-583. [DOI: https://dx.doi.org/10.1002/cphy.c170027] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29687889]
65. Khananshvili, D. The SLC8 Gene Family of Sodium-Calcium Exchangers (NCX)-Structure, Function, and Regulation in Health and Disease. Mol. Asp. Med.; 2013; 34, pp. 220-235. [DOI: https://dx.doi.org/10.1016/j.mam.2012.07.003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23506867]
66. Remigante, A.; Spinelli, S.; Pusch, M.; Sarikas, A.; Morabito, R.; Marino, A.; Dossena, S. Role of SLC4 and SLC26 Solute Carriers during Oxidative Stress. Acta Physiol.; 2022; 235, e13796. [DOI: https://dx.doi.org/10.1111/apha.13796] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35143116]
67. Lam, S.H.; Lui, E.Y.; Li, Z.; Cai, S.; Sung, W.K.; Mathavan, S.; Lam, T.J.; Ip, Y.K. 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 Genom.; 2014; 15, 921. [DOI: https://dx.doi.org/10.1186/1471-2164-15-921] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25342237]
68. Vallaeys, L.; Van Biervliet, S.; De Bruyn, G.; Loeys, B.; Moring, A.S.; Van Deynse, E.; Cornette, L. Congenital Glucose-Galactose Malabsorption: A Novel Deletion within the SLC5A1 Gene. Eur. J. Pediatr.; 2013; 172, pp. 409-411. [DOI: https://dx.doi.org/10.1007/s00431-012-1802-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22843301]
69. Marshall, W.S. Mechanosensitive Signalling in Fish Gill and Other Ion Transporting Epithelia. Acta Physiol.; 2011; 202, pp. 487-499. [DOI: https://dx.doi.org/10.1111/j.1748-1716.2010.02189.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20874805]
70. Feeney, G.P.; Zheng, D.; Kille, P.; Hogstrand, C. The Phylogeny of Teleost ZIP and ZnT Zinc Transporters and Their Tissue Specific Expression and Response to Zinc in Zebrafish. Biochim. Biophys. Acta-Gene Struct. Expr.; 2005; 1732, pp. 88-95. [DOI: https://dx.doi.org/10.1016/j.bbaexp.2005.12.002]
71. Voss, J.G. Effect of Inorganic Cations on Bactericidal Activity of Anionic. J. Bacteriol.; 1963; 86, pp. 207-211. [DOI: https://dx.doi.org/10.1128/jb.86.2.207-211.1963]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Simple Summary
Procambarus clarkii is an economic freshwater aquaculture species which is popular with consumers for its delicious flavor and high protein content. Salinization of freshwater ecosystems is an increasingly pressing global issue that poses a significant threat to aquaculture. Salinity is an important environmental factor directly affecting the metabolism, growth, reproduction, and physiological processes of aquatic animals. In this study, crayfish were subjected to acute low salt (6 ppt) and high salt (18 ppt) stress and investigated using transcriptome sequencing technology. The response of the crayfish to different salinity stresses, especially immunity, metabolism, ion transport, and osmoregulation, was analyzed to illustrate the resistance mechanism of crayfish facing salt stress. The results of this study are intended to deepen our understanding of the mechanisms by which freshwater organisms respond to salinity stress and provide useful references for the healthy culture of crayfish and the utilization of saline soils.
AbstractSalinization of freshwater ecosystems is a pressing global issue. Changes in salinity can exert severe pressure on aquatic animals and jeopardize their survival. Procambarus clarkii is a valuable freshwater aquaculture species that exhibits some degree of salinity tolerance, making it an excellent research model for freshwater aquaculture species facing salinity stress. In the present study, crayfish were exposed to acute low salt (6 ppt) and high salt (18 ppt) conditions. The organisms were continuously monitored at 6, 24, and 72 h using RNA-Seq to investigate the mechanisms of salt stress resistance. Transcriptome analysis revealed that the crayfish responded to salinity stress with numerous differentially expressed genes, and most of different expression genes was observed in high salinity group for 24h. GO and KEGG enrichment analyses indicated that metabolic pathways were the primary response pathways in crayfish under salinity stress. This suggests that crayfish may use metabolic pathways to compensate for energy loss caused by osmotic stress. Furthermore, gene expression analysis revealed the differential expression of immune and antioxidant-related pathway genes under salinity stress, implying that salinity stress induces immune disorders in crayfish. More genes related to cell proliferation, differentiation, and apoptosis, such as the Foxo, Wnt, Hippo, and Notch signaling pathways, responded to high-salinity stress. This suggests that regulating the cellular replication cycle and accelerating apoptosis may be necessary for crayfish to cope with high-salinity stress. Additionally, we identified 36 solute carrier family (SLC) genes related to ion transport, depicting possible ion exchange mechanisms in crayfish under salinity stress. These findings aimed to establish a foundation for understanding crustacean responses to salinity stress and their osmoregulatory mechanisms.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details


1 Key Laboratory of South China Sea Fishery Resources Exploitation and Utilization, Ministry of Agriculture and Rural Affairs, South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou 510300, China;
2 Key Laboratory of South China Sea Fishery Resources Exploitation and Utilization, Ministry of Agriculture and Rural Affairs, South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou 510300, China;
3 Key Laboratory of South China Sea Fishery Resources Exploitation and Utilization, Ministry of Agriculture and Rural Affairs, South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou 510300, China;
4 Key Laboratory of Efficient Utilization and Processing of Marine Fishery Resources of Hainan Province, Sanya Tropical Fisheries Research Institute, Sanya 572018, China;