- 1Guangxi Key Laboratory of Aquatic Genetic Breeding and Healthy Aquaculture, Guangxi Academy of Fishery Sciences, Nanning, China
- 2Guangxi Shrimp and Crab Breeding Engineering Technology Research Center, Guangxi Academy of Fishery Sciences, Nanning, China
- 3Key Lab of Freshwater Animal Breeding, Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction, Ministry of Education, College of Fishery, Huazhong Agriculture University, Wuhan, China
To characterize the cold tolerance mechanism of the Pacific white shrimp (Litopenaeus vannamei), we performed single-cell RNA sequencing (scRNA-seq) of ∼5185 hepatopancreas cells from cold-tolerant (Lv-T) and common (Lv-C) L. vannamei at preferred and low temperatures (28°C and 10°C, respectively). The cells fell into 10 clusters and 4 cell types: embryonic, resorptive, blister-like, and fibrillar. We identified differentially expressed genes between Lv-T and Lv-C, which were mainly associated with the terms “immune system,” “cytoskeleton,” “antioxidant system,” “digestive enzyme,” and “detoxification,” as well as the pathways “metabolic pathways of oxidative phosphorylation,” “metabolism of xenobiotics by cytochrome P450,” “chemical carcinogenesis,” “drug metabolism-cytochrome P450,” and “fatty acid metabolism.” Reconstruction of fibrillar cell trajectories showed that, under low temperature stress, hepatopancreas cells had two distinct fates, cell fate 1 and cell fate 2. Cell fate 1 was mainly involved in signal transduction and sensory organ development. Cell fate 2 was mainly involved in metabolic processes. This study preliminarily clarifies the molecular mechanisms underlying cold tolerance in L. vannamei, which will be useful for the breeding of shrimp with greater cold tolerance.
Introduction
Temperature may directly affect shrimp metabolism, shelling, growth rate, oxygen consumption, and survival rate; temperature may also indirectly affect shrimp growth via environmental factors such as salinity and microalgae growth (Moullac and Haffner, 2000). Temperature also significantly affects prawn blood cell density (Gomez-Jimenez et al., 2000; Cheng and Chen, 2001), phagocytic activity (Steenbergen et al., 1978), and body fluid phenoloxidase activity. Indeed, cold is one of the main environmental factors affecting the growth and survival of the Pacific white shrimp (Litopenaeus vannamei): suitable growth temperatures for L. vannamei are 25–35°C, and this species poorly tolerates cold exposure. At temperatures below 18°C, L. vannamei stop eating (Ponce-Palafox et al., 1997). Low temperatures also limit L. vannamei breeding season duration and breeding range, affecting the economics of shrimp breeding. In order to breed L. vannamei with a greater range of tolerated temperatures, it is important to study the molecular genetic mechanisms underlying low temperature tolerance in this species. Although some genetic studies of cold tolerance in L. vannamei are available, these studies used high-throughput sequencing to identify microRNAs related to cold adaptation (He et al., 2018), to characterize the transcriptomes of the liver and pancreas under low temperature stress (Zhuo et al., 2021), and to investigate proteomic features (Peng et al., 2018; Fan et al., 2019).
Single-cell transcriptome sequencing is useful because even cells derived from the same individual or cell line may differ genomically, transcriptomically, and epigenetically (Wen and Tang, 2016) due to cellular heterogeneity (Kalisky et al., 2011). Conventional multi-cell transcriptomic analyses require large numbers of cells, and the results thus reflect average gene-expression values across the cell population used. As these average values mask differences in expression among cellular subgroups, it is difficult to distinguish the transcriptional characteristics of these subgroups within the cellular population (Lu and Fuchou, 2014; Wang and Navin, 2015). Single-cell transcriptome sequencing (scRNAseq) solves this problem by revealing expression differences among cell populations at the level of the individual cell (Gawad et al., 2016). In particular, scRNAseq can be used to isolate, identify, and describe different cell populations by sequencing the transcriptome libraries prepared from thousands of single cells in parallel (Trapnell, 2015; Moris et al., 2016; Tanay and Regev, 2017). Because scRNAseq analysis allows for the posterior identification of each cell in a cellular population, it is possible to study heterogeneous cell populations to discover unknown cell types or to compare different points in the developmental trajectory of a single cell or cell type (Briggs et al., 2018; Farrell and Wang, 2018; Wagner et al., 2018).
To date, single-cell transcriptome sequencing has rarely been applied to aquatic animals, and the few available studies have tended to focus on model fish or ascidians. For example, single-cell sequencing was used to reveal the heterogeneous effects of bisphenol A on zebrafish embryo development (Chen et al., 2020); to characterize the diversity of newborn neurons in the adult zebrafish brain (Lange and Rost, 2020); to show that homologous genes in human and zebrafish oocytes are highly expressed and belong to similar functional categories (Can et al., 2020); to reconstruct the trajectory of zebrafish embryogenesis and development (Farrell and Wang, 2018); to produce a transcriptome of zebrafish gene expression from the pharyngeal embryo stage to the larval stage (Farnsworth et al., 2020); to create a transcriptome map of the larval brain of Ciona (Sharma et al., 2019); to determine the regulatory network of dopaminergic neurons in embryonic ascidians (Horie et al., 2018); and to analyze head kidney-derived leukocytes in tilapia (Niu et al., 2020).
However, single-cell transcriptome sequencing has not previously been used to analyze the cold response of in L. vannamei. Here, we chose to focus on the hepatopancreas, as this organ plays an important role in immunity, hematopoiesis, metabolism, energy storage, detoxification, and excretion in shrimp (Yepiz-Plascencia et al., 2000; Verri et al., 2001; Shekhar et al., 2013). We therefore aimed to construct a single-cell map of the hepatopancreas of L. vannamei and to analyze the transcriptomes of its highly heterogeneous cell population; to compare singe-cell transcriptomes between cold-tolerant and common L. vannamei; and to characterize the cold tolerance mechanism of L. vannamei.
Method Details
Shrimp Used and Rapid Cooling Stress
Cold-tolerant L. vannamei (Lv-T) and common (Lv-C) L. vannamei were provided by Guangxi Shrimp Genetics and Breeding Center (Fangchenggang City, Guangxi Province, China). The average survival time (30.07 h) of Lv-T shrimp at 10°C are significantly greater than those (20.17 h) of Lv-C shrimp (Supplementary Table S1). All shrimp (average weight ∼32 g) were kept in the seawater pools (5 m × 3 m × 2 m) at a water temperature of 28 ± 0.5°C, salinity of 30‰, and pH of 7.9 ± 0.1 for 1 month to adapt to the experimental environment. After the acclimation period, an air-cooled condenser (FNV-73.9/240, Shengzhou Xinba Refrigeration Equipment Co., Ltd., Shaoxing City, Zhejiang Province, China) was used to reduce the water temperature, at a rate of 1°C per hour, from 28°C to 10°C and then to hold the water temperature at 10°C for 24 h. We randomly selected shrimp from among the cold-tolerant and common shrimp before cooling and after 24 h at 10°C to establish four treatment groups (n = 5 shrimp per group): cold-tolerant shrimp at 10°C (Lv-T10), cold-tolerant shrimp at 28°C (Lv-T28), common shrimp at 10°C (Lv-C10), and common shrimp at 28°C (Lv-C28). The hepatopancreas was removed from each selected shrimp, immediately transferred to MACS Tissue Storage Solution (Miltenyi Biotec, Germany), and held at 2–8°C for later use.
Preparation of Hepatopancreas Single-Cell Suspensions
We cut 100 mg of each stored shrimp hepatopancreas sample into small pieces. We added the pieces to pre-cooled HBSS buffer (4°C) and gently pipetted the pieces in and out five times to wash the blood from the hepatopancreas. After natural sedimentation, the supernatant was discarded. Enzyme solution (2.5 g/L trypsin plus 0.1 g/L EDTA at a ratio of 1:1) was added to the clean tissue and allowed to digest for 10 min. When the tissue had swelled noticeably and appeared bleached, 2 ml of pre-cooled HBSS buffer solution was added and repeatedly pipetted gently to dissociate and release the hepatopancreas cells. After natural sedimentation for about 15 s, the cells were transferred to a new centrifuge tube using a pipette. This step was repeated three times for each sample (using the same centrifuge tube) to recover as many hepatopancreas cells as possible. The resulting cell suspension was filtered using a cell sieve with a pore size of 40 μm. The filtered liquid was centrifuged at 500 g/min for 3 minutes, the supernatant was discarded, and the remaining bottom layer (the cell pellet) was retained. Each cell pellet was resuspended in HBSS buffer, stained with trypan blue, and observed under a microscope to calculate cell concentration and percent viability. Cell suspension with cell viability >90% and cell concentration of 1,000–5000 cell/μL was used for future sequencing.
cDNA Library Preparation and Single Cell RNA-Sequencing
Cell concentrations were adjusted to approximately 1,000 cell/μL using HBSS buffer. Gel beads containing barcode information were bound to the mixture of cells and enzymes, and bound molecules were then wrapped in oil surfactant droplets using the microfluidic “double cross” system to create gel beads-in-emulsions (GEMs). The GEMs flowed into the reservoir and were collected, the gel beads were dissolved to release the barcode sequences, and the cDNA fragments were reverse transcribed. The cDNA sequences were used as templates for PCR amplification, and the PCR amplicons were mixed to construct a sequencing library. The established libraries were sequenced using the Illumina sequencing platform.
Cell Clustering
The raw sequences were cleaned by comparison to the L. vannamei genome (Zhang et al., 2019). Cell Ranger was used to filter, compare, quantify, identify, and reclaim cells; this software also generated a gene expression matrix for each cell. We then used Seurat to perform additional cell filtration and standardization, to classify cell clusters, to identify the genes differentially expressed among subgroups, and to screen marker genes. Cells that met the following conditions were retained for the future analysis: containing 400–1,000 genes; containing 750–3,000 unique molecular identifiers (UMI); and <20% of the UMIs were mitochondrial genes. The gene expression data for each retained cell were normalized against total gene expression levels across all retained cells using the global-scaling normalization method “LogNormalize” in the R software package Seurat 3.0 (https://satijalab.org/seurat/). Specifically, the gene expression data were multiplied by a scaling factor (10,000) and loge transformed as follows: gene expression level = ln (1 + (UMIA ÷ UMITotal) × 10,000), where UMIA was the UMI of gene A in the target cells; and UMITotal was the sum of all UMI quantities in the target cell. We performed a canonical correspondence analysis (CCA) of the normalized expression levels to eliminate the batch effect, and then integrated the expression data. The integrated data were Z-score normalized and a principal component analysis (PCA) was performed to reduce variable dimensionality. Finally, graph theory clustering algorithms were used to cluster cells across all four treatment groups. We then used t-SNE (McDavid et al., 2013) to visualize and explore the cell clusters.
Differentially Expressed Transcripts and Functional Annotation
The likelihood-ratio test (Fonville et al., 2013) was used to identify patterns of differential gene expression among cell clusters. Significantly upregulated genes in each cluster were identified as those meeting all of the following criteria: expressed in >25% of all cells in the experimental (cold) treatment or in the control treatment, and log2FC was ≥0.360,674 (p ≤ 0.01) between the control and experimental treatment groups (i.e., upregulated at least 1.28-fold in response to cold). The five most strongly upregulated genes (i.e., with the greatest logFC values) in each cell cluster were selected as potential marker genes for that cluster. We then constructed heat maps and bubble charts for the marker genes.
We annotated all of the transcripts differentially expressed between the cold (experimental) treatment and the control treatment across all cell clusters using the GO database (http://www.geneontology.org/) (Falcon and Gentleman, 2007). We counted the number of transcripts associated with each term and performed a hypergeometric test to identify GO terms that were significantly enriched in the differentially expressed transcripts (FDR-corrected p ≤ 0.05). We determined the primary biological functions of each of the differentially expressed transcripts based on the significantly enriched GO functions. We then used the hypergeometric test, followed by correction for the false discovery rate (FDR), to identify pathways in the KEGG database (Kanehisa and Goto, 2000) that were significantly enriched in differentially expressed transcripts (Q-value ≤0.05). Based on the pathways significantly enriched in the differentially expressed transcripts, particularly those associated with biochemical metabolisms and signal transduction, we manually grouped the cell clusters into the four major cell types found in the crustacean hepatopancreas tubules: resorptive cells (R-cells), blister-like cells (B-cells), fibrillar cells (F-cells), and embryonic cells (E-cells).
Differentially Expressed Gene Analysis
To explore molecular differences in the cold stress response between cold-tolerant (Lv-T) and common (Lv-C) L. vannamei, we performed scRNA-based transcriptome analyses on samples from the four treatment groups (Lv-T10, Lv-T28, Lv-C10, and Lv-C28). We identified DEGs in four sets of comparisons: Lv-C28-vs-Lv-C10 using Seurat (Stuart et al., 2019); genes were considered significantly differentially expressed when |log2FC| was ≥0.36, the proportion of cells expressing the target gene in a given cluster was ≥0.1, and ≥25% of all cells in that same cluster expressed the target gene. We then used the hurdle model in Model-based Analysis of Single-cell Transcriptomics (MAST) (Finak et al., 2015) to identify the DEGs, based on each set of comparisons, in each cluster. p-values were corrected for multiple tests using the Benjamini-Hochberg method in Seurat; and genes were considered DEGs when the corrected p-value was ≤0.05. We then analyzed the GO and KEGG enrichment of the DEGs as described in Section 2.5.
Pseudotime Analysis
Single-cell gene expression profiles allow the deconstruction of population heterogeneity and the reconstruction of cell differentiation trajectories. As F-cells are important for hepatopancreatic function in L. vannamei cells, we chose to investigate the differentiation trajectory of these cells in response to cold stress using pseudotime (cell-cell trajectory) analysis (Trapnell et al., 2014) in Monocle (Tryselius and Hultmark, 1997; Qiu et al., 2017). In brief, we used Monocle to reduce the single-cell gene expression data matrix to two dimensions. The cells were plotted as a minimum spanning tree, and difference in expression were determined based on biological meaning. We considered the original cell cluster that with the lowest pseudotime in the cell trajectory, and the pseudotime values for all cells were calculated from this. Based on pseudotime values, Monocle models gene expression levels as smooth, non-linear pseudotime functions to visualize changes in gene expression level over pseudotime. Genes with FDR<1e-5 were identified as differentially expressed genes over pseudotime.
Results
Ten Cell Clusters Were Identified Across 5185 Cells From the Four Treatment Groups
Illumina sequencing identified 61,193 genes in 252,817 L. vannamei hepatopancreas cells across the four treatment groups: cold-tolerant shrimp at 10°C (Lv-T10), cold-tolerant shrimp at 28°C (Lv-T28), common shrimp at 10°C (Lv-C10); and common shrimp at 28°C (Lv-C28). Of these cells, 5185 hepatopancreas cells were identified as high quality, based on the number of genes per cell (100–400; Figure 1B), the number of UMIs per cell (750–3,000), and the proportion of UMIs that were mitochondrial genes (<20%). These 5185 cells clustered into 10 subgroups (Figure 1A) using Seurat, a clustering algorithm based on graph theory.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792172/fgene-12-792172-HTML/image_m/fgene-12-792172-g001.jpg)
FIGURE 1. Single-cell sequencing reveals the subtypes of shrimp hepatopancreatic cells. (A) t-SNE visualization showing the 10 single-cell clusters identified in the 5185 hepatopancreas cells. Each dot represents a single cell; clusters are distinguished by color. (B) The number of genes per cell in each of the four treatments. (C) Histogram showing the number of genes upregulated in the cold (experimental) treatment group as compared to the control group in each of the single-cell subgroups shown in panel (A). (D) Heat map showing the expression patterns of the five most strongly upregulated genes in each cell cluster. Each column represents a cell, and each row represents a gene. Brighter yellow cells are more upregulated; brighter pink cells are more downregulated. (E) Bubble chart showing the expression levels of candidate molecular marker genes across all cell subgroups. The size of the bubble represents the ratio of the total expression of each marker gene in the corresponding cluster to its total expression across all clusters. That is, larger bubbles indicate genes that were relatively frequently expressed in the corresponding cluster, while darker bubbles indicate genes that were relatively highly expressed in the corresponding cluster. (F) t-SNE visualization showing the cell clusters classified into the four major types of hepatopancreas cells: resorptive cells (R-cells; orange), blister-like cells (B-cells; blue), fibrillar cells (F-cells; green), and embryonic cells (E-cells; red). Cell clusters were classified as cell types based on the molecular markers and KEGG signaling pathways. (G) Bubble chart showing the expression levels of candidate molecular marker genes across all cell subgroups. The size of the bubble represents the ratio of the total expression of each marker gene in the corresponding cluster to its total expression across all clusters. That is, larger bubbles indicate genes that were relatively frequently expressed in the corresponding cluster, while darker bubbles indicate genes that were relatively highly expressed in the corresponding cluster.
To analyze the transcriptional regulation of each cell cluster, and to screen marker genes unique to each cluster, we first identified the genes in each cluster that were significantly upregulated in response to cold stress. In cell clusters 0–9, we identified 76, 213, 227, 178, 266, 142, 120, 172, 368, and 145 significantly upregulated genes, respectively, in the cold treatment group as compared to the control treatment group (Figure 1C). According to the expression analysis of genes in each cluster, we further selected the five most strongly upregulated genes in each cluster as marker genes. A heatmap (Figure 1D) shows the expression patterns of the marker genes, and an enriched bubble chart (Figure 1E) and a bubble chart (Figure 1G) show the expression levels of marker genes across all cell subgroups. These genes may be useful molecular markers for these cell clusters in future studies.
The 10 Cell Clusters Were Classified Into Four Cell Types
According to pathways and molecular markers, we classified the 10 cell clusters into 4 cell types, including R-cells, B-cells, F-cells, and E-cells (Figure 1F) as follows (Niu et al., 2020): Several pathways related to the fat metabolism, including “fatty acid metabolism” (ko01212), “fatty acid degradation” (ko00071), and the “glyoxylate and dicarboxylate metabolism” (ko00630), were significantly enriched in cell cluster 0; in particular, the glyoxylate and dicarboxylate metabolism is important for the conversion of fat to glycogen. Two molecular function GO terms associated with metal ions, “metalloexopeptidase activity” (GO:0,008,235) and “metallopeptidase activity” (GO:0,008,237), were also significantly enriched in cell cluster 0. As R-cells primarily absorb and store nutrients (e.g., lipid droplets, glycogen, and calcium) and trace elements (e.g., copper, zinc, sulfur, and calcium) (Al-Mohanna and Nott, 1987), but are also involved in fat metabolism (Sousa and Petriella, 2001; Tian et al., 2012), we classified cell cluster 0 as R-cells (Figure 1F). The pathway “glycosphingolipid biosynthesis-ganglio series” (ko00604) was significantly enriched in cell subgroup 6. Glycosphingolipids are complex glycolipids (Chatterjee and Kwiterovich, 1984). As lipid storage is an important function of R-cells (Sousa and Petriella, 2001; Tian et al., 2012), we also classified cell cluster 6 as R-cells (Figure 1F).
Several pathways and GO cellular component terms associated with B-cells, including “phagosome” (ko04145), “endocytosis” (ko04144), “lysosome” (ko04142), “synaptic vesicle cycle” (ko04721), “lysosome” (GO:0,005,764), “secondary lysosome” (GO:0,005,767), “vacuum lumen” (GO:0,005,775), and “autophagosome” (GO:0,005,776) were significantly enriched in cell cluster 5. As mature B-cells contain large amounts of rough endoplasmic reticulum (rER), Golgi apparatus, mitochondria, and lysosomes, and exhibit endocytosis (Vogt, 2019), we classified cell cluster 5 as B-cells (Figure 1F).
The “steroid hormone biosynthesis” (ko00140) pathway was significantly enriched in cell cluster 4. The synthesis of steroid hormones is controlled by several highly selective cytochrome P450 enzymes, in conjunction with some steroid dehydrogenases and reductases, especially aromatase (Sanderson, 2006). In addition, the hepatopancreas detoxifies organic compounds using cytochrome P450 (James and Boyle, 1998; Brouwer and Lee, 2007). We also identified 189 transcripts associated with detoxification in cell clusters 1 and 4. These transcripts were significantly enriched in the pathways “drug metabolism-cytochrome P450” (ko00982), “metabolism of xenobiotics by cytochrome P450” (ko00980), and “longevity regulating pathway-worm” (ko04212) in cell clusters 1 and 4. In cell cluster 1, the “glutathione S-transferase genes” pathway (including the three subtypes LOC113817634, LOC113822637, and LOC113829084) was significantly enriched in the “longevity regulating pathway-worm” (ko04212) pathway. There is evidence that the “glutathione S-transferase gene” is localized in the F-cells of the crab Callinectes sapidus, and this enzyme plays an important role in the detoxification of organic heterologous substances (Keeran and Lee, 1987). Indeed, in several crustaceans the F-cells have been identified as the possible detoxification site for organic foreign bodies (Brignac-Huber et al., 2016; Vogt, 2019). We thus classified cell clusters 1 and 4 as F-cells (Figure 1F).
The “carbohydrate digestion and absorption” (ko04973) and “starch and sucrose metabolism” (ko00500) pathways were significantly enriched in cell cluster 8, while genes related to digestive enzymes, including carboxypeptidase B (CP-B), amylase (α-AMY), and cathepsin 1 (CTSL1), were significantly upregulated in cell subset 9. CTSL1, as well as CTSB, CLEC4M, Mgst1, Mgst3, CYP9E2, and CYP9E2 (Figure 1E), were also upregulated in cell cluster 3. CTSL1 is a typical lysosomal enzyme that is found in the digestive juices of decapod animals (Hu and Leung, 2007). in addition, previous in situ hybridization studies based on chymotrypsin (Van wormhoudt et al., 1995), amylase, chitinase, cellulase, and trypsin (Lehnert and Johnson, 2002) showed that F-cells were the main site of digestive enzyme synthesis. F-cells also synthesize digestive enzymes and blood proteins, and participate in oxygen transport and immune defense (Vogt et al., 1989; Möhrlen et al., 2001; Vogt, 2019). Finally, previous studies have shown that F-cells in Penaeus penicillatus secrete zymogen granules through exocytosis (Al-Mohanna et al., 1985). Thus, as the main function of F-cells is to synthesize and secrete digestive enzymes (Al-Mohanna et al., 1985; Toullec et al., 1992), we classified cell clusters 3, 8, and 9 as F-cells (Figure 1F).
Several signaling pathways related to immunity, including “B cell receptor signaling pathway” (ko04662), “viral carcinogenesis” (ko05203), “human cytomegalovirus infection” (ko05163), and “HTLV-I infection” (ko05166) were significantly enriched in cluster 7, while “B cell receptor signaling pathway” (ko04662), “human cytomegalovirus infection” (ko05163), “kaposi sarcoma-associated herpesvirus infection” (ko05167), “pathways in cancer” (ko05200), “leukocyte transendothelial migration” (ko04670), and “cell adhesion molecules (CAMs)” (ko04514) were significantly enriched in cluster 9. The most abundant protein in hemolymph, while plays a key role in the crustacean immune system, is the oxygen-carrier hemocyanin (>90% of all proteins in hemolymph) (Stöcker et al., 1988). As a previous study used in situ hybridization of shrimp hepatopancreas slices with hemocyanin to show that hemocyanin was synthesized in F-cells (Lehnert and Johnson, 2002; Wang et al., 2007), we classified the cell clusters 7 and 9 as F-cells.
The “hedgehog signaling pathway” (ko04340), which regulates cell fate, including proliferation, apoptosis, migration, and differentiation (Qu et al., 2021; Wutka et al., 2014), was significantly enriched in cell cluster 2. Ccnd2, which was significantly upregulated in cell cluster 2, is a key gene in the “hedgehog signaling” pathway (Li et al., 2020) (Figure 1E). E-cells are undifferentiated embryonic cells from which R-, B-, and F-cells are differentiated (Toullec et al., 1992; Sousa and Petriella, 2000; Tian et al., 2012). Therefore, we identified cell cluster 2 as E-cells.
Differentially Expressed Genes Associated With Cold Tolerance
Across all cell clusters, we identified 1,145 DEGs between common L. vannamei at 28°C and common L. vannamei at 10°C (Lv-C28-vs-Lv-C10), 525 upregulated and 620 downregulated (Figure 2A); 1,209 DEGs between cold-tolerant L. vannamei at 28°C and cold-tolerant L. vannamei at 10°C (Lv-T28-vs-Lv-T10), 616 upregulated and 593 downregulated (Figure 2B); 655 DEGs between cold-tolerant L. vannamei at 10°C and common L. vannamei at 10°C (Lv-C10-vs-Lv-T10), 372 upregulated and 283 downregulated (Figure 2C); and 926 DEGs between common L. vannamei at 28°C and cold-tolerant L. vannamei at 28°C (Lv-C28-vs-Lv-C28), 442 upregulated and 484 downregulated (Figure 2D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792172/fgene-12-792172-HTML/image_m/fgene-12-792172-g002.jpg)
FIGURE 2. Differentially expressed genes (DEGs) between cold-tolerant and common Litopenaeus vannamei at 10°C and 28°C. (A–D) Differentially expressed genes across all cell cultures (Total) and in cell clusters 0–9 for the following comparisons: (A) Lv-C28-vs-Lv-C10 (common at 28°C versus common at 10°C), (B) Lv-T28-vs-Lv-T10 (cold-tolerant at 28°C versus cold-tolerant at 10°C), (C) Lv-C10-vs-Lv-T10 (common at 10°C versus cold-tolerant at 10°C), and (D) Lv-C28-vs-Lv-T28 (common at 28°C versus cold-tolerant at 28°C). (E,F) Venn diagrams showing the (E) upregulated and (F) downregulated genes shared and unique among the four comparisons.
Six genes were co-upregulated in all four comparison groups and 14 were co-downregulated (Figures 2E,F). The shared upregulated genes included Kintoun (LOC113828196), CtBP (LOC113800533), Hnrnph1 (LOC113819995), and La (LOC113809655), while the shared downregulated genes included nrf-6 (LOC113816674), PHYH (LOC113811729), AK (LOC113816366), and Gs2 (LOC113823144). We identified 202 co-upregulated DEGs (Figure 2E) and 183 co-downregulated DEGs (Figure 2F) in cold-tolerant and common shrimp at 28°C as compared to 10°C (i.e., Lv-C28-vs-Lv-C10 and Lv-T28-vs-Lv-T10), including ACBP (LOC113815675), Leng8 (LOC113805896), GAT-2 (LOC113816986), Ankrd2 (LOC113828833), and ND1 (LOC113820286).
NR annotations indicated that, across all four comparisons, many of the DEGs were related to immune defense, the cytoskeleton, the antioxidant system, digestive enzymes, and detoxification, including C-type lectin (CTL), leukocyte receptor (Lr8), hemocyanin (HC), integrin beta-PS-like (ITGB), actin (ACT), titin-like (Ttn), catalase (CAT), glutathione S-transferase (GST), glutathioneperoxidase (GPX), peroxiredoxin 4 (PRDX4), preamylase 1 (Amy1), glutamate carboxypeptidase (GCP), chymotrypsin BII (Prss2), UDP-glucuronosyltransferases (UGTs), cytochrome P450 (CYP450) (Table 1). Other DEGs included low-density lipoprotein receptor 2 (LDLR2), calcium-activated 4A-like (Clca4), taurine transporter (TAUT), heat shock protein (HSC70), zinc proteinase (Zinc), and vital membrane outer layer protein 1-like (VMP1) (Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792172/fgene-12-792172-HTML/image_m/fgene-12-792172-t001.jpg)
TABLE 1. Genes differentially expressed in the hepatopancreas of Litopenaeus vannamei in response to cold stress. (N/A, not significantly differentially expressed).
Pathways Significantly Enriched in the Differentially Expressed Genes Were Associated With Cold-Stress Resistance
The DEGs in the Lv-C28-vs-Lv-C10 dataset were most commonly significantly enriched in the GO biological process terms “metabolic process (GO:0,008,152)”, “single-organism metabolic process” (GO:0,044,710), and “oxoacid metabolic process” (GO:0,043,436); in the GO cellular component term “proton-transporting two-sector atpase complex” (GO:0,016,469); and in the GO molecular function terms “catalytic activity” (GO:0,003,824) and “oxidoreductase activity” (GO:0,016,491; Figure 3A). The DEGs in the Lv-T28-vs-Lv-T10 dataset were most commonly significantly enriched in the GO biological process terms “single-organism metabolic process” (GO:0,044,710), “proton transport (GO:0,015,992), and “hydrogen transport” (GO:0,006,818); in the GO cellular component terms “proton-transporting two-sector atpase complex” (GO:0,016,469) and “proton-transporting two-sector atpase complex, catalytic domain” (GO:0,033,178); and in the GO molecular function terms “catalytic activity” (GO:0,003,824), “oxidoreductase activity” (GO:0,016,491), and “oxidoreductase activity, acting on CH-OH group of donor” (GO:0,016,614; Figure 3B). The DEGs in the Lv-C10-vs-Lv-T10 dataset were most commonly significantly enriched in the GO biological process terms “mRNA metabolic process” (GO:0,016,071), “RNA catabolic process” (GO:0,006,401), and “multicellular organismal macromolecule metabolic process” (GO:0,044,259); in the GO cellular component terms “proton-transporting two-sector atpase complex” (GO:0,016,469) and “endosomal part” (GO:0,044,440); and in the GO molecular function terms “hydrolase activity” (GO:0,016,787), “cation transmembrane transporter activity” (GO:0,008,324), and “peptidase activity” (GO:0,008,233; Figure 3C). The DEGs in the Lv-C28-vs-Lv-T28 dataset were most commonly significantly enriched in the GO biological process terms “metabolic process” (GO:0,008,152), “monovalent inorganic cation transport” (GO:0,015,672), and “proton transport” (GO:0,015,992); in the GO cellular component term “proton-transporting two-sector atpase complex” (GO:0,016,469); and in the GO molecular function terms “oxidoreductase activity” (GO:0,016,491), “ion transmembrane transporter activity” (GO:0,015,075), and “cation transmembrane transporter activity” (GO:0,008,324; Figure 3D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792172/fgene-12-792172-HTML/image_m/fgene-12-792172-g003.jpg)
FIGURE 3. Circle plots showing GO enrichment of the genes differentially expressed between cold-tolerant and common Litopenaeus vannamei at 10°C and 28°C. (A) Lv-C28-vs-Lv-C10 (common at 28°C versus common at 10°C), (B) Lv-T28-vs-Lv-T10 (cold-tolerant at 28°C versus cold-tolerant at 10°C), (C) Lv-C10-vs-Lv-T10 (common at 10°C versus cold-tolerant at 10°C), and (D) Lv-C28-vs-Lv-T28 (common at 28°C versus cold-tolerant at 28°C). The outermost circle of each graph shows the 20 GO terms most commonly significantly enriched in the differentially expressed genes (DEGs). The second circle shows the number of genes associated with each GO term; redder bars correspond to greater Q-values. The third circle shows the ratio of upregulated to downregulated DEGs associated with each GO term as bar graph, with upregulated genes shown in light purple and downregulated genes shown in dark purple; the specific numbers of up- and downregulated genes are shown below bar. The innermost (fourth) circle shows the Rich Factor value for each GO term (i.e., the number of DEGs versus the number of non-DEGs associated with the GO term); each gridline represents 0.1.
Nine KEGG pathways were significantly enriched (Q-value ≤ Q-val in the Lv-C28-vs-Lv-C10, Lv-T28-vs-Lv-T10, and Lv-C28-vs-Lv-T28 DEG datasets, including “oxidative phosphorylation” (ko00190), “metabolic pathways” (ko01100), “biosynthesis of amino acids” (ko01230), “metabolism of xenobiotics by cytochrome P450 " (Ko00980), “chemical carcinogenesis” (ko05204), “drug metabolism-cytochrome P450” (ko00982), and “fatty acid metabolism” (ko01212; Figure 4). These pathways may be related to cold stress. The pathways “pentose and glucuronate interconversions” (ko00040), “steroid hormone biosynthesis” (ko00140), and “bile secretion” (Ko04976) were significantly enriched in the Lv-C10-vs-Lv-T10 and Lv-C28-vs-Lv-T28 DEG datasets after cold stress, but not in the Lv-C28-vs-Lv-C10 DEG dataset. That is, these pathways were significantly enriched in the cold-tolerant shrimp, but not the common shrimp, in response to cold stress. Thus, these pathways may be important for cold-stress resistance in cold-tolerant L. vannamei families. Notably, these three pathways were also significantly enriched in cell cluster 0 (i.e., R-cells; Figure 4). In the DEG datasets Lv-C10-vs-Lv-T10 and Lv-C28-vs-Lv-T28, R-cells were significantly enriched in the “metabolism of xenobiotics by cytochrome P450” (ko00980), “chemical carcinogenesis” (ko05204), “drug metabolism-cytochrome P450” (ko00982), “antigen processing and presentation” (ko04612), and “drug metabolism-other enzymes” (ko00983) pathways (Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792172/fgene-12-792172-HTML/image_m/fgene-12-792172-g004.jpg)
FIGURE 4. Q-value heatmaps of KEGG pathways associated with the genes differentially expressed between cold-tolerant and common Litopenaeus vannamei at 10°C and 28°C across all cell clusters (Total) and in individual cell clusters (0–9). (A) common at 28°C versus common at 10°C, (B) cold-tolerant at 28°C versus cold-tolerant at 10°C, (C) common at 10°C versus cold-tolerant at 10°C, and (D) common at 28°C versus cold-tolerant at 28°C. Q-values (i.e., p-values corrected for the false discovery rate) are shown in cells: cells where Q ≤ 0.05 are shown in red; cells where Q > 0.05 are shown in green; cells not corresponding to differential expression are shown in grey. Corresponding KEGG signaling pathways are shown in the leftmost column of each graph.
Across all DEG datasets, R-cells were significantly enriched in the “antigen processing and presentation pathway” (ko04612) and “PPAR signaling pathway” (ko03320), while F-cells were significantly enriched the “PPAR signaling pathway” (ko03320). B-cells were not significantly enriched in any KEGG pathways in any DEG dataset. was significantly enriched in all comparison groups.
F-Cell Trajectories Differed Between Cold-Tolerant and Common F-Cells in Response to Cold Stress
After F-cell differentiation in response to cold-stress (the pre-branch), 2 cell fates were identified: cell fate 1 and cell fate 2 (Figure 5A). F-cells from clusters 1 and 3 were primarily found at the beginning of the pre-branch, transitioning to cell clusters 7 and 8 just before differentiation. After differentiation, cell fate 1 primarily included cell cluster 4, while cell fate 2 primarily included cell cluster 9 (Figure 5B). The pre-branch primarily included cells from both common and cold-tolerant shrimp at 28°C; cells from the common shrimp at 10°C were mainly found on the cell fate 1 branch, while cells from the cold-tolerant shrimp at 10°C were mainly found on the cell fate 2 branch (Figure 5C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792172/fgene-12-792172-HTML/image_m/fgene-12-792172-g005.jpg)
FIGURE 5. Reconstruction of the differentiation trajectory of F-cells in a quasi-chronological analysis. (A–C) Differentiation trajectory of F-cells with different mappings: (A) pseudotime, (B) cell clusters, and (C) treatment groups. Three cell trajectories are shown in these graphs: the pre-branch, cell fate 1, and cell fate 2. (D) Heatmap showing simulated gene expression dynamics in F-cells during cell state transitions from the pre-branch to cell fate 1 or 2. Pseudotime is shown on the x-axis (with the origin at the center of the axis), and genes are graphed along the y-axis. Gene sets with similar expression patterns are clustered. (E,F) The top 10 GO terms with the smallest Q value were used for mapping. GO terms are shown on the y-axis; the percentages of the number of GO terms to the total number of differences are shown on the x-axis (the darker the color, the smaller the Q value). The digitals on the columns are the GO term numbers and Q values. E is the branch of cell fate 1, and F is the branch of cell fate 2.
A heatmap showing gene expression dynamics over pseudotime showed distinct patterns of gene expression among the three branches (Figure 5D). Several genes exhibited similar expression patterns and cell fates over pseudotime in both common and cold-tolerant shrimp. For example, apolipoprotein D (APOD), cysteine proteinase 4-like (CP4), carboxypeptidase B (CPB), pancreatic triacylglycerol lipase-like (PTL), vitelline membrane outer Layer 1 (VMO1), and zinc proteinase (Zinc) were upregulated, while cytochrome P450 9e2-like (CYP450), hemocyanin (HC), and V-type proton atpase (V-ATPase) were downregulated (Figure 6A). Other genes exhibited different expression patterns in cells with different fates in response to cold stress. For example, astacin-like (AST), chitinase precursor (Chia), carboxypeptidase B (CPB), C-type lectin receptor (CTL), low-density lipoprotein receptor-related protein 2-like (LDLR2), macrophage mannose receptor 1-like (MMR1), and trypsin-1-like (Trypsin-1) were upregulated in cells on the cell fate 2 branch, but downregulated in cells on the cell fate 1 branch (Figure 6B). Conversely, anti-lipopolysaccharide factor (ALF), calcium-activated chloride channel regulator 1 (CLCA1), and headcase protein-like (HECA) were downregulated in cells on the cell fate 2 branch, but upregulated in cells on the cell fate 1 branch (Figure 6B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/792172/fgene-12-792172-HTML/image_m/fgene-12-792172-g006.jpg)
FIGURE 6. Expression patterns of genes in cold-tolerant and wild-type Litopenaeus vannamei in response to cold stress over pseudo-time. Yellow dot represents common at 10°C, blue dot represents common at 28°C, red dot represents cold-tolerant at 10°C, and green dot represents cold-tolerant at 28°C. (A) Gene expression patterns in cells with similar fates after cold stress. Solid black line represents best fit line for gene expression level. (B) Gene expression patterns in cells with different fates after cold stress. Solid black line represents the best fit line for the expression levels of genes in cells associated with cell fate 1, and the dashed line represents the best fit line for the expression levels of genes in cells associated with cell fate 2.
GO enrichment analysis indicated that cell fate 1 was mainly associated with signal pathway transduction and sensory organ development: the GO terms “intracellular signal transduction" (GO:0,035,556), “sensory organ development" (GO:0,007,423), “Ras protein signal transduction" (GO:0,007,265), and “small gtpase mediated signal transduction" (GO:0,007,264) were significantly enriched in cell fate 1 (Figure 5E). Cell fate 2 was mainly associated with metabolic processes: “collagen metabolic process" (GO:0,032,963), “multicellular organism metabolic process" (GO:0,044,236), “multicellular organismal macromolecule metabolic process" (GO:0,044,259), and “protein metabolic process" (GO:0,019,538) were significantly enriched in cell fate 2 (Figure 5F). Thus, pseudotime analysis showed that F-cells had different trajectories in cold-tolerant and common L. vannamei in response to low-temperature stress.
Discussion
Cell Clusters Identification in the Hepatopancreas of L. Vannamei
Although electron microscopic observation of hepatopancreas tubules have shown that the crustacean hepatopancreas is primarily composed E-cells, R-cells, B-cells, and F-cells (Al-Mohanna and Nott, 1989; Toullec et al., 1992; Sousa and Petriella, 2000; Hopkin and Nott, 2009), precise molecular markers for these cells types have yet to be identified. To address this knowledge gap, we comprehensively analyzed gene expression in the shrimp hepatopancreas at the level of a single cell, using scRNA-seq to identify candidate markers for each cell type in the shrimp hepatopancreas.
R-cells play important roles in the absorption and metabolism of nutrients, in energy (lipids and glycogen) and mineral element storage, in lipoprotein synthesis and export, in heavy metal detoxification, and in uric acid excretion (Al-Mohanna and Nott, 1987; Sousa and Petriella, 2001; Tian et al., 2012). Here, the pathways significantly enriched in the putative R-cells, including “glyoxylate and dicarboxylate metabolism,” “cysteine and methionine metabolism,” “metalloexopeptidase activity,” “metalloexopeptidase activity,” and “glycosphingolipid biosynthesis-ganglio series” were consistent with the functions of R-cells. The pathways significantly enriched in the putative F-cells, including “drug metabolism-cytochrome P450” and “metabolism of xenobiology by cytochrome P450” were consistent with the detoxification function of F-cells (James and Boyle, 1998; Brouwer and Lee, 2007; Brignac-Huber et al., 2016). Previous in situ hybridization analyses have also implicated F-cells in immune defense (Wang et al., 2007). Here, several pathways related to viral infection and immune signaling were significantly enriched in the F-cells.
Some previous studies have suggested that B-cells secrete digestive enzymes, similar to F-cells (Barker and Gibson, 1979; Dall et al., 1983), while other studies have indicated that only F-cells have this function; indeed, immunohistochemical and in situ hybridization studies have shown that B-cells cannot synthesize classic digestive enzymes (Vogt et al., 1989). In addition, B-cells participate in the absorption of nutrients from the hepatopancreatic duct and concentrate absorbed substances in large vacuoles, but do not store nutrients (Al-Mohanna and Nott, 1986; Hopkin and Nott, 2009). E-cells are undifferentiated embryonic cells that can differentiate into R-cells, F-cells, and B-cells (Al-Mohanna and Nott, 1989; Toullec et al., 1992; Sousa and Petriella, 2000).
Several Signaling Pathways Related to Cold Tolerance Were Enriched in L. Vannamei
Several pathways, including “oxidative phosphorylation,” “metabolism of xenobiology by cytochrome P450,” “chemical carcinogenesis,” “drug metabolism cytochrome P450,” and “fatty acid metabolism,” were significantly enriched in the cold-tolerant shrimp. Previous studies suggest that these signaling pathways may be important for cold tolerance. For example, oxidative phosphorylation provides ATP to meet increased energy demands under cold stress (Hall et al., 2012; Colinet et al., 2017; Chamberlain and Sheng., 2019), cytochrome P450 participates in detoxification and removal of harmful substances (Nanji and Hiller-Sturmhöfel., 1997; Nebert et al., 2004), and the fatty acid metabolism generates heat to cope with the low-temperature environments (Hui et al., 2014; Lee et al., 2019). Drosophila melanogaster uses the ATP produced by oxidative phosphorylation to meet energy requirements under cold stress (Williams et al., 2018). Indeed, it has been shown that cold tolerance in L. vannamei is associated with oxidative phosphorylation (Fan et al., 2019; Peng et al., 2018). Cytochrome P450 may be upregulated in shrimp in response to cold stress; cytochrome P450 may also participate in detoxification and the removal of harmful substances (Vergara-Amado et al., 2017), as well as the biotransformation of endogenous chemical products (such as reactive oxygen species) produced under chemical stress (Andersson and Förlin, 1992; Zhang et al., 2012). Increased cytochrome P450 expression has been detected in marine oligochaetes (Thalassodrilides sp.) exposed to pollutants and heat stress (Ito et al., 2016). Some functional analyses have shown that cytochrome P450 genes closely related to the fat metabolism (e.g., CYP7A1 and CYP1A1) play important roles in cold stress (Chen et al., 2012). Interestingly, exposure to low temperatures inhibits gene expression and most cellular processes in rats, with the exception of cytochrome P450 enzymes (CYP1A) in liver, which are instead activated by cold (Perepechaeva and Grishanova, 2013). Here, detoxification-associated genes were upregulated at 28°C or 10°C in the cold-tolerant shrimp (LV-T) as compared to the common shrimp (LV-C; Table 1). However, these genes were downregulated in response to low temperatures in both common and cold-tolerant shrimp (Table 1). This suggested that cold-tolerant shrimp may have greater detoxification capacity than the common shrimp, but that detoxification capacity in both lines was decreased by low-temperature stress.
The fatty acid metabolism is critical for cold adaptation (Gombos et al., 1994; Tuo et al., 2014; Zhang et al., 2020). Long-chain acyl-coenzyme A synthetase participates in the first step of the fatty acid metabolism and is one of the key enzymes for fat synthesis and catabolism (Wang et al., 2004). Studies have shown that the activity levels of long-chain acyl-CoA synthase increase in the rainbow trout heart during low temperature adaptation (Hicks et al., 1996; Patey and Driedzic, 2011), while fatty acyl-CoA synthetase-1 (adipose acyl-CoA synthetase-1, ACSL1) guides fatty acids to beta-oxidation, which is necessary for cold and heat generation (Ellis et al., 2010). In addition, stearoyl-CoA desaturase (SCD) is involved in the adaptation of large yellow croaker (Pseudosciaena crocea) to cold stress (Xu et al., 2015). Here, several genes associated with the fatty acid metabolism, including acyl-CoA dehydrogenase, long-chain-acyl-CoA dehydrogenase, and stearoyl-CoA desaturase, were significantly enriched in the cold-tolerant shrimp, indicating that the fatty acid metabolism plays an important role in the cold tolerance of L. vannamei. This was consistent with the results of Fan et al. (2019) and He et al. (2018).
Several immune-related genes, including C-type lectin (CTL), leukocyte receptor cluster member 8 (Lr8), hemocyanin (HC), integrin beta-PS-like (ITGB) and macrophage mannose receptor 1-like (MMR1), were upregulated in the cold-tolerant shrimp as compared to the common shrimp at 28°C or 10°C. In addition, CTL, HC, and MMR1 were upregulated in response to low-temperature stress, suggesting that the immune defense ability of the cold-tolerant shrimp was greater than that of the common shrimp. We speculate that low-temperature stress might affect the expression of immune-related genes, causing immune disorders and damaging the immune system.
Catalase (CAT), glutathione S-transferase (GST), glutathione peroxidase (GPX), and peroxiredoxin (PRDX4) are important antioxidants that can reduce ROS and prevent oxidative damage (Lim et al., 2018). GST, GPX, and PRDX4 were upregulated in cold-tolerant shrimp as compared to common shrimp at 28°C, suggesting that the immune defense abilities of cold-tolerant shrimp were better than those of common shrimp. In both cold-tolerant and common shrimp, these genes were downregulated in response to cold stress, indicating that low temperatures are associated with antioxidant system dysfunction in shrimp.
At both 10°C and 28°C, the taurine transporter (TAUT) was upregulated in the cold-tolerant shrimp as compared to the common shrimp. Taurine is an abundant, sulfur-rich free amino acid that is associated with antioxidant activity, ion transport, and DNA repair (Kuz’mina et al., 2010). Thus, our results suggested that the self-repair abilities of the cold-tolerant shrimp exceeded those of the common shrimp. In both common and cold-tolerant shrimp, heat shock protein 70 (HSC70) was upregulated in response to low temperature stress, suggesting that this protein is involved in the cold stress response mechanism. This was consistent with the results of Fan et al.(Fan et al., 2019).
Reconstruction of F-Cell Differentiation Fate Trajectory Under Low Temperature Stress
F-cells participate in cell detoxification, immune defense, and the secretion of digestive enzymes. However, the differentiation trajectory of F-cells under low temperature stress is not clear. Reconstruction of the differentiation trajectory of F-cells will help to clarify the transcription factors that cause F-cell differentiation in response to cold stress. In response to low temperature stress, CP4, CPB, CTL, PTL, VMO1, and Zinc were upregulated over pseudotime in all cells (Figure 6A). CPB, a an immune-related metallopeptidase (Qian et al., 2019), requires zinc ions (Zn2+) to participate in catalysis; CPB has a variety of functions, including promoting protein digestion and the processing and formation hormones and neuropeptides. Zinc also plays an important role in low temperature stress in Arabidopsis (Nakai et al., 2013) and L. vannamei (Peng et al., 2018). Overexpression of Cystatin B increased resistance to salinity, drought, oxidative stress, and cold stress in transgenic yeast and Arabidopsis (Zhang et al., 2008). CTL is an important immune factor, and studies have suggested that Osmerus mordax CTL is homologous to antifreeze proteins (AFPs) (Achenbach and Ewart, 2002). Low temperature induces CTL expression in Ruditapes philippinarum. In squirrels (Spermophilus tridecemlineatus), PTL is strongly upregulated and has high activity levels at low temperatures; PTL provides a mechanism to release the fatty acids required by this species during hibernation (Squire et al., 2003).
CYP450, HC, and V-ATPase were downregulated in response to cold stress in all cells over pseudotime (Figure 6B). Fan et al. (2013) and Fan et al. (2019) found that, in response to cold stress, shrimp blood cell concentrations decreased and differentially expressed genes associated with hematopoiesis were downregulated. Here, we speculated that the shrimp hematopoietic system was damaged by exposure to cold, leading to the downregulation of HC. In crickets (Gryllus pennsylvanicus), cold acclimation can reduce the production of primary urine by downregulating carbonic anhydrase and V-ATPase in Markov tubules, thereby reducing the loss of hemolymph volume in cold environments (Des Marteaux et al., 2017). AST and VMO1 participate in the adaptive response of Daphnia pulex to environmental temperature changes (Schwerin et al., 2009). LDLR is an important regulator of cell differentiation and proliferation (Tsukamoto et al., 2012; Auderset et al., 2020).
After cold exposure, AST, Chia, CPB, CTL, LDLR2, MMR1, and Trypsin-1 were upregulated in the cell fate 2 branch (primarily associated with cold-tolerant shrimp cells) and slightly downregulated in the cell fate 1 branch (primarily associated with common shrimp cells); the converse relationship was observed for ALF, CLCA1, and HECA (Figure 6B). These genes determine the fate of F-cells under low temperature stress. Thus, the differences in gene expression dynamics between cold-tolerant and common cells may help to identify the transcription factors that drive differences in F-cell differentiation in response to low-temperature stress.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, SRP330719.
Author Contributions
DZ and YZ conceived and directed the project; XC, QL, Qiangyong L, MP, and HW performed the experiments; XC, QY, ZL, ML, CP, and PF analyzed the data; and WZ and CY wrote the manuscript. All of the authors read and commented on the manuscript.
Funding
This work was supported by the Guangxi Science and Technology base and Talent Special Project (AD18294004), Guangxi Science and Technology base and Talent Special Project (AD20297021), National Shrimp and Crab Industry Technology System Construction Project (CARS-48), National Modern Agricultural Industrial Technology System Guangxi Innovation Team Construction Project (nycytxgxcxtd-14–01), and Independent research project of Guangxi Key Laboratory of Aquatic Breeding and Healthy Breeding (17-A-01–09).
Conflict of Interest
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.
Publisher’s Note
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.792172/full#supplementary-material
References
Achenbach, J. C., and Ewart, K. V. (2002). Structural and Functional Characterization of a C-type Lectin-like Antifreeze Protein from Rainbow Smelt (Osmerus mordax). Eur. J. Biochem. 269, 1219–1226. doi:10.1046/j.1432-1033.2002.02761.x
Al-Mohanna, S. Y., and Nott, J. A. (1986). B-cells and Digestion in the Hepatopancreas of Penaeus Semisulcatus (Crustacea: Decapoda). J. Mar. Biol. Ass. 66, 403–414. doi:10.1017/S0025315400043034
Al-Mohanna, S. Y., and Nott, J. A. (1989). Functional Cytology of the Hepatopancreas of Penaeus Semisulcatus (Crustacea: Decapoda) during the Moult Cycle. Mar. Biol. 101, 535–544. doi:10.1007/BF00541656
Al-Mohanna, S. Y., Nott, J. A., and Lane, D. J. W. (1985). Mitotic E- and Secretory F-Cells in the Hepatopancreas of the Shrimp Penaeus Semisulcatus (Crustacea: Decapoda). J. Mar. Biol. Ass. 65, 901–910. doi:10.1017/S0025315400019408
Al-Mohanna, S. Y., and Nott, J. A. (1987). R-cells and the Digestive Cycle in Penaeus Semisulcatus (Crustacea: Decapoda). Mar. Biol. 95, 129–137. doi:10.1007/BF00447494
Andersson, T., and Förlin, L. (1992). Regulation of the Cytochrome P450 Enzyme System in Fish. Aquat. Toxicol. 24, 1–19. doi:10.1016/0166-445X(92)90014-E
Auderset, L., Pitman, K. A., Cullen, C. L., Pepper, R. E., Taylor, B. V., Foa, L., et al. (2020). Low-Density Lipoprotein Receptor-Related Protein 1 (LRP1) Is a Negative Regulator of Oligodendrocyte Progenitor Cell Differentiation in the Adult Mouse Brain. Front. Cel Dev. Biol. 8, 564351. doi:10.3389/fcell.2020.564351
Briggs, J. A., Weinreb, C., Wagner, D. E., Megason, S., Peshkin, L., Kirschner, M. W., et al. (2018). The Dynamics of Gene Expression in Vertebrate Embryogenesis at Single-Cell Resolution. Science 360. doi:10.1126/science.aar5780
Brignac-Huber, L. M., Park, J. W., Reed, J. R., and Backes, W. L. (2016). Cytochrome P450 Organization and Function Are Modulated by Endoplasmic Reticulum Phospholipid Heterogeneity. Drug Metab. Dispos 44, 1859–1866. doi:10.1124/dmd.115.068981
Brouwer, M., and Lee, R. F. (2007). “Responses to Toxic Chemicals at the Molecular, Cellular, Tissue, and Organismal Level,” in The bluecrab, Callinectes sapidus. Editors V. S. Kennedy, and L. E. Cronin. College Park: Maryland Sea Grant College, 565–654.
Can, H., Chanumolu, S. K., Gonzalez-Muñoz, E., Prukudom, S., Otu, H. H., and Cibelli, J. B. (2020). Comparative Analysis of Single-Cell Transcriptomics in Human and Zebrafish Oocytes. BMC Genomics 21, 471. doi:10.1186/s12864-020-06860-z
Chamberlain, K. A.-O., and Sheng, Z. A.-O. (2019). Mechanisms for the Maintenance and Regulation of Axonal Energy Supply. J. Neurosci. Res. 97, 897–913. doi:10.1002/jnr.24411
Chatterjee, S., and Kwiterovich., P. O. (1984). Glycosphingolipids and Plasma Lipoproteins: a Review. Can. J. Biochem. Cel Biol. 62, 385–397. doi:10.1139/o84-054
Chen, L., Wang, Z., Gu, W., Zhang, X.-x., Ren, H., and Wu, B. (2020). Single-Cell Sequencing Reveals Heterogeneity Effects of Bisphenol A on Zebrafish Embryonic Development. Environ. Sci. Technol. 54, 9537–9546. doi:10.1021/acs.est.0c02428
Chen, X., Jiang, R., and Geng, Z. (2012). Cold Stress in Broiler: Global Gene Expression Analyses Suggest a Major Role of CYP Genes in Cold Responses. Mol. Biol. Rep. 39, 425–429. doi:10.1007/s11033-011-0754-x
Cheng, W., and Chen, J.-C. (2001). Effects of Intrinsic and Extrinsic Factors on the Haemocyte Profile of the Prawn, Macrobrachium Rosenbergii. Fish Shellfish Immunol. 11, 53–63. doi:10.1006/fsim.2000.0293
Colinet, H., Renault, D., and Roussel, D. (2017). Cold Acclimation Allows Drosophila Flies to Maintain Mitochondrial Functioning under Cold Stress. Insect Biochem. Mol. Biol. 80, 52–60. doi:10.1016/j.ibmb.2016.11.007
Dall, W., and Moriarty, D. J. W. (1983). “Functional Aspects of Nutrition and Digestion,” in 4 - Functional Aspects of Nutrition and Digestion. Internal Anatomy and Physiological Regulation. Editor L. H. Mantel (Academic Press), 215–261. doi:10.1016/b978-0-12-106405-1.50015-1
Des Marteaux, L. E., McKinnon, A. H., Udaka, H., Toxopeus, J., and Sinclair, B. J. (2017). Effects of Cold-Acclimation on Gene Expression in Fall Field Cricket (Gryllus pennsylvanicus) Ionoregulatory Tissues. BMC Genomics 18, 357. doi:10.1186/s12864-017-3711-9
Ellis, J. M., Li, L. O., Wu, P.-C., Koves, T. R., Ilkayeva, O., Stevens, R. D., et al. (2010). Adipose Acyl-CoA Synthetase-1 Directs Fatty Acids toward β-Oxidation and Is Required for Cold Thermogenesis. Cel Metab. 12, 53–64. doi:10.1016/j.cmet.2010.05.012
Falcon, S., and Gentleman, R. (2007). Using GOstats to Test Gene Lists for GO Term Association. Bioinformatics 23, 257–258. doi:10.1093/bioinformatics/btl567
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
Fan, L., Wang, L., and Wang, Z. (2019). Proteomic Characterization of the Hepatopancreas in the Pacific white Shrimp Litopenaeus Vannamei under Cold Stress: Revealing the Organism Homeostasis Mechanism. Fish Shellfish Immunol. 92, 438–449. doi:10.1016/j.fsi.2019.06.037
Farnsworth, D. R., Saunders, L. M., and Miller, A. C. (2020). A Single-Cell Transcriptome Atlas for Zebrafish Development. Dev. Biol. 459, 100–108. doi:10.1016/j.ydbio.2019.11.008
Farrell, J. A., Wang, Y., Riesenfeld, S. J., Shekhar, K., Regev, A., and Schier, A. F. (2018). Single-cell Reconstruction of Developmental Trajectories during Zebrafish Embryogenesis. Science 360. doi:10.1126/science.aar3131
Finak, G., McDavid, A., Yajima, M., Deng, J., Gersuk, V., Shalek, A. K., et al. (2015). MAST: a Flexible Statistical Framework for Assessing Transcriptional Changes and Characterizing Heterogeneity in Single-Cell RNA Sequencing Data. Genome Biol. 16, 278. doi:10.1186/s13059-015-0844-5
Fonville, J. M., Carter, C. L., Pizarro, L., Steven, R. T., Palmer, A. D., Griffiths, R. L., et al. (2013). Hyperspectral Visualization of Mass Spectrometry Imaging Data. Anal. Chem. 85, 1415–1423. doi:10.1021/ac302330a
Gawad, C., Koh, W., and Quake, S. R. (2016). Single-cell Genome Sequencing: Current State of the Science. Nat. Rev. Genet. 17, 175–188. doi:10.1038/nrg.2015.16
Gombos, Z., Wada, H., and Murata, N. (1994). The Recovery of Photosynthesis from Low-Temperature Photoinhibition Is Accelerated by the Unsaturation of Membrane Lipids: a Mechanism of Chilling Tolerance. Proc. Natl. Acad. Sci. 91, 8787–8791. doi:10.1073/pnas.91.19.8787
Gomez-Jimenez, S., Uglow, R. F., and Gollas-Galvan, T. (2000). The Effects of Cooling and Emersion on Total Haemocyte Count and Phenoloxidase Activity of the Spiny lobsterPanulirus Interruptus. Fish Shellfish Immunol. 10, 631–635. doi:10.1006/fsim.2000.0277
Hall, C. N., Klein-Flügge Mc Fau - Howarth, C., Howarth C Fau - Attwell, D., and Attwell, D. (2012). Oxidative Phosphorylation, Not Glycolysis, powers Presynaptic and Postsynaptic Mechanisms Underlying Brain Information Processing. J. Neurosci. 32, 8940–8951. doi:10.1523/JNEUROSCI.0026-12.2012
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
Hicks, J. M. T., Bailey, J. R., and Driedzic, W. R. (1996). Acclimation to Low Temperature Is Associated with an Increase in Long-Chain Acyl-CoA Synthetase in Rainbow trout (Oncorhynchus mykiss) Heart. Can. J. Zool. 74, 1–7. doi:10.1139/z96-001
Hopkin, S. P., and Nott, J. A. (2009). Studies on the Digestive Cycle of the Shore Crab Carcinus maenas (L.) with Special Reference to the B Cells in the Hepatopancreas. J. Mar. Biol. Ass. 60, 891–907. doi:10.1017/S0025315400041977
Horie, T., Horie, R., Chen, K., Cao, C., Nakagawa, M., Kusakabe, T. G., et al. (2018). Regulatory Cocktail for Dopaminergic Neurons in a Protovertebrate Identified by Whole-Embryo Single-Cell Transcriptomics. Genes Dev. 32, 1297–1302. doi:10.1101/gad.317669.118
Hu, K.-J., and Leung, P.-C. (2007). Food Digestion by Cathepsin L and Digestion-Related Rapid Cell Differentiation in Shrimp Hepatopancreas. Comp. Biochem. Physiol. B: Biochem. Mol. Biol. 146, 69–80. doi:10.1016/j.cbpb.2006.09.010
Hui, M., Liu, Y., Song, C., Li, Y., Shi, G., and Cui, Z. (2014). Transcriptome Changes in Eriocheir Sinensis Megalopae after Desalination Provide Insights into Osmoregulation and Stress Adaption in Larvae. PLoS One 9, e114187. doi:10.1371/journal.pone.0114187
Ito, M., Ito, K., Ohta, K., Hano, T., Onduka, T., and Mochida, K. (2016). Transcription of a Novel P450 Gene Varies with Some Factors (Pollutant Exposure, Temperature, Time, and Body Region) in a marine Oligochaete ( Thalassodrilides sp.). Mar. Pollut. Bull. 109, 344–349. doi:10.1016/j.marpolbul.2016.05.055
James, M. O., and Boyle, S. M. (1998). Cytochromes P450 in Crustacea. Comp. Biochem. Physiol. C: Pharmacol. Toxicol. Endocrinol. 121, 157–172. doi:10.1016/s0742-8413(98)10036-1
Kalisky, T., Blainey, P., and Quake, S. R. (2011). Genomic Analysis at the Single-Cell Level. Annu. Rev. Genet. 45, 431–445. doi:10.1146/annurev-genet-102209-163607
Kanehisa, M., and Goto, S. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28, 27–30. doi:10.1093/nar/28.1.27
Keeran, W. S., and Lee, R. F. (1987). The Purification and Characterization of Glutathione S-Transferase from the Hepatopancreas of the Blue Crab, Callinectes sapidus. Arch. Biochem. Biophys. 255, 233–243. doi:10.1016/0003-9861(87)90390-0
Kuz'mina, V. V., Gavrovskaia, L. K., and Ryzhova, O. V. (2010). Taurine. Effect on Exotrophia and Metabolism in Mammals and Fish. Zh Evol. Biokhim Fiziol 46, 17–23.
Lange, C., Rost, F., Machate, A., Reinhardt, S., Lesche, M., Weber, A., et al. 2020). Single Cell Sequencing of Radial Glia Progeny Reveals Diversity of Newborn Neurons in the Adult Zebrafish Brain. 147, 1855951, doi:10.1242/dev.185595
Le Moullac, G., and Haffner, P. (2000). Environmental Factors Affecting Immune Responses in Crustacea. Aquaculture 191, 121–131. doi:10.1016/s0044-8486(00)00422-1
Lee, J. I., Kim, S. S., and Kang, D. H. (2019). Susceptibility of Escherichia coli O157:H7 Grown at Low Temperatures to the Krypton-Chlorine Excilamp. Sci. Rep. 9, 563. doi:10.1038/s41598-018-37060-1
Lehnert, S. A., and Johnson, S. E. (2002). Expression of Hemocyanin and Digestive Enzyme Messenger RNAs in the Hepatopancreas of the Black Tiger Shrimp Penaeus monodon. Comp. Biochem. Physiol. Part B: Biochem. Mol. Biol. 133, 163–171. doi:10.1016/s1096-4959(02)00123-9
Li, D., Nie, H., Jahan, K., and Yan, X. (2020). Expression Analyses of C-type Lectins (CTLs) in Manila Clam under Cold Stress Provide Insights for its Potential Function in Cold Resistance of Ruditapes Philippinarum. Comp. Biochem. Physiol. C: Toxicol. Pharmacol. 230, 108708. doi:10.1016/j.cbpc.2020.108708
Lim, S., Kwon, M., Joung, E.-J., Shin, T., Oh, C.-W., Choi, J., et al. (2018). Meroterpenoid-Rich Fraction of the Ethanolic Extract from Sargassum Serratifolium Suppressed Oxidative Stress Induced by Tert-Butyl Hydroperoxide in HepG2 Cells. Mar. Drugs 16, 374. doi:10.3390/md16100374
Lu, W., and Fuchou, T. (2014). Recent Progress in Single-Cell RNA-Seq Analysis. Yi Chuan 36, 1069–1076. doi:10.3724/sp.j.1005.2014.1069
McDavid, A., Finak, G., Chattopadyay, P. K., Dominguez, M., Lamoreaux, L., Ma, S. S., et al. (2013). Data Exploration, Quality Control and Testing in Single-Cell qPCR-Based Gene Expression Experiments. Bioinformatics 29, 461–467. doi:10.1093/bioinformatics/bts714
Möhrlen, F., Baus, S., Gruber, A., Rackwitz, H.-R., Schnölzer, M., Vogt, G., et al. (2001). Activation of Pro-astacin. Eur. J. Biochem. 268, 2540–2546. doi:10.1046/j.1432-1327.2001.02136.x
Moris, N., Pina, C., and Arias, A. M. (2016). Transition States and Cell Fate Decisions in Epigenetic Landscapes. Nat. Rev. Genet. 17, 693–703. doi:10.1038/nrg.2016.98
Nakai, Y., Nakahira, Y., Sumida, H., Takebayashi, K., Nagasawa, Y., Yamasaki, K., et al. (2013). Vascular Plant One-Zinc-finger Protein 1/2 Transcription Factors Regulate Abiotic and Biotic Stress Responses in Arabidopsis. Plant J. 73, 761–775. doi:10.1111/tpj.12069
Nanji, A. A., and Hiller-Sturmhöfel, S. (1997). Hiller-Sturmhöfel Apoptosis and Necrosis: Two Types of Cell Death in Alcoholic Liver Disease.21, 325–330.
Nebert, D. W., Dalton Tp Fau - Okey, A. B., and Okey Ab Fau - Gonzalez, F. J. (2004). Gonzalez Role of Aryl Hydrocarbon Receptor-Mediated Induction of the CYP1 Enzymes in Environmental Toxicity and Cancer. J. Biol. Chem. 279, 23847–23850. doi:10.1074/jbc.R400004200
Niu, J., Huang, Y., Liu, X., Zhang, Z., Tang, J., Wang, B., et al. (2020). Single-cell RNA-Seq Reveals Different Subsets of Non-specific Cytotoxic Cells in Teleost. Genomics 112, 5170–5179. doi:10.1016/j.ygeno.2020.09.031
Patey, C. P., and Driedzic, W. R. (1997). Cold Acclimation Increases Activities of Mitochondrial Long-Chain Acyl-CoA Synthetase and Carnitine Acyl-CoA Transferase I in Heart of Rainbow trout (Oncorhynchus mykiss). Can. J. Zool. 75, 324–331. doi:10.1139/z97-041
Peng, J.-X., He, P.-P., Wei, P.-Y., Zhang, B., Zhao, Y.-Z., Li, Q.-Y., et al. (2018). Proteomic Responses under Cold Stress Reveal Unique Cold Tolerance Mechanisms in the Pacific White Shrimp (Litopenaeus Vannamei). Front. Physiol. 9, 1399. doi:10.3389/fphys.2018.01399
Perepechaeva, M. L., and Grishanova, A. (2013). Comparison of the Effects of Cold on Cytochromes Р450 1А in the Rat Lungs versus Liver. Ab 1, 7–12. doi:10.11648/j.ab.20130102.11
Ponce-Palafox, J., Martinez-Palacios, C. A., and Ross, L. G. (1997). The Effects of Salinity and Temperature on the Growth and Survival Rates of Juvenile white Shrimp, Penaeus Vannamei, Boone, 1931. Aquaculture 157, 107–115. doi:10.1016/S0044-8486(97)00148-8
Qian, X., Lai, Y., and Zhu, F. (2019). Molecular Characterization of Carboxypeptidase B-like (CPB) in Scylla Paramamosain and its Role in white Spot Syndrome Virus and Vibrio Alginolyticus Infection. Fish Shellfish Immunol. 94, 434–446. doi:10.1016/j.fsi.2019.09.036
Qiu, X., Mao, Q., Tang, Y., Wang, L., Chawla, R., Pliner, H. A., et al. (2017). Reversed Graph Embedding Resolves Complex Single-Cell Trajectories. Nat. Methods 14, 979–982. doi:10.1038/nmeth.4402
Qu, J. Y., Xiao, Y. T., Zhang, Y. Y., Xie, H. T., and Zhang, M. C. (2021). Hedgehog Signaling Pathway Regulates the Proliferation and Differentiation of Rat Meibomian Gland Epithelial Cells. Invest. Ophthalmol. Vis. Sci. 62–33. doi:10.1167/iovs.62.2.33
Sanderson, J. T. (2006). The Steroid Hormone Biosynthesis Pathway as a Target for Endocrine-Disrupting Chemicals. Toxicol. Sci. 94, 3–21. doi:10.1093/toxsci/kfl051
Schwerin, S., Zeis, B., Lamkemeyer, T., Paul, R. J., Koch, M., Madlung, J., et al. (2009). Acclimatory Responses of the Daphnia pulex Proteome to Environmental Changes. II. Chronic Exposure to Different Temperatures (10 and 20°C) Mainly Affects Protein Metabolism. BMC Physiol. 9, 8. doi:10.1186/1472-6793-9-8
Sharma, S., Wang, W., and Stolfi, A. (2019). Single-cell Transcriptome Profiling of the Ciona Larval Brain. Dev. Biol. 448, 226–236. doi:10.1016/j.ydbio.2018.09.023
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
Sousa, L. G., and Petriella, A. M. (2001). Changes in the Hepatopancreas Histology of Palaemonetes Argentinus (Crustacea, Caridea) during Moult. Biocell 25, 275–281.
Sousa, L. G., and Petriella, A. M. (2000). Histology of the Hepatopancreas of the Freshwater Prawn Palaemonetes Argentinus (Crustacea, Caridea). Biocell 24, 189–195.
Squire, T. L., Lowe, M. E., Bauer, V. W., and Andrews, M. T. (2003). Pancreatic Triacylglycerol Lipase in a Hibernating Mammal. II. Cold-Adapted Function and Differential Expression. Physiol. Genomics 16, 131–140. doi:10.1152/physiolgenomics.00168.2002
Steenbergen, J. F., Steenbergen, S. M., and Schapiro, H. C. (1978). Effects of Temperature on Phagocytosis in Homarus americanus. Aquaculture 14, 23–30. doi:10.1016/0044-8486(78)90137-0
Stöcker, W., Raeder, U., Bijlholt, M. M. C., Wichertjes, T., van Bruggen, E. F. J., and Markl, J. (1988). The Quaternary Structure of Four Crustacean Two-Hexameric Hemocyanins: Immunocorrelation, Stoichiometry, Reassembly and Topology of Individual Subunits. J. Comp. Physiol. B 158, 271–289. doi:10.1007/BF00695326
Stuart, T., Butler, A., Hoffman, P., Hafemeister, C., Papalexi, E., Mauck, W. M., et al. (2019). Comprehensive Integration of Single-Cell Data. Cell 177, 1888–1902. doi:10.1016/j.cell.2019.05.031
Tanay, A., and Regev, A. (2017). Scaling Single-Cell Genomics from Phenomenology to Mechanism. Nature 541, 331–338. doi:10.1038/nature21350
Tian, Z., Kang, X., and Mu, S. (2012). The Molt Stages and the Hepatopancreas Contents of Lipids, Glycogen and Selected Inorganic Elements during the Molt Cycle of the Chinese Mitten Crab Eriocheir Sinensis. Fish. Sci. 78, 67–74. doi:10.1007/s12562-011-0426-8
Toullec, J. Y., Chikhi, M., and van Wormhoudt, A. (1992). In Vitro protein Synthesis and α Amylase Activity in F Cells from Hepatopancreas ofPalaemon Serratus (Crustacea; Decapoda). Experientia 48, 272–277. doi:10.1007/bf01930474
Trapnell, C., Cacchiarelli, D., Grimsby, J., Pokharel, P., Li, S., Morse, M., et al. (2014). The Dynamics and Regulators of Cell Fate Decisions Are Revealed by Pseudotemporal Ordering of Single Cells. Nat. Biotechnol. 32, 381–386. doi:10.1038/nbt.2859
Trapnell, C. (2015). Defining Cell Types and States with Single-Cell Genomics. Genome Res. 25, 1491–1498. doi:10.1101/gr.190595.115
Tryselius, Y., and Hultmark, D. (1997). Cysteine Proteinase 1 (CP1), a Cathepsin L-like Enzyme Expressed in the Drosophila melanogaster Haemocyte Cell Line Mbn-2. Insect Mol. Biol. 6, 173–181. doi:10.1111/j.1365-2583.1997.tb00085.x
Tsukamoto, H., Zhu, N.-L., Wang, J., Asahina, K., and Machida, K. (2012). Morphogens and Hepatic Stellate Cell Fate Regulation in Chronic Liver Disease. J. Gastroenterol. Hepatol. 27 (Suppl. 2), 94–98. doi:10.1111/j.1440-1746.2011.07022.x
Tuo, B., Li, C., Peng, L., Ye, M., Liu, W., Zhong, X., et al. (2014). Analysis of Differentially Expressed Genes in Cold-Exposed Mice to Investigate the Potential Causes of Cold-Induced Hypertension. Exp. Ther. Med. 8, 110–114. doi:10.3892/etm.2014.1703
Van wormhoudt, A., Bourreau, G., and Le Moullac, G. (1995). Amylase Polymorphism in Crustacea Decapoda: Electrophoretic and Immunological Studies. Biochem. Syst. Ecol. 23, 139–149. doi:10.1016/0305-1978(94)00090-4
Vergara-Amado, J., Silva, A. X., Manzi, C., Nespolo, R. F., and Cárdenas, L. (2017). Differential Expression of Stress Candidate Genes for thermal Tolerance in the Sea Urchin Loxechinus Albus. J. Therm. Biol. 68, 104–109. doi:10.1016/j.jtherbio.2017.03.009
Verri, T., Mandal, A., Zilli, L., Bossa, D., Mandal, P. K., Ingrosso, L., et al. (2001). D-glucose Transport in Decapod Crustacean Hepatopancreas. Comp. Biochem. Physiol. A: Mol. Integr. Physiol. 130, 585–606. doi:10.1016/s1095-6433(01)00434-2
Vogt, G. (2019). Functional Cytology of the Hepatopancreas of Decapod Crustaceans. J. Morphol. 280, 1405–1444. doi:10.1002/jmor.21040
Vogt, G., Stcker, W., Storch, V., and Zwilling, R. (1989). Biosynthesis of Astacus Protease, a Digestive Enzyme from Crayfish. Histochemistry 91, 373–381. doi:10.1007/bf00493824
Wagner, D. E., Weinreb, C., Collins, Z. M., Briggs, J. A., Megason, S. G., and Klein, A. M. (2018). Single-cell Mapping of Gene Expression Landscapes and Lineage in the Zebrafish Embryo. Science 360, 981–987. doi:10.1126/science.aar4362
Wang, Y.-C., Chang, P.-S., and Chen, H.-Y. (2007). Tissue Expressions of Nine Genes Important to Immune Defence of the Pacific white Shrimp Litopenaeus Vannamei. Fish Shellfish Immunol. 23, 1161–1177. doi:10.1016/j.fsi.2007.04.004
Wang, Y. L., Guo, Y., Zang Y Fau, G. C., Yaney Gc Fau, G., Fau Getty-Kaushik, L. V. G., Getty-Kaushik, L., et al. (2004). Acyl Coenzyme a Synthetase Regulation: Putative Role in Long-Chain Acyl Coenzyme a Partitioning. Obes. Res. 12, 1781–1788. doi:10.1038/oby.2004.221
Wang, Y., and Navin, N. E. (2015). Advances and Applications of Single-Cell Sequencing Technologies. Mol. Cel 58, 598–609. doi:10.1016/j.molcel.2015.05.005
Wen, L., and Tang, F. (2016). Single-cell Sequencing in Stem Cell Biology. Genome Biol. 17, 71. doi:10.1186/s13059-016-0941-0
Williams, C. M., Rocca, J. R., Edison, A. S., Allison, D. B., Morgan, T. J., and Hahn, D. A. (2018). Cold Adaptation Does Not Alter ATP Homeostasis during Cold Exposure in Drosophila melanogaster. Integr. Zoolog. 13, 471–481. doi:10.1111/1749-4877.12326
Wutka, A., Palagani, V., Barat, S., Chen, X., El Khatib, M., Götze, J., et al. (2014). Capsaicin Treatment Attenuates Cholangiocarcinoma Carcinogenesis, PLoS One 9, e95605. doi:10.1371/journal.pone.0095605
Xu, H., Zhang, D. L., Yu, D. H., Lv, C. H., Luo, H. Y., and Wang, Z. Y. (2015). Molecular Cloning and Expression Analysis of Scd1 Gene from Large Yellow Croaker Larimichthys Crocea under Cold Stress. Gene 568, 100–108. doi:10.1016/j.gene.2015.05.027
Yepiz-Plascencia, G., Galván, T. G., Vargas-Albores, F., and García-Bañuelos, M. (2000). Synthesis of Hemolymph High-Density Lipoprotein β-Glucan Binding Protein by Penaeus Vannamei Shrimp Hepatopancreas. Mar. Biotechnol. 2, 485–492. doi:10.1007/s101260000030
Zhang, G., Fang, X., Guo, X., Li, L., Luo, R., Xu, F., et al. (2012). The Oyster Genome Reveals Stress Adaptation and Complexity of Shell Formation. Nature 490, 49–54. doi:10.1038/nature11413
Zhang, H., Jiang, C., Ren, J., Dong, J., Shi, X., Zhao, X., et al. (2020). An Advanced Lipid Metabolism System Revealed by Transcriptomic and Lipidomic Analyses Plays a Central Role in Peanut Cold Tolerance. Front. Plant Sci. 11, 1110. doi:10.3389/fpls.2020.01110
Zhang, X., Liu, S., and Takano, T. (2008). Two Cysteine Proteinase Inhibitors from Arabidopsis thaliana, AtCYSa and AtCYSb, Increasing the Salt, Drought, Oxidation and Cold Tolerance. Plant Mol. Biol. 68, 131–143. doi:10.1007/s11103-008-9357-x
Zhang, X., Yuan, J., Sun, Y., Li, S., Gao, Y., Yu, Y., et al. (2019). Penaeid Shrimp Genome Provides Insights into Benthic Adaptation and Frequent Molting. Nat. Commun. 10, 356. doi:10.1038/s41467-018-08197-4
Keywords: single-cell RNA sequencing, Litopenaeus vannamei, cold tolerance, mechanisms, pacific white shrimp
Citation: Zhu W, Yang C, Chen X, Liu Q, Li Q, Peng M, Wang H, Chen X, Yang Q, Liao Z, Li M, Pan C, Feng P, Zeng D and Zhao Y (2022) Single-Cell Ribonucleic Acid Sequencing Clarifies Cold Tolerance Mechanisms in the Pacific White Shrimp (Litopenaeus Vannamei). Front. Genet. 12:792172. doi: 10.3389/fgene.2021.792172
Received: 09 October 2021; Accepted: 25 November 2021;
Published: 12 January 2022.
Edited by:
Kui Li, Institute of Animal Sciences (CAAS), ChinaReviewed by:
Xianyun Ren, Yellow Sea Fisheries Research Institute (CAFS), ChinaChenghua Li, Ningbo University, China
Copyright © 2022 Zhu, Yang, Chen, Liu, Li, Peng, Wang, Chen, Yang, Liao, Li, Pan, Feng, Zeng and Zhao. 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.
*Correspondence: Yongzhen Zhao, yongzhenzhao@hotmail.com; Digang Zeng, zengdigang@126.com
†These authors have contributed equally to this work