Proteomic Responses Under Cold Stress Reveal Unique Cold Tolerance Mechanisms in the Pacific White Shrimp (Litopenaeus vannamei)

The Pacific white shrimp (Litopenaeus vannamei), one of the most widely cultured shrimp species in the world, often suffers from cold stress. To understand the molecular mechanism of cold tolerance in Pacific white shrimp, we conducted a proteomic analysis on two contrasting shrimp cultivars, namely, cold-tolerant Guihai2 (GH2) and cold-sensitive Guihai1 (GH1), under normal temperature (28°C), under cold stress (16°C), and during recovery to 28°C. In total, 3,349 proteins were identified, among which 2,736 proteins were quantified. Based on gene ontology annotations, differentially expressed proteins largely belonged to biological processes, cellular components, and molecular functions. KEGG pathway annotations indicated that the main changes were observed in the lysosome, ribosomes, and oxidative phosphorylation. Subcellular localization analysis showed a significant increase in proteins present in cytosol, extracellular regions, and mitochondria. Combining enrichment-based clustering analysis and qRT-PCR analysis, we found that glutathione S-transferase, zinc proteinase, m7GpppX diphosphatase, AP2 transcription complex, and zinc-finger transcription factors played a major role in the cold stress response in Pacific white shrimp. Moreover, structure proteins, including different types of lectin and DAPPUDRAFT, were indispensable for cold stress tolerance of the Pacific white shrimp. Results indicate the molecular mechanisms of the Pacific white shrimp in response to cold stress and provide new insight into breeding new cultivars with increased cold tolerance.


INTRODUCTION
The Pacific white shrimp (Litopenaeus vannamei) is a widely cultured species in subtropical areas such as China as well as many other Asian countries. Many biotic factors affect the growth and yield of Pacific white shrimp, such as early mortality syndrome (EMS) (Chuchird et al., 2015), Taura syndrome virus (TSV) (Zeng et al., 2013), associated bacterial communities (Zheng et al., 2017), and Spiroplasma penaei sp. nov., a bacterial species associated with mortality (Nunan et al., 2005). Abiotic stresses also considerably influence the aquaculture of these organisms. These stresses include low temperatures, which can cause persistent effects on fish muscle (Schnurr et al., 2014), and acute ammonia stress, which can lead to the death of the Pacific white shrimp (Heisterkamp et al., 2016;Lu et al., 2016).
In recent years, farming areas in China have been significantly expanded to meet the increasing demands for the Pacific white shrimp. However, cultivation in the subtropical areas of China, including the Guangdong Province, the Guangxi Zhuang Autonomous Region, the Hainan Province, and Northern China, is adversely and frequently affected by cold stress (i.e., water temperature below 16 • C). Low water temperature causes growth retardation, digestion malfunction, and energy metabolism disorders in the Pacific white shrimp. Breeding of new shrimp cultivars with cold tolerance is, therefore, urgently required (Peng et al., 2015). However, the conventional generation cycle of the Pacific white shrimp is long and the breeding efficiency is extremely low (Castillo-Juarez et al., 2015). Although cold resistance in organisms is usually controlled by quantitative traits associated with multiple physiological and biochemical processes, identifying major genes conferring cold resistance is very important. Zinc finger-containing glycine-rich RNAbinding proteins confer cold tolerance in rice (Kim and Kang, 2006;Chaikam and Karlson, 2008;Park et al., 2009). Recently, transcriptomic, proteomic, and metabolomic studies have been performed to effectively investigate stress tolerance in animals and plants (Chen et al., 2008;Kowalczuk et al., 2018). For example, in the Pacific white shrimp, high-throughput sequencing was used to identify the miRNAs that respond to white spot syndrome virus (WSSV) infections  and in organisms under acute ammonia stress (Lu et al., 2016). To reveal the molecular mechanisms of cold tolerance in Pacific white shrimp, a cold-tolerant cultivar GH2 and a cold-sensitive cultivar GH1 were investigated under normal temperature, low temperature, and in recovery stages using proteomics, bioinformatics, and qRT-PCR techniques. The results provide new insight into the molecular mechanisms of the Pacific white shrimp associated with cold stress to facilitate the breeding of new cultivars with increased cold tolerance.

Pacific White Shrimp Cultivars and Treatment
Two Pacific white shrimp cultivars, Guihai2 (GH2) with cold tolerance and Guihai1 (GH1) with high yields but sensitive to cold stress, were reared at Fangchenggang Aquaculture Base of the Guangxi Academy of Fishery Sciences (E108 • 41 62 , N21 • 61 30 ). GH1 is a high-yield cultivar developed by the Guangxi Academy of Fishery Sciences in 2013. Since 2015, the acclimation conditions for cold-tolerant cultivar breeding have included exposure to 16 • C for 144 h, and these breeding conditions for the cold-resistant cultivar GH2 result in its resistance to low temperatures. GH2 shrimp grew well at 20-22 • C and the low-temperature resistance improved by 3-5 • C, and its adaptability to variable temperatures enhanced significantly. Organisms were adapted to conditions of 28 • C, 32-35% salinity, and a dissolved oxygen level of ≥5 mg L −1 for 1 week prior to experimentation. Organisms used were 3 months old with a weight of 10-15 g. Shrimp were fed three times daily, with a daily ration of approximately 5% of the total body weight of a shrimp. Two groups of shrimp were cooled to specific experimental temperatures and then maintained at those temperatures for sampling at different points. The control group was maintained at 28 • C. To acclimate the experimental groups to a lower temperature, the temperature was lowered to 16 • C, which was maintained for 144 h, with feeding halted only in the experimental groups from this point forward. After 144 h, one group was allowed to recover to 28 • C ( Table 1). Three organisms from each group were sampled, with shrimp sacrificed by destroying the main center nerve and rapidly isolating the hepatopancreas for subsequent protein extraction (He et al., 2018). Each of the three group treatments were repeated twice. Reproducibility analysis of the two repeated trials was conducted by Pearson correlation coefficient methods.

Proteomic Analysis
Shrimp samples were sonicated three times on ice using a highintensity ultrasonic processor (Ningbo Scientz Biotechnology Co., Ltd., Ningbo, China) in lysis buffer [8 M urea, 2 mM (Ethylenedinitrilo)tetraacetic acid (EDTA), 10 mM DL-Dithiothreitol (DTT), and 1% Protease inhibitor cocktail]. The remaining debris was removed by centrifugation at 20,000 g at 4 • C for 10 min. Finally, protein was precipitated with cold 15% TCA for 2 h at −20 • C. After centrifugation at 4 • C for 10 min, the supernatant was discarded, and the remaining precipitate was washed three times with cold acetone. Protein precipitate was redissolved in buffer [8 M urea, 100 mM Triethylammonium bicarbonate (TEAB), pH 8.0], and the protein concentration was determined by using a 2-D Quant Kit (GE Healthcare, Beijing, China) according to the manufacturer's instructions. For trypsin digestion, protein samples were diluted by adding 100 mM TEAB to a urea concentration less than 2 M. Finally, trypsin was added at a 1:50 trypsin-to-protein mass ratio for the first digestion overnight and a 1:100 trypsin-to-protein mass ratio for a second 4 h-digestion. After trypsin digestion, peptides were desalted by using a Strata X C18 SPE column (Phenomenex, Tianjin, China) and vacuum-dried for tandem mass tag (TMT) labeling. Samples were mixed and then fractionated by using high-pH reverse-phase high-performance liquid chromatography (HPLC) using an Agilent 300 Extend C18 column (5 µm particles, 4.6 mm ID, 250 mm length). Peptides were dissolved in 0.1% FA and directly loaded onto a reversed-phase precolumn (Acclaim PepMap 100, Thermo Fisher Scientific, Shanghai, China). Peptide separation was performed using a reversed-phase analytical column (Acclaim PepMap RSLC, Thermo Fisher Scientific, Shanghai, China). The gradient comprised an increase from 6 to 23% solvent B (0.1% FA in 98% ACN) over 26 min, 23 to 35% in 8 min and climbing to 80% in 3 min, and then holding at 80% for the last 3 min, all at a constant flow rate of 400 nL min −1 on an EASY-nLC 1000 ultra performance liquid chromatography (UPLC) system. The peptides were subjected to an NSI source followed by tandem mass spectrometry (MS/MS) in Q Exactive TM plus (Thermo Fisher Scientific, Shanghai, China) coupled online to the UPLC.

Protein Annotation and Subcellular Localization
Gene ontology (GO) annotations were created by searching the UniProt-GOA database 1 . The functional 1 http://www.ebi.ac.uk/GOA/ description of protein domains were annotated by using the InterProScan tool 2 based on protein sequence alignment; the InterPro domain database was also used. Subcellular localization was performed using the WoLF PSORT online software 3 .

GO Enrichment Analysis
Proteins were classified by GO annotation into three categories: biological process, cellular compartment, and molecular function. For each category, a two-tailed Fisher's exact test was employed to test the enrichment of the differentially expressed protein against all the identified proteins. Correction for multiple hypothesis testing was carried out using standard false discovery rate (FDR) control methods. A GO with a corrected p-value < 0.05 was considered to be significant.

Pathway Enrichment Analysis
The Kyoto Encyclopedia of Genes and Genomes (KEGG) database 4 was used to identify the enriched pathways by a two-tailed Fisher's exact test to investigate the enrichment of the differentially expressed protein against all identified proteins. Correction for multiple hypothesis testing was carried out using standard false discovery rate control methods. The pathway with a corrected p-value < 0.05 was considered significant. These pathways were classified

Protein Domain Enrichment Analysis
For each category of protein, InterPro, a database resource that provides a functional analysis of protein sequences by classifying them into families and predicting the presence of domains and important sites 5 was searched and a two-tailed Fisher's exact test was employed to test the enrichment of the differentially expressed protein against all identified proteins. Correction for multiple hypothesis testing was carried out using standard false discovery rate control methods and domains with a corrected p-value < 0.05 were considered significant. 5 https://www.ebi.ac.uk/interpro/

Comparative Clustering Enrichment Analysis
All the protein groups obtained after enrichment were collated, along with their p-values, and then filtered for categories that were at least enriched in one of the clusters with p-value < 0.05. This filtered p-value matrix was transformed by the function x = −log10 (p-value). Finally, these x values were z-transformed for each category. The z scores were then clustered by using one-way hierarchical clustering (Euclidean distance, average linkage clustering) in Genesis Software. Cluster membership was visualized by a heat map using the "heatmap.2" function from the "gplots" R-package. All the data related to this study are available on iProX 6 with id IPX0001223000/PXD009889.

RNA Extraction and cDNA Synthesis
Total RNA was extracted from the two shrimp cultivars (18 samples) using the RP5611 RNA Rapid Extraction Kit (BioTeke Corporation, Beijing, China). The quality and concentration of the extracted RNA were determined by agarose gel electrophoresis and measured by using a spectrophotometer. First-strand cDNA was synthesized from 2 mg of the total RNA with MMLV reverse transcriptase and random hexamer primer (Takara, Dalian, China) according to the manufacturer's instructions.

Quantitative RT-PCR (qRT-PCR) Assays
Transcriptomic analysis of the same samples (data not published) was performed in addition to protein identification. The ID numbers of proteins were the same as that of the transcriptome sequence. After selecting the protein to be verified according to the results of the differential analysis, the transcriptome sequence with the corresponding ID number was retrieved for primer design and verification. Based on proteomic analysis, 18 differentially expressed proteins under cold stress were selected for designing primers and qRT-PCR verifications in the two shrimp cultivars. The primers of the selected genes were designed by using Primer Premier 6 (PREMIER Biosoft, Palo Alto, CA, United States) and synthesized by GENEWIZ Biotechnology (Suzhou, China). The qRT-PCR assays were performed using 2× SYBR Green qPCR ProMix (Takara, Dalian, China) on an ABI 7500 Real-Time PCR System (ABI, Thermo Fisher Scientific, Shanghai, China) following the manufacturer's instructions. Each plate was analyzed independently in triplicate for all the reference and selected genes. The beta-actin gene was used as a reference gene, and the 2 − Ct method (Livak and Schmittgen, 2001) was used to evaluate the relative gene expression levels. The protein accession numbers, description, and primers used in the qRT-PCR tests are listed in Table 2.

RESULTS
We first analyzed the proteome quality of two repeats (Figure 2). The similarity of the two repeats ranged from 0.879 in GH1-R to 0.918 in GH2-28 and GH1-16 based on Pearson correlation coefficient analysis, and these values suggested that the proteome quality was suitable for subsequent analysis.

Proteomic Analysis Revealed Differentially Expressed Proteins Between Cold-Tolerant Cultivar GH2 and Cold-Sensitive Cultivar GH1 Under Cold Stress
In total, 3,349 identified proteins were found in the proteome, among which 2,736 proteins were quantified (Table 3 and Supplementary Tables S1, S2). The number of differential proteins between the two shrimp cultivars under the same temperature were much lower than those under different temperatures. At the normal temperature of 28 • C (GH2-28 vs. GH1-28), there were 71 upregulated and 59 downregulated proteins differentially expressed between the cold-tolerant cultivar GH2 and the cold-sensitive cultivar GH1 (Supplementary Table S3). However, under the low temperature of 16 • C (GH2-16 vs. GH1-16), the number of differentially expressed proteins increased to 274 upregulated and 139 downregulated proteins. With the recovery regime of 16 to 28 • C (GH2-R vs. GH2-R), there were 236 upregulated and 129 downregulated proteins. These results indicated that there were different proteomic profiles in cold-tolerant and cold-sensitive cultivars, especially under the low-temperature treatments and in the recovery period.  The GO Distribution of Differentially Expressed Proteins To further understand the functions and features of the identified and quantified proteins, they were classified into four categories, namely, gene ontology, domain, pathway, and subcellular localization. According to the GO annotation information of the identified proteins, the amount of the differentially expressed proteins for each GO term of level 2 was aggregated ( Table 4 and   Supplementary Table S4). A low temperature caused important biological process changes in the two cultivars. First, the number of proteins for metabolic processes, cellular processes, single-organism processes, localization, and biological regulation significantly increased. Second, the low temperature caused an increase in the response of proteins to stimuli. Third, the numbers of proteins for cellular component organization or biogenesis increased only in the recovery group. There were 130 proteins  differentially expressed in the cold-tolerant cultivar GH2 and the cold-sensitive cultivar GH1 under normal conditions. A low temperature caused an increase in the number of proteins attributed to biological processes, especially the response of proteins to stimuli in the cold-tolerant cultivar GH2. During the recovery process, the growth of shrimp recovered, and a decrease in the number of proteins that responded to stimuli and an increase in the number of proteins attributed to cellular component organization or biogenesis was observed. In addition, 22 and 11 proteins, respectively, of an unknown function, were present in the low-temperature and recovery groups. Significant changes in the upregulated and downregulated proteins were also found according to cellular component and molecular function categories. Not only did the number of proteins associated with the cell, organelle, macromolecular complexes, and the membrane extracellular region categories increase, but the proteins associated with the extracellular regions also increased in the cold-tolerant cultivar GH2 under the low temperature.

Classification of the Identified Proteins Based on Subcellular Location
The amount of differentially expressed proteins in each subcellular location was determined according to the subcellular location annotation of the identified proteins ( Table 5 and Supplementary Table S5). Low-temperature treatments induced the increase of proteins associated with the plasma membrane, cytosol, mitochondria, nuclear cytosol, as well as extracellular and nuclear proteins. Interestingly, new proteins were found in the peroxisome, plasma membrane, and mitochondria, as well as new extracellular and nuclear proteins in both GH2 and GH1 under cold stress conditions.

Functional Enrichment of Differentially Quantified Proteins
After the proteins were assigned to different categories, the quantities were calculated via the −log10 (p-value) method (Figure 3 and Supplementary Table S6). Under normal temperature (GH2-28 vs. GH1-28), the difference in protein content varied from 1.37 for cytoplasm to 4.75 for hydrolase activity. After cold treatments (GH2-16 vs. GH1-16), the difference in protein content varied from 1.56 for monovalent inorganic cation transmembrane transporter activity to 4.85 for small-molecule metabolic processes. During the recovery phase (GH2-R vs. GH1-R), the difference in protein content varied from 1.4 for pyridoxal phosphate binding to 6.03 for intracellular proteins. These results indicate that during cold treatment and the recovery process, the number of proteins and protein contents increased both in GH1 and GH2.

Functional Enrichment-Based Clustering for Comparable Groups
After a GO-based enrichment analysis of all the proteins, KEGG pathway enrichment-based clustering analysis was employed to compare all the changes among GH1 and GH2 under the different treatments (Figure 4 and Supplementary  Table S3). Significant changes were observed in oxidative phosphorylation, glycine, serine, and threonine metabolism, and in cardiac muscle contraction under cold stress (GH2-16 vs. GH1-16). In GH2, a low temperature caused changes in proteins associated with hypertrophic cardiomyopathy (HCM), extracellular matrix (ECM)-receptor interaction, toxoplasmosis, and antigen processing (GH2-28 vs. GH2-16). On the contrary, in GH1, a low temperature caused protein changes in the lysosome, with other glycan degradation proteins, in glycosphingolipid biosynthesis, and in glutathione metabolism (GH1-28 vs. GH1-16). Protein domains were analyzed to further explore specific protein families (Figure 5 and Supplementary Table S3). Significant changes were observed in proteins with 2Fe-2S ferredoxin-type iron-sulfur-binding domains, with beta-grasp domains, in aldehyde oxidase/xanthine dehydrogenases, and with molybdopterin binding under cold stress (GH2-16 vs. sGH1-16). In GH2, a low temperature caused changes in proteins with leucine-rich repeat (LRR) domains, L domain-like is one kind of protein with thioredoxin domains, disulphide isomerases, and laminin EGF domains (GH2-28 vs. GH2-16). On the contrary, in GH1, a low temperature caused changes in glutathione S-transferases (GSTs), glycoside hydrolase superfamily, PA domain, and C-type lectin (GH1-28 vs. GH1-16). These results corresponded with that of biological processes and cellular components.

Validation of Proteomic Data by qRT-PCR Analysis
To validate the proteomic data, we performed the qRT-PCR analysis of 18 genes belonging to four groups: enzymes (i.e., tetraspanin CD63, m7GpppX pyrophosphatase, and zinc proteinase), transcription factors (i.e., TF, AP2, and zinc-finger), lectin family proteins, and DAPPUDRAFT family proteins. The mRNA expression of 18 genes corresponded to 18 proteins that were identified as differentially expressed in GH1 and GH2 under cold stress (Figure 6 and Supplementary Table  S7). The proteomic and qRT-PCR data exhibited the same trends in most proteins such as lectin 1 (CL9159Contig1) and zinc-finger in Ran-binding protein and other domaincontaining proteins (CL23275Contig1). These data indicate that the proteomic data were reliable and could be used for future studies. Furthermore, we summarized the domains of differentially expressed proteins in all the nine compared groups ( Table 6). There were 18 proteins with significant differences between GH2-28 vs. GH1-28. However, under cold stress and under the recovery phase, differentially expressed proteins with significant differences included lectin, those with 2Fe-2S ferredoxin-type iron-sulfur-binding domains, alcohol dehydrogenase, GST, and LRR domains, among others.

DISCUSSION
Pacific white shrimp are native to tropical areas, but they are now widely cultivated in subtropical areas. Pacific white shrimp can be raised for 2-3 cycles each year in southern China, whereas only one growth cycle can occur in northern China due to temperature limitations. Therefore, to extend the breeding cycle of these shrimp in northern China, it is necessary to breed coldresistant cultivars to reduce death rates caused by cold stress. Suppression-subtractive hybridization (Peng et al., 2016) and transcriptomic studies (Long et al., 2013;Chopra et al., 2015;Wang et al., 2015) of shrimp and zebrafish have been carried out to better understand the molecular mechanisms involved in cold tolerance in aquatic organisms. Recently, proteomic techniques were extensively used in research associated with human disease (Zhu et al., 2018), protein o-glycosylation (You et al., 2018), and apicomplexan biology (Yakubu et al., 2018). Little is known about the proteome of Pacific white shrimp (Lu et al., 2016), especially under cold stress. Here, we carried out a comprehensive study of variations in proteins in a cold-tolerant cultivar (GH2) and a cold-sensitive cultivar (GH1) under low-temperature treatments. The hepatopancreas performs some of the same functions that the pancreas and the liver perform in humans; therefore, it is often used as an indicator of organismal health and for the nutritional, metabolic, and disease status in shrimp (Shekhar et al., 2013;Zhang et al., 2014). The hepatopancreas has also been used in detecting cold responsive genes and proteins Peng et al., 2016). Here, the results demonstrated significant differences between the expressed proteins in the cold-tolerant and cold-sensitive cultivars, and these protein functions were largely associated with metabolic processes, cellular processes, single-organism processes, localization, biological regulation, and in response to stimuli. Heat shock proteins have previously been found in Drosophila melanogaster (Burton et al., 1988;Fujikake et al., 2005;Colinet et al., 2010), quail spleen (Ren et al., 2018), and chicken hearts , under cold stress conditions. However, proteomics data reported herein did not indicate that heat shock proteins were involved in the cold response in the hepatopancreas of white shrimp. These may be due to the functions of the heat shock proteins in acute cold response (Leandro et al., 2004). Here, the coldtolerant cultivar GH2 was reared for cold adaption. Based on proteomics and qRT-PCR results, we speculate that the following are the main mechanisms of cold tolerance in Pacific white shrimp.
in white shrimp . Overall, the data suggest that more stress response and detoxing proteins and enzymes were synthesized under cold treatments in the cold-tolerant cultivar GH2.

Transcriptional Regulation Was Involved in Cold Stress Response
Transcription factors were highly accumulated in the coldtolerant cultivar GH2. LRR proteins are known to be involved in cold stress response in plants (Meyer et al., 1999;Yang et al., 2014). In humans, LRR proteins are considered to be related to Parkinson's disease (Kett and Dauer, 2012) and lipid rafts (Hatano et al., 2007). Our results suggested that LRR transcription factors may play roles in lipid membrane protection resulting in cold tolerance in GH2. We found that AP2 and zinc-finger transcription factors were highly expressed in GH2-28 and GH2-16 compared to those of GH1-28 and GH1-16. AP2 interaction with (G/a) (C/t)CGAC motif (Xue, 2002) has been confirmed to participate in cold tolerance in plants (Kang et al., 2011;Du et al., 2016). Zinc-finger transcription factors have been shown to be involved in the response to cold acclimation in catfish (Ju et al., 2002), bees (Xu et al., 2017), and rice (Bai et al., 2015). Cold-responsive genes, such as heat shock proteins (HSPs) and serine/threonine kinases (STKs) and especially rely on the lectin type, were significantly upregulated through transcriptional regulation (Xu et al., 2017).

Lectin and DAPPUDRAFT (Alkaline Phosphatase) Protein Families Played Versatile Roles in Cold Tolerance
On the structural level, the lectin protein family is tightly related to cold tolerance (Moulessehoul et al., 1992;Gronwald , 1998). Interestingly, the function of lectin proteins relies on alpha1,3-galactosyltransferase (Khraltsova et al., 2000), cytoskeleton (Timofeeva et al., 2000), and especially rely on the lectin type (Meng et al., 2017;Hung et al., 2018). In this study, three C-type lectin proteins were differentially expressed under the different temperature treatments (Figure 5), suggesting that each individual lectin played a specific role in cold tolerance in this organism. Interestingly, the protein expression levels observed from the proteomic analysis were similar to those observed in the qRT-PCR analysis (Figure 6). DAPPUDRAFT has been identified as an alkaline phosphatase (EC3.1.3.1) protein family in the water flea (Colbourne et al., 2011). Alkaline phosphatases are key enzymes in sea bream fish that are involved in cold tolerance (Mateus et al., 2017). These enzymes have also shown to have protective roles in liver (Bruinsma et al., 2015;Hjorleifsson and Asgeirsson, 2016), bone (Mateus et al., 2017), and ischemia-reperfusion injuries in rats (Li R. et al., 2014). Our analyses found that the expression levels of several DAPPUDRAFT proteins (240262, 307838, and 206907) were higher in GH2 than in GH1, suggesting that the cold tolerance of GH2 was related to alkaline phosphatase activity, most likely due to the fact that these enzymes can help protect the liver, bones, and blood of shrimp from cold-stress damage.
Interestingly, not all the cold-responsive gene expressions that procede protein synthesis, and this could be due to a lack of energy for modification at post-transcriptional and posttranslational levels, such as glycosylation (Yakubu et al., 2018;You et al., 2018). It is well known that qRT-PCR is an effective way to perform quantitative proteomic analysis (Yan et al., 2006;Fan et al., 2013;Zhang et al., 2016), and there are five types of expression patterns at the qRT-PCR and protein level, namely: (A) upregulated at both the transcriptional and protein level, (B) upregulated at the transcriptional level but downregulated at the protein level, (C) downregulated at the transcriptional level but upregulated at the protein level, (D) no change at the transcriptional level but upregulated at the protein level, and (E) no change at the transcriptional level but downregulated at the protein level (Zhang et al., 2016). Expressions of most of the lectin genes exhibited the same patterns as the protein expression ( Table 6 and Figure 6), suggesting that the key proteins involved in cold-tolerance in GH2 were lectin and phosphatase, among others. Based on the findings from previous studies, it can be speculated that cold stress first leads to the damage of the plasma membrane (Takahashi et al., 2014) that will induce signal cascades by phosphatase (Kristensen et al., 2016), zinc proteinase, and m7GpppX diphosphatase. Through a second messenger (Zieger et al., 2011;Kocharunchitt et al., 2012) or hormone synthesis (Cerny et al., 2014) at the transcriptional level, the AP-2 transcription complex and zinc-finger transcription factors (Koehler et al., 2012) regulate the expression of functional genes, including different lectin and DAPPUDRAFT proteins (Figures 5, 6). After translation and regulation at the posttranslational level (Chen et al., 2014), proteins are synthesized and used for the cold stress response in Pacific white shrimp.

ETHICS STATEMENT
This study was carried out in accordance with the recommendations of the guidelines of the Ministry of Agriculture and Rural Affairs of the China.

AUTHOR CONTRIBUTIONS
J-XP and X-LC analyzed and interpreted the data and drafted the manuscript. J-XP, P-PH, P-YW, BZ, Y-ZZ, and Q-YL set the experimental design and coordinated the study. J-XP, P-PH, and P-YW performed the proteome study. XC, MP, D-GZ, and C-LY assisted with proteomic analysis.

ACKNOWLEDGMENTS
We thank LetPub (www.letpub.com) for their linguistic assistance during the preparation of this manuscript.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys. 2018.01399/full#supplementary-material TABLE S1 | Protein amino sequences in proteomics analysis of GH1 and GH2 under different temperature treatments.