Physiological and Transcriptome Analyses of CaCl2 Treatment to Alleviate Chilling Injury in Pineapple

The post-harvest ripening of pineapples can be effectively postponed by refrigerated storage. Nevertheless, internal browning (IB) frequently appears in pineapples after refrigerated storage during the course of the shelf life at room temperature, which is known as chilling injury (CI). In this study, the chilling injury of pineapple fruit was induced by a low temperature (6 °C) and transferred to normal-temperature storage; the best concentration of 50 μmol/L of CaCl2 was selected by the IB appearance and electrical conductivity. Fruit quality, reactive oxygen species (ROS), antioxidants, and transcription factors were investigated. The physiological data results indicated that pineapples treated with 50 μmol/L of CaCl2 maintained fruit quality, decreased reactive oxygen species (ROS), and enhanced the antioxidant activity of fruits, alleviating internal browning (IB) symptoms in pineapple fruit. The expressions of related genes were also consistent with the physiological changes by the transcriptome data analysis. In addition, we focused on some related metabolic pathways, including phenylpropanoid biosynthesis, MAPK pathway, plant hormone, plant–pathogen interaction, tricarboxylic acid cycle (TAC), and fatty acid biosynthesis. We performed integrative analyses of transcriptome data combined with a series of physiology and experimental analyses on the internal browning of pineapples, which will be of great significance to extending the shelf life of pineapples through molecular biology in the future.


Introduction
The pineapple (Ananas comosus L. Merr.) is one of three major tropical economic crops (among the most significant harvest crops) and a tasty tropical fruit [1]. To ensure long-distance transport trade, the harvested fruits are usually stored in low-temperature conditions. However, it is easy to induce the internal browning of pineapple fruit, which has a negative effect on the commodity value and the price of this fruit [2]. In addition to pineapple, these unfavorable changes (caused by cold storage) are common in many other harvested horticultural fruits, including apples, pears, peaches, and bananas [3][4][5][6],. Internal browning, also termed endogenous browning or black heart, is a physiological disorder of pineapple fruit [7]. To date, the mechanism of harvested fruit browning has not been elucidated under cold storage.
Cold storage for tropical fruits leads to a series of physiological disorders within plant cells, including dehydration, flesh browning, and so on [8,9]. These processes are also closely related to the compositions of various flavor-related metabolites and phytohormones, such as sugars, fatty acids (FAs), amino acids, ethylene (ET), jasmonate (JA), abscisic acid (ABA), salicylic acid (SA), and gibberellic acid (GA), which lead to flavor deterioration and cold-stress induced modulation in plants [3,10]. Several researchers have observed that flavor-related metabolite compositions (sugars) will stabilize the membrane Moreover, 50 μmol/L of CaCl2 can strikingly decrease the relative conduc ple fruit and reduce the membrane damage to fruit storage.
By the IB appearance and electrical conductivity, 50 μmol/L of CaCl centration selection. Therefore, the following physiological indicators an data are only for the 50 μmol/L CaCl2-treated group and the control grou  (B) in pineapple fruit during s and 2 w) after transfer to room temperature (3 and 6 d), and treatment with CaCl μmol/L) or water (CK). Vertical bars represent the standard error of the mean of t Characters at the same storage times represent significant differences (p < 0.05).

Membrane ROS Production and Membrane Lipid Peroxidation
As shown in (Figure 2A,B), the O2 • − production rate and H2O2 conten ited overall increments in fruit DP during room temperature storage, bu production rate, there was no significant difference ( Figure 2A). The H2O control group was more than the 50 μmol/L CaCl2-treated group in the HP The DP control group was significantly greater than the 50 μmol/L CaC after 2 w from 0.57 ± 0.014 umol/g of FW to 0.73 ± 0.023 umol/g of FW. I μmol/L of CaCl2 reduced the ROS damage to fruit cells.
MDA content in both the control group and the 50 μmol/L of the CaC continuously increased, but DP had higher HP 3 and 6 d after the transfer tion to room temperature ( Figure 2C). This proves that pineapple fruit membrane lipid peroxidation during storage (from cold storage to room Figure 1. IB appearance (A), electrical conductivity (B) in pineapple fruit during storage at 6 • C (1 and 2 w) after transfer to room temperature (3 and 6 d), and treatment with CaCl 2 (50, 100, 150 µmol/L) or water (CK). Vertical bars represent the standard error of the mean of triplicate assays. Characters at the same storage times represent significant differences (p < 0.05).

Membrane ROS Production and Membrane Lipid Peroxidation
As shown in (Figure 2A,B), the O 2 • − production rate and H 2 O 2 content in fruit exhibited overall increments in fruit DP during room temperature storage, but in the HP O 2 • − production rate, there was no significant difference ( Figure 2A). The H 2 O 2 content of the control group was more than the 50 µmol/L CaCl 2 -treated group in the HP during storage. The DP control group was significantly greater than the 50 µmol/L CaCl 2 -treated group after 2 w from 0.57 ± 0.014 umol/g of FW to 0.73 ± 0.023 umol/g of FW. In conclusion, 50 µmol/L of CaCl 2 reduced the ROS damage to fruit cells.
MDA content in both the control group and the 50 µmol/L of the CaCl 2 -treated group continuously increased, but DP had higher HP 3 and 6 d after the transfer from refrigeration to room temperature ( Figure 2C). This proves that pineapple fruit is subjected to membrane lipid peroxidation during storage (from cold storage to room temperature).

Ascorbic Acid (ASA) and Glutathione (GSH) Contents
Initially, the ascorbic acid content of the control and the 50 μmol/L of CaC fruit were approximately 2.72 ± 0.160 g Kg −1 and 2.72 ± 0.071 g Kg −1 in refrigeratio 3A). Following the storage time elongation at room temperature, the overall in of the ASA content in HP. Moreover, the 50 μmol/L CaCl2-treated group wa higher than the control group. Comparatively, the DP was lower during most sto periods. However, 50 μmol/L of the CaCl2 group was more than the control grou In Figure 3B, regarding the GSH content, there are no significant changes i in DP, the 50 μmol/L CaCl2-treated group had a higher GSH content than th group; the value was from 0.099 ± 0.0013 g Kg −1 to 0.013 ± 0.0010 g Kg −1 after 1 w 0.089 ± 0.0034 g Kg −1 to 0.082 ± 0.0096 g Kg −1 after 2 w. Based on the change traje control group's HP decreased after 1 and 2 w, but the HP of the 50 μmol/L CaC group experienced a slower change (a 'rise again' phase); this was, respectively, and 2 w + 3 d. It shows that the 50 μmol/L CaCl2-treated group had more stress cold stress. during storage at 6 • C (1 and 2 w) after transfer to room temperature (3 and 6 d), and treatment with CaCl 2 (50 µmol/L) or water (CK). HP: "health part", DP: "disease part". Vertical bars represent the standard error of the mean of triplicate assays. Characters at the same storage times represent significant differences (p < 0.05).

Ascorbic Acid (ASA) and Glutathione (GSH) Contents
Initially, the ascorbic acid content of the control and the 50 µmol/L of CaCl 2 -treated fruit were approximately 2.72 ± 0.160 g Kg −1 and 2.72 ± 0.071 g Kg −1 in refrigeration ( Figure 3A). Following the storage time elongation at room temperature, the overall increments of the ASA content in HP. Moreover, the 50 µmol/L CaCl 2 -treated group was sharply higher than the control group. Comparatively, the DP was lower during most storage time periods. However, 50 µmol/L of the CaCl 2 group was more than the control group.
In Figure 3B, regarding the GSH content, there are no significant changes in HP, but in DP, the 50 µmol/L CaCl 2 -treated group had a higher GSH content than the control group; the value was from 0.099 ± 0.0013 g Kg −1 to 0.013 ± 0.0010 g Kg −1 after 1 w and from 0.089 ± 0.0034 g Kg −1 to 0.082 ± 0.0096 g Kg −1 after 2 w. Based on the change trajectory, the control group's HP decreased after 1 and 2 w, but the HP of the 50 µmol/L CaCl 2 -treated group experienced a slower change (a 'rise again' phase); this was, respectively, 1 w + 6 d and 2 w + 3 d. It shows that the 50 µmol/L CaCl 2 -treated group had more stress ability to cold stress.  . ASA content (A) and GSH content (B) in pineapple fruit during storage at 6 °C (1 and 2 w) after transfer to room temperature (3 and 6 d), and treatment with CaCl2 (50 μmol/L) or water (CK). HP: "health part", DP: "disease part". Vertical bars represent the standard error of the mean of triplicate assays. Characters at the same storage times represent significant differences (p < 0.05).

Phenolic Metabolism and Active Oxidase
As shown in Figure 4A, total phenolic content decreased after an initial increase in the HP of the control group and the 50 μmol/L CaCl2-treated group after 1 w. Moreover, the HP of the 50 μmol/L CaCl2-treated group experienced relative lagging. Respectively, the HP initial values of the control group and the 50 μmol/L CaCl2-treated group were 0.589 ± 0.0406 g Kg −1 and 0.525 ± 0.0277 g Kg −1 (1 w + 0 d), then rose to 0.956 ± 0.0225 g Kg −1 and 0.818 ± 0.0811 g Kg −1 (1 w + 3 d), and finally went down to 0.229 ± 0.0189 g Kg −1 and 0.551 ± 0.0132 g Kg −1 (1 w + 6 d). The total phenolic content continuously decreased in the control group HP after 2 w; the value went from 0.738 ± 0.0143 g Kg −1 to 0.216 ± 0.0100 g Kg −1 . However, the HP of the 50 μmol/L CaCl2-treated group was first up and then down from 0.410 ± 0.0262 g Kg −1 (2 w + 0 d) to 0.836 ± 0.0151 g Kg −1 (2 w + 3 d) to 0.622 ± 0.0060 g Kg −1 (2 w + 6 d). Evidently, the HP of the 50 μmol/L CaCl2-treated group was higher than the HP of the control group (2 w + 3 d and 2 w + 6 d). There is no obvious difference in the fruit DP.
As shown in Figure 4B,C, the PPO activity of the DP fruit was higher than the HP fruit. For HP fruit, there was no difference for both the control group and the 50 μmol/L CaCl2-treated group in 1w; however, the 50 μmol/L CaCl2-treated group was lower at any storage time. The DP fruit (1 w + 3 d and 2 w + 6 d 50 μmol/L of the CaCl2-treated group) was lower than the control group. Separately, the values were 0.395 ± 0.0225 Ug −1 FW and 0.391 ± 0.0218 Ug −1 FW ( Figure 4B). To summarize, PPO activity was higher in DP than in HP. POD activity was HP and far outweighed DP, whether it was the control group or the 50 μmol/L CaCl2-treated group ( Figure 4C). The HP of the control group had similar changing trajectories as the 50 μmol/L CaCl2-treated group HP after 1 w. However, initially, the control group HP was higher than the HP of the 50 μmol/L CaCl2-treated group but after the control group HP was lower than the 50 μmol/L CaCl2-treated group HP, the values (separately) were 7.85 ± 0.280 Ug −1 FW and 17.53 ± 1.075 Ug −1 FW. There was no difference with DP; however, the 50 μmol/L CaCl2-treated group was slightly higher after 1 w. Both the control group and the 50 μmol/L CaCl2-treated group experienced obviously rises. Separately, there were increases from 4.59 ± 0.223 Ug −1 FW to 15.33 ± 0.328 Ug −1 FW and from 12.64 ± 0.858 Ug −1 FW to 25.72 ± 0.170 Ug −1 FW. . ASA content (A) and GSH content (B) in pineapple fruit during storage at 6 • C (1 and 2 w) after transfer to room temperature (3 and 6 d), and treatment with CaCl 2 (50 µmol/L) or water (CK). HP: "health part", DP: "disease part". Vertical bars represent the standard error of the mean of triplicate assays. Characters at the same storage times represent significant differences (p < 0.05).

Phenolic Metabolism and Active Oxidase
As shown in Figure 4A, total phenolic content decreased after an initial increase in the HP of the control group and the 50 µmol/L CaCl 2 -treated group after 1 w. Moreover, the HP of the 50 µmol/L CaCl 2 -treated group experienced relative lagging. Respectively, the HP initial values of the control group and the 50 µmol/L CaCl 2 -treated group were 0.589 ± 0.0406 g Kg −1 and 0.525 ± 0.0277 g Kg −1 (1 w + 0 d), then rose to 0.956 ± 0.0225 g Kg −1 and 0.818 ± 0.0811 g Kg −1 (1 w + 3 d), and finally went down to 0.229 ± 0.0189 g Kg −1 and 0.551 ± 0.0132 g Kg −1 (1 w + 6 d). The total phenolic content continuously decreased in the control group HP after 2 w; the value went from 0.738 ± 0.0143 g Kg −1 to 0.216 ± 0.0100 g Kg −1 . However, the HP of the 50 µmol/L CaCl 2treated group was first up and then down from 0.410 ± 0.0262 g Kg −1 (2 w + 0 d) to 0.836 ± 0.0151 g Kg −1 (2 w + 3 d) to 0.622 ± 0.0060 g Kg −1 (2 w + 6 d). Evidently, the HP of the 50 µmol/L CaCl 2 -treated group was higher than the HP of the control group (2 w + 3 d and 2 w + 6 d). There is no obvious difference in the fruit DP.
As shown in Figure 4B,C, the PPO activity of the DP fruit was higher than the HP fruit. For HP fruit, there was no difference for both the control group and the 50 µmol/L CaCl 2 -treated group in 1w; however, the 50 µmol/L CaCl 2 -treated group was lower at any storage time. The DP fruit (1 w + 3 d and 2 w + 6 d 50 µmol/L of the CaCl 2 -treated group) was lower than the control group. Separately, the values were 0.395 ± 0.0225 Ug −1 FW and 0.391 ± 0.0218 Ug −1 FW ( Figure 4B). To summarize, PPO activity was higher in DP than in HP. POD activity was HP and far outweighed DP, whether it was the control group or the 50 µmol/L CaCl 2 -treated group ( Figure 4C). The HP of the control group had similar changing trajectories as the 50 µmol/L CaCl 2 -treated group HP after 1 w. However, initially, the control group HP was higher than the HP of the 50 µmol/L CaCl 2 -treated group but after the control group HP was lower than the 50 µmol/L CaCl 2 -treated group HP, the values (separately) were 7.85 ± 0.280 Ug −1 FW and 17.53 ± 1.075 Ug −1 FW. There was no difference with DP; however, the 50 µmol/L CaCl 2 -treated group was slightly higher after 1 w. Both the control group and the 50 µmol/L CaCl 2 -treated group experienced obviously rises. Separately, there were increases from 4.59 ± 0.223 Ug −1 FW to 15.33 ± 0.328 Ug −1 FW and from 12.64 ± 0.858 Ug −1 FW to 25.72 ± 0.170 Ug −1 FW. , and POD activity (C) in pineapple fruit during storage at 6 °C (1 and 2 w) after transfer to room temperature (3 and 6 d), and treatment with CaCl2 (50 μmol/L) or water (CK). HP: "health part", DP: "disease part". Vertical bars represent the standard error of the mean of triplicate assays. Characters at the same storage times represent significant differences (p < 0.05).

Antioxidant Enzyme Activities
SOD, CAT, and APX activity were higher in HP than DP for both the control group and the 50 μmol/L CaCl2-treated group ( Figure 5).
Regarding SOD activity, the control group HP continued rising after 1 and 2 w; the values went from 8.96 ± 0.122 Ug −1 FW to 24.90 ± 0.500 Ug −1 FW and from 15.33 ± 1.248 Ug −1 FW to 31.53 ± 0.373 Ug −1 FW ( Figure 5A). However, for the 50 μmol/L CaCl2-treated group, HP first up and then down for 1 w + 3 d and 2 w + 3 d, reaching the highest point separately at 30.33 ± 2.698 Ug −1 FW and 30.87 ± 0.310 Ug −1 FW. HP was higher than DP in both the control group and the 50 μmol/L CaCl2-treated group at any storage time.
Regarding CAT activity, HP for the control fruit (1 w + 0 d) was higher than the 50 μmol/L CaCl2-treated group value, which was 4.08 ± 0.082 Ug −1 FW ( Figure 5B). For the rest of the storage time, the GP of the 50 μmol/L CaCl2-treated group was far more than the control group HP except at 1 w + 3 d. There is no obvious change in DP, but 2 w + 6 d for the 50 μmol/L CaCl2-treated group was higher than the control group; the value was 0.810 ± 0.0120 Ug −1 FW. HP was higher than DP both in the control group and the 50 μmol/L CaCl2-treated group at any storage time.
The APX activity HP 50 μmol/L of the CaCl2-treated group was significantly higher , and POD activity (C) in pineapple fruit during storage at 6 • C (1 and 2 w) after transfer to room temperature (3 and 6 d), and treatment with CaCl 2 (50 µmol/L) or water (CK). HP: "health part", DP: "disease part". Vertical bars represent the standard error of the mean of triplicate assays. Characters at the same storage times represent significant differences (p < 0.05).

Antioxidant Enzyme Activities
SOD, CAT, and APX activity were higher in HP than DP for both the control group and the 50 µmol/L CaCl 2 -treated group ( Figure 5).
Regarding SOD activity, the control group HP continued rising after 1 and 2 w; the values went from 8.96 ± 0.122 Ug −1 FW to 24.90 ± 0.500 Ug −1 FW and from 15.33 ± 1.248 Ug −1 FW to 31.53 ± 0.373 Ug −1 FW ( Figure 5A). However, for the 50 µmol/L CaCl 2 -treated group, HP first up and then down for 1 w + 3 d and 2 w + 3 d, reaching the highest point separately at 30.33 ± 2.698 Ug −1 FW and 30.87 ± 0.310 Ug −1 FW. HP was higher than DP in both the control group and the 50 µmol/L CaCl 2 -treated group at any storage time.
Regarding CAT activity, HP for the control fruit (1 w + 0 d) was higher than the 50 µmol/L CaCl 2 -treated group value, which was 4.08 ± 0.082 Ug −1 FW ( Figure 5B). For the rest of the storage time, the GP of the 50 µmol/L CaCl 2 -treated group was far more than the control group HP except at 1 w + 3 d. There is no obvious change in DP, but 2 w + 6 d for the 50 µmol/L CaCl 2 -treated group was higher than the control group; the value was 0.810 ± 0.0120 Ug −1 FW. HP was higher than DP both in the control group and the 50 µmol/L CaCl 2 -treated group at any storage time. The APX activity HP 50 µmol/L of the CaCl 2 -treated group was significantly higher than the control group, except at 2 w + 3 d and 2 w + 6 d ( Figure 5C). The DP of the 50 µmol/L CaCl 2 -treated group was far more than the control group. The control group HP continuously increased to the highest point: 50.67 ± 0.333 Ug −1 FW. The HP of the 50 µmol/L CaCl 2 -treated group rose roughly, except at 2 w + 3 d to the highest point 46.36 ± 0.337 Ug −1 FW. HP was higher than DP in both the control group and the 50 µmol/L CaCl 2 -treated group at any storage time.
μmol/L CaCl2-treated group rose roughly, except at 2 w + 3 d to the highest point 46.36 ± 0.337 Ug −1 FW. HP was higher than DP in both the control group and the 50 μmol/L CaCl2treated group at any storage time.
Regarding the GR activity, the control group HP increased; however, the HP of the 50 μmol/L CaCl2-treated group initially decreased and then dramatically increased to 145.18 ± 0.000 UKg −1 FW after 1 w ( Figure 5D). The control group HP reduced rapidly; the value went from 125.36 ± 3.286 UKg −1 FW to 1.18 ± 0.005 UKg −1 FW. However, the HP of the 50 μmol/L CaCl2-treated group initially increased to the highest point at 92.90 ± 1.390 UKg −1 FW and finally decreased after 2 w. In addition, 2 w + 3 d was obviously higher than the control fruit. Moreover, the DP of the 50 μmol/L CaCl2-treated group was far more than the control group DP, where at 1 w + 6 d it was 164.90 ± 2.028 UKg −1 FW. Regarding DP, there were no obvious changes after 2 w. HP was higher than DP in both the control group and the 50 μmol/L CaCl2-treated group at any storage time. Regarding the GR activity, the control group HP increased; however, the HP of the 50 µmol/L CaCl 2 -treated group initially decreased and then dramatically increased to 145.18 ± 0.000 UKg −1 FW after 1 w ( Figure 5D). The control group HP reduced rapidly; the value went from 125.36 ± 3.286 UKg −1 FW to 1.18 ± 0.005 UKg −1 FW. However, the HP of the 50 µmol/L CaCl 2 -treated group initially increased to the highest point at 92.90 ± 1.390 UKg −1 FW and finally decreased after 2 w. In addition, 2 w + 3 d was obviously higher than the control fruit. Moreover, the DP of the 50 µmol/L CaCl 2 -treated group was far more than the control group DP, where at 1 w + 6 d it was 164.90 ± 2.028 UKg −1 FW. Regarding DP, there were no obvious changes after 2 w. HP was higher than DP in both the control group and the 50 µmol/L CaCl 2 -treated group at any storage time. Table 1 shows the data processing results obtained from the sequencing of pineapple fruit under refrigeration (1 w) and the transfer to room temperature (3 d). The proportion of bases with a quality of no less than 20 after filtration (Q20) was 96-98%, the proportion of bases with a quality of no less than 30 after filtration (Q30) was 91-93%, and the overall data sequencing error rate was 0.03%, and the base numbers of G and C after filtration accounted for 47-52% of the total base number. Table 1. Transcriptome raw data for each sample, respectively, raw and clean reads, raw and clean bases, error rate, percentages of bases with Phred values greater than 20 in total bases, percentages of bases with Phred values greater than 30 in total bases, percentages of G and C in four bases in clean reads.

Overall Transcriptome Quality Assessment
The principal component analysis (PCA) is shown in Figure 6A; the PCA1 is 89% variance and PCA2 is 5% variance. The CK0d group and TR0d group are indistinct distinctions. In this connection, TR0 d_vs_CK0 d did not conduct the control analysis. The TR3d group and CK3d group were dispersed and each group sample was gathered together.
As shown in Figure 6B, the Pearson correlation coefficient is more than 0.8 after three repetitions. At the same time, the correlation coefficient of the repeated sample is higher than that of the non-repeated sample, which shows that the sample has good biological repeatability, so it can be concluded that the transcriptome sequencing results are reliable. Figure 6C shows each sample's FAPKM distribution condition. for FPKM values (C). Boxplot C: Red represents"CK0d"; Orange represents"TR0d"; Blue represents"CK3d"; Green represents"TR0d".

Screening Results of DEGs
The co-expression Venn diagram ( Figure 7A) shows that the differentially expressed gene (DEG) analyses resulted in CK3 d_vs_CK0 d, TR3 d_vs_TR0 d, and TR3 d_vs_CK3 d-three comparison combinations; the total DEGs are 6094, 4167, and 1809, respectively. There were 451 (5.9%) among the number of genes co-expressed.
The volcano map ( Figure 7B  The Figure 7C histogram shows the CK3 d_vs_CK0 d, TR3 d_vs_TR0 d, and TR3 d_vs_CK3 d DEG numbers of up-and downregulated genes, respectively, as 3092 and3002,2133 and 2,033,986 and 823.

DEGs, GO, and KEGG Pathway Enrichment Analysis
Gene ontology (GO) is a comprehensive database describing gene functions; through the analysis of DEGs, we found the up-and downregulated GO enrichments most significantly in each comparison combination. From Figure 8A, one can see that the CK3 d_vs_CK0 d-downregulated GO enrichment included plant hormone signal transduction, DNA replication, starch and sucrose metabolism, and cyanoamino acid metabolism. TR3 d_vs_TR0 d-downregulated GO enrichment included plant hormone signal transduction and plant circadian rhythm. Both CK3 d_vs_CK0 d and TR3 d_vs_TR0 d-downregulated GO enrichments have plant hormone signal transductions; however, CK3 d_vs_CK0 d counts more than TR3 d_vs_TR0 d. TR3 d_vs_CK3 d-downregulated GO enrichment may include the biosynthesis of amino acids, glycolysis/gluconeogenesis, and carbon metabolism. For upregulated GO enrichment, CK3 d_vs_CK0 d and TR3 d_vs_TR0 d include the biosynthesis of amino acids, carbon metabolism, alpha-linolenic acid metabolism, tropane, piperidine, and pyridine alkaloid biosynthesis. However, TR3 d_vs_CK3 d was not significantly upregulated in GO enrichment.  MAPK signaling pathway only occurs in CK3 d_vs_CK0 d, proving that the control group can more easily produce ROS compared to the treatment group. In addition, regarding phenylpropanoid biosynthesis occurring in TR3 d_vs_TR0 d, the treatment group can produce more antioxidant bioactive compounds. A further explanation is that the treatment group can alleviate the IB of pineapple.  As presented in Figure 8B, the Kyoto Encyclopedia of Genes and Genomes (KEGG) shows three comparison combinations CK3 d_vs_CK0 d, TR3 d_vs_TR0 d, and TR3 d_vs_CK3 d (the enrichment results are shown in the figure below). It was found that the above downregulated and upregulated genes were mainly involved in plant hormone signal transduction, starch and sucrose metabolism, biosynthesis of amino acids, carbon metabolism, phenylpropanoid biosynthesis, glycolysis/gluconeogenesis (KEGG ID: obr00010), and the MAPK signaling pathway plant (KEGG ID: obr04016). MAPK signaling pathway is a stress response pathway of plants to the surrounding environment. The MAPK signaling pathway only occurs in CK3 d_vs_CK0 d, proving that the control group can more easily produce ROS compared to the treatment group. In addition, regarding phenylpropanoid biosynthesis occurring in TR3 d_vs_TR0 d, the treatment group can produce more antioxidant bioactive compounds. A further explanation is that the treatment group can alleviate the IB of pineapple.

Screening and Differential Expressions of Oxidation and Antioxidation-Related Genes
Using a heatmap to further prove the relationship between the occurrence of pineapple IB and oxidation and antioxidation systems. A total of 83 DEG-related genes were screened from transcriptome data, including POD (36) (22), respectively. The APX family (g.109727824, g.109726514), CAT family (g.109729040, g.109707674), SOD family (g.109716880), and GR (g.109727484), respectively, were upregulated after fruit disease. The PPO family (g.109719920, g.109719939) was upregulated. The GRX 4 family (g.109711896) upregulated during the treatment group (3 d). Many genes of the GST family and POD family were upregulated under fruit disease. Perhaps these are related to the occurrence of IB. These phenomena show the pineapple IB is related to oxidation and antioxidation systems. Pineapple fruit will activate its own antioxidant system to eliminate ROS, but the oxidase activity will also increase after browning. These are also consistent with the physiological results ( Figure 9).

Phenylpropanoid Biosynthesis and MAPK Pathway
The plant shikimate pathway is the entry to the biosynthesis of phenylpropanoids. The pathway of shikimic acid involves producing aromatic amino acids and other aromatic compounds (phenylpropanoids) ( Figure 10C). Seventy-four DEGs were detected in this pathway, including cytochrome (CYP), 4-coumarate-CoA ligase (4CL), putrescine hydroxycinnamoyl transferase(PHT), peroxidase (POD), β-glucosidase, cationic peroxidase (Cat-POD), caffeoyl shikimate esterase (CSE), trans-cinnamate 4-monooxygenase (T-C4M), phenylalanine ammonia-lyase-like (PAL), hydroxy cinnamoyl transferase(HCT), cinnamyl alcohol dehydrogenase (CAD), BAHD acyltransferase (BAHD-acylt), caffeoyl-CoA O-methyltransferase (CCoAOMT), aldehyde dehydrogenase (ALD), putrescine hydroxycinnamoyl transferase (PHT), vinorine synthase (VSY), salutaridinol 7-O-acetyltransferase (SOAT), and cinnamoyl-CoA reductase (CCoAR). PAL, 4CL, and HCT upregulation can lead to CAD upregulation and will increase the lignin content. During the heatmap analysis, PAL (g.109727686, g.109707567) was upregulated in the control and treatment groups at 3 d. Moreover, 4CL (g.109720819, g.109711909) and HCT (g.109703812, g.109706886) upregulation led to CAD (g.109715017) upregulation. Upregulation of these enzyme genes will increase lignin synthesis ( Figure 10D). Lignin biosynthesis ( Figure 10E) is a cross-linked phenolic polymer. Therefore, total phenol content will enhance. It can be seen from the physiological results in Figure 4 that the total phenol content at 1 w + 3 d is higher than 1 w + 0 d, which is consistent with the transcriptome data.  The MAPK pathway ( Figure 10A) is mainly related to pathogen infection, pathogen attack, phytohormones, external environmental stimulation (cold, salt, drought, osmotic stress, etc.), ozone, and wounding. Under a series of enzymatic reactions, H 2 O 2 can further cause cell death and the H 2 O 2 product. ABA can increase fruit stress adaptation and CAT can suppress the H 2 O 2 product increase regarding the fruit stress-tolerant response. Furthermore, Ca 2+ can regulate ROS homeostasis. As seen in Figure 10B, the heat map analysis of the MAPK pathway chitinase is highly expressed in HP (g.109713885). At the same time, the expression in the TR group is significantly more than in the CK group, indicating that the disease in the control group is more serious. The WRKY gene family is also related to plant resistance; g.109709450, g.10970389, and g.109719592 are highly expressed in DP. The expression of the CK group was significantly more than that of the TR group, indicating that the expressions of WRKY gene family genes would increase when plants were stimulated by the outside world.
Plant cells are affected by many factors, Figure 10F shows
The plant-pathogen interaction ( Figure 11B) includes the antiviral protein, pathogenic protein, heat shock protein, respiratory burst oxidase, and cysteine protease. Respiratory burst oxidase can produce a little of ROS in the plant body to resist external pressure. Genes (g.109722976, g.109713174) are great quantity expressions in DP, g.109717636 is expressed in HP. Cytochrome oxidase (g.109703608) promotes ATP synthesis in DP ( Figure 11E).
The tricarboxylic acid cycle (TAC) is the second stage of aerobic respiration ( Figure 11F) that shows pyruvate dehydrogenase (g.109721645, g.109728788, g.109719561), promotes acetone to produce acetyl coenzyme A, generates citric acid, and enters the citric acid cycle after fruit morbidity. During gluconeogenesis, oxaloacetic acid generates phosphoenolpyruvate and carbon dioxide under the action of phosphoenolpyruvate carboxylation kinase (g.109710186, g.109714503) in DP. TR group genes express more than the CK group.

Transcription Factor Analysis
Among all of the gene families, we selected the 15 gene families with the largest p portions (Figure 12). The six families with the largest proportions are Pkinase (2 Pkinase Tyr (134), P450 (79), Myb DNA-binding (60), NB AR (59), and PPR (46). B Pkinase and Pkinase Tyr are related to the MAPK pathway. When plants are under verse conditions, MAPKK and MAPK can be activated in turn and jointly regulate cell physiological and pathological processes, mainly receptor proteins, protein kinase, signal transductions. P450 main includes cytochrome, premnaspirodiene oxygen oryzalexin E synthase, abscisic acid 8'-hydroxylase, geraniol 8-hydroxylase, Ent-k renoic acid oxidase, alkane hydroxylase MAH, isoflavone e 2'-hydroxylase, allene ox synthase, flavonoid 3'-monooxygenase, 3,9-dihydroxypterocarpan 6A-monooxygen and trans-cinnamate 4-monooxygenase. NB_AR is mainly related to disease resista proteins. The heatmap of these gene families is shown in Figure 13.  Figure 11C,D,G show that fatty acid biosynthesis, fatty acid biodegradation, and fatty acid metabolism play important roles in the cold resistance of plants. In the process of fatty acid synthesis, acetyl coenzyme A generates malonyl coenzyme A (g.109713239) under the action of acetyl coenzyme A carboxylase (ACC) (g.109722793) in DP ( Figure 11C). Acyl carrier protein (ACP) is a domain in the fatty acid synthase (FAS) peptide chain, which is used to synthesize fatty acids. ACP (g.109723533, g.109726849, g.109714538) is upregulated in DP, but ACP (g.109713115) is in HP. Peroxisomal (g.109716412) β-oxidation mainly oxidizes ultra-long-chain fatty acids ( Figure 11D), Fatty acyl coenzyme A with a shortened chain is formed and then transported to mitochondria for complete oxidation.
Fatty acid oxidation also takes place in other organelles, such as endoplasmic reticulum ω-oxidation and the introduction of hydroxyl groups on the last carbons of fatty acids, followed by oxidation to aldehyde and carboxyl groups to form dicarboxylic acids. The enzymes involved are ω-hydroxylase (g.109721150), alcohol dehydrogenase (g.109707359, g.109721267), and aldehyde dehydrogenase (g.109728684, g.109723285), which are also involved in the metabolic process of ethanol ( Figure 11G).

Transcription Factor Analysis
Among all of the gene families, we selected the 15 gene families with the largest proportions ( Figure 12). The six families with the largest proportions are Pkinase (244), Pkinase Tyr (134), P450 (79), Myb DNA-binding (60), NB AR (59), and PPR (46). Both Pkinase and Pkinase Tyr are related to the MAPK pathway. When plants are under adverse conditions, MAPKK and MAPK can be activated in turn and jointly regulate cellular physiological and pathological processes, mainly receptor proteins, protein kinase, and signal transductions. P450 main includes cytochrome, premnaspirodiene oxygenase, oryzalexin E synthase, abscisic acid 8 -hydroxylase, geraniol 8-hydroxylase, Ent-kaurenoic acid oxidase, alkane hydroxylase MAH, isoflavone e 2 -hydroxylase, allene oxide synthase, flavonoid 3 -monooxygenase, 3,9-dihydroxypterocarpan 6A-monooxygenase, and trans-cinnamate 4-monooxygenase. NB_AR is mainly related to disease resistance proteins. The heatmap of these gene families is shown in Figure 13.  Calmodulin-binding receptor-like cytoplasmic kinase 2 (CRCK2) can against various diseases. From Figure 13A,B, Calmodulin-binding receptor-like cytoplasmic kinase 2 (g.109712591, g.109728905, g.109727024, and g.109725844) is upregulated in DP (3 d), but g.109727204 and g.109718861 are upregulated in HP (0 d). Calcium-dependent protein kinase (g.109727237, g.109714958, g.109715490, g.109715759) is upregulated in DP (3 d). Compared with the CK group, the FPKM of the TR group are lower. It shows that CaCl 2 is effective in fruit disease. Flavonoid 3'-monooxygenase-like (F3'H) belongs to the cytochrome P450 subfamily, which catalyzes NADPH and the O 2 -dependent monooxygenation. Anthocyanins and proanthocyanidins can be produced, which can prevent cells from being damaged by strong UV-B radiation. Flavonoid 3'-hydroxylase can also participate in the acetic malonic acid pathway through flavanone substances to protect plants from animals, plants, pathogenic microorganisms, etc. F3'H (g.109720891, g.109721770) is upregulated in HP (0 d) ( Figure 13C). It shows that the content of anthocyanins and other antioxidant substances decreases after fruit disease. Allene oxide synthase (AOS) generates precursors of the defense hormone jasmonate (JA). AOS (g.109720968) is upregulated in DP (3 d). It is proven that a large number of resistance hormones will be synthesized after fruit disease. The MYB TF family refers to a class of transcription factors containing the MYB domain, which is a peptide of about 51-52 amino acids and contains a great deal of highly conserved amino acid residues and spacer sequences. The first involves tryptophan residues with regular intervals of about 18 amino acids, which participate in the formation of hydrophobic cores in the spatial structure ( Figure 13D). From (Figure 13E), the largest proportion of the resistance protein is upregulated in HP (0 d). However, the TR group DP FPKM is further than the CK group. The pentatricopeptide repeat-containing protein is the largest protein family in plants for plant growth development ( Figure 13F). Calmodulin-binding receptor-like cytoplasmic kinase 2 (CRCK2) can against various diseases. From Figure 13A,B, Calmodulin-binding receptor-like cytoplasmic kinase 2 (g.109712591, g.109728905, g.109727024, and g.109725844) is upregulated in DP (3 d), but g.109727204 and g.109718861 are upregulated in HP (0 d). Calcium-dependent protein kinase (g.109727237, g.109714958, g.109715490, g.109715759) is upregulated in DP (3 d). Compared with the CK group, the FPKM of the TR group are lower. It shows that CaCl2 is effective in fruit disease. Flavonoid 3'-monooxygenase-like (F3'H) belongs to the cytochrome P450 subfamily, which catalyzes NADPH and the O2-dependent monooxygenation. Anthocyanins and proanthocyanidins can be produced, which can prevent cells from being damaged by strong UV-B radiation. Flavonoid 3'-hydroxylase can also participate in the acetic malonic acid pathway through flavanone substances to protect plants from animals, plants, pathogenic microorganisms, etc. F3'H (g.109720891, g.109721770) is upregulated in HP (0 d) ( Figure 13C). It shows that the content of anthocyanins and other antioxidant substances decreases after fruit disease. Allene oxide synthase (AOS) generates precursors of the defense hormone jasmonate (JA). AOS (g.109720968) is upregulated in DP (3 d). It is proven that a large number of resistance hormones will be synthesized after fruit disease. The MYB TF family refers to a class of transcription factors containing the MYB domain, which is a peptide of about 51-52 amino acids and contains a great deal of highly conserved amino acid residues and spacer sequences. The first involves tryptophan residues with regular intervals of about 18 amino acids, which participate in the

Discussion
The pineapple, which is a tropical or sub-tropical fruit, is usually transported by a cold chain process [23]. Refrigeration is a general method used to store fruits and vegetables, but cold stress can induce the chilling injury of pineapple fruit, which principally restricts the postharvest fruit quality and shelf life of pineapples [24]. To date, plenty of physical and chemical tests have shown good progress in delaying pineapple CI. However, it is difficult to bring about massive applications for economic costs, the stability of efficacy, and operation accessibility; therefore, developing more practical and reliable approaches to ameliorate the storability of pineapple fruit is (as usual) required.
Ca, as a second messenger, is important in plant cold stress response. Studies have shown that CaCl 2 treatment alleviated bermudagrass chilling injury by regulating reactive oxygen species (ROS) of products and lessening cell damage [25]. In addition, whether it involves pre-and postharvest foliar application calcium, the effects of delayed aging or ripening and improving physiological disorders are apparent in fruits and vegetables [26]. Ca can improve apple fruit firmness and reduce fruit decay [27]. CaCl 2 can reduce the occurrence of the internal browning of pineapple fruit under low temperatures [28]. Our present study also obtained similar results.
In the present study, pineapples were processed via exogenous CaCl 2 (50, 100, and 150 µmol/L) under chilling temperatures at different times; they were transferred to a normal temperature to simulate the sale of pineapple fruit in actual production. We investigated what affects the chilling temperature, inducing the CI development of postharvest pineapple via exogenous CaCl 2 treatment and no treatment. To screen for the best CaCl 2 concentrations, we measured some characterization indicators and related indexes of cell membrane destruction, such as the incidence rate, IB index, and electrical conductivity. A large amount of research has confirmed that CaCl 2 can alleviate CI (the incidence rate, IB index, and electrical conductivity) [29]. In our present study, after the transfer from refrigeration to room temperature, the incidence rate, IB index, and electrical conductivity increased, indicating that their membranes had become damaged. Previous studies have proved that CaCl 2 treatment can alleviate CI symptoms. Meanwhile, along with the extension of storage time, the lower the temperature, the more serious the chilling injury [30]. In our study, it was not hard to find similar results. It was proved (again) that calcium can enhance cold tolerance in pineapple fruit.
Low-temperature induction of reactive oxygen species (ROS), such as O 2 • − and H 2 O 2 , in turn, caused cell damage and resulted in chilling symptoms [31]. In our study, certain similarities changed between the ROS (O 2 • − and H 2 O 2 ) and MDA of fruit DP, showing the accumulation of lipid peroxide and enhanced membrane permeability after fruit disease [32].
The ascorbate-glutathione (AsA-GSH) cycle plays an important role in the antioxidant pathway, which can directly scavenge the H 2 O 2 produced. Moreover, ASA and GSH are the most abundant low molecular non-enzymatic antioxidants and GR is a key enzyme in the AsA-GSH cycle [33]. Some studies have shown that CaCl 2 treatment alleviated CI by enhancing the AsA-GSH cycle and antioxidant enzyme activities to quench ROS in loquat fruit [34]. In the present study, AsA-GSH content and APX and GR activity distinctly decreased in DP, indicating that antioxidant capacity is not enough to resist external cold stress, leading to the IB of pineapple fruit.
During storage, pineapple fruit CI browning symptoms were a result of a physiological disorder, which is usually due to oxidation of phenolic compounds and accumulation of polyphenol oxidase [35]. In the present study, there was PPO activity contrary to the POD activity of HP and DP; however, the POD activity changing tendency was similar to other antioxidant enzymes, which were enhanced during cold storage and decreased strongly at room temperature [36]. This indicates that POD is a two-way street enzyme, which may serve as both oxidase and antioxidant enzymes.
ROS was reflected in the total antioxidant capacity and APX activity [32]. ROS homeostasis is usually maintained by the balance of ROS production and scavenging in plant organisms. Plant antioxidant systems-to eliminate ROS-including enzymatic and nonenzymatic antioxidative components (to maintain ROS homeostasis) [37]. Antioxidative enzymes, such as SOD, APX, and CAT, play important roles in the ROS -scavenging pathway [37]. In our study, there was a rapid increase in the activities of the enzymatic antioxidative biomarkers (SOD, APX, CAT) by CaCl 2 treatment. Furthermore, the activities of enzymatic antioxidative biomarkers decreased in the DP of fruit. It was confirmed (on various fruits) that CaCl 2 treatment maintained higher activities of enzymatic antioxidative biomarkers [38].
Through the transcriptome data analysis, we concluded that fruit disease was related to oxidation and oxidation resistance. Phenylpropanoid biosynthesis is related to the synthesis of phenolic substances. Total phenols, flavonoids, and anthocyanins are antioxidant bioactive compounds [39]. In this study, the upregulation of genes related to lignin synthesis led to the increase in total phenol content, which is consistent with the physiological results. The MAPK pathway is a stress response of plants to the outside world and is related to ROS. It has been proven that ROS can induce MAPK activation; meanwhile, MAPK cascades can modulate ROS production and responses [40]. By the transcription factor analysis, we can assume that the two largest families are related to MAPK and the plant disease resistance protein. It may be mainly related to the occurrence of the IB of pineapple.

Plant Material and Treatments
Mature green stage pineapple (Ananas comosus L. Merr.) was harvested from a local orchard in Qionghai city, Hainan Province, and transferred to a laboratory at the Hainan University Food Science and Engineering building within 3 h. A total of 600 fruits were randomly selected, with 150 fruits in each group (uniform size and color). All fruits were cleaned, disinfected, and then dried for the next step. Three groups were used, separately immersed in 50, 100, and 150 µmol/L of CaCl 2 solution in 100 L of deionized water for 5 min at 25 • C. One group was dipped in deionized water for 5 min at 25 • C and was served as the control. After drying with natural air, fruits were packed in a plastic box with plastic wrap for three hours to ensure full absorption. Then, all of the fruits were transferred to cold storage (6 ± 1 • C and 85 ± 5% relative humidity (RH)) for subsequent experiments.
After cold storage for 1 w, we divided all processed fruits into two parts-one part was transferred to normal temperature (25 • C ± 1 • C and 85 ± 5% relative humidity (RH)) and the remainder of the fruits did not move from the cold storage for 2 w. The pulps near the cores of the pineapples were taken for experimental study. Fifteen fruits were sampled from each treatment after cold storage (1 w) and were transferred to room temperature at 0, 3, and 6 d. There were three replicates included in each treatment. Similarly, the fruit samples obtained after cold storage (2 w) were transferred to room temperature at 0, 3, and 6 d. According to the statistical results of internal browning, three pineapple fruits with similar statistical conditions were selected from each group, and the flesh "F/C (Flesh/core)" parts were taken, respectively. One part was used to determine the physiological indexes, and the other part was immediately mixed with liquid nitrogen and stored in a −80 • C ultra-low-temperature refrigerator for reserve. Meanwhile, all of the samples (ground to powder under liquid nitrogen) were stored at −80 • C for subsequent experiments.

IB Appearance
The fruit was cut in half and the IB appearance was photograph-recorded during each sampling.

Electrical Conductivity Measurement
Relative conductivity involved the eigenvalue index reflecting membrane permeability. Methods refer to [41] and were slightly modified. Nine fruits were cut in half-half were used for conductivity determination and half for sampling. A hole punch (6 mm in diameter) was used randomly, with 5 little discs picked from the pulp near the heart of the fruit, for a total of 15 little discs. Then they were cleaned 3 times with distilled water and rocked for half an hour in the water bath. The weight and temperature remained unchanged before and after boiling for 20 min. Meanwhile, the initial (E 0 ) and final (E t ) values were measured using a conductivity meter before and after boiling. The formula is as follows: Relative conductivity (%) = E0 Et × 100 %

Assays of ROS Production and Membrane Lipid Peroxidation MDA Contents
Malondialdehyde (MDA) is one of the membrane lipid peroxidation products determined by the thiobarbituric acid assay [41], slightly modified. Accumulation caused damage to the cytoplasmic membrane and organelles of fruit thinning cells. Moreover, 0.1 g of flesh was added to 1 mL of 100 g/L of TCA and centrifuged at 10,000× g at 4 • C for 20 min. We took the supernatant and added 2.0 mL of 0.67%TBA boiled in boiling water for 20 min, cooling the centrifugation again. Then the absorbance value was measured at 450,532,600 nm. The MDA content was calculated by formula.
The ROS production was assessed by the Superoxide Anion Radical (O 2 • − ) generation rate and the hydrogen peroxide (H 2 O 2 ) content was determined by using kits from Suzhou Grace (China).

Determination of Total Phenolic Content
Polyphenols play important roles for plants because of their antioxidant properties. We followed the method in [42], slightly modified. We weighed 1 g of the powdered sample and added it to 80% pre-cooled methanol solution at 4 • C for 20 min after oscillation blending. We took 500 µL of extract, added 4 mL of pre-cooled methanol solution, and 500 µL of Folin phenol, and let it stand for 5 min. Then we added 5 mL of 20% sodium carbonate solution in darkness at 25 • C for 1 h. The total phenol content was calculated according to the gallic acid standard curve.
PPO and POD were determined by a UV-5000PC spectrophotometer, METASH, China. Respectively, absorbances were at 420 and 470 nm. The cuvette was added to 1.6 mL of 50 mM of phosphate buffer (pH 5.5), 0.4 mL of 10 mM pyrocatechol, and then quickly added 0.2 mL of enzyme extract solution, blending the measure absorbance change curve. One unit (U) of PPO activity was defined as the amount of the enzyme that caused a change in A 420 of 0.01 min −1 . The cuvette was added to 2 mL of 25 mM of guaiacol, 0.2 mL of H 2 O 2 , and 0.1 mL of an enzyme extract to measure POD activity. One U of POD activity was defined as the amount of the enzyme that caused an increase in the A 470 of min −1 .

Assessments of Ascorbic Acid (ASA) and Glutathione (GSH) Contents
ASA and GSH play pivotal roles in the AsA-GSH cycle. Glutathione regenerates ASA in the AsA-GSH cycle and regulates H 2 O 2 [33]. The ASA content was determined by the red phenanthroline method [45], slightly modified. We took 0.1 g of pulp and added 1 mL of TCA solution, extracted for 10 min, and centrifuged at 12,000× g at 4 • C for 30 min. Then we took 0.5 mL of the extract, added 50 g/LTCA and 1 mL of anhydrous ethanol in turn, and mixed. Finally, 0.5 mL of 0.4% phosphoric acid-ethanol solution, 1 mL of 5 g/LBPethanol solution, and 0.5 mL of 0.3 g/L FeCl 3 -ethanol solution were successively added. The absorbance value was measured at 534 nm after 60 min of a water bath at 30 • C. The content of the ascorbic acid was calculated by the standard curve.
The content of GSH was determined by the dithio-dinitrobenzoic (DTNB) acid method [46], slightly modified. We added 0.1 g of flesh to 1 mL of precooled 50 g/L trichloroacetic acid solution (containing 5 mmol/EDTA-Na2) and centrifuged at 12,000× g at 4 • C for 20 min. We took two test tubes and added 0.5 mL of supernatant and 1.0 mL of 0.1 mol/L pH 7.7 phosphoric acid buffer solution, respectively. Then, to one tube, we added 0.5 mL of 4 mmol/LDTNB solution; to another tube, we added 0.5 mL of 0.1 mol/L pH 6.8 phosphoric acid buffer solution. Two test tubes were placed at 25 • C (thermal reaction) for 10 min. The content of the GSH was calculated by the standard curve.

Assessment of Antioxidant Enzyme Activities
Enzyme activity was conducted by the kinetic analysis of the absorbance value. The determination of GR [47]: 0.1 g of flesh was added to a 5.0 mL, 4 • C precooled, 0.1 mol/L pH 7.5 phosphoric acid buffer solution (containing 1 mmol/EDTA-Na2) and centrifuged at 12,000× g at 4 • C for 20 min. Then, we added a 0.9 mL 0.1 mol/L pH 7.5 phosphoric acid buffer solution (containing 1 mmol/EDTA), 0.34 mL 5 mmol/L GSSG solution, and 0.067 mL enzyme extract solution. Finally, 13.33 µL of 4 mmol/L NADPH solution was added to initiate the enzyme reaction. We immediately measured the dynamic OD value every 5 s for 35 s at 340 nm. The enzyme activity was calculated according to the formula. The determination of APX refers to [48], slightly modified. Flesh powder (0.1 g) was used to analyze each APX activity. The APX activity was presented by U kg −1 FW.
Catalase activity (CAT) and superoxide dismutase (SOD) were determined by using kits from Suzhou Grace (China).

Transcriptome Data Quality Assessment
The transcriptome data quality were evaluated by analyzing the transcriptome raw data, principal component analysis (PCA), pearson correlation coefficient and boxplot, the transcriptome data quality Respectively, raw data was assessed by Er%, Q20%, Q30%, and GC% in the proportion of the total base number. The principal component analysis (PCA), Pearson correlation coefficient, and FPKM value boxplot that were analyzed show the correlation between samples.

GO and KEGG Enrichment Analysis
GO is a database of gene function information; it is divided into three types: biological process (BP), cellular component (CC), and molecular function (MF), and uses less than 0.05 padj as the threshold of significant enrichment.
KEGG is a database related to metabolic pathways and uses less than 0.05 padj as the threshold of significant enrichment.

Screening of Related DEGs, Analysis of Related Pathways, and Transcription Factor Analysis
We screened the genes related to physiological indicators through GO and KEGG annotation analysis and then analyzed the expressions of the screened genes in each sample by a heatmap. We found the metabolic pathway related to this experiment by consulting, analyzing the genes in the pathway by the heat map, excavating the material change information in the metabolic pathway, and obtaining the conclusions related to this experiment. We analyzed the transcription factors of all genes, counted the top 15 transcription factors with the largest proportions, and made a pie chart. In addition, we made a heat map analysis of the genes of the largest six transcription factor families.

Conclusions
In the present study, pineapple IB was related to ROS homeostasis, which included the organism's relative balance between enzymatic and non-enzymatic antioxidants and oxides. Physiological data and transcriptome data also support this triangulation. The antioxidant capacities (SOD, POD, CAT, APX, and GR) decreased and the oxidation capacity (PPO) increased after fruit browning. Moreover, pineapple fruit treated with CaCl 2 improved this 'out-of-balance'. However, this change is difficult to understand through physiological indicators. Therefore, we analyzed it at the molecular level through the transcriptome. With pineapple IB, some metabolic pathway genes will be regulated, including phenylpropanoid biosynthesis and the MAPK pathway. These transcriptome analyses provide us with many relevant gene networks, which are very helpful for us to solve pineapple IB at the molecular level by genetic engineering.