Printer Friendly

Post-harvest quality model of pineapple guava fruit according to storage and weather conditions of cultivation/Modelo da qualidade pos-colheita do fruto da goiaba serrana baseado em condicoes de armazenamento e climaticas de cultivo.


Pineapple guava or feijoa (Acca sellowiana (O. Berg) Burret, Myrtaceae family) is a subtropical shrub that is adapted to the weather conditions of the Colombian Andean region and lives at altitudes ranging from 1,800 to 2,700 m.a.s.l. (Parra-Coronado & Fischer, 2013). Covering an estimated total production area of 650 ha, the pineapple guava crops in Colombia include such varieties as Apollo, Triumph, the clones 8-4 and Quimba, which facilitate pollination and high-quality fruit production (Quintero, 2012).

Pineapple guava fruits exhibit considerable differences in their physio-chemical characteristics during harvest, even among fruit from the same cultivar, and these differences directly affect their post-harvest quality (Parra-Coronado & Fischer, 2013). The physio-chemical characteristics of fruit are influenced by weather conditions during cultivation and vary depending on altitude, temperature, radiation, precipitation or available soil moisture and relative humidity (Gruda, 2005; Martinez-Vega et al., 2008; Fischer et al., 2012; Galindo et al., 2014).

Variations in quality during the post-harvest period can be quantified using models based on mathematical functions (Correa Filho et al., 2015). Using a good model makes it possible to establish relationships between the components that would not be evident if the components were considered separately; a good model can be used to facilitate experiments that are not possible in real systems (Hernandez et al., 2009).

This study aimed to develop a model of post-harvest quality in pineapple guava fruit that considers the effects of storage variables (time, temperature and relative humidity (RH)) and prevailing weather conditions during the pre-harvest period to predict the physio-chemical characteristics during storage.


Pineapple guava fruits were harvested when they reached physiological maturity in two production areas located within the Department of Cundinamarca, Colombia. Clone 41 ('Quimba') trees were planted in these areas in 2006. The first site is located in the Tenjo municipality, with geographic coordinates 4[degrees] 51' 23" N and 74[degrees] 6' 33" W. This site has an average altitude of 2.580 m.a.s.l., a mean temperature of 12.5 [degrees]C, a RH between 74 and 86%, a bimodal rainfall pattern, and an average annual precipitation of 765 mm. The second site is located in the San Francisco de Sales municipality, with geographic coordinates 4[degrees] 57' 57" N and 74[degrees] 16' 27" W. This site has an average altitude of 1,800 m.a.s.l., an average temperature of 20.6 [degrees]C, a RH between 63 and 97%, a bimodal rain pattern, and an annual average precipitation of 1,493 mm.

The weather conditions of the study sites from fruit growth and development until harvest (Table 1) were obtained for a two-year period (2012 to 2014). Meteorological data were obtained from weather stations at each sampling site, which provided data on the temperature, RH, precipitation and total radiation on an hourly basis.

Because pineapple guava is a perennial crop, 10 trees per elemental plot and 2 plots per site were used, resulting in a total of 40 trees in this study. The data from one plot for each harvest and farm were used to develop the model, and the data from another plot were used to validate the models. Two harvests were conducted per site, and 300 fruits were collected per plot and per harvest. The fruits that were collected had reached physiological maturity and were free of defects and physical damage. They were collected manually and then transported to the post-harvest laboratory of the Engineering Faculty of the Universidad Nacional de Colombia in Bogota city. The fruits were treated with a 1 mL [L.sup.-1] sodium hypochlorite solution. The fruits were separated according to site of origin and plot and stored at 18 [+ or -] 1 [degrees]C (76 [+ or -] 5% RH, over 11 days, 90 fruits) and 5 [+ or -] 1 [degrees]C (87 [+ or -] 5% RH, over 31 days, 210 fruits).

The measured quality parameters were soluble solids (SS), titratable acidity (TA), skin and pulp firmness, weight loss (WL) and skin color (hue angle, [degrees]h). The technical standard Norma Tecnica Colombiana (NTC) 4624 (ICONTEC, 1999b) was applied when measuring SS with an Eclipse refractometer (Bellingham Stanley, Tunbridge Wells, UK) using a scale from 0 to 32 and a precision of 0.2 [degrees]Brix. To determine the TA, the NTC 4623 Technical Standard (ICONTEC, 1999a) was followed: a sample of juice was taken, which was titrated with a solution of 0.1 N NaOH to neutralize the acid. To determine the fruit skin and pulp firmness, a Brookfield CT3-4500 texture analyzer (Brookfield Engineering, Middleboro, MA, USA) was used, and two readings per fruit were taken with a TA39 probe at [+ or -] 0.5% precision. A gravimetric method (analytical balance "Precisa XT220A", capacity: 220 g, precision: 0.0001 g, Precisa instruments, Switzerland) was used to determine the weight variation during storage; for these measurements, five samples were taken of two fruits each from each storage condition.

The skin color was determined with a Minolta CR-400 colorimeter (Konica Minolta, Ramsey, NJ, USA). The quality parameters were evaluated every 2 days for all fruits under the two storage conditions up to day 11. The evaluation continued every 4 days for fruits stored at 5 [degrees]C until day 31. The statistical design was entirely randomized, with five repetitions per trial.

To analyse the variation in quality attributes of the post-harvest pineapple guava fruits, nonlinear polynomial, exponential and sigmoidal models were evaluated as functions of the storage conditions and place of origin (weather conditions) for each of the measured variables during fruit storage. The model with the best fit was selected based on having the highest coefficient of determination ([R.sup.2]) and lowest standard error. The inputs for each of the variables in the models included information on the storage conditions (number of days, temperature and RH) and place of origin (weather conditions of the crop).

The Excel[R] Solver tool was used for parameter estimation. "Solver is an iterative non-linear procedure, and the first estimation is with the initial values assigned to the parameters in the equation, with those values increasing or decreasing until the lowest value of the sum of the square of the deviations is found" (Parra-Coronado et al., 2015).

The IBM-SPSS statistical package v.20 (Chicago, IL, USA) was used to perform a factor analysis by principal components, and a Varimax rotation was applied to display those variables that had the largest effect on the models. A descriptive statistical analysis was also performed.

Each of the models was evaluated using data that were not used to produce the models. Cross validation between the observed and simulated values was performed to demonstrate that the models provided a good fit to the data. The [R.sup.2] value, index of concordance (d) and root of the mean squared error (RMSE) were calculated to determine the performance of each of the models (Parra-Coronado et al., 2015).


An analysis of variance was performed, and it showed that there were significant differences between locations and harvests (Table 2) for the parameter weights, TA and SS. This result indicates that at the time of harvest, those characteristics were strongly influenced by the local weather conditions and the particular harvest conditions (Fischer et al., 2007; Martinez-Vega et al., 2008). Regarding skin firmness, only the second harvest at Tenjo showed significant differences; thus, the effect of weather conditions during crop growth cannot be clearly defined nor used to explain the behavior of this trait at the time of harvest (Kang et al., 2002). Hue angle and pulp firmness did not exhibit significant differences between location and harvest (Table 2), indicating that weather conditions during fruit growth did not affect the values of these physiochemical characteristics at the time of harvest.

The results of an initial multivariate analysis showed that weather variables and storage conditions (temperature and RH) had an effect on the physio-chemical characteristics of pineapple guava during the post-harvest period. Altitude (H) was the weather variable with the largest effect (-0.959).

The correlation between climatic variables indicated that for greater H, there would be higher radiation (Rad) and lower GDD, P and RH. This trend is consistent with the results reported by Fischer et al. (2012), who reported that there is less precipitation and a lower RH at higher altitude (starting above 1,500 m.a.s.l.). Considering that the weather variable with the greatest effect is H and that correlations exist between H and the other weather variables, a second multivariate analysis was performed (Table 3) considering only the weather variable H, which implicitly correlates the remaining weather variables, to use only those variables with the largest impact on the model.

The results of the principal components analysis when only altitude (H) was included as the weather variable (Table 3), corresponding to the second principal component (PC2), explain 26.0% of the variance and yield the highest commonality values. This result indicates that storage conditions have a strong influence on the development of the models. The effects of H and SS are included in the third principal component (PC3), which explains 15.9% of the variance. The storage period (post-harvest days) and pulp firmness are important variables in the development of the models, and they are included in the first PC, which explains 34.6% of the variance.

The variables corresponding to weight loss, TA, hue angle and skin firmness did not influence the models. However, weight loss and skin firmness are important quality parameters, and their respective communalities had high values (greater than 0.7), indicating the variability explained by the set of principal components or variance shared by other variables in the factorial analysis; therefore, these variables were considered for the models of variation for post-harvest physico-chemical characteristics.

The correlations between variables (Table 3) indicate that as the length of the storage period increases, weight loss and SS increase and TA, hue angle and fruit firmness decrease. These results are consistent with reports by Parra-Coronado & Fischer (2013) regarding pineapple guava fruits. Those studies indicate that weight loss and SS increase and TA and firmness decrease during the post-harvest period. Parra-Coronado et al. (2006) reported similar results for pears.

Weight loss during post-harvest ripening is mainly caused by transpiration and respiration, which causes water and substrate loss (Kays & Paull, 2004). The increase in SS during ripening is caused by organic acids turning into sugars or degrading through respiration (Machado et al., 2015). Loss of firmness during ripening is caused by increased enzymatic activity, which degrades pectins (Parra-Coronado et al., 2006).

The correlations also indicated that at higher storage temperature, weight loss and SS are greater and firmness and TA are lower, which is consistent with the results of Machado et al. (2015), who found that the speed of enzymatic reactions increases exponentially with increasing temperature and is responsible for these physico-chemical changes in orange fruits. Amarante et al. (2008) obtained similar results for pineapple guava, and Parra-Coronado et al. (2006) obtained similar results for pears stored at different temperatures.

The correlations between the variables listed in Table 3 indicate that at higher production zone altitudes, there is a higher SS content and lower values for weight loss and skin and pulp firmness during the post-harvest period. The correlations between weather variables (Table 3) indicate that there is greater radiation (Rad) and lower values for average temperature (GDD), P and RH at higher altitudes; these results are consistent with the behavior of SS reported for the Quimba pineapple guava clone (Martinez-Vega et al., 2008), tomato (Gruda, 2005) and other fruits. The results for firmness are consistent with those found for cucumber (Kang et al., 2002).

Based on the results of the multivariate analysis (Table 3), polynomial, exponential, and sigmoidal equations were evaluated as functions of the storage conditions and location (H) for each of the physico-chemical characteristics using data from both localities and the Excel Solver tool. A regression analysis was performed to determine the best-fit equations, which corresponded to a quadratic polynomial model in all cases (Eq. 1). Table 4 shows the coefficients used in the equations for different attributes of pineapple guava.

The variables are coefficients defined as follows:

At = [[a x T + b x (1/RH) + c x t + d x H].sup.2] + (1)

+ [e x T + f x (1/RH) + g x t + h x H] + i

At--attribute; WL--weight loss of the fruit, %; TA--total titratable acidity, % citric acid; SS--total soluble solids, [degrees]Brix; SkiFir--skin firmness, N; PulFir--pulp firmness, N; hue--color of the epidermis of the fruit; hue angle, [degrees]h; a, b, c, d, e, f, g, h, i--coefficients estimated with Solver; T--storage temperature, [degrees]C; RH--relative humidity of storage, %; t--time in storage, days; and, H--altitude of the production zone, m.a.s.l.

According to the criteria for [R.sup.2], the model of fruit weight loss provides excellent prediction, whereas the model for SS provides adequate prediction. The models for skin and pulp firmness indicate that the model explains more than 50% of their variance, whereas the models for TA and skin color have low predictive ability for the post-harvest period.

Figure 1 shows the cross validation of the models, which involved comparing the predicted and observed values during fruit storage (the data used for cross validation were not used for the development of the model). The cross-validation and regression analysis revealed that the model for weight loss has excellent predictive performance; the model for skin firmness indicate that the model explains more than 50% of their variance, whereas the models for SS and pulp firmness explains less than 50% of the variance, although these can be still considered useful. The values of the concordance index and RMSE also indicate the good fit of these models for the post-harvest period. The models for TA and epidermis color (hue angle) do not adequately estimate these parameters, confirming the results from the multivariate analysis.

The mathematical model for weight loss in the pineapple guava fruit during storage (Eq. 1) could be used as a tool for estimating the maximum storage time in specific conditions so that the fruit exhibits acceptable weight loss, which directly influences the economic yield. Equations for the physico-chemical characteristics of pineapple guava fruits are applicable for storage temperatures between 5 and 18 [degrees]C and for pineapple guava fruits cultivated at altitudes between 1,800 and 2,600 m.a.s.l.

The results presented here have helped to formulate equations that can be used to determine variations in the main physico-chemical characteristics of pineapple guava fruits during the post-harvest period according to storage conditions (number of days, temperature and RH) and prevailing weather conditions during crop growth, which is represented by the altitude of fruit production.


1. The storage conditions (number of days, temperature, and RH) and the prevailing weather conditions during growth, which were represented by the altitude of the fruit production site, strongly influenced the physico-chemical characteristics of pineapple guava fruit during post-harvest ripening.

2. The mathematical models that best predict the changes in post-harvest pineapple guava fruit quality are those for weight loss, SS, and skin and pulp firmness.

3. These models are applicable for storage temperatures between 5 and 18 [degrees]C and for pineapple guava fruits grown at altitudes between 1,800 and 2,600 m.a.s.l. Models for titratable acidity and epidermis color (hue angle) do not adequately estimate these parameters.



Amarante, C. V. T do; Steffens, C. A.; Ducroquet, J. P. H. J.; Sasso, A. Qualidade de goiaba-serrana em resposta a temperatura de armazenamento e ao tratamento com 1-metilciclopropeno. Pesquisa Agropecuaria Brasileira, v.43, p.1683-1689, 2008. https://

Correa Filho, L. C.; Andrade, E. T de; Martinazzo, A. P.; D'Andrea, E. M.; Sousa, F. A. de; Figueira, V. G. Cinetica de secagem, contracao volumetrica e analise da difusao liquida do figo (Ficus carica L.). Revista Brasileira de Engenharia Agricola e Ambiental, v.19, p.797-802, 2015. v19n8p797-802

Fischer, G.; Ebert, G.; Ludders, P. Production, seeds and carbohydrate contents of cape gooseberry (Physalisperuviana L.) fruits grown at two contrasting Colombian altitudes. Journal of Applied Botany and Food Quality, v.81, p.29-35, 2007.

Fischer, G.; Ramirez, F.; Almanza-Merchan, P. J. Induccion floral, floracion y desarrollo del fruto. In: Fischer G. (ed.). Manual para el cultivo de frutales en el tropico. Bogota: Produmedios, 2012. Cap.3. p.120-140.

Galindo, A.; Rodriguez, P.; Collado-Gonzalez, J.; Cruz, Z. N.; Torrecillas, E.; Ondono, S.; Corell, M.; Moriana, A.; Torrecillas, A. Rainfall intensifies fruit peel cracking in water stressed pomegranate trees. Agricultural and Forest Meteorology, v.194, p.29-35, 2014.

Gruda, N. Impact of environmental factors on product quality of greenhouse vegetables for fresh consumption. Critical Reviews in Plant Sciences, v.24, p.227-247, 2005. https://doi. org/10.1080/07352680591008628

Hernandez, F.; Soto, F.; Caballero, A. Modelos de simulacion de cultivos. Caracteristicas y usos. Cultivos Tropicales, v.30, p.73-82, 2009.

ICONTEC--Instituto Colombiano de Normas Tecnicas. Norma tecnica colombiana NTC 4623. Productos de frutas y verduras. Determinacion de la acidez titulable. Bogota: Instituto Colombiano de Normas Tecnicas y Certificacion, 1999a. 6p.

ICONTEC--Instituto Colombiano de Normas Tecnicas. Norma tecnica colombiana NTC 4624. Jugos de frutas y hortalizas. Determinacion del contenido de solidos solubles. Metodo refractometrico. Bogota: Instituto Colombiano de Normas Tecnicas y Certificacion, 1999b. 9p.

Kang, H. M.; Park, K. W.; Saltveit, M. E. Elevated growing temperatures during the day improve the postharvest chilling tolerance of greenhouse-grown cucumber (Cucumis sativus) fruit. Postharvest Biology and Technology, v.24, p.49-57, 2002. https://

Kays, S. J.; Paull, R. E. Postharvest biology. Athens: Exxon Press, 2004. 568p.

Machado, F. L. de C.; Cajazeira, J. P.; Costa, J. M. C. da. Color change and quality response of 'lane late' orange submitted to degreening process. Engenharia Agricola, v.35, p.144-153, 2015. https://doi. org/10.1590/1809-4430-Eng.Agric.v35n1p144-153/2015

Martinez-Vega, R. R.; Fischer, G.; Herrera, A.; Chaves, B.; Quintero, O. C. Caracteristicas fisico-quimicas de frutos de feijoa influenciadas por la posicion en el canopi. Revista Colombiana de Ciencias Horticolas, v.2, p.21-32, 2008.

Parra-Coronado, A.; Fischer, G. Maduracion y comportamiento poscosecha de la feijoa (Acca sellowiana (O. Berg) Burret). Una revision. Revista Colombiana de Ciencias Horticolas, v.7, p.98-11, 2013.

Parra-Coronado, A.; Fischer, G.; Chaves-Cordoba, B. Tiempo termico para estados fenologicos reproductivos de la feijoa (Acca sellowiana (O. Berg) Burret). Acta Biologica Colombiana, v.20, p.167-177, 2015.

Parra-Coronado, A.; Hernandez, J. E. H.; Camacho-Tamayo, J. H. Estudio de algunas propiedades fisicas y fisiologicas precosecha de la pera variedad Triunfo de Viena. Revista Brasileira de Fruticultura, v.28, p.55-59, 2006.

Quintero, O. C. Feijoa (Acca sellowiana Berg). In: Fischer, G. (ed.). Manual para el cultivo de frutales en el tropico. Bogota: Produmedios, 2012. Cap.4, p.443-473.

Alfonso Parra-Coronado (1), Gerhard Fischer (2) & Jesus H. Camacho-Tamayo (1)

(1) Universidad Nacional de Colombia/Facultad de Ingenieria/Departamento de Ingenieria Civil y Agricola. Bogota, Colombia. E-mail: (Corresponding author);

(2) Universidad Nacional de Colombia/Facultad de Ciencias Agrarias/Departamento de Agronomia. Bogota, Colombia. E-mail:

Caption: Figure 1. Cross validation of the physico-chemical characteristics of the pineapple guava fruit in the post-harvest period: (A) weight loss model, (B) TA model, (C) SS model, (D) skin firmness model, (E) pulp firmness model, and (F) color model (hue angle)
Table 1. Weather conditions in the studied zones during the
growth  and development of pineapple guava fruits

Zone               Harvest   Days (1)   GDD (2)   T (3)


Tenjo                          180       1,979    12.3
(2,580 m.a.s.l.)      2        180       1,966    12.3
San Francisco         1        155       2,728    18.5
(1,800 m.a.s.l.)      2        155       2,627    18.0

Zone               Harvest   RH (4)   P (5)      Rad (6)
                              (%)     (mm)    (W [m.sup.-2])

Tenjo                         76.4     190        12,303
(2,580 m.a.s.l.)      2       84.3     417         9861
San Francisco         1       86.1     573         7814
(1,800 m.a.s.l.)      2       95.1    1,400       10,021

(1) Days--Calendar days elapsing from anthesis to harvest;
(2) GDD--Growing degree days accumulated from anthesis to harvest;
(3) T--Mean temperature during the study period;
(4) RH--Mean relative humidity during the study period;
(5) P--Accumulated precipitation from anthesis to harvest;
(6) Rad--Accumulated radiation from anthesis to harvest

Table 2. Mean values (1) of the physico-chemical characteristics
of pineapple guava fruit at the time of harvest

Parameter                         Location--Harvest

                         Tenjo -1                   Tenjo -2

Fresh            38.23 [+ or -] 4.23 (bc)    98.93 [+ or -] 12.62 (a)
  weight (g)
SS ([degrees]    13.35 [+ or -] 0.66 (a)     11.73 [+ or -] 0.91 (b)
TA (% citric      1.91 [+ or -] 0.12 (a)     11.68 [+ or -] 0.09 (b)
Hue angle        124.72 [+ or -] 0.75 (a)   123.63 [+ or -] 1.35 (a)
Skin             14.82 [+ or -] 3.51 (a)     10.21 [+ or -] 2.64 (b)
  firmness (N)
Pulp             106.90 [+ or -] 2.13 (a)    16.14 [+ or -] 4.03 (a)
  firmness (N)

Parameter                        Location--Harvest

                     San Francisco-1            San Francisco-2

Fresh             30.53 [+ or -] 4.67 (c)    45.73 [+ or -] 6.83 (b)
  weight (g)
SS ([degrees]     11.19 [+ or -] 0.81 (b)   11.59 [+ or -] 0.66 (b)
TA (% citric     11.58 [+ or -] 0.07 (b)    11.93 [+ or -] 0.06 (a)
Hue angle        121.63 [+ or -] 2.53 (a)   124.16 [+ or -] 1.36 (a)
Skin              16.20 [+ or -] 1.30 (a)    14.18 [+ or -] 1.87 (a)
  firmness (N)
Pulp             15.47 [+ or -] 2.08 (a)    16.12 [+ or -] 1.89 (a)
  firmness (N)

(1) Mean [+ or -] sd; Means followed by different letters for the
same parameter indicate significant differences according to Tukey's
test (p [less than or equal to] 0.05)

Table 3. Multivariate analysis of all of the data from Tenjo
and San Francisco considering altitude (H) as the only weather

Variable            PC1      PC2      PC3      Communality

Post-harvest days   -0.730    0.481    0.063      0.769
H                    0.076    0.047   -0.936      0.884
Storage              0.010   -0.975   -0.039      0.953
Storage RH          -0.010    0.975    0.039      0.953
Weight loss         -0.604   -0.605    0.261      0.800
SS                  -0.264   -0.101   -0.848      0.800
TA                   0.479    0.419    0.069      0.410
Hue Angle            0.696   -0.056    0.295      0.575
Skin firmness        0.573    0.024    0.654      0.756
Pulp firmness        0.772    0.197    0.349      0.757
Eigenvalue           3.457    2.603    1.595
Variance (%)        34.575   26.033   15.945
Cumulative          34.575   60.608   76.553
  variance (%)

* PC1, PC2, PC3: principal components 1, 2 and 3. Values > 0.7
(absolute value) appear in bold; H--Altitude; TA--Titratable acidity

Table 4. Coefficients of the equations for different attributes
of pineapple guava


Coefficients      WL           TA           SS

a               -2.6131     -2.6973      -0.7709
b              1492.2892   2004.4856     302.4858
c               -0.0242     0.000605     -0.0116
d              0.000069    1.8909E-06   -0.000011
e              -55.3924     -30.3962     -7.5536
f              -0.00064     -0.00064     -0.00064
g               0.3782      -0.01418      0.0289
h               -0.0021     -0.00007      0.0041
i              286.4847     86.5955      42.1316
[R.sup.2](*)     0.93         0,32         0.68
St.E(*)          1.29         0.12         0.99


Coefficients    SkiFir     PulFir        hue

a              -1.7104     -2.078      -2.9739
b              711.0051   997.5706    2435.8724
c               0.0176     0.0089      0.0027
d              0.000036   -0.000027   -2.61E-06
e              -34.6905   -45.2278     -21.850
f              -0.00064   -0.00064    -0.00064
g              -0.1591     -0.2026    -0.19716
h              -0.0065     -0.0023     -0.0013
i              214.3106   252.1828    94.41056
[R.sup.2](*)     0.56       0.54        0.24
St.E(*)          2.33       1.48        1.13

WL-weight loss of the fruit, %; TA--total titratable acidity, %
citric acid; SS--total soluble solids, [degrees]Brix; SkiFir--skin
firmness, N; PulFir--pulp firmness, N; hue--color of the epidermis
of the fruit; hue angle, [degrees]h. [R.sup.2](*)--Determination
coefficient; St.E(*).--Standard error
COPYRIGHT 2017 ATECEL--Associacao Tecnico Cientifica Ernesto Luiz de Oliveira Junior
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Parra-Coronado, Alfonso; Fischer, Gerhard; Camacho-Tamayo, Jesus H.
Publication:Revista Brasileira de Engenharia Agricola e Ambiental
Article Type:Ensayo
Date:Sep 1, 2017
Previous Article:Potassium fertilization in the cultivation of colored cotton irrigated with saline water/Adubacao potassica no cultivo do algodoeiro de fibra...
Next Article:Physico-chemical composition of native fruits of the Ceara coast at different development stages/Composicao fisico-quimica de frutas nativas do...

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters