Hanxi Li
Xuena Huang
Aibin Zhan- 1Key Laboratory of Environmental Biotechnology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, China
- 2College of Resources and Environment, University of the Chinese Academy of Sciences, Chinese Academy of Sciences, Beijing, China
Fluctuating environmental changes impose tremendous stresses on sessile organisms in marine ecosystems, in turn, organisms develop complex response mechanisms to keep adaptive homeostasis for survival. Physiological plasticity is one of the primary lines of defense against environmental challenges, and such defense often relies on the antioxidant defense system (ADS). Hence, it is imperative to understand response mechanisms of ADS to fluctuating environments. Invasive species provide excellent models to study how species cope with environmental stresses, as invasive species encounter sudden, and often recurrent, extensive environmental challenges during the whole invasion process. Here, we studied the roles of ADS on rapid response to recurrent cold challenges in a highly invasive tunicate (Ciona robusta) by simulating cold stresses during its invasion process. We assessed antioxidative indicators, including malondialdehyde (MDA), total antioxidant capacity (T-AOC), superoxide dismutase (SOD), catalase (CAT), and glutathione (GSH), as well as transcriptional changes of ADS-related genes to reveal the physiological plasticity under recurring cold stresses. Our results demonstrated that physiological homeostasis relied on the resilience of ADS, which further accordingly tuned antioxidant activity and gene expression to changing environments. The initial cold stress remodeled baselines of ADS to promote the development of stress memory, and subsequent stress memory largely decreased the physiological response to recurrent environmental challenges. All results here suggest that C. robusta could develop stress memory to maintain physiological homeostasis in changing or harsh environments. The results obtained in this study provide new insights into the mechanism of rapid physiological adaption during biological invasions.
Introduction
Nowadays, organisms often experience frequent and intensive environmental challenges derived from multiple interacting stressors, such as global climate change and habitat transitions influenced by increasing human activities (e.g., biological invasions; Pickering et al., 2013; Calosi et al., 2017). Amongst several strategies for mitigating those stresses, physiological plasticity is a crucial mechanism to tolerate fluctuating environmental changes (Somero, 2012; Jimenez et al., 2015; Spicer et al., 2019). It is incontrovertible that the adaptation within one generation is critically important to buffer environmental stresses to ensure survival, which is the prerequisite for subsequent adaptation and evolution of a population (Carlson et al., 2014). Rapid physiological regulation serves as one of the most direct and effective ways to mitigate threats to the survival of a species within one generation (Hofmann and Todgham, 2010; Seebacher et al., 2014). Thus, understanding physiological response dynamics is crucial for deeply investigating the homeostasis maintenance during environmental challenges and further predicting the influence of environmental changes on the survival and geographical distribution shifts of a species (Hoffmann and Sgro, 2011; Nandamuri et al., 2017).
Rapid environmental changes often trigger excessive loads of reactive oxygen species (ROS) by boosting the mitochondrial respiratory activity, reduce the fluidity of the mitochondrial membrane, and ruin the electron transport chain (i.e., causing oxidative stress; Murphy, 2009). Antioxidant defense system (ADS), as the first defense line and one of the most vital physiological mechanisms, can effectively eliminate the imbalance state of redox and protect organismal homeostasis against environmental stresses (Storey, 1996; Jimenez et al., 2015; Calabrese and Mattson, 2017). By activating antioxidants which can catalyze the decomposition of ROS into less toxic molecules within seconds to minutes, organisms keep the content of free-radicals and reactive oxidants at appropriate levels during acute environmental challenges (Sthijns et al., 2016; Oliveira et al., 2018). On the other hand, ADS can change the synthesis rate and content of downstream genes, including most of antioxidants, by altering the expression of transcriptional factors (TF) in ADS-related signaling pathways (Zhang, 2006; Kensler et al., 2007; Shadel and Horvath, 2015). The nuclear factor erythroid-2-related factor 2 (Nrf2)-keap1 signaling pathway is one of the most relevant and extensively studied pathways in resistance to oxidative stresses (Nguyen et al., 2003; Ristow, 2014; Sies et al., 2017). Nrf2 is suppressed by bounding to Kelch-like epichlorohydrin associating protein1 (Keap1) in a normal condition and lose the ability of translocation to the nucleus. While under stresses, Keap1 is repressed and the released Nrf2 is subsequently translocated forward to the nucleus where it accumulates substantially to activate the expression of its target genes (Ma, 2013). By modulating the expression of its target genes involved in antioxidants, such as superoxide dismutase (SOD), glutathione peroxidase (GPx), catalase (CAT), phase II detoxifying enzyme (e.g., glutathione S-transferase, GST; heme oxygenase-1, HO-1) and other antioxidative compounds, the Nrf2-keap1 signaling pathway can decline the content of free-radicals and/or reactive oxidants to a harmless level, thus protecting organisms against various environmental stresses for maintaining redox homeostasis ultimately (Schieber and Chandel, 2014; Davies, 2016; Sparks et al., 2019).
Recent studies have demonstrated that previously encountered abiotic stresses could affect organismal response against later stresses (i.e., stress memory; Crisp et al., 2016). Available evidence showed that organisms could respond more effectively to recurring stimulus through two distinct types of stress memory: type I shows the enhanced response with a faster and/or stronger pattern, and such a type could provide the benefits of potentiation to subsequent challenges. Type II has a feature of increased tolerance with steadier and less intense response, and such a type can provide the benefits of tolerance during recurrent stresses (Avramova, 2015). For example, after being subjected to several dehydration/rehydration cycles, Arabidopsis thaliana and Zea mays (maize) showed the transcriptional stress memory mainly by producing stronger reactions to improve the ability of moisture retention when compared to unstressed plants (Ding et al., 2013; Virlouvet and Fromm, 2015). The red seaweed Bangia fuscopurpurea could alter specific membrane fatty acid composition to own heat memory and increase heat tolerance to protect against future heat stresses (Kishimoto et al., 2019). However, in marine ecosystems, we have very limited knowledge on the mechanism of physiological response based on ADS under fluctuating environments (except Mytilus californianus, Jimenez et al., 2015; Crassostrea gigas, Moreira et al., 2016).
Biological invasions provide an excellent “natural experiment” to study how species cope with environmental stresses, as invasive species encounter sudden, and often repeated, extensive environmental challenges during the invasion process (Zhan et al., 2015; Ni et al., 2018). Particularly for human activity-mediated introductions, invasive species often encounter environmental changes much greater and/or faster than what species can experience under natural conditions such as seasonal fluctuations (Zhan et al., 2015; Chen et al., 2018a). Taking shipping-mediated invasions as an example, environmental conditions can change dramatically in ballast tanks, with temperature suddenly dropped from 20°C to the freezing point, and this type of great changes can be recurrent during the whole shipping voyage (Briski et al., 2013; Chen et al., 2018b). The response to the former stress may cause stress memory via achieving a renewed balance of homeostasis and tolerance to influence subsequent stress response. Thus, invasive species represent good candidates to study stress memory, in turn, the obtained results on stress memory can largely promote our understanding of invasion success, particularly on mechanisms of physiologic plasticity under environmental challenges.
Ascidians are excellent models to investigate how species respond to repeated harsh environments in coastal ecosystems (Zhan et al., 2015; Huang et al., 2017). Survival of repeated acute stresses during the transport and establishment stages is an indispensable ability and prerequisite for successful invasions to new habitats in one generation (Lenz et al., 2011). Many invasive species, especially for sessile organisms such as ascidians, have no chance to moderate stresses via behaviors such as escape from harsh environments, but largely rely on strong adaptive homeostasis through pre-adaptation or physiological plasticity to maintain homeostasis (Richards et al., 2006; Bruce et al., 2007). The highly invasive ascidian Ciona robusta, which is presumably a native to Northern Atlantic, has successfully invaded many coasts globally (Carver et al., 2006; Rius et al., 2014). It has become a popular model animal for evolutionary biology and invasion science (Delsuc et al., 2006; Zhan et al., 2015). C. robusta can tolerate a wide temperature range from −1°C to 35°C (Dybern, 1965). Amazingly, as a simultaneous hermaphrodite, C. robusta can produce eggs under low temperature as low as −1°C (Carver et al., 2003). Hence, the capacity of cold tolerance makes C. robusta a good model to explore mechanisms of physiological response to recurrent cold stresses.
Here we use this model species to study physiological plasticity represented by the regulation of Nrf2-depended ADS at the levels of both enzyme activity and transcriptional synthesis in response to repeated cold stresses. Our study used two rounds of cold stress-recovery to simulate repeated environmental changes during the invasion process. We aim to (1) explore the response profile of ADS of C. robusta under two rounds of cold stresses, and (2) verify the existence of stress memory based on ADS by comparing different responses on the variation of antioxidant indicators and ADS-related genes expression. The results obtained here are expected to promote our understanding of the maintenance of physiological homeostasis in fluctuating environments and provide new insights into the rapid physiological adaption process during biological invasions.
Materials and Methods
Sample Collection and Acclimation
Adult ascidians were collected from scallop cages in the Longwangtang maricultural area of Dalian, Liaoning Province, China (38°48′53″ N, 121°24′06″ E) in September 2018. Ascidians were transported into tanks in the laboratory and acclimated for one week at the temperature of 20 ± 1°C according to the field condition. In the laboratory, all animals were fed twice daily with algae powder mixture of wall-broken Chlorella (Chlorella spp.) and spirulina (Arthrospira platensis). At the end of the acclimation, 12 adults were randomly chosen as the unstressed group.
Cold Stress Treatment and Recovery
We used two cycles of “cold stress-recovery” experiment to simulate the recurring environmental stresses during the invasion process. The whole experiment consisted of two rounds (Round I and II) of challenges, each round included two phases of cold stress (1S and 2S) and recovery to ambient condition (1R, 2R), and every single phase lasted 24 h based on surveys from our former study (Huang et al., 2018). In order to clearly describe the comparison between response within two rounds, we defined samples from the ambient condition as “unstressed” and the final status of the first round (1R-24) as “prestressed.” The temperature of cold stress was set at 5 ± 1°C according to the minimum temperature within the normal tolerance range and the recovery temperature was 20 ± 1°C (i.e., field condition, see Figure 1 for the experimental design).
Figure 1. Experimental design scheme. Ascidians were acclimated for one week and then transferred to the temperature-controlled aquarium at 5 ± 1°C for cold stress treatment (1S). After 24-hour treatment, the temperature was recovered to 20 ± 1°C and lasted for another 24-hour as the first recovery stage (1R). The ascidians that experienced the first round of stress-recovery treatment were then subjected to the new round of stress-recovery experiment, which was exactly the same as the Round I. Red points indicate the time of sampling.
During the first stress (1S), 200 adult ascidians with similar body size were transferred from the acclimation tank at 20 ± 1°C to the temperature-controlled aquarium (80 cm × 45 cm × 40 cm) at 5 ± 1°C for the cold challenge. The density of ascidians was approximately 1.4 ind/L. After 24 h of cold exposure, the temperature was recovered to 20 ± 1°C and lasted for 24 h as the first recovery phrase (1R). The treated ascidians by the first round of stress-recovery were then subjected to a new round of stress-recovery experiment. Six individuals were randomly sampled at different sampling time points including 1, 6, and 24 h at each of the four phases (Figure 1). To obtain an accurate profile of ADS changes during acute cold stresses, we also sampled two additional points at 3-hour (1S-3) and 12-hour (1S-12) of the first phase (Figure 1). Ascidians were dissected on ice-chilled culture dishes and separated into two parts, including 100–150 mg pharynx muscle mass for RNA isolation and the rest of the pharynx muscle for physiological analyses for each individual. After dissection, all samples were flash-frozen in liquid nitrogen and stored at −80°C. Species identification of all samples was conducted according to Zhan et al. (2010) and Smith et al. (2012) to avoid species confusion due to the morphological similarity among Ciona ascidians.
Antioxidant Status Assay
Malondialdehyde (MDA) Content and Total Antioxidant Capacity (T-AOC) Assay
Tissue supernatant was prepared for each bioindicator analysis. Each sample (80 mg pharynx muscle) was individually prepared in 0.9% sodium chloride solution with the proportion ratio of 1: 9 (w/v) and homogenized by tissue homogenizer in ice bath. Smear test was taken after homogenization to confirm the accomplishment of cell disruption. Subsequently, the homogenate was centrifuged at 10,000 g for 20 min under 4°C. Supernatants were immediately used for the determination of the total protein (TP), MDA concentration, T-AOC, special total SOD activity, special CAT activity and GSH concentration. All the following indices for specific enzymatic activity were defined as the concentration/activity unit of 1 mg tissue protein. TP concentration was determined using the Kaumas blue method-based kit (Nanjing Jiancheng, China).
The comprehensive and accurate assessment of antioxidant degrees was measured by both cellular damage products and antioxidant capacity. The cellular damage level was represented by the extent of lipid peroxidation, which was further indicated by the concentration of MDA using the thiobarbituric acid (TBA) method-based kit (Nanjing Jiancheng, China). T-AOC was measured by the assay kit (Nanjing Jiancheng, China) according to the instructions.
Antioxidant Enzyme Activity Assay
Animals rely on robust ADS to maintain the redox balance, including antioxidase system and other non-enzymatic antioxidant molecules. In this study, we selected three types of antioxidants, including total superoxide dismutase (T-SOD), catalase (CAT), and reduced glutathione (GSH). These indicators were tested by relevant assay kits (Nanjing Jiancheng, China). T-SOD was detected by the hydroxylamine method, whereas CAT and GSH were measured by the ultraviolet method and dithiol dinitrobenzoic acid method, respectively. All indicators mentioned above were measured by an ultraviolet spectrophotometer (DR6000, HACH Inc., United States).
Total RNA Isolation and RT-PCR
Total RNA was extracted from the muscle tissue (50–100 mg) by TRIzol reagent (Ambion, United States) with the standard protocol. RNA integrity was tested by agarose gel electrophoresis and RNA concentration was determined by Nano 2000 spectrophotometer (Nanodrop Technologies, United States). A total of 3 μg high-quality RNA was reverse-transcribed to cDNA using the PrimerScriptII cDNA Synthesis Kit (Takara, Japan). Then cDNA was treated with RNase-free DNase I (Promega, United States) to avoid genomic DNA contamination, and 10-fold cDNA dilution was applied in downstream assays.
We assessed the gene expression of two transcription factors (Nrf2 and Keap1) and their target genes including Cu/ZnSOD, MnSOD, CAT, GPX, and GST to monitor the transcriptional change under repeated cold stresses. The β-actin gene was chosen as the reference gene for gene expression normalization (Fujikawa et al., 2010; Huang et al., 2016; Huang et al., 2019). The sequences of all genes were downloaded from Ensembl1 and all primers were designed by Primer Premier 5.02. Primers were tested to ensure specificity and efficiency, and we used the primer pairs with the amplification efficiency of 1.9-2.1 for downstream analyses (Table 1).
Table 1. Primer sequences for the genes used in the real-time PCR.
Real-time PCR (RT-RCR) was conducted on Roche LightCycler® 96 detection system (Roche Applied Science, Germany) by using FastStart Essential DNA Green Master (Roche Applied Science, Germany). A 20 μl reaction mix was composed of 10 μl SYBR Green Master Mix, 8 μl PCR-free water, 0.5 μl of forward and reverse primers and 1 μl cDNA template. PCR program consisted of 40 cycles with 95°C for 10 s, 60°C for 10 s and 72°C for 15 s. Negative controls with non-template were used to eliminate contamination of the reaction system during the whole experiment process. Relative gene expression of treated groups to corresponding unstressed groups was normalized by the 2–ΔΔcq method (Pfaffl, 2001). For each gene, the transcriptional abundance of unstressed group was standardized as 1.0.
Data Analysis
All data including antioxidative indicators and gene expression was described as mean ± standard error. The significance test of all pairwise comparisons between treated and unstressed groups was analyzed using One-way ANOVA after normal test and homogeneity of variance by Kolmogorov-Smirnov test in SPSS (Version 25.0, SPSS Inc., United States). After One-way ANOVA with stress time as the factor for each treatment group, the significance test between unstressed and treated groups was assessed by Least Significant Difference (LSD) post hoc test. Principal Component Analysis (PCA) was also conducted using the IBM SPSS Statistics 25 and all statistical graphs were generated by OriginPro 9.1.
The antioxidant response element (ARE), as the main DNA binding sequence of Nrf2, was predicted within 2,000 bp upstream of the transcriptional start site of each gene with MatInspector tool in Genomatix Software Suite3. We focused on potential Nrf2 target genes recognized by ARE core sequence (5′-G/ATGACNNNGC-3′), and then described the number, distribution, and correlation with ARE for each gene.
Results
Antioxidative Indicators
General Change Patterns of Antioxidative Indicators
We first explored changes of antioxidative indicators by principal component analysis (PCA) to identify the major variation among all samples. Three axes explained 73.8% of the total variation in organismal antioxidant capacity (Figure 2A). The PCA loading plot showed that the ascidians collected in Round I did not separate from those sampled in Round II. The samples collected in two recovery phases were separated from those under cold stresses, but were not separated from samples collected under unstressed condition. These results illustrated that cold challenge was the dominant source of changes on antioxidative indicators, and C. robusta should use different antioxidative defense strategies at the stress and recovery stages.
Figure 2. Principal component analysis (PCA) based on antioxidant indicators (A) and gene expression data (B). Each data point represents a sample plotted using the first three principal components and the coordinate axis represents the first three principal component dimensions and their contribution to the overall variation among individuals. PCA was carried out on the log-transformed mean-centered data matrix.
Response of Organismal Antioxidative Status
Oxidative status of ascidians was assessed by MDA content and T-AOC. MDA content, which commonly indicates lipid peroxidation of the cell membrane, significantly increased at 1S-12 and 1R-6 and then returned to the unstressed level (Figure 3A). However, MDA content did not change in the subsequent round of cold stress. Our results showed a relatively stable status of MDA content without significant changes except for 1S-12, and such results suggest that C. robusta either be able to withstand the cold environment and/or have initiated ADS to maintain membrane integrity.
Figure 3. Organismal antioxidative level of Ciona robusta under unstressed (20°C), cold challenge (5°C) and recovery (20°C). (A) MDA (malondialdehyde) concentration; (B) T-AOC (total antioxidative capacity) activity. Each error bar represents the mean ± SE (N = 6). *p < 0.05, **p < 0.01.
Total antioxidant capacity is a comprehensive indicator of antioxidative ability. During the first cold stress, T-AOC significantly increased at 12 h, and then returned to the control level at the end of the first stress (Figure 3B). In the following first recovery, T-AOC significantly increased at 1R-6 and 1R-24. During the second round of cold stress-recovery, T-AOC significantly increased at most time points except 2S-1 (Figure 3B). However, when compared to the prestressed status, T-AOC did not significantly change during the second cold stress and recovery, suggesting that the prestressed T-AOC should have been optimized or served as a new baseline to cope with the repeated environmental changes (Supplementary Table S1).
Response of Organismal Antioxidative Indicators
Total superoxide dismutase activity only significantly increased at 1S-12 and 1R-6, and then remained at the unstressed level throughout all subsequent treatments (Figure 4A). CAT activity significantly decreased during cold exposures and got resilience at recovery stages. The minimum value of CAT activity in the two rounds dropped to 16 and 18% of the unstressed value at 1S-6 and 2S-24, respectively (Figure 4B). While during the recovery stages, CAT activity slightly recovered and returned to the unstressed level at the end of both rounds (Figure 4B). When comparing the response between two rounds of cold stress and recovery, we found that CAT activity slowly decreased in Round II with a delayed response, suggesting the increased tolerance in Round II.
Figure 4. Organismal antioxidative indicators of Ciona robusta under unstressed (20°C), cold challenge (5°C) and recovery (20°C). (A) T-SOD (total superoxide dismutase) activity; (B) CAT (catalase) activity; (C) reduced GSH (glutathione) concentration; Each error bar represents the mean ± SE (N = 6). In the second round, the significant difference between the unstressed and prestressed was represented by black and red color, respectively. *p < 0.05, **p < 0.01.
The concentration of GSH significantly increased during cold stresses while decreased at recovery stages in both rounds (Figure 4C). The GSH content increased with the stress duration and showed the highest value at the end of each stress (1S-24 and 2S-24) with similar fold change (3-fold). Whereas, the decline of GSH content during recovery stages did not show complete resilience but was still significantly higher than the unstressed level. In Round II, GSH content was maintained at a significantly elevated prestressed level and showed significant increase in 2S-24 relative to the prestressed status (Supplementary Table S1). All these results indicate that C. robusta re-adjusted the status with the elevated GSH content to antioxidative for subsequent cold challenges.
Given that the experimental temperature was sublethal for C. robusta (Dybern, 1965; Marin et al., 1987), the results here confirmed that the used conditions were not chronic stress but mutually independent repeated stresses. Hence, the steady antioxidative levels under the repeated stresses revealed the important role of the physiological plasticity of ADS for homeostasis.
Transcriptional Response of ADS-Related Genes
General Patterns of Transcriptional Changes
The gene expression data was analyzed using PCA to examine the general pattern of transcriptional changes. PCA results showed that three axes explained 87.7% of the total variation of transcriptional changes (Figure 2B). Similarly to the results on physiological changes, samples collected in cold stresses were separated from those in unstressed and recovery phases by PC2 (13.3%). Meanwhile, unstressed ascidians were clustered separately from those in recovery phases by PC3 (9.1%). Given the consideration between two stress phases, samples collected in the second stress were an overlap of those in the first stress with a more concentrated distribution. Such a pattern suggest that the response amplitude of the repeated stress should be narrower than that of the initial response and the recurrent response should be weak (i.e., cold stress memory).
Expression of Downstream Genes of Nrf2-keap1 Signal Pathway
Amongst the TFs (Nrf2, Keap1) and downstream genes (MnSOD, Cu/ZnSOD, CAT, GST, and GPx), the expression range of Round II was narrower than that of Round I, suggesting the decreased response after prestressing and the appearance of cold stress memory (Figures 5, 6A).
Figure 5. Relative expression level in pharynx muscle of Ciona robusta under unstressed (20°C), cold challenge (5°C) and recovery (20°C). (A) superoxide dismutase (SOD); (B) catalase (CAT); (C) glutathione S-transferases (GSTs); (D) glutathione peroxidase (GPX); (E) Nuclear factor (erythroid-derived 2)-like 2 (Nrf2); (F) Kelch-like ECH-associated protein 1 (Keap1). Each bar represents the mean ± SE (N = 6). In Figure 5A, significant difference of MnSOD and Cu/ZnSOD was represented by asterisk (*) and pound sign (#), respectively. In the second round, the significant difference between the unstressed and prestressed was represented by black and red color, respectively. */#p < 0.05, **/##p < 0.01.
Figure 6. Gene expression and correlation between Nrf2 and its downstream genes under recurrent environment stress. (A) Heat map of gene expression, Un = unstressed; (B) Pair-wise Pearson correlation coefficients (r) between transcriptional factors Nrf2 and its downstream genes. (C) The numbers and distribution of identified ARE (antioxidant-responsive element) by core sequence 5′-G/ATGACNNNGC-3′ in downstream genes. *p < 0.05, **p < 0.01.
Two isoforms of SOD, Cu/ZnSOD and MnSOD, showed similar response patterns (Figure 5A). The expression of both isoforms was tripled at 1S-6 and then sharply decreased at 1S-24. However, the two isoforms were merely upregulated to 1.1- and 1.2-fold in the second round of cold challenge, respectively. During recovery phases, the expression of these two SOD isoforms was significantly increased in the first 6 h (1R-1∼1R-6; 2R-1∼2R-6) and then declined at the end of each round. MnSOD was more sensitive and showed a robust response with a broader range than Cu/ZnSOD. Meanwhile, during the first round, CAT expression was significantly upregulated by 7. 5-, 17-, and 7.5-fold at 1S-1, 1S-6, and 1R-6, respectively, which was much higher than 1.5-fold at 2R-1 and 3-fold at 2R-6 (Figure 5B).
Two genes, GPx and GST, were selected to detect GSH metabolism. GPx expression showed 3-fold increase after 6-hour stress and then was significantly decreased to the value lower than that of the unstressed (Figure 5C). While in the second round, GPx was only increased to the unstressed level at 2S-6. The expression of GST was down-regulated during the whole experiment (Figure 5D). Meanwhile, the expression change range of the first round was between 0.1- and 0.8-fold, which were much wider than that of the second round (0.1- to 0.15-fold).
As the crucial molecular switch of Nrf2-keap1 signal pathway, the transcriptional expression of Nrf2 under cold exposure was sharply increased at 1S-1 and 1S-6 by 4- and 3-fold, respectively, and then significantly declined to the lowest value at 1S-24 (Figure 5E). During the first recovery stage, Nrf2 expression slowly recovered from the lowest value up to the similar value of the unstressed. The expression of Nrf2 was significantly decreased at 2S-1 and then continuously rose and stabilized to the prestress level. As the expression of Keap1 can be induced by the accumulation of Nrf2, similarly to Nrf2, the expression of Keap1 was only upregulated at 1S-1, 1S-6 and 2R-6 (Figure 5F).
In addition, given the comparison between the expression in the second round (six timepoints in total) and the unstressed, the expression of ADS-related genes at most timepoints (4/6 – 6/6) showed significant changes. When comparing the expression level of the second round to that of the prestressed, we found that expression of ADS-related genes showed significant changes at a fewer timepoints (1/6 - 5/6; Supplementary Table S2). For instance, expression of GST at all timepoints in the second round showed significant changes when compared to the unstressed value, while did not change when compared to the prestressed value except for 2S-24. Such a pattern demonstrated that the prestressing could remodel the baseline of gene expression, and the expression level during the second round was much closer to that after prestressing, but different from that of the unstressed.
Co-expression Relationship Between Nrf2 and Other Antioxidative Related Genes
By predicting the number of Antioxidant Response Element (ARE, DNA binding sites of Nrf2), we found that all the downstream genes on Nrf2-keap1 pathway had 4-6 ARE sites scattered along the 2000 bp upstream of each sequence (Figure 6C). Nrf2 expression was positively related to the expression level of Keap1, Cu/ZnSOD, CAT, GPx (r > 0.5, p < 0.05) and weakly related to MnSOD expression (r = 0.355, p = 0.032; Supplementary Table S3). This co-expression relationship between the transcription factor and its target genes indicated the potential “memory factor” role of Nrf2 to mediate the stress memory of its target genes in stress memory.
Discussion
Dynamic Response of the Antioxidant Defense System Under Cold Stress
Intense temperature changes such as sudden cold stress could largely disturb physiological conditions and metabolism of a species (Chen et al., 2019). Such adverse effects are highly related to the balance between generation and elimination of reactive oxygen species (ROS), which further triggers the activation of ADS (Sanchez-Nuno et al., 2018). Furthermore, cold-induced antioxidant defense largely relies on the dynamic changes of the activity and expression of ADS-related genes to cope with environmental changes (Wu et al., 2019). SOD is the first line in the elimination of ROS loads, since it can convert O2– to O2 and H2O2, which is subsequently transformed into H2O by CAT and GPx (Ma, 2013). Our results showed the transient activation of T-SOD, which was consistent with the excessive accumulation of MDA at the same time points (Figures 3A, 4A). The results here suggest that the activation of T-SOD should be a protective mechanism to prevent excessive oxidation of polyunsaturated fatty acids (PUFAs) and defend against further damages (Hofmann and Todgham, 2010; Kammer et al., 2011). The stable status in most timepoints might be the result of redox balance between superoxide production and antioxidant defense (Chen et al., 2019). At the transcriptional level, our results showed that the response pattern of different SOD isoforms was consistent (Figure 5A), while the response amplitude of MnSOD was greater than that of Cu/ZnSOD. These results suggest that ROS loads should be eliminated by the synthesis of SOD to protect tissues from oxidative damages, and MnSOD was the major functional isoform. The different response amplitudes of two SOD isoforms might be related to their various cellular distributions (Pietsch, 2017). MnSOD converts H2O2 in the mitochondrial matrix space, which is the main place producing ROS and the most sensitive position suffering from ROS effects (Murphy, 2009; Sies et al., 2017). Although Cu/ZnSOD has been approved as the major isoform (95% of the total SOD) in cells, it functions in intermembrane space and might play a weak role in cleaning ROS (Shadel and Horvath, 2015; Wu et al., 2019).
The catalyzation of CAT and biotransformation of GSH have complementary roles in the decomposition of H2O2 (Lesser, 2006; Sthijns et al., 2016). GSH and H2O2 could be catalyzed to GSSH and H2O by GPx (Zhang et al., 2019). In our study, CAT activity showed significant decrease under cold stresses, while the GSH content and the expression of GPx gene were elevated (Figures 4C, 5B,C). These results suggest that C. robusta might largely rely on GSH and GPx to scavenge H2O2 and lipid hydroperoxides. Additionally, we observed discordant response patterns of CAT on enzyme activity and gene expression regulatory levels (Figures 4B, 5B), which might suggest the potential effects of post-transcriptional and post-translational regulation of CAT in response to recurrent cold challenges. Since the half-life period of CAT mRNA is rather long (about 42 h in human cell and 30 h in saltwater clam Meretrix meretrix), some unidentified redox sensitive proteins can bind to the 5′ UTR of the CAT mRNA to regulate its translational rate and plus at the protein level, as observed in rat PC12 cancer cells exposed to H2O2 (Haque et al., 2012). It has been also reported that human ARPE-19 retinal pigment epithelial cells under several stressors including H2O2 can alter miRNA expression profiles for ROS, of which miR-30b can bind to the 3′ UTR of CAT mRNA, leading to a drastic decrease of CAT protein (Glorieux et al., 2015; Leisegang et al., 2018). Post-translational modifications, such as phosphorylation and ubiquitinylation, can regulate the turnover and degradation of CAT enzyme (Cao et al., 2003). In human cancer cells exposed to oxidative stress (H2O2), tyrosine kinases c-abl (abelson murine leukemia viral oncogene homolog 1) and c-abl-related gene (Arg) were able to phosphorylate human CAT enzyme at residues Tyr231 and Tyr386 to inhibit the CAT at the enzyme activity level. Therefore, although CAT gene expression was upregulated under cold stresses, the enzyme activity of CAT was still decreased possibly due to its translational repression induced by epigenetic changes or later biochemical processes.
The induction of ADS-related genes is regulated by several signaling pathways, in which Nrf2-Keap1 is a master regulatory pathway (Ma, 2013). In the present study, the elevation of Nrf2 expression showed an active response, which was further confirmed by the significant positive correlation between the expression of Nrf2 and its downstream genes (Figure 6B). Nevertheless, the continuous expression of Nrf2 could cause negative effects, such as free radical damages and cell apoptosis (Zheng et al., 2016). The elevation of Keap1 expression could activate the feedback regulation and increase Nrf2 degradation to control Nrf2 abundance (Zhang, 2006). Our results showed coordinated increase in the expression of Nrf2 and Keap1 (Figures 5E,F), supporting the role of keap1 in facilitating the termination of the pathway timely and preventing the overexpression of genes in Nrf2-keap1 signaling pathway. The positive response and feedback regulation of genes on Nrf2-keap1 signaling pathway ensured the ADS to maintain the balance of reactive oxygen free radicals through a series of changes and promote adaptive homeostasis (Kensler et al., 2007). In addition, Nrf2 could orchestrate the regulation of the gene expression involved in innate immunity, which confers protection by modulation of inflammatory defensive response against various environmental stresses (Battino et al., 2018). Cold stress can induce the upregulation of Nrf2 transcripts, which elicit antioxidant response to keep redox balance. Meanwhile, as the regulator of the innate immune response, Nrf2 can regulate the expression of the transcription factor in NF-κB signaling pathway to inhibit the transcriptional activity of NF-κB and protect animals against the adverse environment (Thimmulappa et al., 2006). By simultaneously modulating ADS and innate immune system, Nrf2 could enhance the rapid response for the stress protection and furthermore potentially invoke inflammatory defense to tissue integrity/repair under fluctuating environments.
Our results showed that after the transient increase at the expression level (1S-6 and 2S-6), most cold-induced gene expressions returned to renewed status (i.e., the rebuilt baselines) after a short period of time (i.e., transcriptional resilience; Figure 5). Transcriptional resilience was an indicator to predict stress tolerance (Franssen et al., 2011). For example, the study on the resilient patterns of recurrent bleaching on corals (genus Acropora) showed that the resilience capacity after bleaching strongly contributed to the sensitivity to future bleaching, and the species A. gemmifera with a fast recovery speed resulted in a high level of survival under repeated heat extremes (Thomas et al., 2019). Similarly, the transcriptional resilience in this study adds evidence on the adaptative response against recurrent stresses and should be functional in adaptative homeostasis simultaneously.
Cold Stress Memory Induced by Recurrent Stresses
Both the decreased response amplitude of ADS-related gene expression (Figure 5) and enhanced antioxidant capacity (Figure 3B) illustrated that C. robusta obtained cold stress memory to maintain adaptive homeostasis. Stress memory has been detected in multiple species and was characterized by decreased reactivity (tolerance) or enhanced response (potentiation) to subsequent challenges (Avramova, 2015). Study on murine embryonic fibroblasts (MEF) showed that pre-exposure to H2O2 could increase cell counts, BrdU incorporation and protein carbonyl content in subsequent challenges (Pickering et al., 2013). However, the transcriptional stress memory of A. thaliana was characterized by the expression decrease of osmotic and ionic equilibrium related genes under recurrent drought stresses (Ding et al., 2013). From a broad perspective, stress memory is a mechanism for organisms to store the message from an initial stress and becomes more robustness in response to future stresses, leading to the tolerating, priming, or acclimation affects after prestressing (Bruce et al., 2007).
Studies have proposed several molecular mechanisms underpinning stress memory. One mechanism was proposed to be sustained modulations of key signaling metabolites or transcription factors, which could initiate cascaded effects on downstream genes (Mozgova et al., 2019). For instance, recurring heat stresses induced the sustained expression of miR156, which regulated the transcriptional factor, Squamosa Promoter Binding Protein-Like (SPL), in A. thaliana, affecting endogenous timing of leaf initiation and developmental transitions (Stief et al., 2014). Another possible avenue was proposed to rely on chromatin modifications and epigenetic marks. Flowering Locus C (FLC) was transcriptionally repressed by epigenetic silence of H3K27me3 under cold exposure, and repression was then epigenetically consolidated during subsequent warmer temperature, facilitating a memory of cold challenges (Berry and Dean, 2015).
Nuclear factor erythroid-2-related factor 2 signaling pathway involves the maintenance of redox balance and stress response regulation, and the upregulated Nrf2 could result in transcriptional activation of its downstream genes to protect against various environmental stresses (Ma, 2013). The significant positive correlation between the expression of Nrf2 and ADS-related genes (Figure 6B; Supplementary Table S3), as well as the presence of ARE in the promoter regions of ADS-related genes (Figure 6C), illustrated the initiation of the Nrf2 pathway and implied the potential regulatory relationship between Nrf2 and other ADS-related genes. Considering sustained/accumulated expression levels of transcription factors for stress memory, the sustained changes of Nrf2 expression would potentially be related to the cold stress memory of ADS. For the potential mechanism of cold stress memory in this study, Nrf2 could act as the key transcription factor and potential “memory factor” whose sustained expression could regulate the transcriptional level of ADS-related genes during the subsequent cold exposure. The evidence obtained in this study provides baselines for further investigation of transcription factor patterns to reveal dynamic changes of ADS for comprehensive understanding of stress memory under environmental stresses.
The Effects of Cold Stress Memory
The acquisition of stress memory has several benefits, such as effectively resisting recurring environmental stresses and remaining homeostasis status (Ding et al., 2013; Avramova, 2015; Virlouvet and Fromm, 2015; Mozgova et al., 2019). Studies have approved that cells after prestressing were considerably more tolerant than non-prestressed cells to defend oxidative stresses in later challenges (Wiegant et al., 2011). Coral species experiencing heat stresses in previous years could reduce mortality during later heat shocks (Hughes et al., 2018; Harrison et al., 2019). Because organisms could modulate their internal states, such as energy allocation to activate/inactivate the expression of related genes and to coordinate physiological needs (López-Maury et al., 2008), the stress memory could shrink the energy allocation of ADS-related genes in recurrent stresses. Such a process allows prestressed individuals to cost less than those unstressed in response to future stresses. In the context of trade-off on resource allocation, stress memory could facilitate the expression of genes related to other functions such as growth and reproduction (van Hulten et al., 2006). Thus, stress memory provides an energy-efficient strategy to balance gene functions, which would possibly promote the homeostasis from rapid response to evolutionary adaptation (López-Maury et al., 2008).
Stress memory might promote the adaptive capacity to cope with harsh environments, making positive contributions to enlarging the ecological amplitude of a species. Similar to the function of feedback regulation in Nrf2-Keap1 signaling pathway, stress memory could limit the adverse impacts on overexpressed or persistent activation in the function of a protective effect against harsh environments (Melillo et al., 2018). Persistent adaptive response can be triggered by the duration of stresses, which can further increase the tolerance and make the physiological changes reduced gradually (potentially cause desensitization; Muniz-Castillo et al., 2019).
However, stress memory could lead to continuous resource consumption to maintain the readjusted baseline (i.e., stress memory can also entail costs; van Hulten et al., 2006). It is obvious that in an ambient condition (i.e., no disturbance), the ADS is always under dynamic balance and an organism has the highest fitness (Ristow, 2014). We assume that after prestressing, the redox balance would be no more disrupted until recurrent stresses appear, and the fitness of prestressed individuals would not be higher than that of unstressed ones. Unless the repeated challenges occur, organisms would undertake the costs for maintaining a remodeled baseline and merely gain growth benefits. Therefore, the recurring environmental pressures would be preconditions for an organism to get the opportunity to gain benefits from stress memory. Thus, when judging whether stress memory is the optimal strategy for survival of a species, we should consider the frequency of stress occurrence, as well as the associated benefits and costs (Maleck et al., 2000; van Hulten et al., 2006).
Potential Roles of Cold Stress Memory in Invasions
Our recurrent cold stress/recovery treatments on C. robusta simulated the fluctuating environments during the process of biological invasions. Usually, great and/or fast environmental challenge is an inescapable threat for invasive species. Stress memory might exert fundamental influence in response to future challenges, as the enduring effects of prior experience along with changing environments unfold over time (Hughes et al., 2018). At the early stages of invasions, C. robusta could both complete the geographical shifts and remain adaptive homeostasis by successfully invoking ADS to defense stresses. Additionally, organisms could effectively store “memory (rebuilding the baseline of ADS status)” from experienced stresses and optimize stress response at subsequent invasion stages. Stress memory could promote response by shrinking the amplitude of transcriptional changes, as well as by allocating limited resources to ADS related genes. In this perspective, stress memory can strongly impel to the process of later population establishment, as it probably tunes the resource allocation in favor of other functions such as growth and reproduction. We conclude that stress memory is a strategy to stabilize the status of ADS after multiple stresses, and it is imperative to lucubrate possible mechanisms of stress memory to deeply understand rapid adaptation during biological invasions.
Conclusion
Our results showed C. robusta relied on the resilience of ADS to maintain the adaptive homeostasis against environmental challenges. The resilience process of ADS transiently elevated the enzyme activities and transcription of ADS-related genes to readjusted baselines. The readjustment of ADS status after the initial stress was regarded as an active response to result in the stress memory. Subsequently, stress memory largely decreased the physiological response to recurrent environmental challenges. Overall, stress memory was an optimized strategy to allow C. robusta to store the “information” derived from initial invasion stages and decreased the subsequent responses to promote the homeostasis in the invasion process. Our results provide new insights into the mechanism of rapid physiological adaption under recurrent stresses and maintenance of adaptive homeostasis during biological invasions.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
Author Contributions
All authors conceived the study, designed the experiments, and took part in the interpretation of results and preparation of the manuscript. HL performed the experiments.
Funding
This work was supported by the National Natural Science Foundation of China (31272665 and 31772449).
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.
Acknowledgments
We thank Profs. Zunchun Zhou and Bei Jiang for assistance with ascidian sampling at the aquaculture farm and providing experiment sites for acclimation and treatments.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.00094/full#supplementary-material
Footnotes
- ^ http://asia.ensembl.org/index.html
- ^ http://www.premierbiosoft.com/primerdesign/index.html
- ^ http://www.genomatix.de
References
Avramova, Z. (2015). Transcriptional ‘memory’ of a stress: transient chromatin and memory (epigenetic) marks at stress-response genes. Plant J. 83, 149–159. doi: 10.1111/tpj.12832
Battino, M., Giampieri, F., Pistollato, F., Sureda, A., de Oliveira, M. R., Pittala, V., et al. (2018). Nrf2 as regulator of innate immunity: a molecular Swiss army knife! Biotechnol. Adv. 36, 358–370. doi: 10.10164/j.biotechadv.2017.12.012
Berry, S., and Dean, C. (2015). Environmental perception and epigenetic memory: mechanistic insight through FLC. Plant J. 83, 133–148. doi: 10.1111/tpj.12869
Briski, E., Bailey, S. A., Casas-Monroy, O., DiBacco, C., Kaczmarska, I., Lawrence, J. E., et al. (2013). Taxon- and vector-specific variation in species richness and abundance during the transport stage of biological invasions. Limnol. Oceanogr. 58, 1361–1372. doi: 10.4319/lo.2013.58.4.1361
Bruce, T. J. A., Matthes, M. C., Napier, J. A., and Pickett, J. A. (2007). Stressful “memories” of plants: evidence and possible mechanisms. Plant Sci. 173, 603–608. doi: 10.1016/j.plantsci.2007.09.002
Calabrese, E. J., and Mattson, M. P. (2017). How does hormesis impact biology, toxicology, and medicine? NPJ Aging Mech. Dis. 3:13. doi: 10.1038/s41514-017-0013-z
Calosi, P., Melatunan, S., Turner, L. M., Artioli, Y., Davidson, R. L., Byrne, J. J., et al. (2017). Regional adaptation defines sensitivity to future ocean acidification. Nat. Commun. 8:13994. doi: 10.1038/ncomms13994
Cao, C., Leng, Y., and Kufe, D. (2003). Catalase activity is regulated by c-Abl and Arg in the oxidative stress response. J. Biol. Chem. 278, 29667–29675. doi: 10.1074/jbc.M301292200
Carlson, S. M., Cunningham, C. J., and Westley, P. A. H. (2014). Evolutionary rescue in a changing world. Trends Ecol. Evol. 29, 521–530. doi: 10.1016/j.tree.2014.06.005
Carver, C., Mallet, A., and Vercaemer, B. (2006). Biological Synopsis of the Solitary Tunicate Ciona Intestinalis. Dartmouth: Bedford Institute of Oceanography.
Carver, C. E., Chisholm, A., and Mallet, A. L. (2003). Strategies to mitigate the impact of Ciona intestinalis (L.) biofouling on shellfish production. J. Shellfish Res. 22, 621–631.
Chen, B., Zhang, W., Niu, C., Li, W., Jia, H., and Storey, K. B. (2019). Antioxidant response to acute cold exposure and during recovery in juvenile Chinese soft-shelled turtles (Pelodiscus sinensis). J. Exp. Biol. 222:jeb197863. doi: 10.1242/jeb.197863
Chen, Y., Li, S., Lin, Y., Li, H., and Zhan, A. (2018a). Population genetic patterns of the solitary tunicate, Molgula manhattensis, in invaded Chinese coasts: large-scale homogeneity but fine-scale heterogeneity. Mar. Biodivers. 48, 2137–2149. doi: 10.1007/s12526-017-0743-y
Chen, Y., Shenkar, N., Ni, P., Lin, Y., Li, S., and Zhan, A. (2018b). Rapid microevolution during recent range expansion to harsh environments. BMC Evol. Biol. 18:187. doi: 10.1186/s12862-018-1311-1
Crisp, P. A., Ganguly, D., Eichten, S. R., Borevitz, J. O., and Pogson, B. J. (2016). Reconsidering plant memory: intersections between stress recovery, RNA turnover, and epigenetics. Sci. Adv. 2:e1501340. doi: 10.1126/sciadv.1501340
Davies, K. J. (2016). Adaptive homeostasis. Mol. Aspects Med. 49, 1–7. doi: 10.1016/j.mam.2016.04.007
Delsuc, F., Brinkmann, H., Chourrout, D., and Philippe, H. (2006). Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature 439, 965–968. doi: 10.1038/nature04336
Ding, Y., Liu, N., Virlouvet, L., Riethoven, J. J., Fromm, M., and Avramova, Z. (2013). Four distinct types of dehydration stress memory genes in Arabidopsis thaliana. BMC Plant Biol. 13:229. doi: 10.1186/1471-2229-13-229
Dybern, B. I. (1965). The life cycle of Ciona intestinalis (L.) f. typica in relation to the environmental temperature. Oikos 16, 109–131.
Franssen, S. U., Gu, J., Bergmann, N., Winters, G., Klostermeier, U. C., Rosenstiel, P., et al. (2011). Transcriptomic resilience to global warming in the seagrass Zostera marina, a marine foundation species. Proc. Natl. Acad. Sci. 108, 9276–9281. doi: 10.1073/pnas.1107680108
Fujikawa, T., Munakata, T., Kondo, S.-I., Satoh, N., and Wada, S. (2010). Stress response in the ascidian Ciona intestinalis: transcriptional profiling of genes for the heat shock protein 70 chaperone system under heat stress and endoplasmic reticulum stress. Cell Stress Chaperones 15, 193–204. doi: 10.1007/s12192-009-0133-x
Glorieux, C., Zamocky, M., Sandoval, J. M., Verrax, J., and Calderon, P. B. (2015). Regulation of catalase expression in healthy and cancerous cells. Free Radic. Biol. Med. 87, 84–97. doi: 10.1016/j.freeradbiomed.2015.06.017
Haque, R., Chun, E., Howell, J. C., Sengupta, T., Chen, D., and Kim, H. (2012). MicroRNA-30b-mediated regulation of catalase expression in human ARPE-19 cells. PLoS One 7:e42542. doi: 10.1371/journal.pone.0042542
Harrison, H. B., Alvarez-Noriega, M., Baird, A. H., Heron, S. F., MacDonald, C., and Hughes, T. P. (2019). Back-to-back coral bleaching events on isolated atolls in the Coral Sea. Coral Reefs 38, 713–719. doi: 10.1007/s00338-018-01749-6
Hoffmann, A. A., and Sgro, C. M. (2011). Climate change and evolutionary adaptation. Nature 470, 479–485. doi: 10.1038/nature09670
Hofmann, G. E., and Todgham, A. E. (2010). Living in the now: physiological mechanisms to tolerate a rapidly changing environment. Annu. Rev. Physiol. 72, 127–145. doi: 10.1146/annurev-physiol-021909-135900
Huang, X., Gao, Y., Jiang, B., Zhou, Z., and Zhan, A. (2016). Reference gene selection for quantitative gene expression studies during biological invasions: a test on multiple genes and tissues in a model ascidian Ciona savignyi. Gene 576, 79–87. doi: 10.1016/j.gene.2015.09.066
Huang, X., Li, S., Gao, Y., and Zhan, A. (2018). Genome-wide identification, characterization and expression analyses of heat shock protein-related genes in a highly invasive ascidian Ciona savignyi. Front. Physiol. 9:1043. doi: 10.3389/fphys.2018.01043
Huang, X., Li, S., Ni, P., Gao, Y., Jiang, B., Zhou, Z., et al. (2017). Rapid response to changing environments during biological invasions: DNA methylation perspectives. Mol. Ecol. 26, 6621–6633. doi: 10.1111/mec.14382
Huang, X., Li, S., and Zhan, A. (2019). Genome-wide identification and evaluation of new reference genes for gene expression analysis under temperature and salinity stresses in Ciona savignyi. Front. Genet. 10:71. doi: 10.3389/fgene.2019.00071
Hughes, T. P., Kerry, J. T., Connolly, S. R., Baird, A. H., Eakin, C. M., Heron, S. F., et al. (2018). Ecological memory modifies the cumulative impact of recurrent climate extremes. Nat. Clim. Chang. 9, 40–43. doi: 10.1038/s41558-018-0351-2
Jimenez, A. G., Jayawardene, S., Alves, S., Dallmer, J., and Dowd, W. W. (2015). Micro-scale environmental variation amplifies physiological variation among individual mussels. Proc. Biol. Sci. 282:20152273. doi: 10.1098/rspb.2015.2273
Kammer, A. R., Orczewska, J. I., and O’Brien, K. M. (2011). Oxidative stress is transient and tissue specific during cold acclimation of threespine stickleback. J. Exp. Biol. 214, 1248–1256. doi: 10.1242/jeb.053207
Kensler, T. W., Wakabayash, N., and Biswal, S. (2007). Cell survival responses to environmental stresses via the Keap1-Nrf2-ARE pathway. Annu. Rev. Pharmacol. Toxicol. 47, 89–116. doi: 10.1146/annurev.pharmtox.46.120604.141046
Kishimoto, I., Ariga, I., Itabashi, Y., and Mikami, K. (2019). Heat-stress memory is responsible for acquired thermotolerance in Bangia fuscopurpurea. J. Phycol. 55, 971–975. doi: 10.1111/jpy.12895
Leisegang, M. S., Schroder, K., and Brandes, R. P. (2018). Redox regulation and noncoding RNAs. Antioxid. Redox Sign. 29, 793–812. doi: 10.1089/ars.2017.7276
Lenz, M., da Gama, B. A., Gerner, N. V., Gobin, J., Groner, F., Harry, A., et al. (2011). Non-native marine invertebrates are more tolerant towards environmental stress than taxonomically related native species: results from a globally replicated study. Environ. Res. 111, 943–952. doi: 10.1016/j.envres.2011.05.001
Lesser, M. P. (2006). Oxidative stress in marine environments: biochemistry and physiological ecology. Annu. Rev. Physiol. 68, 253–278. doi: 10.1146/annurev.physiol.68.040104.110001
López-Maury, L., Marguerat, S., and Bähler, J. (2008). Tuning gene expression to changing environments: from rapid responses to evolutionary adaptation. Nat. Rev. Genet. 9, 583–593. doi: 10.1038/nrg2398
Ma, Q. (2013). Role of nrf2 in oxidative stress and toxicity. Annu. Rev. Pharmacol. Toxicol. 53, 401–426. doi: 10.1146/annurev-pharmtox-011112-140320
Maleck, K., Levine, A., Eulgem, T., Morgan, A., Schmid, J., Lawton, K. A., et al. (2000). The transcriptome of Arabidopsis thaliana during systemic acquired resistance. Nat. Genet. 26, 403–410. doi: 10.1038/82521
Marin, M. G., Bressan, M., Beghi, L., and Brunetti, R. (1987). Thermohaline tolerance of ciona intestinalis l. 1767 at different developmental stages. Cah. Biol. Mar. 28, 47–58.
Melillo, D., Marino, R., Italiani, P., and Boraschi, D. (2018). Innate immune memory in invertebrate metazoans: a critical appraisal. Front. Immunol. 9:1915. doi: 10.3389/fimmu.2018.01915
Moreira, A., Figueira, E., Soares, A. M., and Freitas, R. (2016). The effects of arsenic and seawater acidification on antioxidant and biomineralization responses in two closely related Crassostrea species. Sci. Total Environ. 54, 569–581. doi: 10.1016/j.scitotenv.2015.12.029
Mozgova, I., Mikulski, P., Pecinka, A., and Farrona, S. (2019). “Epigenetic mechanisms of abiotic stress response and memory in plants,” in Epigenetics in Plants of Agronomic Importance: Fundamentals and Applications, eds R. Alvarez-Venegas, C. De-la-Peña, and J. Casas-Mollano (Cham: Springer), 1–64. doi: 10.1007/978-3-030-14760-0_1
Muniz-Castillo, A. I., Rivera-Sosa, A., Chollett, I., Eakin, C. M., Andrade-Gomez, L., McField, M., et al. (2019). Three decades of heat stress exposure in caribbean coral reefs: a new regional delineation to enhance conservation. Sci. Rep. 9:11013. doi: 10.1038/s41598-019-47307-0
Murphy, M. P. (2009). How mitochondria produce reactive oxygen species. Biochem. J. 417, 1–13. doi: 10.1042/Bj20081386
Nandamuri, S. P., Yourick, M. R., and Carleton, K. L. (2017). Adult plasticity in African cichlids: rapid changes in opsin expression in response to environmental light differences. Mol. Ecol. 26, 6036–6052. doi: 10.1111/mec.14357
Nguyen, T., Sherratt, P. J., and Pickett, C. B. (2003). Regulatory mechanisms controlling gene expression mediated by the antioxidant response element. Annu. Rev. Pharm. Toxicol. 43, 233–260. doi: 10.1146/annurev.pharmtox.43.100901.140229
Ni, P., Li, S., Lin, Y., Xiong, W., Huang, X., and Zhan, A. (2018). Methylation divergence of invasive Ciona ascidians: significant population structure and local environmental influence. Ecol. Evol. 8, 10272–10287. doi: 10.1002/ece3.4504
Oliveira, M. F., Geihs, M. A., Thiago, F. A. F., Moreira, D. C., and Hermes-Lima, M. (2018). Is “preparation for oxidative stress” a case of physiological conditioning hormesis? Front. Physiol. 9:945. doi: 10.3389/fphys.2018.00945
Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29:e45. doi: 10.1093/nar/29.9.e45
Pickering, A. M., Vojtovich, L., Tower, J., and Davies, K. J. A. (2013). Oxidative stress adaptation with acute, chronic, and repeated stress. Free Radic. Biol. Med. 55, 109–118. doi: 10.1016/j.freeradbiomed.2012.11.001
Pietsch, C. (2017). Zearalenone (ZEN) and Its influence on regulation of gene expression in carp (Cyprinus carpio L.) liver tissue. Toxins 9:E283. doi: 10.3390/toxins9090283
Richards, C. L., Bossdorf, O., Muth, N. Z., Gurevitch, J., and Pigliucci, M. (2006). Jack of all trades, master of some? On the role of phenotypic plasticity in plant invasions. Ecol. Lett. 9, 981–993. doi: 10.1111/j.1461-0248.2006.00950.x
Ristow, M. (2014). Unraveling the truth about antioxidants: mitohormesis explains ROS-induced health benefits. Nat. Med. 20, 709–711. doi: 10.1038/nm.3624
Rius, M., Clusella-Trullas, S., McQuaid, C. D., Navarro, R. A., Griffiths, C. L., Matthee, C. A., et al. (2014). Range expansions across ecoregions: interactions of climate change, physiology and genetic diversity. Glob. Ecol. Biogeogr. 23, 76–88. doi: 10.1111/geb.12105
Sanchez-Nuno, S., Sanahuja, I., Fernandez-Alacid, L., Ordonez-Grande, B., Fontanillas, R., Fernandez-Borras, J., et al. (2018). Redox challenge in a cultured temperate marine species during low temperature and temperature recovery. Front. Physiol. 9:923. doi: 10.3389/fphys.2018.00923
Schieber, M., and Chandel, N. S. (2014). ROS function in redox signaling and oxidative stress. Curr. Biol. 24, R453–R462. doi: 10.1016/j.cub.2014.03.034
Seebacher, F., White, C. R., and Franklin, C. E. (2014). Physiological plasticity increases resilience of ectothermic animals to climate change. Nat. Clim. Chang. 5, 61–66. doi: 10.1038/nclimate2457
Shadel, G. S., and Horvath, T. L. (2015). Mitochondrial ROS signaling in organismal homeostasis. Cell 163, 560–569. doi: 10.1016/j.cell.2015.10.001
Sies, H., Berndt, C., and Jones, D. P. (2017). Oxidative stress. Annu. Rev. Biochem. 86, 715–748. doi: 10.1146/annurev-biochem-061516-0-45037
Smith, K. F., Thia, J., Gemmill, C. E. C., Cary, S. C., and Fidler, A. E. (2012). Barcoding of the cytochrome oxidase I (COI) indicates a recent introduction of Ciona savignyi into New Zealand and provides a rapid method for Ciona species discrimination. Aquat. Invasion. 7, 305–313. doi: 10.3391/ai.2012.7.3.002
Somero, G. N. (2012). The physiology of global change: linking patterns to mechanisms. Ann. Rev. Mar. Sci. 4, 39–61. doi: 10.1146/annurev-marine-120710-100935
Sparks, C., Marnewick, J., Toefy, R., Snyman, R., and Odendaal, J. (2019). Baseline levels of antioxidant activities in Mytilus galloprovincialis along the coast of Cape Town, South Africa. Mar. Pollut. Bull. 140, 287–293. doi: 10.1016/j.marpolbul.2019.01.059
Spicer, J. I., Morley, S. A., and Bozinovic, F. (2019). Physiological diversity, biodiversity patterns and global climate change: testing key hypotheses involving temperature and oxygen. Philos. Trans. R. Soc. B Biol. Sci. 374:20190032. doi: 10.1098/rstb.2019.0032
Sthijns, M. M., Weseler, A. R., Bast, A., and Haenen, G. R. (2016). Time in redox adaptation processes: from evolution to hormesis. Int. J. Mol. Sci. 17:E1649. doi: 10.3390/ijms17101649
Stief, A., Altmann, S., Hoffmann, K., Pant, B. D., Scheible, W. R., and Baurle, I. (2014). Arabidopsis miR156 regulates tolerance to recurring environmental stress through spl transcription factors. Plant Cell 26, 1792–1807. doi: 10.1105/tpc.114.123851
Storey, K. B. (1996). Oxidative stress: animal adaptations in nature. Braz. J. Med. Biol. Res. 29, 1715–1733.
Thimmulappa, R. K., Lee, H., Rangasamy, T., Reddy, S. P., Yamamoto, M., Kensler, T. W., et al. (2006). Nrf2 is a critical regulator of the innate immune response and survival during experimental sepsis. J Clin. Invest. 116, 984–995. doi: 10.1127/Jci25790
Thomas, L., Lopez, E. H., Morikawa, M. K., and Palumbi, S. R. (2019). Transcriptomic resilience, symbiont shuffling, and vulnerability to recurrent bleaching in reef-building corals. Mol. Ecol. 28, 3371–3382. doi: 10.1111/mec.15143
van Hulten, M., Pelser, M., van Loon, L. C., Pieterse, C. M., and Ton, J. (2006). Costs and benefits of priming for defense in arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 103, 5602–5607. doi: 10.1073/pnas.0510213103
Virlouvet, L., and Fromm, M. (2015). Physiological and transcriptional memory in guard cells during repetitive dehydration stress. New Phytol. 205, 596–607. doi: 10.1111/nph.13080
Wiegant, F. A., Prins, H. A., and Van Wijk, R. (2011). Postconditioning hormesis put in perspective: an overview of experimental and clinical studies. Dose Response 9, 209–224. doi: 10.2203/dose-response.10-004.Wiegant
Wu, H., Liu, Y., Zhang, X., Zhang, J., and Ma, E. (2019). Antioxidant defenses at enzymatic and transcriptional levels in response to acute lead administration in Oxya chinensis. Ecotox. Environ. Safe 168, 27–34. doi: 10.1016/j.ecoenv.2018.10.061
Zhan, A., Briski, E., Bock, D. G., Ghabooli, S., and MacIsaac, H. J. (2015). Ascidians as models for studying invasion success. Mar. Biol. 162, 2449–2470. doi: 10.1007/s00227-015-2734-5
Zhan, A., Macisaac, H. J., and Cristescu, M. E. (2010). Invasion genetics of the Ciona intestinalis species complex: from regional endemism to global homogeneity. Mol. Ecol. 19, 4678–4694. doi: 10.1111/j.1365-294X.2010.04837.x
Zhang, D. (2006). Mechanistic studies of the Nrf2-Keap1 signaling pathway. Drug Metab. Rev. 38, 769–789. doi: 10.1080/03602530600971974
Zhang, W., Chen, B., Niu, C., Yuan, L., Jia, H., and Storey, K. B. (2019). Response of the chinese soft-shelled turtle to acute heat stress: insights from the systematic antioxidant defense. Front. Physiol. 10:710. doi: 10.3389/fphys.2019.00710
Zheng, J., Zeng, L., Shen, B., Xu, M., Zhu, A., and Wu, C. (2016). Antioxidant defenses at transcriptional and enzymatic levels and gene expression of Nrf2-Keap1 signaling molecules in response to acute zinc exposure in the spleen of the large yellow croaker Pseudosciaena crocea. Fish Shellfish Immunol. 52, 1–8. doi: 10.1016/j.fsi.2016.02.03-1
Keywords: physiology plasticity, environmental stress, ascidian, antioxidant defense system, stress memory, adaptive homeostasis
Citation: Li H, Huang X and Zhan A (2020) Stress Memory of Recurrent Environmental Challenges in Marine Invasive Species: Ciona robusta as a Case Study. Front. Physiol. 11:94. doi: 10.3389/fphys.2020.00094
Received: 23 October 2019; Accepted: 27 January 2020;
Published: 13 February 2020.
Edited by:
Youji Wang, Shanghai Ocean University, ChinaReviewed by:
Mohamed Abdelmoneim Mostafa Hegazi, Tanta University, EgyptXiaoli Hu, Ocean University of China, China
Copyright © 2020 Li, Huang and Zhan. 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: Aibin Zhan, zhanaibin@hotmail.com; azhan@rcees.ac.cn