Original Research ARTICLE
Air Exposure Affects Physiological Responses, Innate Immunity, Apoptosis and DNA Methylation of Kuruma Shrimp, Marsupenaeus japonicus
- 1State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
- 2Fujian Key Laboratory of Genetics and Breeding of Marine Organisms, Xiamen University, Xiamen, China
Air exposure stress is a common phenomenon for commercial crustacean species in aquaculture and during waterless transportation. However, the antioxidant responses to air exposure discussed in previous studies may be insufficient to present the complexities involved in this process. The comprehensive immune responses, especially considering the immune genes, cell apoptosis, and epigenetic changes, are still unknown. Accordingly, we investigated the multifaceted responses of Marsupenaeus japonicus to air exposure. The results showed that the expression profiles of the apoptosis genes (e.g., IAP, TXNIP, caspase, and caspase-3) and the hypoxia-related genes (e.g., hsp70, hif-1α, and HcY) were all dramatically induced in the hepatopancreas and gills of M. japonicus. Heart rates, T-AOC (total antioxidant capacity) and lactate contents showed time-dependent changes upon air exposure. Air exposure significantly induced apoptosis in hepatopancreas and gills. Compared with the control group, the apoptosis index (AI) of the 12.5 h experimental group increased significantly (p < 0.05) in the hepatopancreas and gills. Most individuals in the experimental group (EG, 12.5 h) had lower methylation ratios than the control group (CG). Air exposure markedly reduced the full-methylation and total-methylation ratios (31.39% for the CG and 26.46% for the EG). This study provided a comprehensive understanding of the antioxidant responses of M. japonicus considering its physiology, innate immunity, apoptosis, and DNA methylation levels, and provided theoretical guidance for waterless transportation.
The Kuruma shrimp, Marsupenaeus japonicus, is an important commercial species that is widely distributed in the Indo-Western Pacific (Tsoi et al., 2014). Compared to other shrimps, M. japonicus have numerous advantages, such as fast growth rates, good reproductive capability, and survival during waterless transportation (Duan et al., 2016a). Air exposure stress is a common phenomenon for commercial crustacean species in aquaculture and during live transport (Lorenzon et al., 2008; Fotedar and Evans, 2011; Romero et al., 2011). Upon air exposure, aquatic animals are subject to water shortage and hypoxia stress, accompanied by the generation of reactive oxygen species (ROS) (Romero et al., 2007; Paital and Chainy, 2010; Paital, 2013; Abasubong et al., 2018). Air exposure affects cellular damage and metabolic capacity for Mytilus galloprovincialis (Andrade et al., 2019) and causes a significant but reparable immunological response in Litopenaeus vannamei (Xu et al., 2019). Excessive ROS have an obvious oxidation effect on lipids, nucleic acids, and proteins that leads to an imbalance of homeostasis and physiological metabolism (Romero et al., 2011; Lennicke et al., 2015; Sies, 2017). Continuous oxidative stress triggers caspase-independent cell death and induces apoptosis or necrosis (Simon et al., 2000; Wang et al., 2012, 2019; Nathan and Cunningham-Bussel, 2013; Holze et al., 2018; Ondricek and Thomas, 2018).
Crustaceans mainly rely on non-specific innate immunity, including physical defense, and cellular and humoral immunity (Sadaaki and Bok Luel, 2005; Li and Xiang, 2013). To eliminate ROS-mediated damages, organisms protect themselves with enzymatic and non-enzymatic antioxidant defenses (Storey, 1996). The antioxidant enzyme systems mainly include superoxide dismutase (SOD), glutathione peroxidase (GPx) and catalase (CAT) (Ighodaro and Akinloye, 2018). Duan et al. (2016a) investigated the effect of air exposure on antioxidant enzyme activities of Kuruma shrimp and found that the activities of SOD and CAT increased significantly. T-AOC (total antioxidant capacity), including enzymatic and non-enzymatic antioxidants, reflects the comprehensive antioxidant capacity of an organism (Jia et al., 2011; Liu H.L. et al., 2015). Air exposure induced significant oxidative and antioxidant responses in L. vannamei (Xu et al., 2018). Other research has detected significant changes in antioxidant enzyme activities in Lithodes santolla (Schvezov et al., 2019), Paralomis granulosa (Romero et al., 2011), Apostichopus japonicus (Huo et al., 2018), Penaeus monodon (Duan et al., 2016b), Babylonia areolate (Liu et al., 2017), and Oryzias melastigma (Cui et al., 2019). Hypoxia exposure had significant effects on the expression profiles of apoptosis and hypoxia-related genes, such as p53, Hif-1α, Hsp70, GSH-Px (glutathione peroxidase) and CAT (catalase) (Xiao, 2015; Sun et al., 2016a; Cai et al., 2018; Ondricek and Thomas, 2018; Camacho-Jiménez et al., 2019). Hypoxia-inducible factors (HIFs) are recognized as master regulators of the cellular response to hypoxia stress by binding HRE (hypoxia-responsive element) (Schofield and Ratcliffe, 2004; Choudhry and Harris, 2018). The role of scHIF-1α in the regulation of the HIF signaling pathway and the immune response was confirmed in mandarin fish (Siniperca chuatsi) (He et al., 2019). Cancer magister responds to hypoxia by increasing Hc concentrations and regulating subunit compositions (Head, 2010). Heat shock protein 70 (Hsp70) plays a vital role in preventing protein aggregation, folding and refolding of proteins and the degradation of unstable and misfolded proteins (Daugaard et al., 2007). The thioredoxin system, including TRX, TRXR, and NADPH, plays a critical role in resisting oxidative stress, which has been suggested to have complex regulatory functions in shrimp (Garcia-Orozco et al., 2012; Lu and Holmgren, 2014; Liu P.-F. et al., 2015; Zuo et al., 2019). DNA methylation, a conservative method of epigenetic regulation in eukaryotes, regulates many biological processes such as transcription, DNA repair and cell differentiation (Zemach et al., 2010; Sahu et al., 2013; Wang et al., 2015). DNA methylation levels are usually influenced by various biotic and abiotic environmental stresses, including oxidative stress (Franco et al., 2008; Murgatroyd et al., 2009; Dowen et al., 2012).
In general, air exposure is a comprehensive stimulation to aquatic animals. However, most previous studies have mainly focused on antioxidant enzyme activities. The comprehensive immune responses, especially immune genes, cell apoptosis, and epigenetic changes, are still unknown. Our objectives in the study were to determine the effects of air exposure on the physiological responses, non-specific immunity, apoptosis and DNA methylation levels of M. japonicus.
Materials and Methods
All healthy M. japonicus individuals (weight: 6.25 ± 1.3 g) were obtained from an aquaculture farm in Zhangzhou (Fujian, China) and were acclimated in environmentally controlled cement pools for 2 weeks (25°C, 29‰ salinity and under continuous aeration). The filtered seawater was renewed every other day and the shrimps were fed twice daily with moderate amounts of commercial pelleted food. Shrimps were fasted for 24 h before the stress treatment.
Air Exposure and Sampling
After the conditioning period, two hundred individuals were transferred into an air-conditioned room at 25°C, and an air-exposure challenge was performed as follows. For the air-exposure treatment, the water on the surfaces of all individuals was sucked up using gauze and all individuals were then placed in six flat-bottomed rectangular tanks (80 cm × 50 cm × 30 cm). Every 20 min, these shrimp were sprayed with seawater to keep the skin moist. During the air exposure period, the hepatopancreas, gills and muscles from nine individuals were sampled at control (unstressed), 2.5, 5, 7.5, 10, and 12.5 h after exposure. Further, we recorded the numbers of dead individuals and counted the cumulative mortality. The criterion determining death was the loss of swimming ability after resubmersion. The heart rates of six individuals were determined by direct counting for 2 min at each sample point. The hepatopancreas and gills were immediately preserved in RNAfixer (BioTeck, China) for RNA extraction. Parts of the gill and hepatopancreas were rinsed with normal saline and then immediately preserved in 4% paraformaldehyde for detecting cell apoptosis. Part of the muscle and hepatopancreas tissues were immediately frozen at −80°C for lactate content determination and methylation-sensitive amplification polymorphism (MSAP) analysis, respectively.
Homology Cloning and Sequence Analysis
Total RNA was extracted using RNAiso Plus (TaKaRa, Japan) following the specification. The purity and integrity were assessed. First-strand cDNA was synthesized using TransScript II First-Strand cDNA SuperMix (TransGen Biotech, China) according to the manufacture’s protocol. Primers were designed based on our transcriptome unigenes and sequences of other species that have been uploaded to NCBI (Table 1). PCR amplification products were ligated into the pMD19-T simple vector and sequenced in both directions. The open reading frames (ORFs) of IAP (inhibitor of apoptosis) and TXNIP (thioredoxin interacting protein) were predicted using the ORF finder1. The amino acid sequence was deduced by the EMBOSS Transeq2. ExPASy3 was used to analyze the protein molecular mass and isoelectric point (pI). The protein domains were predicted by SMART software4. Sequence alignments and phylogenetic trees were performed using MEGA 7.0 software (Kumar et al., 2016).
The Quantitative Real-Time PCR Analysis
The expression profiles of apoptosis and hypoxia-related genes in the hepatopancreas and gills were measured by qRT-PCR, which was implemented on the QuantStudio 6 Flex Real-time PCR System (Applied Biosystems, United States) using SYBR® Premix ExTaq II (2×) (TaKaRa, Japan) following the protocols. All samples were run in triplicate and the M. japonicus elongation factor 1-α (EF1-α) was served as the reference gene (Table 1). Gene expression change was calculated by the 2–ΔΔCt method (Livak and Schmittgen, 2001). Expression data were shown with the mean ± standard deviation.
Total Antioxidant Capacity (T-AOC) and Lactate Content Determination
The total antioxidant capacity (T-AOC) of the hepatopancreas tissues was determined spectrophotometrically at 520 nm according to the instructions of a commercial assay kit (Jiancheng, China). Muscle tissues were used for lactate content determination at 530 nm using a lactic acid assay kit (Jiancheng, China). The determinations were conducted in triplicate for each sample using the Infinite M200 Pro system (Tecan, Switzerland) following the protocols. One unit of T-AOC was defined as the amount of protein per mg needed to increase 0.01 of the absorbance values every minute under the assay conditions (U/mg protein). The lactate content was expressed as mmol/g protein.
DAPI Staining and TUNEL Assay
The hepatopancreas and gills from five individuals in each group (control and 12.5 h) were sampled and preserved in 4% paraformaldehyde for 24 h. The fixed tissues were embedded in paraffin following ethanol dehydration and vitrification in dimethylbenzene. The paraffin sections were dewaxed, repaired and ruptured. The number of apoptotic cells in hepatopancreas and gills was determined using a One-Step TUNEL apoptosis assay kit (Beyotime, China) following the manufacturer’s instructions. Apoptotic cell nuclei and total cell nuclei were quantitated by TUNEL staining (green fluorescence) and DAPI staining (blue fluorescence), respectively. The double-stained sections were examined using an inverted fluorescence microscope (NIKON Eclipse TI-SR, Japan). For each sample, four visual fields were selected randomly and were photographed by an inverted fluorescence microscope. Image-Pro Plus 6.0 software was used to count the number of apoptotic and common cells. The apoptosis index (AI) was defined as the proportion of apoptotic cells in the total cells.
DNA Methylation Levels and Patterns Under Air Exposure
Using the phenol-chloroform method, total genomic DNA was extracted from the hepatopancreas tissues of samples from the control and at 12.5 h after air exposure. The DNA integrity was visualized by agarose gel electrophoresis (1%) and concentration was measured using the Qubit 2.0 fluorometer (Life Technologies, CA, United States). DNA samples with clear bands were diluted to 100 ∼ 200 ng/μL for the downstream analysis. The MSAP analysis was performed according to Xiong et al. (1999) and He et al. (2015) with modifications. The enzyme-digestion products of EcoRI/HpaII and EcoRI/MspI were ligated to adapters with T4 DNA ligase (TaKaRa) and performed pre-amplification and selective amplification with fluorescent-labeled primers (Table 1).
The selective amplification products were detected by capillary electrophoresis (CE) on an ABI 3730XL auto DNA sequencer (Shanghai Sangon, China). The MSAP fragments were analyzed using GeneMapper v3.2 software. Different enzyme digestion combinations of the same DNA sample were denoted as EH/EM. Three kinds of bands were detected and counted, including EH + /EM + (type I), EH + /EM- (type II) and EH-/EM + (type III). The methylation ratio was calculated as the following formula: Hemi-methylation ratio = Type II/(Type I + II + III) and Full-methylation ratio = Type III/(Type I + II + III).
The data were analyzed using SPSS Statistics version 22 (IBM, United States). The results were presented as the means ± SD. The significant difference levels were detected by one-way analysis of variance (ANOVA) followed by Tukey’s test, and a p-value less than 0.05 (<0.05) is statistically significant.
Sequence Analysis of TXNIP and IAP
The open reading frame (ORF) of MjTXNIP was 1044 bp (GenBank accession number: MN265397), encoding a polypeptide of 347 amino-acid residues with a calculated molecular weight of 39.97 kDa and a theoretical pI of 8.61. SMART analysis revealed that the predicted TXNIP protein contained an arrestin C domain located at 184-315 residues with an expected value of 2.44e-23 (Figure 1A). Homology comparisons showed that TXNIP shared amino acid identities with TXNIP from L. vannamei (98.85%), Daphnia magna (75.43%), Limulus polyphemus (65.99%), Drosophila melanogaster (61.73%), Danio rerio (22.67%), and Homo sapiens (21.52%). The open reading frame (ORF) of MjIAP (GenBank accession number: MN265396) was 2136 bp, encoding a polypeptide of 711 amino-acid residues with a calculated molecular weight of 78.48 kDa and a theoretical pI of 5.71. The SMART analysis revealed that the predicted IAP protein contains three BIR domains and a RING domain, which are located at 17-82, 108-172, 257-322, and 664-699 residues, respectively (Figure 1B). Multiple sequence alignment showed that MjIAP shared high amino-acid identity with IAP from P. monodon (85.96%, ABO38431.1), L. vannamei (84.83%, ADH03018.1) and Scylla paramamosain (59.73%, AMW91738.1). Phylogenetic analysis showed that the MjIAP clustered with the IAP protein of other crustaceans with a bootstrapping value of 100.
Figure 1. Amino acid sequence analysis of Marsupenaeus japonicus TXNIP and IAP. (A) Neighbor joining phylogenetic tree and protein domains analysis of TXNIPs. Homo sapiens (NP_006463.3), Danio rerio (AAH67376.1), Drosophila melanogaster (NP_612004.2), Cyprinus carpio (XP_018951332.1), Litopenaeus vannamei (AYE56137.1), Limulus polyphemus (XP_013787968.1), Tribolium castaneum (XP_972691.1), Nasonia vitripennis (XP_001604921.2), Myzus persicae (XP_022178427.1), Mus musculus (NP_001009935.1), Daphnia magna (KZS06009.1), Drosophila melanogaster (NP_612004.2). (B) Neighbor joining phylogenetic tree and protein domains analysis of IAPs. Mus musculus (AAC53531.1), Homo sapiens (AAC50371.1), Danio rerio (AAP04483.1), Bombyx mori (XP_012551311.1), Eurytemora affinis (XP_023329578.1), Frankliniella occidentalis (XP_026285387.1), Hyalella azteca (XP_018007268.1), Mizuhopecten yessoensis (XP_021344767.1), Palaemon carinicauda (AHC55308.1), Scylla paramamosain (AMW91738.1), Penaeus monodon (ABO38431.1), Litopenaeus vannamei (ADH03018.1).
Expression Profiles of MjIAP, MjTXNIP, hsp70, hif-1α, HcY, caspase, and caspase3 Under Air Exposure Stress
The expression profiles of the apoptosis- and the hypoxia-related genes are shown in Figures 2, 3. Overall, these genes showed time-dependent changes after 5 h post-exposure in the hepatopancreas and after 2.5 h in the gill tissues. The mRNA expression levels of MjIAP gradually increased with time during air exposure. The expression levels in the hepatopancreas showed 2. 34-, 4. 3-, and 9.88-fold increases (p < 0.05) at 7.5, 10, and 12.5 h post-exposure, respectively (Figure 2A). In the gills, the mRNA expression levels of MjIAP showed 2.4- and 3.8-fold increases (p < 0.05) at 2.5 and 5 h, respectively (Figure 3A). In the hepatopancreas, the transcriptional level of MjTXNIP at the 10 h time point showed a significant incremental change (p < 0.05) (Figure 2B), while the expression levels in the gills peaked at 12.5 h with a 1.38-fold increase (p < 0.05) (Figure 3B). In the hepatopancreas, the mRNA expression levels of caspase increased 1.79-fold at the 5 h time point (p < 0.05) and peaked at the 7.5 h time point with a 1.91-fold increase (Figure 2C). In the gills, the mRNA levels showed a significant incremental change at the 10 h time point with a 1.44-fold increase (Figure 3C). The caspase-3 levels in the hepatopancreas were significantly up-regulated at 7.5 h post-exposure and peaked at 10 h post-exposure with a 2.51-fold increase (Figure 2D). In the gills, the expression levels were significantly elevated at 5 and 7.5 h with 1.36-fold and 1.94-fold increases, respectively (Figure 3D).
Figure 2. Expression profiles of several apoptosis and hypoxia-related genes in the hepatopancreas. The relative mRNA expression levels of IAP (A), TXNIP (B), Caspase (C), Caspase3 (D), Hsp70 (E), Hif-1α (F), HcY (G). Each bar represents the mean ± SD (n = 6). A significant difference between groups at p < 0.05 (ANOVA) is indicated by “*”.
Figure 3. Expression profiles of several apoptosis and hypoxia-related genes in the gill. The relative mRNA expression levels of IAP (A), TXNIP (B), Caspase (C), Caspase3 (D), Hsp70 (E), Hif-1α (F), HcY (G). Each bar represents the mean ± SD (n = 6). A significant difference between groups at p < 0.05 (ANOVA) is indicated by “*”.
After 5 h post-exposure, the hsp70 transcriptional levels in the hepatopancreas showed significant up-regulation that peaked at the 7.5 h time point, and was 3.14-fold higher when compared with the control group (p < 0.05), these levels then decreased significantly back to the baseline levels (Figure 2E). In the gills, the mRNA expression of hsp70 was significantly elevated at 2.5 and 5 h with 2.13- and 3.5-fold increases (p < 0.05), respectively (Figure 3E). The transcriptional levels of hif-1α in the hepatopancreas increased significantly by 1.87- and 2.83-fold (p < 0.05) at 10 and 12.5 h post- exposure, respectively (Figure 2F). In the gills, the hif-1α transcriptional levels showed a significant incremental change at the 10 h time point with a 1.95-fold increase (p < 0.05) (Figure 3F). The obvious up-regulation of the HcY expression in the hepatopancreas was detected at the 5 h time point with a 1.2-fold increase (Figure 2G). The HcY expression levels in the gills were significantly elevated at 2.5 and 5 h with 1.34- and 1.88-fold increases (p < 0.05), respectively (Figure 3G).
The Cumulative Mortality and Physiological Indexes
The cumulative mortality and physiological indexes showed time-dependent changes upon air exposure (Figure 4). Several individuals showed stress symptoms, and nine individuals died at the 5 h time point (the mortality rate was 4.5%) (Figure 4A). As the exposure time progressed, the cumulative mortality rate was 15% at 7.5 h post-exposure, and then increased markedly to 89.5% at 12.5 h. The heart rates exhibited an initial increase and a subsequent decrease (Figure 4A). The heart rates significantly up-regulated at the 2.5 h time point and then peaked at the 5 h time point at 181 beats/min. The heart rates then quickly decreased until faint beating was observed at the 12.5 h time point. T-AOC significantly increased at 2.5 h post-exposure, which was a 1.3-fold increase (p < 0.05) compared with the control group (Figure 4B). After the peak, T-AOC decreased rapidly to the lowest level at the 12.5 h time point (1.9 U/mg protein). The lactate concentrations remained elevated through the 2.5 and 5 h post-exposure points, which were increased 1.72- and 2.1-fold, respectively, compared with the control group and then decreased rapidly (Figure 4B). At the 12.5 h time point, the lactate concentrations were still higher than those of the control group.
Figure 4. Cumulative mortality rates and physiological indexes of M. japonicus under air exposure. (A) Cumulative mortality rate and heart rates. (B) T-AOC and lactate contents.
Detection of Apoptosis in the Hepatopancreas and Gills via TUNEL
In this study, air exposure significantly induced apoptosis in the hepatopancreas and gill cells (Figure 5). TUNEL-positive cell nuclei appeared green and common cells appeared blue. For each sample, four visual fields were randomly selected for counting the apoptosis index (AI). In the hepatopancreas, the apoptosis index (AI) values of the control and the 12.5 h post-exposure group were 10.5 and 30.5%, respectively. Compared with the control group, the AI for the 12.5 h experimental group increased significantly (p < 0.05). In the gills, the AI for 12.5 h post-exposure group was significantly higher than that for the control group (p < 0.05), at 24.5 and 9.5%, respectively. Compared with the gill tissues, the hepatopancreas showed a higher mean AI at 12.5 h post-exposure. However, no significant changes (p > 0.05) were observed between the two groups.
Figure 5. Apoptosis analysis of M. japonicus hepatopancreas and gill (200×). TUNEL-positive cell nuclei appear green, and parts are indicated by the arrows.
Changes in DNA Methylation Patterns in the Hepatopancreas
Methylation-sensitive amplification polymorphism fragment gel files were transformed into chromatogram files by GeneMapper (Figure 6A and Supplementary Table S1). All fragments were sieved and classified by size using a Python script (Figure 6B). The methylation ratios of ten individuals from the control group and from 12.5 h after exposure are presented in Figure 7A. Most individuals in the control group (CG) had higher methylation ratios than those in the experimental group (EG). A total of 4635 fragments for the CG and 4524 fragments for the EG were obtained (Figure 7B). The numbers of fragments of types I, II, and III were 3180, 586, and 869 for CG and were 3327, 545, and 652 for EG, respectively. The computational results showed a significant difference in the full-methylation levels (p < 0.05). However, there was no significant difference between the hemi-methylation levels. For both CG and EG, the full-methylation ratios were higher than the hemi-methylation ratio (Figure 7B). The average methylation ratio was 31.39% for CG and 26.46% for EG.
Figure 6. Methylation patterns detected in M. japonicus hepatopancreas tissues. (A) A chromatogram file from the GeneMapper software. The height peak indicates the quantity of amplification and the horizontal scale represents the molecular weight of the polymorphic fragments. (B) Polymorphic fragments were sieved and classified by size using a Python script. T1 = EH+/EM+, T2 = EH+/EM–, and T3 = EH–/EM+.
Figure 7. Methylation levels of individuals and groups. (A) The methylation ratio of ten individuals from control and 12.5 h post-exposure. (B) The methylation level of M. japonicus hepatopancreas under air exposure. I = EH+/EM+, II = EH+/EM–, III = EH–/EM+, T_Met = II + III, T_bands = I + II + III. H_Met = Half methylation, F_Met = Full methylation, T_Met = Total methylation. Each bar represents the mean ± SD (n = 10). A significant difference between groups at p < 0.05 (ANOVA) is indicated by different letters.
As a swimming and burrowing species, M. japonicus is usually subjected to aerial exposure during the processes of packing and transportation (Hewitt and Duncan, 2001; Abe et al., 2007; Fotedar and Evans, 2011). In our study, some individuals showed stress symptoms and nine individuals died at 5 h post- exposure. The cumulative mortality rate rapidly increased after 7.5 h post-exposure. At the 12.5 h time point, the remaining individuals were only sufficient for sampling. The aerial stress tolerance was consistent with the results of Duan et al. (2016a). Heart rate, a direct measure of the physiological state, is correlated with oxygen consumption and the metabolic rate (Green, 2011; Claireaux and Chabot, 2016; Lund et al., 2017). In the early stages, the heart rate significantly increased for oxygen supply and ATP production. Because the almost transparent heart is covered with a carapace, there are still no reliable tools to record heart rates. The lower heart rates were a sign of a decline in overall bodily function. During the experiment, we monitored the gill ventilation frequency, which reflected individual vitality levels. The gill ventilation frequencies exhibited an initial increase and a subsequent decrease, and could not be observed directly at the final stage.
During air exposure, T-AOC showed time-dependent changes. In the early stages of exposure, the antioxidant defense system increased T-AOC. As the exposure time progressed, the increased activity of the antioxidant enzymes was not sufficient to effectively remove ROS, however, the body needs to constantly consume antioxidants to reduce oxidative damage. Liu H.L. et al. (2015) showed that T-AOC in the hepatopancreas of L. vannamei decreased after 10 min of air exposure and then increased. During hypoxia, crustaceans utilize anaerobic glycolysis to satisfy their energy requirements, thus accumulating lactate (Soñanez-Organis et al., 2010; Dunbar et al., 2017; Valère-Rivet et al., 2019). For the early stages of air exposure, the lactate concentrations remained elevated through 2.5 and 5 h post-exposure, which were consistent with the results of Abe et al. (Abe et al., 2007; Aparicio-Simón et al., 2018). Subsequently, the comprehensive metabolic system of the shrimp was severely affected, and part of the lactic acid was converted into glucose to maintain homeostasis and the energy supply, thus reducing the lactic acid content.
The duration of air exposure significantly influenced the expression profiles of the hypoxia-related genes in the hepatopancreas and gills. Hepatopancreas, a multifunctional organ, integrate metabolism and immune functions. Gills, a crustacean respiratory organ, participate in ion transport, acid-base balance, ammonia excretion, and heavy metal accumulation (Henry et al., 2012). Overall, gill tissues may be more sensitive to desiccation stress than the hepatopancreas tissues, which makes sense given their direct exposure to the environment. In this study, the Hif-1α expression levels in the hepatopancreas and gill tissues significantly increased at the 10 and 7.5 h time points, which is in accordance with most research results (Sun et al., 2016b; Okamura et al., 2018). In the hepatopancreas, the transcriptional level of Hif-1α peaked at the 12.5 h time point, while the expression levels in the gills exhibited an initial increase and a subsequent decrease. As an important organ of respiration, the microvascular cavity at the edge of the gill filaments lobule is filled with a large amount of hemocyanin for oxygen transport. The HcY expression increased gradually within 2.5 h of air exposure in the gills, but was not obvious in the hepatopancreas. With a change in dissolved oxygen, crustaceans increase the hemolymph pH or change the composition of the Hc molecules to enhance the hemocyanin oxygen affinity (Coates and Nairn, 2014). At the initial stage of air exposure, the up-regulation of lactate concentrations might affect pH values, which weakened the binding capacity between hemocyanin and oxygen (Whiteley et al., 1997). Further studies are needed to investigate hemocyanin concentrations during air exposure.
Continuous oxidative stress triggers caspase-independent cell death and induces apoptosis or necrosis (Simon et al., 2000). Several stress conditions, such as high temperatures, hypoxia exposure, and ultraviolet (UV) irradiation, have induced extensive apoptosis in the larval embryos of the Japanese flounder (Yabu et al., 2003; Williams et al., 2017). In this study, air exposure significantly affected the percentage of TUNEL-positive cells. Compared with the control group, the apoptosis index (AI) for the 12.5 h experimental group increased significantly (p < 0.05) in the hepatopancreas and gills. A higher apoptosis index was detected in the hepatopancreas than in the gills, which was possibly because to the hepatopancreas is the major metabolic center of ROSs production in crustaceans (Arun and Subramanian, 1998; Moore et al., 2007). The MjIAP showed different expression patterns in the hepatopancreas and gills. An earlier upward trend was presented in gills at 2.5 h, which might suppress the early apoptosis of the gill cells. Further studies are needed to investigate the histological symptoms and ultrastructure of the hepatopancreas and gills.
As the major factors for apoptosis, the caspases are a family of structurally related cysteine proteases and contain inflammatory mediator (1,4,5,11), activator (2,8,9,10) and executioner (3,6,7) subfamilies (Shalini et al., 2015; Julien and Wells, 2017). Chronic exposure of croaker to hypoxia increased caspase-3/7 activity, as well as the apoptosis of ovarian follicle cells (Ondricek and Thomas, 2018). Hypoxia ischemia can result in a significant increase in the caspase-3 expression levels and neuronal apoptosis in the brains of neonatal mice (Deng et al., 2019). In this study, the caspase and caspase-3 expression levels increased initially and then decreased, and the caspase gene was more sensitive than the caspase-3 gene. The caspase gene showed sequence similarity to the caspase-8 gene of other species (Wang et al., 2008). The caspase protein was essential for virus-induced apoptosis in M. japonicus. As an activator, caspase-8 activation can induce direct cleavage of downstream caspase-3/7 (Kaufmann et al., 2012; Feltham et al., 2017). Caspase-3/8 activities will be detected in further studies. As an endogenous inhibitor of caspases, IAP does not bind to caspase-8 but inhibits its substrate, caspase-3, to execute anti-apoptotic roles (Wei et al., 2008). The MjIAP expression levels in the hepatopancreas gradually increased, reaching an increase of 9.88-fold (p < 0.05) at 12.5 h post- exposure, which may induce a decrease in caspase and caspase-3 expression in the final stage. TXNIP, an inhibitor of the major ROS scavenger TRX, was positively associated with antibacterial responses but was negatively correlated with antiviral responses in shrimp (Zuo et al., 2019). In the hepatopancreas and gills, the transcriptional level of MjTXNIP exhibited similar variation trends.
DNA methylation is an important epigenetic modification, which may contribute to environmentally induced phenotypic variations by regulating gene transcription (Angers et al., 2010; Roberts and Gavery, 2012). Zhang et al. (2017) showed that the total methylation levels of Crassostrea gigas gills and adductors significantly increased at 12 h post-exposure. The methylation levels were different for different tissues of the same individual or for different states of the same tissue, even in the offspring of full-sibling families. The high methylation diversity among individuals may be related to the high genome heterozygosity of shrimp (Yuan et al., 2018; Zhang et al., 2019). The control group and the experimental group all had one abnormal individual, but this did not affect the overall outcome. The air exposure led to a lower methylation ratio in the experimental group than in the control group. There are varying degrees of methylation in invertebrates, which are often below the rate for vertebrates (Gavery and Roberts, 2010; Roberts and Gavery, 2012). The methylation rates in this study were higher than those of Fenneropenaeus chinensis (He et al., 2015). DNA methylation has an obvious tissue-specific pattern (Sun et al., 2014; Zhao et al., 2015). Except for the hepatopancreas, the DNA methylation changes in other tissues still need further verification. It is well known that DNA methylation in promoter regions normally represses gene expression (Suzuki and Bird, 2008; Maunakea et al., 2010; Jones, 2012; Zhou et al., 2019). Wan et al. (2015) determined the existence of two different modes of DNA methylation-dependent gene regulation. The DNA methylation of IAP, TXNIP, caspase, and caspase-3 genes may play a key role in the differential expression patterns, which needs further verification.
This study provided a relatively comprehensive understanding of the desiccation and oxidative stress responses of M. japonicus in terms of physiology, innate immunity, cell apoptosis, and DNA methylation levels, and facilitated further investigation into the molecular mechanisms of shrimp under oxidative stress. In the coming centuries, aquatic animals will face more severe challenges, including ocean warming and the expansion of the death zone (Keeling et al., 2010; Altieri and Gedan, 2015; Vinagre et al., 2019). This research can serve as a case study for future research on hypoxia adaptation and the antioxidant mechanisms of other marine groups.
Data Availability Statement
The data of MjTXNIP (MN265397) and MjIAP (MN265396) are publicly accessible in the NCBI (https://www.ncbi.nlm.nih.gov/).
All experimental procedures were conducted in conformity with institutional guidelines for the care and use of laboratory animals in Xiamen University, Xiamen, China.
PW, YM, JW, and YS designed the experiments. PW and ZL conducted computational and statistical analyses. PW and YM drafted the manuscript. All authors read and approved the final manuscript.
This study was supported by the China’s Agricultural Research System (Grant No. CARS-48); the Fujian Provincial Department of Science and Technology Project (Grant No. 2016NZ0001-4); the Marine Economy Innovation and Development Project of Xiamen (Grant No. 16CZY009SF05); Special funds for Ocean and Fishery Structure Adjustment in Fujian Province in 2019 (Fujian Financial Index No. 140); and Science and Technology Plan Project of Ningbo (2019B10011).
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.
The authors are grateful to the Dongshan Swire Marine Station of Xiamen University (D-SMART) for providing the laboratory space.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.00223/full#supplementary-material
- ^ https://www.ncbi.nlm.nih.gov/orffinder/
- ^ https://www.ebi.ac.uk/Tools/st/emboss_transeq/
- ^ https://www.expasy.org/
- ^ http://smart.embl.de/
Abasubong, K. P., Liu, W.-B., Zhang, D.-D., Yuan, X.-Y., Xia, S.-L., Xu, C., et al. (2018). Fishmeal replacement by rice protein concentrate with xylooligosaccharides supplement benefits the growth performance, antioxidant capability and immune responses against Aeromonas hydrophila in blunt snout bream (Megalobrama amblycephala). Fish Shellfish Immunol. 78, 177–186. doi: 10.1016/j.fsi.2018.04.044
Abe, H., Hirai, S., and Okada, S. (2007). Metabolic responses and arginine kinase expression under hypoxic stress of the kuruma prawn Marsupenaeus japonicus. Compar. Biochem. Physiol. Part A 146, 40–46. doi: 10.1016/j.cbpa.2006.08.027
Andrade, M., De Marchi, L., Soares, A. M., Rocha, R. J., Figueira, E., and Freitas, R. (2019). Are the effects induced by increased temperature enhanced in Mytilus galloprovincialis submitted to air exposure? Sci. Total Environ. 647, 431–440. doi: 10.1016/j.scitotenv.2018.07.293
Angers, B., Castonguay, E., and Massicotte, R. (2010). Environmentally induced phenotypes and DNA methylation: how to deal with unpredictable conditions until the next generation and after. Mol. Ecol. 19, 1283–1295. doi: 10.1111/j.1365-294X.2010.04580.x
Aparicio-Simón, B., Piñón, M., Racotta, R., and Racotta, I. S. (2018). Neuroendocrine and metabolic responses of Pacific whiteleg shrimp Penaeus vannamei exposed to hypoxia stress. Latin Am. J. Aquat. Res. 46, 364–376. doi: 10.3856/vol46-issue2-fulltext-12
Arun, S., and Subramanian, P. (1998). Antioxidant enzymes in freshwater prawn Macrobrachium malcolmsonii during embryonic and larval development. Compar. Biochem. Physiol. Part B Biochem. Mol. Biol. 121, 273–277. doi: 10.1016/s0305-0491(98)10100-1
Cai, X., Zhang, D., Wang, J., Liu, X., Ouyang, G., and Xiao, W. (2018). Deletion of the fih gene encoding an inhibitor of hypoxia-inducible factors increases hypoxia tolerance in zebrafish. J. Biol. Chem. 293, 15370–15380. doi: 10.1074/jbc.RA118.003004
Camacho-Jiménez, L., Leyva-Carrillo, L., Peregrino-Uriarte, A. B., Duarte-Gutiérrez, J. L., Tresguerres, M., and Yepiz-Plascencia, G. (2019). Regulation of glyceraldehyde-3-phosphate dehydrogenase by hypoxia inducible factor 1 in the white shrimp Litopenaeus vannamei during hypoxia and reoxygenation. Compar. Biochem. Physiol. Part A 235, 56–65. doi: 10.1016/j.cbpa.2019.05.006
Cui, Q., Chen, F.-Y., Chen, H.-Y., Peng, H., and Wang, K.-J. (2019). Benzo [a] pyrene (BaP) exposure generates persistent reactive oxygen species (ROS) to inhibit the NF-κB pathway in medaka (Oryzias melastigma). Environ. Pollut. 251, 502–509. doi: 10.1016/j.envpol.2019.04.063
Daugaard, M., Rohde, M., and Jäättelä, M. (2007). The heat shock protein 70 family: highly homologous proteins with overlapping and distinct functions. FEBS Lett. 581, 3702–3710. doi: 10.1016/j.febslet.2007.05.039
Deng, C., Li, J., Li, L., Sun, F., and Xie, J. (2019). Effects of hypoxia ischemia on caspase-3 expression and neuronal apoptosis in the brain of neonatal mice. Exp. Therap. Med. 17, 4517–4521. doi: 10.3892/etm.2019.7487
Dowen, R. H., Pelizzola, M., Schmitz, R. J., Lister, R., Dowen, J. M., Nery, J. R., et al. (2012). Widespread dynamic DNA methylation in response to biotic stress. Proc. Natl. Acad. Sci. U.S.A. 109, E2183–E2191. doi: 10.1073/pnas.1209329109
Duan, Y., Zhang, J., Dong, H., Wang, Y., Liu, Q., and Li, H. (2016a). Effect of desiccation and resubmersion on the oxidative stress response of the kuruma shrimp Marsupenaeus japonicus. Fish Shellfish Immunol. 49, 91–99. doi: 10.1016/j.fsi.2015.12.018
Duan, Y., Zhang, Y., Dong, H., and Zhang, J. (2016b). Effect of desiccation on oxidative stress and antioxidant response of the black tiger shrimp Penaeus monodon. Fish Shellfish Immunol. 58, 10–17. doi: 10.1016/j.fsi.2016.09.004
Dunbar, S. G., Shives, J., and Boskovic, D. S. (2017). Lactate accumulation in the intertidal hermit crab, Pagurus samuelis, in response to burial-induced hypoxia. Crustacean Res. 46, 121–132. doi: 10.18353/crustacea.46.0_121
Garcia-Orozco, K. D., Sanchez-Paz, A., Aispuro-Hernandez, E., Gomez-Jimenez, S., Lopez-Zavala, A., Araujo-Bernal, S., et al. (2012). Gene expression and protein levels of thioredoxin in the gills from the whiteleg shrimp (Litopenaeus vannamei) infected with two different viruses: the WSSV or IHHNV. Fish Shellfish Immunol. 32, 1141–1147. doi: 10.1016/j.fsi.2012.03.020
Gavery, M. R., and Roberts, S. B. (2010). DNA methylation patterns provide insight into epigenetic regulation in the Pacific oyster (Crassostrea gigas). BMC Geno. 11:483. doi: 10.1186/1471-2164-11-483
He, J., Yu, Y., Qin, X. W., Zeng, R. Y., Wang, Y. Y., Li, Z. M., et al. (2019). Identification and functional analysis of the Mandarin fish (Siniperca chuatsi) hypoxiainducible factor-1α involved in the immune response. Fish Shellfish Immunol. 92, 141–150. doi: 10.1016/j.fsi.2019.04.298
He, Y., Du, Y., Li, J., Liu, P., Wang, Q., and Li, Z. (2015). Analysis of DNA methylation in different tissues of Fenneropenaeus chinensis from the wild population and Huanghai No. 1. Acta Oceanol. Sin. 34, 175–180. doi: 10.1007/s13131-015-0765-x
Henry, R. P., Lucu, C., Onken, H., and Weihrauch, D. (2012). Multiple functions of the crustacean gill: osmotic/ionic regulation, acid-base balance, ammonia excretion, and bioaccumulation of toxic metals. Front. Physiol. 3:431. doi: 10.3389/fphys.2012.00431
Hewitt, D., and Duncan, P. (2001). Effect of high water temperature on the survival, moulting and food consumption of Penaeus (Marsupenaeus) japonicus (Bate, 1888). Aquacult. Res. 32, 305–313. doi: 10.1046/j.1365-2109.2001.00560.x
Holze, C., Michaudel, C., Mackowiak, C., Haas, D. A., Benda, C., Hubel, P., et al. (2018). Oxeiptosis, a ROS-induced caspase-independent apoptosis-like cell-death pathway. Nat. Immunol. 19, 130–140. doi: 10.1038/s41590-017-0013-y
Huo, D., Sun, L., Ru, X., Zhang, L., Lin, C., Liu, S., et al. (2018). Impact of hypoxia stress on the physiological responses of sea cucumber Apostichopus japonicus: respiration, digestion, immunity and oxidative damage. PeerJ 6:e4651. doi: 10.7717/peerj.4651
Ighodaro, O., and Akinloye, O. (2018). First line defence antioxidants-superoxide dismutase (SOD), catalase (CAT) and glutathione peroxidase (GPX): their fundamental role in the entire antioxidant defence grid. Alexand. J. Med. 54, 287–293. doi: 10.1016/j.ajme.2017.09.001
Jia, F.-X., Dou, W., Hu, F., and Wang, J.-J. (2011). Effects of thermal stress on lipid peroxidation and antioxidant enzyme activities of oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae). Florida Entomol. 94, 956–964.
Lennicke, C., Rahn, J., Lichtenfels, R., Wessjohann, L. A., and Seliger, B. (2015). Hydrogen peroxide–production, fate and role in redox signaling of tumor cells. Cell Commun. Signal. 13:39. doi: 10.1186/s12964-015-0118-6
Liu, H.-L., Yang, S.-P., Guo, W.-J., Tan, Z.-H., and Chan, S. F. (2017). Effects of air exposure and re-submersion on oxidative stress of marine gastropod, Babylonia areolata. Israeli J. Aquacult. Bamidgeh 69, 1–12.
Liu, H. L., Yang, S. P., Wang, C. G., Chan, S. M., Wang, W. X., Feng, Z. H., et al. (2015). Effect of air exposure and resubmersion on the behavior and oxidative stress of Pacific white shrimp Litopenaeus vannamei. North Am. J. Aquacult. 77, 43–49. doi: 10.1080/15222055.2014.955157
Liu, P.-F., Liu, Q.-H., Wu, Y., and Huang, J. (2015). Thioredoxin of Litopenaeus vannamei facilitated white spot syndrome virus infection. J. Invert. Pathol. 129, 57–62. doi: 10.1016/j.jip.2015.05.009
Lorenzon, S., Giulianini, P. G., Libralato, S., Martinis, M., and Ferrero, E. (2008). Stress effect of two different transport systems on the physiological profiles of the crab Cancer pagurus. Aquaculture 278, 156–163. doi: 10.1016/j.aquaculture.2008.03.011
Lund, M., Dahle, M. K., Timmerhaus, G., Alarcon, M., Powell, M., Aspehaug, V., et al. (2017). Hypoxia tolerance and responses to hypoxic stress during heart and skeletal muscle inflammation in Atlantic salmon (Salmo salar). PLoS One 12:e0181109. doi: 10.1371/journal.pone.0181109
Maunakea, A. K., Nagarajan, R. P., Bilenky, M., Ballinger, T. J., D’souza, C., Fouse, S. D., et al. (2010). Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature 466:253. doi: 10.1038/nature09165
Moore, M. N., Viarengo, A., Donkin, P., and Hawkins, A. J. (2007). Autophagic and lysosomal reactions to stress in the hepatopancreas of blue mussels. Aquat. Toxicol. 84, 80–91. doi: 10.1016/j.aquatox.2007.06.007
Murgatroyd, C., Patchev, A. V., Wu, Y., Micale, V., Bockmühl, Y., Fischer, D., et al. (2009). Dynamic DNA methylation programs persistent adverse effects of early-life stress. Nat. Neurosci. 12:1559. doi: 10.1038/nn.2436
Okamura, Y., Mekata, T., Elshopakey, G. E., and Itami, T. (2018). Molecular characterization and gene expression analysis of hypoxia-inducible factor and its inhibitory factors in kuruma shrimp Marsupenaeus japonicus. Fish Shellfish Immunol. 79, 168–174. doi: 10.1016/j.fsi.2018.05.015
Ondricek, K., and Thomas, P. (2018). Effects of hypoxia exposure on apoptosis and expression of membrane steroid receptors, ZIP9, mPRα, and GPER in Atlantic croaker ovaries. Compar. Biochem. Physiol. Part A 224, 84–92. doi: 10.1016/j.cbpa.2018.07.002
Paital, B. (2013). Antioxidant and oxidative stress parameters in brain of Heteropneustes fossilis under air exposure condition; role of mitochondrial electron transport chain. Ecotoxico. Environ. Saf. 95, 69–77. doi: 10.1016/j.ecoenv.2013.05.016
Paital, B., and Chainy, G. (2010). Antioxidant defenses and oxidative stress parameters in tissues of mud crab (Scylla serrata) with reference to changing salinity. Compar. Biochem. Physiol. Part C 151, 142–151. doi: 10.1016/j.cbpc.2009.09.007
Romero, M. C., Ansaldo, M., and Lovrich, G. A. (2007). Effect of aerial exposure on the antioxidant status in the subantarctic stone crab Paralomis granulosa (Decapoda: Anomura). Compar. Biochem. Physiol. Part C 146, 54–59. doi: 10.1016/j.cbpc.2006.06.009
Romero, M. C., Tapella, F., Sotelano, M. P., Ansaldo, M., and Lovrich, G. A. (2011). Oxidative stress in the subantarctic false king crab Paralomis granulosa during air exposure and subsequent re-submersion. Aquaculture 319, 205–210. doi: 10.1016/j.aquaculture.2011.06.041
Sahu, P. P., Pandey, G., Sharma, N., Puranik, S., Muthamilarasan, M., and Prasad, M. (2013). Epigenetic mechanisms of plant stress responses and adaptation. Plant Cell Rep. 32, 1151–1159. doi: 10.1007/s00299-013-1462-x
Schvezov, N., Lovrich, G. A., Tapella, F., Gowland-Sainz, M., and Romero, M. C. (2019). Effect of the temperature of air exposure on the oxidative stress status of commercial male southern king crab Lithodes santolla. Fish. Res. 212, 188–195. doi: 10.1016/j.fishres.2018.12.020
Soñanez-Organis, J. G., Racotta, I. S., and Yepiz-Plascencia, G. (2010). Silencing of the hypoxia inducible factor 1–HIF-1-obliterates the effects of hypoxia on glucose and lactate concentrations in a tissue-specific manner in the shrimp Litopenaeus vannamei. J. Exp. Mar. Biol. Ecol. 393, 51–58. doi: 10.1016/j.jembe.2010.06.031
Sun, S., Gu, Z., Fu, H., Zhu, J., Ge, X., and Xuan, F. (2016a). Molecular cloning, characterization, and expression analysis of p53 from the oriental river prawn, Macrobrachium nipponense, in response to hypoxia. Fish Shellfish Immunol. 54, 68–76. doi: 10.1016/j.fsi.2016.03.167
Sun, S., Xuan, F., Fu, H., Ge, X., Zhu, J., Qiao, H., et al. (2016b). Molecular characterization and mRNA expression of hypoxia inducible factor-1 and cognate inhibiting factor in Macrobrachium nipponense in response to hypoxia. Compar. Biochem. Physiol. Part B 196, 48–56. doi: 10.1016/j.cbpb.2016.02.002
Sun, Y., Hou, R., Fu, X., Sun, C., Wang, S., Wang, C., et al. (2014). Genome-wide analysis of DNA methylation in five tissues of Zhikong scallop, Chlamys farreri. PLoS One 9:e86232. doi: 10.1371/journal.pone.0086232
Tsoi, K. H., Ma, K. Y., Wu, T., Fennessy, S. T., Chu, K. H., and Chan, T. Y. (2014). Verification of the cryptic species Penaeus pulchricaudatus in the commercially important kuruma shrimp P. japonicus (Decapoda: Penaeidae) using molecular taxonomy. Invert. Syst. 28, 476–490.
Valère-Rivet, M. G., Boskovic, D. S., Estevez, D., and Dunbar, S. G. (2019). Changes in hemolymph lactate and ammonia in the hermit crab Pagurus samuelis (Stimpson, 1857)(Decapoda: Anomura: Paguridae) during shallow burial. J. Crustacean Biol. 39, 172–180. doi: 10.1093/jcbiol/ruy114
Vinagre, C., Dias, M., Cereja, R., Abreu-Afonso, F., Flores, A. A., and Mendonça, V. (2019). Upper thermal limits and warming safety margins of coastal marine species–Indicator baseline for future reference. Ecol. Indic. 102, 644–649. doi: 10.1016/j.ecolind.2019.03.030
Wan, J., Oliver, V. F., Wang, G., Zhu, H., Zack, D. J., Merbs, S. L., et al. (2015). Characterization of tissue-specific differential DNA methylation suggests distinct modes of positive and negative gene expression regulation. BMC Genom. 16:49. doi: 10.1186/s12864-015-1271-4
Wang, S., Lv, J., Zhang, L., Dou, J., Sun, Y., Li, X., et al. (2015). MethylRAD: a simple and scalable method for genome-wide DNA methylation profiling using methylation-dependent restriction enzymes. Open Biol. 5:150130. doi: 10.1098/rsob.150130
Wang, S., Zhang, Q., Zheng, S., Chen, M., Zhao, F., and Xu, S. (2019). Atrazine exposure triggers common carp neutrophil apoptosis via the CYP450s/ROS pathway. Fish Shellfish Immunol. 84, 551–557. doi: 10.1016/j.fsi.2018.10.029
Wang, Z., Jiang, H., Chen, S., Du, F., and Wang, X. (2012). The mitochondrial phosphatase PGAM5 functions at the convergence point of multiple necrotic death pathways. Cell 148, 228–243. doi: 10.1016/j.cell.2011.11.030
Williams, T. A., Bergstrome, J. C., Scott, J., and Bernier, N. J. (2017). CRF and urocortin 3 protect the heart from hypoxia/reoxygenation-induced apoptosis in zebrafish. Am. J. Physiol. Regul. Integr. Compar. Physiol. 313, R91–R100. doi: 10.1152/ajpregu.00045.2017
Xiong, L., Xu, C., Maroof, M. S., and Zhang, Q. (1999). Patterns of cytosine methylation in an elite rice hybrid and its parental lines, detected by a methylation-sensitive amplification polymorphism technique. Mol. Gen. Genet. MGG 261, 439–446. doi: 10.1007/s004380050986
Xu, Z., Guan, W., Xie, D., Lu, W., Ren, X., Yuan, J., et al. (2019). Evaluation of immunological response in shrimp Penaeus vannamei submitted to low temperature and air exposure. Dev. Compar. Immunol. 100:103413. doi: 10.1016/j.dci.2019.103413
Xu, Z., Regenstein, J. M., Xie, D., Lu, W., Ren, X., Yuan, J., et al. (2018). The oxidative stress and antioxidant responses of Litopenaeus vannamei to low temperature and air exposure. Fish Shellfish Immunol. 72, 564–571. doi: 10.1016/j.fsi.2017.11.016
Yuan, J., Zhang, X., Liu, C., Yu, Y., Wei, J., Li, F., et al. (2018). Genomic resources and comparative analyses of two economical penaeid shrimp species, Marsupenaeus japonicus and Penaeus monodon. Mar. Genom. 39, 22–25. doi: 10.1016/j.margen.2017.12.006
Zhang, X., Yuan, J., Sun, Y., Li, S., Gao, Y., Yu, Y., et al. (2019). Penaeid shrimp genome provides insights into benthic adaptation and frequent molting. Nat. Commun. 10:356. doi: 10.1038/s41467-018-08197-4
Zhao, Y., Chen, M., Storey, K. B., Sun, L., and Yang, H. (2015). DNA methylation levels analysis in four tissues of sea cucumber Apostichopus japonicus based on fluorescence-labeled methylation-sensitive amplified polymorphism (F-MSAP) during aestivation. Compar. Biochem. Physiol. Part B 181, 26–32. doi: 10.1016/j.cbpb.2014.11.001
Zhou, R., De Koning, D. J., Mccormack, H., Wilson, P., and Dunn, I. (2019). Short tandem repeats and methylation in the promoter region affect expression of cystathionine beta-synthase gene in the laying hen. Gene 710, 367–374. doi: 10.1016/j.gene.2019.05.049
Keywords: Marsupenaeus japonicus, air exposure, oxidative stress, non-specific immunity, apoptosis, DNA methylation
Citation: Wang P, Wang J, Su Y, Liu Z and Mao Y (2020) Air Exposure Affects Physiological Responses, Innate Immunity, Apoptosis and DNA Methylation of Kuruma Shrimp, Marsupenaeus japonicus. Front. Physiol. 11:223. doi: 10.3389/fphys.2020.00223
Received: 11 September 2019; Accepted: 26 February 2020;
Published: 12 March 2020.
Edited by:Ignacio Ruiz-Jarabo, University of Cádiz, Spain
Reviewed by:Maria Giulia Lionetto, University of Salento, Italy
Brad Buckley, Portland State University, United States
Copyright © 2020 Wang, Wang, Su, Liu and Mao. 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: Yong Mao, email@example.com