Next Article in Journal
Effects of Extraction Conditions on Banana Peel Polyphenol Oxidase Activity and Insights into Inactivation Kinetics Using Thermal and Cold Plasma Treatment
Next Article in Special Issue
Combination of Stable Isotopes and Fatty Acid Composition for Geographical Origin Discrimination of One Argan Oil Vintage
Previous Article in Journal
Change Regularity of Taste and the Performance of Endogenous Proteases in Shrimp (Penaens vannamei) Head during Autolysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Preliminary Study on Tracing the Origin and Exploring the Relations between Growing Conditions and Isotopic and Elemental Fingerprints of Organic and Conventional Cavendish Bananas (Musa spp.)

1
Food Quality & Design Group, Wageningen University & Research, P.O. Box 17, 6700 AA Wageningen, The Netherlands
2
Wageningen Food Safety Research, Wageningen University & Research, P.O. Box 230, 6700 AE Wageningen, The Netherlands
*
Author to whom correspondence should be addressed.
Submission received: 6 April 2021 / Revised: 29 April 2021 / Accepted: 5 May 2021 / Published: 8 May 2021
(This article belongs to the Special Issue Food Origin Analysis with Isotope Fingerprints)

Abstract

:
The stable isotopic ratios and elemental compositions of 120 banana samples, Musa spp. (AAA Group, Cavendish Subgroup) cultivar Williams, collected from six countries (Colombia, Costa Rica, Dominica Republic, Ecuador, Panama, Peru), were determined by isotope ratio mass spectrometry and inductively coupled plasma mass spectrometry. Growing conditions like altitude, temperature, rainfall and production system (organic or conventional cultivation) were obtained from the sampling farms. Principal component analysis (PCA) revealed separation of the farms based on geographical origin and production system. The results showed a significant difference in the stable isotopic ratios (δ13C, δ15N, δ18O) and elemental compositions (Al, Ba, Cu, Fe, Mn, Mo, Ni, Rb) of the pulp and peel samples. Furthermore, δ15N was found to be a good marker for organically produced bananas. A correlation analysis was conducted to show the linkage of growing conditions and compositional attributes. The δ13C of pulp and peel were mainly negatively correlated with the rainfall, while δ18O was moderately positively (R values ~0.5) correlated with altitude and temperature. A moderate correlation was also found between temperature and elements such as Ba, Fe, Mn, Ni and Sr in the pulp and peel samples. The PCA results and correlation analysis suggested that the differences of banana compositions were combined effects of geographical factors and production systems. Ultimately, the findings contribute towards understanding the compositional differences of bananas due to different growing conditions and production systems linked to a defined origin; thereby offering a tool to support the traceability of commercial fruits.

Graphical Abstract

1. Introduction

Banana (Musa spp.) is one of the most representative fruits which are well-liked by consumers for their flavor and nutritional value. Generally, the cooking banana, also called plantain, is used as a staple food for millions of people in many countries, especially in developing economies, because it is an affordable source of energy [1]. Another kind of popular banana is called the dessert banana, which is used as an everyday fruit because of its soft texture, sweet taste and bioactive components [2]. The dessert banana is typically used as a food companion with products such as breakfast cereals, ice cream, and other desserts apart from just its raw consumption. In the current international market, the Cavendish variety is the most commercialized bananas and consumed world-wide [3]. Cavendish bananas are rich in carbohydrates, vitamins, minerals and dietary fiber, indicating that they have anti-oxidative and anti-ageing activities [4].
In 2019 alone, around 21 million tons of bananas, excluding plantains, were exported from Ecuador, Colombia, Costa Rica, etc. Latin America and the Caribbean are the largest exporting regions, and the European Union (EU), United States and China are the largest global importers [5]. As the quality of the banana is highly dependent on climate and cultivation system, the label of origins is considered an important parameter for quality and price [6]. The production of organic bananas requires stricter management, such as forbidding the use of pesticides, chemical fertilizers, and pesticides [7]. In the case of consumer preferences, the geographical origin and production system play a role in the purchase intent of consumers [8]. However, there are only a few papers about the authentication of the geographical origin and production system of bananas.
In recent years, several papers have been published concerning the geographical origin and organic production of agricultural products [9,10,11]. This is because food fraud and public health incidents occur frequently around the world. Food fraud, mainly known as economically motivated adulteration, is an ongoing problem [12]. Although banana is not an extremely vulnerable target, authenticity research of banana still has important value [13]. The banana is a staple food in many developing countries to ensure household food security. Meanwhile, the income from the banana industry is also important for approximately 400 million workers globally [5]. In recent years, the production of bananas is facing many challenges, such as infectious diseases Fusarium wilt, mislabeled organic bananas, sustainable production with less pesticides, etc. A reliable tracing system is helpful in banana chains to protect smallholder farmers, the retailer from unfair competition, and make sure the consumers get what they pay for [14]. Previous studies have reported that the combined application of stable isotopes and elemental compositions could effectively identify the geographical origin of a variety of agricultural products such as grape wine [15], pork [16], honey [17], tea [18], etc. Nevertheless, the measurement of stable isotopes could also reveal if the foods were produced in organic cultivations. The combination of δ13C, δ15N and δ18O could provide a sound solution for rice origin using chemometrics [19]. For bell peppers, the linking of δ18O in the sample and local water could indicate the geographical origin [20]. Previous research also confirmed that elemental analysis of Mn, Cr, Sr, Ag and Co could differentiate wines from different Romanian vineyards [21]. In terms of bananas, there are limited papers about the geographical origin of bananas. One relevant study explored the impact of farming type, variety and geographical origin on bacterial community of banana, and was aimed at discriminating organic bananas from conventional ones [22]. Most of the published research focused on the discrimination of geographical origin and production system, while only a few research studies aimed to explore the underlining reasons that foods have various composition differences according to different origins and production systems. There are already studies which reported that the lithological properties of origin could change the element contents of garlic [23], while latitudes had significant influences on the depleted δ2H and δ18O values of poultry [24]. Our previous study revealed that the isotopic ratios and elemental profiles of bananas were highly related to climate, topography and soil conditions in a single country [25]. To the best of our knowledge, only limited works have studied the relation of growing conditions and the composition of bananas in different countries. In fact, most of the related studies focuses on ways to alter the growing conditions to increase the yield of bananas.
In this preliminary study, the main stable isotope ratios and elemental compositions of banana samples from different origins were characterized with the aim to establish correlations among these fingerprints and the growing conditions of bananas. Banana samples, including pulp and peel, were collected from ten farms located in six different countries in South America. Besides the related growing conditions such as altitude, rainfall and temperatures, some farms cultivated bananas according to an organic and some to a conventional production practices. The stable isotopic ratios and elements of pulp and peel from different countries were firstly characterized through isotope ratio mass spectrometry (IR-MS) and inductively coupled plasma mass spectrometry (ICP-MS). Furthermore, the network between stable isotopic and elemental fingerprints of organic/conventional bananas and their growing conditions were correlated. Due to the limitations in exploring the effects of external (environmental) factors such as harvest year and season together with the use of limited samples numbers, the paper was aimed to provide the feasibility study for the relations between growing conditions and isotopic and elemental fingerprints of bananas.

2. Materials and Methods

2.1. Sample Collection and Growing Conditions

Banana (Musa spp.) samples (AAA Group, Cavendish Subgroup) cultivar (cv.) Williams were collected in a 3-month period from February to April 2018 from ten farms located in six different countries namely, Colombia, Costa Rica, Dominica Republic, Ecuador, Panama, Peru. The locations of the farms, based on their Global Positioning System (GPS) coordinates, are shown in Figure 1. For each farm, twelve bunches from the top position of the banana trees were randomly selected. The bunches were then placed in clean polyethene bags, labelled with the sampling site location and transported in low temperature and dark conditions to Wageningen University & Research. Upon receival of the bunches, two banana fingers were randomly selected from each bunch and combined as one research sample. The sampled bananas were then separated into pulp and peel, cut into slices, freeze-dried and crushed to fine powder. Ultimately, 120 powdered banana pulp and peel samples were collected and stored at −20 °C until IR-MS and ICP-MS analyses. The growing conditions that include the geographical factors and production system data were also obtained from local farms. The overview of the collected banana samples, altitude, monthly mean temperature, annual rainfall and production system data is shown in Table 1.

2.2. Determining the Stable Isotope Ratios

The three kinds of stable isotopes (Carbon, Nitrogen and Oxygen) of bananas were performed according to the protocol by Erasmus et al. with slight changes [26]. The banana pulp and peel samples were analyzed for δ13C and δ15N using an isotope ratio mass spectrometer (Flash 2000 in combination with Delta V advantage, Thermo Scientific, Waltham MA USA). The results were expressed as stable isotope ratios (13C/12C or 15N/14N) against Vienna Pee Dee Belimnite (VPDB) for δ13C and atmospheric AIR for δ15N. The oxygen isotope ratios were performed by a SerCon high-temperature elemental analyzer interfaced with a SerCon 20-20 IR-MS with SysCon electronics (Sercon, Cheshire, UK). The final delta unit (δ18O) was expressed as 18O/16O relative to international standards VSMOW2 (Vienna Standard Mean Ocean Water). For δ13C and δ15N, the samples were weighed into tin capsule, while 4 × 6 mm silver cups were used for δ18O. To reduce the uncertainty and ensure the quality of the measurement, the international atomic energy agency (IAEA) standard samples (IAEA-600 for Carbon, Nitrogen and IAEA-601 for Oxygen) were used. Every 12 samples were interspersed with a standard sample for calibration. All analyses were conducted in triplicate.

2.3. Determining the Elemental Compositions

The 11 kinds of elemental compositions (Al, Ba, Cr, Cu, Fe, Mn, Mo, Ni, Rb, Sr and Zn) of banana pulp and peel were determined by ICP-MS [27]. The freeze-dried and fine homogeneous samples (pulp/peel) were digested using a MARS 6 microwave (CEM Corporation, Matthews NC, USA). About 0.8 g sample was accurately weighed into Teflon vessels. For digestion, 3 mL 70% Nitric acid (HNO3) and 1 mL 30% Hydrogen peroxide (H2O2) were added to each vessel. Besides, another vessel only within HNO3 and H2O2 was selected as a blank group. After digestion, the samples were transferred to tubes and filled up to 10 mL with Milli-Q water. Standard solutions were prepared from 1000 mg/L ICP-MS stock solutions (ICP-MS calibration standard, ULTRA Scientific, North Kingstown, Rhode Island). The multi-elemental analyses were performed using a NexION 300D ICP-MS (Perkin Elmer, Waltham, MA, USA). The analytical performance was verified by processing certified reference materials Lichen 482 (Sigma-Aldrich, Buchs, Switzerland), tuning solution A and B (Merck, Darmstadt, Germany). The concentrations of elemental composition in samples were determined using external calibration curves and the rhodium was used as an internal standard.

2.4. Statistical Analysis

The compositional data were statistically analyzed using a one-way analysis of variance (ANOVA) and Tukey’s test (multiple-range) to assess the impact of the geographical factors and production system of bananas, as a single categorical factor, on the different stable isotopic ratios and elemental compositions. Means with p values below 0.05 were significantly different [28]. The multivariate characterization of the samples concerning stable isotopic ratios and elements, and growing conditions were performed by multiple variable regression and principal component analysis (PCA) [29]. The stable isotopic ratios and elemental compositions were pre-processed using a log transformation before PCA analysis was performed. Pearson correlation coefficients (R) were calculated to evaluate the correlation between growing conditions and banana compositions [30]. Data analyses and the visualizations were conducted using R 4.0.2 (R Foundation for Statistical Computing, Vienna, Austria).

3. Results and Discussion

3.1. Summary of the Growing Conditions of Banana Farms

The geographical factors (altitude, monthly mean temperature and annual rainfall) and production system (organic and conventional cultivation) were regarded as growing conditions of the banana farms. In total, ten farms from six different countries were selected as sampling sites due to their differences in growing conditions between the different farms (Figure 1). Table 1 shows the six conventional farms (CO1, CR1, CR2, CR3, EC2, PA1) and four organic farms (DR1, DR2, EC1, PE1). The farms were specifically selected to explore the effects of geographical factors and production system. For example, the farms with a conventional production system, but from different countries, could be used to compare the effects of geographical origins. Whereas the farms with different production systems, but from the same country, allow one to determine the effect of production systems. Figure 2 indicates the differences in geographical factors such as altitude, temperature and rainfall due to the different locations. The ten farms occupied different spatial locations on the scatterplot composed of altitude, temperature and rainfall. For example, the farm CR2 from Costa Rica had the highest annual rainfall of 5014 mm/year and an altitude of 47 m. Therefore, it grouped separately from Farm PE1 with a rainfall of 200 mm/year and an altitude of 40 m. However, farms EC2, DR1 and DR2 grouped closely together as they share similar geographical factors.

3.2. Differences of Stable Isotopes between Bananas from Different Growing Conditions

For the banana pulp samples, as shown in Table 2, the three stable isotopes ratios (δ13C, δ15N and δ18O) of bananas from different countries in Latin America were compared. For all ten farms, the δ13C of pulp samples ranged from −23.8‰ (CR2) to −22.6‰ (PE1). The δ13C value for PE1 was significantly different to the carbon isotope ratio of CR3, EC2 and PA1. For the δ15N value, the results of DR2 (6.3 ± 1.3‰) was remarkedly higher than the ratios for the other farms. In fact, the δ15N values for the conventional farms (CR1, CR3, EC2 and PA1) were significantly lower than most of the organic farms such as DR1, DR2, EC1 and PE1. The stable isotope ratios of oxygen were the highest for CR1 (31.9 ± 0.7‰), while EC2 had the lowest value at 28.4 ± 0.5‰.
The stable isotope ratios of the peel samples were consistent with that of the pulp samples. Table 3 shows that the values of δ13C fluctuated from −25.6‰ (DR2) to −23.6‰ (DR1 and PE1). The variation among the δ15N values of the samples was greater than that of the δ13C values. The nitrogen ratios of the organic farms were significantly higher than the conventional farms; such as 8.0 ± 0.4‰ for DR2 and 0.8 ± 0.6‰ for PA1. For the δ18O result, the highest value was reported for DR1 (29.9 ± 0.3‰), while the lowest was 25.2 ± 0.5‰ (EC2). The values from different farms changed according to related growing conditions (p < 0.05). With regard to method precision, the standard deviations of repeated pulp samples (n = 12) of 10 groups were ranged from 0.2–0.8‰ for δ13C, 0.1–1.5‰ for δ15N and 0.3–1.1‰ for δ18O. In the peel samples, the standard deviations of repeated measurements (n = 6) were ranged from 0.2–0.8‰ for δ13C, 0.2–0.6‰ for δ15N and 0.3–1.3‰ for δ18O.
There are several studies reporting that the δ15N value could be used to determine if a food product were produced organically [31]. It is well-known that the main type of nitrogen isotope in synthesized fertilizers is isotope 14N, while the proportion of the heavier isotope 15N is much higher than 14N in organic fertilizers [32]. This is because the conventional fertilizers are chemically synthesized by N2 in the air under high temperature and high pressure, the content of 15N in chemical nitrogen fertilizer is equivalent to that in the atmosphere. However, organic fertilizers are derived from animal matter, animal manure, and vegetable matter, etc., therefore, stable isotope 15N is continuously enriched in plant and animals due to bioconcentration effects. Due to the enrichment of 15N, the abundance of 15N in plants or animals is significantly higher than that in air. Therefore, the abundance of 15N in organic fertilizers is significantly higher than that in chemical fertilizers [33]. In line with the literature, the pulp and peel samples from organic farms such as DR1, DR2, EC1 and PE1 had generally higher δ15N (15N/14N) values than the samples from some of the other conventional farms. The highest δ15N was 6.3‰ and 8.0‰ in pulp and peel samples of farm DR2, respectively. Most of the conventional farms were reported having δ15N values lower than 1.0‰ except for CO1 with a nitrogen isotope ratio of 2.9‰ for pulp and 3.8‰ for the peel samples. This could be due to the fact that some farms choose to apply chemical and organic fertilizers simultaneously to improve soil vitality and agricultural sustainability [34]. Therefore, the application of organic fertilizer could also improve the proportion of the stable isotope 15N in the conventionally produced bananas of farm CO1.

3.3. Differences of Elemental Compositions between Bananas from Different Growing Conditions

More than 20 types of minerals were detected in the bananas by ICP-MS. In total, 10 elements in pulp and 11 elements in peel samples were reported higher than the limit of detection and hence, only these were selected to compare the effects of geographical factors (altitude, monthly mean temperature and annual rainfall) and production system (organic and conventional cultivation). With regard to method precision, the relative standard deviations of pulp samples among all groups for elements, Fe were less than 10%, Cu, Ni, Sr and Zn were less than 20% and Al, Ba, Mn, Mo, Rb were less than 30%. In the peel samples, slightly higher standard deviations were reported, with Ba, Cu, Fe, Mo, Sr, Zn less than 20%, and Al, Cr, Mn, Ni, Rb were ranged from 20% to 40%. The concentrations of elements for the banana pulp and peel samples obtained from the different countries are listed in Table 2 and Table 3 and Tables S1 and S2 in supplementary materials. For the pulp samples, 10 elements (Al, Ba, Cu, Fe, Mn, Mo, Ni, Rb, Sr and Zn) were detected in the samples. The mean concentration of the elements was significantly different among farms. The concentration of Mn was reported as 21.9 mg/kg for CR1, which was higher (p < 0.05) than the concentration for DR1 (5.2 mg/kg). The pulp samples from the organic farms, DR1 and DR2, contained a significantly higher concentration of Zn. For the peel samples, 11 elements (Al, Ba, Cr, Cu, Fe, Mn, Mo, Ni, Rb, Sr and Zn) were obtained and were significantly different according to geographical factors and the production system. The elemental profiles of the banana peel samples were quite different to that of the pulp samples; with the concentration of Al, Fe, Mn, Sr, and Zn being higher in the peel samples. The organic farm, PE1, had a notably lower concentration of Fe when compared to other conventional farms (p < 0.05). Most of the farms have significantly different concentrations of elements in the peel samples, but there are still a few farms that did not show any significant differences in some elements (p > 0.05). As seen in Table 2 and Table 3, the concentrations of the elements of the pulp and peel samples were significantly different among most of the farms. The results were consistent with the findings of a recent paper about the elemental compositions of bananas [35]. Hardisson et al. detected Fe, Cu, Zn, Mn and B in the banana pulp of Musa AAA bananas, the same variety used in the current study [36], which is in line with the current results. As reported in comparative research of eight banana cultivars, the trace elements such as Zn, Fe and Mn were also detected in peel samples [37].
The current study systematically described and compared the trace elements of banana pulp and peel from different growing conditions. The ANOVA results indicates significant differences of elements for most of pulp and peel samples. For example, for the pulp samples, the farm DR2 showed the lowest concentration of Al at 5.4 mg/kg, significantly lower than farm DR1 at 23.5 mg/kg, even though both farms are in the same country. Therefore, it is vital to also consider the geographical (or external) factors as opposed to just the “country border” when studies on the geographical origin of foods are performed. For example, the two geographically distant farms, DR2 (organic) and EC2 (conventional), have similar geographical factors in relation to altitude, temperature and rainfall; the scatter plot (Figure 2) showed that the space location of both farms were almost coincident. Likewise, the concentration of most of the elements such as Ba, Cr, Cu, Fe, Mn, Mo, Ni, Rb and Zn of the peel samples were not significantly different, although these two farms have different production systems. Hence, there are likely other external factors that have a greater influence on the variation in elemental compositions of bananas.

3.4. Principal Component Analysis of Bananas between the Different Growing Condition

To explore the effects of growing conditions on bananas compositions, the IR-MS and ICP-MS data sets were subjected to multivariate data analysis. The PCA was used to show the distribution and groupings of the pulp and peel samples (Figure 3). The PCA plot shows that most of the farms could be separated according to differences in the stable isotopic ratios and elemental compositions of the samples. For the pulp and peel samples, PC1 and PC2 explained 39.6% and 41.4%, respectively, of the total variance. Some farms from the same countries, like CR1, CR2, and CR3 from Costa Rica and DR1, DR2 from the Dominican Republic, overlapped in the PCA plots. Figure 3a,b show the farm distribution according to different geographical factors and production system, respectively. DR1 and DR2 were clearly separated from the other farms (Figure 3a). The differences in the concentrations of Mo, N, Sr, Rb and Ni mainly contributed to the farms’ distributions at geographical level. The effects of production system showed that the variables δ15N, Mo and Sr were characteristics of organic farms, while Mn, Ni, Al, Ba and Fe were characteristics of conventional farms for the pulp samples (Figure 3c).
The stable isotope ratios and concentration of elements of the peel samples resulted in a different distribution of farms compared to the pulp samples. For instance, the concentration of Sr and Cu show a greater contribution to the discrimination of the peel samples based on geographical origin and production system. As seen in Figure 3e, the peel samples from DR1 and DR2 were grouped separately from samples from Central and South America mainland, with only a slight overlap seen between the different farms. Whereas Figure 3a,c show that the internal variation between peel samples within each farms is lower than their variation in pulp samples, Figure 3c shows that the variables δ15N, Mo, Sr and Cr were characteristics of organic farms, while Ni, Ba, Mn, Cu, Al, Rb and Zn were characteristics of conventional farms for the peel samples.
There is considerable evidence that the application of stable isotopes and elements is of great interest for tracing the geographical origin and production system of many food products, such as tea [38], maca [39], orange [40], etc. The use of multi-isotopes (i.e., δ13C, δ15N) and multi-elements (e.g., Na, Fe, Zn and K) have been reported as reliable tools in authenticity testing of organically grown rice [41]. In the present study, the combination of stable isotopic ratios and elemental compositions were employed for the visualization of multifactorial analysis. The PCA loadings and the contribution of variances are shown in Figure 3c,d. A clear separation was identified between organic and conventional production for the pulp and peel samples. The most important variables driving the separation between groups were δ13C, Mo, N, Fe, Ba, Mn and Ni for the pulp samples, and δ18O, δ15N, Se, Mo, Cu, Al, Mn and Ba for the peel samples. The different contribution of elements in the pulp and peel sample could be because the distribution of minerals in the different parts of the fruits were different. For the results on the geographical origin of the bananas (Figure 3e), the peel samples indicated clearer grouping compared to pulp samples, which could be caused by the lower internal variation between peel samples in the same farm compared with the variation in the pulp samples. The PCA separation was consistent with the results of the ANOVA test of the peel samples, while the above-mentioned elements and stable ratios were also significantly different among the 10 farms (Figure 3d). As seen from Table 2, the farms which showed big differences in the scatter plot (Figure 2), using altitude, temperature and rainfall, also grouped separately in the PCA plots. The relation between the separation in “geographical factors” and “PCA groups” could possibly indicate the effects of growing conditions on the stable isotopic ratios and elemental compositions of bananas. Figure 3 showed that organic farms DR1 and DR2 were separated from other organic farms (EC1 and PE1) in both of pulp and peel results, which could be due to the differences in geographical factors (or other external factors) such as altitude, monthly mean temperature and annual rainfall. For example, the monthly mean temperature of DR1 and DR2 were higher than EC1 and PE1, also the rainfall of DR farms was higher than PE1 and lower than EC1.
At the same time, some values of stable isotopes and elements of organic farms were similar when compared with conventional farms. PCA plots also showed that conventional farms and organic farms were completely separated, such as EC1 and EC2 (Figure 3). This phenomenon could be because organic fertilizers were also used in conventional farms for sustainable banana farming [42], which could explain why some organic farms, such as PE1, were plotted closely to conventional farm such as EC2 in PCA results. Overall, the differences in stable isotope ratios and elemental compositions of bananas from selected countries were the result of the combined effects of the geographical factors and production systems.

3.5. Correlation Analysis of Banana Compositions and Their Growing Conditions

To protect the integrity of food and the interests of consumers, most of the related works of literature attempt to discriminate the origin of foods based on the geographical distances such as the “national or provincial boundaries”, whereas only a few papers have been published concerning the underlying mechanisms linked to geographical origin. However, both macromolecular nutrients and micro-metabolites of foods are influenced by the environmental conditions of the origin, planting methods, temperature, soil conditions, etc. [43,44]. Therefore, exploring the correlations between the authentic growing conditions and the composition of foods could providing more in-depth evidence for its traceability. The banana is a suitable fruit to use for correlation research between growing conditions and compositional attributes, as most banana plantations are only located in tropical and subtropical regions, while being sold worldwide. The fresh banana is the most important fruit traded internationally. Verifying the geographical origin of bananas is important to protect the consumers’ pocket and farmers’ benefits [5]. At the same time, identifying and tracing organic bananas (i.e., production system) could further ensure the sustainable supply of bananas and protect consumers’ right to know what they are buying.
To further explain the relationships between the stable isotopes and elements of different sampling farms and the local growing conditions, the compositions and their growing conditions were correlated using Pearson’s correlation coefficients (R). The relationship between growing conditions and banana compositions were shown using the network in Figure 4. The distance between the attributes indicated the strength of the correlation coefficient, while the blue or red colors were used to show if they were positive or negative, respectively.
For banana pulp, a slight negative relationship was found in δ13C and rainfall (R = −0.32, p = 0.18 × 10−3), which is in agreement with the several studies where δ13C of fruits negatively correlate with the rainfall during sugar accumulation in fruit body [45]. δ18O of wine was reported to be highly correlated with daily temperatures and local precipitation [46], and the network analysis of banana pulp showed that δ18O was slightly positively correlated with altitude (R = 0.33, p = 0. 14 × 10−3) and rainfall (R = 0.33, p = 0. 12 × 10−3) (Figure 4a). Table 1 and Table 2 also showed that farm CR1 had the higher altitude and rainfall values than CO1, at the same time, significant higher δ18O were detected in pulp samples from CR1. As for stable isotope N, only a weak positive correlation with temperature and δ15N (R = 0.32, p < 0.01) were observed. On the other hand, the element Fe indicated a strong correlation with rainfall (R = 0.74, p = 0.01), a significant and moderate correlation was found between altitude and two elements: Mn (R = 0.66, p = 0.01) and Ni (R = 0.53, p < 0.01). A negative correlation with Fe (R = −0.38, p < 0.01) and Mn (R = −0.36, p < 0.01) due to the effects of temperatures was also found in network analysis. Generally, the differences among sampling locations were mainly caused by individual geographical features such as the altitude, temperature and rainfall. The growth of bananas requires sufficient humidity, light and moisture, otherwise the yield and quality of bananas would be decreased [47].
For banana peel, Figure 4b shows that most of the compositional attributes of peel samples were associated with geographical factors. The differences of δ13C from different farms could be associated with temperature conditions (R = −0.41, p = 0.90 × 10−3), while the variation of δ15N was slightly correlated with altitude (R = −0.36, p = 0.42 × 10−2) and rainfall (R = −0.42, p = 0.77 × 10−3). It could suggest that other factors likely also contribute to the variations of stable isotope N from different farms. Most of the research indicates that organic production tends to increase the values of δ15N more than conventional cultivation [48]. The PCA plots also confirmed that organic farms could have higher values of δ15N. The slight negative relationship was found between rainfall and δ13C (R = −0.32, p = 0.18 × 10−3). Similarly, the elements Mn (R = 0.65, p < 0.01) and Ni (R = 0.41, p = 0.94 × 10−3) showed significantly positive correlations with altitude. A slight positive relationship was found in rainfall and Ba (R = 0.44, p = 0.04 × 10−3).
The significant correlation between growing conditions of the farms and the compositional attributes (stable isotopic ratios and elemental concentrations) was an obvious indication of the relationship between growing conditions and banana compositions. The underlying reason may be that geographical factors (e.g., altitude, temperature, rainfall) affected the absorption capacity of banana roots to elements, thereby promoting or inhibiting the concentration of elements in pulp and peel. Heat stress could induce various physiological injuries to plants such as root growth inhibition, therefore the absorption of elements from soil to roots will be inhibited as well [49]. Subsequently, a negative correlation between temperature and most of the elements of banana pulp and peel were observed in this study. Especially for farm DR1 and DR2 with higher monthly mean temperature than most farms, the elemental contents like Ba, Cu, Fe, Mn, Ni and Cu of pulp samples were also lower than other farms correspondingly. In the effects of altitude on fruit compositions, Crespo et al. reported that chemical composition such as vitamin C, organic acids and anthocyanins of strawberry changed due to increasing altitude [50]. The possible reason could be that the microclimate in a high-altitude region changed a lot in air temperature, precipitation and evaporation compared with a low-altitude region [51]. As found in this study, the altitude also correlated with Mn, Mo, Ni positively or negatively of pulp and peel samples. However, not all stable isotope and elements correlated with the growing conditions, which means more geographical factors such as soil type, sunshine data, etc., should be collected to strengthen the findings.

4. Conclusions

A comprehensive application of IR-MS and ICP-MS was employed to identify the isotopes and mineral elements fingerprints of Cavendish bananas from different origins. The ANOVA results demonstrated that the stable isotopic ratios and elemental concentrations changed significantly according to different farm locations. Based on the obtained PCA plots, the selected farms could be separated according to different geographical factors and production systems. This revealed that the variation among banana compositions were potentially related to their growing conditions. Moreover, the δ15N values of pulp and peel were significantly higher in the organic system than the conventional system. On the other hand, relations exist between the geographical factors and the concentration of elements in bananas; particularly, Mn and Ni were significantly positively correlated with altitude. The underlying reason could be that the interaction of roots and mineral elements and the banana metabolism is affected by the growing conditions. Overall, the combination of stable isotopes and elemental compositions could be used to explore the effects of local growing conditions concerning origin. However, this feasibility study was limited to one single season and a low number of samples. It is important that more comprehensive discussions on the effects of environmental factors are strengthened in future work. In the following research, a soil map based on correlation of banana composition and growing conditions will be generated to predict the compositional features of bananas without huge chemical database.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/foods10051021/s1, Table S1: The minimum and maximum values of the stable isotopes ratios (13C/12C, 15N/14N and 18O/16O) and elemental compositions (mg/kg) of banana pulp samples (n = 12). Table S2: The minimum and maximum values of the stable isotopes ratios (13C/12C, 15N/14N and 18O/16O) and elemental compositions (mg/kg) of banana pulp samples (n = 6).

Author Contributions

Conceptualization, Z.W.; data curation, Z.W.; formal analysis, Z.W. and S.W.E.; funding acquisition, S.M.v.R.; investigation, Z.W.; methodology, Z.W.; resources, S.M.v.R.; software, Z.W. and S.W.E.; visualization, Z.W. and S.W.E.; writing—original draft, Z.W.; writing—review and editing, S.W.E. and S.M.v.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Dutch Top sector Agri & Food project, grant number AF-16008 and the first author Zhijun Wang received a PhD scholarship from the China Scholarship Council (grant number 201706820017).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not available.

Acknowledgments

We would like to thank Erwin Brouwer from Wageningen Food Safety Research for obtaining banana samples. We also hope to thank Greet van Bemmel from Wageningen Food Safety Research and Michel Breuer from Animal Nutrition group in Wageningen University for their outstanding technical guidance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Oyeyinka, B.O.; Afolayan, A.J. Comparative Evaluation of the Nutritive, Mineral, and Antinutritive Composition of Musa sinensis L. (Banana) and Musa paradisiaca L. (Plantain) Fruit Compartments. Plants 2019, 8, 598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Aurore, G.; Parfait, B.; Fahrasmane, L. Bananas, raw materials for making processed food products. Trends Food Sci. Technol. 2009, 20, 78–91. [Google Scholar] [CrossRef]
  3. Dale, J.; James, A.; Paul, J.-Y.; Khanna, H.; Smith, M.; Peraza-Echeverria, S.; Garcia-Bastidas, F.; Kema, G.; Waterhouse, P.; Mengersen, K.; et al. Transgenic Cavendish bananas with resistance to Fusarium wilt tropical race 4. Nat. Commun. 2017, 8, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Pereira, A.; Maraschin, M. Banana (Musa spp.) from peel to pulp: Ethnopharmacology, source of bioactive compounds and its relevance for human health. J. Ethnopharmacol. 2015, 160, 149–163. [Google Scholar] [CrossRef]
  5. FAO. Banana Market Review 2019; FAO: Roma, Italy, 2020. [Google Scholar]
  6. Aprile, M.C.; Caputo, V.; Nayga, R.M.N., Jr. Consumers’ valuation of food quality labels: The case of the European geographic indication and organic farming labels. Int. J. Consum. Stud. 2012, 36, 158–165. [Google Scholar] [CrossRef]
  7. Bellamy, A.S. Banana Production Systems: Identification of Alternative Systems for More Sustainable Production. Ambio 2013, 42, 334–343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Santeramo, F.G.; Lamonaca, E. Evaluation of geographical label in consumers’ decision-making process: A systematic review and meta-analysis. Food Res. Int. 2020, 131, 108995. [Google Scholar] [CrossRef] [Green Version]
  9. Sun, S.; Guo, B.; Wei, Y. Origin assignment by multi-element stable isotopes of lamb tissues. Food Chem. 2016, 213, 675–681. [Google Scholar] [CrossRef]
  10. Shen, M.; Zhao, S.; Zhang, F.; Huang, M.; Xie, J. Characterization and authentication of olive, camellia and other vegetable oils by combination of chromatographic and chemometric techniques: Role of fatty acids, tocopherols, sterols and squalene. Eur. Food Res. Technol. 2021, 247, 411–426. [Google Scholar] [CrossRef]
  11. Katerinopoulou, K.; Kontogeorgos, A.; Salmas, C.E.; Patakas, A.; Ladavos, A. Geographical Origin Authentication of Agri-Food Products: A Review. Foods 2020, 9, 489. [Google Scholar] [CrossRef]
  12. Robson, K.; Dean, M.; Haughey, S.; Elliott, C. A comprehensive review of food fraud terminologies and food fraud mitigation guides. Food Control 2021, 120, 107516. [Google Scholar] [CrossRef]
  13. Van Ruth, S.M.; Luning, P.A.; Silvis, I.C.J.; Yang, Y.; Huisman, W. Differences in fraud vulnerability in various food supply chains and their tiers. Food Control 2018, 84, 375–381. [Google Scholar] [CrossRef]
  14. Van Ruth, S.M.; Huisman, W.; Luning, P.A. Food fraud vulnerability and its key factors. Trends Food Sci. Technol. 2017, 67, 70–75. [Google Scholar] [CrossRef]
  15. Bronzi, B.; Brilli, C.; Beone, G.M.; Fontanella, M.C.; Ballabio, D.; Todeschini, R.; Consonni, V.; Grisoni, F.; Parri, F.; Buscema, M. Geographical identification of Chianti red wine based on ICP-MS element composition. Food Chem. 2020, 315, 126248. [Google Scholar] [CrossRef]
  16. Zhao, Y.; Tu, T.; Tang, X.; Zhao, S.; Qie, M.; Chen, A.; Yang, S. Authentication of organic pork and identification of geographical origins of pork in four regions of China by combined analysis of stable isotopes and multi-elements. Meat Sci. 2020, 165, 108129. [Google Scholar] [CrossRef]
  17. Magdas, D.A.; Guyon, F.; Puscas, R.; Vigouroux, A.; Gaillard, L.; Dehelean, A.; Feher, I.; Cristea, G. Applications of emerging stable isotopes and elemental markers for geographical and varietal recognition of Romanian and French honeys. Food Chem. 2021, 334, 127599. [Google Scholar] [CrossRef]
  18. Baskali-Bouregaa, N.; Milliand, M.-L.; Mauffrey, S.; Chabert, E.; Forrestier, M.; Gilon, N. Tea geographical origin explained by LIBS elemental profile combined to isotopic information. Talanta 2020, 211, 120674. [Google Scholar] [CrossRef]
  19. Wang, J.; Chen, T.; Zhang, W.; Zhao, Y.; Yang, S.; Chen, A. Tracing the geographical origin of rice by stable isotopic analyses combined with chemometrics. Food Chem. 2020, 313, 126093. [Google Scholar] [CrossRef]
  20. De Rijke, E.; Schoorl, J.C.; Cerli, C.; Vonhof, H.B.; Verdegaal, S.J.A.; Vivó-Truyols, G.; Lopatka, M.; Dekter, R.; Bakker, D.; Sjerps, M.J.; et al. The use of δ2H and δ18O isotopic analyses combined with chemometrics as a traceability tool for the geographical origin of bell peppers. Food Chem. 2016, 204, 122–128. [Google Scholar] [CrossRef] [Green Version]
  21. Geana, I.; Iordache, A.; Ionete, R.; Marinescu, A.; Ranca, A.; Culea, M. Geographical origin identification of Romanian wines by ICP-MS elemental analysis. Food Chem. 2013, 138, 1125–1134. [Google Scholar] [CrossRef]
  22. Bigot, C.; Bugaud, C.; Camilo, J.; Kapitan, A.; Montet, D.; Meile, J.-C. Impact of farming type, variety and geographical origin on bananas bacterial community. Food Control 2020, 109, 106925. [Google Scholar] [CrossRef]
  23. Choi, S.-H.; Bong, Y.-S.; Park, J.H.; Lee, K.-S. Geographical origin identification of garlic cultivated in Korea using isotopic and multi-elemental analyses. Food Control 2020, 111, 107064. [Google Scholar] [CrossRef]
  24. Rees, G.; Kelly, S.D.; Cairns, P.; Ueckermann, H.; Hoelzl, S.; Rossmann, A.; Scotter, M.J. Verifying the geographical origin of poultry: The application of stable isotope and trace element (SITE) analysis. Food Control 2016, 67, 144–154. [Google Scholar] [CrossRef]
  25. Wang, Z.; Erasmus, S.W.; Dekker, P.; Guo, B.; Stoorvogel, J.J.; Van Ruth, S.M. Linking growing conditions to stable isotope ratios and elemental compositions of Costa Rican bananas (Musa spp.). Food Res. Int. 2020, 129, 108882. [Google Scholar] [CrossRef] [PubMed]
  26. Erasmus, S.W.; Muller, M.; van der Rijst, M.; Hoffman, L.C. Stable isotope ratio analysis: A potential analytical tool for the authentication of South African lamb meat. Food Chem. 2016, 192, 997–1005. [Google Scholar] [CrossRef]
  27. Herwig, N.; Stephan, K.; Panne, U.; Pritzkow, W.; Vogl, J. Multi-element screening in milk and feed by SF-ICP-MS. Food Chem. 2011, 124, 1223–1230. [Google Scholar] [CrossRef]
  28. Yan, J.; Oey, S.B.; van Leeuwen, S.P.J.; van Ruth, S.M. Discrimination of processing grades of olive oil and other vegetable oils by monochloropropanediol esters and glycidyl esters. Food Chem. 2018, 248, 93–100. [Google Scholar] [CrossRef]
  29. Vera, D.N.; Jiménez-Carvelo, A.M.; Cuadros-Rodríguez, L.; Ruisánchez, I.; Callao, M.P. Authentication of the geographical origin of extra-virgin olive oil of the Arbequina cultivar by chromatographic fingerprinting and chemometrics. Talanta 2019, 203, 194–202. [Google Scholar] [CrossRef] [PubMed]
  30. Liu, N.; Pustjens, A.M.; Erasmus, S.W.; Yang, Y.; Hettinga, K.; van Ruth, S.M. Dairy farming system markers: The correlation of forage and milk fatty acid profiles from organic, pasture and conventional systems in the Netherlands. Food Chem. 2020, 314, 126153. [Google Scholar] [CrossRef]
  31. Laursen, K.H.; Schjoerring, J.K.; Kelly, S.D.; Husted, S. Authentication of organically grown plants—advantages and limitations of atomic spectroscopy for multi-element and stable isotope analysis. TrAC Trends Anal. Chem. 2014, 59, 73–82. [Google Scholar] [CrossRef]
  32. Laursen, K.H.; Mihailova, A.; Kelly, S.D.; Epov, V.N.; Bérail, S.; Schjoerring, J.K.; Donard, O.F.X.; Larsen, E.H.; Pedentchouk, N.; Marca-Bell, A.D.; et al. Is it really organic? – Multi-isotopic analysis as a tool to discriminate between organic and conventional plants. Food Chem. 2013, 141, 2812–2820. [Google Scholar] [CrossRef] [PubMed]
  33. Camin, F.; Bontempo, L.; Perini, M.; Piasentier, E. Stable Isotope Ratio Analysis for Assessing the Authenticity of Food of Animal Origin. Compr. Rev. Food Sci. Food Saf. 2016, 15, 868–877. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Mustaffa, M.M.; Kumar, V. Banana production and productivity enhancement through spatial, water and nutrient manage-ment. J. Hortic. Sci. 2012, 7, 1–28. [Google Scholar]
  35. Emaga, T.H.; Andrianaivo, R.H.; Wathelet, B.; Tchango, J.T.; Paquot, M. Effects of the stage of maturation and varieties on the chemical composition of banana and plantain peels. Food Chem. 2007, 103, 590–600. [Google Scholar] [CrossRef]
  36. Hardisson, A.; Rubio, C.; Baez, A.; Martin, M.; Alvarez, R.; Diaz, E. Mineral composition of the banana (Musa acuminata) from the island of Tenerife. Food Chem. 2001, 73, 153–161. [Google Scholar] [CrossRef]
  37. Sulaiman, S.F.; Yusoff, N.A.M.; Eldeen, I.M.; Seow, E.M.; Sajak, A.A.B.; Ooi, K.L. Correlation between total phenolic and mineral contents with antioxidant activity of eight Malaysian bananas (Musa sp.). J. Food Compos. Anal. 2011, 24, 1–10. [Google Scholar] [CrossRef]
  38. Peng, C.-Y.; Zhang, Y.-L.; Song, W.; Lv, Y.-N.; Xu, Q.; Zheng, P.; Zhang, Z.-Z.; Wan, X.-C.; Hou, R.-Y.; Cai, H.-M.; et al. Using stable isotope signatures to delineate the geographic point-of-origin of Keemun black tea. J. Sci. Food Agric. 2019, 99, 2596–2601. [Google Scholar] [CrossRef] [PubMed]
  39. He, Y.; Sun, Q.; Zhang, X.; Bao, X.; Wang, Y.; Rasheed, M.; Guo, B. Authentication of the geographical origin of Maca (Lepidium meyenii Walp.) at different regional scales using the stable isotope ratio and mineral elemental fingerprints. Food Chem. 2020, 311, 126058. [Google Scholar] [CrossRef]
  40. Rapisarda, P.; Camin, F.; Fabroni, S.; Perini, M.; Torrisi, B.; Intrigliolo, F. Influence of Different Organic Fertilizers on Quality Parameters and the δ15N, δ13C, δ2H, δ34S, and δ18O Values of Orange Fruit (Citrus sinensis L. Osbeck). J. Agric. Food Chem. 2010, 58, 3502–3506. [Google Scholar] [CrossRef]
  41. Liu, Z.; Yuan, Y.; Xie, T.; Zhang, Y.; Shao, S.; Nie, J.; Xia, W.; Rogers, K.M.; Zhang, W. Long-Term Agricultural Effects on the Authentication Accuracy of Organic, Green, and Conventional Rice Using Isotopic and Elemental Chemometric Analyses. J. Agric. Food Chem. 2020, 68, 1213–1225. [Google Scholar] [CrossRef]
  42. Gonçalves, A.L.; Kernaghan, J.R. Banana production methods. 2014. Available online: https://www.naturskyddsforeningen.se/sites/default/files/dokument-media/banana_report_final_version.pdf (accessed on 4 July 2014).
  43. Sim, I.; Suh, D.H.; Singh, D.; Do, S.-G.; Moon, K.H.; Lee, J.H.; Ku, K.-M.; Lee, C.H. Unraveling Metabolic Variation for Blueberry and Chokeberry Cultivars Harvested from Different Geo-Climatic Regions in Korea. J. Agric. Food Chem. 2017, 65, 9031–9040. [Google Scholar] [CrossRef] [PubMed]
  44. Cao, T.-T.; Liu, J.; Zhang, X.Y.; Wei, X.P.; Qi, Y.; Zhang, B.; Liu, H.; Xiao, P.-G. Metabolomics characterization of different geographical origins of Flos Carthami using UPLC-QTOF/MS and their correlation with climate factors. Anal. Methods 2019, 11, 2700–2711. [Google Scholar] [CrossRef]
  45. Camin, F.; Dordevic, N.; Wehrens, R.; Neteler, M.; Delucchi, L.; Postma, G.; Buydens, L. Climatic and geographical dependence of the H, C and O stable isotope ratios of Italian wine. Anal. Chim. Acta 2015, 853, 384–390. [Google Scholar] [CrossRef]
  46. Martin, G.J.; Martin, M.L. Climatic significance of isotope ratios. Phytochem. Rev. 2003, 2, 179–190. [Google Scholar] [CrossRef]
  47. Coltro, L.; Karaski, T.U. Environmental indicators of banana production in Brazil: Cavendish and Prata varieties. J. Clean. Prod. 2019, 207, 363–378. [Google Scholar] [CrossRef]
  48. Nishida, M.; Sato, Y. Characteristics of the Relationship between Natural 15N Abundances in Organic Rice and Soil. Plant Prod. Sci. 2015, 18, 180–186. [Google Scholar] [CrossRef] [Green Version]
  49. Bita, C.E.; Gerats, T. Plant tolerance to high temperature in a changing environment: Scientific fundamentals and production of heat stress-tolerant crops. Front. Plant Sci. 2013, 4, 273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Crespo, P.; Bordonaba, J.G.; Terry, L.A.; Carlen, C. Characterisation of major taste and health-related compounds of four strawberry genotypes grown at different Swiss production sites. Food Chem. 2010, 122, 16–24. [Google Scholar] [CrossRef]
  51. Bennie, J.; Huntley, B.; Wiltshire, A.; Hill, M.O.; Baxter, R. Slope, aspect and climate: Spatially explicit and implicit models of topographic microclimate in chalk grassland. Ecol. Model. 2008, 216, 47–59. [Google Scholar] [CrossRef]
Figure 1. The sampling sites of the collected bananas. CO: Colombia; CR: Costa Rica; DR: Dominican Republic; EC: Ecuador; PA: Panama; PE: Peru; the numbers refer to the individual farms in a specific country.
Figure 1. The sampling sites of the collected bananas. CO: Colombia; CR: Costa Rica; DR: Dominican Republic; EC: Ecuador; PA: Panama; PE: Peru; the numbers refer to the individual farms in a specific country.
Foods 10 01021 g001
Figure 2. Differences in geographical factors of banana farms from Colombia (CO), Costa Rica (CR), Dominican Republic (DR), Ecuador (EC), Panama (PA), and Peru (PE). The numbers refer to individual farms in a country (please refer to the web version of this article for the interpretation of the colors used in the figure).
Figure 2. Differences in geographical factors of banana farms from Colombia (CO), Costa Rica (CR), Dominican Republic (DR), Ecuador (EC), Panama (PA), and Peru (PE). The numbers refer to individual farms in a country (please refer to the web version of this article for the interpretation of the colors used in the figure).
Foods 10 01021 g002
Figure 3. The principal component analysis (PCA) of banana pulp (n = 12) and peel (n = 6) samples according to stable isotopic ratios and elemental compositions. (a) PCA plot showing scores for the pulp samples from the different geographical origins and loadings for the variables. Explained variance: PC1 = 22.8%, PC2 = 16.8%; (b) PCA plot showing scores for the pulp samples from the different production systems and loadings for the variables. Explained variance: PC1 = 22.8%, PC2 = 16.8%; (c) the loading plot of stable isotopic ratios and elemental compositions for banana pulp samples; (d) the loading plot of stable isotopic ratios and elemental compositions for banana peel samples; (e) PCA plot showing scores for the peel samples from the different geographical origins and loadings for the variables. Explained variance: PC1 = 21.3%, PC2 = 20.1%; (f) PCA plot showing scores for the peel samples from the different production systems and loadings for the variables. Explained variance: PC1 = 21.3%, PC2 = 20.1%.CO: Colombia; CR: Costa Rica; DR: Dominican Republic; EC: Ecuador; PA: Panama; PE: Peru (please refer to the web version of this article for the interpretation of the colors used in the figure).
Figure 3. The principal component analysis (PCA) of banana pulp (n = 12) and peel (n = 6) samples according to stable isotopic ratios and elemental compositions. (a) PCA plot showing scores for the pulp samples from the different geographical origins and loadings for the variables. Explained variance: PC1 = 22.8%, PC2 = 16.8%; (b) PCA plot showing scores for the pulp samples from the different production systems and loadings for the variables. Explained variance: PC1 = 22.8%, PC2 = 16.8%; (c) the loading plot of stable isotopic ratios and elemental compositions for banana pulp samples; (d) the loading plot of stable isotopic ratios and elemental compositions for banana peel samples; (e) PCA plot showing scores for the peel samples from the different geographical origins and loadings for the variables. Explained variance: PC1 = 21.3%, PC2 = 20.1%; (f) PCA plot showing scores for the peel samples from the different production systems and loadings for the variables. Explained variance: PC1 = 21.3%, PC2 = 20.1%.CO: Colombia; CR: Costa Rica; DR: Dominican Republic; EC: Ecuador; PA: Panama; PE: Peru (please refer to the web version of this article for the interpretation of the colors used in the figure).
Foods 10 01021 g003
Figure 4. The correlation network of stable isotopes, elements and growing conditions of banana pulp (a) and peel (b) (refer to the web version of this article for the interpretation of the colors used in the figure).
Figure 4. The correlation network of stable isotopes, elements and growing conditions of banana pulp (a) and peel (b) (refer to the web version of this article for the interpretation of the colors used in the figure).
Foods 10 01021 g004
Table 1. The banana pulp and peel samples collected from different countries and related growing conditions.
Table 1. The banana pulp and peel samples collected from different countries and related growing conditions.
CountryFarmPulp (n)Peel (n)Altitude (m)Monthly Mean Temperature (°C)Annual Rainfall (mm/Year)Production System
ColombiaCO11266623.21837Conventional
Costa RicaCR112672623.42857Conventional
CR21264724.45014Conventional
CR31262426.31837Conventional
Dominica RepublicDR11266526.7925Organic
DR21262726.7925Organic
EcuadorEC11263222.91511Organic
EC21262226.5843Conventional
PanamaPA11262019.73679Conventional
PeruPE11264024.1200Organic
n: number of samples.
Table 2. The analysis of variance (ANOVA) of the stable isotopes ratios (13C/12C, 15N/14N and 18O/16O) and elemental compositions (mg/kg) of banana pulp samples (n = 12).
Table 2. The analysis of variance (ANOVA) of the stable isotopes ratios (13C/12C, 15N/14N and 18O/16O) and elemental compositions (mg/kg) of banana pulp samples (n = 12).
Farmδ13Cδ15Nδ18OAlBaCuFeMnMoNiRbSrZn
CO1−23.0 ab ± 0.82.9 b ± 0.830.9 bc ± 0.50.8 bcd ± 0.21.6 a ± 0.23.6 cd ± 0.511.5 b ± 1.315.7 b ± 3.80.1 cd ± 0.10.3 a ± 0.11.9 c ± 0.32.6 b ± 0.46.1 bcd ± 1.1
CR1−23.3 ac ± 0.5−0.9 f ± 0.931.9 a ± 0.70.9 ac ± 0.30.5 b ± 0.13.9 bcd ± 0.610.7 b ± 0.921.9 a ± 5.50.0 d ± 0.00.3 a ± 0.11.3 c ± 0.32.0 bc ± 0.55.2 de ± 0.8
CR2−23.8 ac ± 0.61.1 cd ± 1.531.4 ab ± 0.61.2 a ± 0.51.6 a ± 0.34.5 ab ± 0.412.7 a ± 0.512.9 bc ± 1.90.2 b ± 0.00.1 cd ± 0.09.0 b ± 1.12.5 b ± 0.65.9 bce ± 1.2
CR3−23.5 bc ± 0.4−0.5 ef ± 0.631.2 ac ± 0.40.9 bc ± 0.10.4 bc ± 0.03.7 cd ± 0.210.7 b ± 0.69.2 ce ± 2.90.2 b ± 0.10.2 cd ± 0.07.0 b ± 2.01.9 bc ± 0.25.9 bce ± 0.5
DR1−23.3 ac ± 0.51.9 cd ± 0.430.5 c ± 0.30.6 cd ± 0.20.3 cd ± 0.04.2 ac ± 0.78.8 c ± 0.75.2 f ± 1.70.3 a ± 0.10.1 cd ± 0.01.9 c ± 0.83.3 a ± 0.67.0 ab ± 0.6
DR2−23.2 ac ± 0.66.3 a ± 1.330.9 bc ± 1.10.6 cd ± 0.30.5 bc ± 0.24.5 bd ± 0.89.3 c ± 0.97.6 ef ± 2.50.2 b ± 0.10.1 d ± 0.02.9 c ± 1.83.5 a ± 0.67.1 a ± 1.1
EC1−23.1 ac ± 0.61.5 ce ± 0.629.1 d ± 0.70.7 cd ± 0.20.2 d ± 0.12.8 e ± 0.28.6 c ± 0.55.6 ef ± 0.80.1 d ± 0.00.2 bc ± 0.07.3 b ± 2.61.2 de ± 0.46.0 bce ± 0.4
EC2−23.6 bc ± 0.2−0.4ef ± 0.228.4 d ± 0.50.4 d ± 0.10.3 bd ± 0.03.6 cd ± 0.48.8 c ± 0.55.3 ef ± 1.20.2 bc ± 0.00.2 b ± 0.01.8 c ± 0.41.0 e ± 0.16.0 bcd ± 0.3
PA1−23.4 bc ± 0.20.0 def ± 0.131.0 ac ± 0.40.8 bc ± 0.20.5 bd ± 0.14.7 a ± 0.412.7 a ± 0.811.6 cd ± 4.20.2 b ± 0.00.1 d ± 0.03.4 c ± 1.52.0 bc ± 0.46.8 ac ± 0.5
PE1−22.6 a ± 0.31.7 c ± 0.331.5 ab ± 0.81.1 ab ± 0.20.4 bc ± 0.13.4 de ± 0.28.7 c ± 0.78.0 def ± 2.00.2 bc ± 0.00.1 d ± 0.012.1 a ± 5.21.7 cd ± 0.25.0 e ± 0.3
a–f Different mean values in the same column with different superscript letters are significantly different (p < 0.05) according to Tukey’s significant difference test; the Cr was not presented because of values lower than limit of detection.
Table 3. The analysis of variance (ANOVA) of stable isotopes ratios (13C/12C, 15N/14N and 18O/16O) and elemental compositions (mg/kg) of banana peel samples (n = 6).
Table 3. The analysis of variance (ANOVA) of stable isotopes ratios (13C/12C, 15N/14N and 18O/16O) and elemental compositions (mg/kg) of banana peel samples (n = 6).
Farmδ13Cδ15Nδ18OAlBaCrCuFeMnMoNiRbSrZn
CO1−25.0 bd ± 0.83.8 b ± 0.325.7 cd ± 0.324.3 bc ± 2.68.7 b ± 1.10.1 ab ± 0.05.1 ab ± 0.624.8 a ± 9.343.2 bcd ± 7.80.2 bc ± 0.10.4 a ± 0.13.1 d ± 0.722.0 c ± 3.517.8 ab ± 1.7
CR1−24.7 bd ± 0.5−1.1 f ± 0.528.0 b ± 1.023.5 cd ± 7.53.6 cd ± 0.60.1 ab ± 0.14.0 cd ± 0.319.6 ac ± 1.482.2 a ± 15.60.1 d ± 0.00.3 ab ± 0.16.0 d ± 4.720.8 cd ± 3.117.8 ab ± 2.0
CR2−25.0 bd ± 0.30.5 df ± 0.327.9 b ± 1.324.7 bc ± 5.910.6 a ± 1.50.1 ab ± 0.05.8 a ± 0.721.8 ab ± 1.261.7 ab ± 26.90.2 bc ± 0.00.2 b ± 0.024.3 bc ± 7.222.3 c ± 8.515.7 ab ± 1.9
CR3−24.6 bc ± 0.3−0.3 ef ± 0.228.2 b ± 0.929.1 bc ± 13.63.4 cd ± 1.40.0 b ± 0.04.5 bc ± 0.816.2 bcd ± 4.735.0 ce ± 11.30.2 bc ± 0.00.2 c ± 0.116.4 bd ± 3.921.8 c ± 4.416.5 ab ± 1.6
DR1−23.6 a ± 0.52.6 c ± 0.429.8 a ± 0.323.5 cd ± 4.82.4 cd ± 0.30.1 ab ± 0.03.7 cd ± 0.619.9 ac ± 2.518.0 ce ± 6.90.5 a ± 0.10.2 c ± 0.08.8 cd ± 5.341.5 a ± 4.619.0 a ± 2.5
DR2−25.6 d ± 0.58.0 a ± 0.427.2 bc ± 1.35.4 e ± 1.43.0 cd ± 0.10.1 a ± 0.03.4 d ± 0.420.5 ac ± 2.431.7 ce ± 2.50.3 b ± 0.00.2 c ± 0.011.8 cd ± 5.333.7 b ± 1.915.3 ab ± 2.7
EC1−24.4 ab ± 0.51.5 cde ± 0.526.7 bd ± 0.711.1 de ± 2.31.8 d ± 0.40.1 ab ± 0.03.6 cd ± 0.416.5 bcd ± 1.028.0 ce ± 3.60.1 cd ± 0.00.3 bc ± 0.028.7 b ± 5.213.3 de ± 2.415.5 ab ± 1.4
EC2−24.9 bd ± 0.30.5 df ± 0.225.2 d ± 0.521.3 cd ± 4.32.1 cd ± 0.10.1 ab ± 0.04.2 bd ± 0.316.3 bcd ± 1.421.2 de ± 7.00.2 bc ± 0.00.3 bc ± 0.110.8 cd ± 7.610.5 e ± 2.114.3 b ± 1.2
PA1−25.3 cd ± 0.30.8 cde ± 0.628.4 ab ± 0.819.8 cd ± 1.72.6 cd ± 0.60.0 b ± 0.05.6 a ± 0.513.7 cd ± 0.845.8 bc ± 18.70.2 bc ± 0.0n.d.11.9 cd ± 6.023.5 c ± 3.819.0 a ± 1.9
PE1−23.6 a ± 0.21.9 cd ± 0.527.4 b ± 0.837.3 b ± 11.72.9 cd ± 0.20.1 ab ± 0.04.4 bc ± 0.512.2 d ± 1.234.0 ce ± 4.50.2 b ± 0.00.2 c ± 0.147.7 a ± 20.517.0 ce ± 1.818.2 a ± 1.9
a–f Different mean values in the same column with different superscript letters are significantly different (p < 0.05) according to Tukey’s significant difference test; n.d.: none detected.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Erasmus, S.W.; van Ruth, S.M. Preliminary Study on Tracing the Origin and Exploring the Relations between Growing Conditions and Isotopic and Elemental Fingerprints of Organic and Conventional Cavendish Bananas (Musa spp.). Foods 2021, 10, 1021. https://0-doi-org.brum.beds.ac.uk/10.3390/foods10051021

AMA Style

Wang Z, Erasmus SW, van Ruth SM. Preliminary Study on Tracing the Origin and Exploring the Relations between Growing Conditions and Isotopic and Elemental Fingerprints of Organic and Conventional Cavendish Bananas (Musa spp.). Foods. 2021; 10(5):1021. https://0-doi-org.brum.beds.ac.uk/10.3390/foods10051021

Chicago/Turabian Style

Wang, Zhijun, Sara W. Erasmus, and Saskia M. van Ruth. 2021. "Preliminary Study on Tracing the Origin and Exploring the Relations between Growing Conditions and Isotopic and Elemental Fingerprints of Organic and Conventional Cavendish Bananas (Musa spp.)" Foods 10, no. 5: 1021. https://0-doi-org.brum.beds.ac.uk/10.3390/foods10051021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop