Noni fruit’s water spot appearance on the second day of harvest: a trade-off between resistance and energy

Noni (Morinda citrifolia Linn.) is a tropical tree that bears climacteric fruit. Previous observations and research have shown that the second day (2 d) after harvest is the most important demarcation point when the fruit has the same appearance as the freshly picked fruit (0 d); however, they are beginning to become water spot appearance. We performed a conjoint analysis of metabolome and transcriptome data for noni fruit of 0 d and 2 d to reveal what happened to the fruit at the molecular level. Genes and metabolites were annotated to KEGG pathways and the co-annotated KEGG pathways were used as a statistical analysis. We found 25 pathways that were significantly altered at both metabolic and transcriptional levels, including a total of 285 differentially expressed genes (DEGs) and 11 differential metabolites through an integrative analysis of transcriptomics and metabolomics. The energy metabolism and pathways originating from phenylalanine were disturbed the most. The upregulated resistance metabolites and genes implied the increase of resistance and energy consumption in the postharvest noni fruit. Most genes involved in glycolysis were downregulated, further limiting the available energy. This lack of energy led noni fruit to water spot appearance, a prelude to softening. The metabolites and genes related to the resistance and energy interacted and restricted each other to keep noni fruit seemingly hard within two days after harvest, but actually the softening was already unstoppable. This study provides a new insight into the relationship between the metabolites and genes of noni fruit, as well as a foundation for further clarification of the post-ripening mechanism in noni fruit.

digestion, reduction in allergy and asthma symptoms, lower pain, fewer infections, and better sleep [7]. More and more studies have provided information on phytochemical compositions of noni fruit, finding variations according to regions, picking time, storage life, and processing methods [8]. Noni fruit is climacteric and easily ripens after harvest, which is not conducive to long-term storage and long-distance transportation and marketing, limiting its wider promotion to some extent. As a result, it is of great significance to study the softening mechanism of noni fruit after harvest.
Noni fruit is fleshy, greenish-white to yellowish-white and hard in its initial stage of maturity. At full maturity, the fruit becomes gray-white and soft, when the fruit is at an advanced stage of ripening, it has a unique odor (some people think it is unpleasant, but others find it acceptable) [9]. A single noni fruit can contain up to 260 seeds [10], which can form a framework which prevents the fruit from collapsing quickly, even if it softens. Softening as an indication of ripening is, therefore, often disregarded. Before noni fruit becomes gray-white, water spot appearance around the inlaid seeds appears on the fruit surface at the beginning of the initiation of softening. Some reported mechanisms of softening in postharvest fruit involve cell wall metabolism, starch metabolism, energy metabolism, hormone regulation, and pathogen infection [11][12][13]. The biggest challenge of postharvest fruit is the proliferation of fungi and bacteria, because of their ability to grow under storage conditions, thereby, damaging fruit directly [14]. When attacked by pathogens or fungi, postharvest fruit can activate their defense systems to improve their resistance [15]. By activating resistance mechanisms, plants prioritize the use of energy for their survival against a threat, which drives changes in the allocation of energy resources within the fruit.
Although some research progress has been made on the softening mechanism in some fruits [12][13][14], there are very few reports on how the fruit enters the softening stage, or the molecular mechanism regulating fruit softening in harvested noni. In a previous study, Wu et al., using liquid chromatography-mass spectrometry (LC-MS), showed that the metabolites within noni fruit vary at different ripening stages (0, 2, 4, and 6 days after harvest) when the fruit is in postharvest storage [16]. The 2nd day has been identified as the key point when noni fruit appears similar to that of 0-day fruit, but in the subsequent day, it quickly softens and smells. Furthermore, RNA-seq was performed using 0-and 2-day postharvest noni fruit, and conducted verification of the qRT-PCR. The results formed the bases for the present research.

Plant material
The noni plants were grown on a farm in Yuanjiang, Yunnan province, China, with an average elevation of approximately 350 m, average temperature of 25 °C, and an average annual precipitation of 560 mm. When the noni fruit is greenish-white in color, slightly hard, and has a relatively smooth surface, it has reached maturity and is at optimal harvesting time (Fig. 1a). When the fruit is cut in half vertically, the flesh is fresh and white (Fig. 1b). This is the 0-day (0 d) sample used in the experiment. After 2 days (48 h) of storage under 25 °C and 50-70% humidity, the color and the hardness of noni fruit seem similar with that of 0 d fruit, but on close observation, some water spot appearance can be found near the seed insertion on the surface (Fig. 1c). Cut vertically, the fruit shows a slightly blackened flesh, which is somewhat water spot appearance (Fig. 1d). This is the 2nd day (2 d) sample. After this point, the water spot increases rapidly and the fruit begins to smell until the whole fruit is soft (Figs. 1E, 1F), and it is difficult to extract high-quality RNA. The second day is the critical point of noni fruit softening, and the samples of 2 d vs. 0 d were, therefore, used for the conjoint transcriptome and metabolome analysis.
A biological replicate consisted of five randomly selected healthy noni plants, and three to five fruit were collected from each plant. Three biological replicates were conducted for each sample. The freshly picked fruit was sliced, quick-frozen with liquid nitrogen, and stored at -80 °C. Every replicate was divided into two sections, one was used for metabolomics analysis, while the other was used to extract RNA for RNA-seq.

Data resource
The RNA was extracted by TRIzol ® Reagent (Invitrogen), and mRNA was separated using Magnetic Separation Rack (Invitrogen). The platform was the Illumina Hiseq 2000. The sequencing method was Single end (SE), and the average fragments obtained were 1032 bp. The transcriptome dataset was submitted to the database of National Center for Biotechnology Information (NCBI), with an accession number of PRJNA503490. A total of three noni fruit replicates at 0 days (0 d) and three noni fruit replicates at 2nd day after harvest (2 d) were included in PRJNA503490. In transcriptomics analysis, the software were SeqPrep (https ://githu b.com/jstjo hn/ SeqPr ep) and Sickle (https ://githu b.com/najos hi/sickl e) for quality control of sequencing data, Trinity (http:// trini tyrna seq.sourc eforg e.net/, version number: trinityrnaseq-r2013-02-25) for sequence assembly and ORF prediction, and RSEM (http://dewey lab.biost at.wisc. edu/rsem/) for gene expression analysis. The metabolite profilings of 0 d and 2 d fruit using the liquid chromatography-mass spectrometry (LC-MS) method were detailed and the metabolome dataset was obtained from the study of Wu et al. [16].

Correlation analysis between transcriptome and metabolome
To understand the relationship between genes and metabolites in the noni fruit of 2 d relative to 0 d, a correlation analysis was carried out for gene expression levels and metabolites using Spearman's Rank Correlation Coefficient in SPSS 24.0 (SPSS Inc., USA). The genes and metabolites were selected when the correlation was ≥ 0.9 or ≤ −0.9. The correlation analysis was performed between DEGs (differentially expressed genes) in the transcriptome and differential metabolites in the metabolome. The genes and metabolites were annotated to KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways and the co-annotated KEGG pathways were used as a statistical analysis. The same statistical analysis was conducted on significantly different genes and metabolites.

pathways were significantly varied at both metabolic and transcriptional levels
There were 25 pathways significantly varied at both metabolic and transcriptional levels (Additional file 1: Table S1), including a total of 285 DEGs and 11 differential metabolites. To better understand these pathways, they were classified into primary and secondary metabolic pathways. Among them, the primary metabolic pathways included three pathways related to carbohydrate and energy metabolism and seven amino acid-related pathways. Secondary metabolic pathways included three pathways originating from phenylalanine, three terpene-related pathways, and a glycerolipid metabolism pathway. As for other pathways, they only contained a few downregulated genes or the FPKM (Fragments Per Kilobase Million) had a value of 0 or close to 0 in one biological replicate, and no or very few metabolites with significant differences were detected; therefore, these other pathways were not analyzed in the following correlation analysis.

Disordered energy metabolism, enrichment of DEGs in galactose metabolism, and downregulation of D-Tagatose
There are three pathways related to carbohydrate and energy metabolism: carbon fixation in photosynthetic organisms, carbon metabolism, and galactose metabolism. Among these pathways, galactose metabolism enriched the most differential genes. The thirteen DEGs showed significantly upregulated and nine DEGs showed significantly downregulated (Table 1).
Central carbon metabolism, consisting of glycolysis, pentose phosphate pathway (PPP), and the tricarboxylic acid cycle (TCA), functions in the energy and cofactor regeneration as well as biosynthetic precursor supply in higher plants [17]. The DEGs in Table 1 could be targeted in the associated pathways to better understand energy metabolism. Glycolysis provides energy in the form of ATP and reducing power [18]. The genes of PGK (phosphoglycerate kinase), PPDK (pyruvate phosphate dikinase), PKC (chloroplastic pyruvate kinase), PKP (plastid pyruvate kinase), and HXK (hexokinase) were all involved in glycolysis and were all downregulated except for one copy of HXK. PGK, a major enzyme in glycolysis and the first ATP-generating step of this pathway, was downregulated by 5.56-fold (Table 1). PGK can participate in photosynthetic and glycolytic/gluconeogenic reactions at the same compartment and/or at the same time, so the function of it is of crucial importance [19]. Pyruvate kinase (PK) is responsible for catalyzing in the last step of glycolysis, including both cytosolic (PKC) and plastid (PKP) isozymes [18]. PKP is vital for providing pyruvate and ATP to multiple plastidic biosynthetic pathways, whereas PKC functions more complicatedly because it acts in a range of biosynthetic pathways and respiration [20], such as resistance to disease [21] and fruit ripening [22]. PPDK catalyzes pyruvate to phosphoenolpyruvate (PEP), consuming ATP to produce AMP in the process. HXKs can regulate glycolysis and supply substrates to several pathways and also act as glucose sensors regulating the gene expression of central cellular components and secondary metabolism [23]. In summary, the universal downregulation of these genes in glycolysis by more than fivefold on the 2nd day after harvesting of noni fruit is speculated to cause a weakening of glycolysis, which would lead to a decrease in energy availability, and then directly lead to postharvest senescence of the fruit.
In addition, one gene in TAC and one gene in PPP were significantly upregulated. MDH (malate dehydrogenase) is a key enzyme of TAC, catalyzing oxaloacetate to produce malic acid with NAD + [24]. MDH was upregulated by 4.14-fold in the 2 d sample relative to the 0 d sample (Table 1). G6PDH (glucose-6-phosphate 1-dehydrogenase) was upregulated by 3.05-fold (Table 1) and is a key enzyme limiting the PPP by reducing NAD + to NADH [25]. Thus, TAC and PPP were also affected on the 2nd Table 1 DEGs related to carbohydrate and energy metabolism 0d and 2d means fresh-picked noni fruit and the fruit of start-to-soften (store at room temperature for 2 day), respectively. And the numbers 1, 2, and 3 mean the number of replicates. The same as Table 2 No. 0d-1 0d-2 0d-3 2d-1 2d-2 2d-3 FC Description day after harvest in noni fruit, which leads to changes in the energy state via the consumption of NAD + . In galactose metabolism, there were the DEGs associated with sugar metabolism, including genes encoding RFS (galactinol-sucrose galactosyltransferase) (five DEGs) and bifunctional UGE (UDP-glucose 4-epimerase) and UXE (UDP-xylose 4-epimerase) (two DEGs), S6PDH (NADP-dependent D-sorbitol-6-phosphate dehydrogenase), and Gal (galactosidase). The five RFS genes were all upregulated by 3.37--11.20-fold in the 2 d fruit relative to the 0 d fruit (Table 1). RFS is the key enzyme catalyzing sucrose and galactinol into raffinose and is the second step of the biosynthesis of RFOs (raffinose family of oligosaccharides) [26]. Raffinose is a widespread oligosaccharide in plants and participates in the transport and storage of carbon [27]. In addition, raffinose and its precursor galactinol accumulate in plants under various abiotic stresses, so RFS could be regarded as a protector regulating gene expression by controlling reactive oxygen species (ROS) levels [28]. S6PDH, namely aldose 6-phosphate reductase, is a vital enzyme in sorbitol biosynthesis [29]. Its expression is closely related to photosynthesis, source-sink intensity and distribution of photosynthetic products, and has an obvious relationship with the plant's resistance to abiotic and biotic stresses [30]. The consistent upregulation of RFS and S6PDH genes in noni fruit are likely to increase the defense resistance and oxidation resistance of the fruit after harvesting, which consumes energy.
In response to the postharvest noni fruit, besides RFS and S6PDH, other resistance-related genes were found that involved the upregulated genes of G6PDH and AKR (aldo-keto reductase) as well as the downregulated genes of CYSC (bifunctional L-3-cyanoalanine synthase/cysteine synthase) and GolS (galactinol synthase). G6PDH is reported in response to abiotic stresses and pathogenesis and reduces the level of ROS in plants [31]. AKRs are a superfamily of enzymes participating in redox reactions, and due to their lack of substrate specificity, the genes may be involved in widespread toxicant metabolism [32]. It was reported that a deficit in tissue energy was a key factor in the aging of postharvest fruits [33]. The upregulation of these genes is beneficial for the response to biotic and abiotic stress after harvest in noni, but the upregulation of resistance genes is an energyconsuming process, and the continuous upregulation of resistance genes may lead to energy dissipation in postharvest tissues at the later stage of storage.
CYSC, a kind of cyanoalanine synthase, could degrade the toxicity of cyanide and scavenge ROS in plants [34]. GolS genes are good targets for improving the abiotic stress tolerance of higher plants [35]. The downregulated expressions of these two genes is unfavorable for detoxification and response to stress in noni fruit after harvest, but it may save energy to some extent. In a word, from the expression characteristics of these resistance genes, the genes in postharvest noni fruit might not only relate to biotic and abiotic stress but also might be involved in energy consumption, which constitutes a trade-off between the resistance response and energy deficit in the harvested noni fruit.
α-Gal was upregulated by 8.16-fold, while β-Gal was downregulated by 6.14-fold in the 2 d samples relative to the 0 d samples (Table 1). α-Gal is rich in many fruits during early softening stages [36], degrading cell galactose, changing the interaction between pectin and other polymers, decreasing cell wall integrity, and producing galactose [37]. It also functions in other responses, including ethylene release and stress responses [38]. β-Gal affects fruit softening by expanding and softening cell walls [39]. It is speculated that although the downregulation of β-Gal hinders the softening of noni fruit to some extent, the upregulation of α-Gal makes the softening inevitable, and α-Gal may also start the ethylene pathway and resistance reaction with the consumption of energy.
D-Tagatose, a peculiar hexokinase, a natural isomer of D-galactose, and the product of galactose metabolism [40], was downregulated by 2.44-fold (Additional file 1: Table S2). In the biosynthesis pathway of D-tagatose, galactose is a raw material, catalyzed by galactitol dehydrogenase (GDH), and requires NAD + as a co-factor [41]. GolS could catalyze UDP galactose to synthesize galactinol, and RFS could consume galactinol [42]. The downregulation of GolS gene by 15.79-fold (Table 1), as mentioned above, could result in the decrease of galactinol. Moreover, all five RFS genes were significantly upregulated, indicating the consumption of galactinol, which led to a further decrease of galactinol (Fig. 2). Not only was galactinol decreased, but the GDH gene that could synthesize D-Tagatose was not detected in the 2 d samples relative to the 0 d samples, which eventually led to the downregulation of D-tagatose. NAD + also decreased after harvest of noni fruit; thus, the lack of NAD + may be an important reason for the downregulation of D-tagatose. In addition, D-tagatose is antibacterial [43], so its reduction could reduce the resistance to pathogens, leading to the softening of noni fruit. To summarize, the downregulation of GolS, the upregulation of five RFSs, and a decrease in NAD + collectively caused the downregulation of D-tagatose. These changes in the noni fruit could reduce resistance reactions and result in fruit softening; however, on the other hand, it may also be a way to save energy after harvest in noni fruit.

Seven amino acid-related pathways, with the upregulation of Alanine and downregulation of Leucine and 2-oxobutanoate
Seven amino acid-related pathways changed in the noni fruit from the 0 to 2nd day after harvesting (Additional file 1: Table S1). The genes involved in these pathways underwent complex changes, including ten upregulated genes, such as homocysteine S-methyltransferase (HMT), and twelve downregulated genes (Additional file 1: Table S3). During senescence, amino acids metabolism can counteract the reduction in electron supply from the TCA cycle [44]. The decrease of amino acids during ripening after harvest is usually related to the pathways releasing aroma volatiles and ethylene [45].
Amino acids may also be acted in fruit responses to different stresses and physiological disorders [46]. Numerous amino acid pathways were altered to answer the noni fruit of the second day, and likely impacted energy metabolism, secondary pathways, resistance response, and even ethylene production in harvested noni fruit.
Alanine (Ala) was upregulated by 5.59-fold in the 2ndday noni fruit (Additional file 1: Table S2). Ala existed in alanine, aspartate and glutamate metabolism, cysteine and methionine metabolism, biosynthesis of amino acids, as well as protein digestion and absorption. Production of Ala is a vital adaptive process because Ala acts as an osmoprotectant under stress conditions and balances a quick decline in carbohydrates; moreover, the accumulation of Ala seems to be a response to hypoxia in higher plants [47]. Therefore, Ala is one of the amino acids that play a key role in noni fruit tolerance to abiotic stresses during postharvest senescence. The upregulation of Ala is speculated to mainly affect energy metabolism and resistance reaction in the postharvested noni fruit.
L-Leucine (Leu) participating in valine, leucine, and isoleucine degradation, biosynthesis of amino acids as well as protein digestion and absorption to reply to the postharvest noni fruit was downregulated by 15.81fold (Additional file 1: Table S2). Leu often decreases in many fruits, such as apples, bananas, strawberries, and tomatoes, to release volatiles, which is the property of branched-chain amino acids like Leu [48][49][50]. Volatiles often play vital roles in plants defense performance determination [51]. Thus, the decline of Leu indicates the improvement of the resistance response after harvest in noni fruit, which usually takes energy.
2-Oxobutanoate was downregulated in the 2 d samples relative to the 0 d samples (Additional file 1: Table S2). In glycine, serine and threonine metabolism, there are two threonine synthase genes: the one upregulated by 26.68-fold, and the other by 7.83-fold (Additional file 1: Table S3), suggesting threonine should accumulate within the fruit. Threonine can be catabolized to 2-oxobutanoate by threonine ammonialyase [52]. Because of the downregulation of 2-oxobutanoate and there was no significant difference in threonine ammonia-lyase gene expression in the 2 d samples relative to the 0 d samples, it was speculated that the gene did not function. In cysteine and methionine metabolism, homocysteine S-methyltransferase (HMT) can catalyze the conversion of L-homocysteine (Hcy) to L-methionine [53]. In the 2 d samples relative to the 0 d samples, three genes were encoding HMT with an upregulation of 13.69-fold, 18.37-fold, and 12.77-fold (Additional file 1: Table S3), indicating that a large amount of methionine could be accumulated. However, methionine was not found in the metabolome in the 2nd-day noni fruit, indicating that it was further metabolized. Methionine can be catalyzed into S-adenosyl-L-methionine (SAM), which is the vital methyl donor in cells and used for the biosynthesis of 2-oxobutanoate and ethylene [54]. 2-Oxobutanoate and S-adenosylhomocysteinase (SAHase) gene were both downregulated, so the branch might direct toward ethylene synthesis (Fig. 3). Ethylene can trigger ripening, particularly in climacteric fruits [55]. However, ethylene was not detected in our results, which could be due to the use of LC-MS rather than GC-MC (Gas chromatography-mass spectrometry), or because the related intermediate metabolites of ethylene were not detected. Intermediate metabolites of ethylene may not have been detected because large-scale synthesis of ethylene had not yet occurred on the second day after noni harvest, which is consistent with the noni fruit looking similar on day 2 vs. day 0, with no significant softening. In short, based on our result of transcriptome and metabolomics, 2-oxobutanoate was downregulated not only because the threonine ammonia-lyase gene did not work but also because the synthesis route was biased to the ethylene pathway to have provoked ripening.

In pathways originating from Phe, β-glu gene is the most important, with upregulation of coumarin and vanillin
Three pathways originated from phenylalanine (Phe): phenylpropanoid biosynthesis, phenylalanine metabolism, and flavonoid biosynthesis. Phenylpropanoid biosynthesis was disturbed the most. Table 2 lists the DEGs Fig. 3 The diagrammatic drawing of 2-oxobutanoate and ethylene in cysteine and methionine metabolism. The upregulated gene in dark red solid rectangles and downregulated gene in blue solid rectangles; the upregulated metabolites in pink solid oval and downregulated metabolite in green solid oval. Ethylene was marked in yellow that are significantly up-and down-regulated. There were 11 β-glu genes showed significant changes, especially the upregulations of three β-glu genes by 805.58fold, 21.67-fold, and 31.97-fold (Table 2), suggesting their importance. β-Glu belongs to the glucohydrolase enzyme family 1 and functions in detoxification, ripening, and defense [56]. The deposition of active β-glu in the cell wall during fruit ripening and its ability to hydrolase complex glycans could be indicative of a role in degrading cell wall components to trigger fruit soften. The alterations in β-glu activity and its products of enzymatic hydrolysis are firmly associated with product quality [57]. Under some unfavorable conditions, β-glu could cleave the glucoside to synthesize glucose and coumaric acid, which is spontaneously cyclized to coumarin [58]. In the 2nd-day noni fruit, the upregulation of coumarin by 1.69-fold (Additional file 1: Table S2) was likely due to the large expressions of β-glu genes (Fig. 4).
Coumarin is a compound with a simple structural unit and is found in numerous plants, including in noni. It can interact with diverse enzymes and receptors in living organisms [59]. Coumarin has many important properties, including acting as a plant growth regulator, and providing antioxidant, antibacterial, antifungal, analgesic, anti-inflammatory, and antimutagenic properties [60], and it could also be used as biological pesticide [61]. The upregulation of coumarin in noni fruit after harvest on the second day was positive for the disease resistance and storage tolerance.
Our results showed the genes of PAL (phenylalanine ammonia-lyase), COMT (caffeic acid 3-O-methyltransferase), CAD (cinnamyl alcohol dehydrogenase), and POX (peroxidase) were changed a lot in phenylalanine metabolism, for example, PAL upregulated by 29.63fold and POX upregulated by 7.29-fold (Table 2), indicating the pathway was disturbed most. The pathway is closely linked to fruit ripening processes [62], and these key enzymes are also associated with lignin synthesis [63]. Therefore, the ripening and softness of noni fruit after harvest might be heavily influenced by phenylalanine metabolism.
Vanillin is one of the secondary metabolites of phenylalanine metabolism, and numerous key enzymes in the pathway possibly act in the catalytic reactions to synthesize vanillin, including PAL [64] and COMT [65]. The significant changes in genes of PAL, COMT, and CCoAOMT in the upstream of vanillin suggested they might participate in the biosynthetic pathway of vanillin (Fig. 4). Vanillin, the established active and representative ingredient in noni fruit [66], was upregulated by 1.77-fold (Additional file 1: Table S2). It possesses an intensely sweet and very tenacious creamy vanillalike odor [67], and could be used as a nutraceutical, antioxidant, antimicrobial and ripening agent [68,69]. Because of the antioxidant and antimicrobial properties, vanillin could be used as a food preservative [70]. Consequently, the increase of vanillin in postharvest noni fruit helped in the fruit's resistance to bacteria and oxidation and promoted ripening. Moreover, vanillin is found in a conjugated form as β-D-glucoside, which could be catalyzed by β-glu to free vanillin [71]. In the 2nd-day noni fruit, numerous β-glu genes were identified but not the accumulation of vanillin glucoside, so it was speculated that some β-glu genes took part in the biological process that dissociated vanillin glucoside. In other words, the accumulation of vanillin was partly due to β-glu. In addition, it was reported that vanillin could be produced from esters of coniferyl alcohol by lipoxygenase [72], and based on the downregulation of coniferyl alcohol in the metabolome in the 2 d noni samples relative to the 0 d samples, coniferyl alcohol might have been consumed to synthesize vanillin. In brief, the upregulation of vanillin in the postharvest noni fruit was owing to three reasons; firstly, the upregulation of PAL genes resulted in a lot of substrate; secondly, beta-glucosidase caused vanillin glucoside to convert to vanillin; thirdly, a part of vanillin might be formed from coniferyl alcohol. Vanillin in noni fruit was significantly upregulated on the second day after harvest and provides the material base for the increase of defense.
Contrary to the increase of coumarin and vanillin, coniferyl alcohol and ( −)-epigallocatechin (EGC) decreased (Additional file 1: Table S2). Coniferyl alcohol could be synthesized from coniferyl aldehyde catalyzed by CAD [73] with the downregulation by 2.58-fold in the 2 d samples relative to the 0 d samples ( Table 2). The CAD gene was downregulated by 7.07-fold (Table 2), which should be the direct cause of the downregulation of coniferyl alcohol (Fig. 4). The other cause could be that coniferyl alcohol was transformed into vanillin, as mentioned above. Considering that coniferyl alcohol is one Fig. 4 The diagrammatic drawing showing the key genes and metabolites in the postharvest noni fruit in the pathways originating from phenylalanine. The upregulated gene in dark red solid rectangles and downregulated gene in blue solid rectangles; the upregulated metabolites in pink solid oval and downregulated metabolite in green solid oval of the lignin precursors [74], and the coniferyl alcohol was positively correlated with hardness in fruit [63], the downregulation of coniferyl alcohol may be one of the causes for noni fruit softening after harvest.
In flavonoid biosynthesis, only one metabolite and one gene were significantly altered, downregulating by 4.30-fold (Additional file 1: Table S2) in EGC and by 19.71-fold (Table 2) in the acetyl-CoA-benzyl alcohol acetyltransferase (BEAT) gene. BEAT could catalyze p-coumaroyl-CoA in flavonoid biosynthesis [75]. Flavonoid biosynthesis starts with general phenylpropanoid metabolism and the precursor of EGC is p-coumaroyl-CoA [76]. The downregulation of the BEAT gene might be because p-coumaroyl-CoA was consumed and synthesized into vanillin since p-coumaroyl-CoA is also a precursor substrate for vanillin. In other words, when the substrate was competed for the pathway to vanillin won over the pathway to BEAT (Fig. 4). As a secondary metabolite synthesized downstream in flavonoid biosynthesis, EGC must be downregulated due to its inactive enzymatic reaction upstream. It was reported that EGC has stronger antioxidant and antibacterial effects [77]; therefore, its decrease might result in the decrease of antioxidant and antibacterial effects, which would negatively impact the fruit's resistance to pathogens. However, from the other angle, the decline of resistance might be favorable for energy retention.
In summary, coumarin and vanillin showed upregulation, which is effective at increasing resistance but consumed energy, while coniferyl alcohol and EGC showed downregulation, which may be an effort to save energy.
HPPD can catalyze the conversion of 4-hydroxyphenylpyruvate to homogentisate [80], and the homogentisate might be one of the substrates that go into terpenoid backbone biosynthesis. Given that there was no upregulated gene found in the upstream of isoprene in the terpenoid backbone biosynthesis, we speculate that the upregulated expression of the HPPD gene would produce homogentisate, which would start the terpenoid backbone biosynthesis and further accumulate isoprene (Fig. 5).
Isoprene was upregulated by 1.68-fold in the noni fruit on the second day after harvesting (Additional file 1: Table S2). It is a volatile organic compound with a petroleum-like odor. Its biosynthesis occurs through the methylerythritol 4-phosphate (MEP) pathway in terpenoid backbone biosynthesis, and only parts of plants could make isoprene through MEP pathway [81]. Nearly all the carbon producing isoprene comes directly from photosynthetic intermediates [82], and the emission of isoprene is a character of having a high metabolism in higher plants [83]. Its emission hardly reaches 2-5% of the carbon fixed by photosynthesis in healthy plants, but it often increases under stressful conditions [84]. Therefore, the upregulation of isoprene may indicate that the noni fruit on the second day of harvest was in vigorous metabolism and was consuming energy in abundance. Furthermore, isoprene is very active due to its diene with conjugated double bonds; consequently, the increase of isoprene could increase the fruit's resistance to abiotic and biotic stresses [85,86]. In summary, the upregulation of isoprene is beneficial to post harvested noni fruit responding to stress conditions; however, it also accelerates the consumption of photosynthetic products, which might be a double-edged sword.
2,3-Epoxymenaquinone was upregulated by 3.96-fold (Additional file 1: Table S2). It, is also named vitamin K2 2,3-epoxide, is a vitamin K metabolite occurring in the ubiquinone and other terpenoid-quinone biosynthesis. In the pathway, vitamin K2 (menaquinone) is in the upstream of 2,3-epoxymenaquinone. The upregulation of 2,3-epoxymenaquinone meant that vitamin K2 was consumed. There are two main forms of vitamin K, vitamin K1 found in plant foods and vitamin K2 found in animal and fermented foods [87]. Therefore, the upregulation of 2,3-epoxymenaquinone might suggest that the harvested noni fruit was beginning to ferment, which also caused the fruit flesh water spot appearance and even to be soft.
Zeatin biosynthesis could be regarded as a branch of terpenes, and there was one downregulated metabolite, isopentenyl adenosine (iPA), and one upregulated gene, zeatin O-glucosyltransferase (ZOG), in the pathway (Fig. 5). The multiple of downregulated iPA was the largest among the metabolites mentioned in this study, reaching 125.33 times (Additional file 1: Table S2). IPA, a member of the cytokinin family, is a product of the MEP. And isoprene is also located in MEP. Dimethylallylpp is the common precursor of iPA and isoprene, which might lead to the downregulation of iPA due to substrate competition. Not only that, but iPA is in the upstream of ZOG, and ZOG was upregulated, meaning it needed to consume iPA; therefore, iPA must be further downregulated. ZOG can catalyze zeatin to O-glucosylzeatin. The formation of O-glucosylzeatin is reversible by betaglucosidase [88], and several beta-glucosidase genes were massively upregulated as mentioned in phenylpropanoid biosynthesis. This reversibility implied the activity of cytokinin-related pathways. In summary, the upregulation of isoprene and the ZOG gene together resulted in the downregulation of iPA by a large amount.
Cytokinins inhibit fruit ripening and act in the resistance [89]. As the fruit ripens, the amount of cytokinin increases [42]. In reaction to the postharvest noni fruit, although the iPA decreased, the downstream metabolite trans-zeatin-O-glucoside synthesized by ZOG might increase. In addition, the upregulated ICME genes in terpenoid backbone biosynthesis were positive regulators of ABA signaling, suggesting the possibility of accumulation of ABA. All of these results indicated that cytokinins were involved in the post-ripening process of noni fruit. Although cytokinins can reduce the post-ripening process, not all cytokinins work at the same time, which could consume more energy and negatively affect fruit quality.

Glycerol as substrate could be oxidized to DHA in glycerolipid metabolism
Glycerolipid metabolism was disturbed in postharvest noni fruit. The metabolite dihydroxyacetone (DHA) was upregulated by 4.51-fold (Additional file 1: Table S2),  Table S5). DHA, also known as glycerone, is the simplest ketone form of ketoses and is an important intermediate in carbohydrate metabolism in higher plants during glycolysis [90]. Therefore, the upregulation of DHA is obviously closely related to energy consumption.
Glycerol as a substrate can be biosynthesized to DHA and also be catalyzed by GPAT4 to 1-acyl-sn-glycerol 3-phosphate (Additional file 1: Fig. S1). GPAT is the ratelimiting enzyme in glycerolipid synthesis and plays a pivotal role [91], and 1-acyl-sn-glycerol 3-phosphate can be further catalyzed by AGPAT5, which was reduced in the 2nd-day noni fruit, meaning that 1-acyl-sn-glycerol 3-phosphate might accumulate. Glycerolipid metabolism may participate in stress response signaling in response to the phospholipid precursor [92], so it might increase the resistance of postharvest noni fruit more or less by energy expenditure.

Conclusions
A hypothesis model was proposed by the conjoint analysis of transcriptomic and metabolomic profiles in noni fruit on the second day vs. the initial day of harvesting (Fig. 6). This model indicated that the primary metabolic pathways, including amino acid-related pathways as well as carbohydrate and energy metabolism, were disturbed. The affected secondary metabolic pathways included pathways originating from phenylalanine, terpene-related pathways, zeatin biosynthesis, as well as glycerolipid metabolism. Numerous resistance metabolites, including Ala, coumarin, vanillin, and isoprene, and resistance genes, including β-glu, PAL, POX, G6PDH, S6PDH, and AKR, were upregulated, which implied the increase of resistance in the postharvest noni fruit. As a result, a great deal of energy was consumed. Not all resistant metabolites were upregulated, such as EGC and D-tagatose, which may be trivial compared with the number that were upregulated, but could indicate a struggle between consumption and saving energy. In addition, PGK, Fig. 6 The proposed model explaining the response to the harvest in the noni fruit of 2d relative to 0d. The bold lines above or below represent pathways, and the names of the pathways are labeled above or below. The solid rectangles in blue and red are significantly altered metabolites, upregulated in red, and downregulated in blue. The yellow solid arrow simply indicates the main cause of the upregulation or downregulation of this metabolite. The colorful boxes with rounded corners in the center show several major changes, including resistance, energy, soften, and hormone. The relationship between metabolites and changes is linked by arrows PPDK, PKC, PKP, and HXK were all downregulated in glycolysis, further limiting energy. The lack of energy led noni fruit to water spot appearance, a prelude to softening, and the downregulation of coniferyl alcohol and D-tagatose directly led to softening. Upregulation of 2,3-epoxymenaquinone suggested the initiation of fermentation, which would also lead to softening. The metabolites and genes related to the resistance and energy interacted and restricted each other to keep noni fruit seemingly hard within two days after harvest, but actually the softening was already unstoppable. Moreover, the cytokinin and ethylene pathways were also subtly altering in preparation for the following acceleration of softening. Our results lay a foundation for the clarification of noni's post-ripening mechanism and could scientifically guide the optimal storage time of noni fruit.
Additional file 1: Fig. S1. The diagrammatic drawing showing the key genes and the upregulation of glycerone (DHA) in the postharvest noni fruit in glycerolipid metabolism. The upregulated genes in dark red solid rectangles and downregulated genes in blue solid rectangles; the upregulated metabolite in red solid oval. Table S1 The differential metabolites in all three comparisons. Table S2. Metabolites located in the pathways significantly altered at both metabolic and transcriptional levels. Table S3. DEGs in amino acid-related pathways. Table S4. DEGs in terpene-related pathways. Table S5. DEGs in glycerolipid metabolism.