- Research article
- Open Access
Comparative proteomic and metabolomic profiling of citrus fruit with enhancement of disease resistance by postharvest heat treatment
BMC Plant Biology volume 13, Article number: 44 (2013)
From field harvest to the consumer’s table, fresh citrus fruit spends a considerable amount of time in shipment and storage. During these processes, physiological disorders and pathological diseases are the main causes of fruit loss. Heat treatment (HT) has been widely used to maintain fruit quality during postharvest storage; however, limited molecular information related to this treatment is currently available at a systemic biological level.
Mature ‘Kamei’ Satsuma mandarin (Citrus unshiu Marc.) fruits were selected for exploring the disease resistance mechanisms induced by HT during postharvest storage. Proteomic analyses based on two-dimensional gel electrophoresis (2-DE), and metabolomic research based on gas chromatography coupled to mass spectrometry (GC-MS), and liquid chromatography quadrupole time-of-flight mass spectrometry (LC-QToF-MS) were conducted. The results show resistance associated proteins were up-regulated in heat treated pericarp, such as beta-1, 3-glucanase, Class III chitinase, 17.7 kDa heat shock protein and low molecular weight heat-shock protein. Also, redox metabolism enzymes were down-regulated in heat treated pericarp, including isoflavone reductase, oxidoreductase and superoxide dismutase. Primary metabolic profiling revealed organic acids and amino acids were down-regulated in heat treated pericarp; but significant accumulation of metabolites, including tetradecanoic acid, oleic acid, ornithine, 2-keto-d-gluconic acid, succinic acid, turanose, sucrose, galactose, myo-inositol, glucose and fructose were detected. Noticeably, H2O2 content decreased, while, lignin content increased in heat treated pericarp compared to the control, which might increase fruit resistibility in response to external stress. Also, flavonoids, substances which are well-known to be effective in reducing external stress, were up-regulated in heat treated pericarp.
This study provides a broad picture of differential accumulation of proteins and metabolites in postharvest citrus fruit, and gives new insights into HT improved fruit disease resistance during subsequent storage of ‘Kamei’ Satsuma mandarin. Interpretation of the data for the proteins and metabolites revealed reactive oxygen species (ROS) and lignin play important roles in heat treatment induced fruit resistance to pathogens and physiological disorders.
Fresh citrus fruits, from harvest to human consumption, require a lengthy period for shipping, storing and marketing. During that time, the main damage to fruit comes from biological diseases and physical damage. Storage pre-treatments, such as biotic or abiotic typological treatments, can be used to inhibit disease and extend the shelf-life. Among the various types of treatment available, heat treatment has previously been reported to induce disease resistance in apples, citrus, grapes, peaches, pears and mangos against a multitude of pathogens [1–3].
In the fruit industry, hot water treatment is an alternative physical treatment, which has been used in citrus fruit storage since 1922 . More recently, hot air and hot water systems have been widely adopted in several countries to inhibit postharvest decay of fruits and to alleviate chilling injuries during storage. Heat treatment (HT) has the advantage of minimizing surface infection in pericarp without leaving any chemical residues . HT can prevent the development of plant pathogens and disease by heat-killing disease organisms or by inhibiting pathogen mycelial growth, by removing surface insects and by structurally altering epicuticular wax . HT can also enhance fruit resistance to chilling injury in cold-sensitive cultivars, and help the fruit retain quality during cold storage and help maintain shelf life . Aside from controlling pathogenic organisms, HT may also reduce disease by changing the physical characteristics of the interior pericarp. For example, HT not only delayed fruit ripening and coloration, inhibited fruit softening, and slowed ethylene production , it also increased fruit hardness, as well as increasing the content of polyphenolics, flavonoids and dopamine [6–8]. However, less is understood about the regulatory genes which are active in the response to microbial attack, and there is no exhaustive metabolomics profile detailing the changes in heat treated pericarp.
In the past few decades, researchers have given increasing attention to the study of gene expression in heat treated fruits. Sapinitskaya et al. reported HT activates stress genes and lipid modification genes in response to chilling stress in grapefruit, such as glutathione-dependent formaldehyde dehydrogenase, dehydrin, heat shock protein (HSP), ascorbate peroxidase, superoxide dismutase, fatty acid desaturase and lipid transfer protein. Some scientists discovered HT induces phenylalanine ammonia lyase and accumulates phenylpropanoid compounds, which play a vital role in heat-induced fruit resistibility . Dixon et al. studied the functions of phenylpropanoid compounds in plants involved in local and systemic signalling for defence-gene induction. However, little information is available related to which particular class of phenylpropanoid compounds plays a core role in defensive functions, and even less is known about the regulatory genes which orchestrate the rapid induction of phenylpropanoid defences in response to HT. Proteome analysis reveals the induction of HSPs, allergen proteins, dehydrin, and other proteins involved in the stress response. Repression of polyphenol oxidase caused by HT may constitute the molecular basis for the protection against chilling stress in peach fruit . However, few proteins involved in HT were reported to have induced fruit resistance to deterioration, especially in citrus fruit.
Also, proteomics analysis and metabolic profiling approaches make possible comparative parallel analyses of global changes at the proteome and metabolome levels, facilitating an understanding of the relationships between changes in specific proteins and subsequent alterations in metabolites in response to HT. In the present study, fruits of ‘Kamei’ Satsuma mandarin (Citrus unshiu Marc.) received HT. The comparative proteomics analysis and metabolic profiling (the primary metabolome and the secondary metabolome) analysis were performed based on two-dimensional electrophoresis (2-DE), gas chromatography and mass spectrometry (GC-MS) and high-performance liquid chromatography/electrospray ionization-time of flight-mass spectrometry (HPLC-QTOF-MS). Comprehensive physiological and biochemical analysis were completed. Results will provide additional theoretical evidence to improvement of disease resistance in heat treated pericarp.
Satsuma mandarin is one of the main citrus crops in East Asia, including China, South Korea and Japan. However, this cultivar has poor storability mainly because it is susceptible to biotic stress and physical damage. HT improved fruit resistance to stress and damage effectively and inhibited postharvest decay of fruits. In present study, HT fruits were dipped in a 52°C warm water bath for 2 min, and control fruits were dipped in a 25°C water bath for 2 min. After treatments, fruits were stored in a ventilated warehouse (storage temperature: 12–16°C; relative humidity: 90–95%). Samples were collected at 2 h, and 1, 2, 3, 4, 6, 9, 12, 16, 20, 24, 28, 32, 38, 44, 50 and 60 days after treatments (DAT) for further analysis.
HT increased fruit weight loss slightly, but did not affect fruit quality
Fruit quality directly determines its commercial value. Parameters related to fruit quality, including weight loss, respiration rate, and soluble solids content, were evaluated to investigate the influence of HT on fruit quality. Weight loss increased significantly in heat treated pericarp compared to control pericarp during the storage period (Figure 1A). However, no obvious changes were observed in the respiration rate and soluble solids content between HT and control fruits during the entire storage period (Figures 1B, C). Since soluble solids content is a well-known standard in measuring fruit internal quality and respiration acts as an important standard for fruit storability, the minor changes in soluble solids content and respiration indicated that heat treatment did not affect fruit flesh quality.
HT suppression of Penicilliumgermination on pericarp
Artificial inoculation tests were conducted with the blue mould Penicillium italicum to assess the effects of HT on the fruit’s response to biotic stresses. Few spores had germinated on the peel of HT fruits 4 d after inoculation, but about 90% of the control fruits showed typical symptoms of blue mould development at the inoculation sites (Figure 2A). As for the rate of disease incidence, the incidence of disease on control fruit was significantly higher than on HT fruit from 4–7 d after inoculation (Figure 2B); lesion diameters were significantly higher in control fruit than HT fruit from 4–7 d after inoculation (Figure 2C). Based on the above experiments, HT effectively inhibited the germination of Penicillium on pericarp.
Differentially accumulated proteins in HT citrus pericarp
To explore the differential expression of proteins between HT and control at the same period, the pericarp of ‘Kamei’ Satsuma mandarin fruit were analysed with samples collected at 1, 6, 12 and 32 DAT. After 2-DE analysing, more than 600 protein spots were reproducibly detected (Figure 3). The protein relative abundance was obtained with the PDQuest 2-D software version 7.4. Protein spots were scored only when they were reproducibly observed in four independent replicates. Data were normalized as a percentage of the total density of all spots on the corresponding gel. Comparative analysis was carried out using 1.5-fold comparative analysis and one-way analysis of variance (ANOVA). Result showed there were only 50 differentially displayed protein spots whose abundance was altered at least 1.5-fold when comparing the HT and control during the same period of storage. Of these, 37 were identified using Applied Biosystems 4800 matrix-assisted laser desorption/ionization-time-of-flight tandem mass spectrometry (MALDI-TOF MS) based on tryptic peptide sequences (Figure 3, Table 1). Based on previously reported functions of these proteins, they were classified into five categories (Table 1): glycolysis and TCA cycle (nine proteins, 24%), redox regulation (nine proteins, 24%), stress response (seven proteins, 19%), other metabolism (six proteins, 16%), and protein folding (six proteins, 16%).
Differentially accumulated metabolites in HT citrus pericarp
Differentially accumulated metabolites were studied to investigate the changes in metabolic composition between HT and control pericarp. The biochemical changes were determined by GC-MS and HPLC-QTOF-MS. A parallel analysis of non-polar and polar extracts using methanol extraction enabled the detection of a large number of compounds of different classes.
A total of 62 identified metabolites were monitored in the same sample sets using GC-MS, including mainly alcohols, amino acids, sugars, organic acids, and fatty acids. After conducting an analysis of significance with ANOVA (P <0.05), 45 detected metabolites, mostly organic acids, sugars and amino acids, were found to be significantly up/down-regulated in heat treated pericarp when compared to control pericarp with the same storage period (Table 2).
GC-MS metabolite profiling (primary metabolites) was complemented using HPLC-QTOF-MS analysis of mostly secondary metabolites. Statistical filtering was applied to the mass signals data to estimate the number of up/down-regulated metabolites in heat treated pericarp. After MS/MS analysis, a total of 58 metabolites were found to be significantly regulated (2-fold; P < 0.05) in heat treated pericarp when compared to control pericarp during the same period of storage, and 27 of these metabolites were identified (Table 3). The metabolites detected by the HPLC-QTOF-MS technology contained mainly polyphenols and flavonoids (Table 3).
HT up-regulated stress response proteins in citrus pericarp
In disease response, HT induced Beta-1, 3-glucanase (H25) and Class III chitinase (H20) accumulated in pericarp (Table 1), where these proteins are involved in reducing stress resistance in plants. In the stress response, HT up-regulated 17.7 kDa heat shock protein (H22) and low molecular weight heat-shock protein (H24) accumulated in pericarp, especially at early storage stage (Table 1). Also, other proteins accumulated in heat treated pericarp in higher amounts when compared to control pericarp during the same period of storage, which might be indirectly involved in the fruit’s response to stress, including the pyruvate dehydrogenase E1 beta subunit isoform 2 (H7), chloroplast stromal ascorbate peroxidase (H18), putative ATP synthase beta subunit (H36), putative inorganic pyrophosphatase, isoflavone reductase related protein (H10), and cytoplasmic malate dehydrogenase (H4, Table 1).
HT down-regulated reactive oxygen species metabolism-related proteins in citrus pericarp
HT down-regulated reactive oxygen species metabolism-related proteins, especially during later periods of storage (Table 1), including oxidoreductase, zinc-binding dehydrogenase family protein (H11), putative NAD(P)H oxidoreductase, isoflavone reductase (H12), aldo/keto reductase family protein (H13), copper/zinc superoxide dismutase (H16), and putative NAD(P)H oxidoreductase, isoflavone reductase (H17). Also, other proteins were down-regulated in heat treated pericarp (relative to control), such as triosephosphate isomerase, cytosolic (H5), abscisic stress ripening-like protein (H19), and thiazole biosynthetic enzyme, chloroplast precursor (H30). In addition, some proteins were up-regulated in heat treated pericarp (relative to control) only in the middle of the experiment, including, phosphoglycerate kinase, putative (H2), pyruvate dehydrogenase E1 beta subunit isoform 2 (H7), aldo/keto reductase family protein (H13), mitochondrial processing peptidase (H15), putative inorganic pyrophosphatase (H29). Putative ATP synthase beta subunit (H36) was up-regulated at the end but down-regulated at the beginning in heat treated pericarps (relative to control).
HT increased contents of sugars and flavonoids in citrus pericarp
Primary metabolite profiling of the HT and control pericarp revealed prominent changes in a group of metabolites during the entire storage phases (Table 2), mainly including sugars and fatty acids. Interestingly, the contents of many kinds of sugars were increased in heat treated pericarp (relative to control) at 2 h after HT, but decreased to the same level with control at later storage phases, such as galactose, fructose, glucose, sucrose and 4-keto-glucose. Some primary metabolites contents were increased in heat treated pericarp (relative to control) during whole storage, including oleic acid, tetradecanoic acid and ornithine (Table 2).
Secondary metabolite profiling of the HT and control pericarp revealed that a group of metabolites were increased in heat treated pericarp (relative to control) during some periods of storage, including quercetin-dihexose-deoxyhexose, hesperetin, hydroxybenzoic acid-hexose, naringenin, naringenin chalcone-hexose, vanillic acid, diosmin, rutin and isosakuranetin, which all might play a vital role in the response to stress.
HT decreased contents of organic acids, amino acids and phenylpropanoids
Primary metabolic profiling revealed organic acids, fatty acids and amino acids were down-regulated in heat treated pericarp when compared to control pericarp, especially at 2 h after heat treatment. Levels of five types of organic acids (citric acid, malic acid, 4-N-Trimethylsilylmethylaminobutyric acid, oxalic acid and GABA) had decreased in heat treated pericarp compared to the control at 2 h after HT (Table 2), but levels of ten types of organic acid increased in heat treated pericarp (relative to the control) at 32 DAT (Table 2). The content of ten types of amino acids were lower in heat treated pericarp than in the control at 2 h after HT, while only five of them were lower in heat treated pericarp than in the control at 32 DAT, although two of them were higher in heat treated pericarp than in the control at 32 DAT (Table 2). In a similar way, compared to organic acids, levels of four types of fatty acids decreased in heat treated pericarp compared to the control at 2 h after heat treatment; levels of three of them increased in heat treated pericarp compared to the control at 32 DAT (Table 2). There is common occurrence here; levels of many types of organic acids, amino acids and fatty acids decreased in heat treated pericarp compared to the control at 2 h after HT, but levels of many of them increased in heat treated pericarp compared to the control at 32 DAT. Compared to increased sugar concentrations, it appears there is a sudden transformation of organic, amino and fatty acids into sugars which is driven by HT; also, levels of many types of organic, amino and fatty acids increased in heat treated pericarp compared to the control during the later stages of storage.
Further metabolic characterization of chemicals which induced resistance in heat treated pericarp was performed using an HPLC-QTOF-MS. Contents of ferulic acid, sinapic acid, cinnamic acid and caffeic acid, the precursors of lignin synthesis, decreased in heat treated pericarp when compared to the control pericarp, especially at 2 h and 12 DAT. Moreover, the contents of other secondary metabolites decreased in heat treated pericarp, especially at 2 h or 12 DAT, such as hydroxylated naringenin-hexose, sinensetin, p-coumaroyl quinic acid, feruloylquinic acid, pantothenic acid-hexose, quercetin-O-dihexoside, eriodictyol-O-dihexoside, and jasmonic acid.
HT decreased H2O2content in citrus pericarp
The content of H2O2 decreased in heat treated pericarp during the storage phase while the content of H2O2 increased in the control pericarp (Figure 4A). This suggests HT decreased the content of H2O2 in pericarp, and the content of H2O2 was held at a lower level in heat treated pericarp when compared to control pericarp during storage.
Superoxide dismutase (SOD) and peroxidase (POD) play important roles in deoxidization of H2O2. In this study, the activity of POD showed no significant change between HT and the control (Figure 4C). The activity of SOD was higher in heat treated pericarp relative to the control at 6 DAT, but was lower in heat treated pericarp relative to the control at 12 and 32 DAT (Figure 4D). So, the activity of SOD fluctuated during the entire storage period in heat treated pericarp compared to the control (Figure 4D).
Concerning redox regulation, most of the proteins were down-regulated in heat treated pericarp when compared to control pericarp, especially in the later stages of storage (Table 1). This suggests the redox reaction was weaker in heat treated pericarp during storage.
HT increased fruit firmness and the content of lignin in citrus pericarp
Fruit hardness gradually decreased during storage (Figure 5A); however, HT increased fruit hardness, especially at 1 DAT (Figure 5A), which led to the heat treated fruit remaining firmer than the control fruits during the entire storage period. Lignin is an important component responsible for fruit hardness. So, samples from 2 h after treatment and 1, 6 and 32 DAT were used to detect lignin content changes over time. Results show HT significantly induced lignin accumulation in pericarp (Figure 5B). Also, the activity of phenylalanine ammonia-lyase (PAL), a critical enzyme in lignin synthesis, initially increased in heat treated pericarp at 1 DAT but subsequently fell into the same level as in the control (Figure 4B). Furthermore, levels of ferulic acid, sinapic acid, cinnamic acid and caffeic acid, the precursors of lignin synthesis, declined in heat treated pericarp when compared to control pericarp (Table 3). This suggests a rapid transformation occurred from ferulic acid, sinapic acid, cinnamic acid and caffeic acid to lignin in heat treated pericarp after HT.
HT improved fruit pathogen resistance without negatively impacting fruit quality
Physical treatments are used frequently to keep fruit fresh, including microporous membranes, pre-storage treatment (sweating), ultraviolet illumination, radiofrequency treatment, heat treatments (heat therapy), and other storage techniques. Postharvest heat treatments such as hot water treatment and hot air treatment provide a type of quarantine, inhibit the development of pathogens, enhance fruit resistance to chilling injury in cold-sensitive cultivars, and help preserve fruit quality during storage and extend its shelf life . In the present study, HT effectively inhibited the development of pathogens. Also, HT (52°C, 2 min) caused a statistically significant reduction the incidence of chilling injury [9, 11].Also, postharvest HT had no obvious negative effects on fruit flesh quality during subsequent storage. This suggests HT increased the ability of fruit to tolerate storage stress without affecting their commercial qualities, and shows strong potential application in the citrus industry.
HT induced stress response proteins
HT can effectively inhibit fruit chilling injury during cold storage . After HT, some proteins were up-regulated, while others were down-regulated. The up-regulated proteins might be responsible for heat induced fruit resistance. Interesting, some proteins were up-regulated at 1 d after HT, and decreased to the previous level of control during the storage period. Those kinds of proteins might be just involved in response to HT, but not play a vital role in fruit disease resistance. The proteins, which are up-regulated during whole storage period might be directly involved in fruit disease resistance.
In the present study, heat shock protein (HSP) 17.7, putative HSP 70 and low molecular weight HSP were up-regulated in heat treated pericarp compared to the control during entire storage period (Table 1). HSPs are encoded by a multigene family and are located in different subcellular compartments . These molecular chaperones play a crucial role in protecting plants against stress by re-establishing normal protein configuration and thus helping to maintain cellular homeostasis . Abiotic stresses usually causes protein dysfunction; so maintaining proteins in their functional configurations and preventing the aggregation of non-native proteins is particularly important for cell survival under stress . HSPs are responsible for protein folding, assembly, translocation and degradation during cellular processes. HSPs also stabilize proteins and membranes, and can assist in protein refolding under stressful conditions . Our results suggest HSPs play a crucial role in reducing the susceptibility of HT fruit to biotic and abiotic stress.
HT significantly inhibited pathogen development and fruit decay during storage . Chitinase and β-1, 3-glucanase were up-regulated in heat treated fruit in this study; these enzymes have long been thought to be active in fruit antifungal defences. Pathogenesis-related (PR) proteins have been classified into five groups and two of these groups contain proteins with β-1, 3-glucanase and chitinase . Chitinase catalyses the hydrolysis of chitin, a β-1, 4-linked polymer of N-acetyl D-glucosamine and is a major component of cell walls of most phytopathogenic fungi . Crude protein extracts of chitinase and β-1, 3-glucanase were observed to effectively inhibit fungal activity . Also, β-1, 3-glucanase and chitinase are thought to contribute to the level of resistance provided and are known to inhibit mycelial growth of a wide range of fungal pathogens . This suggests β-1, 3-glucanase and chitinase were directly involved in HT induced pathogen resistance in fruits.
HT increased the content of stress response metabolites
Primary metabolites are directly involved in normal growth, development, and reproduction. In the present study, the content of ornithine, oleic acid and tetradecanoic acid increased in heat treated pericarp compared to the control during storage, and levels of seven types of sugars increased in heat treated pericarp compared to the control at 2 h after HT (Table 2). In plants, ornithine is required for the synthesis of polyamines and alkaloids, which contribute to oxidative stress tolerance in plants subjected to severe water stress .
Flavonoids are important defensive compounds which protect fruit against the development of pathogens, acting either as phytoanticipins or phytoalexins . Flavonoids were observed to be up-regulated in heat treated pericarp (compared to the control), including quercetin-dihexose-deoxyhexose, hesperetin, naringenin, naringenin chalcone-hexose, hydroxybenzoic acid-hexose, isosakuranetin, diosmin, and rutin, which may play a vital role in response to stress . Rutin is a bioflavonoid with strong antioxidant activity, which had an antimicrobial effect on pathogens, as evidenced by reduced conidial germination and appressorium formation of the pathogen . Quercetin protected cells from oxidative stress-induced cell death by blocking both the cyclooxygenase and lipooxygenase pathways . Quercetin has a protective function against oxidative stress by reducing oxidative injury to the cells during H2O2 treatment . Also, plant extracts rich in quercetin have been applied as an interesting resource for functional food products and other foods used to reduce the risks of age-related macular degeneration . Hesperidin acts as an antioxidant, based on in vitro studies, and hesperidin also provides strong cellular antioxidant protection against the damaging effects induced by stress . Diosmin, hesperidin and naringin are flavonoid glycosides that occur naturally in citrus fruits; they exert a variety of properties such as antioxidant and free radical scavenging . Interestingly, naringenin and rutin are not only involved in response to stress, but also were assumed to be signalling molecules .
HT induced fruit pathogen resistance and increased lignin content
The functions of phenylpropanoid compounds in plant defence range from preformed or inducible physical and chemical barriers against infection . Defensive functions are not restricted to a particular class of phenylpropanoid compounds, but are found in the lignification of monolignols to lignin . In the present study, HT decreased the level of precursors needed in lignin synthesis (ferulic acid, sinapic acid, cinnamic acid and caffeic acid), but increased lignin content (Table 3, Figure 5). Lignin is one of the major components of secondary cell walls, providing cells with mechanical support and isolation from the environment. In plants, cells can sense cell wall polymer structure and this activates trans-plasma membrane sensor proteins which might trigger a cascade of signal transduction events interwoven with other pathways that might modulate cellular functions by the activation or inhibition of specific transcription factors as well as by affecting posttranslational control of gene expression and protein function . If a cell wall is damaged, cells activate a series of responses to the invading disease organisms, which might negatively impact fruit quality. The accumulation of lignin induced by HT creates a thickening of cell walls which forms an effective physical barrier to pathogens, delaying the invasion of disease organisms. Also, HT increased fruit firmness; this might be a result of HT induced lignin thickened fruit cell walls.
HT reduced the content of H2O2
HT is an effective method for improving citrus storability. In storage, fruits suffer a variety of biotic and abiotic environmental stresses. All these stresses could lead to the generation of reactive oxygen species (ROS). To avoid the accumulation of ROS, which may lead to cell death, plants have up-regulated enzymatic and non-enzymatic antioxidants to scavenge ROS and alleviate their negative effects . In this study, HT down-regulated H2O2 content in pericarp, which was consistent with the effect of HT on yeast cells. HT yeast cells accumulated less ROS than untreated cells in response to stresses .
Also, isoflavone reductase (IFR) might have an important role as an antioxidant . IFR is an enzyme involved in the production of isoflavone phytoalexins and it accumulates in response to pathogenic attack, fungal elicitor and biotic stress . Isoflavonoids are important secondary metabolites for essential physiological processes. Kim et al. reported the isoflavone reductase-like gene may act as a down-regulator by producing antioxidant chemicals to prevent the elevation of ROS. In this study, isoflavone reductase, a related protein (H10), was up-regulated in heat treated pericarp when compared to control pericarp (Table 1), which might be another factor leading to lower H2O2 levels and stronger resistance to pathogens in heat treated pericarp when compared to control pericarp.
Ascorbate peroxidase scavenges H2O2 in specific organelles of the cell . Also, Cu/Zn-SOD is the main cellular isoform of SODs which play a key role in the antioxidant defence system through the disproportionation of H2O2. Both of these can prevent ROS caused damages to the cellular membrane and act as a primary defence mechanism when organisms were exposed to oxidative stress. These proteins are considered the main enzymatic and non-enzymatic systems for protecting cells against oxidative damage, and are responsible for disease resistance. In response to biotic stress, these proteins also act as the primary defence mechanism against oxidative stresses caused by pathogens, therefore preventing damage of ROS to cellular membranes. However, the activities of POD and SOD were not elevated in Heat treated pericarp when compared to control pericarp, which might be caused by the lower H2O2 content in heat treated pericarp. Also, other proteins involved in scavenging ROS were down-regulated in heat treated pericarp when compared to control pericarp, including oxidoreductase (H11), oxidoreductase (H12, 17), aldo/keto reductase (H13), 2-oxoacid dehydrogenase (H14), superoxide dismutase (H16). This suggests HT induced a level of resistance to pathogens attacking fruits by causing a decline in fruit H2O2 content and by down-regulated scavenging ROS enzymes.
Role of sugars, organic acids and amino acids in improving storability
In this study, the profiling of primary metabolics showed HT induced the accumulation of sugars at early stages of storage (Table 2). Soluble sugars, especially sucrose, glucose, and fructose, accumulated when fruits were subjected to stress, and accumulated soluble sugars play an obviously vital role in the response of fruit to a number of stressors . Soluble sugars do not only act as nutrients and metabolites, but also act as signal substance in the modification of redox homeostasis . However, less information was available on how other sugars are involved in the stress response, especially arabinose, mannose, galactose, and ribose.
Moreover, primary metabolic data shows HT reduced organic acids and amino acids content especially at 2 h after treatment in this study. Although organic acids and amino acids play an obvious central role in energy and metabolism at the cellular and whole-organism levels, there seems to be a rapid conversion from organic acids and amino acids to sugars in heat treated pericarp. However, little is known about the decreased levels of organic acids and amino acids seen during a stress response, and further study will be required to validate this response.
This study provides a broad picture of differential accumulation of proteins and metabolites, and gives a new insight into how HT increases fruit stress resistance during subsequent storage of citrus fruit. In metabolic substrates, the up-accumulation of secondary metabolites played vital roles in increasing the ability of fruit to deal with stress. Flavonoids are directly involved in the response to external stress. The increased lignin and decreased H2O2 contents involved in HT probably increased fruit resistibility in response to external stress.
‘Kamei’ Satsuma mandarin (Citrus unshiu Marc.) fruits were harvested at commercial maturity (fruit colour has turned completely orange and the level of soluble solids approaches 9%) from an orchard in the city of Yichang, Hubei Province, China, in October 2009. Fruits of a uniform size and colour, free of visible injury or blemishes, were selected for the experiments. Three hundred kilograms of fruit were divided into two equal parts, 150 kg were used for the HT and the other 150 kg used for a control. After treatments, fruits were stored in a ventilated warehouse (storage temperature: 12–16°C; relative humidity: 90–95%). Samples were collected at 2 h, and 1, 2, 3, 4, 6, 9, 12, 16, 20, 24, 28, 32, 38, 44, 50 and 60 DAT. The pericarp was sampled along the equatorial plane of each fruit; 30 fruits were selected as one sample and ground to powder in liquid nitrogen, then stored at –80°C for further analysis.
Heat treatment (HT)
Heat treated fruits were dipped in a 52°C warm water bath for 2 min, and control fruits were dipped in a 25°C water bath for 2 min. After treatments, fruits were air dried and stored in a ventilated warehouse prior to sampling.
Fruit quality determination
Fruits were checked for weight loss (%), respiration rate (mg kg–1 h–1) and total soluble solids (TSS,%) at 1, 2, 3, 4, 6, 9, 12, 16, 20, 24, 28, 32, 38, 44, 50 and 60 DAT. The weight was measured using an electronic hydrostatic balance (Model: MP31001, Shanghai Selon Scientific Instrument, Co., Ltd, Shanghai, China) with an accuracy of ± 0.01 g. The respiration rate of fruit was measured with an infrared gas analyser (Model: GXH-305H, Junfang Science & Technology Institute of Physical and Chemical Research, Beijing, China). TSS were determined with a refractometer (Model: Pocket PAL-1, Atago Inc., Toyko, Japan) according to the manufacturer’s instructions.
Antifungal assay of HT in wounded citrus fruits
Blue mould (Penicillium italicum) was chosen as an indicator to evaluate the resistance of treated fruit to fungal infections. 150 fruits from each treatment were used for fungal inoculation. A uniform lesion (3 mm deep, 4 mm wide) was made at the equator of the fruit using a sterile nail. Aliquots of 10 μl suspension of P. italicum at 1×105 spore ml–1 were inoculated into each wound site. After fungal inoculation, fruits were stored in a storage chamber (relative humidity, 95%; temperature, 25°C), and disease incidence was detected at 1, 2, 3, 4, 5, 6 and 7 d after inoculation. Disease incidence rates and lesion diameters were recorded based on the following equations:
Two-dimensional gel electrophoresis and MAILDI-TOF/TOF
Samples at 1, 6, 12 and 32 DAT were used for total protein extraction based on the method described by Isaacson et al.. The precipitate was then nitrogen gas-dried and solubilized in the lysis buffer containing 8 M urea, 0.2% (w/v) Bio-Lyte, 4% CHAPS, 65 mM DTT. The protein concentration was determined with a Bio-Rad protein assay kit (Bio-Rad Laboratories Inc., Hercules, CA, USA) based on the Bradford method using BSA as the standard. Two-dimensional gel electrophoresis and gel staining was carried out based on the method of Yun et al.. The stained gels were imaged with a UVP imager (Bio-Rad) in the ‘trans white’ mode using PDQuest 2-D analysis software version 7.4 (Bio-Rad). Protein spots matching and differential protein spots analysis were completed based on the method of Yun et al.. Isoelectric points and molecular weights of the proteins were determined by comparison with the markers.
In-gel digestion and protein identification were performed using a 4800 Proteomics Analyzer MALDI TOF/TOF (Applied Biosystems, Foster City, CA) as described by Yun et al.. MS/MS data were then submitted to the MASCOT program for protein identification against the green plants database. The following are the search parameters: trypsin enzyme, one missed cleavage, fixed modifications of carbamidomethyl (C), fragment mass tolerance of ± 0.5 Da, variable modifications of oxidation (Met), peptide tolerance of 100 ppm, and peptide charge of 1+. Only results of peptides with MS/MS P < 0.05 confidence were accepted.
The primary metabolic profiling
Samples collected at 2 h and 1, 6 and 32 DAT were used for the differential primary metabolic profiling analysis. A 300 mg sample was extracted in 2,700 μl of methanol as described by Roessner-Tunali et al.. 300 μl of 0.2 mg ml–1 ribitol in water was added as a quantification internal standard. A derivatization reaction was performed based on the protocol of Zhang et al. with little modification. Extracts were incubated in 50 μl of 20 mg ml–1 methoxyamine hydrochloride in pyridine for 30 min at 50°C, followed by a 40 min treatment at 60°C using 50 μl BSTFA (containing 1% TMCS). Each sample (1 μl) was injected into the gas chromatograph system through a fused-silica capillary column (30 m × 0.25 mm i.d., 0.25 μm) DB-5 MS stationary phase. The injector temperature was 250°C, with a carrier gas flow rate of 1.0 ml min–1. The column temperature was held at 100°C for 1 min; increased to 184°C at a rate of 3°C min–1, increased to 190°C at a rate of 0.5°C min–1, increased to 280°C at 15°C min–1. The flow rate of the carrier helium (99.999%) gas was 1 ml min–1. The following were MS operating parameters: ionization voltage, 70 eV (electron impact ionization); ion source temperature, 200°C; interface temperature, 250°C. TIC (total ioncurrent) spectra were recorded in the mass range of 45–600 atomic mass units in scanning mode.
The secondary metabolic profiling
Samples at 2 h and 1, 6 and 32 DAT were used for differential secondary metabolic profiling analysis using HPLC-MS. 300 mg dried powder was extracted with methanol, the mixture was filtered through a 0.22–mm polytetrafluoroethylene membrane filter.
Metabolite profiling was performed using a QTOF 6520 mass spectrometer (Agilent Technologies, Palo Alto, CA, USA) coupled to a 1200 series Rapid Resolution HPLC system as described by Page et al.. 2 μl of sample extract was loaded onto a Zorbax Eclipse Plus C18 1.8 μm, 2.1 × 100 mm reverse-phase analytical column (Agilent Technologies). Mobile phase A comprised 0.1% formic acid in water and mobile phase B was acetonitrile. The following gradient was used: 0 min–10% B; 20 min–95% B; 22 min–95% B; 22.1 min–10% B; 30 min–10% B. The flow rate was 0.3 ml min–1 and the column temperature was held at 35°C for the duration. The source conditions for electrospray ionization were as follows: gas temperature was 350°C with a drying gas flow rate of 10 l min–1 and a nebulizer pressure of 40 psig. The capillary voltage was 3.5 kV in both positive and negative ion mode. The fragmentor voltage was 135 V and the skimmer was 65 V.
Assay of phenylalanine ammonia-lyase (PAL), superoxide dismutase (SOD) and peroxidase (POD) activities
Samples collected at 1, 6, 12 and 32 DAT were used for enzyme activity determination. A 5 g sample was homogenized in 25 ml of 0.05 M sodium borate buffer (pH 8.8 for PAL, containing 5 mM β-mercaptoethanol) and 25 ml of ice-cold 100 mM sodium phosphate (pH 7.8 for POD and SOD, containing 0.5 g polyvinyl polypyrrolidone). PAL activity was determined by measuring the absorbance of cinnamic acid at 290 nm over a period of 30 min at 30°C . POD activity was analysed using guaiacol as the substrate . SOD activity was determined using the xanthine-xanthine oxidase method, and nitro blue tetrazolium was used as the indicator of superoxide radical production .
Determination the content of H2O2in the peel
Samples collected at 2 h and 1, 6 and 32 DAT were used to measure hydrogen peroxide levels based on the methods of Patrick and Wagner  with some modifications. Briefly, peel tissues (0.5 g) were ground in liquid nitrogen and homogenized in 5 ml of physiological saline. Samples were then subjected to an ultrasonic bath for 15 min and centrifuged at 6,000 g for 20 min. The supernatants (50 μl) were subsequently added into 0.5 ml of solution I and warmed in a 37°C water bath for 10 min before 0.5 ml of solution II was added. The absorbance of the mixture was measured at 405 nm and the content of H2O2 in each sample was calculated by comparison with a standard calibration curve.
Determination of fruit firmness
Samples collected at 2 h and 1, 6 and 32 DAT were used for the determination of fruit firmness. Eighteen fruits were measured and every fruit was measured on three different spots around the equatorial plane. Firmness was measured with a fruit sclerometer (GY-B, Jilin, China) with 4 mm diameter and expressed as N cm–2 and the average value per fruit was used for further analysis, after removing three largest and three smallest measurements; a total twelve replicates were comparative statistically analysed per sample using Student’s t-test (P < 0.05).
Determination of lignin in pericarp
Samples collected at 2 h and 1, 6 and 32 DAT were used to measure fruit lignin content. Lignin levels were determined using the method described by Lei et al..
The experimental design was completely randomized with more than three replications. Four individual replicates were completed in 2-DE, GC-MS, HPLC-MS, H2O2 content, lignin content and enzyme activity determination. Twelve replicates were completed to measure weight loss, respiration rate, total soluble solids and fruit firmness. Data for each sample was statistically analysed using Student’s t-test (P < 0.05).
Two-dimensional gel electrophoresis
Gas chromatography coupled to mass spectrometry
Heat shock protein
Liquid chromatography quadrupole time-of-flight mass spectrometry
- MALDI-TOF MS:
Matrix-assisted laser desorption/ionization-time-of-flight tandem mass spectrometry
National center for biotechnology information
Reactive oxygen species
Total soluble solids.
Gonzalez-Candelas L, Alamar S, Sanchez-Torres P, Zacarias L, Marcos JF: A transcriptomic approach highlights induction of secondary metabolism in citrus fruit in response to Penicillium digitatum infection. BMC Plant Biol. 2010, 10: 194. 10.1186/1471-2229-10-194.
Tu K, Shao XF, Tu SC, Su J, Zhao Y: Effects of heat treatment on wound healing in gala and Red fuji apple fruits. J Agr Food Chem. 2010, 58: 4303-4309. 10.1021/jf904273m.
Lara MV, Borsani J, Budde CO, Lauxmann MA, Lombardo VA, Murray R, Andreo CS, Drincovich MF: Biochemical and proteomic analysis of ‘Dixiland’ peach fruit (Prunus persica) upon heat treatment. J Exp Bot. 2009, 60: 4315-4333. 10.1093/jxb/erp267.
Gould GW: History of science - spores - lewis B perry memorial lecture 2005. J Appl Microbiol. 2006, 101: 507-513. 10.1111/j.1365-2672.2006.02888.x.
Hong SI, Lee HH, Kim D: Effects of hot water treatment on the storage stability of satsuma mandarin as a postharvest decay control. Postharvest Biol Tec. 2007, 43: 271-279. 10.1016/j.postharvbio.2006.09.008.
Zhang ZK, Bi Y, Ge YH, Wang JJ, Deng JJ, Xie DF, Wang Y: Multiple pre-harvest treatments with acibenzolar-S-methyl reduce latent infection and induce resistance in muskmelon fruit. Sci Hortic-Amsterdam. 2011, 130: 126-132. 10.1016/j.scienta.2011.06.024.
Schirra M, D’Aquino S, Cabras P, Angioni A: Control of postharvest diseases of fruit by heat and fungicides: efficacy, residue levels, and residue persistence. A Review. J Agr Food Chem. 2011, 59: 8531-8542. 10.1021/jf201899t.
Perotti VE, Del Vecchio HA, Sansevich A, Meier G, Bello F, Cocco M, Garrán SM, Anderson C, Vázquez D, Podestá FE: Proteomic, metabalomic, and biochemical analysis of heat treated Valencia oranges during storage. Postharvest Biol Tec. 2011, 62: 97-114. 10.1016/j.postharvbio.2011.05.015.
Sapitnitskaya M, Maul P, McCollum GT, Guy CL, Weiss B, Samach A, Porat R: Postharvest heat and conditioning treatments activate different molecular responses and reduce chilling injuries in grapefruit. J Exp Bot. 2006, 57: 2943-2953. 10.1093/jxb/erl055.
Dixon RA, Achnine L, Kota P, Liu CJ, Reddy MSS, Wang LJ: The phenylpropanoid pathway and plant defence - a genomics perspective. Mol Plant Pathol. 2002, 3: 371-390. 10.1046/j.1364-3703.2002.00131.x.
Ghasernnezhad M, Marsh K, Shilton R, Babalar M, Woolf A: Effect of hot water treatments on chilling injury and heat damage in ‘satsuma’ mandarins: Antioxidant enzymes and vacuolar ATPase, and pyrophosphatase. Postharvest Biol Tec. 2008, 48: 364-371. 10.1016/j.postharvbio.2007.09.014.
Vigneault C, Lu J, Charles MT, Goyette B, Raghavan GSV: Effect of heat treatment uniformity on tomato ripening and chilling injury. Postharvest Biol Tec. 2010, 56: 155-162. 10.1016/j.postharvbio.2010.01.005.
Vierling E: The roles of heat shock proteins in plants. Annu Rev Plant Biol. 1991, 42: 579-620. 10.1146/annurev.pp.42.060191.003051.
Wang WX, Vinocur B, Shoseyov O, Altman A: Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9: 244-252. 10.1016/j.tplants.2004.03.006.
Basha E, Jones C, Wysocki V, Vierling E: Mechanistic differences between two conserved classes of small heat shock proteins found in the plant cytosol. J Biol Chem. 2010, 285: 11489-11497. 10.1074/jbc.M109.074088.
Hu WZ, Jiang AL, Jin LM, Liu CH, Tian MX, Wang YY: Effect of heat treatment on quality, thermal and pasting properties of sweet potato starch during yearlong storage. J Sci Food Agr. 2011, 91: 1499-1504. 10.1002/jsfa.4340.
Schraudner M, Ernst D, Langebartels C, Sandermann H: Biochemical plant responses to ozone: III. Activation of the defence-related proteins β-1, 3-glucanase and chitinase in tobacco leaves. Plant Physiol. 1992, 99: 1321. 10.1104/pp.99.4.1321.
Chalavi V, Tabaeizadeh Z, Thibodeau P: Enhanced resistance to Verticillium dahliae in transgenic strawberry plants expressing a Lycopersicon chilense chitinase gene. J Am Soc Hortic Sci. 2003, 128: 747-753.
Mauch F, Mauch-Mani B, Boller T: Antifungal hydrolases in pea tissue: II. Inhibition of fungal growth by combinations of chitinase and β-1, 3-glucanase. Plant Physiol. 1988, 88: 936-942. 10.1104/pp.88.3.936.
Magnin-Robert M, Trotel-Aziz P, Quantinet D, Biagianti S, Aziz A: Biological control of Botrytis cinerea by selected grapevine-associated bacteria and stimulation of chitinase and beta-1,3 glucanase activities under field conditions. Eur J Plant Pathol. 2007, 118: 43-57. 10.1007/s10658-007-9111-2.
Kalamaki MS, Alexandrou D, Lazari D, Merkouropoulos G, Fotopoulos V, Pateraki I, Aggelis A, Carrillo-Lopez A, Rubio-Cabetas MJ, Kanellis AK: Over-expression of a tomato N-acetyl-L-glutamate synthase gene (SlNAGS1) in Arabidopsis thaliana results in high ornithine levels and increased tolerance in salt and drought stresses. J Exp Bot. 2009, 60: 1859-1871. 10.1093/jxb/erp072.
Radhakrishna Shetty R, Jensen JD: Silicon-induced changes in antifungal phenolic acids, flavonoids and Key phenylpropanoid pathway genes during the interaction between miniature roses and the biotrophic pathogen podosphaera pannosa. Plant Physiol. 2011, 157: 2194-2205. 10.1104/pp.111.185215.
Treutter D: Significance of flavonoids in plant resistance: a review. Environ Chem Lett. 2006, 4: 147-157. 10.1007/s10311-006-0068-8.
Nijveldt RJ, van Nood E, van Hoorn DEC, Boelens PG, van Norren K, van Leeuwen PAM: Flavonoids: a review of probable mechanisms of action and potential applications. Am J Clin Nutr. 2001, 74: 418-425.
Saviranta NMM, Veeroos L, Granlund LJ, Hassinen VH, Kaarniranta K, Karjalainen RO: Plant flavonol quercetin and isoflavone biochanin A differentially induce protection against oxidative stress and inflammation in ARPE-19 cells. Food Res Int. 2011, 44: 109-113. 10.1016/j.foodres.2010.10.056.
Wilmsen PK, Spada DS, Salvador M: Antioxidant activity of the flavonoid hesperidin in chemical and biological systems. J Agr Food Chem. 2005, 53: 4757-4761. 10.1021/jf0502000.
Kanaze FI, Gabrieli C, Kokkalou E, Georgarakis M, Niopas I: Simultaneous reversed-phase high-performance liquid chromatographic method for the determination of diosmin, hesperidin and naringin in different citrus fruit juices and pharmaceutical formulations. J Pharmaceut Biomed. 2003, 33: 243-249. 10.1016/S0731-7085(03)00289-9.
Denness L, McKenna JF, Segonzac C, Wormit A, Madhou P, Bennett M, Mansfield J, Zipfel C, Hamann T: Cell wall damage-induced lignin biosynthesis is regulated by a reactive oxygen species- and jasmonic acid-dependent process in arabidopsis. Plant Physiol. 2011, 156: 1364-1374. 10.1104/pp.111.175737.
Seifert GJ, Blaukopf C: Irritable walls: the plant extracellular matrix and signaling. Plant Physiol. 2010, 153: 467-478. 10.1104/pp.110.153940.
Apel K, Hirt H: Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol. 2004, 55: 373-399. 10.1146/annurev.arplant.55.031903.141701.
Liu J, Wisniewski M, Droby S, Tian SP, Hershkovitz V, Tworkoski T: Effect of heat shock treatment on stress tolerance and biocontrol efficacy of Metschnikowia fructicola. FEMS Microbiol Ecol. 2011, 76: 145-155. 10.1111/j.1574-6941.2010.01037.x.
Kim SG, Kim ST, Wang YM, Kim SK, Lee CH, Kim KK, Kim JK, Lee SY, Kang KY: Overexpression of rice isoflavone reductase-like gene (OsIRL) confers tolerance to reactive oxygen species. Physiol Plantarum. 2010, 138: 1-9. 10.1111/j.1399-3054.2009.01290.x.
Kim ST, Cho KS, Yu S, Kim SG, Hong JC, Han CD, Bae DW, Myung AE, Kang KY: Proteomic analysis of differentially expressed proteins induced by rice blast fungus and elicitor in suspension-cultured rice cells. Proteomics. 2003, 3: 2368-2378. 10.1002/pmic.200300577.
Zhang HX, Xia Y, Wang GP, Shen ZG: Excess copper induces accumulation of hydrogen peroxide and increases lipid peroxidation and total activity of copper-zinc superoxide dismutase in roots of Elsholtzia haichowensis. Planta. 2008, 227: 465-475.
Rosa M, Prado C, Podazza G, Interdonato R, González JA, Hilal M, Prado FE: Soluble sugars—metabolism, sensing and abiotic stress: a complex network in the life of plants. Plant Signal Behav. 2009, 4: 388. 10.4161/psb.4.5.8294.
Rosenwasser S, Rot I, Sollner E, Meyer AJ, Smith Y, Leviatan N, Fluhr R, Friedman H: Organelles contribute differentially to reactive oxygen species-related events during extended darkness. Plant Physiol. 2011, 156: 185. 10.1104/pp.110.169797.
Isaacson T, Damasceno CMB, Saravanan RS, He Y, Catala C, Saladie M, Rose JKC: Sample extraction techniques for enhanced proteomic analysis of plant tissues. Nat Protoc. 2006, 1: 769-774. 10.1038/nprot.2006.102.
Yun Z, Li WY, Pan ZY, Xu J, Cheng YJ, Deng XX: Comparative proteomics analysis of differentially accumulated proteins in juice sacs of ponkan (Citrus reticulata) fruit during postharvest cold storage. Postharvest Biology and Technology. 2010, 56: 189-201. 10.1016/j.postharvbio.2010.01.002.
Roessner-Tunali U, Hegemann B, Lytovchenko A, Carrari F, Bruedigam C, Granot D, Fernie AR: Metabolic profiling of transgenic tomato plants overexpressing hexokinase reveals that the influence of hexose phosphorylation diminishes during fruit development. Plant Physiol. 2003, 133: 84-99. 10.1104/pp.103.023572.
Zhang JJ, Wang X, Yu O, Tang JJ, Gu XG, Wan XC, Fang CB: Metabolic profiling of strawberry (Fragariaxananassa Duch.) during fruit development and maturation. J Exp Bot. 2011, 62: 1103-1118. 10.1093/jxb/erq343.
Page M, Sultana N, Paszkiewicz K, Florance H, Smirnoff N: The influence of ascorbate on anthocyanin accumulation during high light acclimation in Arabidopsis thaliana: further evidence for redox control of anthocyanin synthesis. Cell & Environment: Plant; 2011.
Ballester AR, Lafuente MT, Gonzalez-Candelas L: Spatial study of antioxidant enzymes, peroxidase and phenylalanine ammonia-lyase in the citrus fruit-Penicillium digitatum interaction. Postharvest Biol Tec. 2006, 39: 115-124. 10.1016/j.postharvbio.2005.10.002.
Meng XH, Han J, Wang Q, Tian SP: Changes in physiology and quality of peach fruits treated by methyl jasmonate under low temperature stress. Food Chem. 2009, 114: 1028-1035. 10.1016/j.foodchem.2008.09.109.
Sala JM, Lafuente MT: Antioxidant enzymes activities and rindstaining in ‘Navelina’ oranges as affected by storage relative humidity and ethylene conditioning. Postharvest Biol Tec. 2004, 31: 277-285. 10.1016/j.postharvbio.2003.10.002.
Patrick WA, Wagner HB: Determination of hydrogen peroxide in small concentrations. Anal Chem. 1949, 21: 1279-1280. 10.1021/ac60034a038.
Lei Y, Liu Y-Z, Gu Q-Q, Yang X-Y, Deng X-X, Chen J-Y: Comparison of cell wall metabolism in the pulp of three cultivars of ‘Nanfeng’ tangerine differing in mastication trait. J Sci Food Agr. 2011, 92: 496-502.
We thank Prof. Hanhui Kuang of Huazhong Agricultural University for help in revising the manuscript. This work was supported by Huazhong Agricultural University Scientific & Technological Self-innovation Foundation, the Natural Science Foundation of China (NSFC No. 31271968), The Important Project of Ministry of Education (311029), and National Modern Agriculture (Citrus) Technology Systems of China (No. CARS-27).
The authors declare that they have no competing interests.
ZY, HJG and YJC conceived the study, designed all the experiments and performed silicon and biochemical analyses. SZL carried out primary metabolite detection and statistical analysis. SJ carried out secondary metabolite detection and statistical analysis. ZY, HJG, TL, QX, JX, YJC and XXD interpreted the experimental dada and participated in writing the manuscript. All the authors read and approved the final manuscript.
Ze Yun, Huijun Gao contributed equally to this work.