Skip to main content

Comparative proteomic analysis of Ulva prolifera response to high temperature stress



Ulva prolifera belongs to green macroalgae and is the dominant species of green tide. It is distributed worldwide and is therefore subject to high-temperature stress during the growth process. However, the adaptation mechanisms of the response of U. prolifera to high temperatures have not been clearly investigated yet.


In this study, isobaric tags for relative and absolute quantitation (iTRAQ) labelling was applied in combination with the liquid chromatography-tandem mass spectrometry (LC-MS/MS) to conduct comparative proteomic analysis of the response of U. prolifera to high-temperature stress and to elucidate the involvement of this response in adaptation mechanisms. Differentially expressed proteins (DEPs) of U. prolifera under high temperature (denote UpHT) compared with the control (UpC) were identified. Bioinformatic analyses including GO analysis, pathway analysis, and pathway enrichment analysis was performed to analyse the key metabolic pathways that underlie the thermal tolerance mechanism through protein networks. Quantitative real-time PCR and western blot were performed to validate selected proteins.


In the present study, 1223 DEPs were identified under high temperature compared with the control, which included 790 up-regulated and 433 down-regulated proteins. The high-temperature stimulus mainly induced the expression of glutathione S-transferase, heat shock protein, ascorbate peroxidase, manganese superoxide dismutase, ubiquitin-related protein, lhcSR, rubisco activase, serine/threonine protein kinase 2, adenylate kinase, Ca2+-dependent protein kinase (CDPK), disease resistance protein EDS1, metacaspase type II, NDPK2a, 26S proteasome regulatory subunit, ubiquinone oxidoreductase, ATP synthase subunit, SnRK2s, and cytochrome P450. The down-regulated proteins were photosynthesis-related proteins, glutathione reductase, catalase-peroxidase, thioredoxin, thioredoxin peroxidase, PP2C, and carbon fixation-related proteins. Furthermore, biological index analysis indicated that protein content and SOD activity decreased; the value of Fv/Fm dropped to the lowest point after culture for 96 h. However, APX activity and MDA content increased under high temperature.


The present study implied an increase in proteins that were associated with the stress response, oxidative phosphorylation, the cytokinin signal transduction pathway, the abscisic acid signal transduction pathway, and the glutathione metabolism pathway. Proteins that were associated with photosynthesis, carbon fixation in photosynthesis organisms, and the photosynthesis antenna protein pathway were decreased. These pathways played a pivotal role in high temperature regulation. These novel proteins provide a good starting point for further research into their functions using genetic or other approaches. These findings significantly improve the understanding of the molecular mechanisms involved in the tolerance of algae to high-temperature stress.


Ulva prolifera (Chlorophyta) belongs to Chlorophyta, Ulvales, Ulvaceae and is a type of adaptable large algae. It is rich in protein, carbohydrates, crude fibre, vitamins, amino acids, unsaturated fatty acids and other bioactive substances. U. prolifera was developed as a food material, for refining raw materials, and as aquaculture feed. The natural reproductive ability of U. prolifera is strong, and it is the dominant species in the Yellow Sea and East China Sea. The breakout of U. prolifera in Qingdao since June 2008 has destroyed the natural environment and has a great influence on tourism and the aquatic industry. In addition, U. prolifera had been abundant in Xiangshan of the East China Sea from January to March. However, it disappeared in April. Therefore, change in the environment, especially temperature, had a great influence on the growth of U. prolifera. To adapt to these environmental changes, U. prolifera has developed mechanisms to adapt to different types of stresses including high temperatures, cold, hypersalinity and ultraviolet radiation [1, 2]. High temperature is a major environmental factor that limits U. prolifera yield. It can affect photosynthesis, respiration, water balance, membrane stability, hormone levels, and primary and secondary metabolites [3]. Direct injuries that are associated with high temperature include protein denaturation, aggregation, and increased fluidity of membrane lipids. Indirect or slower heat injuries include inactivation of enzymes in chloroplasts and mitochondria, inhibition of protein synthesis, protein degradation and loss of membrane integrity [4,5,6]. Photosynthesis is a process that is very sensitive to heat stress. The inhibition of photosystem II (PSII) leads to an invariable change in chlorophyll fluorescence [7]. To counter the effects of heat stress on cellular metabolism, plants respond to temperature changes by reprogramming their transcriptome, proteome, metabolome and lipidome. Such changes are aimed at establishing a new steady-state balance of metabolic processes that can enable the organism to function, survive and even reproduce at a higher temperature [5]. Previous studies have focused on physiological changes including photosynthesis, respiration, cell-membrane stability, hormone changes and induced antioxidant systems and heat-shock-protein expression in U. prolifera at high temperatures [8,9,10]. The effect of heat stress on U. prolifera appears to be associated with multiple processes and mechanisms including stress-related genes, transcription factors, and metabolism [11]. However, the transcription patterns are not always directly consistent with protein expression levels. To the best of our knowledge, the transcriptome expression analysis of the responses of U. prolifera to abiotic stresses such as low light or low temperature have been reported [1]. However, little attention has been given to the proteomic changed analysis of the response of U. prolifera to high temperatures.

Recent advances in the use of iTRAQ labelling of peptides and proteins for the relative quantification of an entire organism or cell proteome under different experimental conditions offer promising tools for the discovery of biomarkers at the protein level. Proteomic responses to heat stress have been widely studied in many species, such as Aspergillus flavus [12], Phaeodactylum tricornutum [13], and Pyropia yezoensis [14]. In addition, proteomic analysis also has been used to identify the metabolic responses of Synechocystis PCC 6803 to biofuel butanol, ethanol and hexane treatments [15,16,17]. However, it has been rarely reported in U. prolifera.

We describe here a procedure for proteomic analysis of the U. prolifera using iTRAQ labelling to identify differentially expressed proteins upon exposure to high temperature. The goal of the study is to use the molecular information in U. prolifera and provide a better understanding of the thermotolerance mechanism in U. prolifera and other macroalgae of the same habitats.


Macroalgae collection and preparation

Macroalgae used for this study were collected from the intertidal zone of the Xiangshan port of the East China Sea, Ningbo, China (121.82424 E, 29.552086 N) during low tide. U. prolifera samples were extensively rinsed with filtered (0.2 μm) seawater and placed into an enamel tray to further remove debris and epiphytes. Furthermore, the U. prolifera samples were disinfected with 0.2% KI for 10 min and flushed with sterile seawater. Next, U. prolifera samples that were inoculated in Provasoli medium [18] with 200 μg/mL ampicillin were placed in a biochemical incubator that was set to the appropriate temperature (25 °C) and light intensity (100 μmol s− 1 m− 2 on a 12 h: 12 h light dark cycle). To obtain relatively sterile materials, the U. prolifera samples were disinfected and transferred to fresh Provasoli medium every 7 days and continuously cultured for 10 days at 25 °C.

Experimental design

U. prolifera samples were transferred to a high-temperature 35 °C growth chamber and subjected to the same irradiance and photoperiod (denoted UpHT) as the treatment group, with three biological repeats. In addition, the U. prolifera samples were cultured at 25 °C as the control group (denoted UpC), with three biological repeats. The gametophyte thalli were harvested at 3 h, snap frozen using liquid nitrogen and maintained at − 70 °C until protein extraction. In addition, U. prolifera were collected at 3 h, 12 h, 24 h, 48 h and 96 h for the determination of physiological indexes.

Protein extraction

Proteins were extracted using the following method: the U. prolifera samples were disrupted in lysis buffer (including 0.1 M/L Tris-HCl, 1.4 M/L NaCl, 0.02 M/L Na2EDTA, 2% CTAB, 0.1% DIECA, 2% PVP K-30, 0.2% β-mercaptoethanol with pH adjusted to 8.0) with enzyme inhibitors (phenylmethyl sulphonyl fluoride, PMSF) using a tissue grinder and were sonicated on ice. The expected proteins were extracted after centrifugation at 25,000 g for 20 mins. The supernatants were carefully removed and mixed with an appropriate volume of cold acetone and stored at − 20 °C overnight. The mixture was centrifuged again. The pellets were dissolved using lysis buffer. Next, 10 mM DTT (dithiothreitol) was added to the solution, and the solution was maintained at 56 °C for 1 h to reduce the disulphide bonds of peptides. Next, 55 mM IAM (iodacetamide) was added to the solution, and the solution was maintained in a darkroom for 45 min. An appropriate volume of chilled acetone was added into the solution, which was maintained at − 20 °C for 2 h. The pellet was dissolved with lysis buffer to obtain the protein solution. Final protein concentrations were determined using BCA kits according to the manufacturer’s instructions.

Protein digestion and iTRAQ labelling

Each 100 μg of protein was digested in trypsin solution (1:10) and incubated at 37 °C for 4 h. The digested peptides were labelled using iTRAQ reagents according to the manufacturer’s instructions. iTRAQ Reagent 6-Plex kit (Applied Biosystems, Foster City, CA, USA) was used for iTRAQ labelling according to the manufacturer’s protocol with some modifications. The iTRAQ labelling reagents 113, 114 and 115 were used to label three biological replicates from the control group; similarly, 116, 117 and 118 were used to label three biological replicates from UpHT.

Peptide fractionation via reverse phase chromatography

Prominence high-performance liquid chromatography (HPLC) system (LC-20AB; Shimadzu, Kyoto, Japan) that was connected to a reversed-phase column was used to perform the first dimension of peptide separation. The iTRAQ-labelled sample was dried using a speedvac and resuspended in 2 mL of 5% acetonitrile, and the pH was adjusted to approximately pH 9.8. The resuspended sample solution was injected into a Gemini C18 column (4.6 × 250 mm, 5-μm particles). The following separation gradient with buffer B (98% ACN, 0.1% FA) was used at a flow rate of 1000 μL/min: 5% of mobile phase B in 10 min, 5 to 35% of mobile phase B in 40 min, and 35 to 95% of mobile phase B in 1 min. The system was then maintained in 95% of mobile phase B in 3 min, which was decreased to 5% within 1 min before equilibrating with 5% solvent B for 10 min. Elution was monitored by measuring the absorbance at 214 nm, and fractions are collected every 1 min. The peptides were subjected to nanoelectrospray ionization followed by tandem mass spectrometry (MS/MS) in a Q EXACTIVE (Thermo Fisher Scientific, San Jose, CA) that was coupled online to the HPLC.

Identification and quantification of peptides and proteins

The MS/MS data were searched against a protein sequence database that was derived from the transcriptome database of U. prolifera and partion for other green algae proteins using the ProteinPilot™ software 4.2 (SCIEX) for peptide identification and quantification. The MS/MS spectra obtained were searched using the following user-defined search parameters. The MS/MS spectra were searched against a decoy database to estimate the false discovery rate (FDR) for peptide identification. For iTRAQ studies, 95% protein confidence level was used as the identification criterion.

Proteomics data processing

The raw MS/MS data were converted into “.mgf” files using ProteinPilot software (AB Sciex). Mascot version 2.3.0 (Matrix Sciences, London, UK) was used to search against the transcriptome database of U. prolifera and partion for other green algae proteins. To identify false positives, raw spectra from the actual database were compared with a generated database of random sequences. Only peptides with significant scores at the 95% confidence level were considered reliable and used for protein identification. Protein quantitative ratios were weighted and normalized relative to the median ratio in Mascot. Only proteins with significant quantitative ratios between the two treatments (p < 0.05) and with fold changes > 1.2 or < 0.83 were considered differentially expressed. Functional annotation and classification of proteins for the DEPs was conducted using the Blast2GO program ( [19]. KEGG pathway analysis was performed using ( The clustering of the heat map was conducted using Mev4.9. Enrichment analysis was conducted using the singular enrichment analysis (SEA) tool in the agriGO toolkit [20]. Subcellular localizations of proteins were determined using TargetP [21].

RNA isolation and qRT-PCR

The isolated RNA sequencing samples were also used to perform real-time quantitative (qRT-PCR) analysis. Extraction of total RNA from three biological replicates in UpHT and three biological replicates in UpC was performed using plant RNA Kit (Omega). The total RNA was reverse transcribed to cDNA using TransScript® All-in-One First-Strand cDNA Synthesis SuperMix for qPCR (one-step gDNA removal) (Trans) according to the manufacturer’s protocol. The purity and integrity of the RNA was verified via absorbance measurements at 260 nm using a Nanodrop ND-1000 spectrophotometer (NanoDrop Technologies Inc., USA) and via electrophoretic separation using a Bioanalyzer and the RNA Nano kit (Agilent Technologies). From the DEGs, 8 related genes were selected to verify the reliability of the transcriptome by real-time fluorescence quantitative PCR. Quantitative real-time PCR was performed on an ABI Quant Studio6 Flex (USA) using SYBR Green PCR kits (Trans, China) according to the manufacturers’ instructions. Each cycle of the thermal amplification followed the universal protocol according to the manufacturer’s instructions: 94 °C for 30 s, followed by 40 cycles with 15 s at 94 °C and 1 min at 55 °C. For all qRT-PCR results, β-actin was used as the housekeeping gene, for which sequences of the specific primers have been previously published. The gene-specific primers used to amplify each gene were as follows: β-actin forward primer was 5’-AGGATGCATACGTTGGTGAA-3’, and β-Actin reverse primer was 5’-TTGTGGTGCCAAATCTTCTC-3’ [8].

Western blot

The protein samples were electrophoresed via 12% SDS-PAGE and transferred onto nitrocellulose membranes (Millipore, Germany). After blocking with 5% skim milk, the membranes were incubated with the rabbit polyclonal anti-HSP70 (Agrisera AS0837, Sweden,1: 3000), Rabbit polyclonal Rubisco large subunit,form I and form II antibody (Agrisera AS03037, Sweden 1:10000) at 4 °C overnight, respectively. Subsequently, the membranes were washed with PBS five times and incubated with goat anti-rabbit IgG(H + L) secondary antibody (Thermo Fisher Scientific, USA) at a ratio of 1:5000 at room temperature for 2 h. Finally, the membranes were developed using an enhanced chemiluminescence kit (Santa Cruz, USA) [22].

PPI network analysis

The PPI data of U. prolifera were downloaded from the STRING database. Each interaction has a combined score, which represents the reliability of the interaction between the proteins. The PPI interactions with a combined score (0: lowest confidence; 1: highest confidence) larger than 0.4 were used for further network analysis. All differentially expressed proteins were mapped onto the PPI network, and the Cytoscape tool was used to visualize the network.

Physiological and biochemical indicators

It was proposed that physiological adaptations of U. prolifera may enable it to survive the harsh intertidal environment and contribute to subsequent blooms. To investigate the effects of high temperature on antioxidant activities and photosynthesis indexes, the U. prolifera samples were quickly frozen using liquid nitrogen and stored at − 80 °C. For antioxidative enzyme extraction, 1.0 g of U. prolifera was homogenized in 5.0 mL of extraction buffer containing 1 mM EDTA, 0.05% Triton-X-100, 2% PVP, and 1 mM ascorbate in 50 mM phosphate buffer, pH 7.8. This mixture was centrifuged at 12,000 g for 20 min at 4 °C [11, 23]. The resulting supernatant was stored at − 20 °C for the assay of the following antioxidant enzymes.

Superoxide dismutase activity determination

Total superoxide dismutase activity was determined via the inhibition of the photochemical reduction of the chloride nitroblue tetrazolium (NBT) at 560 nm. The enzyme activity was expressed as unit U/g fresh weight, and one SOD unit was defined as the quantity required to inhibit the photoreduction rate of NBT by 50% [24].

Ascorbate peroxidase (APX) activity determination

Ascorbate peroxidase (APX) activity was determined by measuring the decrease in absorbance of ascorbic acid (AsA) at 290 nm. The reaction mixture contained 2.6 mL of PBS (pH 7.5, containing 0.1 mmol/L EDTA and 0.5 mmol/L AsA), 0.1 mL of sample, and 0.3 mL of 2 mmol/L H2O2 [24, 25].

Estimation of lipid peroxidation (MDA content)

Fresh samples (500 mg each) were homogenized in 10 mL of 0.1% trichloroacetic acid (TCA). The homogenate was centrifuged at 15,000 g for 5 min. Next, 2 mL aliquot of supernatant was taken and 4 mL of 0.5% thiobarbituric acid (TBA) in 20% TCA was added to it. The mixture was heated at 95 °C for 30 min and quickly cooled in an ice bath. After centrifugation at 10,000 g for 10 min to remove the suspended turbidity, the absorbance of the supernatant was recorded at 532 nm absorbance on a UV-visible spectrophotometer (Chemito Spectrascan, UV 2600). The value of non-specific absorption at 600 nm was subtracted. The MDA content was calculated using its absorption coefficient of 155 mmol/ cm [26, 27].

Chlorophyll fluorescence parameters measurements

Chlorophyll fluorescence parameters of the maximum photochemical efficiency (Fv/Fm) were measured using Water-PAM. Detection of chlorophyll fluorescence was made using a pulse modulation fluorometer (JUNIOR-PAM, Walz, Germany). At least three algal samples were used for each measurement of chlorophyll fluorescence, and the algae were acclimated to darkness for 10 min before being analysed. The maximum quantum yield of the photosystem (PS) II of U. prolifera was estimated as FV/Fm; the photochemical quenching coefficient (qP) and nonphotochemical quenching coefficient (NPQ) were also determined. The rapid light curves (RLCs) consisted of the fluorescence response to eight different and increasing actinic irradiance levels over the range of 0  820 μmol photons m− 2 s− 1. The parameters of the RLCs were calculated following the formula described by Jassby and Platt: rETR (relative electron transport rate) = rETRmax × tanh( × I/ rETRmax), where rETRmax is the saturated maximum rETR, tanh is the hyperbolic tangent function, is the initial slope of the RLC (the efficiency of the electron transport), and I is the incident irradiance [28].

Results and discussion

Annotation of proteome data

Protein response to heat stress in U. prolifera were revealed by iTRAQ analysis. A total of 283,344 spectra were generated; 15,546 peptides and 4449 proteins were identified with the cut-off of Mascot Percolator Q value<= 0.01 [29]. Differentially expressed proteins (DEPs) were identified upon an expression ratio > 1.20 or < 0.83 and p < 0.05 under UpHT compared to UpC [30,31,32]. During high-temperature stress, a total of 1223 proteins were found to be significantly changed, of which, 790 DEPs were up-regulated and 433 DEPs were down-regulated (Fig. 1, Table 1), which indicated that almost 27.50% of the identified proteins changed their abundance significantly in response to the heat stress.

Fig. 1
figure 1

Identification of DEPs. a Volcano of differentially expressed proteins(DEPs)under UpHT comparing to UpC. This plot depicts volcano plot of log2 fold-change (x-axis) versus -log10 Qvalue (y-axis, representing the probability that the protein is differentially expressed). Qvalue<=0.05, Foldchange> 1.2 and < 0.83 are set as the significant threshold for differentially expression. Dots in green indicated down-regulated proteins, Dots in red indicated up-regulated proteins, black indicated no significantly change. b the number of DEPs (UpHT vs UpC)

Table 1 Partion DEPs under UpHT vs UpC conditions

GO function analysis of DEPs

GO analysis was conducted on DEPs using the Blast2GO program; 469, 481 and 473 DEPs were successfully assigned to the biological processes, molecular function, and cellular components GO categories, respectively. The distribution of GO terms (the second level) for the three categories is shown in Table 2. For the biological process category, the most abundant groups were metabolic process (42.20%), cellular process (39.40%), and translation (24.90%). In addition, 126 DEPs belonged to response to stimulus (26.86%), of which 86 DEPs were up-regulated and 40 DEPs were down-regulated. For the cellular component category, cells (63.59%), cell parts (63.59%), and organelles (51.79%) were the most abundant groups. In the molecular function category, most of the DEPs were classified into catalytic activity (58.07%) and binding functions (45.70%); antioxidant activity only accounted for 0.72% of the DEPs (Table 2). Furthermore, subcellular localization of the 790 up-regulated characterized proteins showed that 54 proteins (6.84%) were located in the chloroplast, 139 proteins (17.60%) were assigned to the mitochondria, 42 proteins (5.32%) belonged to the secretory pathway, and 512 proteins (12%) were classified as belonging to other locations. Forty-three of the DEPs had unknown locations (Fig. 2). Furthermore, subcellular localization analysis of the 433 down-regulated proteins showed that quite many chloroplast proteins and mitochondria proteins are related to the thermotolerance of U. prolifera.

Table 2 GO analysis of DEPs (UpHT vs UpC) (top15)
Fig. 2
figure 2

Subcellular localization of the differentiallyexpressed proteins under heat stress. a Subcellular localization of the up-regulated proteins b Subcellular localization of the down-regulated proteins Note: C: Chloroplast, i.e. the sequence contains cTP, a chloroplasttransit peptide; M: Mitochondrion, i.e. the sequence contains mTP, amitochondrial targeting peptide; S: Secretory pathway, i.e. the sequencecontains SP, a signal peptide; −: other location

COG function class of DEPs

In the COG function analysis, the largest category was “translation, ribosomal structure and biogenesis” (161 DEPs up-regulated, 13 DEPs down-regulated), followed by “post-translational modification, protein turnover, molecular chaperone” (79DEPs up-regulated and 20 DEPs down-regulated) and “general function prediction only” (39 DEPs up-regulated and 24 DEPs down-regulated) (Fig. 3).

Fig. 3
figure 3

COG function analysis DEPs of UpHT vs UpC. Colour in red indicated up-regulated proteins number, blue indicated down-regulated proteins number

KEGG pathway analysis of DEPs

KEGG pathway analysis of DEPs were conducted online at; 557 up-regulated proteins were mapped to 97 KEGG pathways, and the top 15 pathways are listed in Table 3. The pathway enrichment analysis showed that the most significantly enriched pathways among the 557 up-regulated proteins were ribosome (134, 24.05%), oxidative phosphorylation (32, 5.75%), proteasome (14, 2.51%), glutathione metabolism (15, 2.69%), SNARE interactions in vesicular transport (5, 0.89%) and linoleic acid metabolism (3, 0.53%) (Fig. 4a). Two hundred and thirty-three down-regulated proteins were enriched in photosynthesis (15, 6.44%), carbon fixation in photosynthesis organisms (15, 6.44%), photosynthesis antenna proteins (8, 8.43%), metabolic pathways (91, 39.06%) and lysine biosynthesis (3, 1.2%) (Fig. 4b). At present, glutathione S-transferase, heat shock proteins, MnSOD and ubiquitin-related proteins were up-regulated. The results indicated that upon high-temperature stimulus, U. prolifera induces defence mechanisms, and the expression of heat shock proteins and antioxidant-associated proteins is increased.

Table 3 KEGG pathway of DEPs(top15)(UpHT vs UpC)
Fig. 4
figure 4

Pathway enrichment statistical scatter plot of DEPs. a UpHT vs UpC up-regulated proteins b UpHT vs UpC down-regulated proteins. The vertical axis represents the name of the pathway; the horizontal axis represents the pathway corresponding Rich factor. The colour indicated p value. Rich factor refers to the ratio of the number of differentially expressed genes in the pathway and the number of all annotated genes in the pathway. Higher Rich factors indicate greater degrees of enrichment. Q values are often completed after multiple hypothesis testing with corrected P value values ranging from 0 to 0.05. The closer they are to zero, the more significant the enrichment

Expression analysis of DEPs in response to stimulus

Abiotic stresses such as high temperature, low temperature, and drought induced plant defence mechanisms, including the expression of antioxidant enzymes to regulate their adaptation [33]. In higher plants, the signal response to high-temperature stress involves a reduction in the synthesis of normal proteins and is accompanied by accelerated expression of heat-responsive genes and HSPs. Reactive oxygen species (ROS) refer to free radicals, including hydrogen peroxide (H2O2), singlet oxygen (O2), and hydroxyl radicals (OH). Excessive amounts of ROS can damage macromolecules and cell membranes. To remove toxic ROS, algae have developed a number of antioxidant systems that serve as protective mechanisms, namely, superoxide dismutases (SOD), catalases, peroxidases, thioredoxin, and glutathione [34]. In the present study, high-temperature stress can lead to protein misfolding and inactivation and desiccation of plant cells. Heat shock proteins are important molecular chaperones of the Hsp family, which prevents the aggregation of misfolded proteins. Significant up-regulation was observed for the heat shock protein family, which included 90 (CL5666_Contig1_All, 1.5-fold), 90C (Unigene40313_All, 1.78-fold), 70 (gi|219764948|emb|CAW63927_1|, 1.4-fold), ClpB chaperone, Hsp100 family (CL517_Contig2_All, 1.64-fold). Therefore, it could play a role in protecting the protein structure early in the stress response. DNA repair protein Rad50 (CL6208_Contig1_All, 1.22-fold) were identified and significantly up-regulated under high temperature. It means that the DNA repair protein Rad50 played an important role in DNA damage induced by high temperature. We also determined that antioxidant proteins exhibited highly dynamic changes in response to the high-temperature stimulus. Reactive oxygen or nitrogen species are generated in the plant cell during the extreme stress condition, which produces toxic compounds after reacting with the organic molecules. The glutathione-S-transferase (GST) enzymes play a significant role in detoxifying these toxins and help in their excretion or sequestration; the enzymes catalyse reduction of glutathione (GSH; a tripeptide Glu-Cys-Gly) by electrophilic and hydrophobic toxic molecules that are generated under stress to convert them to non-toxic and soluble conjugates [35]. The expression of glutathione S-transferase (Unigene15839_All,1.45-fold), MnSOD (CL6691_Contig2_All, 1.42-fold), ascorbate peroxidase (CL1454_Contig2_All,1.28-fold) and polyubiquitin (Unigene6725_All, 1.54-fold) were up-regulated. In addition, the expression of catalase-peroxidase (CL5101_Contig1_All, 0.66-fold), glutathione reductase (CL514_Contig1_All, 0.77-fold), and thioredoxin-dependent peroxidase (Unigene9950_All, 0.77-fold) were down-regulated.

In addition, the 26S proteasome regulatory subunit is a component of the hexameric ring of AAA-ATPases that forms the base of the 19S regulatory particle (RP). This subunit plays specific roles in plant proteasomes by helping to promote the assembly of the RP with the 20S core protease (CP) and gate the CP to prevent indiscriminate degradation of cytosolic and nuclear proteins. In plants, this subunit plays an important role in diverse processes that include shoot and root apical meristem maintenance, cell size regulation, and stress responses [36]. In the present study, the 26S proteasome regulatory subunit (CL1786_Contig3_All, CL2484_Contig2_All, 1.3 or 1.32-fold) was induced under high-temperature stress. The result is consistent with the observation that the expression of the 26S proteasome subunit RPN10 is upregulated by salt stress in Dunaliella viridis [37]. The result indicated that the 26S proteasome regulatory subunit played a key role in stress tolerance. Programmed cell death (PCD) is an important biological phenomenon. Algae can activate PCD when exposed to external pathogens and other stresses. Metacaspases often play an important role in plant PCD, and a large number of experiments have shown that metacaspases are involved in plant PCD. In PCD, type II metacaspase was found to be transported from the cytoplasm to the nucleus, which promoted chromosome degradation. In the present study, mov34-domain-containing protein (Unigene14048_All, 1.33-fold) and metacaspase type II (CL6485_Contig1_All, 1.25-fold) were up-regulated. The results indicated that the high-temperature stimulus accelerated the accumulation of ROS, inhibited the growth of U. prolifera, induced the expression of partial stress related proteins, and alleviated the stimulation. This process is consistent with the results that are reported for the brown alga Ectocarpus siliculosus [38].

Expression analysis of DEPs that are involved in the oxidative phosphorylation pathway

The oxidative phosphorylation (OXPHOS) process comprises an electron-transfer chain (ETC) that is driven by substrate oxidation and is coupled to the synthesis of ATP through an electrochemical transmembrane gradient. Oxidative phosphorylation is a vital part of metabolism, it also produces reactive oxygen species such as superoxide and hydrogen peroxide, which leads to propagation of free radicals, damages cells and contributes to disease, and possibly ageing. The previous study showed that low temperature, UVB radiation and phoxim-induced increase in oxidative phosphorylation gene [39,40,41]. In the present study, 34DEPs were mapped to the oxidative phosphorylation pathway (Figs. 5 and 6). The expression of 32 DEPs were found to be increased, and 2 DEPs were down-regulated. NADH:ubiquinone oxidoreductase (CL6587_Contig2_All), succinate dehydrogenase (CL2133_Contig2_All), cytochrome C oxidase subunit (CL1977_Contig2_All), and ATP synthesis subunit (gi|807046094|gb|AKC35152_1) were all significantly up-regulated by 1.36 fold, 1.59 fold, 1.23 fold, and 1.83 fold, respectively. The results are consistent with previous results. These results indicated that high-temperature induction enhanced the activity of the ETC complex in the mitochondria. The imbalance between the activity changes of ETC may lead to the accumulation of reactive oxygen species.

Fig. 5
figure 5

The oxidative phosphorylation pathway

Fig. 6
figure 6

The cluster of oxidative phosphorylation pathway(UpHT vs UpC). Red indicated the expression of protein up-regulated, green indicated the expression of protein down-regulated

Expression analysis of DEPs involved in phytohormone and signal transduction pathways

Studies have investigated the expression of the hormone-related proteins that were changed under abiotic stress conditions [42]. In the present study, expression pattern analysis showed that the high-temperature stimulus changed the expression of key enzymes involved in phytohormone signal transduction and hormone biosynthesis. Approximately 12 DEPs were identified in phytohormone signal transduction (Fig. 7). In carotenoid biosynthesis, cytochrome P450 (Unigene18531_All) was significantly up-regulated by 1.23 fold. The expression of violaxanthin de-epoxidase-related protein (CL1646_Contig1_Alland CL1850_Contig11_All), prolycopene isomerase (Unigene205_All) and zeaxanthin epoxidase (CL4385_Contig1_All) were down-regulated by 0.83 fold, 0.74 fold, 0.81 fold, and 0.82 fold compared to the control. In brassinosteroid (BR) biosynthesis, molecular and biochemical studies revealed that two different cytochrome P450s, CYP90B2 and CYP724B1/D11, are the rate-limiting step of brassinosteroid biosynthesis. Three significant DEPs that were mapped to the BR synthesis pathway the expression of cytochrome P450s 90B/724B (gi|158273526|gb|EDO99315_1|, 0.69-fold) were decreased. The expression of 90C1D1 (CL3922_Contig1_All, 1.29-fold) and CYP735A (Unigene18531_All, 1.23-fold) were up-regulated. CYP735A encodes cytokinin hydroxylases that catalyse iP-nucleotide converts to trans-Zeatin. The signalling transduction of the cytokinin pathway and abscisic acid (ABA) pathway played an important role in the growth and development of plants and the regulation of stress resistance. AHP (Histidine-containing phosphotransfer protein) mainly functions as two-component phosphorelay mediators between cytokinin sensor histidine kinases and response regulators (B-type ARRs) and plays an important role in promoting cytokinin signal transduction via the multistep His-to-Asp phosphorelay. Biochemical and molecular genetic studies in Arabidopsis have identified protein phosphatase 2C (PP2C) enzymes as negative regulators in plant signal transduction processes [43]. In the absence of ABA, PP2C mainly inhibited the active state of SnRK2s protein. The inhibition of the active state of SnRK2s protein result in a decrease in downstream transcriptional regulation. ABA relieved the inhibition effect of PP2C on SnRK2s, activated the SnRK2s protein, and initiated downstream signal transduction. In the present study, the expression of the AHP protein (CL7329_Contig1_All, 1.25-fold) and SnRK2s (CL4735_Contig1_All, 1.27-fold) protein was up-regulated. However, the expression of PP2C (CL5191_Contig2_All, 0.8-fold) was down-regulated, which indicated that the signal transduction pathway of cytokinin and ABA was induced at a high temperature. The result is consistent with previous rzesearch [42, 43]. Plant growth is co-ordinately regulated by environmental and hormonal signals. Brassinosteroid (BR), cytokinin and ABA played essential roles in growth regulation via light and temperature; however, the interactions between hormone signal transduction and environmental signals remain poorly understood at the molecular level. In summary, the results indicated that phytohormone signal regulation played a key role in the high-temperature stimulus and enhanced cytokinin and ABA signal transduction pathway. However, the regulation mechanism should be further considered.

Fig. 7
figure 7

The cluster of plant hormone biosynthesis and signal transduction related proteins(UpHT vs UpC). Red indicated the expression of protein up-regulated, green indicated the expression of protein down-regulated

Expression analysis of DEPs involved in photosynthesis-related proteins

Photosynthesis is very sensitive to high-temperature stress. Changes in environmental temperature are primarily reflected by photosynthesis, which triggers a response aimed at attaining the best possible performance under the new conditions. Several studies have shown that high-temperature stress can significantly inhibit photosynthesis [44, 45]. In green algae, light-harvesting chlorophyll protein complexes are the major light-harvesting complexes. These complexes include photosystem II (light-harvesting complex II and photosynthetic reaction centre) and photosystem I (LHCI and photosynthetic reaction centre). The light-harvesting complex is a complex of subunit proteins that may be part of a larger supercomplex of the photosystem, and the function is to collect more of the incoming light. In the previous study, LHCI of photosystem I Chlamydomonas reinhardtii is composed of nine genes. Three product complexes (Lhca1, 3 and 7) showed emission peaks in the range of 680–690 nm, thus accounting for the blue contribution observed in the native preparations. Three other complexes (Lhca5, 6 and 8) exhibit emission in the 695–700 nm range. This emission component is not clearly resolved in the fluorescence spectra of the native LHCI preparations, possibly due to energy transfer to the red forms of a neighbouring complex or due to a change in the absorption properties of these complexes, because of protein–protein interactions, when embedded in the PSI–LHCI supercomplex, as previously proposed [46]. Three more complexes, namely, Lhca2, 4 and 9, exhibited emission above 707 nm with the red-most form (715 nm) associated with Lhca2, which indicates that Lhca2 is responsible for the red-most emission in vivo. Up-regulation of Lhca4 and Lhca9 was observed under iron-deficiency conditions, and it correlates with a shift to the blue end of the emission spectrum (from 710 to 704 nm) and with a relative increase in intensity of the emission above 700 nm, compared to that for the emission below 700 nm [47]. This indicates an increase in the content of red forms that emit approximately 705 nm in the antenna, in agreement with the fluorescence characteristics of Lhca4 and Lhca9. Under the same conditions, down-regulation of Lhca5 was also detected, leadin to the suggestion that Lhca4 can substitute for Lhca5, which possibly provides a better sink for energy dissipation. The finding that these antennas contain the lowest energy forms supports the hypothesis. A similar effect was recently observed in higher plants, where Lhca5 could substitute for Lhca4 [48, 49]. In the present study, the LHCI subunits lhca1 (CL4203_Contig1_All, 0.72-fold) and lhca3 (CL654_Contig2_All, 0.8-fold) were down-regulated, and lhca2 (LhcSR, gi|315319013|gb|ADU04518_1|, 1.3-fold) was up-regulated. It may be possible that the decrease in the emission of blue light results in energy transfer to the red forms. Lhca2 is responsible for the red-most emission. The up-regulation of lhca2 indicated that the lowest energy forms were native organized. Lhca1 and lhca3 were unstable and underwent decomposition under high temperature. The results were consistent with that obtained by Naumann et al. [47]. In addition, the two most abundant light-harvesting complex II are Lhcb1 and Lhcb2, which make up light-harvesting complex (LHC) II trimers. They are also involved in facilitating state transitions, a process during which energy balance between photosystem (PS) II and I is achieved. In the present study, the LHCII subunits lhcb1 and lhcb2 were also down-regulated. Down-regulation of LHCII would significantly influence light absorbance and energy transfer in U. prolifera upon a high–temperature stimulus (Fig. 8). The results were consistent with the observation that salt stress induces a decrease in excitation energy transfer from phycobilisomes to photosystem II but an increase to photosystem I in the cyanobacterium Spirulina platensis [50]. Eighteen DEPs were mapped to the photosynthesis pathway, and the expression of 14 DEPs were down-regulated. Ferredoxin-NADP+ reductase (FNR) is one of the important enzymes involving in many biochemical and physiological metabolism processes, such as electron transfer, carbon dioxide fixation, nitrogen assimilation and antioxidation. Ferredoxin-NADP+reductase (Unigene9738_All, 0.82-fold and CL7087_Contig1_All, 0.81-fold) and photosystem I and photosystem II subunit proteins were identified and found to be decreased.

Fig. 8
figure 8

The expression of photosynthesis-related proteins (UpHT vs UpC). Red indicated the expression of protein up-regulated, green indicated the expression of protein down-regulated

However, the expression of the chloroplast ATP synthase subunit that was encoded by 5 unigenes (gi|807,046,094|gb|AKC35152_1|, 1.82-fold) was up-regulated.

This result indicated that the high temperature inhibited photosynthesis of U. prolifera and induced stress-related LhcSR protein expression, which was consistent with the expression of LhcSR in U. linza in a study that was conducted by Dong et al. [51]. The electrons from PSI and PSII might transfer to oxygen, and thereby lead to substantial generation of reactive oxygen species. Enrichment of the oxidative phosphorylation pathway under cold treatment suggested that high–temperature induced ROS production might occur. However, this hypothesis still lacks experimental evidence.

Expression analysis of DEPs involved in the carbon-fixation pathway

According to CO2 assimilation mechanisms, plants can be classified into three types: C3, C4, and CAM. C4 plants have acquired a series of high-performance photosynthetic genes during evolution, which confer more efficient photosynthesis than that in C3 plants under adverse conditions such as high light intensity, temperature, and low CO2 concentration. In the present study, 19 DEPs were mapped to carbon fixation in photosynthetic organisms (Fig. 9). The expression of 14 DEPs was down-regulated, and that of 5 DEPs was up-regulated. In C4-dicarboxylic acid cycle, the key regulated enzymes aspartate aminotransferase (Unigene42677_All, 0.71-fold), pyruvate orthophosphate dikinase (Unigene19826_All, 0.73-fold), and NADP malic enzyme (CL5095_Contig3_All, 0.81-fold) were down-regulated. At the same time, pyruvate orthophosphate dikinase (Unigene19826_All, 0.73-fold) and NADP malic enzyme were also found to be down-regulated in CAM. Moreover, the expressions of plastid transketolase (Unigene20747_All, 0.43-fold), fructose-1,6-bisphosphatase (CL438_Contig1_All, 0.81-fold), glyceraldehyde-3-phosphate dehydrogenase (CL7324_Contig1_All, 0.72-fold) and phosphoglycerate kinase (CL1803_Contig1_All, 0.81-fold) were down-regulated. Ribulose-1,5-bisphosphatecarboxylase/oxygenase large subunit was found to be encoded by 7 DEPs, of which, 4 DEPs were up-regulated, and 3 DEPs were down-regulated. The result demonstrated that the expression of the key enzyme in the carbon-fixation pathway was down-regulated.

Fig. 9
figure 9

The cluster of carbon fixation related proteins (UpHT vs UpC). Red indicated the expression of protein up-regulated, green indicated the expression of protein down-regulated

The activation of rubisco in vivo requires the presence of the regulatory protein rubisco activase. This enzyme facilitates the release of sugar phosphate inhibitors from rubisco catalytic sites thereby influencing carbamylation. Rubisco activase was used to explore the role of the enzyme in C4 photosynthesis at a high temperature. It is concluded that during short-term treatment at 40 °C, rubisco activase content is not the only factor that modulates rubisco carbamylation during C4 photosynthesis [52]. In the present study, rubisco activase (CL649_Contig6_All, 1.32-fold) was up-regulated; however, it is not the only factor. At the same time, soluble carbonic anhydrase precursor encoded by CL7120_Contig2_All was down-regulated, and CL6328_Contig1_All was up-regulated. In other words, to protect the energy balance, U. prolifera decreased photosynthesis and carbon fixation of metabolic activity under high-temperature conditions. However, increase in the expression of carbonic anhydrase precursor and rubisco activase play an important role in the acceleration to photosynthesis balance adjustment.

Expression analysis of DEPs involved in the glutathione metabolism pathway

A total of 18 DEPs were mapped to the glutathione metabolism pathway. The expression of key enzymes in glutathione metabolism pathway were observed to significantly change. Glutathione S-transferase (GST) catalyses the nucleophilic fusion of reduced glutathione (GSH; a tripeptide Glu-Cys-Gly) with electrophilic and hydrophobic toxic molecules, which were generated under stress, to convert them to non-toxic and soluble conjugates. GSTs play an important role in various metabolic pathways, and they are involved in detoxification of oxidative lipid peroxide metabolites, hormone metabolism, stress resistance, protection of cell oxidative stress injury, cell signal transduction and disease resistance [53,54,55]. Ascorbate peroxidase (APXs) catalyses the H2O2-dependent oxidation of ascorbate in plants, algae and certain cyanobacteria. Glutathione reductase (GR) catalyses the reduction of glutathione disulphide (GSSG) to the sulphydryl form glutathione (GSH), which is a critical molecule in resisting oxidative stress and maintaining the reducing environment inside the cell [56,57,58]. In the present study, APX (CL1454_Contig2_All, 1.28-fold), glutathione S-transferase (CL2243_Contig1_All, 1.53-fold), glucose-6-phosphate 1-dehydrogenase 2(CL2243_Contig1_All, 1.53-fold) and 6-phosphogluconate dehydrogenase (gi|158274793|gb|EDP00573_1|, 1.32-fold), decarboxylating(gi|158274793|gb|EDP00573_1|, 1.32-fold) were found to significantly increase. However, the expression of GR(CL514_Contig1_All, 0.77-fold) was down-regulated (Fig. 10). It is indicated that high temperature stimulus accumulated more ROS, induced the expression of resistance-related protein glutathione S-transferase and APX.

Fig. 10
figure 10

The cluster of glutathione metabolism related proteins (UpHT vs UpC). Red indicated the expression of protein up-regulated, green indicated the expression of protein down-regulated

Expression analysis of Ca2+-binding protein

Ca2+-dependent signalling processes are beneficial to plant perception, and they respond to diverse environmental stressors, such as osmotic stress and high temperature [59,60,61]. In our study, the expression of Hsp90 (Unigene40313_All, 1.78-fold), CDPK (Unigene26962_All, 1.35-fold) and EDS1 (gi|158,282,407|gb|EDP08159_1|, 1.4-fold) was up-regulated in U. prolifera under the high-temperature stimulus. In addition, CDPK can activate H2O2 accumulation through regulating NADPH oxidase. The results indicated that the high-temperature stimulus accelerated ROS signalling, CDPK and EDS1 expression. In conclusion, the Ca2+-signal also played an important role in high-temperature stress in U. prolifera.

Expression analysis of DEPs involved in other types of metabolisms

Linoleic acid (LA), a carboxylic acid, belongs to one of the two families of essential fatty acids. It is converted by various lipoxygenases and cyclooxygenases and the cytochrome P450 enzyme. All of these LA products exhibit bioactivity, and they are implicated in human physiology and pathology as indicated in the cited linkages. In the present study, 4 DEPs were mapped to linoleic acid metabolism. The expression of the alpha-beta hydrolase superfamily (gi|158279131|gb|EDP04893_1|, 2.05-fold), cytochrome P450 (Unigene18531_All, 1.23-fold), lipoxygenase (Unigene17329_All, 3.21-fold), and 12-oxophytodienoic acid reductase (CL2323_Contig1_All, 1.47-fold) was up-regulated. In addition, glycan, taurine, and hypotaurine played a key role in plant response to stress. The expression of glycan-biosynthesis-related proteins changed significantly under high-temperature stimulus. The expression of dolichyl-phosphate beta-glucosyl transferase (ALG5, Unigene5913_All, 1.29-fold) was up-regulated. However, alpha-1,2-mannosyltransferase (ALG9, CL1555_Contig1_All, 0.78-fold) was down-regulated. Glycosyl transferase is encoded by two DEPs; one is down-regulated, and the other is up-regulated. Moreover, it has been suggested that taurine, hypotaurine and their metabolic precursors (cysteic acid, cysteamine and cysteine sulphinic acid) might act as antioxidants through scavenging ·OH, HOCl and H2O2. It is concluded that cysteamine and hypotaurine are far more likely to act as antioxidants in vivo than taurine, provided that they are present in a sufficient concentration at sites of oxidant generation [62]. In our study, the expression of cystathionine beta-synthase (Unigene42486_All, 1.47-fold) was observed to be up-regulated. This result indicates that the high-temperature stimulus induced cysteamine biosynthesis and increased the antioxidant levels.

qPCR to verify the reliability of the proteome

The expression of 8 unigenes were investigated via qPCR. The results were consistent with the available proteome data (Additional file 1: Table S1). Quantitative real time polymerase chain reaction results indicated that the unigenes that encode heat shock protein 90, APX, MnSOD, glutathione reductase, glutathione S-transferase, RbcL, heat shock protein 70 and ATP synthase were all significantly up-regulated by 3.152, 2.502, 2.106, 0.420, 2.320, 2.522, 2.242 and 3.252-fold, respectively.

Validation of differentially expressed proteins identified by western blot

Two proteins, Hsp70 and RbcL identified DEPs with marked differences in expression determined by iTRAQ based quantitative analysis were selected to be verified by western blot analysis (Fig. 11). Hsp70 and RbcL protein was significantly up-regulated in UpHT groups as compared with control group (p < 0.05). The results which were found by western blot is consistent with the findings in iTRAQ analysis. It is well established that Hsp70 and RbcL could play a role in protecting the protein structure and photosynthesis early in the stress response.

Fig. 11
figure 11

Western blotting.showing the changes of Hsp70 and RbcL level in UpHT compare with the control UpC. Note: 1.UpC(Control), 2. UpC(Control), 3.UpHT (High temperature), 4. UpHT (High temperature)

PPI network

In plants, proteins do not function in isolation within cells but as part of a network. In this study, a PPI network was generated to highlight the interactions and relationships between different proteins. To obtain interaction between proteins, we constructed a PPI network that was based on data downloaded from the STRING database. For the UpHT vs UpC protein, we chose 102 proteins that were involved in response to stimulus, photosynthesis, carbon fixation in photosynthetic organisms, oxidative phosphorylation, glutathione metabolism, plant signal transduction and plant-pathogen interaction to construct the PPI. Based on the differential expression patterns at the protein levels in UpHT compared to UpC (Fig. 12). In the present study, there are 37 unigenes nodes and 194 interactions in the PPI network. The PPI network indicated that heat shock protein 90C (CL2081_Contig1_All, 1.49-fold), heat shock protein 90C (gi|158270891|gb|EDO96722_1|, 1.44-fold), heat shock protein 70 (gi|219764948|emb|CAW63927_1|, 1.4-fold), ClpB chaperone, Hsp100 family (CL517_Contig2_All and Unigene15739_All, 1.64-fold and 1.46-fold) and heat shock protein 70B (gi|297592060|gb|ADI46845_1|, 1.82-fold) are involved in interactions with proteins from other functional groups and played important signal regulated. In addition, soluble inorganic pyrophosphatase 2 (CL408_Contig1_All, 1.51-fold), plastid transketolase (Unigene20747_All, 0.43-fold), F1F0 ATP synthase gamma subunit (Unigene20819_All, 1.69-fold), ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit (gi|323,320,892|gb|ADX36432_1|, 0.8-fold), serine/threonine-protein kinase CTR1 (gi|300257667|gb|EFJ41912_1|, 2.66-fold), component of cytosolic 80S ribosome and 40S small subunit (CL2825_Contig3_All, 1.3-fold), flagellar associated protein (CL4919_Contig2_All, 1.53) as the key regulation factor and construct PPI under high temperature.

Fig. 12
figure 12

PPI network of U. prolifera response to high temperature stimuls. Red indicated the expression of protein up-regulated,blue indicated the expression of protein down-regulated. CL2081_Contig1_All: heat shock protein 90C; gi|158270891|gb|EDO96722_1|:heat shock protein 90C; gi|219764948|emb|CAW63927_1|:heat-shock inducible Hsp70;Unigene20819_All:F1F0 ATP synthase gamma subunit; CL517_Contig2_All:ClpB chaperone, Hsp100 family; gi|1001185371|gb|AML80572_1|:NADH-ubiquinone oxidoreductase 49 kDa subunit; CL4445_Contig2_All:prohibitin; gi|315319013|gb|ADU04518_1|:LhcSR; gi|300264644|gb|EFJ48839_1|:Stress-induced-phosphoprotein; CL408_Contig1_All:soluble inorganic pyrophosphatase 2; gi|300257667|gb|EFJ41912_1|:Serine/threonine-protein kinase CTR1;CL2825_Contig3_All:component of cytosolic 80S ribosome and 40S small subunit;Unigene15739_All:ClpB chaperone, Hsp100 family; CL4169_Contig1_All:nucleolar protein, component of C/D snoRNPs; Unigene20747_All:plastid transketolase; Unigene24720_All:DNA binding helix-turn helix protein; CL2133_Contig2_All:iron-sulfur subunit of mitochondrial succinate dehydrogenase;gi|297,592,060|gb|ADI46845_1|:heat shock protein 70B;CL4919_Contig2_All:flagellar associated protein;gi|300260825|gb|EFJ45042_1|:component of cytosolic 80S ribosome and 60S large subunit; CL7329_Contig1_All:histidine-aspartic acid phosphotransferase 1;CL4832_Contig3_All:ATP-dependent Clp protease regulatory subunit ClpC

A pathway model of high temperature stress responses in U. prolifera

Using our results and previous studies, we propose a putative synergistic response network for U. prolifera that responds to high-temperature stress. As shown in Fig. 13, the reaction of U. prolifera to high-temperature stress rapidly accumulates ROS, which leads to a series of metabolic changes. First, high-temperature stress induces stress signals to generate differential accumulation of signal transduction components. High temperature causes the cell ROS concentration to rapidly increase, and high concentrations of ROS result in an increasing release of Ca2+ into the cytoplasm. When the concentration of Ca2+ reaches a particular threshold, CDPK is activated. As a result, plants are able to produce appropriate responses to the stress conditions. In an emergency response, HSP, as the key molecule, is significantly up-regulated to protect proteins from misfolding. To remove the accumulated ROS, the antioxidant system is activated, and proteins that regulate the redox balance, including MnSOD, GSTs, APX and other response factors, are up-regulated. These proteins can break down excessive intracellular ROS and restore the normal intracellular redox environment. Proteins with important functions in energy metabolism are continuously upregulated to generate ATP, which allows plants to perform vital metabolic functions and to combat stress. Oxidative phosphorylation is the major approach to generate ATP, accompanied by the accumulation of ROS. The high-temperature stimulus induced upregulation of electron transport chain complex protein: NADH:ubiquinone oxidoreductase, succinate dehydrogenase, cytochrome C oxidase subunit, and ATP synthesis subunit in mitochondria. The imbalance between the activity changes of ETC may lead to accumulation of reactive oxygen species. The possible reason was that changes in the quinone pool redox state are responsible for the apparent inorganic phosphate activation of complex III. Complex III is responsible for higher ROS production during physiological working conditions relative to complex I [63]. At the same time, protein synthesis increased, and cytosolic 80S ribosome and 40S small subunit, nucleolar protein, and component of C/D snoRNPs were all up-regulated. However, the effect of high temperature on proteins involved in photosynthesis and carbon fixation are downregulated, and the result was consistent with the proteome response of Brachypodium distachyon to drought stress [64].

Fig. 13
figure 13

A pathway model of high temperature stress responses in U. prolifera. MAPK: mitogen-activated protein kinase; CNGC:cyclic-nucleotide gated channel;CDPK:Calcium-dependent protein kinases;ROS: reactive oxygen species;APX: ascorbate peroxidase;GST:glutathione S-transferase;MnSOD: manganese superoxide dismutase

Biological index and photosynthesis index

Antioxidant enzymes played an important role in preventing the oxidation of biological molecules and protecting plant tissues from damage caused by ROS that was induced by multiple environmental stresses. SODs catalysed O2− + 2H+ → O2 + H2O2. In addition, APX and CAT catalysed H2O2 → O2 + H2O. In the present study, the effects of high temperature were investigated on the activity of antioxidant enzymes, the chlorophyll content, and the chlorophyll fluorescence parameter of U. prolifera. The result indicated that under a high-temperature stimulus, the protein content decreased during the all process compared to the control group (Additional file 1: Figure S1). SOD activity was decreased under high temperature compared to the control group. However, the APX activity increased in the later period (Additional file 1: Figure S2, Fig. 4). The MDA content increased over the time and reached the highest value after 96 h, which was significantly different from the control group (Additional file 1: Figure S3). Furthermore, the value of Fv/Fm dropped to the lowest point after 96 h of culture (Additional file 1: Figure S5). Our results provide important information for the selection of heat-resistant U. prolifera.


A total of 1223 differentially expressed proteins were identified in response to heat stress. The up-regulated proteins were enriched in ribosome, oxidative phosphorylation, proteasome, glutathione metabolism, SNARE interactions in vesicular transport and linoleic acid metabolism pathway. The down-regulated proteins were in photosynthesis, carbon fixation, photosynthesis antenna proteins, and metabolic pathways. The up-regulated proteins were mainly glutathione S-transferases, heat shock proteins, manganese superoxide dismutase, ubiquitin-related protein, lhcSR, rubisco activase, serine/threonine protein kinase 2, adenylate kinase, Ca2+-dependent protein kinase (CDPK), disease resistance protein EDS1, mov34-domain-containing protein, metacaspase type II, NDPK2a, 26S proteasome regulatory subunit, ubiquinone oxidoreductase, ATP synthase subunit and cytochrome C oxidase subunit AHP, SnRK2s, and cytochrome P450. The down-regulated proteins were photosynthesis-related proteins, glutathione reductase, thioredoxin, thioredoxin peroxidase and carbon-fixation-related proteins. Overall, stress response, oxidative phosphorylation, cytokinin signal transduction pathway, abscisic acid signal transduction pathway, and Ca2+ signal transduction played key roles in high-temperature regulation. Moreover, high temperature might lead to photosynthesis and inhibition of carbon fixation in U. prolifera. These novel proteins provide a good starting point for further research into their functions using genetic or other approaches. These findings significantly improve the understanding of the molecular mechanisms that are involved in the tolerance of algae to high-temperature stress.


  1. Li YX, Zhang XW, Dong X. Differential gene expression in Ulva prolifera under low light and low temperature conditions. Curr Genet. 2012;58(4):235–44.

    Article  CAS  Google Scholar 

  2. Wang Y, Wang Y, Zhu L. Comparative studies on the ecophysio-logical differences of two green tide macroalgae under controlled laboratory conditions. PLoS One. 2012;7(8):1–10.

    Article  Google Scholar 

  3. Wahid A, Gelani S, Ashraf M, Foolad MR. Heat tolerance in plants: an overview. Environ Exp Bot. 2007;61(3):199–223.

    Article  Google Scholar 

  4. Howarth CJ. Genetic improvements of tolerance to high temperature. In: Ashraf M, PJC H, editors. Abiotic stresses: plant resistance through breeding and molecular approaches. New York: Haworth Press Inc.; 2005.

    Google Scholar 

  5. Bokszczanin KL, Frargkostefanakis S. Perspectives on deciphering mechanisms underlying plantheat stress response and thermotolerance. Front Plant Sci. 2013;4(315):315–20.

    PubMed  PubMed Central  Google Scholar 

  6. Liu GT, Ma L, Duan W, et al. Differential proteomic analysis of grapevine leaves by iTRAQ reveals responses to heat stress and subsequent recovery. BMC Plant Biol. 2014;14(1):110–7.

    Article  Google Scholar 

  7. Stirbet A. On the relation between the kautsky effect (chlorophyll a fluorescence induction) and photosystem II: basics and applications of the OJIP fluorescence transient. J Photoch Photobio. 2011;104(1–2):236–57.

    Article  CAS  Google Scholar 

  8. Fan MH, Sun X, Xu NJ, Liao Z, Wang RX. cDNA cloning, characterization and expression analysis of manganese superoxide dismutase in Ulva prolifera. J Appl Phycol. 2016;28(2):1391–401.

    Article  CAS  Google Scholar 

  9. Pang QT, Li F, Liu XQ, Wang JT. Compare the growth of Enteromorphaprolifera under different nutrient conditions. Environ Sci. 2013;34(9):3398–404.

    Google Scholar 

  10. Zhu M, Liu ZP, Xu JT, Mao YX, Yao DR. The release, adhesion and germination of spores of Enteromorpha prolifera and its adaptation to different desiccation conditions. Mar Sci. 2011;35(7):1–6.

    Google Scholar 

  11. Fan MH, Sun X, Liao Z, Li YH, Wang JX, Fan YP, Cui DL, Li P, Miao ZL. Integration of deep transcriptome and proteome analyses of salicylic acid regulation high temperature stress in Ulva prolifera. Sci Rep. 2017;7(1):11052–71.

    Article  Google Scholar 

  12. Bai YH, Wang S, Zhong H, Yang Q, Zhang F, Zhuang ZH, Yuan J, Nie XY, Wang SH. Integrative analyses reveal transcriptome-proteome correlation in biological pathways and secondary metabolism clusters in A. flavusin response to temperature. Scientific Reports. 2015;5:1–13. 5:14582.

    Article  CAS  Google Scholar 

  13. Bai XC, Song H, Lavoie M, Zhu K, Su YY, Ye HQ, Chen S, Fu ZW, Qian HF. Proteomic analyses bring new insights into the effect of a dark stress on lipid biosynthesis in Phaeodactylum tricornutum. Sci Rep. 2016;6:1–10.

    Article  CAS  Google Scholar 

  14. Sun PP, Mao YX, Li GY, Cao M, Kong FN, Wang L, Bi GQ. Comparative transcriptome profiling of Pyropia yezoensis (Ueda) M.S. Hwang & H.G. Choi in response to temperature stresses. BMC Genomics. 2015;16:463–79.

    Article  Google Scholar 

  15. Tia X, Chen L, Wang J, Qiao J, Zhang W. Quantitative proteomics reveals dynamicresponses of Synechocystis sp. PCC 6803 to next-generation biofuel butanol. J Proteome. 2013;78(1):326–45.

    Google Scholar 

  16. Qiao J, Wang J, Chen L, Tian X, Huang S, Ren X, Zhang W. Quantitative iTRAQ LC-MS/MS proteomics reveals metabolic responses to biofuel ethanol in cyanobacterial Synechocystis sp.PCC 6803. J Proteome Res. 2012;11(11):5286–300.

    Article  CAS  Google Scholar 

  17. Liu J, Chen L, Wang JX, Qiao J, Zhang W. Proteomic analysis reveals resistance mechanism against biofuel hexane in Synechocystis sp. PCC 6803. Biotechnol Biofuels. 2012;5(1):1–17.

    Article  Google Scholar 

  18. Provasoli L. Media and prospect for the cultivation of marine algae, in: Watanabe A, Hattori A (eds) Cultures and collection of algae. Jap Soc Plant Physiol 1968; 63–75.

  19. Conesa A, Götz S, García-Gómez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.

    Article  CAS  Google Scholar 

  20. Du Z, Zhou X, Ling Y, Zhang Z, Su Z. AgriGO: a GO analysis tool kit for the agricultural community. Nucleic Acids Res. 2010;38:W64–70.

    Article  CAS  Google Scholar 

  21. Emanuelsson O, Nielsen H, Brunak S, Heijne G. Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J Mol Biol. 2000;300(4):1005–16.

    Article  CAS  Google Scholar 

  22. Zhang Y, Wang YZ, Li SJ, Zhang XT, Li WH, Luo SX, Sun ZY, Nie RJ. ITRAQ-based quantitative proteomic analysis of processed Euphorbia lathyris L. for reducing the intestinal toxicity. Proteome Sci. 2018;16:8.

    Article  Google Scholar 

  23. Hussain S, Liu GQ, Liu DF, Ahmed M, Hussain N, Teng YW. Study on the expression of dehydrin genes and activities of antioxidative enzymes in oral buds of two sand pear (Pyrus pyrifolia Nakai) cultivars requiring different chilling hours for bud break. Turk J Agric For. 2015;39(6):930–9.

    Article  CAS  Google Scholar 

  24. Abd El Baky HH, Nofal OA, El Baroty GS. Enhancement of antioxidant enzymes activities, drought stress tolerances and quality of potato plants as response to algal foliar application. Recent Pat Food Nutr Agric. 2016;8(1):70–7.

  25. Zhang XY, Hu CG, Yao JL. Tetraploidization of diploid dioscorea results in activation of the antioxidant defense system and increased heat tolerance. J Plant Physiol. 2010;167(2):88–94.

    Article  CAS  Google Scholar 

  26. Shi JY, Gao LP, Zuo JH, Wang Q, Wang Q, Fan LL. Exogenous sodium nitroprusside treatment of broccoli florets extends shelf life, enhances antioxidant enzyme activity, and inhibits chlorophyll-degradation. Postharvest Biol Tec. 2016;116:98–104.

    Article  CAS  Google Scholar 

  27. Khare T, Desai D, Kumar V. Effect of MgCl2 stress on germination, plant growth, chlorophyll content, proline content and lipid peroxidation in sorghum cultivars. J Stress Physiol Biochem. 2012;8(4):169–78.

    Google Scholar 

  28. Wu HY. Effect of different light qualities on growth, pigment content, chlorophyll fluorescence, and antioxidant enzyme activity in the red alga Pyropia haitanensis (Bangiales Rhodophyta). Biomed Res Int. 2016;7383918:1–8.

    Google Scholar 

  29. Brosch M, Yu L, Hubbard T, Choudhary J. Accurate and sensitive peptide identification with mascot percolator. J Proteome Res. 2009;8(6):3176–81.

    Article  CAS  Google Scholar 

  30. Nguyen THN, Brechenmacher L, Aldrich JT, Clauss TR, Gritsenko MA, Hixson K, Libault M, Tanaka K, Yang F, Yao QM, Ljiljana PT, Xu D, Nguyen HT, Stacey G. Quantitative phosphorproteomic analysis ofsoybean root hairs inoculated with radyrhizobium japonicum. Mol Cell Proteomics. 2012;11(11):1140–55.

    Article  Google Scholar 

  31. Tolin S, Arrigoni G, Trentin AR, Veljovic-Jovanovic S, Pivato M. Biochemical and quantitative proteomics investigations in Arabidopsis ggt1 mutant leaves reveal a role for the gammaglutamyl cycle in plant's adaptation to environment. Proteomics. 2013;13(12–13):2031–45.

    Article  CAS  Google Scholar 

  32. Tolin S, Arrigoni G, Moscatiello R, Masi A, Navazio L, Sablok G, Squartini A. Quantitative analysis of the naringenin-inducible proteome in Rhizobium leguminosarum by isobaric tagging and mass spectrometry. Proteomics. 2013;13(12–13):1961–72.

    Article  CAS  Google Scholar 

  33. Kim HJ, Chiang YH, Kieber JJ, Schaller GE. SCF(KMD) controls cytokinin signaling by regulating the degradation of type-B response regulators. P Natl Acad Sci USA. 2013;110(24):10028–33.

    Article  CAS  Google Scholar 

  34. Wang JX, Milton S, Hu Q. Cloning and expression of isoenzymes of superoxide dismutase in Haematococcus pluvialis (Chlorophyceae) under oxidative stress. J Appl Phycol. 2011;23(6):995–1003.

    Article  CAS  Google Scholar 

  35. Tiwari V, Patel MK, Chaturvedi AK, Mishra A, Jha B. Functional characterization of the tau class glutathione-s-transferases gene (SbGSTU) promoter of Salicornia brachiata under salinity and osmotic stress. PLoS One. 2016;11(2):e0148494 1–20.

    Article  Google Scholar 

  36. Lee KH, Marshall RS, Slivicke LM, Vierstra RD. Genetic analyses of the Arabidopsis 26S proteasome regulatory particle reveal its importance during light stress and a specific role for the N-terminus of RPT2 in development. Plant Signal Behav. 2012;7(8):973–8.

    Article  CAS  Google Scholar 

  37. Sun XB, Meng XZ, Xu ZK, Song RT. Expression of the 26S proteasome subunit RPN10 is upregulated by salt stress in Dunaliella viridis. J Plant Physiol. 2010;167(12):1003–8.

    Article  CAS  Google Scholar 

  38. Dittami SM, Scornet D, Petit JL, Ségurens B, Silva CD, Corre E, Dondrup M, Glatting KH, König R, Sterck L, Rouzé P, de Peer YV, Cock JM, Boyen C, Tonon T. Global expression analysis of the brown alga Ectocarpus siliculosus (Phaeophyceae) reveals large scale reprogramming of the transcriptome in response to abiotic stress. Genome Biol 2009;10(6):R66.1-R66.20.

    Article  Google Scholar 

  39. Tian DQ, Pan XY, Yu YM, Wang WY, Zhang F, Ge YY, Shen XL, Shen FQ, Liu XJ. Denovo characterization of the Anthurium transcriptome and analysis of its digital gene expression under cold stress. BMC Genomics. 2013;14:827–40.

    Article  CAS  Google Scholar 

  40. Tsai ML, Chang KY, Chiang CS, Shu WY, Weng TC, Chen CR, Huang CL, Lin HK, Hsu IC. UVB radiation induces persistent activation of ribosome and oxidative phosphorylation pathways. Radiat Res. 2009;171:716–24.

    Article  CAS  Google Scholar 

  41. Li FC, Xu KZ, Ni M, Wang BB, Gu ZY, Shen WD, Li B. Effect of oxidative phosphorylation signaling pathway on silkworm midgut following exposure to phoxim. Environ Toxicol. 2015;32(1):167–77.

    Article  Google Scholar 

  42. Reguera M, Peleg Z, Abdel-Tawab YM, Tumimbang EB, Delatorre CA, Blumwald E. Stress-induced cytokinin synthesis increases drought tolerance through the coordinated regulation of carbon and nitrogen assimilation in rice. Plant Physiol. 2013;63(4):1609–22.

    Article  Google Scholar 

  43. Rodriguez PL. Protein phosphatase 2C (PP2C) function in higher plants. Plant Mol Biol. 1998;38(6):919–27.

    Article  CAS  Google Scholar 

  44. Salvucci ME, Crafts-Brandner SJ. Inhibition of photosynthesis by heat stress: the activation state of rubisco as a limiting factor in photosynthesis. Plant Physiol. 2004;120(2):179–86.

    Article  CAS  Google Scholar 

  45. Xu Y, Chen CS, Ji DH, Hang N, Xie C. Proteomic profile analysis of Pyropia haitanensis in response to high-temperature stress. J Appl Phycol. 2014;26(1):607–18.

    Article  CAS  Google Scholar 

  46. Bassi R, Soen SY, Frank G, Zuber H, Rochaix JD. Characterization of chlorophyll a/b proteins of photosystem I from Chlamydomonas reinhardtii. J Biol Chem. 1992;267:25714–21.

    CAS  PubMed  Google Scholar 

  47. Naumann B, Stauber EJ, Busch A, Sommer F, Hippler M. N-terminal processing of Lhca3 is a key step in remodeling of the photosystem I-light-harvesting complex under iron deficiency in Chlamydomonas reinhardtii. J Biol Chem. 2005;280:20431–41.

    Article  CAS  Google Scholar 

  48. Mozzo M, Mantelli M, Passarini F, Caffarri S, Croce R, Bassi R. Functional analysis of photosystem I light-harvesting complexes (Lhca) gene products of Chlamydomonas reinhardtii. Biochim Biophys Acta. 2010;1797:212–21.

    Article  CAS  Google Scholar 

  49. Wientjes E, van Stokkum IH, van Amerongen H, Croce R. Excitation-energy transfer dynamics of higher plant photosystem I light-harvesting complexes. Biophys J. 2011;100(5):1372–80.

    Article  CAS  Google Scholar 

  50. Zhang T, Gong HM, Wen XG, Lu CM. Salt stress induces a decrease in excitation energy transfer from phycobilisomes to photosystem II but an increase to photosystem I in the cyanobacterium Spirulina platensis. J Plant Physiol. 2010;167(12):951–8.

    Article  CAS  Google Scholar 

  51. Dong MT, Zhang XW, Zhuang ZM. Characterization of the LhcSR gene under light and temperature stress in the green alga Ulva linza. Plant Mol Biol Rep. 2012;30(1):10–6.

    Article  CAS  Google Scholar 

  52. Hendrickson L, Sharwood R, Ludwig M, Whitney SM, Badger MR, von Caemmerer S. The effects of rubisco activase on C4 photosynthesis and metabolism at high temperature. J Exp Bot. 2008;59(7):1789–98.

    Article  CAS  Google Scholar 

  53. Nishinaka T, IchijoY IM, Kimura M, Katsuyama M, Iwata K, Miura T, Terada T, Yabe-Nishimura C. Curcumin activates human glutathione S-transferase P1 expression through antioxidant response element. Toxicol Lett. 2007;170(3):238–47.

    Article  CAS  Google Scholar 

  54. Lyu K, Gu L, Li BP, Lu YC, Wu CC, Guan HY, Yang Z. Stress-responsive expression of a glutathione S-transferase (delta)gene in waterflea daphnia magna challenged by microcystin-producing and microcystin free microcystis aeruginosa. Harmful Algae. 2016;56:1–8.

    Article  Google Scholar 

  55. Pflugmacher S, Schwarz S, Pachur HJ, Steinberg CEW. Effects of tributyltin chloride (tbtcll) on detoxication enzymes in aquatic plants. Environ Toxicol2015;15(15):225–233.

    Google Scholar 

  56. Deponte M. Glutathione catalysis and the reaction mechanisms of glutathione-dependent enzymes. Biochim Biophys Acta. 2013;1830(5):3217–66.

    Article  CAS  Google Scholar 

  57. Meister A. Glutathione metabolism and its selective modification. J Biol Chem. 1988;263(33):17205–8.

    CAS  PubMed  Google Scholar 

  58. Vanacker H, Carver TLW, Foyer CH. Early H2O2 accumulation in mesophyll cells leads to induction of glutathione during the hyper sensitive response in the barley-powdery mildew interaction. Plant Physiol. 2000;123(4):1289–300.

    Article  CAS  Google Scholar 

  59. Nawaz Z, Kakar KU, Saand MA, Shu QY. Cyclic nucleotide-gated ion channel gene family in rice, identification, characterization and experimental analysis of expression response to planthormones, biotic and abiotic stresses. BMC Genomics. 2014;15(1):853–60.

    Article  Google Scholar 

  60. Jammes F, Hu HC, Villiers F, Bouten R, Kwak JM. Calciumpermeable channels in plant cells. FEBS J. 2011;278(22):4262–76.

    Article  CAS  Google Scholar 

  61. Qi Z, Verma R, Gehring C, Yamaguchi Y, Zhao YC, Ryan CA, Berkowitzb GA. Ca2+ signaling by plant Arabidopsis thaliana pep peptides depends on AtPepR1, a receptor with guanylyl cyclase activity, and cGMP-activated Ca2+ channels. P Natl Acad Sci USA. 2010;107(49):21193–8.

    Article  CAS  Google Scholar 

  62. Lodish H, Berk A, Zipursky S. Molecular Cell Biology, vol. 5. 7th ed. New York: Freeman WH and CO; 2013. p. 675.

    Google Scholar 

  63. Bazil JN, Beard DA, Vinnakota KC. Catalytic coupling of oxidative phosphorylation, ATP demand, and reactive oxygen species generation. Biophys J. 2016;110:962–7.

    Article  CAS  Google Scholar 

  64. Bian YW, Deng X, Yan X, Zhou JX, Yuan LL, Yan YM. Integrated proteomic analysis of Brachypodium distachyon roots and leaves reveals a synergistic network in the response to drought stress and recovery. Sci Rep. 2017;7(46183):1–15.

    Google Scholar 

Download references


Not applicable.


This work was sponsored by Zhejiang ocean university doctoral research startup funding project (Q1701), University student innovation training program(xj2018011), and the National Natural Science Foundation of China (41606129, 31671009, 41276122). This research was also sponsored by the K.C.Wong Magna Fund in Ningbo University.

Availability of data and materials

Please contact author for data requests.

Author information

Authors and Affiliations



NJX and MHF designed the study. MHF performed proteomics experiments. MHF, XS and ZL participated in the design of the study and performed the statistical analysis. MHF wrote the manuscript with the help of NJX and ZL and JXW. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Meihua Fan or Nianjun Xu.

Ethics declarations

Ethics approval and consent to participate

Samples of U. Prolifera were collected from the intertidal zone of the eastern gulf, Xiangshan, Ningbo, China (29°32.803 N, 121°45.599E). The Xiangshan municipal government is responsible for this area. The sampling sites were not located within protected areas; however, permission for collection was still obtained. Samples were handled in accordance with guidelines on experimental sample processing for scientific purposes set by of Zhejiang Ocean University, Zhoushan, China.

Consent for publication

Not applicable.

Competing interests

All the authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Effects of high temperature on protein content of U. prolifera. Figure S2. Effects of high temperature on superoxide dismutase activity of U. prolifera. Figure S3. Effects of high temperature on MDA contents of U. prolifera. Figure S4. Effects of high temperature on APX activity of U. prolifera. Figure S5. Effects of high temperature on the maximum quantum yield of U. prolifera. Table S1. qPCR validation of the proteome data. (DOCX 402 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fan, M., Sun, X., Liao, Z. et al. Comparative proteomic analysis of Ulva prolifera response to high temperature stress. Proteome Sci 16, 17 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: