ORIGINAL RESEARCH article

Front. Nutr., 01 December 2022

Sec. Nutrition and Food Science Technology

Volume 9 - 2022 | https://doi.org/10.3389/fnut.2022.1066043

Transcriptomic and physiological analysis reveals the possible mechanism of ultrasound inhibiting strawberry (Fragaria × ananassa Duch.) postharvest softening

  • Beijing Advanced Innovation Center for Food Nutrition and Human Health (BTBU), School of Food and Health, Beijing Higher Institution Engineering Research Center of Food Additives and Ingredients, Beijing Technology and Business University, Beijing, China

Abstract

Ultrasound effectively inhibited strawberry softening but the mechanism was not clear. In this study, physical data including firmness, soluble pectin (SP) contents, pectin esterase (PE), polygalacturonase (PG) activity and transcriptome sequencing data were analyzed to explore the mechanism of strawberry response to ultrasonic treatment. After 24 days storage, the firmness reduction rate and soluble contents (SP) increased rate of the strawberry treated with ultrasound (25 kHz, 0.15 W/cm2) for 3 min decreased 41.70 and 63.12% compared with the control, respectively. While the PG and PE enzyme activities of ultrasound-treated strawberries were significantly lower than control after storage for 18 days. A total of 1,905 diferentially expressed genes (DEGs) were identified between ultrasound-treated and control, with 714 genes upregulated and 1,254 genes downregulated, including 56 genes in reactive oxygen species (ROS), auxin (AUX), ethylene (ETH) and jasmonic acid (JA) signaling pathways. At 0 h, 15 genes including LOX, JMT, ARP, SKP, SAUR, IAA, ARF, and LAX were significantly upregulated compared with the control group, which means reactive oxygen specie, auxin, ethylene and jasmonic acid-mediated signaling pathway respond to ultrasound immediately. ERF109, ERF110, and ACS1_2_6 downregulated before 2 days storage indicated ethylene signaling pathway was inhibited, while after 2 days, 9 genes including ERF027, ERF109, and ERF110 were significantly upregulated indicating that the response of the ethylene signaling pathway was lagging. Therefore, in strawberry ultrasound enhanced ROS scavenging and activated JA biosynthesis, which acts as a signal for delaying the activation of ET signaling pathway, thus suppressing the activity of pectin-degrading enzymes PE and PG, and ultimately inhibiting postharvest softening.

Introduction

Strawberry (Fragaria × ananassa Duch.) is widely appreciated due to its fragrance, deliciousness and read high nutritional value, which is a rich source of vitamins, minerals, amino acids, and biologically active ingredients such as phenolic compounds, flavonoids, and anthocyanins (). However, postharvest strawberries are highly perishable during transport and storage, because they have thin skin and juicy pulp, thus quite vulnerable to mechanical injury, as well as to fungal pathogens spoilage (). As a pulp fruit, the firmness is one of the most essential quality attributes, which is in relation to fruit quality, resistance to postharvest diseases, shelf life, transportation capacity, and especially the acceptability of consumers ().

Fruit softening occurs naturally during ripening, primarily as a consequence of polysaccharides (pectin, hemicellulose and cellulose) metabolism of the cell wall and middle lamella (). Numerous researches have concluded that cellulose-hemicellulose network crosslinked with pectin, via non-covalent interactions, contributes to cell wall strength and tissue integrity (, , ). Among the three carbohydrate polymers, pectin consists of 1,4-linked α-D-galacturonic acid residues predominantly, mostly present in the homogalacturonan (HG) and rhamnogalacturonan I (RG-I) (), hemicelluloses are composed of xylans, xyloglucans as well as mannans, and cellulose is a homopolymer of β-(1,4)-linked glucose residues (, ). The disassembly of the cell wall structure requires the orderly participation of multiple cell wall modifying enzymes and proteins. In strawberries, softening is characterized by increasing pectin solubilization, i.e., an increase in the quantity of pectin loosely bound to the cell wall (). Pectin degrading enzymes, including pectinesterase (PE, EC 3.1.1.11), polygalacturonase (PG, EC 3.2.1.15) and pectate lyase (PL, EC 4.2.2.2), display different activities in pectin metabolism (, ). PE is a hydrolase that demethoxylates homogalacturonans in pectin, provide substrates for PG to hydrolyze demethylated polygalacturonic acid to galacturonic acid, resulting in the increase of soluble pectin (SP)content (, ), and PL can furtherly cleave β(1-4) linkages between galacturonosyl residues by β-elimination (, ). In metabolism of other cell wall polysaccharides, β-galactosidase (β-Gal, EC 3.2.1.23) removes galactosyl residues in side chains of cell wall polysaccharides (pectin and xyloglucan), cellulase (Cx, EC 3.2.1.4) can degrade hemicellulose and cellulose in cell wall (), β-xylosidase (β-Xyl) (EC 3.2.1.37) hydrolizes xylooligosaccharides releasing xylose units (). Besides, endo-1,4-beta-glucanase (EGase, EC 3.2.1.4), xyloglucan endo-transglucosylase/hydrolase (XTH, EC 2.4.1.207, and EC 3.2.1.151) participate in hemicellulose metabolism helping cell wall disassembly (, ). In addition, expansins (EXPs) may play an indirect role in pectin degradation. These EXP proteins without catalytic activity can promote non-enzymatic cell wall loosening by binding cellulose and hemicellulose and disrupting non-covalent interactions ().

As mentioned above, fruit postharvest softening is a complex event that requires the orchestrated participation of a wide variety of cell wall degrading proteins. In this process, multiple signals such as reactive oxygen species (ROS), phytohormones including ethylene (ET), abscisic acid (ABA), jasmonic acid (JA), auxin, brassinosteriods and cytokinins, etc., and a set of transcription factors (TFs) are the organizers (, , ). Low level of ROS can serve as a signal molecule in response to abiotic stress in the biological process of plants (, ), and hydrogen peroxide (H2O2) is an important form of ROS. In recent years, there is increasing evidence that ABA plays an important role in regulating fruit ripening process and participating in softening by promoting pectin degradation, including in strawberry (). ET is considered another important regulator in fruit ripening and softening, accomplished by the regulation of a series of ethylene signal transductions, as well as the effect on gene expression of cell wall-degradation enzymes (). The results of Villarreal et al. () confirmed that ethylene participates in the regulation of strawberry fruit softening. JA, considered as a new type of phytohormone, is reported to regulate the coloring, softening and hormone-related gene expression, thus promoting fruit softening and ripening in strawberry (). Moreover, recent studies have indicated the involvements of auxins, brassinosteriods and cytokinins in fruit ripening and senescence (). In addition, TFs also play essential roles in gene regulation during fruit ripening ().

Power ultrasound (20–100 kHz), propagating in a medium in the form of mechanical vibration (, ), has been evaluated as an innovative nonthermal preservation technology for fruits and vegetables (). Ultrasound technology is considered to be attractive, because it can offer significant advantages over other preservation technologies, such as high reproducibility, low cost, ease of treatment, environmentally friendly, and the removal of residues (, ). Previous studies have shown that ultrasound treatment can maintain quality characteristics including texture, color, flavor and antioxidant compounds, inactivate microorganisms, influence enzyme activity in fruits and vegetables (). Cavitation effect caused by ultrasound can provide positive or negative effects on enzyme activities, mechanically, by altering the molecular structure of enzymes and their corresponding with substrates (, ). On the other hand, the free radicals produced by hydrolysis can interact with amino acid residues and affect the activity of enzymes (, ). There has been some evidence of the effect of ultrasound on gene expression in plant, for instance, changes of mRNA transcription in potato exposed to ultrasound were observed (). Noteworthily, a prior study in Arabidopsis exposed to sonication discovered transcriptomic changes, and changes in reactive oxygen species (ROS) scavenging ability, as well as in phytohormonal metabolism ().

To the best of our knowledge, however, few reports are available on the mechanism involved in ultrasound treatment inhibiting postharvest softening of strawberries. In this study, comprehensive analysis of physiological and transcriptomics were utilized to evaluate the effect of ultrasound on strawberry softening during 4°C storage. Furtherly, a broader systemic perspective on possible physiological and molecular regulatory effect of ultrasound treatment on strawberries softening can be obtained.

Materials and methods

Strawberry fruit materials

Strawberry fruits (Fragaria × ananassa Duch. Cv. “Benihoppe”), an octoploid cultivar were harvested from Da Zi Ran farm (Beijing, China), and were immediately transported to laboratory. Well-formed fresh fruits with similar maturity, size and color were selected and stored at 5°C for a night before the experiments.

Ultrasound treatment and storage

The selected strawberry fruits were completely immersed in a sonicator bath (JXD-02, Beijing Jinxing Ultrasonic Equipment Technology Co., Ltd., China) filled with deionized water for ultrasonic treatment (Frequency 28 kHz, Power 0.15 W/cm2, Time 3 min), the surface of water in the bath was kept at the same level and water temperature (15 ± 1°C) was maintained using a low-temperature thermostatic water bath (DC-2006, Ningbo Scientz Biotechnology Co., Ningbo, China). The ultrasound parameter was derived from the orthogonal design experiment in our previous research. Control strawberries samples were subjected to the same conditions and completely immersed in deionized water for the same time as the U treatment. After the treatment, all fruits were dried with absorbent paper and then placed in polyethylene terephthalate (PET) trays (15 fruit per tray) immediately. Later, packaged fruits were stored at 5 ± 1°C and 75–85% RH for 24 d.

Experimental design

Strawberry fruits after treatments were randomly sampled at different storage time (0 h, 12 h, 1 d, 2 d, 3 d, 6 d, 9 d, 12 d, 15 d, 18 d, 21 d, 24 d) for measurements of hardness, soluble pectin content (SP), PE and PG activity. Among these physiological indicators, fruit samples were homogenized and stored at −20°C after quick-frozen in liquid nitrogen for the determination of PE and PG activity, all the other indicators were measured using fresh samples. Two trays of strawberries were used to measure ethylene production. Signal substance H2O2 content was measured at 0 h, 3 h, 6 h, 9 h, 12 h, 1 d, 2 d, 3 d, 6 d, 9 d, 12 d, 15 d, 18 d, 21 d, 24 d. Due to the decay of strawberries at 24 d, the ethylene production was no longer measured, and the other time points for ethylene measurement was the same as the H2O2 content.

Strawberry samples in U and CK group were randomly selected at 0 h, 3 h, 6 h, 9 h, 12 h, 1 d, 2 d, 3 d, 6 d, 9 d during storage, and then the samples were pulped and frozen with liquid nitrogen to store at −80°C for transcriptome sequencing.

Three biological replicates were set in physiological and transcriptome sequencing experiments.

Firmness and soluble pectin content

Firmness was measured by a digital sclerometer (GY-4, Hangzhou Lvbo Instrument Co., Ltd., China) with 7.9 mm diameter probe, and the actual firmness was twice the shown maximum penetration force (N). Two opposite sides of fruit along the main axis were measured, and five strawberries from each treatment were randomly selected.

Extraction and assay of soluble pectin content was performed according to Cao et al. (). Strawberry fruit (1.0 g) was prepared as alcohol-insoluble material after repeated washing with 95% ethanol for 4 times. Then, the washed material was incubated in 50°C water for 30 min to obtain the soluble pectin extract. The soluble pectin content was determined with carbazole methods colorimetrically using a microplate spectrophotometer (Infinite 200 PRO, Tecan, Switzerland). Galacturonic acid was used as a calibration standard, and soluble pectin content was expressed as galacturonic acid equivalents (GaE).

Measurement of PE and PG activity

Extraction and measurement of PE (E 3.1.1.11) was determined by the method of Vicente et al. () with modifications. Frozen strawberry pulp (5 g) was homogenized with 15 ml of 1 M NaCl solution containing 1% (w/v) polyvinylpolypyrrolidone (PVPP). The mixture was stored at 4°C for 2 h, and then centrifuged at 4°C and 10,000 × g for 10 min. The pH of supernatant was adjusted to 7.5 with 0.1 M NaOH, the resulting enzymatic extract was used for assaying PE activity. Substrate mixture containing 2.8 mL 1% (w/v) partially methylated (≥85% esterified) citrus pectin solution pH = 7.5, 750 μL 30 mM potassium phosphate buffer (PH = 7.5) containing 0.01% (w/v) bromothymol blue. Reaction-mixture after adding 500 μL PE extract was mixed immediately and used for spectrophotometric measuration (UV spectrophotometer, UVmini-1240, Shimadzu Corporation, Japan) every 20 s during 3 min at 620 nm. Calibration curve was built using different concentration of galacturonic acid (pH = 7.5) as standard. PE activity was expressed as μmol of demethylgalacturonic acid (GalA) produced per second and per kilogram of fruit (μmol·s−1·kg−1).

The methodology to measure PG (E 3.2.1.15) activity was adapted from Karakurt and Huber (). Frozen strawberry pulp (10 g) was homogenized with 20 mL of 95% ethanol, placed at 4°C for 10 min, and then centrifuged at 8,000 rpm for 10 min. The precipitate was added with 10 mL 80% ethanol, placed at 4°C for 10 min, and then centrifuged at 8,000 rpm for 5 min. Then, the precipitate was extracted with 5 mL of buffer (50 mM sodium acetate/acetic acid pH 5.5) containing 1.8 M NaCl for 20 min and then centrifuged at 8,000 rpm for 5 min. The supernatant was adjusted to 10 mL, which was the PG extract for activity determination. All the steps were done at 4°C. Reaction mixtures for the determination of PG activity, using polygalacturonic acid as substrate (Shanghai Yuanye Biotechnology Co., Ltd., China.), consisted of 250 μL of the enzymic extract and 250 μL of 10 g/L polygalacturonic acid in 500 μL of 50 mM, pH 5.5 acetic acid-sodium acetate buffer, were incubated for 1 h at 37°C. Enzymic product galacturonic acid was determined by the 3, 5-dinitrosalicylic acid (DNS) method (), and absorbance at 540 nm was measured on microplate spectrophotometer (Infinite 200 PRO, Tecan, Switzerland). Calibration curve was built using different concentration of galacturonic acid (pH = 7.5) as standard. PG activity was expressed as μg of D-(+)-galacturonic acid produced per hour and per gram of fruit (μg·h−1·g−1).

Determination of H2O2 and ethylene production

H2O2 was assayed using the method of Brennan and Frenkel () with slight modifications. In brief, pulped strawberry (4 g) mixed with 6 mL acetone, and then centrifuged at 8,000 rpm 4°C for 5 min. The supernatant was made up to 10 mL with acetone, which was the H2O2 extract. One milliliter H2O2 extract was reacted with 0.1 mL 5% titanium sulfate-sulfuric acid (2 M) and 0.2 mL concentrated ammonia. The precipitate after centrifugation was washed by acetone and finally dissolved in 1 mL 2 M H2SO4. The absorbance at 415 nm was recorded and the H2O2 content results was expressed as μmol g−1 FW.

Ethylene production was measured as described by Mditshwa et al. () with slight modifications. Briefly, five fruits per replicate were used to measure ethylene production. Fruits was weighed using an electronic balance with an accuracy of 0.01 g, and thereafter enclosed in a 400 mL airtight glass jar for 1 h at 20°C. For each ethylene measurement, three jars also were used. SCS56 ethylene analyzer (Storage Control Systems, Inc., USA) was used for ethylene measurements and the results were expressed as μL·kg−1·h−1.

RNA extraction and transcriptome sequencing

Total RNA from strawberry samples was extracted using OmniPlant RNA Kit (DNase I) (ComWin Biotech Co., Ltd., China) according to the manufacturer's instruction. Transcriptome sequencing was performed by Shanghai Majorbio Biopharm Technology Co., Ltd. (Shanghai, China). The cDNA library was constructed using the TruSeq™ RNA sample preparation kit (Illumina Inc.). Briefly, 1 μg of total qualified RNA was used for poly-A based mRNA enrichment selection using Oligo (dT) magnetic beads, and then fragmentation buffer was added to randomly interrupt the enriched complete mRNA, small fragments of about 300 bp were separated by bead screening. Under the action of reverse transcriptase and random hexamers (six-base random primers), reverse transcription to synthesize one-strand cDNA was performed using mRNA as template, and then followed by double stranded cDNA synthesis. Double stranded cDNA was end repaired adding “A” tails to the 3′ end for adaptor ligation. Finally, transcriptome sequencing was performed on the Illumina Novaseq 6000 platform with 150 bp paired-end sequencing length. Subsequent bioinformatics analysis was based on clean, high-quality reads obtained by quality control. All raw sequence read data were submitted in the NCBI Sequence Read Archive (SRA, http://www.ncbi.nlm.nih.gov/Traces/sra) database with the accession number PRJNA613956.

Bioinformatics analysis

The clean data were mapped to the octoploid strawberry reference genome (https://www.rosaceae.org/species/Fragaria_x_ananassa/genome_v1.0.a1) using HISAT2 (, ). All genes were annotated in GO, KEGG, COG, NR, Swiss-Prot, and Pfam databases. Gene expression quantitative analysis was carried out using RSEM with TPM (Transcripts Per Million reads) as the quantitative indicator (). DESeq2 software was used to analyze the differences between different samples and treatment groups (). Differentially expressed genes (DEGs) were identified when |log2 fold change| ≥ 1 and p-adjust < 0.05 using the Benjamin-Hochberg multiple hypothesis test correction method (). KEGG enrichment analysis of the DEGs were performed using Goatools software base on R script, KEGG pathway was considered to be significantly enriched in the gene set when p-adjust < 0.05 in Fisher test.

qRT-PCR validation

In order to confirm the results of the transcriptome sequencing, 9 genes related to ROS, ethylene, auxin signaling pathways, MAPK signal transduction pathway, and fruit softening metabolic pathways, were selected for qRT-PCR (quantitative real-time PCR) analysis. Primer sequence were designed by Primer Premier 5 and their specificity were validated. Synthetic primers were purchased from Synbio Technology Co., Ltd. (China) and shown in Supplementary Table S1. The reference gene FaActin primers were F: 5′-GGTGACGAGGCTCAATCCAA-3′, R: 5′-GGGCAACACGAAGCTCATTG-3′.

cDNA was synthesized from quantified total RNA by HiScript® QRT SuperMix for qPCR (+gDNA wiper) (Vazyme Biotech Co., Ltd., Nanjing, China) according to the manufacturer's instruction. Synthesized cDNA (50 ng/μL) which were frozen at −20°C, were diluted 10-fold used as the template for qRT-PCR. qRT-PCR was performed using SuperReal PreMix Plus (SYBR Green) (Tiangen Biotech Co., Ltd., China) in CFX96 Touch™ Real-Time PCR Detection System (Bio-Rad Laboratories, USA). The 10 μL reaction mixture contained 0.3 μL of forward primer (5 μM), 0.3 μL reverse primer (5 μM), 5 ng cDNA template, 5 μL 2 × SuperReal Premix Plus (SYBR Green), and 3.4 μL water. The qRT-PCR conditions were as follows: 15 min at 95°C, followed by 40 cycles of 10 s at 95°C and 20 s at 55°C, and 20 s at 72°C. The melting curves were measured from 65 to 95°C at 0.5°C increments. Relative expression levels of target gene were normalized to internal control FaActin, with CK_0h sample as the benchmark, and were calculated by 2−ΔΔCt method (). For each gene analyzed, a template-free negative control was included in each run, and three independent biological replicates were made for each gene in each sample.

Statistical analysis

Bioinformatics analysis and graphical presentations of RNA-seq data was completed with online data processing and drawing tools in website (http://www.majorbio.com/). Physiological data was statistically analyzed by independent t-tests using SPSS 25 (IBM, Armonk, NY, USA), and differences at P < 0.05 level were considered statistically significant. Figures were made by Origin 2019 (Microcal Software Inc., Northampton, MA, USA). The data in the figures were presented as means ± standard errors (SE).

Results

Firmness and soluble pectin content after ultrasound treatment

Fruit firmness is considered to be one of the most important indicators of fruit quality, and it decreases gradually along with fruit ripening and softening during postharvest storage. There was a general decrease in fruit firmness of both the control and ultrasound-treated strawberries during cold storage at 24 d (Figure 1). By the end of storage at day 24, the firmness of the control fruits decreased from 17.36 to 10.52 N, a decrease of 39.40%, while the firmness of the ultrasound-treated fruits decreased only by 22.97%, from 16.98 to 13.08 N. Relative to the control, the firmness of the ultrasound treated strawberries was maintained at a significant (p < 0.05) higher level after 18 d. The results showed that the decline in strawberry fruit firmness was inhibited by ultrasound treatment during storage, which was in accordance with the study of Cao et al. (), Alexandre et al. (), and Gani et al. (), finding that ultrasound can effectively maintain the textural quality of strawberry.

Figure 1

The increase of soluble pectin (SP) content in plant cell walls is one of the main reasons for fruit firmness decrease (). According to Figure 2B, SP content increased with postharvest storage time in two groups. SP content of the control fruits increased from 2.63 to 4.63 mg·g−1, an increase of 76.05%, whereas the SP content of the ultrasound- treated fruits increased by 28.05 %, from 3.03 to 3.88 mg·g−1. Compared with the control, SP content in ultrasound-treated strawberries was maintained at a significant (p < 0.05) lower level generally after 9 d. The results indicated that ultrasound inhibited the production of SP in late storage period, thereby maintained fruit firmness and retarded strawberry fruit softening effectively.

Figure 2

Activities of PE and PG after ultrasound treatment

PE can catalyze the deesterification of polygalacturonic acid in plant cell walls into pectic acid, forming a substrate that facilitates the action of PG, PL and other pectinase, promote the degradation of pectin (). As shown in Figure 3A, PE activity in ultrasound-treated strawberries showed the same trend as the control. PE activity decreased during the first 9 d, increased rapidly from days 9-12, and then decreased again after day 12. There was no significant difference in PE activity between U and CK group during the first 6 d, and significant differences appeared since 9 d. The PE activity of strawberries in U and CK group reached the highest level at 15 and 12 d respectively, which were 4.17 (U) and 4.51 (CK) μmol·s−1·kg−1, decreased by 7.54% compared to the control. Moreover, the PE activity of strawberries in U group was significantly lower than CK (P < 0.05) after day 18.

Figure 3

PG is one of the main hydrolytic enzymes of plant cell walls, which can further hydrolyze the pectic acid in the fruit cell wall to generate oligomeric galacturonic acid or galacturonic acid, leading to fruit softening. In the two groups, PG activity increased on 0–2 d and then decreased during the first 6 d of storage (Figure 3B) approximately, reached a peak after day 9 and decreased gradually thereafter. There was no significant difference in PG activity between U and CK group during the first 6 d. The PG activity of strawberries in U and CK group peaked on days 15 and 12 respectively, which were 490.20 (U) and 531.41 (CK) μg·h−1·g−1. Compared with the control, the highest PG activity in ultrasound treated strawberries decreased by 7.75%, and PG activity was also maintained at a significantly lower level since 18 d. These results were similar to the performance of PE activity. In the present study, the time when the PE and PG activity reached the highest level was delayed by ultrasound treatment, and their activities in late storage period were also suppressed by ultrasound.

H2O2 production and ethylene production after ultrasound treatment

H2O2 is an important reactive oxygen species (ROS) component, which acts as an important signal substance to regulate the quality of fruits and vegetables during postharvest storage. H2O2 content was monitored through time for both the control and ultrasound treated strawberries (Figure 4A). It can be seen that the H2O2 content of ultrasound treated strawberries was significantly lower than that in control at 0 h, 3 h, indicating that ultrasound significantly inhibited H2O2 production in the early storage period. From 6 h to 15 d, the H2O2 content almost had no obvious change in two groups, and had no significant difference in strawberries of two groups except on day 2 and 6. Since day 18, the H2O2 content increased rapidly along with storage time, and the H2O2 content of the ultrasound group was significantly lower than that of control, indicating that ultrasound reduced the H2O2 content of the strawberry in the late storage period. The results exhibit that ultrasound affects the ROS metabolism of strawberry fruits mainly by reducing the H2O2 production in the early and late stages of storage.

Figure 4

Ethylene (ET) is a gaseous hormone that regulates fruit maturation and senescence. Studies have found that ethylene can be used as a signal substance to regulate the texture, color and aroma of non-climacteric fruits. As shown in Figure 4B, during cold storage period after 0 h, the ethylene production of strawberry fruits showed an upward trend. During the first 6 d, the ethylene production of strawberries in the U group was generally significantly higher than that in CK. No significant difference was observed between U and CK group during 9–18 d, while the ethylene production of U group was significantly (p < 0.05) lower than CK at 21 d finally.

Transcriptome sequencing and DEGs analysis

A total of 250.05 Gb clean data was obtained from 60 samples, for each sample, above 6.27 Gb clean data (with Q30 > 93.68%) was obtained, and about 85.59–92.18% clean data was mapped to a reference Fragaria × ananassa genome (v1.0.a1). In total, 85, 634 expressed genes were detected, including 82, 205 known genes and 3, 429 new genes. Differentially expressed genes (DEGs) were identified by comparing the transcripts of the U group samples with CK at the same storage time. A total of 1,905 DEGs were identified at all the 10 time intervals, of these, 714 genes were significantly up-regulated and 1,254 genes were significantly down-regulated. The number of DEGs in samples varies with storage time greatly (Table 1; Figure 5). It was found that 498, 739, and 769 genes differentially expressed at 0, 6, and 12 h after ultrasound treatment, respectively. While the number of DEGs decreased gradually after 12 h, which means that the majority of DEGs (up to 769/1,905) respond to ultrasound treatment in the early stage of storage.

Table 1

Storage timeDEGsUp-regulated DEGsDown-regulated DEGs
0 h498118380
3 h624
6 h739375364
9 h261214
12 h769210559
1 d541242
2 d351322
3 d431
6 d743
9 d26818
Total DEGs1,9057141,254

Numbers of DEGs at different storage time after ultrasound treatment.

Figure 5

DEGs GO annotation and KEGG enrichment analysis

All the 1,905 DEGs were subjected to GO functional annotation analysis to determine the biological process, cellular component and molecular function of genes (Figure 6). Among these, 609 genes and 495 genes were associated with metabolic processes and cellular processes in biological process, respectively; 294 genes and 275 genes were associated with membrane part and cell part in cellular component, respectively; 752 genes and 682 genes were associated with binding function and catalytic activity in molecular function, respectively. In addition, there were 258 genes involved in the biological regulation process, and 80 genes owned transcription regulator activity.

Figure 6

For both up-regulated and down-regulated genes of all 1,905 DEGs, KEGG pathway enrichment analysis was performed (Figure 7). Of the upregulated genes, largest quantity of genes was involved in the plant-pathogen interaction and glycerophospholipid metabolism pathway; of the downregulated genes, largest quantity of genes was also involved in the plant-pathogen interaction metabolism pathway, and followed pathways were plant hormone signal transduction and MAPK signaling pathway. The KEGG pathway enrichment analysis showed that ultrasound treatment had a significant effect on the interaction between strawberry fruit and pathogens. The MAPK signaling pathway can amplify signals through gradual phosphorylation and participate in plant defense responses under stress. The phytohormone signal transduction pathway, MAPK signaling pathway, starch and sucrose metabolic pathway, and galactose metabolic pathway were significantly down-regulated during the first 9 d storage, indicating that the stress response and signal transduction process was milder and metabolism level was lower in ultrasound treated strawberries than in control, which helping to maintain fruit quality.

Figure 7

Expression profile of genes in ROS, hormone metabolism and signaling pathways and transcription factor family

A total of 13 genes related to ROS metabolism was investigated (Figure 8). Compared with the control group, four genes differently expressed significantly at 0 h, of which three genes, and encoding LOX were upregulated, and one gene encoding RBOH (NADPH oxidase) () was downregulated in the ultrasound-treated group. ROS scavenging of genes encoding LOX was upregulated at 0 h, which was the earliest gene to promote ROS scavenging. Down-regulation of genes encoding RBOH which can catalyze the production of superoxide anion and H2O2 showing that ROS production was inhibited at 0 h. At 6 h after ultrasound treatment, six DEGs were observed. Four of which were upregulated, including one gene encoding PPO, one gene encoding 4CL, and two genes encoding POD. Another two genes were downregulated, one gene encoding POD and the other gene encoding HCT. Genes encoding PPO and POD which related to ROS scavenging were partially upregulated, and gene encoding 4CL which can promote phenol synthesis to reduce ROS was upregulated, which mean that ROS production continued to be inhibited and ROS scavenging was enhanced at 6 h. At 12 h, two genes were downregulated, one gene encoding POD and the other gene encoding RBOH. These results indicated that the production of ROS was suppressed and the elimination of ROS was enhanced continuously during the first 12 h after ultrasound treatment.

Figure 8

Nine genes involved in JA metabolism and signaling were differentially expressed during storage (Figure 8). At 0 h, three genes coding JAZ protein which can inhibit JA signaling showed downregulation, resulted in the increasing release of JA signal. At 6 h, the expression of three genes were remarkably regulated. The transcripts of two JMT genes which coding JA degrading enzyme jasmonate O-methyltransferase (JMT), were dramatically increased, whereas transcript abundance of IAR3 which coding IAA-amino acid hydrolase displayed decline, indicating the inhibition of JA synthesis and enhancement of JA degradation to certain degree at 6 h. At 12 h, two JAZ genes exhibited significant downregulation, suggested that the release of JA signal was enhanced again at which time.

There were 18 ethylene related genes displayed markedly differences in expression during storage (Figure 8). At 0 h, six ethylene response factor genes were significantly down-regulated which including three genes of ERF109 and three genes of ACS1_2_6. The ethylene signal pathway in the ultrasound treatment group began to respond at 0 h, the control group began to up-regulate from 12 h, and the ultrasound group began to significantly up-regulate at 2 days. The ethylene signal pathway response after ultrasound was prolonged and lagged. Ethylene is synthesized from S-aenosyl methionine (SAM) via the intermediate 1-aminocyclopropane-1-carboxylic acid (ACC), in the biosynthesis process, SAM is converted to ACC by ACC synthase (ACS) and ACC to ethylene by ACC oxidase (ACO). Additionally, ERFs are critical downstream components in ethylene signaling pathway (). It was reported that the accumulation of endogenous ethylene can suppress ACS expression, and the transcript abundance of genes coding ethylene synthesis-related enzyme regulated their enzymic activities, influencing the production of ethylene. Thus, the downregulation of two ACS genes and upregulation of one ACO gene, was coincident with the performance of higher ethylene production in ultrasound treated strawberries within first 12 h of storage.

Transcription factors (TFs), are DNA-binding proteins that can specifically bind to the promoter region of eukaryotic genes or interact with certain auxiliary regulators, thereby regulating the specific expression of plant genes (). According to the differences of DNA binding domains, TFs are primarily classified into several families, including MYB, AP2/ERF, WRKY, NAC, MADS-box, Homeobox, GRAS, bZIP, bHLH, etc. (, ). WRKY TFs compose one of the largest TF families in plants (). Differential expression analysis showed that seven WRKY TFs displayed dramatic shifts in their mRNA levels (Figure 8). A WRKY33 showed reduced transcription at 0 h. At 12 h after ultrasound treatment, six WRKY TFs including two WRKY29 genes and four WRKY33 genes were significantly down-regulated. Prior researches have revealed the vital regulatory function of WRKY TFs in response to biotic or abiotic stress (). The differential expression of WRKY TFs indicated that WRKY29 and WRKY33 may play an important regulatory role in response to ultrasound treatment within 12 h.

qRT-PCR analysis

To validate the transcriptome sequencing data, nine candidate genes were randomly selected for qRT-PCR verification (Figure 9). The genes included one gene related to ROS metabolism pathway, three genes involved in auxin and ethylene signaling pathway, one gene related to MAPK signaling pathway, involved in fruit softening, and two genes involved in physiological metabolism of strawberry. For the nine genes, the results of change with storage time, and difference between two groups in gene expression obtained by qRT-PCR, were consistent with the results obtained from transcriptome sequencing analysis.

Figure 9

Discussion

Effect of ultrasound treatment on physiological characteristics and transcriptome process related to strawberry fruit postharvest softening

Softening is a vital factor that causes the deterioration of strawberry fruit quality. In this study, the firmness of postharvest strawberry fruit decreased gradually with storage time. However, ultrasound treatment significantly inhibited the firmness loss in the strawberry fruit during low temperature storage, accompanied with decreased SP content, PE and PG activity after 9 d of storage (Figure 2). Similar results were also observed in previous studies on ultrasound treated strawberry fruit (, , , ), and other fruit species, such as tomato (, ), cherry (), and jujube (). The inactivation effect of ultrasound on PE and PG has been studied to some extent, Raviyan et al. () found that the inactivation of PE in tomato increased with cavitation intensity, and similar study on PG indicated that high-intensity and prolonged ultrasound (22 kHz, >8.1 W·ml−1, >35 min) could induce the PG inactivation process (). Explanations for PE, PG inactivation could be that free radicals and intense shear forces generated from extreme ultrasound conditions led to enzyme structural changes, and meanwhile changed the reactions with substrates and their surroundings (, ). Application of relatively low-intensity ultrasound (28 kHz, 0.15 W/cm2, 3 min) in this work, similarly, suppressed the activity levels of PG and PE. This result might be possible that inactivation effect on PE and PG of different ultrasound intensity varies from fruit species, and is influenced by regulation of gene expression and plant hormone in living plant.

ROS and ET, two kinds of important signal molecules, were investigated in this work. Our results indicated that ultrasound-treated strawberries showed lower H2O2 levels than the control within the initial 3 h and after day 18 of storage, although the H2O2 level was higher on days 2 and 6 (Figure 4A). As an important component of ROS signal, we speculated that ultrasound activated ROS metabolism to reduce H2O2, subsequently regulated the downstream response to ultrasound. Moreover, ultrasound significantly inhibited the production of ethylene, thereby delaying the softening and senescence of strawberry fruit (Figure 4B). Notably, ethylene production was higher in strawberry fruit treated with ultrasound during the first 6 d, this could be explained that ethylene production is one of the most common and marked reactions to abiotic stress abiotic stress in plant biology (). The result of ethylene production was similar to previous report claiming that higher level of ethylene was observed in ultrasound-treated carrots ().

According to transcriptomic results, changes in gene expression of several plant signaling pathway and TFs family were detected. In ROS metabolism and signaling pathway, gene expression related to ROS production was down-regulated, meanwhile gene expression linked to ROS scavenging was up-regulated at 0, 6, and 12 h, which indicated the inhibition of excessive ROS production and enhancement of ROS scavenging during the first 12 h (Figure 8). The results of ROS gene expression were aligned with the physiological response that H2O2 content of ultrasound-treated strawberries was lower within the initial 3 h of storage. The upregulation of ROS scavenging genes in ultrasound treated strawberry, was similar to the finding that DEGs encoding ROS scavenging enzyme GST were up-regulated in potato in vitro exposed to piezoelectric ultrasound (). The differential expression of JA related genes suggests that there was an apparent response to ultrasound in JA metabolism and signaling pathway during the first 12 h, the release of JA was activated at 0 and 12 h while suppressed momently at 6 h. As shown in Figure 8, the general up-regulation of ethylene-related genes started at 12 h in the control strawberries, while started by 2 d in ultrasound-treated strawberries, which may cause to the relatively higher ethylene production of ultrasound-treated strawberries between days 2 and 6. These results suggest that the ultrasound treatment resulted in an earlier increase in ethylene production, which is associated closely with the regulation that ET metabolism and signaling pathway was activated until day 2. Interestingly, changes in expression profiles of auxin-related genes were also detected during the first 12 h, but the response was too complicated to explained now. In addition, we also found the differential expression of TFs genes including WRKY and ERF. WRKY33 and WRKY29 genes underwent strong downregulation in the initial 12 h, suggesting that WRKY33 and WRKY29 may play an important role in regulating downstream signaling and metabolic genes.

Physiological mechanisms of ultrasound treatment inhibiting strawberry fruit postharvest softening

In fact, as reported in numerous studies, there are extremely complicated interactions in phytohormone, TFs and cell wall degrading enzymes.

The earlier change in ROS (H2O2) content can be considered as an immediate response took place right after ultrasound application (). Acting as a signaling molecule, H2O2 may activate JA biosynthesis and signaling pathway (). H2O2 and JA signals may work together for activating the subsequent responses, such as increased ethylene production as late response in ultrasound-treated carrot (). There have been many studies demonstrated that ET can affect the gene expression of cell wall-degradation related enzymes in strawberry. With the application of ethylene, Villarreal et al. () found that FaPG1, FaGal1, and FaGal2 were up-regulated, whereas FaPME1, FaXyl1, and FaXTH1 were down-regulated, meanwhile FaPLa did not show any change at all in strawberry. Others reported the downregulation of FaPG2 and FaPLa, the upregulation of FaPE1, and no changes in FaPG1 in ethylene treated strawberry (). In present study, the ethylene production of ultrasound-treated strawberries was higher than the control during 2–6 d, although the differential expression of PE and PG, the enzyme activity of PE and PG were markedly lower than the control after 9 d. This could be caused by the influence of ultrasound on the structure of the existing PE PG enzymes in strawberries.

On the other hand, TFs including WRKY, ERF, NAC, and MYB family genes can regulate phytohormone biosynthesis genes and signal transduction (). ERFs, act as a key regulator, integrate the ET, ABA, JA, and redox signaling in plant, and conversely, JA and ABA have also been reported to participate in the regulation of ERFs under abiotic stresses (, ). Moreover, WRKY TFs are reported to be involved in JA mediated resistance (). Eulgem and Somssich () reported that WRKY33 can be activated by MAPK cascade pathway in Arabidopsis thaliana. However, the KEGG enrichment analysis revealed that MAPK signaling pathway was downregulated by ultrasound, which could be the reason why WRKY33 genes were downregulated within the first 12 h. Additionally, extensive studies have shown that ERF genes are involved in the regulation of cell wall-modifying genes and fruit softening. In papaya, CpERF9 acts as a transcriptional repressor of CpPME1/2 and CpPG5 (). In peach, PpeERF2 can suppress the expression of PpePG1 by binding to its promotors ().

Based on the results presented in present work and prior researches, the mechanism ultrasound treatment inhibiting strawberry fruits postharvest softening was shown as Figure 10. After ultrasound treatment, the elimination of ROS (H2O2) and the release of JA signal was enhanced in the early stage (within 12 h). Under the effect of ROS and JA signal, ethylene signaling pathway was activated later (1.5 d later, at 2 d), resulting in longer release time of ET signal during the early middle stage (2–6 d), which is a likely important process for decreased PE and PG capacity in late storage time (after 9 d), inhibiting strawberry fruit postharvest softening ultimately (18–24 d). Herein, the transcription factors WRKY29 and WRKY33 may also play an important regulatory role. For further study, the hypothetical mechanism remains to be elucidated, through a chemical biology approach by using activators or inhibitors of the targeted signal molecules.

Figure 10

Conclusion

Postharvest softening in strawberry fruits was effectively inhibited by ultrasound treatment. Ultrasound maintained the firmness and reduced the production of soluble pectin by retarding and suppressing the activities of pectin-degrading enzymes such as PE and PG, and therefore helping to keep textural quality of strawberry. Moreover, the content of signal substances hydrogen peroxide (H2O2) and ethylene had response to ultrasound treatment in early and late storage stage. Transcriptomic analysis revealed that ultrasound significantly changed the expression of genes related to ROS, ET, JA as well as other hormone signaling pathways and WRKY family including RBOH, LOX, PPO, POD, ACS1_2_6, ACO, ERF, JAZ, IAA, SAUR, ARF, together with WRKY29 and WRKY33 etc. A hypothetical mechanism model for the postharvest softening inhibition in strawberry by ultrasound treatment was proposed based on physiological and transcriptomic results of present study, as well as evidence from previous reports.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The data presented in the study are deposited in the NCBI repository, accession numbers BioProject ID PRJNA895742 and SRA ID SRP405228.

Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.1066043/full#supplementary-material

    Abbreviations

  • ACO

    1-aminocyclopropane-1-carboxylate acid oxidase

  • ACS

    1-aminocyclopropane-1-carboxylate acid synthase

  • ARF

    auxin response factor

  • cDNA

    Complementary DNA

  • DEGs

    differentially expressed genes

  • FC

    fold change

  • GO

    Gene Ontology

  • JAZ

    JASMONATE ZIM-domain protein

  • KEGG

    Kyoto Encyclopedia of Genes and Genomes

  • LOX

    Lipoxygenase

  • MAPK

    Mitogen activated protein kinase

  • NR

    Non-Redundant Protein Sequence Database

  • NCBI

    National Center for Biotechnology Information

  • PE

    pectin methylesterase

  • PG

    polygalacturonase

  • qRT-PCR

    real-time quantitative reverse transcription polymerase chain reaction

  • ROS

    Reactive oxygen species

  • TF

    Transcription factor

  • RNA-seq

    RNA sequencing

  • NCED

    9-cis-epoxy carotenoid dioxygenase.

References

Appendix A

Table A1

Gene IDGene nameSequences
FxaC_27g39130HSPFTCTTCGACCCATTTTCCCTC
RCAGCTCAACCTTCACTTCTTCC
FxaC_6g03070HydFAGTTTCCCTTTTGGTTGTGC
RCGTCCCTGGATACTCTCTGTG
FxaC_4g12140PGFCCAGTAACTTCACGCTTTTTCTC
RTCCTAACCGCTCTCTCCCTC
FxaC_11g01590PPOFCGAGGAGCAAGAATGAAAAGG
RGCGGCAACAACACCAAGC
FxaC_1g11640CHIFTCGGCAGAACAATGCGGTAGA
RTTGGGGTTGGTGATGAGGGAG
FxaC_1g17390MAP3K3FCCATCTACAACCTATGCCTCTA
RTGGGAATCATCACTGAACTATC
FxaC_26g07610ERF027FGATACGGCCCTCAATTTCCCC
RTCCACATTTTCAACCCCACCC
FxaC_8g31020ARFFTAACCAACAGGAACAGAGAATAATG
RCAAAAGAGGAGGATGTGAAATACCA
FxaC_9g45470ARPFTGGAGAGACGAGGAGCAAGC
RACGACGAACCCTATTCACAAAAC

Primers used for qRT-PCR in strawberry fruits.

The raw Illumina mRNA-seq data were submitted to NCBI and the processed data were deposited under GEO ID GSE123176, GSE135294, BioProject ID PRJNA507769, PRJNA558262, PRJNA895742 and SRA ID SRP171630, SRP217133.

Summary

Keywords

ultrasound, strawberry, transcriptome sequencing, postharvest softening, mechanism

Citation

Zhang J, Jiang H, Li Y, Wang S, Wang B, Xiao J and Cao Y (2022) Transcriptomic and physiological analysis reveals the possible mechanism of ultrasound inhibiting strawberry (Fragaria × ananassa Duch.) postharvest softening. Front. Nutr. 9:1066043. doi: 10.3389/fnut.2022.1066043

Received

10 October 2022

Accepted

16 November 2022

Published

01 December 2022

Volume

9 - 2022

Edited by

Yongsheng Chen, Jinan University, China

Reviewed by

Qiya Yang, Jiangsu University, China; Yue Wang, Zhejiang University, China

Updates

Copyright

*Correspondence: Shaojia Wang

This article was submitted to Nutrition and Food Science Technology, a section of the journal Frontiers in Nutrition

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics