<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mar. Sci.</journal-id>
<journal-title>Frontiers in Marine Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mar. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-7745</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmars.2022.974501</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Marine Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Transcriptomic analysis of adaptive mechanisms in response to inland saline-alkaline water in the mud crab, <italic>Scylla paramamosain</italic>
</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Liang</surname>
<given-names>Guoling</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Qin</surname>
<given-names>Kangxiang</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Yuhao</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Niu</surname>
<given-names>Mingming</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Wang</surname>
<given-names>Huan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/795226"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Wang</surname>
<given-names>Chunlin</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/555673"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Mu</surname>
<given-names>Changkao</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/509762"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Lizhi</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Fangfang</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Su</surname>
<given-names>Qiong</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhu</surname>
<given-names>Ritong</given-names>
</name>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>School of Marine Science, Ningbo University</institution>, <addr-line>Ningbo</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Key Laboratory of Aquacultural Biotechnology Ministry of Education, Ningbo University</institution>, <addr-line>Ningbo</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Aquatic Technology Promotion Station, Sanmen Rural Bureau</institution>, <addr-line>Taizhou</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Guangxi Engineering Technology Research Center of Breeding of New Mariculture Varieties, Guangxi Institute of Oceanology Co., Ltd.</institution>, <addr-line>Beihai</addr-line>, <country>China</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>Aquatic Technology Promotion Station, Yanjin Rural Bureau</institution>, <addr-line>Xinxiang</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Hongyu Ma, Shantou University, China</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Ardavan Farhadi, Hainan University, China; Biswaranjan Paital, Orissa University of Agriculture and Technology, India</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Huan Wang, <email xlink:href="mailto:wanghuan1@nbu.edu.cn">wanghuan1@nbu.edu.cn</email>; Chunlin Wang, <email xlink:href="mailto:wangchunlin@nbu.edu.cn">wangchunlin@nbu.edu.cn</email>
</p>
</fn>
<fn fn-type="other" id="fn002">
<p>This article was submitted to Aquatic Physiology, a section of the journal Frontiers in Marine Science</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>02</day>
<month>08</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>9</volume>
<elocation-id>974501</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>06</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>11</day>
<month>07</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Liang, Qin, Chen, Niu, Wang, Wang, Mu, Chen, Wang, Su and Zhu</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Liang, Qin, Chen, Niu, Wang, Wang, Mu, Chen, Wang, Su and Zhu</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>
<italic>Scylla paramamosain</italic> is an important marine aquaculture crustacean on the southeastern coast of China. Given the difficulties of overcrowded coastal aquaculture space and insufficient production, inland saline-alkaline water culture has the potential to alleviate this bottleneck. This study separated 600 crabs into four groups: normal salinity (12-18 ppt, NS), coastal low salinity (3-5 ppt, CS), inland low salinity saline-alkaline (1.5 ppt, IS), and acute low salinity (23 ppt down to 3 ppt, AS), followed by a transcriptomic analysis of the gills. CS-vs-NS, IS-vs-NS, and AS-vs-NS obtained 1154, 1012, and 707 DEGs, respectively. GO analysis showed that the DEGs of the three comparison groups were mainly involved in cellular process, metabolic process, biological regulation, organelle, membrane, extracellular region, binding, and catalytic activity. The findings demonstrate that a significant number of genes are engaged in controlling metabolic processes in the osmoregulation process, and that cell membrane catalysis and different enzymes play a vital part in the environmental adaption process. KEGG enrichment analysis revealed that IS possesses a considerable number of signaling pathways that play important roles in osmoregulation when compared to NS. The cAMP signaling pathway increased the expression of CaM and Na<sup>+</sup>/K<sup>+</sup>-ATPase. These findings show that cations like Ca<sup>2+</sup>, Na<sup>+</sup>, and K<sup>+</sup> are critical for <italic>S. paramamosain</italic> to adapt to inland low salinity saline-alkaline water, and that the cAMP signaling pathway regulates their levels. This research provides a theoretical foundation for improving the saline-alkaline aquaculture technology of <italic>S. paramamosain</italic>.</p>
</abstract>
<kwd-group>
<kwd>
<italic>scylla paramamosain</italic>
</kwd>
<kwd>inland saline-alkaline water</kwd>
<kwd>transcriptomic analysis</kwd>
<kwd>low salinity</kwd>
<kwd>adaptive mechanisms</kwd>
</kwd-group>
<counts>
<fig-count count="6"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="48"/>
<page-count count="14"/>
<word-count count="5871"/>
</counts>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>Soil salinization is a global problem that affects both developing and developed countries. It affects around 20% of the world&#x2019;s agricultural land and nearly half of the world&#x2019;s irrigated lands (<xref ref-type="bibr" rid="B19">Kumari et&#xa0;al., 2021</xref>). Salinization not only affects the productivity of the land and threatens the sustainability of agriculture, but it can also lead to wetland degradation, resulting in imbalances in ecosystems that affect biodiversity (<xref ref-type="bibr" rid="B13">Herbert et&#xa0;al., 2015</xref>). This is particularly important where irrigation waters are compromised by salinity as this directly affects food production. Salinity thus intersects with major global concerns, including food security, desertification, and biodiversity protection. (<xref ref-type="bibr" rid="B3">Bal et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B12">Harper et&#xa0;al., 2021</xref>). Developing strategies to make use of saline land will be crucial for addressing the problem of meeting the challenge of providing food security for the projected global population of 9.3 billion people by 2050 (<xref ref-type="bibr" rid="B22">Liu and Wang, 2021</xref>). The advent of aquaculture technology has opened up new options for the conversion of underutilized saline land natural resources into productive resources. The growing demand for aquatic products, however, has pushed aquaculturists to favor inland aquaculture because of the restricted yield of caught fisheries, as well as the expensive cost of coastal land, tight regulations, and environmental concerns. Furthermore, inland ground saline water is unsuitable for traditional agriculture in many regions of the world, piquing aquaculturists&#x2019; interest in inland saline aquaculture (<xref ref-type="bibr" rid="B23">Li et&#xa0;al., 2017</xref>).</p>
<p>There are approximately 99.13 million hectares of saline-alkaline land areas in China, which are primarily situated in the north, northeast, and northwest portions of the country, yet, these limited land resources could be transformed into great potential for fisheries development (<xref ref-type="bibr" rid="B44">Zhang and Wang, 2021</xref>). In China&#x2019;s Henan, Shandong, Ningxia, etc., there have been fishery development and utilization of saline-alkaline water resources, which has created a new way for the development and utilization of saline-alkaline water resources (<xref ref-type="bibr" rid="B7">Chen et&#xa0;al., 2020</xref>). Inland saline-alkaline waters are rich in resources, but conventional breeding species struggle to survive due to the high salinity, alkalinity, pH values, and complex ionic composition of such waters, severely limiting the development and utilization of saline-alkaline waters. Land-based aquaculture employing inland saline-alkaline groundwater is referred to as inland saline-alkaline waters culture. The ionic composition of inland saline-alkaline waters differs from that of coastal saltwater, and identifying species that are tolerant of these changes is one of the most difficult aspects of the inland saline-alkaline water culture. <italic>Litopenaeus Vannamei</italic> (<xref ref-type="bibr" rid="B33">Roy et&#xa0;al., 2010</xref>) and Nile tilapia (<xref ref-type="bibr" rid="B46">Zhao et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B47">Zhao et&#xa0;al., 2020</xref>) are currently the most popular marine species for inland saline-alkaline water farming.</p>
<p>Mud crab, <italic>Scylla paramamosain</italic> is distributed along the coast of the South China Sea (<xref ref-type="bibr" rid="B8">Chung and Lin, 2006</xref>). Because of its abundance, fast growth rate and high market value, the species is an important aquaculture crab in China. mud crab inhabits waters of salinity around 5-33 ppt (<xref ref-type="bibr" rid="B27">Niu et&#xa0;al., 2020</xref>). Even though the mud crab adapted to survive in a wide range of salinity, changes in salinity beyond its range of osmoregulation capacity of it can lead to individual death, resulting in severe production losses (<xref ref-type="bibr" rid="B42">Xu et&#xa0;al., 2017</xref>). Changes in salinity had a significant effect on the growth and moltin of mud crab (<xref ref-type="bibr" rid="B30">Qi et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B16">Ji et&#xa0;al., 2022</xref>). <xref ref-type="bibr" rid="B45">Zhang et&#xa0;al. (2020)</xref> previously reported transcriptome analysis of <italic>S. paramamosain</italic> under salinity stress. Expression patterns under low and high salinity conditions has been compared with normal salinity to absorb the salinity stress and recognize osmoregulation-related genes of the mud crab, indicated that the mud crab might regulate some mechanisms such as metabolism, immunity responses, and osmoregulation to adapt to the alteration of the low salinity environment. Salinity also affects the nutritional quality of mud crab, it reared at high salinity had better nutritional value than that at low salinity (<xref ref-type="bibr" rid="B41">Wu et&#xa0;al., 2019</xref>). The adaptability of mud crab to environmental salinity is mainly reflected in the ability to regulate osmotic pressure and ion concentration. The gills are the most essential organ in the osmoregulation process (<xref ref-type="bibr" rid="B32">Romano and Zeng, 2012</xref>). The osmoregulation of the posterior gills allows the mud crab to maintain normal physiological function when the water salinity changes (<xref ref-type="bibr" rid="B38">Wang et&#xa0;al., 2018a</xref>; <xref ref-type="bibr" rid="B39">Wang et&#xa0;al., 2018b</xref>).</p>
<p>In 2017, our research team began investigating the aquaculture mode of <italic>S. paramamosain</italic> in the saline-alkaline land of the Yellow River basin, and in 2018 they successfully established the inland low salinity saline-alkaline water (1.5-2ppt) aquaculture technology system, with yield per acre reaching the coastal medium level, marking the first time that <italic>S. paramamosain</italic> has been successfully cultured in inland saline-alkaline. However, the molecular mechanism of the adaptation of <italic>S. paramamosain</italic> to inland low-salinity saline-alkaline water has yet to be discovered. Therefore, in the present study, <italic>S. paramamosain</italic> was cultured in normal salinity (12-18ppt), coastal low salinity (3-5ppt), and inland low salinity saline-alkaline water (1.5-2ppt), and acute low salinity (23ppt dropped to 3ppt), respectively. Then, the molecular mechanism of gill responses to inland low salinity saline-alkaline water stress by transcriptomic analysis, to further improve inland saline-alkaline water aquaculture technology.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<title>Materials and methods</title>
<sec id="s2_1">
<title>Animals and experimental conditions</title>
<p>A total of 600 crabs weighing 30 &#xb1; 5g were randomly selected from Sanmen Bay (Zhejiang, China), and each group of 50 crabs (n=3) was transferred to the following 4 experimental ponds with the same domestication conditions: normal salinity group (NS), coastal low salinity group (CS), inland low salinity saline-alkaline group (IS) and acute low salinity group (AS). NS, located in Sanmen Bay on the eastern coast of Zhejiang, China, is one of the main producing areas of mud crabs, with a salinity of 12-13ppt (as control). CS, located in Hangzhou Bay, northeastern Zhejiang, China, is affected by river inflow and precipitation, the seawater salinity is low and unstable, with a salinity of 3-5ppt. IS, located in Yanjin on the north bank of the Yellow River in Henan, China, the breeding water is saline-alkaline water, and the mud crab is successfully adapted to inland low salinity saline-alkaline water culture after domestication, with a salinity of 1.5-2ppt. AS, located at Ningbo University, Zhejiang, China, normal salinity (12-13ppt) acclimated crabs were directly transferred to diluted 3ppt seawater. Crabs were held in a 30 m<sup>3</sup> cement pond with disinfected seawater and fed once daily in the afternoon (5:00 pm) with razor clam <italic>Sinonovacula constricta</italic> at approximately 10% body weight. The experiment was conducted at ambient temperature and natural photoperiod. Water quality parameters were monitored every morning. The parameters were as followed: water temperature 25-30&#xb0;C, dissolved oxygen 5.24-7.52 mg L<sup>-1</sup>, ammonia nitrogen concentration 0.16-1.58 mg L<sup>-1</sup>, nitrite concentration 0.009-0.102 mg L<sup>-1</sup>, pH 7.8-8.6, alkalinity 115-160 mg/L CaCO<sub>3</sub>.</p>
<p>NS, CS and IS crabs were long-term cultured and sampled on day 90, respectively. According to a previous study (<xref ref-type="bibr" rid="B38">Wang et&#xa0;al., 2018a</xref>), AS was sampled at 120 h. On each sampling day, crabs were collected from the respective pond with traps and transferred to the laboratory for sampling. the crabs were soaked in iced water until anesthetized and then dissected promptly. The posterior gills of each individual were separated immediately with sterile forceps, frozen in liquid nitrogen, and stored at &#x2212;80&#xb0;C for further analysis of transcriptomic. For the collection of the posterior gills, three replicates were set up at the treatment.</p>
</sec>
<sec id="s2_2">
<title>RNA extraction and RNA-seq</title>
<p>RNA was extracted from crabs in NS, CS, IS and AS groups, respectively. Total RNA was extracted using the mirVana miRNA Isolation Kit (Ambion) following the manufacturer&#x2019;s protocol. RNA integrity was evaluated using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The samples with RNA Integrity Number (RIN) &#x2265; 7 were subjected to the subsequent analysis. The libraries were constructed using TruSeq Stranded mRNA LTSample Prep Kit (Illumina, San Diego, CA, USA) according to the manufacturer&#x2019;s instructions. Then these libraries were sequenced on the Illumina sequencing platform (HiSeqTM 2500) and 125bp/150bp paired-end reads were generated.</p>
</sec>
<sec id="s2_3">
<title>Analysis of sequencing data</title>
<p>Raw data (raw reads) were processed using Trimmomatic (<xref ref-type="bibr" rid="B6">Bolger et&#xa0;al., 2014</xref>). The reads containing ploy-N and the low-quality reads were removed to obtain the clean reads. Then the clean reads were mapped to the reference genome using hisat2 (<xref ref-type="bibr" rid="B18">Kim et&#xa0;al., 2015</xref>).</p>
</sec>
<sec id="s2_4">
<title>Analysis of differentially expressed transcript (DEGs)</title>
<p>FPKM (<xref ref-type="bibr" rid="B18">Kim et&#xa0;al., 2015</xref>) and read counts value of each transcript (protein_coding) was calculated using bowtie2 (<xref ref-type="bibr" rid="B20">Langmead and Salzberg, 2012</xref>) and eXpress (<xref ref-type="bibr" rid="B31">Roberts and Pachter, 2013</xref>). DEGs were identified using the DESeq (<xref ref-type="bibr" rid="B2">Anders and Huber, 2013</xref>) functions estimateSizeFactors and nbinomTest. p value &lt; 0.05 and foldchange &gt;2 or foldchange &lt; 0.5 was set as the threshold for significantly differential expression.</p>
</sec>
<sec id="s2_5">
<title>Cluster analysis, GO, and KEGG enrichment</title>
<p>Hierarchical cluster analysis of DEGs was performed to explore transcripts expression patterns. GO enrichment and KEGG (<xref ref-type="bibr" rid="B17">Kanehisa et&#xa0;al., 2008</xref>) pathway enrichment analysis of DEGs were respectively performed using R based on the hypergeometric distribution.</p>
</sec>
<sec id="s2_6">
<title>Quantitative real-time PCR for RNA-seq results validation</title>
<p>Nine genes were screened out for RNA-seq validation using quantitative real-time PCR (qRT-PCR) with TB Green<sup>&#xae;</sup> Premix Ex Taq<sup>&#x2122;</sup> II (Takara Biomedical Technology Co., Ltd, Beijing, China) according to the manufacturer&#x2019;s instructions.</p>
</sec>
</sec>
<sec id="s3">
<title>Results</title>
<sec id="s3_1">
<title>Sequencing analysis</title>
<p>A total of 12 samples were sequenced in this study, and a total of 80.69 G of clean data was obtained through the Illumina sequencing platform. The effective data volume of each sample is distributed at 6.18-7.02 G, the Q30 base is distributed at 93.0-93.67%, and the average GC content is 48.12% (<xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>). The clean reads were aligned to the reference genome, and the genome alignment of each sample was obtained. Based on the comparison results, protein-coding gene expression analysis was performed. DEGs were screened according to the expression levels of protein-coding genes in different samples, and there were three differential groups (CS-vs-NS, IS-vs-NS, AS-vs-NS).</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Sequencing sample data and quality control.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Sample</th>
<th valign="top" align="center">RawReads</th>
<th valign="top" align="center">RawBases</th>
<th valign="top" align="center">CleanReads</th>
<th valign="top" align="center">CleanBases</th>
<th valign="top" align="center">ValidBases</th>
<th valign="top" align="center">Q30</th>
<th valign="top" align="center">GC</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">CS_G1</td>
<td valign="top" align="center">46.24M</td>
<td valign="top" align="center">6.94G</td>
<td valign="top" align="center">45.46M</td>
<td valign="top" align="center">6.37G</td>
<td valign="top" align="center">91.87%</td>
<td valign="top" align="center">93.34%</td>
<td valign="top" align="center">47.53%</td>
</tr>
<tr>
<td valign="top" align="left">CS_G2</td>
<td valign="top" align="center">49.05M</td>
<td valign="top" align="center">7.36G</td>
<td valign="top" align="center">48.10M</td>
<td valign="top" align="center">6.70G</td>
<td valign="top" align="center">91.08%</td>
<td valign="top" align="center">93.34%</td>
<td valign="top" align="center">49.79%</td>
</tr>
<tr>
<td valign="top" align="left">CS_G3</td>
<td valign="top" align="center">49.51M</td>
<td valign="top" align="center">7.43G</td>
<td valign="top" align="center">48.67M</td>
<td valign="top" align="center">6.79G</td>
<td valign="top" align="center">91.47%</td>
<td valign="top" align="center">93.39%</td>
<td valign="top" align="center">48.23%</td>
</tr>
<tr>
<td valign="top" align="left">NS_G1</td>
<td valign="top" align="center">50.37M</td>
<td valign="top" align="center">7.56G</td>
<td valign="top" align="center">49.50M</td>
<td valign="top" align="center">6.93G</td>
<td valign="top" align="center">91.65%</td>
<td valign="top" align="center">93.24%</td>
<td valign="top" align="center">47.75%</td>
</tr>
<tr>
<td valign="top" align="left">NS_G2</td>
<td valign="top" align="center">49.62M</td>
<td valign="top" align="center">7.44G</td>
<td valign="top" align="center">48.75M</td>
<td valign="top" align="center">6.82G</td>
<td valign="top" align="center">91.56%</td>
<td valign="top" align="center">93.31%</td>
<td valign="top" align="center">47.67%</td>
</tr>
<tr>
<td valign="top" align="left">NS_G3</td>
<td valign="top" align="center">47.90M</td>
<td valign="top" align="center">7.18G</td>
<td valign="top" align="center">47.09M</td>
<td valign="top" align="center">6.64G</td>
<td valign="top" align="center">92.35%</td>
<td valign="top" align="center">93.25%</td>
<td valign="top" align="center">46.79%</td>
</tr>
<tr>
<td valign="top" align="left">IS_G1</td>
<td valign="top" align="center">50.05M</td>
<td valign="top" align="center">7.51G</td>
<td valign="top" align="center">49.12M</td>
<td valign="top" align="center">6.89G</td>
<td valign="top" align="center">91.72%</td>
<td valign="top" align="center">93.00%</td>
<td valign="top" align="center">47.07%</td>
</tr>
<tr>
<td valign="top" align="left">IS_G2</td>
<td valign="top" align="center">50.87M</td>
<td valign="top" align="center">7.63G</td>
<td valign="top" align="center">49.98M</td>
<td valign="top" align="center">7.02G</td>
<td valign="top" align="center">91.98%</td>
<td valign="top" align="center">93.36%</td>
<td valign="top" align="center">49.52%</td>
</tr>
<tr>
<td valign="top" align="left">IS_G3</td>
<td valign="top" align="center">44.98M</td>
<td valign="top" align="center">6.75G</td>
<td valign="top" align="center">44.22M</td>
<td valign="top" align="center">6.18G</td>
<td valign="top" align="center">91.62%</td>
<td valign="top" align="center">93.22%</td>
<td valign="top" align="center">47.71%</td>
</tr>
<tr>
<td valign="top" align="left">AS_G1</td>
<td valign="top" align="center">48.16M</td>
<td valign="top" align="center">7.22G</td>
<td valign="top" align="center">47.39M</td>
<td valign="top" align="center">6.66G</td>
<td valign="top" align="center">92.16%</td>
<td valign="top" align="center">93.47%</td>
<td valign="top" align="center">47.43%</td>
</tr>
<tr>
<td valign="top" align="left">AS_G2</td>
<td valign="top" align="center">49.59M</td>
<td valign="top" align="center">7.44G</td>
<td valign="top" align="center">48.83M</td>
<td valign="top" align="center">6.82G</td>
<td valign="top" align="center">91.73%</td>
<td valign="top" align="center">93.67%</td>
<td valign="top" align="center">49.00%</td>
</tr>
<tr>
<td valign="top" align="left">AS_G3</td>
<td valign="top" align="center">50.21M</td>
<td valign="top" align="center">7.53G</td>
<td valign="top" align="center">49.40M</td>
<td valign="top" align="center">6.87G</td>
<td valign="top" align="center">91.23%</td>
<td valign="top" align="center">93.64%</td>
<td valign="top" align="center">48.93%</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3_2">
<title>DEGs expression analysis</title>
<p>In terms of the degree of aggregation of the samples, the distribution distance between IS and the NS, AS, and CS groups is relatively large, but the NS, AS, and CS samples are more clustered (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1A</bold></xref>). When combined with the sample correlation analysis, the gene expression of NS, CS, and AS is similar, but there is a big difference with IS (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1B</bold></xref>). To elucidate gene expression patterns at different salinities, the numbers of DEGs were compared among the four treatment groups. Compared with NS, CS, IS, and AS obtained 1154 DEGs (689 up-regulated genes and 465 down-regulated genes), 1012 DEGs (580 up-regulated genes and 432 down-regulated genes), and 707 DEGs (340 up-regulated genes and 367down-regulated genes) (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1C</bold></xref>), respectively. Interestingly, 117 common DEGs were found in the three comparison groups (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1D</bold></xref>). Of the 117 DEGs, 11 were found without any description, 12 were not annotated in GO and KEGG, 28 were described as &#x201c;hypothetical&#x201d;, and 6 were described as uncharacterized protein, so these 57 DEGs were excluded. Finally, 60 DEGs with GO annotations were screened and plotted as a heatmap.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Differential gene expression analysis. <bold>(A)</bold> PCA diagram. <bold>(B)</bold> Cluster analysis results. <bold>(C)</bold> Bar Chart. <bold>(D)</bold> VennGraph. In <bold>(A)</bold>, dots represent samples, colors represent groups. In <bold>(B)</bold>, the horizontal axis represents the sample name, the vertical axis represents the corresponding sample name, and the color represents the size of the correlation coefficient. Foldchange&gt; = 2.00 and adjusted p-value &lt;= 0.05.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-09-974501-g001.tif"/>
</fig>
<p>The heatmap shows the expression of common DEGs between NS and CS, IS, and AS. Among the 60 common DEGs, 22 were up-regulated and 34 were down-regulated. In addition, cuticle protein AM1159 (CAFS_SP_G_142498.path1) and arthrodial cuticle protein AMP13.4 (CAFS_SP_G_46684.path1) down-regulated in CS-vs-NS and AS-vs-NS, but up-regulated in IS-vs-NS. Y+L amino acid transporter 2 (CAFS_SP_G_140820.path1) up-regulated in CS-vs-NS and IS-vs-NS, but down-regulated in AS-vs-NS. Cytochrome P450 CYP2B (CAFS_SP_G_17938.path1) down-regulated in CS-vs-NS and IS-vs-NS and up-regulated in AS-vs-NS (<xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2A</bold></xref>, <xref ref-type="table" rid="T2"><bold>Table&#xa0;2</bold></xref>). In IS-vs-NS, serine proteinase stubble (CAFS_SP_G_108715.path1, 4.45-fold) was the most up-regulated, followed by neurotrypsin-like (CAFS_SP_G_30842.path1, 3.74-fold) and crustacean calcium-binding protein 23 (CAFS_SP_G_104586.path1, 3.67-fold) (<xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2B</bold></xref>, <xref ref-type="table" rid="T2"><bold>Table&#xa0;2</bold></xref>). The expression of these common DEGs indicated that they played an important role in the adaptation of <italic>S. paramamosain</italic> to low salinity.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Clustered heatmap of common differentially expressed genes (DEGs). <bold>(A)</bold> Common DEGs expression between the three comparison groups. <bold>(B)</bold> DEGs expression of IS-vs-NS. Rows represent DEGs and each column represents one sample. The colors represent the expression level, and red and blue indicate higher and lower expression levels, respectively.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-09-974501-g002.tif"/>
</fig>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Differentially expressed genes annotation.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Gene ID</th>
<th valign="top" align="center">Description</th>
<th valign="top" align="center">log2FoldChange (CS-vs-NS)</th>
<th valign="top" align="center">log2FoldChange (IS-vs-NS)</th>
<th valign="top" align="center">log2FoldChange (AS-vs-NS)</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">CAFS_SP_G_108715.path1</td>
<td valign="top" align="left">Serine proteinase stubble [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">7.18</td>
<td valign="top" align="center">4.45</td>
<td valign="top" align="center">3.57</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_30842.path1</td>
<td valign="top" align="left">Neurotrypsin-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">5.17</td>
<td valign="top" align="center">3.74</td>
<td valign="top" align="center">2.88</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_117315.path1</td>
<td valign="top" align="left">Putative family 31 glucosidase KIAA1161 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">5.13</td>
<td valign="top" align="center">3.53</td>
<td valign="top" align="center">1.40</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_104586.path1</td>
<td valign="top" align="left">Crustacean calcium-binding protein 23 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">4.73</td>
<td valign="top" align="center">3.67</td>
<td valign="top" align="center">2.49</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_19217.path1</td>
<td valign="top" align="left">Aminopeptidase N-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">4.46</td>
<td valign="top" align="center">2.32</td>
<td valign="top" align="center">2.01</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_19196.path1</td>
<td valign="top" align="left">Aminopeptidase N [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">4.34</td>
<td valign="top" align="center">2.62</td>
<td valign="top" align="center">1.80</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_84247.path1</td>
<td valign="top" align="left">Ubiquitin carboxyl-terminal hydrolase 7 isoform X5 [<italic>Trichoplusia ni</italic>]</td>
<td valign="top" align="center">4.00</td>
<td valign="top" align="center">3.62</td>
<td valign="top" align="center">5.75</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_43189.path1</td>
<td valign="top" align="left">AN1-type zinc finger protein 4 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">3.71</td>
<td valign="top" align="center">1.97</td>
<td valign="top" align="center">1.63</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_19227.path1</td>
<td valign="top" align="left">Aminopeptidase N-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">3.45</td>
<td valign="top" align="center">1.80</td>
<td valign="top" align="center">1.75</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_27052.path1</td>
<td valign="top" align="left">Sarcosine dehydrogenase, mitochondrial [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">3.15</td>
<td valign="top" align="center">1.63</td>
<td valign="top" align="center">1.29</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_140096.path1</td>
<td valign="top" align="left">Branched-chain-amino-acid aminotransferase, cytosolic-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">3.06</td>
<td valign="top" align="center">1.61</td>
<td valign="top" align="center">1.09</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_32249.path1</td>
<td valign="top" align="left">Ubiquitin carboxyl-terminal hydrolase 7 isoform X2 [<italic>Diachasma alloeum</italic>]</td>
<td valign="top" align="center">2.64</td>
<td valign="top" align="center">2.25</td>
<td valign="top" align="center">5.06</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_41437.path1</td>
<td valign="top" align="left">Friend leukemia integration 1 transcription factor [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">2.37</td>
<td valign="top" align="center">2.20</td>
<td valign="top" align="center">2.90</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_144669.path1</td>
<td valign="top" align="left">DNA repair protein RAD51 4 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">2.35</td>
<td valign="top" align="center">2.47</td>
<td valign="top" align="center">1.72</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_13241.path1</td>
<td valign="top" align="left">Phospholipase A2 isozymes PA3A/PA3B/PA5 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">2.27</td>
<td valign="top" align="center">2.60</td>
<td valign="top" align="center">1.44</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_75070.path1</td>
<td valign="top" align="left">Putative urea transporter 1-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">2.16</td>
<td valign="top" align="center">1.93</td>
<td valign="top" align="center">1.62</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_11821.path1</td>
<td valign="top" align="left">Sulfotransferase 1 family member D1-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">2.16</td>
<td valign="top" align="center">1.20</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_31331.path1</td>
<td valign="top" align="left">Angiopoietin-4 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">1.77</td>
<td valign="top" align="center">1.74</td>
<td valign="top" align="center">3.34</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_23895.path1</td>
<td valign="top" align="left">LOW QUALITY PROTEIN: probable tRNA N6-adenosine Threonylcarbamoyltransferase, Mitochondrial [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">1.55</td>
<td valign="top" align="center">1.53</td>
<td valign="top" align="center">1.53</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_136269.path1</td>
<td valign="top" align="left">Estrogen sulfotransferase [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">1.35</td>
<td valign="top" align="center">2.15</td>
<td valign="top" align="center">2.02</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_28351.path1</td>
<td valign="top" align="left">Fatty acid synthase, partial [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">1.25</td>
<td valign="top" align="center">3.08</td>
<td valign="top" align="center">2.72</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_119084.path1</td>
<td valign="top" align="left">Translation initiation factor IF-2, mitochondrial-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">1.14</td>
<td valign="top" align="center">1.10</td>
<td valign="top" align="center">1.35</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_140820.path1</td>
<td valign="top" align="left">PREDICTED: Y+L amino acid transporter 2 [<italic>Dinoponera quadriceps</italic>]</td>
<td valign="top" align="center">1.13</td>
<td valign="top" align="center">1.08</td>
<td valign="top" align="center">-1.47</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_8028.path1</td>
<td valign="top" align="left">Netrin-1 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.02</td>
<td valign="top" align="center">-1.19</td>
<td valign="top" align="center">-1.02</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_15411.path1</td>
<td valign="top" align="left">laminin subunit gamma-1-like isoform X1 [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-1.08</td>
<td valign="top" align="center">-1.26</td>
<td valign="top" align="center">-1.06</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_55274.path1</td>
<td valign="top" align="left">Peptidase S9 prolyl oligopeptidase catalytic domain [<italic>Trinorchestia longiramus</italic>]</td>
<td valign="top" align="center">-1.12</td>
<td valign="top" align="center">-1.58</td>
<td valign="top" align="center">-1.23</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_9506.path1</td>
<td valign="top" align="left">Zinc finger protein [<italic>Armadillidium vulgare</italic>]</td>
<td valign="top" align="center">-1.14</td>
<td valign="top" align="center">-1.37</td>
<td valign="top" align="center">-1.14</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_71311.path1</td>
<td valign="top" align="left">Colorectal mutant cancer protein [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.17</td>
<td valign="top" align="center">-1.61</td>
<td valign="top" align="center">-1.34</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_86701.path1</td>
<td valign="top" align="left">Multiple epidermal growth factor-like domains protein 6 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.19</td>
<td valign="top" align="center">-2.24</td>
<td valign="top" align="center">-2.19</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_53922.path1</td>
<td valign="top" align="left">IDLSRF-like peptide [<italic>Cephus cinctus</italic>]</td>
<td valign="top" align="center">-1.21</td>
<td valign="top" align="center">-2.20</td>
<td valign="top" align="center">-1.36</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_117797.path1</td>
<td valign="top" align="left">Multiple epidermal growth factor-like domains protein 6 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.25</td>
<td valign="top" align="center">-2.62</td>
<td valign="top" align="center">-1.59</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_58422.path1</td>
<td valign="top" align="left">Putative metalloprotease TIKI1-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-1.25</td>
<td valign="top" align="center">-1.40</td>
<td valign="top" align="center">-1.32</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_60795.path1</td>
<td valign="top" align="left">Protein sprint [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.30</td>
<td valign="top" align="center">-1.23</td>
<td valign="top" align="center">-1.54</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_75487.path1</td>
<td valign="top" align="left">PREDICTED: NACHT and WD repeat domain-containing Protein 2-like [<italic>Hyalella azteca</italic>]</td>
<td valign="top" align="center">-1.31</td>
<td valign="top" align="center">-1.69</td>
<td valign="top" align="center">-1.38</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_24383.path1</td>
<td valign="top" align="left">Amino acid/polyamine transporter I [<italic>Trinorchestia longiramus</italic>]</td>
<td valign="top" align="center">-1.38</td>
<td valign="top" align="center">-1.02</td>
<td valign="top" align="center">-1.43</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_33009.path1</td>
<td valign="top" align="left">Cubilin [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-1.49</td>
<td valign="top" align="center">-1.68</td>
<td valign="top" align="center">-1.94</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_64043.path1</td>
<td valign="top" align="left">Gamma-butyrobetaine dioxygenase-like isoform X2 [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-1.50</td>
<td valign="top" align="center">-1.62</td>
<td valign="top" align="center">-1.45</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_85829.path1</td>
<td valign="top" align="left">Aquaporin-3-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-1.55</td>
<td valign="top" align="center">-2.91</td>
<td valign="top" align="center">-3.35</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_24480.path1</td>
<td valign="top" align="left">WD repeat-containing protein 47, partial [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-1.62</td>
<td valign="top" align="center">-1.57</td>
<td valign="top" align="center">-1.06</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_4650.path1</td>
<td valign="top" align="left">MD2 domain protein general type [<italic>Hemigrapsus sanguineus</italic>]</td>
<td valign="top" align="center">-1.71</td>
<td valign="top" align="center">-1.95</td>
<td valign="top" align="center">-1.47</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_112.path1</td>
<td valign="top" align="left">Glycogen-binding subunit 76A [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.74</td>
<td valign="top" align="center">-2.39</td>
<td valign="top" align="center">-1.50</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_3661.path1</td>
<td valign="top" align="left">DNA-directed RNA polymerase III subunit RPC4 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.78</td>
<td valign="top" align="center">-2.50</td>
<td valign="top" align="center">-1.59</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_91846.path1</td>
<td valign="top" align="left">Leucine&#x2013;tRNA ligase, cytoplasmic [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.80</td>
<td valign="top" align="center">-2.12</td>
<td valign="top" align="center">-1.51</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_66662.path1</td>
<td valign="top" align="left">Tyrosine-protein kinase RYK-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-1.90</td>
<td valign="top" align="center">-1.09</td>
<td valign="top" align="center">-2.50</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_17938.path1</td>
<td valign="top" align="left">Cytochrome P450 CYP2B [<italic>Eriocheir sinensis</italic>]</td>
<td valign="top" align="center">-1.92</td>
<td valign="top" align="center">-1.29</td>
<td valign="top" align="center">1.24</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_54393.path1</td>
<td valign="top" align="left">Protein kinase DC2 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.96</td>
<td valign="top" align="center">-1.96</td>
<td valign="top" align="center">-1.62</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_150074.path1</td>
<td valign="top" align="left">Serine protease inhibitor dipetalogastin [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-1.98</td>
<td valign="top" align="center">-1.53</td>
<td valign="top" align="center">-1.62</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_71739.path1</td>
<td valign="top" align="left">Potassium voltage-gated channel subfamily H member 1 ether-a-go-go, partial [<italic>Homarus americanus</italic>]</td>
<td valign="top" align="center">-2.00</td>
<td valign="top" align="center">-2.45</td>
<td valign="top" align="center">-1.93</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_143285.path1</td>
<td valign="top" align="left">TNFSF, partial [<italic>Procambarus clarkii</italic>]</td>
<td valign="top" align="center">-2.03</td>
<td valign="top" align="center">-1.82</td>
<td valign="top" align="center">-1.85</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_66540.path1</td>
<td valign="top" align="left">Uricase-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-2.35</td>
<td valign="top" align="center">-1.81</td>
<td valign="top" align="center">-1.86</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_59896.path1</td>
<td valign="top" align="left">Glutamine synthetase [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-2.46</td>
<td valign="top" align="center">-1.41</td>
<td valign="top" align="center">-1.49</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_46684.path1</td>
<td valign="top" align="left">Arthrodial cuticle protein AMP13.4 [<italic>Callinectes sapidus</italic>]</td>
<td valign="top" align="center">-2.69</td>
<td valign="top" align="center">1.62</td>
<td valign="top" align="center">-2.84</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_20930.path1</td>
<td valign="top" align="left">Protein turtle-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-2.75</td>
<td valign="top" align="center">-2.51</td>
<td valign="top" align="center">-1.26</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_142498.path1</td>
<td valign="top" align="left">Cuticle protein AM1159 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-2.76</td>
<td valign="top" align="center">3.49</td>
<td valign="top" align="center">-3.59</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_117186.path1</td>
<td valign="top" align="left">Multiple epidermal growth factor-like domains protein 6 [<italic>Portunus trituberculatus]</italic>
</td>
<td valign="top" align="center">-2.99</td>
<td valign="top" align="center">-2.95</td>
<td valign="top" align="center">-2.49</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_47398.path1</td>
<td valign="top" align="left">Neuropeptide F receptor [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-3.35</td>
<td valign="top" align="center">-5.31</td>
<td valign="top" align="center">-3.48</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_68360.path1</td>
<td valign="top" align="left">Nicotinamide phosphoribosyltransferase-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-3.54</td>
<td valign="top" align="center">-2.05</td>
<td valign="top" align="center">-2.01</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_84998.path1</td>
<td valign="top" align="left">Homeotic protein empty spiracles-like [<italic>Penaeus vannamei</italic>]</td>
<td valign="top" align="center">-4.07</td>
<td valign="top" align="center">-2.82</td>
<td valign="top" align="center">-3.56</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_126450.path1</td>
<td valign="top" align="left">Zinc transporter ZIP1 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-4.52</td>
<td valign="top" align="center">-3.36</td>
<td valign="top" align="center">-3.30</td>
</tr>
<tr>
<td valign="top" align="left">CAFS_SP_G_96094.path1</td>
<td valign="top" align="left">Zinc transporter ZIP1 [<italic>Portunus trituberculatus</italic>]</td>
<td valign="top" align="center">-4.78</td>
<td valign="top" align="center">-3.43</td>
<td valign="top" align="center">-3.76</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3_3">
<title>GO functional annotation</title>
<p>To better understand how different low salinity affects the osmoregulation process of <italic>S. paramamosain</italic>, all DEGs (CS-vs-NS, IS-vs-NS, AS-vs-NS) were subjected to GO analysis to determine their main functions. GO annotations were divided into three types: biological process, cellular component and molecular function.</p>
<p>The main significant processes of DEGs in CS-vs-NS group, IS-vs-NS group and AS-vs-NS group involved involved cellular process, metabolic process, biological regulation, organelle, membrane, extracellular region, binding, and catalytic activity (<xref ref-type="supplementary-material" rid="SF1"><bold>Figure S1</bold></xref>). This indicates that there are many genes involved in regulating metabolic activities in the osmoregulation process, and the catalysis of the cell membrane and various enzymes plays an important role in the process of environmental adaptation.</p>
</sec>
<sec id="s3_4">
<title>KEGG enrichment analysis of DEGs</title>
<p>We performed comparative analysis through the KEGG pathway to identify the enrichment pathways for DEGs. Differential genes are classified into six categories according to their functions, including metabolism, genetic information processing, environmental information processing, cellular processes, organismal systems, and human diseases.</p>
<p>In the CS-vs-NS group, the pathways significantly enriched in the KEGG Enrichment top 20 are mainly related to Energy metabolism, Amino acid metabolism, Carbohydrate metabolism, Carbohydrate metabolism and Lipid metabolism, such as Oxidative phosphorylation, Valine, leucine and isoleucine degradation, propanoate metabolism, citrate cycle (TCA cycle), fatty acid degradation, pyruvate metabolism, glyoxylate and dicarboxylate metabolism, glycine, serine and threonine metabolism, glycolysis/gluconeogenesis, glycerolipid metabolism and fatty acid biosynthesis (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3A</bold></xref>). In the IS-vs-NS group, the top 20 pathways showed that eight pathways were directly related to the expression regulation of the Na<sup>+</sup>/K<sup>+</sup>-ATPase: cardiac muscle contraction, thyroid hormone synthesis, pancreatic secretion, salivary secretion, protein digestion and absorption, proximal tubule bicarbonate reclamation, carbohydrate digestion and absorption, and bile secretion (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3B</bold></xref>). In the AS-vs-NS group, the pathways significantly enriched mainly include aminoacyl-tRNA biosynthesis, Nitrogen metabolism, renin-angiotensin system, antifolate resistance, ECM-receptor interaction, ABC transporters, AGE-RAGE signaling pathway in diabetes, alanine, aspartate and glutamate complications metabolism, cytosolic DNA-sensing pathway, and PI3K-Akt signaling pathway, etc. (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3C</bold></xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Pathway functional enrichment of differentially expressed genes (DEGs). <bold>(A)</bold> CS-vs-NS group. <bold>(B)</bold> IS-vs-NS group. <bold>(C)</bold> AS-vs-NS group. X axis represents enrichment score. Y axis represents pathway name. The color indicates the q value, the lower p-value indicates the more significant enrichment. Point size indicates DEG number (The larger dots refer to a larger amount).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-09-974501-g003.tif"/>
</fig>
<p>In addition, it was found that the DEGs involved in signal transduction were the most among the three comparison groups (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4</bold></xref>). In the IS-vs-NS group, the analysis of the DEGs involved in Signal transduction showed that four significantly signaling pathways were directly related to the expression regulation of the calmodulin (CaM): cGMP-PKG signaling pathway, calcium signaling pathway, cAMP signaling pathway, apelin signaling pathway (<xref ref-type="supplementary-material" rid="SF2"><bold>Figure S2</bold></xref>). Interestingly, the cAMP signaling pathway is also involved in the regulation of Na<sup>+</sup>/K<sup>+</sup>-ATPase and NMDAR, and its expression levels were significantly up-regulated by 2.115-fold and 6.064-fold, respectively (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5</bold></xref>).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>KEGG pathways classification. <bold>(A)</bold> CS-vs-NS group. <bold>(B)</bold> IS-vs-NS group. <bold>(C)</bold> AS-vs-NS group. There were six branches for KEGG pathways: Cellular Processes, Environmental Information Processing, Genetic Information Processing, Human Disease (For animals only), Metabolism and Organismal Systems.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-09-974501-g004.tif"/>
</fig>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Pathway analysis of cAMP signaling pathway (ko04024) in IS-vs-NS.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-09-974501-g005.tif"/>
</fig>
</sec>
<sec id="s3_5">
<title>Validity of DEGs in transcriptomic data</title>
<p>To validate the RNA sequencing, nine genes were selected and quantified using qPCR, namely Uricase-like (UL), TNFSF, partial (TNFSF), Homeotic protein empty spiracles-like (HPES), Amino acid/polyamine transporter I (AATI), Cytochrome P450 CYP2B(CP450), Aminopeptidase N-like (ANL), Estrogen sulfotransferase (ES), WD repeat-containing protein 47 (WDRCP), Laminin subunit gamma-1-like isoform X1 (LSGIX). The relative expression tendencies of the nine genes were in accordance with the RNA-seq (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6</bold></xref>), validating the results of the mRNA sequencing analysis.</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Candidate gene expression levels revealed by qRT-PCR and RNA-seq. NS represents normal salinity of 12-13ppt, CS represents coastal low salinity of 3-5ppt, IS represents inland low salinity saline-alkaline water of 1.5-2ppt, AS represents acute low salinity, 3ppt from normal salinity (12-13ppt) directly, which dropped. The different letters above the line chart indicate significant differences.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-09-974501-g006.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>Discussion</title>
<p>
<italic>S. paramamosain</italic> is an important marine aquaculture crustacean (<xref ref-type="bibr" rid="B37">Tang et&#xa0;al., 2020</xref>). In 2017, our research team began investigating the aquaculture mode of <italic>S. paramamosain</italic> in the saline-alkaline water of the Yellow River basin, and in 2018 we successfully established the inland low salinity saline-alkaline water (1.5-2ppt) breeding technology system, with yield per acre reaching the coastal medium level, marking the first time that <italic>S. paramamosain</italic> has been successfully cultured in inland saline-alkaline. However, the molecular mechanism by which <italic>S. paramamosain</italic> responds to saline-alkaline water stress with reduced salinity has yet to be discovered. As a result, we investigated the transcriptomic analysis of <italic>S. paramamosain</italic> gills cultivated in normal salinity, coastal low salinity, inland low salinity saline-alkaline water, and acute low salinity, respectively. To analyze the molecular mechanism of <italic>S. paramamosain</italic> responds to inland low salinity saline-alkaline water stress at the transcriptional level, and to provide theoretical reference for improving inland saline-alkaline water aquaculture technology.</p>
<p>Under low salinity seawater stress, <italic>S. paramamosain</italic> can acquire adaptive capacity through osmoregulation and reach an adaptive state after 120 h (<xref ref-type="bibr" rid="B38">Wang et&#xa0;al., 2018a</xref>). The osmoregulation of <italic>S. paramamosain</italic> is to compensate and transport ions through the posterior gills to regulate hemolymph ion concentration (<xref ref-type="bibr" rid="B43">Yao et&#xa0;al., 2020</xref>), which determines hemolymph osmotic pressure (<xref ref-type="bibr" rid="B21">Lignot et&#xa0;al., 2000</xref>), so the ionic composition and concentration of the water environment will affect its osmoregulation capacity. However, the ionic composition of saline-alkaline water differs greatly from that of seawater, and the posterior gill response mechanism may differ. The results of PCA analysis revealed that there were clear distribution differences between the IS group and the NS, AS, and CS groups, but the NS, AS, and CS groups with the same seawater background but different salinity were more clustered. When combined with the sample correlation analysis, the gene expression of NS, CS, and AS is similar, but there is a significant difference when compared to IS. The high degree of sample dispersion may be due to the differences in the osmoregulation capacity of <italic>S. paramamosain</italic> individuals under the degree of salinity stress. This indicates that the molecular mechanism of <italic>S. paramamosain</italic> adapting to saline-alkaline water may be different from the osmoregulation mechanism of adapting to low salinity seawater stress, but the osmoregulation mechanism of seawater with different salinity has a certain similarity.</p>
<p>Osmoregulation in crustaceans is an energy-consuming process (<xref ref-type="bibr" rid="B43">Yao et&#xa0;al., 2020</xref>). Free amino acids are used as energy substrates for physiological stress compensation (<xref ref-type="bibr" rid="B48">Zhou et&#xa0;al., 2011</xref>). L-serine levels in the hemolymph and muscle of <italic>Litopenaeus vannamei</italic> increased dramatically with decreasing salinity in earlier research, and osmoregulation by consuming L-serine as an energy source was demonstrated (<xref ref-type="bibr" rid="B34">Shinji et&#xa0;al., 2012</xref>). Serine is one of the major amino acids in Serine proteinase&#x2019;s catalytic active site (<xref ref-type="bibr" rid="B35">Shi et&#xa0;al., 2008</xref>). We discovered that in a low salinity environment, the expression of serine proteinase stubble was dramatically elevated. CS group was the most up-regulated (7.18-fold) compared to NS group, followed by IS group (4.45-fold) and AS group (3.45-fold) (3.57-fold). This finding may indicate serine proteinase stubble of the gills is the energy source for osmoregulation of <italic>S. paramamosain</italic>.</p>
<p>In crustaceans, gill tissue has been found to modulate ion transport, acclimatization to various salinities, and significant active areas for osmoregulation (<xref ref-type="bibr" rid="B26">Nikapitiya et&#xa0;al., 2015</xref>). Gill filament is the most fundamental functional unit of a gill, and its basic structure is made up of cuticle, subcuticular space, and gill epithelial cells, with cuticle being the most important for osmoregulation (<xref ref-type="bibr" rid="B4">Barra et&#xa0;al., 1983</xref>). The major component of the cuticle, cuticle protein, interacts with chitin to maintain the cuticle&#x2019;s structural integrity (<xref ref-type="bibr" rid="B1">Andersen, 1999</xref>). Cuticle protein is extremely permeable to the ions Na<sup>+</sup> and Cl<sup>-</sup>. The major ions that make up hemolymph&#x2019;s osmotic pressure are Na<sup>+</sup> and Cl<sup>-</sup>, which together determine the osmotic pressure of hemolymph (<xref ref-type="bibr" rid="B21">Lignot et&#xa0;al., 2000</xref>).</p>
<p>Cuticle protein AM1159 and cuticle protein AMP13.4 were down-regulated by 2.69-fold and 2.76-fold in CS group, down-regulated by 2.84-fold and 3.59-fold in AS group, but up-regulated by 3.49-fold and 1.62-fold in IS group, respectively, compared to NS group. The findings imply that cuticle protein is significant in the adaptation of <italic>S. paramamosain</italic> to inland low salinity saline-alkaline water, but its molecular mechanism needs to be investigated further.</p>
<p>Calcium signaling is important for cell growth, death, and metabolism (<xref ref-type="bibr" rid="B40">Wang et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B9">Dubois et&#xa0;al., 2016</xref>). Calmodulin (CaM) modulates cellular processes such as stress response, apoptosis, and calcium homeostasis by rapidly increasing the concentration of Ca<sup>2+</sup> on the cell membrane when cells are stimulated (<xref ref-type="bibr" rid="B36">Tadross et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B5">Berchtold and Villalobo, 2014</xref>; <xref ref-type="bibr" rid="B40">Wang et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B14">Huang et&#xa0;al., 2020</xref>). Additionally, Ca<sup>2+</sup> signal transduction may be a key signaling route for gills to cope with osmotic stress (<xref ref-type="bibr" rid="B25">Marshall et&#xa0;al., 2000</xref>; <xref ref-type="bibr" rid="B10">Fiol et&#xa0;al., 2006</xref>). CaM is involved in the regulation of Na<sup>+</sup> and Cl<sup>-</sup> transport in the posterior gills of <italic>Eriocheir sinensis</italic>, resulting in osmoregulation, according to <xref ref-type="bibr" rid="B28">P&#xe9;queux and Gilles (1992)</xref>. Na<sup>+</sup> and Cl<sup>&#x2014;</sup> transport is the most important osmoregulation method in crustaceans, which is driven by the Na<sup>+</sup>/K<sup>+</sup>-ATPase (<xref ref-type="bibr" rid="B14">Huang et&#xa0;al., 2020</xref>). The Na<sup>+</sup>/K<sup>+</sup>-ATPase is a membrane protein that connects the exchange of two external K<sup>+</sup> for three internal Na<sup>+</sup> linked to the hydrolysis of a single ATP molecule (<xref ref-type="bibr" rid="B29">Post et&#xa0;al., 1972</xref>). <xref ref-type="bibr" rid="B15">Jia and Liu (2018)</xref> found that the activity of the Na<sup>+</sup>/K<sup>+</sup>-ATPase can be controlled by cAMP. Injecting cAMP into <italic>Leptograpsus variegatus</italic> causes the gill filament Na<sup>+</sup>/K<sup>+</sup>-ATPase to activate quickly (<xref ref-type="bibr" rid="B24">Lucu et&#xa0;al., 2000</xref>). IS up-regulated CaM expression in comparison to NS <italic>via</italic> numerous signaling pathways, including the cGMP-PKG signaling system, calcium signaling pathway, cAMP signaling pathway, and apelin signaling pathway. In the cAMP signaling pathway, the expressions of Na<sup>+</sup>/K<sup>+</sup>-ATPase and NMDAR were dramatically up-regulated. It is no surprise that the osmoregulation related genes including Na<sup>+</sup>/K<sup>+</sup>-ATPase and NMDAR were upregulated under the low salinity condition. Under salinity stress, the genes related to ion transport will play a very important role in <italic>S. paramamosain.</italic> <xref ref-type="bibr" rid="B45">Zhang et&#xa0;al. (2020)</xref> analyzed transcriptomic changes in ambient salinity challenges of different salinities, indicated that under low salinity stress, <italic>S. paramamosain</italic> might cause induce changes in genes related to ion transport, such as the significant up regulation of Na<sup>+</sup>/K<sup>+</sup>-ATPase. NMDAR ion channels are highly permeable to Ca<sup>2+</sup>, Na<sup>+</sup> and K<sup>+</sup> ions (<xref ref-type="bibr" rid="B11">Hansen et&#xa0;al., 2018</xref>). Based on the above findings, it is hypothesized that the process of adapting to inland low salinity saline-alkaline water in <italic>S. paramamosain</italic> regulates the expression changes of Na<sup>+</sup>/K<sup>+</sup>-ATPase and NMDAR <italic>via</italic> the cAMP signaling pathway, improving the transport efficiency of Ca<sup>2+</sup>, Na<sup>+</sup>, and K<sup>+</sup>, and finally achieving osmotic pressure balance inside and outside the cell membrane, but the molecular mechanism needs to be investigated further. At the same time, this study also shows that cations such as Ca<sup>2+</sup>, Na<sup>+</sup>, and K<sup>+</sup> are extremely important for <italic>S. paramamosain</italic> to adapt to inland low salinity saline water.</p>
</sec>
<sec id="s5">
<title>Conclusion</title>
<p>In conclusion, we analyzed the differences in gene expression between normal salinity, coastal low salinity, inland saline-alkaline low salinity, and acute low salinity in the gills of <italic>S. paramamosain</italic>. CS-vs-NS, IS-vs-NS, and AS-vs-NS obtained 1154 (689 up-regulated genes and 465 down-regulated genes), 1012 (580 up-regulated genes and 432 down-regulated genes), and 707 (340 up-regulated genes and 367 down-regulated genes) DEGs, respectively. KEGG enrichment analysis showed that IS has numerous signaling pathways involved in the regulation of CaM and Na<sup>+</sup>/K<sup>+</sup>-ATPase when compared to NS. The results suggest that the adaptation of <italic>S. paramamosain</italic> to inland low salinity saline-alkaline water regulates the expression changes of Na<sup>+</sup>/K<sup>+</sup>-ATPase and NMDAR <italic>via</italic> the cAMP signaling pathway, improving Ca<sup>2+</sup>, Na<sup>+</sup>, and K<sup>+</sup> transport efficiency and achieving osmotic pressure balance inside and outside the cell membrane, but the molecular mechanism needs to be investigated further. At the same time, our research reveals that cations like Ca<sup>2+</sup>, Na<sup>+</sup>, and K<sup>+</sup> are critical for <italic>S. paramamosain</italic> to adapt to inland low salinity saline-alkaline water. In this study, the molecular mechanism of <italic>S. paramamosain</italic> adapting to inland low salinity saline-alkaline water was investigated for the first time, and numerous potential genes and KEGG pathways involved in regulation were revealed, to provide a scientific basis for guiding the aquaculture technology of crustaceans aquaculture.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The data presented in the study are deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) repository, accession number PRJNA851163.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>HW conceived and designed the study. GL performed the cultivation of experimental animals. GL, KQ, YC, MN, CW, CM, LC, FW, QS, and RZ performed and analyzed all other experiments. GL and HW wrote the manuscript with support from all authors. All authors read and approved the final manuscript.</p>
</sec>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>This study was sponsored by the National Key R&amp;D Program of China (2020YFD0900203), the Sanmen County Agricultural and Rural Bureau, Zhejiang Province (2020-GK088-1), Key Scientific and Technological Grant of Zhejiang for Breeding New Agricultural Varieties (2021C02069-6), Zhejiang Agricultural Science and Technology Cooperation Project (2021SNLF029).</p>
</sec>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>Author QS is employed by Guangxi Institute of Oceanology Co., Ltd.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<sec id="s11" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmars.2022.974501/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmars.2022.974501/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Image_1.tif" id="SF1" mimetype="image/tiff">
<label>Supplementary Figure&#xa0;1</label>
<caption>
<p>GO annotation of differentially expressed genes (DEGs). <bold>(A)</bold> CS-vs-NS group. <bold>(B)</bold> IS-vs-NS group. <bold>(C)</bold> AS-vs-NS group. DEGs were assigned to second-tier GO categories associated with three parent terms: biological process, cellular component, and molecular function.</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Image_2.tif" id="SF2" mimetype="image/tiff">
<label>Supplementary Figure&#xa0;2</label>
<caption>
<p>The signaling pathway involved in regulating CaM in IS-vs-NS. <bold>(A)</bold> Calcium signaling pathway, ko04020. <bold>(B)</bold> cGMP-PKG signaling pathway, ko04022. <bold>(C)</bold> Apelin signaling pathway, ko04371.</p>
</caption>
</supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Andersen</surname> <given-names>S. O.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Exoskeletal proteins from the crab, cancer pagurus</article-title>. <source>Comp. Biochem. Physiol. A Mol. Integr. Physiol.</source> <volume>123</volume> (<issue>2</issue>), <fpage>203</fpage>&#x2013;<lpage>211</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S1095-6433(99)00051-3</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Anders</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Huber</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>2013</year>). <source>Differential expression of RNA-seq data at the gene level &#x2013; the DESeq package</source> Vol. <volume>1843</volume> (<publisher-name>EMBL</publisher-name>), <fpage>398</fpage>&#x2013;<lpage>435</lpage>.</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bal</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Panda</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Pati</surname> <given-names>S. G.</given-names>
</name>
<name>
<surname>Das</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Agrawal</surname> <given-names>P. K.</given-names>
</name>
<name>
<surname>Paital</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Modulation of physiological oxidative stress and antioxidant status by abiotic factors especially salinity in aquatic organisms</article-title>. <source>Comp. Biochem. Physiol. C Toxicol. Pharmacol.</source> <volume>241</volume>, <fpage>108971</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cbpc.2020.108971</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barra</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Pequeux</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Humbert</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>1983</year>). <article-title>A morphological study on gills of a crab acclimated to fresh water</article-title>. <source>Tissue Cell.</source> <volume>15</volume> (<issue>4</issue>), <fpage>583</fpage>&#x2013;<lpage>596</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0040-8166(83)90009-5</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Berchtold</surname> <given-names>M. W.</given-names>
</name>
<name>
<surname>Villalobo</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>The many faces of calmodulin in cell proliferation, programmed cell death, autophagy, and cancer</article-title>. <source>Biochim. Biophys. Acta</source> <volume>1843</volume> (<issue>2</issue>), <fpage>398</fpage>&#x2013;<lpage>435</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.bbamcr.2013.10.021</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bolger</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Lohse</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Usadel</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Trimmomatic: a flexible trimmer for illumina sequence data</article-title>. <source>Bioinformatics</source> <volume>30</volume> (<issue>15</issue>), <fpage>2114</fpage>&#x2013;<lpage>2120</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btu170</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>X. Z.</given-names>
</name>
<name>
<surname>Lai</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Mo</surname> <given-names>Z. L.</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>H. Y.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>Z. X.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Saline alkali water green breeding technology mode</article-title>. <source>China Fish.</source> <volume>9</volume>, <fpage>61</fpage>&#x2013;<lpage>63</lpage>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chung</surname> <given-names>K. F.</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>H. C.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Osmoregulation and Na,K-ATPase expression in osmoregulatory organs of Scylla paramamosain</article-title>. <source>Comp. Biochem. Physiol. A Mol. Integr. Physiol.</source> <volume>144</volume> (<issue>1</issue>), <fpage>48</fpage>&#x2013;<lpage>57</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cbpa.2006.02.003</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dubois</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Prevarskaya</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Vanden Abeele</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The calcium-signaling toolkit: Updates needed</article-title>. <source>Biochim. Biophys. Acta</source> <volume>1863</volume> (<issue>6 Pt B</issue>), <fpage>1337</fpage>&#x2013;<lpage>1343</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.bbamcr.2015.11.033</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fiol</surname> <given-names>D. F.</given-names>
</name>
<name>
<surname>Chan</surname> <given-names>S. Y.</given-names>
</name>
<name>
<surname>K&#xfc;ltz</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Identification and pathway analysis of immediate hyperosmotic stress responsive molecular mechanisms in tilapia (<italic>Oreochromis mossambicus</italic>) gill</article-title>. <source>Comp. Biochem. Physiol. Part D Genomics Proteomics</source> <volume>1</volume> (<issue>3</issue>), <fpage>344</fpage>&#x2013;<lpage>356</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cbd.2006.08.002</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hansen</surname> <given-names>K. B.</given-names>
</name>
<name>
<surname>Yi</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Perszyk</surname> <given-names>R. E.</given-names>
</name>
<name>
<surname>Furukawa</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Wollmuth</surname> <given-names>L. P.</given-names>
</name>
<name>
<surname>Gibb</surname> <given-names>A. J.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Structure, function, and allosteric modulation of NMDA receptors</article-title>. <source>J. Gen. Physiol.</source> <volume>150</volume> (<issue>8</issue>), <fpage>1081</fpage>&#x2013;<lpage>1105</lpage>. doi: <pub-id pub-id-type="doi">10.1085/jgp.201812032</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Harper</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Dell</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Ruprecht</surname> <given-names>J. K.</given-names>
</name>
<name>
<surname>Sochacki</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Smettem</surname> <given-names>K. R. J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Salinity and the reclamation of salinized lands</article-title>. <source>Soils Landscape Restor.</source> <volume>7</volume>, <fpage>193</fpage>&#x2013;<lpage>208</lpage>. doi: <pub-id pub-id-type="doi">10.1016/B978-0-12-813193-0.00007-2</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Herbert</surname> <given-names>E. R.</given-names>
</name>
<name>
<surname>Boon</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Burgin</surname> <given-names>A. J.</given-names>
</name>
<name>
<surname>Neubauer</surname> <given-names>S. C.</given-names>
</name>
<name>
<surname>Franklin</surname> <given-names>R. B.</given-names>
</name>
<name>
<surname>Ard&#xf3;n</surname> <given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>A global perspective on wetland salinization: ecological consequences of a growing threat to freshwater wetlands</article-title>. <source>Ecosphere</source> <volume>6</volume>, <issue>10</issue>, <fpage>1</fpage>&#x2013;<lpage>43</lpage>. Wiley-Blackwell. doi: <pub-id pub-id-type="doi">10.1890/ES14-00534.1</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>Q. H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Acute salinity stress in gill Na<sup>+</sup>/K<sup>+</sup>-ATPase in the swimming crab, <italic>Portunus trituberculatus</italic> (Brachyura, portunidae)</article-title>. <source>Crustaceana</source> <volume>93</volume> (<issue>6</issue>), <fpage>587</fpage>&#x2013;<lpage>609</lpage>. doi: <pub-id pub-id-type="doi">10.1163/15685403-bja10012</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jia</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>X.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Expression of Na<sup>+</sup>/K<sup>+</sup>-ATPase was affected by salinity change in pacific abalone haliotis discus hannai</article-title>. <source>Front. Physiol.</source> <volume>9</volume>, <elocation-id>1244</elocation-id>. doi: <pub-id pub-id-type="doi">10.3389/fphys.2018.01244</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ji</surname> <given-names>Z. L.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>R. H.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C. L.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Z. Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>W. J.</given-names>
</name>
<name>
<surname>Mu</surname> <given-names>C. K.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Effect of long-term low-salinity culture on the survival, growth, and nutrient composition of mud crab Scylla paramamosain</article-title>. <source>J. Ocean Univ. China Vol.</source> <volume>21</volume>, <fpage>179</fpage>&#x2013;<lpage>185</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s11802-022-4812-x</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kanehisa</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Araki</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Goto</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Hattori</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Hirakawa</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Itoh</surname> <given-names>M.</given-names>
</name>
<etal/>
</person-group>. (<year>2008</year>). <article-title>KEGG for linking genomes to life and the environment</article-title>. <source>Nucleic Acids Res.</source> <volume>36</volume>, <fpage>480</fpage>&#x2013;<lpage>484</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkm882D480-4</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Langmead</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Salzberg</surname> <given-names>S. L.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>HISAT: a fast spliced aligner with low memory requirements</article-title>. <source>Nat. Methods</source> <volume>12</volume> (<issue>4</issue>), <fpage>357</fpage>&#x2013;<lpage>360</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.3317</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kumari</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Harikrishna</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Surasani</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Balange</surname> <given-names>A. K.</given-names>
</name>
<name>
<surname>Rani</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Growth, biochemical indices and carcass quality of red tilapia reared in zero water discharge based biofloc system in various salinities using inland saline ground water</article-title>. <source>Aquaculture</source> <volume>540</volume> (<issue>9</issue>), <fpage>736730</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.aquaculture.2021.736730</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Langmead</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Salzberg</surname> <given-names>S. L.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Fast gapped-read alignment with bowtie 2</article-title>. <source>Nat. Methods</source> <volume>9</volume> (<issue>4</issue>), <fpage>357</fpage>&#x2013;<lpage>359</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.1923</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lignot</surname> <given-names>J. H.</given-names>
</name>
<name>
<surname>Spanings-Pierrot</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Charmantier</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Osmoregulatory capacity as a tool in monitoring the physiological condition and the effect of stress in crustaceans - ScienceDirect</article-title>. <source>Aquaculture</source> <volume>191</volume> (<issue>1&#x2013;3</issue>), <fpage>209</fpage>&#x2013;<lpage>245</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0044-8486(00)00429-4</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Protection of halophytes and their uses for cultivation of saline-alkali soil in China</article-title>. <source>Biology (Basel)</source> <volume>10</volume> (<issue>5</issue>), <fpage>353</fpage>. doi: <pub-id pub-id-type="doi">10.3390/biology10050353</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X. D.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Qin</surname> <given-names>J. G.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>L. Q.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Physiological change and nutritional requirement of pacific white shrimp <italic>Litopenaeus vannamei</italic> at low salinity</article-title>. <source>Rev. Aquac.</source> <volume>9</volume>, <fpage>57</fpage>&#x2013;<lpage>75</lpage>. doi: <pub-id pub-id-type="doi">10.1111/raq.12104</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lucu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Devescovi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Skaramuca</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Kozul</surname> <given-names>V. V.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Gill Na, K-ATPase in the spiny lobster <italic>Palinurus elephas</italic> and other marine osmoconformers: Adaptiveness of enzymes from osmoconformity to hyperregulation</article-title>. <source>J. Exp. Mar. Biol. Ecol.</source> <volume>246</volume> (<issue>2</issue>), <fpage>163</fpage>&#x2013;<lpage>178</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0022-0981(99)00179-3</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Marshall</surname> <given-names>W. S.</given-names>
</name>
<name>
<surname>Bryson</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Luby</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Control of epithelial cl(-) secretion by basolateral osmolality in the euryhaline teleost fundulus heteroclitus</article-title>. <source>J. Exp. Biol.</source> <volume>203</volume> (<issue>Pt 12</issue>), <fpage>1897</fpage>&#x2013;<lpage>1905</lpage>. doi: <pub-id pub-id-type="doi">10.1242/jeb.203.12.1897</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nikapitiya</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Kim</surname> <given-names>W. S.</given-names>
</name>
<name>
<surname>Park</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kim</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>M. O.</given-names>
</name>
<name>
<surname>Kwak</surname> <given-names>I. S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Chitinase gene responses and tissue sensitivity in an intertidal mud crab (<italic>Macrophthalmus japonicus</italic>) following low or high salinity stress</article-title>. <source>Cell Stress Chaperones</source> <volume>20</volume> (<issue>3</issue>), <fpage>517</fpage>&#x2013;<lpage>526</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s12192-015-0576-1</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Niu</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>X. L.</given-names>
</name>
<name>
<surname>Ip</surname> <given-names>J. C. H.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>K. Y.</given-names>
</name>
<name>
<surname>Chu</surname> <given-names>K. H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Multi-omic approach provides insights into osmoregulation and osmoconformation of the crab Scylla paramamosain</article-title>. <source>Sci. Rep.</source> <volume>10</volume> (<issue>1</issue>), <fpage>21771</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-020-78351-w</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>P&#xe9;queux</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Gilles</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>1992</year>). <article-title>Calmodulin as a modulator of NaCl transport in the posterior salt-transporting gills of the Chinese crab eriocheir sinensis</article-title>. <source>Mar. Biol.</source> <volume>113</volume> (<issue>1</issue>), <fpage>65</fpage>&#x2013;<lpage>69</lpage>. doi: <pub-id pub-id-type="doi">10.1007/BF00367639</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Post</surname> <given-names>R. L.</given-names>
</name>
<name>
<surname>Hegyvary</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Kume</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>1972</year>). <article-title>Activation by adenosine triphosphate in the phosphorylation kinetics of sodium and potassium ion transport adenosine triphosphatase</article-title>. <source>J. Biol. Chem.</source> <volume>247</volume> (<issue>20</issue>), <fpage>6530</fpage>&#x2013;<lpage>6540</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0021-9258(19)44725-X</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qi</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Gu</surname> <given-names>X. L.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>K. J.</given-names>
</name>
<name>
<surname>Qiao</surname> <given-names>Z. G.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Effect of salinity on the survival,growth and Na<sup>+</sup>/K<sup>+</sup>-ATPase activity of early juvenile mud crabs, Scylla paramamosain</article-title>. <source>Mar. Sci.</source> <volume>37</volume> (<issue>2</issue>), <fpage>56</fpage>&#x2013;<lpage>60</lpage>.</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roberts</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Pachter</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Streaming fragment assignment for real-time analysis of sequencing experiments</article-title>. <source>Nat. Methods</source> <volume>10</volume> (<issue>1</issue>), <fpage>71</fpage>&#x2013;<lpage>73</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.2251</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Romano</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Zeng</surname> <given-names>C. S.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Osmoregulation in decapod crustaceans: implications to aquaculture productivity, methods for potential improvement and interactions with elevated ammonia exposure</article-title>. <source>Aquaculture</source> <volume>334&#x2013;337</volume>, <fpage>12</fpage>&#x2013;<lpage>23</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.aquaculture.2011.12.035</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roy</surname> <given-names>L. A.</given-names>
</name>
<name>
<surname>Davis</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Saoud</surname> <given-names>I. P.</given-names>
</name>
<name>
<surname>Boyd</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Pine</surname> <given-names>H. J.</given-names>
</name>
<name>
<surname>Boyd</surname> <given-names>C. E.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Shrimp culture in inland low salinity waters</article-title>. <source>Rev. Aquac.</source> <volume>2</volume> (<issue>4</issue>), <fpage>191</fpage>&#x2013;<lpage>208</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1753-5131.2010.01036.x</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shinji</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Okutsu</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Jayasankar</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Jasmani</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Wilder</surname> <given-names>M. N.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Metabolism of amino acids during hyposmotic adaptation in the whiteleg shrimp, litopenaeus vannamei</article-title>. <source>Amino Acids</source> <volume>43</volume> (<issue>5</issue>), <fpage>1945</fpage>&#x2013;<lpage>1954</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00726-012-1266-2</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname> <given-names>X. Z.</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>X. F.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J. X.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Molecular cloning and expression analysis of chymotrypsin-like serine protease from the Chinese shrimp, fenneropenaeus chinensis</article-title>. <source>Fish Shellfish Immunol.</source> <volume>25</volume> (<issue>5</issue>), <fpage>589</fpage>&#x2013;<lpage>597</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fsi.2008.07.011</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tadross</surname> <given-names>M. R.</given-names>
</name>
<name>
<surname>Dick</surname> <given-names>I. E.</given-names>
</name>
<name>
<surname>Yue</surname> <given-names>D. T.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Mechanism of local and global Ca<sup>2+</sup> sensing by calmodulin in complex with a Ca<sup>2+</sup> channel</article-title>. <source>Cell</source> <volume>133</volume> (<issue>7</issue>), <fpage>1228</fpage>&#x2013;<lpage>1240</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cell.2008.05.025</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C. L.</given-names>
</name>
<name>
<surname>Mu</surname> <given-names>C. K.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>H. L.</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>H. Z.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Temperature potentially induced distinctive flavor of mud crab <italic>Scylla paramamosain</italic> mediated by gut microbiota</article-title>. <source>Sci. Rep.</source> <volume>10</volume> (<issue>1</issue>), <fpage>3720</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-020-60685-0</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Tang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>H. L.</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>J. K.</given-names>
</name>
<name>
<surname>Mu</surname> <given-names>C. K.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C. L.</given-names>
</name>
</person-group> (<year>2018</year>a). <article-title>Transcriptomic analysis of adaptive mechanisms in response to sudden salinity drop in the mud crab, <italic>Scylla paramamosain</italic>
</article-title>. <source>BMC Genomics</source> <volume>19</volume> (<issue>1</issue>), <fpage>421</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12864-018-4803-x</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>H. L.</given-names>
</name>
<name>
<surname>Tang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>J. K.</given-names>
</name>
<name>
<surname>Mu</surname> <given-names>C. K.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C. L.</given-names>
</name>
</person-group> (<year>2018</year>b). <article-title>Identification and characterization of miRNAs in the gills of the mud crab (<italic>Scylla paramamosain</italic>) in response to a sudden drop in salinity</article-title>. <source>BMC Genomics</source> <volume>19</volume> (<issue>1</issue>), <fpage>609</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12864-018-4981-6</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Cadmium induces apoptosis in freshwater crab <italic>Sinopotamon henanense</italic> through activating calcium signal transduction pathway</article-title>. <source>PloS One</source> <volume>10</volume> (<issue>12</issue>), <fpage>e0144392</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0144392</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname> <given-names>Q. Y.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Fang</surname> <given-names>S. B.</given-names>
</name>
<name>
<surname>Xie</surname> <given-names>Z. F.</given-names>
</name>
<name>
<surname>Guan</surname> <given-names>M. Y.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>S. K.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Different biochemical composition and nutritional value attribute to salinity and rearing period in male and female mud crab Scylla paramamosain</article-title>. <source>Aquaculture</source> <volume>513</volume>, <fpage>734417</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.aquaculture.2019.734417</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname> <given-names>B. P.</given-names>
</name>
<name>
<surname>Tu</surname> <given-names>D. D.</given-names>
</name>
<name>
<surname>Yan</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Shu</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Shao</surname> <given-names>Q. J.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Molecular characterization of a cDNA encoding Na<sup>+</sup>/K<sup>+</sup>/2Cl<sup>-</sup> cotransporter in the gill of mud crab (<italic>Scylla paramamosain</italic>) during the molt cycle: Implication of its function in osmoregulation</article-title>. <source>Comp. Biochem. Physiol. A Mol. Integr. Physiol.</source> <volume>203</volume>, <fpage>115</fpage>&#x2013;<lpage>125</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cbpa.2016.08.019</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yao</surname> <given-names>H. Z.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Tang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C. L.</given-names>
</name>
<name>
<surname>Mu</surname> <given-names>C. K.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Metabolic mechanism of the mud crab (<italic>Scylla paramamosain</italic>) adapting to salinity sudden drop based on GC-MS technology</article-title>. <source>Aquac. Rep.</source> <volume>18</volume>, <fpage>100533</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.aqrep.2020.100533</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>B. Q.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Study on the harm of saline alkali land and its improvement technology in China</article-title>. <source>IOP Conf. Ser.: Earth Environ. Sci.</source> <volume>692</volume>, <fpage>042053</fpage>. doi: <pub-id pub-id-type="doi">10.1088/1755-1315/692/4/042053</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>Q. Y.</given-names>
</name>
<name>
<surname>Fang</surname> <given-names>S. B.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>S. K.</given-names>
</name>
<name>
<surname>Zheng</surname> <given-names>H. P.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y. L.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>mRNA profile provides novel insights into stress adaptation in mud crab megalopa, <italic>Scylla paramamosain</italic> after salinity stress</article-title>. <source>BMC Genomics</source> <volume>21</volume> (<issue>1</issue>), <fpage>559</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12864-020-06965-5</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Thammaratsuntorn</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>J. H.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Comparative transcriptome analysis of Nile tilapia (<italic>Oreochromis niloticus</italic>) in response to alkalinity stress</article-title>. <source>Genet. Mol. Res.</source> <volume>14</volume> (<issue>4</issue>), <fpage>17916</fpage>&#x2013;<lpage>17926</lpage>. doi: <pub-id pub-id-type="doi">10.4238/2015.December.22.16</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>C. S.</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Song</surname> <given-names>L. Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>J. L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Transcriptome changes for Nile tilapia (<italic>Oreochromis niloticus</italic>) in response to alkalinity stress</article-title>. <source>Comp. Biochem. Physiol. Part D Genomics Proteomics</source> <volume>33</volume>, <fpage>100651</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cbd.2019.100651</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Xian</surname> <given-names>J. A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Variation of free amino acid and carbohydrate concentrations in white shrimp, <italic>Litopenaeus vannamei</italic>: Effects of continuous cold stress</article-title>. <source>Aquaculture</source> <volume>317</volume> (<issue>1-4</issue>), <fpage>182</fpage>&#x2013;<lpage>186</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.aquaculture.2011.04.033</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>