Yaoyao Zhan
Yaqing Chang- Key Laboratory of Mariculture & Stock Enhancement in North China’s Sea, Ministry of Agriculture and Rural Affairs, Dalian Ocean University, Dalian, China
To explore the impact of chronic heat stress on commercial echinoderms, the present study assessed the effects of chronic high temperature on the growth, survival, feeding, and differential gene expression in the sea urchin Strongylocentrotus intermedius cultured in northern Yellow Sea in China. One suitable seawater condition (20°C) and one laboratory-controlled high temperature condition (25°C) were set up. After 28 days incubation, our results showed that: (1) The specific growth, survival, and ingestion rates of S. intermedius reared under high temperature (25°C) decreased compared to those reared under optimal temperature (20°C) conditions; (2) comparative transcriptome analysis identified 2,125 differentially expressed genes (DEGs) in S. intermedius reared under high temperature (25°C) compared to those subjected to optimal temperature condition (20°C), which included 1,015 upregulated and 1,100 downregulated genes. The accuracy of the transcriptome profiles was verified by quantitative real-time PCR (qRT-PCR). Further Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyses revealed that these DEGs mainly enriched the functional categories of ribosome, protein processing in endoplasmic reticulum, and prion diseases. A total of 732 temperature-induced expressed genes, such as ATP5, heat shock protein 70, and heat shock protein 90, were identified as candidates that were closely correlated with heat resistance in S. intermedius. Differentially expressed transcription factors (TFs), such as AP-1, Fos, CREB, and ZNF, were also identified as potential regulators that regulate the molecular network that was associated with responses to heat stress in sea urchins. Observations in the present study provide additional information that improves our understanding of the molecular mechanism of temperate echinoid species in response to heat stress, as well as theoretical basis for the molecular-assisted breeding of heat-resistant sea urchins.
Introduction
Seawater temperature has been proven to be a major environmental factor affecting echinoderms from the biomacromolecular to ecological levels. With climate-induced ocean warming, extensive efforts have been made in studying the impact of elevated seawater temperatures on echinoderms. Laboratory-based studies have demonstrated that increased seawater temperature affects early development, survival, growth, metabolism, immunity, behavior, and gene expression profiles in echinoderms. For example, with the elevation of near-future seawater temperature, fertilization and early development in the sea urchin Heliocidaris erythrogramma would be compromised (Byrne et al., 2009). Increased seawater temperatures can affect and reduce both the specific growth rate (SGR) and contents of highly unsaturated fatty acids (HUFAs) in juvenile sea cucumber (Yu et al., 2016). It has been demonstrated that the parental effect of long acclimatization could increase thermal tolerance in juvenile sea cucumber Apostichopus japonicus (Wang et al., 2015). The existence of species-specific innate immune response variations was investigated in the tropical subtidal sea urchin Lytechinus variegatus and the intertidal sea urchin Echinometra lucunter while coping with rising sea temperatures (Branco et al., 2013). In addition, negative effects of elevated seawater temperature on covering and righting behaviors were observed in the sea urchins L. variegatus and Strongylocentrotus intermedius (Brothers and Mcclintock, 2015; Zhang et al., 2017). Comparative transcriptome study indicated alterations in gene expression profiles under mild, chronic increases in temperature stress in embryos of the sea urchin Strongylocentrotus purpuratus (Runcie et al., 2012). A recent study also showed that the response of juvenile sea urchin Loxechinus albus to acute increases in sea temperature is an integrated differential gene regulatory network that includes heat-shock, membrane potential, and detoxification (Vergara-Amado et al., 2017).
The temperate edible sea urchin S. intermedius is naturally distributed along the intertidal and subtidal rocky bottom of Hokkaido, Japan, the Korean Peninsula, and Russian Far East (Chang et al., 2004; Lawrence, 2013). This species has an average lifespan of 8–10 years, and the sexual maturity age is 1.5–2 years. The thermal tolerance of this species is from -1 to 23°C (Chang et al., 2004), and the suitable sea temperature range for the growth of this species is 15–20°C. In 1989, this species was introduced from Japan to north China by the Dalian Ocean University, and artificial breeding was subsequently performed. To date, S. intermedius has been the predominant commercial valuable sea urchin species that has been widely cultivated along the coastal areas of the north Yellow Sea in China (Chang et al., 2004). Due to global ocean warming, sea water temperatures in the north Yellow Sea in China have often been higher than 25°C (the lethal limit of S. intermedius) in the summer in the past few years (Zeng et al., 2006), resulting in the massive death of cultured S. intermedius. The sustainable development of S. intermedius farming and industry, therefore, is under serious threat. Our previous study demonstrated the existence of genotype by temperature interactions (GEI) in the survival rate (SR) in the selection of S. intermedius (Chang et al., 2016); however, the response to heat-stress, especially the corresponding gene expression mechanism in S. intermedius remain unclear.
In the present study, we investigated the impact of high water temperatures on the growth, survival, and feeding of S. intermedius. Then, we identified candidate genes that were closely correlated to heat tolerance in S. intermedius by comparative transcriptome analysis between suitable (20°C, as control) and high temperature (25°C) seawater conditions in S. intermedius. A gene regulatory network related to heat-tolerance was also predicted by setting up relationships between candidate genes and differential expressed transcription factors (TFs). The findings of this study enrich our knowledge of the molecular responses of sea urchins to heat stress, as well as provide candidate genes that can serve as molecular markers that could be potentially used in the selection of heat tolerance-resistant breeding of S. intermedius.
Materials and Methods
Sea Urchins and Treatments
A total of 360 S. intermedius (average test diameter: 10 ± 0.1 mm) were transported from Dalian Haibao Fisheries Company to the Key Laboratory of Mariculture & Stock Enhancement in the North China’s Sea, Ministry of Agriculture and Rural Affairs at the Dalian Ocean University in August 2015. All of the sea urchins were kept in ∼60-L recirculating sea water tanks; each tank was fitted with an automatic temperature control and monitoring system (Dalian Huixin Titanium Equipment Development Co., Ltd., Liaoning, China). Seawater was sand filtered and continuously aerated. The animals were kept under natural light. All of the specimens were fed kelp (Laminaria japonica). Sea urchins were acclimated to default laboratory conditions [18 ± 0.5°C and 31.22 ± 0.14 (practical salinity units) PSU] for 1 week prior to experimentation. The experiments were conducted between November 2015 and January 2016.
All of the sea urchins were dried with a paper towel and weighed on a digital balance (0.01 g sensitivity; AL204; Mettler Toledo, Shanghai, China) to obtain initial mass (W1). We then randomly divided the sea urchins into three groups of 60 specimens each (three replicates for each temperature). Each group was housed in a separate tank. To reach the desired temperature (20 and 25°C), we removed half of the seawater from each tank every day, and replaced it with seawater at a different temperature. We changed the temperature of the new seawater such that the temperature of the entire tank did not increase by more than 1°C per day; this was based on a previous study on S. intermedius (Lawrence et al., 2009) and on field survey data of the coastal waters of the Yellow Sea (Zhang et al., 2016).
We monitored the temperature in each tank using an automatic temperature control and with a water quality monitor (A329 Portable Meter; Thermo Scientific Orion Star, Beijing, China).
Growth, Survival, and Ingestion in Each Treatment
Before the experiment, S. intermedius individuals in each treatment were dried with a paper towel and weighed on a digital balance (0.01 g sensitivity; AL204; Mettler Toledo, Shanghai, China) to obtain initial average mass (W0). The specific growth rate (SGR), survival rate (S), average food consumption (FC) of individual, and daily feeding rate (FR) were calculated using the following formulae:
where Wt is the average body weight (g) of live S. intermedius on day t; t is the duration of experiment; N0 is the initial number of live S. intermedius; Nt is the number of live S. intermedius on day t; TBt is the bait supplied on days t; and RBt is the total amount of remaining bait on day t (Qin et al., 2011; Chang et al., 2016).
Sample Collection for RNA-Seq
As for individual test diameter of around 10 mm, most of the sea urchins sampled in this study were too small to develop their gonads. In addition, it is generally difficult to dissect or obtain enough tissues such as tube feet, coelomocytes, and perioral membranes for RNA-seq library construction and subsequent validation.
The intestines are important organs for nutrient intake and stress defense in sea urchins, and these are the only sources for tissues that can be sampled under a dissection microscope in this study. Therefore, we opted to utilize the intestines for perform transcriptome analysis in this study. At the end of the experimental period, the number of living S. intermedius individuals in three high temperature group (25°C) replicate tanks was 40 (replicate 1#), 31 (replicate 2#), and 34 (replicate 3#). To ensure that there are three independent samples for the transcriptome validation experiment, we randomly selected 20 S. intermedius specimens from each replicate tank for intestinal RNA extraction. For RNA sequencing library construction, the intestines of 20 living S. intermedius specimens from replicate 1# of high temperature group (25°C) were carefully removed and pooled (designated as Si_TT2_1), and 20 intestines of living S. intermedius specimens from replicate 2# and replicate 3# (10 of each) of high temperature group (25°C) were carefully removed and pooled (designated as Si_TT2_2). Two optimal temperature (20°C) sample pools (as control) were constructed using the procedure employed for high temperature (25°C) sample pool construction. We named two optimal temperature (20°C) sample pools as Si_TT0_1 and Si_TT0_2. All of the pooled samples were stored at -80°C until RNA extraction.
RNA Extraction and Sequencing
Total RNA was extracted from each pooled sample by using TRIzol (Ambion, United States) following the manufacturer’s instructions. Total RNA quantity and integrity were assessed by 1% agarose gel electrophoresis and the RNA Nano 6000 assay kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, United States).
A total amount of 2 μg high quality RNA per pooled sample was used for RNA sample preparation. All high-quality RNA samples were sent to BGI Co., Ltd. (Beijing, China). Sequencing libraries were constructed using the NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, United States) according to the manufacturer’s recommendations, and index codes were added to sequences to distinguish one sample from another. The Agilent Bioanalyzer 2100 system was employed to assess the quality of each RNA library. The index-coded samples were clustered on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumina) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina Hiseq platform and paired-end reads (Hiseq 4000, 101PE).
Transcriptome Assembly and Annotation
Clean data (clean reads) were obtained from the raw data (raw reads) by removing adaptors, reads with poly-N, and low-quality reads. All of the clean reads were then submitted to the National Center for Biotechnology Information (NCBI) Short Read Archive (SRA) Sequence Database (Accession Number PRJNA508827). The percentage of bases with a Phred value > 20 (Q20), the percentage of bases with a Phred value > 30 (Q30), and the content of base G and C (GC-content) were calculated. The high-quality clean data were used for subsequent analyses. Clean reads were assembled into transcriptome as reference using Trinity (Grabherr et al., 2011) with min_kmer_cov set to 2 by default, and all of the other parameters set to default. The clean data were mapped back onto the assembled transcriptome, and read count for each gene was obtained from the mapping results by RSEM (v1.2.12). BUSCO v.3.0.2 was used to assess the completeness of the gene assembly (Simão et al., 2015).
Transcriptome annotation was performed using Basic Local Alignment Search Tool (BLAST) searches against the NCBI non-redundant (Nr) databases, NCBI nucleotide sequences (Nt), Swiss-Prot, InterPro, Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Swiss-Prot. We employed Blast2GO with NR annotation for Gene Ontology (GO) annotation, and InterProScan5 for InterPro annotation.
Single Nucleotide Polymorphism (SNP) and Simple Sequence Repeat (SSR) Identification
Single nucleotide polymorphism (SNPs) and SSRs in the transcriptome level were identified using GATK3 software (v3.4) (Mckenna et al., 2010) and MISA (microsatellite)1, respectively. The parameters for SNP identification were an MQ < 20.0 and QD < 2.0. SSR identification criteria in the MISA script were 1–12, 2–6, 3–5, 4–5, 5–4, and 6–4.
Differentially Expressed Gene (DEG) Analysis
Gene expression levels were calculated as previously described (Li and Dewey, 2011) with RSEM (v1.2.12). Differential expression analysis between optimal temperature (20°C; Si_TT0) and high temperature (25°C; Si_TT2) was performed using the NOISeq R package (v3.1). NOISeq provides statistical routines for determining differential expression in digital gene expression data using a model based on a noise distribution model (Tarazona et al., 2011). The software information: fold-change ≥ 2.0 and probability ≥ 0.8. Since genes with an adjusted | fold-change|≥ 2.0 and probability ≥ 0.8 found by NOISeq were assigned as differentially expressed.
GO and KEGG Pathway Enrichment Analyses
Differentially expressed genes were classified based on GO and KEGG functional annotation, GO, and pathway functional enrichment was performed using phyper, as implemented in the R package (v3.1). P-values were calculated using the hypergeometric test:
To ensure relatively precise results, we calculated the false discovery rate (FDR) for each p-value. In general, the terms in which FDR was not larger than 0.001 were defined as significant. Enriched cluster analysis of candidate DEGs was performed using the R package (v3.1).
qRT-PCR Validation
Annotated DEGs were validated using quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). Twelve DEG candidates were randomly selected (including eight upregulated and four downregulated DEGs) from 2,125 DEGs. The DEG input RNA was used as template for cDNA synthesis. cDNA was synthesized using PrimeScriptTM RT reagent Kit (TaKaRa, Japan). The cytochrome b (Cytb) gene was used as internal control (Yang et al., 2010). Primers used in qRT-PCR analyses were designed by Primer Premier 5.0 (Table 1). qRT-PCR was performed in a total volume of 16 μL, which consisted of 2 μL of the cDNA template, 8 μL of 2× SYBR Green Master mix (TaKaRa, Japan), 0.3 μL of ROX reference dye II, 4.5 μL of PCR-grade water, and 0.6 μL (10 mM) of each primer. The running program was set as follows: 95°C for 30 s; followed by 40 cycles of 95°C for 5 s and annealing temperature 56°C for 32 s. At the end of reaction, PCR melting curve analysis was conducted to confirm single PCR products. The relative expression level of each candidate DEG was determined using the comparative 2-ΔΔCt method (Livak and Schmittgen, 2001). The concrete formula was as follows:
Table 1. Primers used in verification of RNA-Seq results by qRT-PCR.
Protein–Protein Interaction (PPI) Analysis
We used Blastx v2.2.28 (Zhang et al., 2000) with an e-value cutoff of 1e-10 to align the S. intermedius DEG sequences with the protein sequences from the sea urchin S. purpuratus. A PPI network was built using STRING2 (Szklarczyk et al., 2011) based on the PPI network of S. purpuratus. We used Cytoscape v3.5.1 (Shannon et al., 2003) to visualize the PPI network.
Data Analysis
All of the data were expressed as the mean ± standard deviation (SD). All of the statistical analyses were performed with SPSS 16.0 (IBM, Shanghai, China). We first confirmed that our data were normally distributed and homogeneous with the Shapiro–Wilk test and with Levene’s test. We then compared differences in survival rate, SGR among treatments with one-way ANOVA (factor: temperature). We considered p < 0.05 as statistically significant and p < 0.01 as statistically extremely significant. Significant differences between pairs of treatments were identified with Duncan’s multiple range tests.
Results
The Impact of High Seawater Temperature on Growth, Survival, and Feeding of S. intermedius
During the 30-day incubation, differences in S. intermedius growth, survival, and feeding among all of the treatments were compared and analyzed statistically. S. intermedius reared at both 20 and 25°C exhibited an increase in SGR (Figure 1A), but the SGRs of S. intermedius reared at 25°C were significantly lower compared to those reared at 20°C. During the first 14 days of incubation, no significant difference in survival was observed between 20 and 25°C treatments, whereas survival rates decreased after 14 days of incubation at 25°C (Figure 1B). Among all of the experimental groups, food consumption increased in a time-dependent manner, whereas that in S. intermedius reared at 25°C was relatively lower compared to those reared at 20°C (Figure 1C). In addition, significantly reduced feeding rates were observed at 25°C compared to 20°C (Figure 1D).
Figure 1. The effect of high temperature stress on SGR, SR, FC, and FR in S. intermedius. (A) The specific growth rates of S. intermedius cultured at 20 and 25°C. (B) The survival rates of S. intermedius between 20°C (control) and 25°C. (C) S. intermedius food consumption at 20 and 25°C. (D) S. intermedius feeding rates at 20 and 25°C. Values are expressed as the mean ± SD (n = 3). ∗Significant differences at p < 0.05 vs. control; ∗∗extremely significant differences at p < 0.01 vs. control.
RNA Sequencing, Transcriptome Assembly, and Annotation
Four RNA sequencing (RNA-Seq) libraries (Si_TT0_1, Si_TT0_2, Si_TT2_1, and Si_TT2_2) from the intestines of S. intermedius cultured at 20 and 25°C were constructed and subsequently sequenced on an Illumina Hiseq 4000 platform. Approximately 45.24–45.25 million raw reads were obtained from each pooled sample (Figure 2A). After trimming, 44.61–44.73 million clean reads were obtained from each sample (Figure 2A). The Q20 range of all of the samples was 96.85–97.02%, and the Q30 range of all samples was 92.56–92.88%. The average GC content was 37.16 ± 0.14% (Figure 2A).
Figure 2. Summary statistics of the S. intermedius transcriptomes. (A) Summary statistics of the RNA-Seq data and transcriptome assembly results. (B) Pearson correlation coefficients between samples used in the current study. (C) Venn diagram showing NR, COG, KEGG, Swiss-Prot, and InterPro.
The transcriptome of S. intermedius was de novo assembled using the Trinity software with min_kmer_cov set to 2 by default, and all of the other parameters set to default. Completeness of the de novo assembly was assessed with BUSCO using a eukaryotic database3 of 303 genes. The assembled genes were deemed 82.5% complete by BUSCO (61.4% as single genes and 21.1% as duplicated genes). In the intestines of S. intermedius cultured at 20°C (Si_TT0), a total of 188,430 transcripts were obtained, with an average length of 636 bp and N50 length of 1,089 bp, and approximately 102,528 unigenes with a mean length of 721 bp and N50 length of 1,179 bp were also generated (Figure 2A). As for S. intermedius cultured at 25°C (Si_TT2), the intestine transcriptome analysis indicated that a total of 218,307 transcripts were obtained, with an average length of 678.5 bp and N50 length of 1201.5 bp, and approximately 113,170 unigenes with a mean length of 777 bp and N50 length of 1313.5 bp were generated as well (Figure 2A). In summary, approximately 46.33% of the transcripts in Si_TT0 and approximately 53.67% of the transcripts in Si_TT2 were successfully mapped back to the de novo transcriptome assembly, respectively (Figure 2A). Pearson’s correlation coefficients of FPKM distribution among the two biological replicates indicated the reproducibility of RNA-seq data (see Figure 2B). Annotation of assembled transcriptome was performed for comprehensive functional annotation of each unigene. After alignment, 65,349 unigenes were annotated in at least one of seven databases (Figure 2C).
Identification of SNPs and SSRs
A total of 890,154 SNPs and 20,239 SSRs were identified from the assembled transcriptome. In SNPs, percent transition (Ts) and percent transversion (Tv) were 58.5 and 41.5%, respectively (Figure 3A). Ts was higher than Tv among the four libraries, and the Ts/Tv ratio was 117:83. The most common transitions were A–G (262,240, 29.46%) and C–T (258,526, 29.04%), and the predominant transversion type was A–T (128,089, 14.39%). The AG/CT repeat was the most abundant type of motif, the proportion gradient of repeat motifs in decreasing order was as follows: dimers (8,754, 43.25%) > monomers (5,573, 27.54%) > trimers (4,925, 24.33%) > pentamers (451, 2.22%) > quadmers (380, 1.88%) > hexamers (156, 0.77%) (Figure 3B).
Figure 3. Frequency of identified (A) single nucleotide polymorphisms and (B) simple sequence repeats among the transcriptome libraries Si_TT0 (n = 2) and Si_TT2 (n = 2).
Analysis of Differentially Expressed Genes (DEGs)
After removing the duplicated unigenes, approximately 59,846 genes expressed in the intestines of S. intermedius cultured at 20°C (Si_TT0) and 82,676 genes expressed in the intestines of S. intermedius cultured at 25°C (Si_TT2) were identified by Trinity with default parameters. DEGs between Si_TT0 (as control) and Si_TT2 were identified with a fold change ≥ 2.00 and probability ≥ 0.8. A total of 2,125 DEGs were identified, which included 1,015 upregulated and 1,110 downregulated unigenes as compared to those expressed in Si_TT0. Of the 2,125 DEGs, 732 genes were expressed specifically in Si_TT2, and 919 genes were expressed specifically in Si_TT0. qRT-PCR analysis indicated that the expression trends of 12 randomly selected DEGs (8 upregulated and 4 downregulated) were correlated well with those obtained in RNA-seq analysis, indicating the reliability and accuracy of the RNA-seq data obtained in this study (see Figure 4).
Figure 4. Verification of RNA-Seq by qRT-PCR. (A–L) The expression trends of 12 randomly selected DEGs. Values are expressed as the mean ± SD (n = 3). ∗Significant differences at p < 0.05 vs. control. ∗∗Extremely significant differences at p < 0.01 vs. control.
All of the identified DEGs were then annotated with GO terms. Approximately 79, 78, and 78 DEGs were clustered as categories related to metabolic process, cell, and cell parts, respectively (see Figure 5A). A total of 2,125 identified DEGs were enriched in 247 pathways by KEGG analysis. The ribosome pathways were the most enriched pathways, and followed by protein processing in endoplasmic reticulum pathways, prion disease pathways, and protein export pathways (see Figure 5B). Several DEGs related to growth (Supplementary Table S1), energy metabolism (Supplementary Table S2), heat shock responses (Supplementary Table S3), and immune response (Supplementary Table S4) were identified.
Figure 5. Functional annotation of DEGs of Si_TT2 vs. Si_TT0 in S. intermedius. (A) The most enriched GO terms of DEGs of S. intermedius Si_TT0and S. intermedius Si_TT2. The Y-axis represents the categories of annotated DEGs, and the X-axis represents the number of DEGs. (B) The top 20 enriched KEGG terms of DEGs of Si_TT2 vs. Si_TT0 in S. intermedius. The Y-axis represents the KEGG pathway, and the X-axis represents the enrichment factor. Dot size indicates the number of DEGs in the pathway. Dot colors corresponds to different Q-values.
For further elucidate interaction relationships, prediction analysis was conducted. After examining PPI networks, 173 proteins were found to be extremely well connected. In the 25°C treatments, 127 proteins were upregulated, 45 proteins were downregulated compared to the 20°C treatments (Figure 6A). These well-connected proteins included E3 ubiquitin-protein ligase MIB2 (GenBank Accession Number gi| 390364157), with 58 connections (Figure 6B); and heat shock protein 90 (GenBank Acc. No. gi| 390340697), with 29 connections (Figure 6C).
Figure 6. Differentially expressed genes interactive network prediction. (A) Interactive network prediction of all DEGs. (B) Interactive network prediction of MIB2. (C) Interactive network prediction of heat shock protein 90.
Transcription factors are also key regulators involved in heat-stress response in organisms. We identified 458 potential TFs from 25 TF families in S. intermedius cultured at 20 and 25°C. The identified TFs most commonly belonged to the following families: Cys2His2 protein (C2H2; 231 TFs; 50.44%), LIM domain-containing protein (LIM; 44 TFs; 9.61%), basic-helix-loop-helix protein (bHLH; 32 TFs; 6.99%), CCCH-type zinc finger protein (C3H; 22 TFs; 4.80%), and the basic leucine zipper protein (bZIP; 8 TFs; 1.75%) (Figure 7A). In the 25°C treatments, 40 TFs were upregulated and 11 TFs were downregulated compared to the 20°C treatments (Figure 7B and Supplementary Table S5).
Figure 7. Transcription factor families identified in S. intermedius Si_TT0 and S. intermedius Si_TT2 transcriptomes. (A) The distribution of all of the identified transcription factors. The number of each transcription factor is shown. (B) The distribution of differentially expressed transcription factors. The number and the expression of each differentially expressed transcription factor are shown.
Discussion
In the present study, we first investigated the impact of high water temperature on the growth, survival, and feeding behavior of S. intermedius. As we expected, high seawater temperature imparted significantly negative effects on the survival, growth, and feeding of S. intermedius. These results support the opinion that the growth, survival and feeding of sea urchins are sensitive to temperature (Siikavuopio et al., 2008; Pearce et al., 2010; Onitsuka et al., 2013). Moreover, we found that the lower SGRs of S. intermedius suffered from high seawater temperature stress were due to a decrease in feeding rates in this study, which agrees with the findings of Yu et al. (2016), who reported that high seawater temperature reduces the feeding and SGRs of sea cucumber A. japonicus. The observed decrease in food consumption and feeding rate not only explains the slow growth of S. intermedius, but also reflects to some extent the decrease in energy budget of S. intermedius under high temperature stress. This observation is also consistent with the results of Watts et al. (2011).
Organisms undergo metabolic adjustments in the presence of environmental stimuli. Moreover, cellar stress responses (CSRs) can also reflect the responses of an individual to environmental fluctuations (such as chronic high seawater temperature stress in this study) at the cellular and tissue levels. We therefore subsequently investigated the molecular mechanisms underlying high temperature-driven decrease in growth, survival, and feeding in S. intermedius by constructing four high-quality RNA-seq libraries and performing comparative transcriptome analysis to identify gene candidates associated with high-temperature responses.
In terms of the decreased energy budget, our data indicated that the expression of ATP5 [δ or the oligomycin sensitivity-conferring protein (OSCP)] mRNA significantly decreased in S. intermedius cultured at 25°C compared to those incubated at 20°C. ATP5 is one of subunits of ATP synthase, and silencing the expression of ATP5 can affect the activity of ATP synthase and further block the electron transfer chain in organisms. Our data suggest that high temperature stress can block ATP synthesis and affect the energy charge by decreasing atp5 expression in sea urchins. A decrease in ATP5 expression can alter total cellular ATP levels and impair growth in plants (Robison et al., 2009), which bears some resemblance to our observations. However, we could not find any studies involving echinoderm ATP5 to compare our results with. In addition to a reduction in energy, another other common strategy for organisms to adapt environmental stress (e.g., warming in this study) involves decreasing energetically expensive metabolic processes to extend their duration of tolerance (Pörtner and Farrell, 2008). Histones are dynamic proteins that can undergo multiple types of post-translational modifications and regulate gene expression depending on the metabolic state of the cell. Padilla-Gamino et al., 2013 reported that the downregulation histone-encoding genes is a principal transcriptional response accompanying metabolic depression in S. purpuratus larvae cultured at higher temperature (18°C), and that histones possibly act as metabolic sensors. We also observed a reduction in the expression of histone-encoding genes, such as Unigene17018_All (fold change: -2.42), Unigene9687_All (fold change: -7.70), and Unigene13915_All (fold change: -8.64). These findings support the hypothesis that sea urchins exposed to higher temperature consume energy supplies more rapidly than those cultured at optimal temperatures by reducing key metabolites for histone-modifying enzymes. In addition, another strategy for organisms to adapt environmental stress is minimizing body size or delaying growth. In this study, we found a significant reduction in expression of transforming growth factor beta (TGF-β) transcripts. TGF-β acts as a cytokine by imparting immunoregulatory effects, including lymphocyte proliferation, cytokine responsiveness, or cytokine expression (Ruscetti and Palladino, 1991). In echinoderms, TGF-β plays an important role in the symmetrical growth of sea urchin embryos and the biomineralization of larval skeletogenesis (Zito et al., 2003). In this study, we postulate that the reduced growth under high temperature stress in S. intermedius might be due to alterations in TGF signals. Further studies should be conducted to clarify the mechanisms underlying how high-temperature influences the TGF signal pathway.
Heat-shock proteins (HSPs) are molecular chaperones with multiple functions, including stress resistance and adaption to environmental changes in various species (Currie, 2011). The upregulation of Hsps is one of ubiquitous mechanisms of marine organism in coping with thermal stress (Osovitz and Hofmann, 2005). Our data also identified altered mRNA expression levels of some Hsps. Runcie et al. (2012) demonstrated that higher temperatures (18°C vs. 12°C) can cause mild embryonic developmental stress and increase both hsp70 and hsp90 expression in embryos of the sea urchin S. purpuratus. Vergara-Amado et al. (2017) reported that transient warmer temperature treatments (18°C vs. 14°C) induces the up-regulation of hsp70 and hsp90 in juveniles of the sea urchin L. albus. As expected, we found the upregulation of hsp70 and hsp90 in adult S. intermedius after application of chronic high-temperature stress. Hsp70 and Hsp90 have both been implicated in the proteasomal degradation of chaperoned client proteins (Kiang and Tsokos, 1998; Pratt, 1998). Hsp90, in particular, has been characterized as the driver for Hsp-mediated proteasomal degradation. Hsp90 is a protective chaperone when in complex with p50/immunophilin, p23, and ATP, but drives client proteins to poly-ubiquitination and proteasomal degradation when ATP-depleted and bound in complex with Hsp70 (Doong et al., 2003). Interestingly, the present study observed that mind bomb-2 (mib2) was downregulated under high temperature (25°C) conditions as compared to that observed using the optimal temperature (20°C). MIB2 has E3 ubiquitin-protein ligase activity and can promote the ubiquitination and endocytosis of its client protein (Koo et al., 2005). PPI prediction indicated the upregulation of Hsp90 and the downregulation of E3 ubiquitin-protein ligase MIB2 under high temperature (25°C) conditions as compared to that observed using the optimal temperature (20°C) (Figure 6B). This observation suggests that high temperature-induced proteasomal degradation or apoptosis might mainly rely on the Hsp-mediated proteasomal degradation pathway rather than that of being MIB-2-mediated. In addition, the cytokine-like function of Hsp70 what been well documented in several studies (Zhang et al., 2011; Yang et al., 2016; Ying et al., 2016). Combined with the findings of the present study, these observations support the hypothesis that Hsp70 and Hsp90 regulate not only proteasomal degradation but also immune responses when sea urchins are subjected to high-temperature stress (Ying et al., 2016). Since enhanced hsp70 and hsp90 transcripts can be detected from embryos to adults in sea urchins during high temperature stress regardless of whether such exposure is chronic or transient, we therefore hypothesize that the hsp70 and hsp90 genes be included in selective breeding or assisted breeding of high temperature-resistant sea urchins.
The significant upregulation of the glutathione S-transferases (GST) gene was also observed in the present study. GST is a phase II detoxification isozyme that catalyzes the conjugation of glutathione with both xenobiotics and endogenous substrates. GST activity has long been utilized as a bioindicator of environmental contamination in coastal regions (Cunha et al., 2005). Field studies have shown that the maximal GST activity in sea urchins and mussels can be observed in the summer (Moreira and Guilhermino, 2005). However, no study has measured echinoderm GST levels under thermal-stress, and thus we were unable to clarify the relationship between GST activities and high temperature stress in sea urchins.
Transcription factors and cis-acting elements are conserved mechanisms that regulate gene transcription (Murray et al., 1988). Among differentially expressed TFs, we observed the significantly upregulation of mRNA expression of some multiple function TFs in sea urchins in response to high temperature stress, which include activator protein-1 (AP-1), FOS, cAMP-response element-binding protein (CREB), and the zinc finger (ZNF) proteins. AP-1 and FOS are members of the basic leucine zipper protein (bZIP) family. Activated AP-1 has been demonstrated to be a stress-responsive TF and plays a key role in responding to environmental stimuli by regulating various immune signal transduction pathways, such as the Toll-like receptor (TLR), tumor necrosis factor alpha (TNF-α), and mitogen-activated protein kinase (MAPK) pathways, in marine organisms (Qu et al., 2015; Zhan et al., 2018). Fos proteins are a key part of the AP-1 complex and can regulate a wide range of biological process (Hirayama et al., 2005). Production of many immune-related molecules (antioxidant enzymes, chemokines, and interleukin) require Fos expression (Rimoldi et al., 2009). Additionally, it has been shown that Fos cooperates with Notch to regulate cell fate specification of intermediate precursors during Caenorhabditis elegans development (Oommen and Newman, 2007). CREB is one of multi-function TF regulating various signal transduction pathways. In Drosophila, CREB and Hsp70 can additively suppress polyglutamine-mediated toxicity (Iijimaando et al., 2005). In cultured rat primary hippocampal neuron cells, CREB has been demonstrated to be activated by the Hsp90/Akt signal pathway (Cen et al., 2006). Combined with our comparative expression data on hsp70 and hsp90, we postulate that the CREB-mediated gene regulation network of higher thermal tolerance in S. intermedius might be closely correlated with HSPs. Further studies should be conducted to confirm this hypothesis. ZNFs can be found in all eukaryotes and act as TFs that plays critical role in responding to environmental stimuli such as biometals (Villalpando et al., 2017). In plants, ZNFs have been identified as heat response-related gene candidates (Wang et al., 2010; Yan et al., 2016; Fang et al., 2017). The present study has shown that ZNFs are involved in echinoderm heat responses and more studies should be conducted to elucidate the ZNF regulated genes and signal transduction pathways during chronic heat stress in echinoderms.
Moreover, the adjustment of gene structure is the ultimate mechanism for organisms to adapt to long-term stress and maintain population size. Our comparative transcriptomic data also indicated altered SNPs and SSRs. These results will facilitate studies on the genetic structure, population geography, and ecology of sea urchins. Further screening of these SNPs and SSRs may assist in the identification of more valuable heat-resistant genetic markers that may be utilized in the selective breeding of heat-resistant sea urchins.
However, significant expression alterations of some transient thermal-tolerance candidate genes such as cytochrome P450 (CYP450), Na+/K+ ATPase (Vergara-Amado et al., 2017) were not observed in the present study. This suggests that high temperature-induced molecular responses in the same species depend on individual developmental stages and the duration of stress.
Conclusion
The present study has demonstrated that chronically high seawater temperatures negatively influence the growth, survival, and feeding of the sea urchin S. intermedius. Gene candidates (e.g., HSPs, cytokines, and TFs) that were closely correlated with thermal resistance and adaptation were identified by comparative transcriptomic analysis. In summary, our results provide insight into the genes and regulatory networks involved in chronic thermal stress in S. intermedius, as well as enriching transcriptomic and genetic resources for sea urchins and other invertebrates. More studies on the molecular events involved in the thermal resistance and adaption mechanisms for S. intermedius should be conducted to better understand the impact of chronically high temperature stress on sea urchin physiology and ecology.
Data Availability
The datasets generated for this study can be found in NCBI, PRJNA508827.
Author Contributions
YC and YZ conceived and designed the experiments. JL, JS, WZ, YL, DC, and WH performed the experiments. JL, YZ, and WZ analyzed the data. YZ and JL wrote the manuscript. All authors read and approved the manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (31672652) and the grant for Chinese Outstanding Talents in Agricultural Scientific Research (for YC) supported this study.
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00301/full#supplementary-material
Footnotes
- ^ http://pgrc.ipk-gatersleben.de/misa/misa.html
- ^ http://string-db.org/
- ^ http://busco.ezlab.org/datasets/eukaryota_odb9.tar.gz
References
Branco, P. C., Borges, J. C., Santos, M. F., and Jensch Junior, B. E. (2013). The impact of rising sea temperature on innate immune parameters in the tropical subtidal sea urchin Lytechinus variegatus, and the intertidal sea urchin Echinometra lucunter. Mar. Environ. Res. 92, 95–101. doi: 10.1016/j.marenvres.2013.09.005
Brothers, C. J., and Mcclintock, J. B. (2015). The effects of climate-induced elevated seawater temperature on the covering behavior, righting response, and Aristotle’s lantern reflex of the sea urchin Lytechinus variegatus. J. Exp. Mar. Biol. Ecol. 467, 33–38. doi: 10.1016/j.jembe.2015.02.019
Byrne, M., Ho, M., Selvakumaraswamy, P., Nguyen, H. D., Dworjanyn, S. A., and Davis, A. R. (2009). Temperature, but not pH, compromises sea urchin fertilization and early development under near-future climate change scenarios. Proc. R. Soc. B Biol. Sci. 276, 1883–1888. doi: 10.1098/rspb.2008.1935
Cen, X., Nitta, A., Ohya, S., Zhao, Y., Ozawa, N., Mouri, A., et al. (2006). An analog of a dipeptide-like structure of FK506 increases glial cell line-derived neurotrophic factor expression through camp response element-binding protein activated by heat shock protein 90/AKT signaling pathway. J. Neurosci. Off. J. Soc. Neurosci. 26:3335. doi: 10.1523/jneurosci.5010-05.2006
Chang, Y. Q., Ding, J., Song, J., and Yang, W. (2004). Biological Research and Aquaculture of Sea Cucumber and Sea Urchin. Beijing: Ocean Press. 217–218.
Chang, Y. Q., Tian, X. F., Han, F. J., Chen, S., Zhou, M., et al. (2016). Family growth and survival response to two simulated water temperature environments in the sea urchin Strongylocentrotus intermedius. Int. J. Mol. Sci. 17, 1356. doi: 10.3390/ijms17091356
Cunha, I., García, L. M., and Guilhermino, L. (2005). Sea-urchin (Paracentrotus lividus) glutathione S-transferases and cholinesterase activities as biomarkers of environmental contamination. J. Environ. Monitor. 7, 288–294. doi: 10.1039/b414773a
Currie, S. (2011). Temperature heat shock proteins and temperature. Encycl. Fish Physiol. 3, 1732–1737. doi: 10.1016/B978-0-12-374553-8.00196-9
Doong, H., Rizzo, K., Fang, S., Kulpa, V., Weissman, A. M., and Kohn, E. C. (2003). CAIR-1/BAG-3 abrogates heat shock protein-70 chaperone complex-mediated protein degradation accumulation of poly-ubiquitinated hsp90 client proteins. J. Biol. Chem. 278, 28490–28500. doi: 10.1074/jbc.m209682200
Fang, L., Tong, J., Dong, Y., Xu, D., Mao, J., and Zhou, Y. (2017). De novo RNA sequencing transcriptome of Rhododendron obtusum identified the early heat response genes involved in the transcriptional regulation of photosynthesis. PLoS One 12:e0186376. doi: 10.1371/journal.pone.0186376
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat. Biotechnol. 29:644. doi: 10.1038/nbt.1883
Hirayama, J., Cardone, L., Doi, M., and Sassone-Corsi, P. (2005). Common pathways in circadian and cell cycle clocks: light-dependent activation of FOS/AP-1 in zebrafish controls CRY-1a and WEE-1. Proc. Natl. Acad. Sci. U.S.A. 102, 10194–10199. doi: 10.1073/pnas.0502610102
Iijimaando, K., Wu, P., Drier, E. A., Iijima, K., and Yin, J. C. P. (2005). cAMP-response element-binding protein and heat-shock protein 70 additively suppress polyglutamine-mediated toxicity in Drosophila. Proc. Natl. Acad. Sci. U.S.A. 102, 10261–10266. doi: 10.1073/pnas.0503937102
Kiang, J. G., and Tsokos, G. C. (1998). Heat shock protein 70 kDa: molecular biology, biochemistry, and physiology. Pharmacol. Therap. 80, 183–201. doi: 10.1016/s0163-7258(98)00028-x
Koo, B. K., Yoon, K. J., Yoo, K. W., Lim, H. S., Song, R., So, J. H., et al. (2005). Mind bomb-2 IS an e3 ligase for notch ligand. J. Biol. Chem. 280, 22335–22342. doi: 10.1074/jbc.m501631200
Lawrence, J. M., Cao, X. B., Chang, Y. Q., Wang, P., Yu, Y., Lawrence, A. L., et al. (2009). Temperature effect on feed consumption, absorption, and assimilation efficiencies and production of the sea urchin Strongylocentrotus intermedius. J. Shellf. Res. 28, 389–395. doi: 10.2983/035.028.0223
Li, B., and Dewey, C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12:323. doi: 10.1186/1471-2105-12-323
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Mckenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The genome analysis toolkit: a mapreduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. doi: 10.1101/gr.107524.110.20
Moreira, S. M., and Guilhermino, L. (2005). The use of Mytilus galloprovincialis acetylcholinesterase and glutathione S-transferases activities as biomarkers of environmental contamination along the northwest portuguese coast. Environ. Monitor. Assess. 105, 309–325. doi: 10.1007/s10661-005-3854-z
Murray, R. K., Granner, D. K., Mayes, P. A., and Rodwell, V. W. (1988). Harper’s Biochemistry. New York, NY: Prentice Hall.
Onitsuka, T., Unuma, T., and Kawamura, T. (2013). Effects of diet and temperature on post-settlement growth and survival of the short-spined sea urchin Strongylocentrotus intermedius. Aquacult. Res. 45, 1–15. doi: 10.1111/j.1365-2109.2012.03198.x
Oommen, K. S., and Newman, A. P. (2007). Co-regulation by Notch and Fos is required for cell fate specification of intermediate precursors during C. elegans uterine development. Development 134:3999. doi: 10.1242/dev.002741
Osovitz, C. J., and Hofmann, G. E. (2005). Thermal history-dependent expression of the hsp70 gene in purple sea urchins: biogeographic patterns and the effect of temperature acclimation. J. Exp. Mar. Biol. Ecol. 327, 134–143. doi: 10.1016/j.jembe.2005.06.011
Padilla-Gamino, J. L., Kelly, M. W., Evans, T. G., and Hofmann, G. E. (2013). Temperature and CO2 additively regulate physiology, morphology and genomic responses of larval sea urchins, Strongylocentrotus purpuratus. Proc. Biol. Sci. 280:20130155. doi: 10.1098/rspb.2013.0155
Pearce, C. M., Williams, S. W., Yuan, F., Castell, J. D., and Robinson, S. M. C. (2010). Effect of temperature on somatic growth and survivorship of early post-settled green sea urchins, Strongylocentrotus droebachiensis (müller). Aquacult. Res. 36, 600–609. doi: 10.1111/j.1365-2109.2005.01264.x
Pörtner, H. O., and Farrell, A. P. (2008). Physiology and climate change. Science 322, 690–692. doi: 10.1126/science.1163156
Pratt, W. B. (1998). The hsp90-based chaperone system: involvement in signal transduction from a variety of hormone and growth factor receptors. Proc. Soc. Exp. Biol. Med. Soc. Exp. Biol. Med. 217, 420. doi: 10.3181/00379727-217-44252
Qin, Y. J., Li, X., Wu, L. X., and Zhou, Y. B. (2011). Effects of starvation and refeeding on metabolism and growth in sea urchin, Strongylocentrotus intermedius. J. Dalian Ocean Univers. 26, 521–525.
Qu, F., Xiang, Z., Wang, F., Zhang, Y., Tong, Y., Li, J., et al. (2015). A novel molluscan fos gene with immune defense function identified in the Hong Kong oyster, Crassostrea hongkongensis. Dev. Compar. Immunol. 51, 194–201. doi: 10.1016/j.dci.2015.03.012
Rimoldi, S., Terova, G., Brambilla, F., Bernardini, G., Gornati, R., and Saroglia, M. (2009). Molecular characterization and expression analysis of Na+/H+ exchanger (NHE)-1 and c-Fos genes in sea bass (Dicentrarchus labrax, L) exposed to acute and chronic hypercapnia. J. Exp. Mar. Biol. Ecol. 375, 32–40. doi: 10.1016/j.jembe.2009.05.002
Robison, M. M., Ling, X., Smid, M. P. L., Zarei, A., and Wolyn, D. J. (2009). Antisense expression of mitochondrial ATP synthase subunits oscp (ATP5) and (ATP3) alters leaf morphology, metabolism and gene expression in Arabidopsis. Plant Cell Physiol. 50, 1840–1850. doi: 10.1093/pcp/pcp125
Runcie, D. E., Garfield, D. A., Babbitt, C. C., Wygoda, J. A., Mukherjee, S., and Wray, G. A. (2012). Genetics of gene expression responses to temperature stress in a sea urchin gene network. Mol. Ecol. 21, 4547–4562. doi: 10.1111/j.1365-294X.2012.05717.x
Ruscetti, F. W., and Palladino, M. A. (1991). Transforming growth factor-beta and the immune system. Prog. Growth Factor Res. 3, 159–175. doi: 10.1016/S0955-2235(05)80006-7
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.3372/wi.43.43119
Siikavuopio, S. I., Mortensen, A., and Christiansen, J. S. (2008). Effects of body weight and temperature on feed intake, gonad growth and oxygen consumption in green sea urchin, Strongylocentrotus droebachiensis. Aquaculture 281, 77–82. doi: 10.1016/j.aquaculture.2008.05.033
Simão, F. A., Waterhouse, R. M., Panagiotis, I., Kriventseva, E. V., and Zdobnov, E. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212. doi: 10.1093/bioinformatics/btv351
Szklarczyk, D., Franceschini, A., Kuhn, M., Simonovic, M., Roth, A., Minguez, P., et al. (2011). The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 39, 561–568. doi: 10.1093/nar/gkq973
Tarazona, S., Garcíaalcalde, F., Dopazo, J., Ferrer, A., and Conesa, A. (2011). Differential expression in RNA-seq: a matter of depth. Genome Res. 21, 2213–2223. doi: 10.1101/gr.124321.111
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. Ther. Biol. 68(Pt A). 104–109. doi: 10.1016/j.jtherbio.2017.03.009
Villalpando, J. L., Arreola, R., Puente-Rivera, J., Azuara-Liceaga, E., Valdés, J., López-Canovas, L., et al. (2017). TvZNF1 is a C2H2 zinc finger protein of Trichomonas vaginalis. Biometals 30, 1–12. doi: 10.1007/s10534-017-0053-7
Wang, L. J., Fan, L., Loescher, W., Duan, W., Liu, G. J., Cheng, J. S., et al. (2010). Salicylic acid alleviates decreases in photosynthesis under heat stress and accelerates recovery in grapevine leaves. BMC Plant Biol. 10:34. doi: 10.1186/1471-2229-10-34
Wang, Q. L., Yu, S., and Dong, Y. (2015). Parental effect of long acclimatization on thermal tolerance of juvenile sea cucumber, Apostichopus japonicus. PLoS One 10:e0143372. doi: 10.1371/journal.pone.0143372
Watts, S. A., Hofer, S. C., Desmond, R. A., Lawrence, A. L., and Lawrence, J. M. (2011). The effect of temperature on feeding and growth characteristics of the sea urchin Lytechinus variegatus fed a formulated feed. J. Exp. Mar. Biol. Ecol. 397, 188–195. doi: 10.1016/j.jembe.2010.10.007
Yan, J., Yu, L., Xuan, J., Lu, Y., Lu, S., and Zhu, W. (2016). De novo transcriptome sequencing and gene expression profiling of spinach (Spinacia oleracea L.) leaves under heat stress. Sci. Rep. 6:19473. doi: 10.1038/srep19473
Yang, A. F., Zhou, Z. C., Ying, D., Bei, J., Wang, X. Y., Zhong, C., et al. (2010). Stability comparison of cytb and β-actin genes expression in sea cucumber (Apostichopus japonicus). J. Agric. Sci. Technol. 12, 79–84.
Yang, Y., Yu, H., Li, H., Wang, A., and Yu, H. Y. (2016). Effect of high temperature on immune response of grass carp (Ctenopharyngodon idellus) by transcriptome analysis. Fish Shellf. Immunol. 58, 89–95. doi: 10.1016/j.fsi.2016.09.014
Ying, Y., Hui, Y., Hua, L., and Wang, A. (2016). Transcriptome profiling of grass carp (Ctenopharyngodon idellus) infected with aeromonas hydrophila. Fish Shellf. Immunol. 51, 329–336. doi: 10.1016/j.fsi.2016.02.035
Yu, H., Zhang, C., Gao, Q., Dong, S., Ye, Z., and Tian, X. (2016). Impact of water temperature on the growth and fatty acid profiles of juvenile sea cucumber Apostichopus japonicus (Selenka). J. Ther. Biol. 60, 155–161. doi: 10.1016/j.jtherbio.2016.07.011
Zeng, G. E., Lian, S. M., Cheng, X. H., Hua, Z. L., and Qi, Y. Q. (2006). EOF analysis of intra-seasonal variabilities of SST in the east china sea and yellow sea. Adv. Mar. Sci. 24, 146–155.
Zhan, Y. Y., Wang, Y., Li, K. Q., Song, J., and Chang, Y. Q. (2018). A novel p38 MAPK gene in the sea cucumber Apostichopus japonicus (Ajp38) is associated with the immune response to pathogenic challenge. Fish Shellf. Immunol. 81, 250–259. doi: 10.1016/j.fsi.2018.07.001
Zhang, A., Zhou, X., Wang, X., and Zhou, H. (2011). Characterization of two heat shock proteins (Hsp70/Hsc70) from grass carp (Ctenopharyngodon idella): evidence for their differential gene expression, protein synthesis and secretion in LPS-challenged peripheral blood lymphocytes. Compar. Biochem. Physiol. Part B Biochem. Mol. Biol. 159, 109–114. doi: 10.1016/j.cbpb.2011.02.009
Zhang, H., Jiang, Z. T., and Liu, E. X. (2016). Seasonal variation characteristics of temperature in the Yellow Sea and the East China Sea. China Water Transport. 16, 167–170.
Zhang, L. S., Zhang, L. L., Shi, D. T., Wei, J., Chang, Y. Q., and Zhao, C. (2017). Effects of long-term elevated temperature on covering, sheltering and righting behaviors of the sea urchin Strongylocentrotus intermedius. PeerJ 5:e3122. doi: 10.7717/peerj.3122
Zhang, Z., Schwartz, S., Wagner, L., and Miller, W. (2000). A greedy algorithm for aligning DNA sequences. J. Comput. Biol. 7, 203–214. doi: 10.1089/10665270050081478
Keywords: heat stress, Strongylocentrotus intermedius, growth, survival, feeding, transcriptome
Citation: Zhan Y, Li J, Sun J, Zhang W, Li Y, Cui D, Hu W and Chang Y (2019) The Impact of Chronic Heat Stress on the Growth, Survival, Feeding, and Differential Gene Expression in the Sea Urchin Strongylocentrotus intermedius. Front. Genet. 10:301. doi: 10.3389/fgene.2019.00301
Received: 12 December 2018; Accepted: 19 March 2019;
Published: 04 April 2019.
Edited by:
Peng Xu, Xiamen University, ChinaReviewed by:
Yang Yu, Institute of Oceanology (CAS), ChinaJie Mei, Huazhong Agricultural University, China
Copyright © 2019 Zhan, Li, Sun, Zhang, Li, Cui, Hu and Chang. 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: Yaqing Chang, yqkeylab@hotmail.com