1. Introduction
Apple (Malus × domestica Borkh) is one of the most widely cultivated fruits in the temperate climate zone, with an annual world harvest of around 96 million tonnes [1]. In the Czech Republic, more than 100,000 tonnes of apples were harvested in production orchards in 2023 (based on data from the Central Institute for Supervising and Testing in Agriculture, 2023). Although they are more expensive compared to imported apples, the regional fruits are favored by consumers as they believe that (thanks to the widespread integral farming practices of local producers) fewer pesticides have been used in their production. Unfortunately, under these conditions, dishonest traders tend to falsely declare information when placing this popular fruit on the market. Alongside the intentional false declaration of geographical origin, the other fraudulent practice is the mislabeling of the particular apple cultivar. In detecting such economically motivated frauds, the development of methods to detect them, the aim of which is both to take preventive measures against fraudulent practices and to protect consumers, is a challenging task [2,3].
To date, a variety of instrumental techniques have been published to authenticate the geographical origin and/or cultivar of apples. This includes near-infrared spectroscopy (NIR) [4], fluorescence spectroscopy [5], head-space solid-phase microextraction coupled with gas chromatography and mass spectrometry (SPME-GC-MS) [6,7], isotope ratio mass spectrometry (IR-MS), either separately [2,8,9] or in combination with elemental analysis by inductively coupled plasma mass spectrometry (ICP-MS) [10], and electronic nose and electronic tongue [3]. Further information on these authentication studies is summarized in Table 1. However, as shown here, some of them were performed with a limited number of samples or did not sufficiently cover factors of natural variability (different apple cultivars, growing locations, farming systems, different harvest years, etc.), and in some cases, the description of statistical methods for data processing was insufficient. Under these conditions, generic applicability of the results could be rather difficult.
Various parts and/or processed forms of apples were used within the studies listed in Table 1; nevertheless, only in one of them [4] was the authentication based on the data collected from the apple surface. However, the NIR technique used in that particular study did not allow for the identification of characteristic marker compounds that could be suitable for authentication via targeted analysis. In this context, it is worth noting that the cuticular layer is an interesting matrix in the search for authenticity markers. The cuticular layer contains a number of secondary metabolites whose profile is not only characteristic of the respective cultivar (genotype) but can also be influenced by various external factors in the respective growing location, such as local weather conditions, application of pesticides and growth regulators, diseases and pests. The cuticle, the outer protective layer of the fruit, consists of structural polymers coated with a layer of wax [11]. While the intracuticular waxes are directly incorporated into the cutin, the epicuticular waxes cover a surface of the cutin polymers [12]. The cuticular layer is a complex mixture of secondary metabolites, which include long-chain hydrocarbons and their derivatives, such as carboxylic acids, alcohols, aldehydes and ketones, esters, etc. The other group of typical metabolites are various triterpenoids [13,14,15].
The composition of cuticular waxes has been analyzed in detail in several studies using different analytical methods. For the wax extraction, several authors have used cuticle membranes enzymatically isolated by pectinase and cellulase [14,16]. On the other hand, simple methods based only on rinsing the apple surface with a solvent (usually chloroform, but also dichloromethane, petroleum ether, hexane, etc.) have also been reported [11,17,18,19]. The most common method for analyzing the isolated compounds is GC-MS [11,14,16,18], but this method requires derivatization. Reversed phase liquid chromatography coupled with mass spectrometry (RPLC-MS) has been increasingly used as an alternative, using both atmospheric pressure chemical ionization (APCI) and electrospray ionization (ESI) [20,21,22,23,24,25]. The advantage of LC-ESI-MS is the possibility of the simultaneous detection of different lipid classes, including difficult-to-ionize neutral wax esters when additives are used. The use of high-resolution (HR) MS allows for identification based on molecular or adduct ion exact masse, and the application of collision energy (tandem mass spectrometry, MS/MS) is suitable to obtain information on molecular structure [20].
The aim of the present study was to investigate the potential of UHPL-HRMS-based metabolomic fingerprinting of cuticular layer components, followed by advanced statistics, for the classification of apple cultivars “Gala”, “Golden Delicious”, “Idared” and “Jonagold”, grown either in the Czech Republic or in Poland. To our knowledge, no other published study has used such an approach.
2. Materials and Methods
2.1. Samples
Apple (Malus domestica) samples analyzed in this study were collected by the Research and Breeding Institute of Pomology Holovousy Ltd. (Holovousy, Czech Republic). A total of 274 authentic apple samples of known geographical origin and cultivar characterization were provided over a period of 3 harvest years (2020–2022). The cultivars available in this study were “Gala”, “Golden Delicious”, “Idared” and “Jonagold”, originating either from the Czech Republic (CZE) or Poland (POL). A total of 35 samples of the “Gala” cultivar (16 CZE, 19 POL), 37 samples of the “Golden Delicious” cultivar (19 CZE, 18 POL), 32 samples of the “Idared” cultivar (17 CZE, 15 POL) and 33 samples of the “Jonagold” cultivar (13 CZE, 19 POL) were available. After delivery to the laboratory, the samples were stored at 4 °C for a maximum of 3 days before further processing.
2.2. Chemicals
Analytical-grade methanol (MeOH), dichloromethane (DCM), ethyl acetate (EtAC), methyl-tert-butyl-ether (MTBE) and isopropanol (iPrOH) were purchased from Merck (Darmstadt, Germany). Deionized water (dH2O) was obtained from a Milli-Q Integral system (Millipore supplied by Merck (Darmstadt, Germany)). The mobile phase modifiers (ammonium formate, formic acid) were purchased from Sigma-Aldrich (Darmstadt, Germany).
2.3. Methods
2.3.1. Sample Preparation
To isolate the metabolites present in the apple cuticular layer, the whole apple was carefully placed into a 1000 mL glass beaker containing 400 mL of solvent mixture. After covering the beaker with parafilm, its content was gently shaken for 10 min. In the experiments undertaken to evaluate the extraction efficiency, the following solvents/mixtures thereof were tested: hexane–EtAc (1:1, v:v), DCM, DCM–MeOH (2:1, v:v), MTBE–MeOH (10:3, v:v) and DCM–MeOH (1:1, v:v). The last one was selected as optimal and then used throughout the study.
To obtain maximum yield of apple surface components, this procedure was repeated twice, always washing three representative apples. The combined extracts of each sample were evaporated stepwise to dryness and the residue was stored at −80 °C. Prior to analysis, the residue was reconstituted in a calculated amount of DCM:MeOH (1:1, v:v) to obtain a standardized concentration of the extracted material of 33.33 mg/mL. The solution was filtered using syringe filters (pore size 0.22 μm) and a 1.5 mL aliquot of each extract was then transferred to glass vials for LC-MS analysis. The quality control (QC) sample was prepared as a pool of aliquots of all prepared standardized apple extracts.
2.3.2. UHPLC-HRMS/MS Non-Target Screening
For the metabolomic analysis of apple extracts, the UHPLC-HRMS/MS technique was employed. A high-performance liquid chromatograph Dionex UltiMate 3000 RS (Thermo Fisher Scientific, Waltham, MA, USA) coupled with quadrupole-time-of-flight TripleTOFTM 6600 mass spectrometer (Sciex, Concord, ON, Canada) was used for this purpose. For sample components separation, Acquity UPLC BEH C18 column (2.1 × 100 mm, 1.7 μm) (Waters, Milford, MA, USA) was employed. The mobile phase consisted of (A) 5 mM ammonium formate in a mixture of dH2O:MeOH (95:5, v:v) with 0.1% formic acid and (B) 5 mM ammonium formate in a mixture of iPrOH:MeOH:dH2O (65:30:5, v:v) with 0.1% formic acid. The following gradient was used for both positive and negative ionization modes: 0 min (70% A), 2 min (50% A), 7 min (20% A), 13 min (0% A), 20 min (0% A), 20.1 min (70% A) and 22 min (70% A), with a constant flow rate of 0,4 mL/min. The column temperature was kept at 60 °C; the temperature of the autosampler at 5 °C and the sample injection volume was 1 μL.
The mass spectrometer was operated in both positive (ESI+) and negative (ESI-) mode with the following ion source settings: nebulizing gas pressure, 55 psi; drying gas pressure, 55 psi; capillary voltage, +4500 V (for ESI+)/−4000 V (for ESI-); ion temperature, 500 °C. Both MS and MS/MS data were acquired using full-scan and information-dependent acquisition (IDA) methods. The mass range in MS mode was set to 100–1200 m/z, and that in MS/MS mode to 50–1200 m/z. The collision energy was 35 V with the spread of ±15 V. Mass spectrometer calibration was performed regularly after every 10 samples based on APCI calibration solution (Sciex, Concord, Canada).
Samples were injected in a randomized order; QC samples were injected during the entire analytical run (after 10 previous sample injections). Blank samples (extraction solvent mixture) were injected at the beginning of the sequence to capture background features.
2.3.3. Data Processing
The UHPLC-HRMS/MS data obtained were processed using the open-source software MS-Dial (version 4.8) [26]. In the first step, the data were converted into the specific *.ibf format. The peak picking parameters were set as follows for both data acquired in ESI+ and ESI− mode: a minimum signal intensity threshold (peak height) of 10,000; a mass accuracy of 0.01 Da for MS data and 0.025 Da for MS/MS data. For data alignment, retention time tolerance of 0.05 min (ESI+) and 0.3 min (ESI−), along with m/z tolerance of 0.015 Da, were used. In ESI+ mode [M+H]+, [M+Na]+ and [M+NH4]+ adducts were considered, and in ESI- mode, [M-H]− and [M+HCOO]− adducts were considered.
The exported data matrices, consisting of all detected features characterized by m/z and retention time, were filtered according to the relative standard deviation (RSD) of the signal intensity (peak area) in the QC samples, with a maximum RSD threshold of 20% for both ESI+ and ESI−. Furthermore, all features with a signal-to-noise ratio (SNR) below 3 were filtered out to obtain the final data matrices.
2.3.4. Statistical Analysis
Prior to statistical analysis, the data were pre-processed to avoid possible misinterpretation of data variability. In this study, the total area sum normalization and logarithmic transformation were performed prior to any univariate statistics, followed by Pareto scaling in the case of multivariate model building.
Within the chemometric processing, the aim was to create models for the classification of apple samples using both univariate and multivariate statistical tools. Principal component analysis (PCA) was used to overview the data. Diagnostics features were selected using t-test/analysis of variance (ANOVA), fold change and receiver operating characteristics (ROC) methods. The combination of used methods, where each of them evaluates the feature significance based on a different algorithm, enables the selection of relevant markers. Based on the selected feature subset, both partial least square discriminant analysis (PLS-DA) and orthogonal partial least square discriminant analysis (OPLS-DA) were applied. The developed classification models were validated using 7-fold internal cross validation and characterized by the described variance (R2X and R2Y), the predicted variance (Q2Y), the root mean square error of estimation (RMSEE) and permutation tests for R2Y and Q2Y.
All statistical analyses were performed using SIMCA® (Sartorius, Göttingen, Germany), Metaboanalyst (metaboanalyst.ca, accessed on 16 January 2024) and using custom built R scripts.
2.3.5. Marker Identification
All significant features (significance is described by the results of the univariate statistical analysis or variable importance on the projection (VIP) score from PLS-DA/OPLS-DA) used for classification models building were subject of structure identification. For these features, *.mat files (containing both MS and MS/MS spectral information) were exported from MS-Dial and imported into the open-source software SIRIUS 4 [27,28], which also integrates CSI:FingerID [29] and CANOPUS [30,31]. Together, these three tools suggest possible molecular formulae, potential structures and compound classes for a given feature, which is compared with online spectral databases (BioCyc, HMDB, COCONUT) [32,33,34]. Tentatively identified markers were characterized by elemental formula, mass error and compound name. In addition, a confidence level of markers identification was classified according to the approaches used in previous studies [35,36,37] for identification of compounds based on LC-MS metabolomic data. The confidence levels range from Level 4 (unknown reproducible signal defined by m/z, retention time and MS spectrum), Level 3 (known compound class with many isomers possibilities), Level 2 (annotated compound based on matched MS/MS spectra and library) and Level 1 (identified compound confirmed with analytical standard) to Level 0 (identified compound including full stereochemistry).
3. Results and Discussion
This study, aimed at the authentication of apple cultivars and their geographical origin, was based on the assumption that characteristic metabolites, authenticity markers, could be identified via the statistical processing of the HPLC-HRMS/MS metabolic fingerprints of cuticular layer extracts. As described in the introduction, several studies have performed well in cultivar or geographical origin classification (Table 1); nevertheless, none of them focused on the analysis of the cuticular layer extracts, which, according to our working hypothesis, have a high application potential for authentication. In the paragraphs below, the steps taken to test this working hypothesis are described.
3.1. Selection of Extraction Solvent/Mixture
The first step was to find an extraction solvent that would enable the reproducible extraction (not necessarily quantitative) of the widest possible range of substances from the apple cuticle. The tested solvents/solvent mixtures, differing in their selectivity, involved hexane–EtAc (1:1, v:v), MTBE–MeOH (10:3, v:v), DCM, DCM–MeOH (2:1, v:v) and DCM–MeOH (1:1, v:v). The comparison of the total ion chromatograms of the tested solvents is shown in Supplementary Materials (Figure S1). The suitability of the extraction solvent was assessed by the distribution of chromatographic peaks in terms of their retention times and the total number of features detected via the reversed-phase UHPLC-HRMS method, which is commonly used by the authors of metabolomic studies. Since wax esters, which are known to occur in large amounts in the apple cuticular layer [13], ionize poorly in ESI- mode, only ESI+ was used for the experiments. A similar approach was used in other studies analyzing a similar matrix [20,22]. Based on the above criteria, the best solvent mixture was DCM–MeOH (1:1, v:v). A total of 11,581 features corresponding to compounds with a wide range of polarities were detected. Only a slightly lower number of features were detected in DCM and DCM–MeOH (2:1, v:v) extracts, while in the chromatograms obtained with the other solvent mixtures, hexane–EtAc (1:1, v:v) and MTBE–MeOH (10:3, v:v), the more polar metabolites, eluted at lower retention times, were not sufficiently represented.
3.2. UHPLC-HRMS/MS Analysis
A number of analytical strategies have been applied for the investigation of apples’ authenticity, either geographically or by variety. As shown in Table 1, in most cases, non-target screening performed by various instrumental techniques, such as NIR, GC-MS or IR-MS, was employed for analyses of various apple parts/forms. In contrast to these approaches, the UHPLC-HRMS/MS technique used in this study allows for not only the acquisition of the metabolomic fingerprints of the respective sample but also the identification of diagnostic markers (without pre-analytical derivatization) that can be used for authentication based on a simpler target screening applicable under routine conditions. The complexity of the apple cuticular layer extract is documented by the total ion chromatograms (TIC) in Figure 1A,B. As expected, wide range of lipophilic compounds was isolated; their main groups are indicated in the figures.
3.3. Chemometric Analysis
Processing of all raw UHPLC-HRMS/MS data with MS-Dial software resulted in the detection of 96,072 features in positive ionization mode and 21,040 features in negative ionization mode. Both of the aligned data matrices obtained were then filtered according to the criteria mentioned in Materials and Methods (Section 2.3.3) based on RSD in the QC sample (20%) and the SNR in the blank sample (≥3). As a result of feature filtering, the final data matrices contained 16,044 features and 2132 features, respectively, in ESI+ and ESI- mode, respectively. These data were used for further processing.
3.3.1. Data Overview
Data visualization via PCA revealed a clustering of the samples based on cultivar (Figure 2A). However, the separation of apples from the Czech Republic and Poland by geographical origin was not pronounced (Figure 2B). The data shown in these figures were obtained for the ESI+ mode; similar trends are documented in the PCA Score plots for the ESI− mode (see Supplementary Materials Figure S2).
In order to obtain reliable models for the classification of both apple cultivars and geographical origin, two alternative methods of supervised learning were combined. PLS-DA was used for cultivar classification, while a systematic strategy was used for geographical origin classification, as the influence of growing location on the metabolome of the cuticular layer is less pronounced. Individual binary OPLS-DA models were created for each apple cultivar. The use of this specific chemometric method allowed us to reduce the effect of collinearity and model overfitting [38].
3.3.2. Apple Cultivars Classification
Prior to building the PLS-DA classification models, an ANOVA false discovery rate (FDR) p-value threshold of 0.05 was used to filter out features that were unimportant for cultivar differentiation. This process resulted in a subset of 14,551 and 1678 significant features for the ESI+ and ESI- modes, respectively, which were used for building PLS-DA models (all four groups of cultivars included). For the model based on ESI+ data, the validation parameters were R2Y = 0.790 and Q2Y = 0.756; for the ESI- data, similar qualitative results were achieved, with R2Y = 0.770 and Q2Y = 0.756. Permutation tests (n = 100) were performed to validate the developed models; for both R2Y and Q2Y, the p-value was below 0.01, indicating valid models [39]. In general, a VIP score > 1 is considered as a threshold value for significant features [40]; nevertheless, in the study here presented, this threshold was increased up to 1.5, which resulted in exclusion of less significant metabolites. This way, only the 197 most significant markers (ESI+ and ESI− combined) were selected for further processing. Table 2 provides an overview of 13 metabolites whose molecular structure could be identified using the criteria specified in Materials and Methods (Section 2.3.5).
Flavonoids, wax esters and other lipids together with triterpenoids were the three classes of compounds identified as markers. The first group included isorhamnetin (methylated metabolite of quercetin), isorhamnetin rhamnoside and luteolin malonyl glucoside, secondary metabolites reported to be present at low levels in apple skin [41,42]; all were present at higher levels in “Idared” and “Jonagold” apple cultivars. The second group of significant markers (according to the VIP score) was identified as wax esters, a lipid subclass defined by the LIPID MAPS structure database (LMSD) [43]. Although many studies have indicated the presence of wax esters in the apple cuticular layer [12,44,45,46], the (tentative) identification of individual representatives relied exclusively on gas chromatography coupled to a flame ionization detector or mass-spectrometry with a simple mass analyzer, where the sample preparation typically involved the hydrolysis of ester bonds and a derivatization step aimed at increasing the volatility of the released fatty acids and alcohols under these conditions; some of the information about the wax structure is lost. Contrary to that approach, in another earlier comprehensive study [47], GC-MS analysis of whole molecules was performed; a database involving electron ionization mass spectra of 154 wax ester standards (various straight-chain and methyl-branched saturated and unsaturated species) was created. As regards LC-MS, several authors have investigated the mass spectra of wax esters obtained via this technique in more detail. In two older studies [48,49], atmospheric pressure ionization (APCI) and two types of mass analyzers (ion trap and Orbitrap) were used for this purpose. The dominant signals in the mass spectra were protonated molecular ions [M+H]+. Also in our study, where ESI was used instead of APCI, the most intensive ions in the wax ester spectra were [M+H]+. On the other hand, some papers report the use of ESI (and ammonium formate was a component of a mobile phase), but the most intensive adducts of wax esters were [M+NH4]+ [50,51]. In our study employing the Q-ToF mass analyzer, the fragmentation spectra of the [M+H]+ ion were further investigated and compared with those reported by Chen et al. [20], who focused on a systematic investigation of collision-induced dissociation (CID) patterns of different wax ester standards using HPLC-ESI-Q-ToF-MS. Based on this information, the fragmentation spectra of some wax esters present in our samples could be interpreted; an example of one of identified markers as wax ester (28:1/18:1) is shown in Figure 3.
The molecular formula C46H88O2 was calculated for the precursor ion m/z 673.6829 based on its exact mass and isotopic envelope. In the MS/MS spectrum, the base ion m/z 283.2659 was identified as an octadecenoic acid fragment [RCOOH2]+, which was formed by breaking the ester bond in the wax molecule. The fragments ions m/z 256.2564 and m/z 247.2445 were residues of oleic acid, namely, [RCO]+ and [RCO-H2O]+, respectively, which is in line with the spectra interpretation introduced by the aforementioned study [20]. The calculated difference between the precursor ion [M+H]+ (m/z 673.6829) and the fatty acyl fragment (m/z 283.2659) 390.4170 corresponded to the elemental composition C28H54, i.e., the loss of the octacosenol alkyl moiety. The fragmentation spectra for other wax esters with identified fatty acid and fatty alcohol moieties are shown in File S1. The description of the wax esters’ structure (number of carbon atoms in the alcohol; number of double bonds/number of carbons in the fatty acid; number of double bonds) corresponds to the nomenclature used in LMSD for this lipid subclass [52]. Another identified marker belonging to lipids was the fatty acid ester of hydroxy fatty acid (FAHFA, 18:1/16:1). These FAHFAs (previously known as indicators of inflammation in human samples) have only recently been quantified in some foods, including apple skin [53]. The other identified markers representing the class of acylglycerols were DAG (46:7), DAG (28:2) and TAG (60:2).
The last group of identified markers of apple cultivar were triterpenic acids, a group of compounds that was the only one analyzed in the apple cuticle by several authors using the UHPLC-HRMS/MS technique [22,23,25]. Caffeoylbetulinic acid, one of the tentatively identified triterpenic acids, was, together with caffeoyloleanolic acid, recognized earlier as a typical metabolite occurring in the skins of russeted apples [25]. It should be noted that susceptibility to russeting is rather cultivar-specific; some apples, such as Idared, develop this defect only very rarely.
It must be emphasized that despite a relatively small number of markers used for cultivars classification, the performance characteristics of the PSL-DA model were acceptable with R2Y = 0.691, Q2Y = 0.666 and recognition ability of 91% [39]. The boxplots illustrating markers can be found in Supplementary Materials (Figure S3).
3.3.3. Classification of Apple Geographical Origin
The comparison of the PCA cosre plots (Figure 2 and Figure S1) documents a bigger impact of the apple cultivar on the fingerprint of the cuticular-layer metabolome than that of the geographical origin. For this reason, the generated data were investigated more in depth and separately for each variety. The input file was a filtered data matrix with 16,044 and 2132 features, obtained in ESI+ and ESI− mode, respectively. Both t-test and ROC were then applied to the data and the fold change was calculated. An OPLS-DA model was created based on the subset of features that met the t-test FDR p-value < 0.05 and where the area under curve (AUC) value was higher than 0.75 [54]. Apart from the model validation performed, which was performed using seven-round cross-validation [55], the validity of the model was also proved over several harvest seasons. The predictive ability of the model was calculated by inserting samples from the 2020 and 2022 harvest seasons into the model created from samples harvested in 2021. The performance characteristics of all created OPLS-DA models for all cultivars are summarized in Table 3.
As shown in Table 3, all OPLS-DA models, except for the one classifying the geographical origin of “Jonagold” apples (ESI− data), performed satisfactorily [39]. The highest number of features differentiating between apples from Poland and the Czech Republic was found for the “Golden Delicious” cultivar. On the other hand, the worst performance of the developed classification models was found for the “Jonagold” cultivar, which is known to be susceptible to mutations (according to the experts from the Czech Research and Breeding Institute of Pomology Holovousy, there are 23 known mutations in “Jonagold” compared to 4 known mutations in “Golden Delicious”, 3 in “Idared” and 8 in “Gala”). Rather high variability of apple metabolite patterns in the case of “Jonagold” cultivar (even within the respective country) was obviously associated with its high susceptibility to mutations. For this reason, it was difficult to identify reliable diagnostic markers. When searching for experiences with the authentication of “Jonagold” in other studies, we found that Chinese authors [6] included this variety in their sample set, together with “Starkrimson”, “Qinguan”, “Gala”, “Golden Delicious” and “Fuji”. However, a different approach, the fingerprinting of volatiles in apple juices via the SPME-HS-GC-MS technique, was used by authors. Interestingly, contrary to other cultivars, the authentication of which was 100% successful, in the case of Jonagold, it was only 89%.
All features selected for geographical origin classification (distinguishing between Czech and Polish apple samples), i.e., those with ROC AUC value > 0.75 and t-test FDR p-value < 0.05, were subjected to the identification process. A total of, 40, 72 and 6 markers were identified for “Gala” (Table S1), “Golden Delicious” (Table S2) and “Idared” 6 (Table S3) cultivars, respectively. With regard to the facts mentioned above, no marker was identified for “Jonagold” cultivar.
The most frequent markers for geographical origin were wax esters, 23 of them for “Gala” and 28 for “Golden Delicious”. In addition to fatty acids, hydroxy fatty acids were also bound in wax ester molecules. The example of the fragmentation spectrum of such a marker (m/z 631.5975; retention time: 13.02 min), which was identified as a hydroxy wax ester (24:0/18:3-O), is shown in Figure 4. The obtained spectra were similar to those of the wax esters discussed in the previous section, but the calculated molecular formula (in this case, C42H78O3) contained three oxygen atoms instead of two. In the fragmentation spectra of these compounds, alike to that shown in Figure 5, there was a visible neutral loss of H2O molecules from both precursor and fragment ions corresponding to the hydroxy fatty acyl [56]. A total of 27 hydroxy wax esters were identified as markers of geographical origin in the analyzed apples. For those, in which fatty acid composition was tentatively identified, the corresponding fragmentation spectra are summarized in File S1.
In general, the most common fatty acids bound in waxes were hexadecanoic acid, hexadecenoic acid (16:0, 16:1) and unsaturated fatty acids containing 18 carbons (18:1, 18:2, 18:3, 18:4) bound to C22-C28 aliphatic alcohol. The heatmaps in Figure 5 show that in both “Gala” and “Golden Delicious” cultivars, the wax esters, as well as hydroxy wax esters, were upregulated in the apples harvested in the Czech Republic (exceptions were shorter chain wax esters, i.e., 30:3, 30:4, 30:5). In “Gala”, the differences were more pronounced than in “Golden Delicious”; the fold change of wax esters ranged from 1.21 to 2.37 and from 1.11 to 3,17, respectively. The higher intensities of wax esters signals (i.e., their higher amounts) observed in the samples from the Czech Republic could be due to the higher average altitude (and thus the different climatic conditions) of the orchards where the apples were collected (450 m in the Czech Republic versus 173 m in Poland), which has been previously reported as a parameter associated with thicker cuticles [57].
The next group of identified geographical origin markers were triterpenic acids, specifically ursane-type triterpenic acids. As already mentioned, these compounds, together with their derivatives (oxo, dihydroxy, oxohydroxy), have been identified in the apple cuticle layer by several authors [22,23].
In the case of our data, eight derivatives of ursolic acid were identified among the markers for the cultivar “Gala”. Interestingly, all ursenoic acids showed increased intensity in Czech apple samples. The higher intensity was statistically significant based on t-test p-value < 0.05 (except for oxohydroxy ursenoic acid), with fold changes ranging from 1.3 to 4. In “Golden Delicious”, in addition to the five ursolic acid derivatives, their precursors lupeone, uvaol, amyrin and hydroxybetulin were also found. Among markers for cultivar “Golden Delicious”, there was no clear trend. As an example, the boxplots of all ursenoic acids selected as markers are shown in the Supplementary Materials (Figures S4 and S5).
As relatively polar metabolites such as sorbitol, fucose, heptulose and mannose (trace amounts were contained in analyzed extracts), alike some flavonoids, represented by phloretin and chlorogenic acid common apple flavonoids [58,59], were also on the list of geographic origin markers.
4. Conclusions
In this study, the UHPLC-HRMS/MS-based metabolomic fingerprinting of cuticular layer extracts followed by advanced data processing was proved to be an applicable strategy for the authentication of apples cultivars and their geographic origin. The results of this study, within which 274 apple samples of four cultivars harvested either in the Czech Republic or Poland in three seasons were analyzed, can be summarized as follows:
PCA showed a more pronounced cultivar impact on the metabolites occurring in the apple cuticle compared to that of geographical origin.
The created PLS-DA models enabled reliable apple cultivar classification; 13 markers encompassing mainly waxes and triterpenoids were identified,
The created OPLS-DA models enabled the safe classification of geographical origins of “Gala”, “Golden Delicious” and “Idared” cultivars; however, for “Jonagold”, it was unsuccessful.
Wax esters, including those with bound hydroxy fatty acids (reported for the first time in apple cuticular wax), represented a significant group of identified markers, the amount of which in “Golden Delicious” and “Gala” cultivars was higher (upregulated) in samples from the Czech Republic compared those from Poland.
Overall, our findings underscore the potential of apple cuticular layer analysis to be used as a robust tool for apple authentication, offering insights into both cultivar and geographical origin distinctions. Moreover, revealing novel compounds that enhance our understanding of apple wax composition is facilitated by this approach.
Conceptualization, K.B., J.S., V.K. (Vladimir Kocourek) and J.H.; methodology, K.B., L.U. and P.V.; software, K.B.; formal analysis, K.B.; investigation, K.B., L.U. and P.V.; resources, J.S.; data curation, K.B.; writing—original draft preparation, K.B. and J.H.; writing—review and editing, J.H., V.K. (Vladimir Kocourek) and V.K. (Vit Kosek); supervision, J.H., V.K. (Vit Kosek), V.K. (Vladimir Kocourek) and J.H.; project administration, J.S., V.K. (Vladimir Kocourek) and J.H. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
The original contributions presented in the study are included in the article and
The authors acknowledge Michaela Rektorisova and Klara Navratilova for their support of this research.
The authors declare no conflicts of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Total ion chromatogram of apple cuticle layer extract, QC sample: (A) ESI+ mode and (B) ESI− mode.
Enlarge this image.Figure 2. PCA Score plots of the complete dataset of ESI+ features. Samples are colored according to (A) their cultivar and (B) their geographical origin.
Enlarge this image.Figure 4. Fragmentation spectra of the marker identified as hydroxy wax ester (24:0/18:3-O).
Enlarge this image.Figure 5. Heatmaps showing fold changes (red and blue colours represent positive and negative fold changes, respectively, darker shadows indicate higher absolute values) of wax ester (WE) and hydroxy wax ester (OHWE) intensities corresponding to the median value of the respective compound in apple samples (on the horizontal axis) harvested either in the Czech Republic (blue) or Poland (yellow).
The overview of studies dealing with the authentication of the origin and/or cultivar of apples.
| Analytical Method | Description of Apple Samples | Classification Factor | Number of Samples | Number of Classes to Be Distinguished within the Sample Set | Classification Method | Performance of Classification | Reference |
|---|---|---|---|---|---|---|---|
| NIR | Surface of whole apple fruits | Cultivar | 300 | 3 (Fuji, Red Star, Gala) | NN, SVM, ELM | Calibration set 98% (ELM) Prediction set 97% (ELM) | [ |
| Geographical origin | 2 (grown in different Chinese provinces) | ||||||
| Fluorescent spectroscopy | Apple juice (squeezed with a juice extractor) | Cultivar | 89 | 2 (grown in different Chinese provinces) | PLS | Calibration set 100% | [ |
| SPME-GC-MS | Apple juice (squeezed with a juicer) | Cultivar | 50 | 6 (Starkrimson, Qinguan, Gala, Jonagold, Golden Delicioius, Fuji) | LDA, SLDA | Predicition set 100% (SLDA) | [ |
| Geographical origin | 5 (grown in different counties within Chinese province) | Predicition set 90% (SLDA) | |||||
| SPME-GC-MS | Apple juice (squeezed with hand press) | Cultivar | 4 (3 kg of apples per sample) | 4 (Rijo, Verde, Ribeiro, Azedo) | PLS-DA, HCA | Vague description of model performance | [ |
| Geographical origin | 2 (different civil parishes of Madeira) | ||||||
| IR-MS + conventional methods | Pulp, juice | Cultivar | 19 | 6 (Topaz, Idared, Golden Delicious, Goldrush, Gala, Gloster) | LDA | Insufficient description of models | [ |
| Geographical origin | 4 (different regions of Slovenia) | ||||||
| Agricultural practice | 2 (way of farming organic, conventional) | ||||||
| IR-MS | Whole apples, peel, pulp, seed | Cultivar | 128 | 4 (Cripps Pink, Gala, Golden Delicious, Granny Smith) | LDA | 71% correctly classified samples | [ |
| Geographical origin | 4 (grown in different districts of northerm Italy) | 99% (LOOCV) | |||||
| IR-MS | Peel, petiole, pulp, seed | Geographical origin | 48 | 2 (grown in different districts of northern Italy) | LDA | Limited information on classification models performance | [ |
| IR-MS, ICP-MS | Apple juice (concentrated to sugar content 65.0°Brix) | Geographical origin | 135 | 6 (grown in different Chinese provinces) | LDA, PLS-DA | Only description of sample clustering in PLS-DA model without information about model validation | [ |
| Electronic nose, electronic tongue | Apple juice (centrifugal juicer) | Cultivar | 126 | 10 (Fuji, Jonagold, Corolla, Gala, Red Delicous, Red Chief Delicious, Cattle Apple, Ralls Janet, Ourin, Tail, Golden Delicous) | LDA, PLS-DA, SVM | 100% (prediction ability of PLS-DA) 100% (accuracy testing rate of SVM) | [ |
| Geographical origin | 7 (grown in different Chinese provinces) |
ELM—extreme learning machine; HCA—hierarchical cluster analysis; LDA—linear discriminant analysis; LOOCV—leave-one-out cross validation; NN—neural networks; PLS—partial least square; PLS-DA—partial least square discriminant analysis; SLDA—stepwise linear discriminant analysis; SVM—support vector machine.
Identified significant metabolites used for apple cultivar classification (compounds sorted in descending order according to the PLS-DA VIP score).
| Marker Ion (m/z) | Retention Time [min] | Adduct Type | Elemental Formula | Mass Error [ppm] | Tentative Identification | PLSDA VIP Score | Confidence Level |
|---|---|---|---|---|---|---|---|
| 701.7138 | 14.03 | [M+H]+ | C48H92O2 | −5.4 | Wax ester (30:1/18:1) | 3.1 | 2 |
| 317.064 | 2.06 | [M+H]+ | C16H12O7 | −6.7 | Isorhamnetine | 2.9 | 3 |
| 673.6829 | 13.69 | [M+H]+ | C46H88O2 | −5 | Wax ester (28:1/18:1) | 2.9 | 2 |
| 461.1111 | 2.14 | [M-H]− | C22H22O11 | 5.9 | Isorhamnetin rhamnoside | 2.8 | 3 |
| 671.6652 | 13.43 | [M+H]+ | C46H86O2 | -8 | Wax esters (46:3) | 2.8 | 3 |
| 699.691 | 14.64 | [M+Na]+ | C46H92O2 | −12.1 | Wax esters (46:0) | 2.7 | 3 |
| 979.8971 | 14.81 | [M+Na]+ | C63H120O5 | −6.4 | TAG (60:2) | 2.4 | 3 |
| 509.4234 | 12.24 | [M+H]+ | C31H56O5 | 5.6 | DAG (28:2) | 2.1 | 3 |
| 533.0917 | 1.33 | [M-H]− | C24H22O14 | 2.6 | Luteolin-O-malonyl glucoside | 2.1 | 3 |
| 663.3906 | 5.99 | [M+HCOO]− | C39H54O6 | 1.3 | Caffeoylbetulinic acid | 2 | 3 |
| 535.4747 | 12.87 | [M-H]− | C34H64O4 | 3.8 | FAHFA (18:1/16:0) | 1.9 | 2 |
| 549.3436 | 3.7 | [M+HCOO]− | C30H48O6 | 1.6 | Triterpenic acid | 1.6 | 3 |
| 749.6105 | 13.21 | [M-H]− | C49H82O5 | 2.8 | DAG (46:7) | 1.6 | 3 |
DAG—diacylglycerol; FAHFA—fatty acid ester of hydroxy fatty acid; PLSDA VIP—variable importance in projection of PLSDA model; TAG—triacylglycerol
Parameters of the OPLS-DA models for the classification of apple samples according to the geographical origin (Poland vs. Czech Republic).
| OPLS-DA Model Parameters | ESI+ | ESI− | ||||||
|---|---|---|---|---|---|---|---|---|
| Gala | Golden Delicious | Idared | Jonagold | Gala | Golden Delicious | Idared | Jonagold | |
| number of features | 506 | 1048 | 156 | 11 | 13 | 44 | 24 | 9 |
| R2X | 0.783 | 0.596 | 0.570 | 0.946 | 0.667 | 0.567 | 0.850 | 0.921 |
| R2Y | 0.735 | 0.635 | 0.886 | 0.561 | 0.639 | 0.738 | 0.646 | 0.480 |
| Q2Y | 0.624 | 0.554 | 0.809 | 0.501 | 0.543 | 0.686 | 0.574 | 0.436 |
| RMSEE | 0.265 | 0.309 | 0.175 | 0.335 | 0.307 | 0.261 | 0.308 | 0.362 |
| p-value of permutation for R2Y | <0.01 | <0.01 | <0.01 | <0.01 | <0.01 | <0.01 | <0.01 | <0.01 |
| p-value of permutation for Q2Y | <0.01 | <0.01 | <0.01 | <0.01 | <0.01 | <0.01 | <0.01 | <0.01 |
| validity of the model over time | 82% | 65% | 85% | 88% | 78% | 77% | 78% | 88% |
R2X—fraction of X variation described by the model; R2Y—fraction of Y variation described by the model; Q2Y—fraction of Y variation predicted by model according to the cross validation; RMSEE—root mean square error of estimation.
Supplementary Materials
The following supporting information can be downloaded at
References
1. Shanbandeh, M. Apple Production Worldwide in 2022. Available online: https://www.statista.com/statistics/237605/production-of-apples-worldwide-by-region-2007/ (accessed on 29 February 2024).
2. Bat, K.; Vidrih, R.; Nečemer, M.; Mozetič Vodopivec, B.; Mulič, I.; Kump, P.; Ogrinc, N. Characterization of Slovenian apples with respect to their botanical and geographical origin and agricultural production practice. Food Technol. Biotechnol.; 2012; 50, pp. 107-116.
3. Wu, H.; Yue, T.; Yuan, Y. Authenticity Tracing of Apples According to Variety and Geographical Origin Based on Electronic Nose and Electronic Tongue. Food Anal. Methods; 2018; 11, pp. 522-532. [DOI: https://dx.doi.org/10.1007/s12161-017-1023-y]
4. Li, C.; Li, L.; Wu, Y.; Lu, M.; Yang, Y.; Li, L. Apple Variety Identification Using Near-Infrared Spectroscopy. J. Spectrosc.; 2018; 2018, e6935197. [DOI: https://dx.doi.org/10.1155/2018/6935197]
5. Bian, H.; Sheng, L.; Yao, H.; Ji, R.; Yu, Y.; Chen, R.; Wei, D.; Han, Y. Application of fluorescence spectroscopy in classifying apple juice according to the variety. Optik; 2021; 231, 166361. [DOI: https://dx.doi.org/10.1016/j.ijleo.2021.166361]
6. Guo, J.; Yue, T.; Yuan, Y. Feature selection and recognition from nonspecific volatile profiles for discrimination of apple juices according to variety and geographical origin. J. Food Sci.; 2012; 77, pp. C1090-C1096. [DOI: https://dx.doi.org/10.1111/j.1750-3841.2012.02914.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23009695]
7. Medina, S.; Perestrelo, R.; Santos, R.; Pereira, R.; Câmara, J.S. Differential volatile organic compounds signatures of apple juices from Madeira Island according to variety and geographical origin. Microchem. J.; 2019; 150, 104094. [DOI: https://dx.doi.org/10.1016/j.microc.2019.104094]
8. Mimmo, T.; Camin, F.; Bontempo, L.; Capici, C.; Tagliavini, M.; Cesco, S.; Scampicchio, M. Traceability of different apple varieties by multivariate analysis of isotope ratio mass spectrometry data. Rapid Commun. Mass Spectrom.; 2015; 29, pp. 1984-1990. [DOI: https://dx.doi.org/10.1002/rcm.7306] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26443397]
9. Brombin, V.; Mistri, E.; Bianchini, G. Multi stable isotope ratio analysis for the traceability of northern Italian apples. Food Chem. X; 2022; 16, 100514. [DOI: https://dx.doi.org/10.1016/j.fochx.2022.100514] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36519099]
10. Liu, X.; Zhao, Y.; Mu, J.; Zhang, J.; Zhang, A. Determination of geographical origin of concentrated apple juice through analysis of stable isotopic and mineral elemental fingerprints: Preliminary results. J. Sci. Food Agric.; 2021; 101, pp. 3795-3803. [DOI: https://dx.doi.org/10.1002/jsfa.11012]
11. Belding, R.D.; Blankenship, S.M.; Young, E.; Leidy, R.B. Composition and Variability of Epicuticular Waxes in Apple Cultivars. J. Am. Soc. Hort. Sci.; 1998; 123, pp. 348-356. [DOI: https://dx.doi.org/10.21273/JASHS.123.3.348]
12. Yang, Y.; Zhou, B.; Zhang, J.; Wang, C.; Liu, C.; Liu, Y.; Zhu, X.; Ren, X. Relationships between cuticular waxes and skin greasiness of apples during storage. Postharv. Biol. Technol.; 2017; 131, pp. 55-67. [DOI: https://dx.doi.org/10.1016/j.postharvbio.2017.05.006]
13. Wu, W.; Jiang, B.; Liu, R.; Han, Y.; Fang, X.; Mu, H.; Farag, M.; Simal-Gandara, J.; Prieto, M.A.; Chen, H. et al. Structures and functions of cuticular wax in postharvest fruit and its regulation: A comprehensive review with future perspectives. Engineering; 2023; 23, pp. 118-129. [DOI: https://dx.doi.org/10.1016/j.eng.2022.12.006]
14. Leide, J.; Xavier de Souza, A.; Papp, I.; Riederer, M. Specific characteristics of the apple fruit cuticle: Investigation of early and late season cultivars ‘Prima’ and ‘Florina’ (Malus domestica Borkh.). Sci. Hortic.; 2018; 229, pp. 137-147. [DOI: https://dx.doi.org/10.1016/j.scienta.2017.10.042]
15. Shepherd, T.; Wynne Griffiths, D. The effects of stress on plant cuticular waxes. New Phytol.; 2006; 171, pp. 469-499. [DOI: https://dx.doi.org/10.1111/j.1469-8137.2006.01826.x]
16. Fernandez-Moreno, J.-P.; Malitsky, S.; Lashbrooke, J.; Biswal, A.K.; Racovita, R.C.; Mellerowicz, E.J.; Jetter, R.; Orzaez, D.; Aharoni, A.; Granell, A. An efficient method for medium throughput screening of cuticular wax composition in different plant species. Metabolomics; 2016; 12, 73. [DOI: https://dx.doi.org/10.1007/s11306-016-0982-0]
17. Lara, I.; Belge, B.; Goulao, L.F. A Focus on the Biosynthesis and Composition of Cuticle in Fruits. J. Agric. Food Chem.; 2015; 63, pp. 4005-4019. [DOI: https://dx.doi.org/10.1021/acs.jafc.5b00013] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25850334]
18. Li, F.; Min, D.; Song, B.; Shao, S.; Zhang, X. Ethylene effects on apple fruit cuticular wax composition and content during cold storage. Postharv. Biol. Technol.; 2017; 134, pp. 98-105. [DOI: https://dx.doi.org/10.1016/j.postharvbio.2017.08.011]
19. Li, F.; Min, D.; Ren, C.; Dong, L.; Shu, P.; Cui, X.; Zhang, X. Ethylene altered fruit cuticular wax, the expression of cuticular wax synthesis-related genes and fruit quality during cold storage of apple (Malus domestica Borkh. C.v. Starkrimson) fruit. Postharv. Biol. Technol.; 2019; 149, pp. 58-65. [DOI: https://dx.doi.org/10.1016/j.postharvbio.2018.11.016]
20. Chen, J.; Green, K.B.; Nichols, K.K. Characterization of Wax Esters by Electrospray Ionization Tandem Mass Spectrometry: Double Bond Effect and Unusual Product Ions. Lipids; 2015; 50, pp. 821-836. [DOI: https://dx.doi.org/10.1007/s11745-015-4044-6]
21. Tada, A.; Jin, Z.-L.; Sugimoto, N.; Sato, K.; Yamazaki, T.; Tanamoto, K. Analysis of the constituents in jojoba wax used as a food additive by LC/MS/MS. Shokuhin Eiseigaku Zasshi. J. Food Hyg. Soc. Jpn; 2005; 46, pp. 198-204. [DOI: https://dx.doi.org/10.3358/shokueishi.46.198]
22. McGhie, T.K.; Hudault, S.; Lunken, R.C.M.; Christeller, J.T. Apple peels, from seven cultivars, have lipase-inhibitory activity and contain numerous ursenoic acids as identified by LC-ESI-QTOF-HRMS. J. Agric. Food Chem.; 2012; 60, pp. 482-491. [DOI: https://dx.doi.org/10.1021/jf203970j] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22148752]
23. Poirier, B.C.; Buchanan, D.A.; Rudell, D.R.; Mattheis, J.P. Differential Partitioning of Triterpenes and Triterpene Esters in Apple Peel. J. Agric. Food Chem.; 2018; 66, pp. 1800-1806. [DOI: https://dx.doi.org/10.1021/acs.jafc.7b04509] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29356521]
24. Andre, C.M.; Greenwood, J.M.; Walker, E.G.; Rassam, M.; Sullivan, M.; Evers, D.; Perry, N.B.; Laing, W.A. Anti-inflammatory procyanidins and triterpenes in 109 apple varieties. J. Agric. Food Chem.; 2012; 60, pp. 10546-10554. [DOI: https://dx.doi.org/10.1021/jf302809k] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23013475]
25. Andre, C.M.; Larsen, L.; Burgess, E.J.; Jensen, D.J.; Cooney, J.M.; Evers, D.; Zhang, J.; Perry, N.B.; Laing, W.A. Unusual immuno-modulatory triterpene-caffeates in the skins of russeted varieties of apples and pears. J. Agric. Food Chem.; 2013; 61, pp. 2773-2779. [DOI: https://dx.doi.org/10.1021/jf305190e] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23418665]
26. Tsugawa, H.; Kind, T.; Nakabayashi, R.; Yukihira, D.; Tanaka, W.; Cajka, T.; Saito, K.; Fiehn, O.; Arita, M. Hydrogen Rearrangement Rules: Computational MS/MS Fragmentation and Structure Elucidation Using MS-FINDER Software. Anal. Chem.; 2016; 88, pp. 7946-7958. [DOI: https://dx.doi.org/10.1021/acs.analchem.6b00770] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27419259]
27. Dührkop, K.; Fleischauer, M.; Ludwig, M.; Aksenov, A.A.; Melnik, A.V.; Meusel, M.; Dorrestein, P.C.; Rousu, J.; Böcker, S. SIRIUS 4: A rapid tool for turning tandem mass spectra into metabolite structure information. Nat. Methods; 2019; 16, 4. [DOI: https://dx.doi.org/10.1038/s41592-019-0344-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30886413]
28. Dührkop, K.; Nothias, L.-F.; Fleischauer, M.; Reher, R.; Ludwig, M.; Hoffmann, M.A.; Petras, D.; Gerwick, W.H.; Rousu, J.; Dorrestein, P.C. et al. Systematic classification of unknown metabolites using high-resolution fragmentation mass spectra. Nat. Biotechnol.; 2021; 39, 4. [DOI: https://dx.doi.org/10.1038/s41587-020-0740-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33230292]
29. Dührkop, K.; Shen, H.; Meusel, M.; Rousu, J.; Böcker, S. Searching molecular structure databases with tandem mass spectra using CSI:FingerID. Proc. Natl. Acad. Sci. USA; 2015; 112, pp. 12580-12585. [DOI: https://dx.doi.org/10.1073/pnas.1509788112] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26392543]
30. Kim, H.W.; Wang, M.; Leber, C.A.; Nothias, L.-F.; Reher, R.; Kang, K.B.; van der Hooft, J.J.J.; Dorrestein, P.C.; Gerwick, W.H.; Cottrell, G.W. NPClassifier: A Deep Neural Network-Based Structural Classification Tool for Natural Products. J. Nat. Prod.; 2021; 84, pp. 2795-2807. [DOI: https://dx.doi.org/10.1021/acs.jnatprod.1c00399]
31. Djoumbou Feunang, Y.; Eisner, R.; Knox, C.; Chepelev, L.; Hastings, J.; Owen, G.; Fahy, E.; Steinbeck, C.; Subramanian, S.; Bolton, E. et al. ClassyFire: Automated chemical classification with a comprehensive, computable taxonomy. J. Cheminform.; 2016; 8, 61. [DOI: https://dx.doi.org/10.1186/s13321-016-0174-y]
32. Karp, P.D.; Billington, R.; Caspi, R.; Fulcher, C.A.; Latendresse, M.; Kothari, A.; Keseler, I.M.; Krummenacker, M.; Midford, P.E.; Ong, Q. et al. The BioCyc collection of microbial genomes and metabolic pathways. Brief. Bioinform.; 2019; 20, pp. 1085-1093. [DOI: https://dx.doi.org/10.1093/bib/bbx085] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29447345]
33. Sorokina, M.; Merseburger, P.; Rajan, K.; Yirik, M.A.; Steinbeck, C. COCONUT online: Collection of open natural products database. J. Cheminform.; 2021; 13, 2. [DOI: https://dx.doi.org/10.1186/s13321-020-00478-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33423696]
34. Wishart, D.S.; Guo, A.C.; Oler, E.; Wang, F.; Anjum, A.; Peters, H.; Dizon, R.; Sayeeda, Z.; Tian, S.; Lee, B.L. et al. HMDB 5.0: The Human Metabolome Database for 2022. Nucleic Acids Res.; 2022; 50, pp. D622-D631. [DOI: https://dx.doi.org/10.1093/nar/gkab1062] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34986597]
35. Reisdorph, N.A.; Walmsley, S.; Reisdorph, R. A Perspective and Framework for Developing Sample Type Specific Databases for LC/MS-Based Clinical Metabolomics. Metabolites; 2019; 10, 8. [DOI: https://dx.doi.org/10.3390/metabo10010008] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31877765]
36. Blaženović, I.; Kind, T.; Ji, J.; Fiehn, O. Software Tools and Approaches for Compound Identification of LC-MS/MS Data in Metabolomics. Metabolites; 2018; 8, 31. [DOI: https://dx.doi.org/10.3390/metabo8020031]
37. Sumner, L.W.; Amberg, A.; Barrett, D.; Beale, M.H.; Beger, R.; Daykin, C.A.; Fan, T.W.-M.; Fiehn, O.; Goodacre, R.; Griffin, J.L. et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics; 2007; 3, pp. 211-221. [DOI: https://dx.doi.org/10.1007/s11306-007-0082-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24039616]
38. Bylesjö, M.; Rantalainen, M.; Cloarec, O.; Nicholson, J.K.; Holmes, E.; Trygg, J. OPLS discriminant analysis: Combining the strengths of PLS-DA and SIMCA classification. J. Chemom.; 2006; 20, pp. 341-351. [DOI: https://dx.doi.org/10.1002/cem.1006]
39. Triba, M.N.; Le Moyec, L.; Amathieu, R.; Goossens, C.; Bouchemal, N.; Nahon, P.; Rutledge, D.N.; Savarin, P. PLS/OPLS models in metabolomics: The impact of permutation of dataset rows on the K-fold cross-validation quality parameters. Mol. BioSyst.; 2015; 11, pp. 13-19. [DOI: https://dx.doi.org/10.1039/c4mb00414k]
40. Cocchi, M.; Biancolillo, A.; Marini, F. Chapter Ten—Chemometric Methods for Classification and Feature Selection. Comprehensive Analytical Chemistry; Jaumot, J.; Bedia, C.; Tauler, R. Data Analysis for Omic Sciences: Methods and Applications Elsevier: Amsterdam, The Netherlands, 2018; Volume 82, pp. 265-299. [DOI: https://dx.doi.org/10.1016/bs.coac.2018.08.006]
41. Pennington, J.A.T. Food Composition Databases for Bioactive Food Components. J. Food Comp. Anal.; 2002; 15, pp. 419-434. [DOI: https://dx.doi.org/10.1006/jfca.2002.1073]
42. Ramírez-Ambrosi, M.; López-Márquez, D.M.; Abad-García, B.; Dapena, E.; Berrueta, L.Á.; Gallo, B. Comparative study of phenolic profile of fruit and juice samples of a progeny of ‘Meana’ × ‘Florina’ from an Asturian cider apple breeding program. Eur. Food Res. Technol.; 2015; 241, pp. 769-784. [DOI: https://dx.doi.org/10.1007/s00217-015-2502-2]
43. Sud, M.; Fahy, E.; Cotter, D.; Brown, A.; Dennis, E.A.; Glass, C.K.; Merrill Jr, A.H.; Murphy, R.C.; Raetz, C.R.H.; Russell, D.W. et al. LMSD: LIPID MAPS structure database. Nucleic Acids Res.; 2007; 35, pp. D527-D532. [DOI: https://dx.doi.org/10.1093/nar/gkl838] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17098933]
44. Veraverbeke, E.A.; Lammertyn, J.; Saevels, S.; Nicolaï, B.M. Changes in chemical wax composition of three different apple (Malus domestica Borkh.) cultivars during storage. Postharv. Biol. Technol.; 2001; 23, pp. 197-208. [DOI: https://dx.doi.org/10.1016/S0925-5214(01)00128-4]
45. Klein, B.; Falk, R.B.; Thewes, F.R.; Anese, R.d.O.; dos Santos, I.D.; Ribeiro, S.R.; Donadel, J.Z.; Brackmann, A.; Barin, J.S.; Cichoski, A.J. et al. Dynamic controlled atmosphere: Effects on the chemical composition of cuticular wax of ‘Cripps Pink’ apples after long-term storage. Postharv. Biol. Technol.; 2020; 164, 111170. [DOI: https://dx.doi.org/10.1016/j.postharvbio.2020.111170]
46. Morice, I.M.; Shorland, F.B. Composition of the surface waxes of apple fruits and changes during storage. J. Sci. Food Agric.; 1973; 24, pp. 1331-1339. [DOI: https://dx.doi.org/10.1002/jsfa.2740241104] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/4763678]
47. Urbanova, K.; Vrkoslav, V.; Valterova, I.; Haková, M.; Cvacka, J. Structural characterization of wax esters by electron ionization mass spectrometry. J. Lipid Res.; 2011; 53, pp. 204-213. [DOI: https://dx.doi.org/10.1194/jlr.d020834] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22058425]
48. Vrkoslav, V.; Urbanová, K.; Cvacka, J. Analysis of wax ester molecular species by high performance liquid chromatography/atmospheric pressure chemical ionisation mass spectrometry. J. Chromatogr. A; 2010; 1217, pp. 4184-4194. [DOI: https://dx.doi.org/10.1016/j.chroma.2009.12.048] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20079497]
49. Vrkoslav, V.; Urbanova, K.; Haková, M.; Cvacka, J. Analysis of wax esters by silver-ion high-performance liquid chromatography–tandem mass spectrometry. J. Chromatogr. A; 2013; 1302, pp. 105-110. [DOI: https://dx.doi.org/10.1016/j.chroma.2013.06.031] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23831001]
50. Sekosan, G. LC-MS Identification of Wax Esters in Cloudy Canola Oil. LCGC Suppl.; 2014; 12, pp. 16-21.
51. Iven, T.; Herrfurth, C.; Hornung, E.; Heilmann, M.; Hofvander, P.; Stymne, S.; Feussner, I. Wax ester profiling of seed oil by nano-electrospray ionization tandem mass spectrometry. Plant Methods; 2013; 9, 24. [DOI: https://dx.doi.org/10.1186/1746-4811-9-24]
52. Liebisch, G.; Fahy, E.; Aoki, J.; Dennis, E.A.; Durand, T.; Ejsing, C.S.; Fedorova, M.; Feussner, I.; Griffiths, W.J.; Köfeler, H. et al. Update on LIPID MAPS classification, nomenclature, and shorthand notation for MS-derived lipid structures. J. Lipid Res.; 2020; 61, pp. 1539-1555. [DOI: https://dx.doi.org/10.1194/jlr.s120001025]
53. Liberati-Čizmek, A.-M.; Biluš, M.; Brkić, A.L.; Barić, I.C.; Bakula, M.; Hozić, A.; Cindrić, M. Analysis of Fatty Acid Esters of Hydroxyl Fatty Acid in Selected Plant Food. Plant Foods Hum. Nutr.; 2019; 74, pp. 235-240. [DOI: https://dx.doi.org/10.1007/s11130-019-00728-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30993530]
54. Mandrekar, J.N. Receiver Operating Characteristic Curve in Diagnostic Test Assessment. J. Thorac. Oncol.; 2010; 5, pp. 1315-1316. [DOI: https://dx.doi.org/10.1097/jto.0b013e3181ec173d] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20736804]
55. Westerhuis, J.A.; Hoefsloot, H.C.J.; Smit, S.; Vis, D.J.; Smilde, A.K.; van Velzen, E.J.J.; van Duijnhoven, J.P.M.; van Dorsten, F.A. Assessment of PLSDA cross validation. Metabolomics; 2008; 4, pp. 81-89. [DOI: https://dx.doi.org/10.1007/s11306-007-0099-6]
56. Spickett, C.M.; Pitt, A.R. Oxidative Lipidomics Coming of Age: Advances in Analysis of Oxidized Phospholipids in Physiology and Pathology. Antioxid. Redox Signal.; 2015; 22, pp. 1646-1666. [DOI: https://dx.doi.org/10.1089/ars.2014.6098] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25694038]
57. Tessmer, M.A.; Antoniolli, L.R.; Appezzato-da-Glória, B. Cuticle of “Gala” and “Galaxy” apples cultivars under different environmental conditions. Braz. Arch. Biol. Technol.; 2012; 55, pp. 709-714. [DOI: https://dx.doi.org/10.1590/S1516-89132012000500010]
58. Lommen, A.; Godejohann, M.; Venema, D.P.; Hollman, P.C.; Spraul, M. Application of directly coupled HPLC-NMR-MS to the identification and confirmation of quercetin glycosides and phloretin glycosides in apple peel. Anal. Chem.; 2000; 72, pp. 1793-1797. [DOI: https://dx.doi.org/10.1021/ac9912303]
59. Awad, M.A.; de Jager, A.; van Westing, L.M. Flavonoid and chlorogenic acid levels in apple fruit: Characterisation of variation. Sci. Hort.; 2000; 83, pp. 249-263. [DOI: https://dx.doi.org/10.1016/S0304-4238(99)00124-7]
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Sedlak, Jiri 2
; Leos Uttl 1
; Kosek, Vit 1
; Vackova, Petra 1 ; Kocourek, Vladimir 1
; Hajslova, Jana 1 1 Department of Food Analysis and Nutrition, University of Chemistry and Technology, Technicka 3, 16628 Prague 6, Czech Republic;
2 Reserach and Breeding Institute of Pomology Holovousy, Holovousy 129, 50801 Holovousy, Czech Republic;