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 RNA-binding 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 (Zeng et al., 2015) 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.

Materials and Methods

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⁻¹ 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.

TABLE 1

TABLE 1. Samples name of proteome analysis.

All the identified proteins were compared among the two cultivars under the different temperature treatments. The experimental procedure is shown in Figure 1. The two shrimp cultivars were treated at three temperature levels (Table 1). Differentially expressed proteins were cross-compared in 12 groups, namely: 1. GH2-28 vs. GH1-28, 2. GH2-28 vs. GH2-16, 3. GH2-R vs. GH2-28, 4. GH2-16 vs. GH1-16, 5. GH2-16 vs. GH2-R, 6. GH2-R vs. GH1-R, 7. GH1-16 vs. GH1-28, 8. GH1-R vs. GH1-28, 9. GH1-16 vs. GH1-R, 10. GH2GH1-16 vs. GH2GH1-28 (hereafter abbreviated as 16 vs. 28), 11. GH2GH1-R vs. GH2GH1-28 (hereafter abbreviated as R vs. 28), and 12. GH2GH1-R vs. GH2GH1-16 (hereafter abbreviated as R vs. 16).

FIGURE 1

FIGURE 1. Schematic diagram of technical route of treatment, proteomic analysis and verification.

Proteomic Analysis

Shrimp samples were sonicated three times on ice using a high-intensity 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⁻¹ 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.

Bioinformatic Analysis

Protein Annotation and Subcellular Localization

Gene ontology (GO) annotations were created by searching the UniProt-GOA database¹. The functional description of protein domains were annotated by using the InterProScan tool² based on protein sequence alignment; the InterPro domain database was also used. Subcellular localization was performed using the WoLF PSORT online software³.

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⁴ 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 into hierarchical categories according to the KEGG website.

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⁵ 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.

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⁶ 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.

TABLE 2

TABLE 2. Protein accession, description and qRT-PCR primer pairs.

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.

FIGURE 2

FIGURE 2. Reproducibility analysis of two repeated trials by Pearson correlation coefficient. For each replicate, relative expressions were calculated and Pearson’s correlation coefficient was calculated.

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.

TABLE 3

TABLE 3. Numbers of differentially expressed proteins to be quantified between different groups.

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.

TABLE 4

TABLE 4. The GO distribution of all up-regulated and down-regulated proteins.

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.

TABLE 5

TABLE 5. The subcellular location of different proteins of nine groups.

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.

FIGURE 3

FIGURE 3. GO-based enrichment analysis of all the proteins.

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).

FIGURE 4

FIGURE 4. KEGG pathway enrichment-based clustering analysis of all the identified proteins.

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.

FIGURE 5

FIGURE 5. Protein domain of proteins identified in each treatment group.

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 domain-containing 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.

FIGURE 6

FIGURE 6. Validation of proteomic data by qRT-PCR analysis. The protein accession number and full names were: SCL8Contig75, tetraspanin-like protein CD63; CL3798Contig3, zinc proteinase; GH1_16_2-c14615_g2_i2, m7GpppX diphosphatase; SCL1Contig301, AP-2 complex subunit alpha; CL23275Contig1, Zn-finger in Ran-binding protein and other domain-containing proteins; CL9159Contig1, lectin 1; SCL12Contig31, lectin D1; CL638Contig3, C-type lectin; CL4347Contig1, lectin 2; CL823Contig3, lectin D2; CL10443Contig1, hypothetical protein DAPPUDRAFT_240263; SCL1Contig877, hypothetical protein DAPPUDRAFT_307838; CL6694Contig3, hypothetical protein BRAFLDRAFT_206907; GH1_28_1-c16133_g2_i1, hypothetical protein DAPPUDRAFT_321849; CL9739Contig1, hypothetical protein DAPPUDRAFT_54086; CL10847Contig1, hypothetical protein DAPPUDRAFT_215063; CL5480Contig1, hypothetical protein DAPPUDRAFT_31637, partial; CL4762Contig4, hypothetical protein DAPPUDRAFT_303198.

TABLE 6

TABLE 6. The subcellular location of different proteins of nine groups.

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 cold-resistant 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 (Fan et al., 2016; 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 (Zhao et al., 2013), 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 cold-tolerant 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.

Functional Enzymes Participate in Cold Tolerance

In Panagrolaimus davidi, cold stress can induce the expression of trehalose-6-phosphate phosphatase, superoxide dismutase, and glutathione peroxidase (Thorne et al., 2014). Similar results were observed in Camellia sinensis during cold acclimation (Wang et al., 2013) and in Antarctic notothenioid fish (Chen et al., 2008). Typically, a cold signal will first lead to damages to the plasma membrane. Next, mitogen-activated protein kinases (MAPKs) and serine/threonine kinases are phosphorylated and enter the nucleus where various transcription factors, enzymes, and other proteins that modulate cellular activities are phosphorylated (Fujiwara and Denlinger, 2007; Wang et al., 2017). Our results indicate that free radical stress scavenging enzymes played key roles in white shrimp under cold stress in these experiments, consistent with studies in other species. For example, cold treatment has been shown to cause significant changes in GST (Tierbach et al., 2018), LRR domains, tetraspanin CD63 (Li G. et al., 2014; Chaiyadet et al., 2017), m7GpppX pyrophosphatase (Malys et al., 2004), zinc proteinase (Loewy et al., 1993), and different types of lectin protein families (Figures 5>, 6). GSTs play important roles in insecticide/drug resistance and stress response in insects (Sookrung et al., 2018), rats (Wang et al., 2018; Yang et al., 2018), and the plant species Arabidopsis thaliana (Banday and Nandi, 2018). In addition, similar results were found in 2D-electrophoresis proteomics studies, i.e., GST Mu 3-like was upregulated more than 1.5-fold under cold stress in white shrimp (Fan et al., 2016). 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 cold-tolerant 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 et al., 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 post-translational 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 post-translational 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.

Funding

This study was supported by the Natural Science Foundation of China (Grant No. 31460690), the Special Funds for Ba Gui Scholars Project Construction (Grant No. BGXZ-NMBDX-05), the China Agriculture Research System (Grant No. CARS-48), and the Fundamental Research Funds for Non-profit Research Institutes under Guangxi Zhuang Autonomous Region (Grant No. GXIF-2014-002).

Conflict of Interest Statement

The 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.

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.

TABLE S2 | Protein accession and description of 3,349 identified proteins in proteomics analysis of GH1 and GH2 under different temperature treatments.

TABLE S3 | Differentially expressed up-regulated and down-regulated proteins of GH1 and GH2 under different temperature treatments.

TABLE S4 | GO annotation information of identified proteins of GH1 and GH2 under different temperature treatments.

TABLE S5 | Subcellular location of differentially expressed proteins of GH1 and GH2 under different temperature treatments.

TABLE S6 | Proteins GO Catalogs and quantities of GH1 and GH2 under different temperature treatments.

TABLE S7 | 18 genes qRT-PCR data and corresponding 18 proteins quantities.

Footnotes

1. ^ http://www.ebi.ac.uk/GOA/
2. ^ http://www.ebi.ac.uk/interpro/search/sequence-search
3. ^ https://wolfpsort.hgc.jp/
4. ^ https://www.genome.jp/kegg/
5. ^ https://www.ebi.ac.uk/interpro/
6. ^ http://www.iprox.org/

References

Bai, B., Wu, J., Sheng, W. T., Zhou, B., Zhou, L. J., Zhuang, W., et al. (2015). Comparative analysis of anther transcriptome profiles of two different rice male sterile lines genotypes under cold stress. Int. J. Mol. Sci. 16, 11398–11416. doi: 10.3390/ijms160511398

PubMed Abstract | CrossRef Full Text | Google Scholar

Banday, Z. Z., and Nandi, A. K. (2018). Arabidopsis thaliana GLUTATHIONE-S-TRANSFERASE THETA 2 interacts with RSI1/FLD to activate systemic acquired resistance. Mol. Plant Pathol. 19, 464–475. doi: 10.1111/mpp.12538

PubMed Abstract | CrossRef Full Text | Google Scholar

Bruinsma, B. G., Wu, W., Ozer, S., Farmer, A., Markmann, J. F., Yeh, H., et al. (2015). Warm ischemic injury is reflected in the release of injury markers during cold preservation of the human liver. PLoS One 10:e0123421. doi: 10.1371/journal.pone.0123421

PubMed Abstract | CrossRef Full Text | Google Scholar

Burton, V., Mitchell, H. K., Young, P., and Petersen, N. S. (1988). Heat shock protection against cold stress of Drosophila melanogaster. Mol. Cell. Biol. 8, 3550–3552. doi: 10.1128/MCB.8.8.3550

PubMed Abstract | CrossRef Full Text | Google Scholar

Castillo-Juarez, H., Campos-Montes, G. R., Caballero-Zamora, A., and Montaldo, H. H. (2015). Genetic improvement of Pacific white shrimp [Penaeus (Litopenaeus) vannamei]: perspectives for genomic selection. Front. Genet. 6:93. doi: 10.3389/fgene.2015.00093

PubMed Abstract | CrossRef Full Text | Google Scholar

Cerny, M., Jedelsky, P. L., Novak, J., Schlosser, A., and Brzobohaty, B. (2014). Cytokinin modulates proteomic, transcriptomic and growth responses to temperature shocks in Arabidopsis. Plant Cell Environ. 37, 1641–1655. doi: 10.1111/pce.12270

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaikam, V., and Karlson, D. (2008). Functional characterization of two cold shock domain proteins from Oryza sativa. Plant Cell Environ. 31, 995–1006. doi: 10.1111/j.1365-3040.2008.01811.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaiyadet, S., Krueajampa, W., Hipkaeo, W., Plosan, Y., Piratae, S., Sotillo, J., et al. (2017). Suppression of mRNAs encoding CD63 family tetraspanins from the carcinogenic liver fluke Opisthorchis viverrini results in distinct tegument phenotypes. Sci. Rep. 7:14342. doi: 10.1038/s41598-017-13527-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Gao, X., Wang, B., Chen, F., Wu, N., and Zhang, Y. (2014). Proteomic approach to reveal the proteins associated with encystment of the ciliate Euplotes encysticus. PLoS One 9:e97362. doi: 10.1371/journal.pone.0097362

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Cheng, C. H., Zhang, J., Cao, L., Chen, L., Zhou, L., et al. (2008). Transcriptomic and genomic evolution under constant cold in Antarctic notothenioid fish. Proc. Natl. Acad. Sci. U.S.A. 105, 12944–12949. doi: 10.1073/pnas.0802432105

PubMed Abstract | CrossRef Full Text | Google Scholar

Chopra, R., Burow, G., Hayes, C., Emendack, Y., Xin, Z., and Burke, J. (2015). Transcriptome profiling and validation of gene based single nucleotide polymorphisms (SNPs) in sorghum genotypes with contrasting responses to cold stress. BMC Genomics 16:1040. doi: 10.1186/s12864-015-2268-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Chuchird, N., Rorkwiree, P., and Rairat, T. (2015). Effect of dietary formic acid and astaxanthin on the survival and growth of Pacific white shrimp (Litopenaeus vannamei) and their resistance to Vibrio parahaemolyticus. Springerplus 4:440. doi: 10.1186/s40064-015-1234-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Colbourne, J. K., Pfrender, M. E., Gilbert, D., Thomas, W. K., Tucker, A., and Oakley, T. H. et al. (2011). The ecoresponsive genome of Daphnia pulex. Science 331, 555–561. doi: 10.1126/science.1197761

PubMed Abstract | CrossRef Full Text | Google Scholar

Colinet, H., Lee, S. F., and Hoffmann, A. (2010). Temporal expression of heat shock genes during cold stress and recovery from chill coma in adult Drosophila melanogaster. FEBS J. 277, 174–185. doi: 10.1111/j.1742-4658.2009.07470.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, C., Hu, K., Xian, S., Liu, C., Fan, J., Tu, J., et al. (2016). Dynamic transcriptome analysis reveals AP2/ERF transcription factors responsible for cold stress in rapeseed (Brassica napus L.). Mol. Genet. Genomics 291, 1053–1067. doi: 10.1007/s00438-015-1161-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, L., Wang, A., Miao, Y., Liao, S., Ye, C., and Lin, Q. (2016). Comparative proteomic identification of the hepatopancreas response to cold stress in white shrimp, Litopenaeus vannamei. Aquaculture 454, 27–34. doi: 10.1016/j.aquaculture.2015.10.016

CrossRef Full Text | Google Scholar

Fan, L., Wang, A., and Wu, Y. (2013). Comparative proteomic identification of the hemocyte response to cold stress in white shrimp, Litopenaeus vannamei. J. Proteomics 80, 196–206. doi: 10.1016/j.jprot.2012.12.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujikake, N., Nagai, Y., Popiel, H. A., Kano, H., Yamaguchi, M., and Toda, T. (2005). Alternative splicing regulates the transcriptional activity of Drosophila heat shock transcription factor in response to heat/cold stress. FEBS Lett. 579, 3842–3848. doi: 10.1016/j.febslet.2005.05.074

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujiwara, Y., and Denlinger, D. L. (2007). p38 MAPK is a likely component of the signal transduction pathway triggering rapid cold hardening in the flesh fly Sarcophaga crassipalpis. J. Exp. Biol. 210, 3295–3300. doi: 10.1242/jeb.006536

PubMed Abstract | CrossRef Full Text | Google Scholar

Gronwald, W., Loewen, M. C., Lix, B., Daugulis, A. J., Sonnichsen, F. D., Davies, P. L., et al. (1998). The solution structure of type II antifreeze protein reveals a new member of the lectin family. Biochemistry 37, 4712–4721. doi: 10.1021/bi972788c

PubMed Abstract | CrossRef Full Text | Google Scholar

Hatano, T., Kubo, S., Imai, S., Maeda, M., Ishikawa, K., Mizuno, Y., et al. (2007). Leucine-rich repeat kinase 2 associates with lipid rafts. Hum. Mol. Genet. 16, 678–690. doi: 10.1093/hmg/ddm013

PubMed Abstract | CrossRef Full Text | Google Scholar

He, P., Wei, P., Zhang, B., Zhao, Y., Li, Q., Chen, X., et al. (2018). Identification of microRNAs involved in cold adaptation of Litopenaeus vannamei by high-throughput sequencing. Gene 677, 24–31. doi: 10.1016/j.gene.2018.07.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Heisterkamp, I. M., Schramm, A., De Beer, D., and Stief, P. (2016). Direct nitrous oxide emission from the aquacultured pacific white shrimp (Litopenaeus vannamei). Appl. Environ. Microbiol. 82, 4028–4034. doi: 10.1128/AEM.00396-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Hjorleifsson, J. G., and Asgeirsson, B. (2016). Cold-active alkaline phosphatase is irreversibly transformed into an inactive dimer by low urea concentrations. Biochim. Biophys. Acta 1864, 755–765. doi: 10.1016/j.bbapap.2016.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Hung, L. D., Ly, B. M., Hao, V. T., Trung, D. T., Trang, V. T. D., Trinh, P. T. H., et al. (2018). Purification, characterization and biological effect of lectin from the marine sponge Stylissa flexibilis (Levi, 1961). Comp. Biochem. Physiol. B Biochem. Mol. Biol. 216, 32–38. doi: 10.1016/j.cbpb.2017.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ju, Z., Dunham, R. A., and Liu, Z. (2002). Differential gene expression in the brain of channel catfish (Ictalurus punctatus) in response to cold acclimation. Mol. Genet. Genomics 268, 87–95. doi: 10.1007/s00438-002-0727-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, H. G., Kim, J., Kim, B., Jeong, H., Choi, S. H., Kim, E. K., et al. (2011). Overexpression of FTL1/DDF1, an AP2 transcription factor, enhances tolerance to cold, drought, and heat stresses in Arabidopsis thaliana. Plant Sci. 180, 634–641. doi: 10.1016/j.plantsci.2011.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Kett, L. R., and Dauer, W. T. (2012). Leucine-rich repeat kinase 2 for beginners: six key questions. Cold Spring Harb. Perspect. Med. 2:a009407. doi: 10.1101/cshperspect.a009407

PubMed Abstract | CrossRef Full Text | Google Scholar

Khraltsova, L. S., Sablina, M. A., Melikhova, T. D., Joziasse, D. H., Kaltner, H., Gabius, H. J., et al. (2000). An enzyme-linked lectin assay for alpha1,3-galactosyltransferase. Anal. Biochem. 280, 250–257. doi: 10.1006/abio.2000.4504

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y. O., and Kang, H. (2006). The role of a zinc finger-containing glycine-rich RNA-binding protein during the cold adaptation process in Arabidopsis thaliana. Plant Cell Physiol. 47, 793–798. doi: 10.1093/pcp/pcj047

PubMed Abstract | CrossRef Full Text | Google Scholar

Kocharunchitt, C., King, T., Gobius, K., Bowman, J. P., and Ross, T. (2012). Integrated transcriptomic and proteomic analysis of the physiological response of Escherichia coli O157:H7 Sakai to steady-state conditions of cold and water activity stress. Mol. Cell. Proteomics 11: M111.009019. doi: 10.1074/mcp.M111.009019

PubMed Abstract | CrossRef Full Text | Google Scholar

Koehler, G., Wilson, R. C., Goodpaster, J. V., Sonsteby, A., Lai, X., Witzmann, F. A., et al. (2012). Proteomic study of low-temperature responses in strawberry cultivars (Fragaria x ananassa) that differ in cold tolerance. Plant Physiol. 159, 1787–1805. doi: 10.1104/pp.112.198267

PubMed Abstract | CrossRef Full Text | Google Scholar

Kowalczuk, L., Matet, A., Dor, M., Bararpour, N., Daruich, A., Dirani, A., et al. (2018). Proteome and metabolome of subretinal fluid in central serous chorioretinopathy and rhegmatogenous retinal detachment: a pilot case study. Transl. Vis. Sci. Technol. 7:3. doi: 10.1167/tvst.7.1.3

PubMed Abstract | CrossRef Full Text | Google Scholar

Kristensen, T. N., Kjeldal, H., Schou, M. F., and Nielsen, J. L. (2016). Proteomic data reveal a physiological basis for costs and benefits associated with thermal acclimation. J. Exp. Biol. 219, 969–976. doi: 10.1242/jeb.132696

PubMed Abstract | CrossRef Full Text | Google Scholar

Leandro, N. S., Gonzales, E., Ferro, J. A., Ferro, M. I., Givisiez, P. E., and Macari, M. (2004). Expression of heat shock protein in broiler embryo tissues after acute cold or heat stress. Mol. Reprod. Dev. 67, 172–177. doi: 10.1002/mrd.10397

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., Endsley, M. A., Somasunderam, A., Gbota, S. L., Mbaka, M. I., Murray, J. L., et al. (2014). The dual role of tetraspanin CD63 in HIV-1 replication. Virol. J. 11:23. doi: 10.1186/1743-422X-11-23

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R., Mi, Y., Tan, G., Zhang, W., Li, G., and Sun, X. (2014). A novel in situ model of liver cold ischemia-reperfusion in rats. J. Surg. Res. 192, 195–199. doi: 10.1016/j.jss.2014.05.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2^−ΔΔC_T method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Loewy, A. G., Santer, U. V., Wieczorek, M., Blodgett, J. K., Jones, S. W., and Cheronis, J. C. (1993). Purification and characterization of a novel zinc-proteinase from cultures of Aeromonas hydrophila. J. Biol. Chem. 268, 9071–9078.

PubMed Abstract | Google Scholar

Long, Y., Song, G., Yan, J., He, X., Li, Q., and Cui, Z. (2013). Transcriptomic characterization of cold acclimation in larval zebrafish. BMC Genomics 14:612. doi: 10.1186/1471-2164-14-612

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Malys, N., Carroll, K., Miyan, J., Tollervey, D., and Mccarthy, J. E. (2004). The ‘scavenger’ m7GpppX pyrophosphatase activity of Dcs1 modulates nutrient-induced responses in yeast. Nucleic Acids Res. 32, 3590–3600. doi: 10.1093/nar/gkh687

PubMed Abstract | CrossRef Full Text | Google Scholar

Mateus, A. P., Costa, R., Gisbert, E., Pinto, P. I. S., Andree, K. B., Estevez, A., et al. (2017). Thermal imprinting modifies bone homeostasis in cold-challenged sea bream (Sparus aurata). J. Exp. Biol. 220, 3442–3454. doi: 10.1242/jeb.156174

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, Q., Zhang, J. H., Zhang, H., Zhou, G. L., Ni, R. R., Zhao, Y. N., et al. (2017). Comparative analysis of C-type lectin domain proteins in the ghost moth, Thitarodes xiaojinensis (Lepidoptera: Hepialidae). Insect Sci. doi: 10.1111/1744-7917.12564 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, K., Keil, M., and Naldrett, M. J. (1999). A leucine-rich repeat protein of carrot that exhibits antifreeze activity. FEBS Lett. 447, 171–178. doi: 10.1016/S0014-5793(99)00280-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Moulessehoul, S., Hirsch, M., and Pouliquen, Y. (1992). Lectin-binding sites in quick-frozen, freeze-substituted extracellular matrix of rabbit corneal stroma: a comparative electron microscopic study with different embedding procedures. Exp. Eye Res. 54, 307–312. doi: 10.1016/S0014-4835(05)80220-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Nunan, L. M., Lightner, D. V., Oduori, M. A., and Gasparich, G. E. (2005). Spiroplasma penaei sp. nov., associated with mortalities in Penaeus vannamei, Pacific white shrimp. Int. J. Syst. Evol. Microbiol. 55, 2317–2322. doi: 10.1099/ijs.0.63555-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S. J., Kwak, K. J., Oh, T. R., Kim, Y. O., and Kang, H. (2009). Cold shock domain proteins affect seed germination and growth of Arabidopsis thaliana under abiotic stress conditions. Plant Cell Physiol. 50, 869–878. doi: 10.1093/pcp/pcp037

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, J., Wei, P., Chen, X., Zeng, D., and Chen, X. (2016). Identification of cold responsive genes in Pacific white shrimp (Litopenaeus vannamei) by suppression subtractive hybridization. Gene 575, 667–674. doi: 10.1016/j.gene.2015.09.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, J., Wei, P., Zhang, B., Zhao, Y., Zeng, D., Chen, X., et al. (2015). Gonadal transcriptomic analysis and differentially expressed genes in the testis and ovary of the Pacific white shrimp (Litopenaeus vannamei). BMC Genomics 16:1006. doi: 10.1186/s12864-015-2219-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, J., Liu, C., Zhao, D., and Fu, J. (2018). The role of heat shock protein 70 in oxidant stress and inflammatory injury in quail spleen induced by cold stress. Environ. Sci. Pollut. Res. Int. 25, 21011–21023. doi: 10.1007/s11356-018-2142-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Schnurr, M. E., Yin, Y., and Scott, G. R. (2014). Temperature during embryonic development has persistent effects on metabolic enzymes in the muscle of zebrafish. J. Exp. Biol. 217, 1370–1380. doi: 10.1242/jeb.094037

PubMed Abstract | CrossRef Full Text | Google Scholar

Shekhar, M. S., Kiruthika, J., and Ponniah, A. G. (2013). Identification and expression analysis of differentially expressed genes from shrimp (Penaeus monodon) in response to low salinity stress. Fish Shellfish Immunol. 35, 1957–1968. doi: 10.1016/j.fsi.2013.09.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Sookrung, N., Reamtong, O., Poolphol, R., Indrawattana, N., Seesuay, W., Saelim, N., et al. (2018). Glutathione S-transferase (GST) of American cockroach, periplaneta Americana: classes, isoforms, and allergenicity. Sci. Rep. 8:484. doi: 10.1038/s41598-017-18759-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Takahashi, D., Nakayama, T., Miki, Y., Kawamura, Y., and Uemura, M. (2014). Proteomic approaches to identify cold-regulated plasma membrane proteins. Methods Mol. Biol. 1166, 159–170. doi: 10.1007/978-1-4939-0844-8_13

PubMed Abstract | CrossRef Full Text | Google Scholar

Thorne, M. A., Kagoshima, H., Clark, M. S., Marshall, C. J., and Wharton, D. A. (2014). Molecular analysis of the cold tolerant Antarctic nematode, Panagrolaimus davidi. PLoS One 9:e104526. doi: 10.1371/journal.pone.0104526

PubMed Abstract | CrossRef Full Text | Google Scholar

Tierbach, A., Groh, K. J., Schonenberger, R., Schirmer, K., and Suter, M. J. (2018). Glutathione S-transferase protein expression in different life stages of zebrafish (Danio rerio). Toxicol. Sci. 162, 702–712. doi: 10.1093/toxsci/kfx293

PubMed Abstract | CrossRef Full Text | Google Scholar

Timofeeva, O., Khokhlova, L., Belyaeva, N., Chulkova, Y., and Garaeva, L. (2000). Cytoskeleton-induced alterations of the lectin activity in winter wheat under cold hardening and abscisic acid (ABA). Cell Biol. Int. 24, 375–381. doi: 10.1006/cbir.1999.0520

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, D. Z., Jin, Y. N., Ding, X. H., Wang, W. J., Zhai, S. S., Bai, L. P., et al. (2017). Gene regulation and signal transduction in the ICE-CBF-COR signaling pathway during cold stress in plants. Biochemistry 82, 1103–1117. doi: 10.1134/S0006297917100030

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Zhang, X., and Liu, J. H. (2015). Deep sequencing-based characterization of transcriptome of trifoliate orange (Poncirus trifoliata (L.) Raf.) in response to cold stress. BMC Genomics 16:555. doi: 10.1186/s12864-015-1629-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. C., Zhao, Q. Y., Ma, C. L., Zhang, Z. H., Cao, H. L., Kong, Y. M., et al. (2013). Global transcriptome profiles of Camellia sinensis during cold acclimation. BMC Genomics 14:415. doi: 10.1186/1471-2164-14-415

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Wu, S., Liu, C., Lu, X., and Chen, Z. (2018). Herba Gelsemii elegantis is detoxified by ramulus et folium Mussaenda pubescentis extract by modulating hepatic cytochrome P450 and glutathione S-transferase enzymes in rats. Exp. Ther. Med. 15, 226–234. doi: 10.3892/etm.2017.5351

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, K., Niu, Q., Zhao, H., Du, Y., and Jiang, Y. (2017). Transcriptomic analysis to uncover genes affecting cold resistance in the Chinese honey bee (Apis cerana cerana). PLoS One 12:e0179922. doi: 10.1371/journal.pone.0179922

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, G. P. (2002). An AP2 domain transcription factor HvCBF1 activates expression of cold-responsive genes in barley through interaction with a (G/a)(C/t)CGAC motif. Biochim. Biophys. Acta 1577, 63–72. doi: 10.1016/S0167-4781(02)00410-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Yakubu, R. R., Weiss, L. M., and Silmon De Monerri, N. C. (2018). Post-translational modifications as key regulators of apicomplexan biology: insights from proteome-wide studies. Mol. Microbiol. 107, 1–23. doi: 10.1111/mmi.13867

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, S. P., Zhang, Q. Y., Tang, Z. C., Su, W. A., and Sun, W. N. (2006). Comparative proteomic analysis provides new insights into chilling stress responses in rice. Mol. Cell. Proteomics 5, 484–496. doi: 10.1074/mcp.M500251-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J. F., Liu, Y. R., Huang, C. C., and Ueng, Y. F. (2018). The time-dependent effects of St John’s wort on cytochrome P450, uridine diphosphate-glucuronosyltransferase, glutathione S-transferase, and NAD(P)H-quinone oxidoreductase in mice. J. Food Drug Anal. 26, 422–431. doi: 10.1016/j.jfda.2017.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Wu, K., Gao, P., Liu, X., Li, G., and Wu, Z. (2014). GsLRPK, a novel cold-activated leucine-rich repeat receptor-like protein kinase from Glycine soja, is a positive regulator to cold stress tolerance. Plant Sci. 21, 19–28. doi: 10.1016/j.plantsci.2013.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

You, X., Qin, H., and Ye, M. (2018). Recent advances in methods for the analysis of protein o-glycosylation at proteome level. J. Sep. Sci. 41, 248–261. doi: 10.1002/jssc.201700834

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, D., Chen, X., Xie, D., Zhao, Y., Yang, C., Li, Y., et al. (2013). Transcriptome analysis of Pacific white shrimp (Litopenaeus vannamei) hepatopancreas in response to Taura syndrome Virus (TSV) experimental infection. PLoS One 8:e57515. doi: 10.1371/journal.pone.0057515

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, D. G., Chen, X. L., Xie, D. X., Zhao, Y. Z., Yang, Q., Wang, H., et al. (2015). Identification of highly expressed host microRNAs that respond to white spot syndrome virus infection in the Pacific white shrimp Litopenaeus vannamei (Penaeidae). Genet. Mol. Res. 14, 4818–4828. doi: 10.4238/2015.May.11.14

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W., Zhang, H., Ning, L., Li, B., and Bao, M. (2016). Quantitative proteomic analysis provides novel insights into cold stress responses in petunia seedlings. Front. Plant Sci. 7:136. doi: 10.3389/fpls.2016.00136

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X.-W., Wang, X.-W., Huang, Y., Hui, K.-M., Shi, Y.-R., Wang, W., et al. (2014). Cloning and characterization of two different ficolins from the giant freshwater prawn Macrobrachium rosenbergii. Dev. Comp. Immunol. 44, 359–369. doi: 10.1016/j.dci.2014.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, F. Q., Zhang, Z. W., Wang, C., Zhang, B., Yao, H. D., Li, S., et al. (2013). The role of heat shock proteins in inflammatory injury induced by cold stress in chicken hearts. Cell Stress Chaperones 18, 773–783. doi: 10.1007/s12192-013-0429-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y., Yu, M., Liu, J., Qiao, Y., Wang, L., Li, Z., et al. (2017). Bacterial community associated with healthy and diseased pacific white shrimp (Litopenaeus vannamei) larvae and rearing water across different growth stages. Front. Microbiol. 8:1362. doi: 10.3389/fmicb.2017.01362

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Q., Wu, N., Liu, G., Zhou, Y., Liu, S., Chen, J., et al. (2018). Comparative analysis of serum proteome in congenital scoliosis patients with TBX6 haploinsufficiency - a first report pointing to lipid metabolism. J. Cell Mol. Med. 22, 533–545. doi: 10.1111/jcmm.13341

PubMed Abstract | CrossRef Full Text | Google Scholar

Zieger, M. A., Gupta, M. P., and Wang, M. (2011). Proteomic analysis of endothelial cold-adaptation. BMC Genomics 12:630. doi: 10.1186/1471-2164-12-630

PubMed Abstract | CrossRef Full Text | Google Scholar