Abstract
Coffee bean chemical compositions has been extensively studied. However, there is a small amount of research on other parts of the coffee plant, including leaves. Fourier transform infrared (FTIR) spectral profiles of Coffea arabica L. cv. IAPAR 59 leaf extracts from a simplex-centroid design were studied by principal component analysis (PCA) to evaluate the effect of solvent extractor on its metabolites. PCA indicated that the extractor solvents containing ethanol were the most suitable for this study. FTIR spectra in conjunction with orthogonal signal correction and partial least squares-discrimination analysis (OSC-PLS-DA) were used to classify and discriminate the leaves of irrigated and non-irrigated plants by bands related to carbohydrates, amino acids and lipids. Leaves receiving different intensities of solar radiation were also discriminated by bands corresponding to caffeine, carbohydrates and lipids. FTIR spectral profile analyzed with chemometric tools showed to be a useful, powerful and simple procedure to discriminate coffee leaves collected from different microclimate conditions.
Keywords:
Coffea arabica leaves; simplex-centroid design; FTIR spectroscopy; chemometric analysis; environmental conditions
Introduction
The coffee plant is a woody, perennial, evergreen and dicotyledon species that belongs to the Rubiaceae family.11 Davis, A. P.; Govaerts, R.; Bridson, D. M.; Stoffelen, P.; Bot. J. Linn. Soc. 2006, 152, 465.,22 de Lima, R. B.; Santos, T. B.; Vieira, L. G. E.; Ferrarese, M. L. L.; Ferrarese-Filho, O.; Donatti, L.; Boeger, M. R. T.; Petkowicz, C. L. O.; Carbohydr. Polym. 2013, 93, 135. Among several species of the genus Coffea, two are considered economically important, i.e. Coffea canephora Pierre (Robusta coffee) and Coffea arabica L. (Arabica coffee).33 Komes, D.; Vojvodić, A. In Processing and Impact on Antioxidants in Beverages; Preedy, V. R., ed.; Academic Press, Elsevier Inc.: Waltham, USA, 2014, ch. 9. The latter species represents 70% of world coffee production22 de Lima, R. B.; Santos, T. B.; Vieira, L. G. E.; Ferrarese, M. L. L.; Ferrarese-Filho, O.; Donatti, L.; Boeger, M. R. T.; Petkowicz, C. L. O.; Carbohydr. Polym. 2013, 93, 135. due to finer flavor and aroma, and therefore is more consumed than Robusta.44 Barbin, D. F.; Felicio, A. L. S. M.; Sun, D. W.; Nixdorf, S. L.; Hirooka, E. Y.; Food Res. Int. 2014, 61, 23.
Coffee is one of the most consumed beverages in the world55 Marquetti, I.; Link, J. V.; Lemes, A. L. G.; Scholz, M. B. S.; Valderrama, P.; Bona, E.; Comput. Electron. Agric. 2016, 121, 313.,66 Wang, N.; Lim, L. T.; J. Agric. Food Chem. 2012, 60, 5446. and an important raw material of the international trade.44 Barbin, D. F.; Felicio, A. L. S. M.; Sun, D. W.; Nixdorf, S. L.; Hirooka, E. Y.; Food Res. Int.
2014, 61, 23. Thus, due to its importance, the chemical composition of coffee beans has been extensively studied. However, there is relatively little research carried out on metabolites in other coffee plant parts, including the leaves.77 Campa, C.; Mondolot, L.; Rakotondravao, A.; Bidel, L. P. R.; Gargadennec, A.; Couturon, E.; La Fisca, P.; Rakotomalala, J. J.; Jay-Allemand, C.; Davis, A. P.; Ann. Bot.
2012, 110, 595.
8 Clifford, M. N.; Kirkpatrick, J.; Kuhnert, N.; Roozendaal, H.; Salgado, P. R.; Food Chem.
2008, 106, 379.
9 Delaroza, F.; Rakocevic, M.; Malta, G. B.; Bruns, R. E.; Scarminio, I. S.; J. Braz. Chem. Soc.
2014, 25, 1929.
10 Domingues Júnior, A. P.; Shimizu, M. M.; Moura, J. C. M. S.; Catharino, R. R.; Ramos, R. A.; Ribeiro, R. V.; Mazzafera, P.; Photochem. Photobiol.
2012, 88, 928.-1111 Mondolot, L.; la Fisca, P.; Buatois, B.; Talansier, E.; de Kochko, A.; Campa, C.; Ann. Bot.
2006, 98, 33. Phenolic compounds found in the coffee leaves77 Campa, C.; Mondolot, L.; Rakotondravao, A.; Bidel, L. P. R.; Gargadennec, A.; Couturon, E.; La Fisca, P.; Rakotomalala, J. J.; Jay-Allemand, C.; Davis, A. P.; Ann. Bot.
2012, 110, 595. have been shown to be potentially beneficial for health, although the influence of consumption of these compounds on the human body requires further research. Other metabolite groups are found in coffee leaves, such as alkaloids,1212 Zheng, X. Q.; Ashihara, H.; Plant Sci. 2004, 166, 807. several carbohydrates1313 Brum, C. N. F.; Melo, E. F.; Barquero, L. O. B.; Alves, J. D.; Chalfun-Júnior, A.; Coffee Sci.
2013, 8, 140.,1414 Lima, R. B.; Santos, T. B.; Vieira, L. G. E.; Ferrarese, M. L. L.; Farrarese-Filho, O.; Donatti, L.; Boeger, M. R. T.; Petkowicz, C. L. O.; Carbohydr. Polym.
2014, 112, 686. and lipids99 Delaroza, F.; Rakocevic, M.; Malta, G. B.; Bruns, R. E.; Scarminio, I. S.; J. Braz. Chem. Soc.
2014, 25, 1929. among others.
Plant metabolites are susceptible to environmental changes and their responses can be evaluated by metabolomics analysis.1515 Zhao, J.; Hu, C.; Zeng, J.; Zhao, Y.; Zhang, J.; Chang, Y.; Li, L.; Zhao, C.; Lu, X.; Xu, G.; Metabolomics
2014, 10, 805. In periods of drought or excessive sun exposure, environmental stress can affect physiological processes in plants owing to their necessity to adapt to the new situation. In excess, these conditions can affect metabolic activities such as photosynthesis and growth rate modifying carbohydrate levels and protein synthesis.1616 Chaves, A. R. M.; Ten-Caten, A.; Pinheiro, H. A.; Ribeiro, A.; DaMatta, F. M.; Trees
2008, 22, 351.
17 Da Matta, F. M.; Ramalho, J. D. C.; Braz. J. Plant Physiol. 2006, 18, 55.
18 Jones, M. M.; Turner, N. C.; Plant Physiol. 1978, 61, 122.
19 Hsiao, T. C.; Acevedo, E.; Fereres, E.; Henderson, D. W.; Philos. Trans. R. Soc., B
1976, 273, 479.-2020 Praxedes, S. C.; DaMatta, F. M.; Loureiro, M. E.; Ferrão, M. A. G.; Cordeiro, A. T.; Environ. Exp. Bot.
2006, 56, 263. In coffee plants, these changes may influence the metabolic quantities and consequently, the coffee bean quality.99 Delaroza, F.; Rakocevic, M.; Malta, G. B.; Bruns, R. E.; Scarminio, I. S.; J. Braz. Chem. Soc.
2014, 25, 1929.,2121 Gobbo-Neto, L.; Lopes, N. P.; Quim. Nova 2007, 30, 374. Information about chemical modifications in coffee leaves subjected to environmental stress can be obtained by metabolomic analysis. It consists of sampling, sample preparation, use of instrumental analytical techniques, processing and data interpretation.2222 Dunn, W. B.; Ellis, D. I.; TrAC, Trends Anal. Chem. 2005, 24, 285. A more global metabolomics approach is called metabolic fingerprinting that can be used for purposes such as the quality control of medicinal plants as well as their characterization and classification.2323 Ivanišević, J.; Thomas, O. P.; Lejeusne, C.; Chevaldonné, P.; Pérez, T.; Metabolomics
2011, 7, 289.
Metabolic fingerprinting of an herbal sample is a characteristic profile of chromatographic or spectroscopic origin, which characterizes its composition.2424 Viaene, J.; Goodarzi, M.; Dejaegher, B.; Tistaert, C.; Hoang Le Tuan, A.; Nguyen Hoia, N.; Chau Van, M.; Quetin-Leclercq, J.; Vander Heyden, Y.; Anal. Chim. Acta 2015, 877, 41. Spectroscopy is an increasingly growing technique due to its rapidity, simplicity, and safety, as well as its ability to measure multiple attributes simultaneously without arduous sample preparation.44 Barbin, D. F.; Felicio, A. L. S. M.; Sun, D. W.; Nixdorf, S. L.; Hirooka, E. Y.; Food Res. Int. 2014, 61, 23. Fourier transform infrared (FTIR) spectroscopy is a valuable fingerprinting tool, owing to its ability to simultaneously measure proteins, lipids, polysaccharides and other metabolic components.2525 Ellis, D. I.; Harrigan, G. G.; Goodacre, R. In Metabolic Profiling: Its Role in Biomarker Discovery and Gene Function Analysis; Harrigan, G. G.; Goodacre, R., eds.; Kluwer Academic Publishers: London, UK, 2003, ch. 7. Furthermore, recent publications show that FTIR spectroscopy has been applied to signature plant material due to its ability to identify environmental effects within plant material, even in living plants.2626 Butler, H. J.; McAinsh, M. R.; Adams, S.; Martin, F. L.; Anal. Methods 2015, 7, 4059.,2727 Ord, J.; Butler, H. J.; McAinsh, M. R.; Martim, F. L.; Anal. Methods 2016, 141, 2896. This spectra methodology combined with chemometric methods can provide specific information about different parameters simultaneously in a direct, reliable and rapid way.2828 Gallardo-Velázquez, T.; Osorio-Revilla, G.; Zuñiga-de Loa, M.; Rivera-Espinoza, Y.; Food Res. Int. 2009, 42, 313.
Chemometrics reduces the large amount of data produced by automated instruments and enables us to extract information and identify interesting patterns or features.2929 Hendriks, M. M. W. B.; Cruz-Juarez, L.; de Bont, D.; Hall, R. D.; Anal. Chim. Acta 2005, 545, 53. Because of the complexity of FTIR spectral data, the chemometric analysis is used to reduce the dimensionality of spectral data to aid the extraction of useful information.3030 Wang, N.; Fu, Y.; Lim, L. T.; J. Agric. Food Chem. 2011, 59, 3220. Principal component analysis (PCA)3131 Wold, S.; Esbensen, K.; Geladi, P.; Chemom. Intell. Lab. Syst. 1987, 2, 37. is a multivariate method which can be used for exploratory data analysis and allows the transformation and visualization of complex data sets providing a new and hopefully simpler perspective, from which more relevant environmental information can be easily perceived.3232 García-Reiriz, A. G.; Olivieri, A. C.; Teixidó, E.; Ginebreda, A.; Tauler, R.; Environ. Sci.: Processes Impacts 2014, 16, 124. Another widely used chemometric method, partial least squares discriminant analysis (PLS-DA), is a linear classification method that combines the properties of partial least squares regression with the discrimination power of a classification technique.3333 Ballabio, D.; Consonni, V.; Anal. Methods 2013, 5, 3790.
The hypothesis of this study was that the water and light availability will leave specific metabolic fingerprints in Arabica coffee leaves. Spectra in the infrared region (FTIR) were employed to analyze C. arabica leaf extracts, prepared using a simplex-centroid design with ethanol, dichloromethane and hexane, in combination with chemometric tools. The first step was performed by studying the effect of these solvents on metabolic extraction from coffee leaves through FTIR spectral profiles and PCA. Then, the FTIR spectra of these extracts were investigated using orthogonal signal correction and partial least squares-discrimination analysis (OSC-PLS-DA) to study metabolites in leaves formed on plants cultivated with and without irrigation and collected from plant layers exposed to different levels of sun irradiation.
Experimental
Plant material
C. Arabica plants, cultivar IAPAR 59, were cultivated in the experimental area of the Agronomic Institute of Paraná, Londrina (23º 18' S, 51º 09' W), Paraná, Brazil. The coffee trees were planted in 2010 in high density and an arrangement at 2.5 × 0.5 m with idealized 1.25 m22 de Lima, R. B.; Santos, T. B.; Vieira, L. G. E.; Ferrarese, M. L. L.; Ferrarese-Filho, O.; Donatti, L.; Boeger, M. R. T.; Petkowicz, C. L. O.; Carbohydr. Polym. 2013, 93, 135. available for the development of each plant. This cultivar is shown to be well adapted to high density planting designs.3434 Matsunaga, F. T.; Tosti, J. B.; Androcioli-Filho, A.; Brancher, J. D.; Costes, E.; Rakocevic, M.; SpringerPlus 2016, 5, 2015. Leaves were collected in October 2012 from non-irrigated and drop by drop irrigated plants. In order to check plant architectural influence and light irradiation effects, the harvest was stratified into an inferior layer (40 cm, self-shaded) and a superior one (> 80 cm, high light exposure).
The leaf samples were dried in circulating air at ambient temperature. The leaves were distributed into trays and every 24 h turned over to allow complete drying. This drying process lasted 15 to 20 days depending on the number of leaves in each tray. Samples were crushed in a domestic blender, passed through a plastic sieve, packed in plastic bags subjected to vacuum and then stored in a freezer/cooler at -18 ºC.
Extract preparation
A simplex-centroid design with three components was used to obtain the leaf extracts. The ethanol (e), dichloromethane (d) and hexane (h) pure solvents, their three (1:1) binary mixtures and their (1:1:1) ternary mixture and three axial ternary mixtures, (4:1:1), (1:4:1) and (1:1:4) were investigated (Table 1). These mixtures were prepared in random order including six replicates at the (1:1:1) center mixture. Solvent selection was based on diversity considering Snyder's solvent selectivity triangle.3535 Snyder, L. R.; Carr, P. W.; Rutan, S. C.; J. Chromatogr. A 1993, 656, 537.
The extractions were performed with 2.0 g of dried and crushed C. arabica leaves and 60 mL of extracting solvent for 24 h, which were subsequently filtered to separate the solution from the coffee leaves. This procedure was repeated three more times, for a total of four repetitions. Thus, the total volume of solvent mixture added to leaves was 240 mL for each point of simplex-centroid design. These solvents were evaporated in a rotary evaporator, kept under forced ventilation, and later the extracts were lyophilized. For leaf extraction, all organic solvents were of analytical grade with dichloromethane and hexane obtained from Anidrol (São Paulo, Brazil) and ethanol from Êxodo Científica (Hortolândia, Brazil).
FTIR fingerprint measurements
Spectra in the infrared region (4000 to 945 cm-1) of crude extracts were obtained with a Thermo Scientific Nicolet iS10 FT-IR spectrometer, using the ATR (attenuated total reflectance) accessory with a Ge window. Spectra were obtained with 32 scans and 4 cm-1 resolution.
FTIR data treatment, model and validation
The analysis of the spectra was performed with PLS toolbox 5.8.1 (Eigenvector Research Inc., Wenatchee, WA, USA) from Matlab R2007a (Mathworks Inc. Natick, MA, USA).
First, PCA was used for an initial exploration of the solvent effects on the C. arabica leaf sample extraction procedure by FTIR data set. PCA is an unsupervised multivariate statistical method which calculates new orthogonal axes (principal components) from the linear combinations of the original data variables, in the direction of explaining maximum data variance.3131 Wold, S.; Esbensen, K.; Geladi, P.; Chemom. Intell. Lab. Syst. 1987, 2, 37. The FTIR spectra were preprocessed using the Savitzky-Golay second-derivative and were mean-centered. Then, the spectra were classified by PLS-DA.3636 Barker, M.; Rayens, W. J.; J. Chemom. 2003, 17, 166. This chemometric tool is considered a supervised method and is based on partial least squares regression (PLSR).3737 Geladi, P.; Kowalski, B. R.; Anal. Chim. Acta 1986, 185, 1. In a very efficient way PLS-DA estimates the best linear combinations of the independent original X-values (latent variables, LV), which optimally correlate with the observed changes of the dependent variable, y (set of binary variables describing the categories of X).3838 Farrés, M.; Piña, B.; Tauler, R.; Metabolomics 2015, 11, 210. The computed score plots give an idea of similarity among the samples whereas the loading plots show the importance of each variable in the modeling.3939 Carlos, C.; Maretto, D. A.; Poppi, R. J.; Sato, M. I. Z.; Ottoboni, L. M. M.; Microchem. J. 2011, 99, 15.
To obtain the PLS-DA model, the second derivative was applied to the spectra (as in the PCA procedure) and the OSC (orthogonal signal correction) algorithm4040 Wold, S.; Antti, H.; Lindgren, F.; Öhman, J.; Chemom. Intell. Lab. Syst. 1998, 44, 175. was utilized to eliminate unnecessary information from the model. The benefit of using the OSC-filter in PLS-DA, i.e. OSC-PLS-DA, is its ability to separate predictive from non-predictive (orthogonal) variations.4141 Bylesjö, M.; Rantalainen, M.; Cloarec, O.; Nicholson, J. K.; Holmes E.; Trygg, J.; J. Chemom. 2006, 20, 341. In this procedure, the X matrix is corrected by the subtraction of the variation which is orthogonal to y and assigns a correspondence class to each sample.3939 Carlos, C.; Maretto, D. A.; Poppi, R. J.; Sato, M. I. Z.; Ottoboni, L. M. M.; Microchem. J. 2011, 99, 15. These data were also mean-centered.
The OSC-PLS-DA model was built considering the X matrix (FTIR spectra), while the y vector was associated with two different classes. The determination of the correct number of latent variables is fundamental to build PLS-DA models. This is commonly performed by cross-validation applied to calibration of samples that are then verified with a validation set of samples.4242 de Almeida, M. R.; Correa, D. N.; Rocha, W. F. C.; Scafi, F. J. O.; Poppi, R. J.; Microchem. J. 2013, 109, 170. Thus, the leaf extracts were divided into validation and calibration subsets with the Kennard-Stone algorithm.4343 Kennard, R. W.; Stone, L. A.; Technometrics 1969, 11, 137. The validation set was used to evaluate the rate of correct assignment of the OSC-PLS-DA model. Furthermore, a threshold value is calculated from the predicted values, and used for sample assignment to the modeled classes.4444 Alves, F. C. G. B. S.; Valderrama, P.; Anal. Methods 2015, 7, 9702. The threshold value for class separation is based on Bayes' theorem.4545 de Almeida, M. R.; Fidelis, C. H. V.; Barata, L. E. S.; Poppi, R. J.; Talanta 2013, 117, 305.,4646 Pérez, N. F.; Ferré, J.; Boqué, R.; Chemom. Intell. Lab. Syst. 2009, 95, 122.
Results and Discussion
Exploratory PCA was performed on the spectra of the simplex-centroid extracts. Spectral data were arranged in a matrix of 1525 rows and 60 columns. Each row represented a wavenumber variable (4000 to 945 cm-1) and each column a sample extract (15 simplex design mixtures for four environmental conditions). The first three components explained 73.38% of the total data variance. The most informative score plot involved the first and third principal components. PC1 differentiated extracts obtained with and without ethanol and is shown in Figure 1a. Extracts prepared with hexane and dichloromethane as pure solvents and its binary mixture are on the more positive side of PC1 and extracts prepared containing ethanol are on the negative side of this principal component. The PC1 loading plot is presented in Figure 1b. It was analyzed in conjunction with the average spectra of the ethanol and non-ethanol extracts given in Figure 2.
(a) PC1 vs. PC3 scores plot, according to solvent mixture compositions and (b) PC1 loading plot of FTIR spectra from C. arabica cv. IAPAR 59 crude extracts of leaves from (●) irrigated and (■) non-irrigated plants.
Average FTIR spectra of crude leaf extracts obtained from the simplex-centroid mixture design. C. arabica cv. IAPAR 59 leaves were collected from irrigated and non-irrigated plants and extracts were prepared (a) with ethanol and (b) without ethanol in the extractor solvent compositions.
All extracts (containing ethanol or not) showed bands in the 3000-2800 cm-1 region. However, the bands at 2960 and 1376 cm-1 contribute to separate extracts containing ethanol, with loading values < -0.1. The band at 2960 cm-1 is characteristic of C-H asymmetric stretching of the CH3. The original spectral plot shows this band overlapped with two other bands at 2918 and 2850 cm-1 referring to C-H symmetric stretching of the CH3 and CH2, respectively, indicating the extraction of lipids/proteins.3030 Wang, N.; Fu, Y.; Lim, L. T.; J. Agric. Food Chem. 2011, 59, 3220.,4747 Movasaghi, Z.; Rehman, S.; Rehman, I. U.; Appl. Spectrosc. Rev.
2008, 43, 134.,4848 Zavoi, S.; Fetea, F.; Ranga, F.; Pop, R. M.; Baciu, A.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca
2011, 39, 82. A band at 1376 cm-1 can be associated to the ionone ring of carotene or due to the C-H symmetrical bending of the methyl of the aliphatic bonds.4949 Kushwaha, K.; Saxena, J.; Tripathi, B. K.; Agarwal, M. K.; J. Biosci. Biotechnol. 2014,3, 253. Other regions with important discrimination between extracts with and without ethanol were around 1676 and 1092 cm-1, due to large negative values and around 1697, 1073 and 1027 cm-1, with large positive loading values (Figure 1b). The bands between 1780 to 1600 cm-1 correspond to the carbonyl (C=O) stretching vibration of organic compounds,66 Wang, N.; Lim, L. T.; J. Agric. Food Chem. 2012, 60, 5446. which can be attributed to diverse compounds such as carotenoids, chlorogenic acids, alkaloids, polysaccharides, hemicellulose, among others, depending on the functional group.66 Wang, N.; Lim, L. T.; J. Agric. Food Chem. 2012, 60, 5446.,99 Delaroza, F.; Rakocevic, M.; Malta, G. B.; Bruns, R. E.; Scarminio, I. S.; J. Braz. Chem. Soc.
2014, 25, 1929.,5050 Craig, A. P.; Franca, A. S.; Oliveira, L. S.; Food Chem. 2012, 132, 1368.
51 Pop, R. M.; Buzoianu, A. D.; Rati, I. V.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca
2014, 42, 545.
52 Ribeiro, J. S.; Salva, T. J.; Ferreira, M. M. C.; J. Food Qual.
2010, 33, 212.-5353 Terrile, A. E.; Marcheafave, G. G.; Oliveira, G. S.; Rakocevic, M.; Bruns, R. E.; Scarminio, I. S.; J. Braz. Chem. Soc. 2016, 7, 1254. The bands between 1300 to 1000 cm-1 are assigned to C-O bond of esters, ethers, alcohols, among others.4848 Zavoi, S.; Fetea, F.; Ranga, F.; Pop, R. M.; Baciu, A.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca
2011, 39, 82.,5151 Pop, R. M.; Buzoianu, A. D.; Rati, I. V.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca
2014, 42, 545.,5454 Kozarski, M.; Klaus, A.; Vunduk, J.; Zizak, Z.; Niksic, M.; Jakovljevic, D.; Vrvic, M. M.; Van Griensven, L. J. L. D.; Food Funct. 2015, 6, 1875.,5555 Schulz, H.; Baranska, M.; Vib. Spectrosc. 2006, 43, 13.
Crude extracts obtained by solvents without ethanol possibly contain lipid compounds, while extracts obtained with solvents containing ethanol extract these compounds as well as other metabolites containing more functional groups. This last group of extracts also contained the most similar spectral fingerprints. Thus, it was decided to apply PLS-DA only to this data. PLS-DA is recommended for attempting to find the best correlation between the X and y when the variability within groups is greater than the variability among groups.3636 Barker, M.; Rayens, W. J.; J. Chemom. 2003, 17, 166. First, the PLS-DA was applied to discriminate the coffee leaves subjected to two water availability conditions, for which the variable y was categorized as class 0 for leaves from irrigated and 1 for those ones from non-irrigated plants. The second derivative was applied to the spectra to perform baseline adjustments (second order polynomial with a window width of 15 points) and the OSC algorithm was used to eliminate unnecessary information in the classification model.
The spectral data set was divided, using the Kennard-Stone algorithm, into two subsets: calibration and validation. The calibration subset contained 19 spectra from crude extracts of irrigated leaves and 15 from non-irrigated ones, which represent 70% of the extracts obtained with ethanol. For the validation subset, five spectra from crude extracts of irrigated leaves and 9 from non-irrigated ones were used, representing 30% of the extracts.
To support the model dimensionality selection, venetian blinds cross-validation was employed, based on the lowest value of the root mean square error of cross-validation (RMSECV). Thus, three latent variables were needed to develop the classification model representing 61.17% of explained variance in X and 75.90% of variance in y, with RMSECV equal to 0.35 and mean error of 0.24 for the test subset root mean square error of calibration (RMSEC).
The score plot for the best OSC-PLS-DA model is presented in Figure 3a, showing the separation between samples submitted to two different water supply conditions. In this case, the extracts of irrigated plants were discriminated, being mostly in the negative regions of LV1 and LV2, while most leaf sample points of non-irrigated plants fall in the positive region. The loading plot, Figure 3b, shows the wavelengths which most contribute to the discrimination of two water availability conditions. The bands at 1665, 1550, and 1073 cm-1 contributed to classify extracts of leaves from non-irrigated plants, because they had positive loadings for LV1 and LV2. For leaves from irrigated plants, the bands 2918; 2850; 1665; 1091; 1076; and 1023 cm-1 were most important, because they had negative loadings for both latent variables.
(a) LV1 vs. LV2 scores plot from OSC-PLS-DA model; (b) loading plots of LV1 (-) and LV2 (...) built from FTIR spectra for crude leaf extracts with second-derivative preprocess and OSC filter. Extracts of leaves from (●) irrigated and (■) non-irrigated plants of C.arabica cv. IAPAR 59 were analyzed.
The irrigated C. arabica leaves were discriminated by bands at 2918 and 2850 cm-1, which may correspond to lipids such as carotenoid.5656 Lóránd, T.; Deli, J.; Molnár, P.; Tóth, G.; Helv. Chim. Acta 2002, 85, 1691. Given this, these two bands were already analyzed by PCA and were associated with carotenoids due to the presence of the band around 1376 cm-1 that was related to ionone ring of carotene.4949 Kushwaha, K.; Saxena, J.; Tripathi, B. K.; Agarwal, M. K.; J. Biosci. Biotechnol. 2014,3, 253. All photosynthetic organisms contain carotenoids, which are essential for photoprotection, usually function also as accessory pigments, and in many cases serve as key regulatory molecules.5757 Sandman, G.; Arch. Biochem. Biophys. 2009, 483, 169. Carotenoids are essential structural components of the photosynthetic antenna and reaction center complexes. Therefore, carotenoids serve as accessory pigments by harvesting radiant light in a region of the spectrum not covered by the chlorophylls.5858 Bartley, G. E.; Scolnik, P. A.; Plant Cell 1995, 7, 1027. They exert their photoprotective action through the rapid uptake of energy absorbed by chlorophylls, excited by light due to excitation transfer or photochemistry.5959 Wollman, F. A.; EMBO J. 2001, 20, 3623. If the excited state of chlorophyll is not rapidly transferred, it can react with molecular oxygen to form radicals that damage many cellular components, especially the photosynthetic membranes. When energy absorbed by the pigments is large and cannot be stored, they can be easily damaged. The excited state of carotenoids does not have enough energy to form radicals, decaying back to the state base while losing its energy as heat.6060 Szabò, I.; Bergantino, E.; Giacometti, G. M.; EMBO Rep. 2005, 6, 629.
Under water stress, amino acid levels may elevate due to a possible activity increase of protease enzymes to break proteins reservations and to osmotically adjust at the stressful environment.6161 Lechinoski, A.; de Freitas, J. M. N.; Castro, D. S.; Lobato, A. K. S.; de Oliveira Neto, C. F.; Cunha, R. L. M.; Rev. Bras. Bioc. 2007, 5, 927. Therefore, the band around 1550 cm-1 might be assigned as the amino group of amino acids.6262 Pavia, D. L.; Lampman, G. M.; Kriz, G. S.; Introduction to Spectroscopy: A Guide of Students of Organic Chemistry, 3rd ed.; Thomson Learning, Books/Cole: Seneca, Mexico, 2011. Under irrigation or not, extracts from coffee leaves presented a band around 1665 cm-1 as important variable. This band is characteristic of the presence of the carbonyl group, for which this region may indicate the presence of several compounds such as proteins,5555 Schulz, H.; Baranska, M.; Vib. Spectrosc. 2006, 43, 13. caffeine,66 Wang, N.; Lim, L. T.; J. Agric. Food Chem. 2012, 60, 5446.,3030 Wang, N.; Fu, Y.; Lim, L. T.; J. Agric. Food Chem. 2011, 59, 3220. ketone between two aromatic rings6262 Pavia, D. L.; Lampman, G. M.; Kriz, G. S.; Introduction to Spectroscopy: A Guide of Students of Organic Chemistry, 3rd ed.; Thomson Learning, Books/Cole: Seneca, Mexico, 2011. as in xanthone structures found in the form of mangiferin detected in coffee leaves,77 Campa, C.; Mondolot, L.; Rakotondravao, A.; Bidel, L. P. R.; Gargadennec, A.; Couturon, E.; La Fisca, P.; Rakotomalala, J. J.; Jay-Allemand, C.; Davis, A. P.; Ann. Bot.
2012, 110, 595. among others. Bands indicated by loadings between 1100 and 1000 cm-1 are possibly C-O stretches of carbohydrates.4848 Zavoi, S.; Fetea, F.; Ranga, F.; Pop, R. M.; Baciu, A.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca
2011, 39, 82.,5151 Pop, R. M.; Buzoianu, A. D.; Rati, I. V.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca
2014, 42, 545.,5454 Kozarski, M.; Klaus, A.; Vunduk, J.; Zizak, Z.; Niksic, M.; Jakovljevic, D.; Vrvic, M. M.; Van Griensven, L. J. L. D.; Food Funct. 2015, 6, 1875.,6363 Černá, M.; Barros, A. S.; Nunes, A.; Rocha, S. M.; Delgadillo, I.; Čopíková, J.; Coimbra, M. A.; Carbohydr. Polym. 2003, 51, 383.,6464 Kačuráková, M.; Mathlouthi, M.; Carbohydr. Res. 1996, 284, 145. Depending on the water deficit intensity, changes in carbon partitioning in the leaf and the plant as a whole may occur6565 Chaves, M. M.; J. Exp. Bot.
1991, 42, 1. and modify the leaf's carbohydrate levels.2020 Praxedes, S. C.; DaMatta, F. M.; Loureiro, M. E.; Ferrão, M. A. G.; Cordeiro, A. T.; Environ. Exp. Bot.
2006, 56, 263.,6565 Chaves, M. M.; J. Exp. Bot.
1991, 42, 1.
66 Chaves Filho, J. T.; Stacciarini-Seraphin, E.; Rev. Bras. Bot.
2011, 24, 199.-6767 Pinheiro, C.; Chaves, M. M.; Ricardo, C. P.; J. Exp. Bot. 2001, 52, 1063. There are reports in the literature about the decrease in starch content under water stress while other carbohydrates increase, representing the osmotic adjustment strategy6868 Gzik, A.; Environ. Exp. Bot.
1996, 36, 29.
69 Lawlor, D. W.; Cornic, G.; Plant, Cell Environ. 2002, 25, 275.-7070 Pelleschi, S.; Rocher, J. P.; Prioul, J. L.; Plant, Cell Environ.
1997, 20, 493. and suggesting that the bands of leaves from non-irrigated plants are the derivative products of carbohydrate.
The distribution of the estimated class values, for both the calibration and validation datasets of the authentication model, are presented in Figure 4. For validation, the test subset was composed of 5 irrigated and 9 non-irrigated extracts. A separation among irrigated and non-irrigated classes occurred. The threshold value was estimated using Bayes' theorem and it is shown as the horizontal line with values of 0.5556 separating the irrigated class and 0.4444 the non-irrigated one. Sensitivity and specificity could be also determined.
OSC-PLS-DA classification for the FTIR spectra of crude leaf extracts from (a) irrigated and (b) non-irrigated C. arabica cv. IAPAR 59 plants considering (●) extract calibration of leaves from irrigated plants; (■) extract calibration of leaves from non-irrigated plants; (○) extract validation of leaves from irrigated plants; (□) extract validation of leaves from non-irrigated plants.
Sensitivity is the model's ability to classify the validation samples belonging to a particular class. If the model classifies all samples in a given class correctly, the sensitivity attributed to this class is equal to 1.4444 Alves, F. C. G. B. S.; Valderrama, P.; Anal. Methods 2015, 7, 9702. For the model, the sensitivities were 1.00 and 0.89 for irrigated and non-irrigated classes, respectively. Based on these results, the model classified correctly all the irrigated extracts and made only one incorrect prediction to the non-irrigated class. Specificity is related to the incorrect prediction of validation samples of other classes in a particular class. Thus, if the model does not present error in predicting a sample class, this model presents a specificity equal to 1.4444 Alves, F. C. G. B. S.; Valderrama, P.; Anal. Methods 2015, 7, 9702. The irrigated extracts present specificity equal to 0.89 due to the erroneous prediction of one non-irrigated extract as belonging to the irrigated class. The non-irrigated class presents specificity equal to 1.00, meaning that no irrigated extract was classified in the non-irrigated class.
Development of a predicted model capable of discriminating stratified leaves into an inferior layer (40 cm, self-shaded) and superior layer (> 80 cm, high light exposed) separately for each irrigation mode was attempted. The model prediction and classification performance for non-irrigated layers was not satisfactory. This result could be explained by the fact of the total leaf area of non-irrigated plants was drastically affected when compared with areas of plants grown under irrigation, allowing a more homogeneous solar radiation access to all plant layers. However, the model prediction and classification performance for irrigated layers was satisfactory. Spectral datasets of crude extracts from irrigated layers were composed of: 17 calibration (9 harvested from the inferior layer and 8 from the superior layer) and 7 validation (3 harvested from the inferior layer and 4 from the superior layer), selected using the Kennard-Stone algorithm.4343 Kennard, R. W.; Stone, L. A.; Technometrics 1969, 11, 137. The FTIR spectra were preprocessed using the second-derivative (second order polynomial with a seventeen point window) and mean-centered before applying the OSC algorithm to support the selection of model dimensionality. Random cross-validation was employed and the RMSECV (0.2371) was calculated to establish the best model. Two latent variables were needed to develop a classification model, representing 61.60% of the variance explained in X and 89.69% of the variance in the y (inferior layer samples were categorized in class 0 and superior layer samples were into class 1).
The distribution of the estimated class values, for calibration and validation datasets from irrigated layers are shown in Figure 5. In this case, the sensitivities are 1.00 and 0.75 for inferior and superior layers' classes, respectively. Based on these results, the model classified all the extracts of leaves from the inferior layer while having one incorrect prediction for those from the superior layer. Extracts of leaves harvested from the inferior canopy layer presented a specificity equal to 0.75 due to the prediction of one extract of superior layer in the inferior layer class. The superior layer class presents specificity equal to 1.00 meaning no extract from the inferior canopy layer was classified in the superior layer class.
OSC-PLS-DA classification for FTIR spectra of (a) inferior and (b) superior layer from C. arabica cv. IAPAR 59 crude leaf extracts considering (●) extract calibration of leaves from inferior layer; (■) extract calibration of leaves from superior layer; (○) extract validation of leaves from inferior layer; (□) extract validation of leaves from superior layer.
The separation between leaves from inferior and superior canopy layers can be observed into LV1 score plot (Figure 6a). Extracts from superior layer were discriminated by the positive part of LV1 and extracts from inferior layer leaves were discriminated by the negative one. The loadings plot (Figure 6b) shows the most intense FTIR bands which contain the major contribution to discriminate between sample groups. In this case, the bands at 1699, 1657, 1075, and 1031 cm-1 contributed to classify the extracts from the superior layer, because they have positive loadings for LV1. For extracts from the inferior layer, the band at 1679 cm-1 is the most important due to the negative loadings value.
(a) LV1 scores plot of the OSC-PLS-DA model and (b) loading plots of LV1 built with preprocessed and OSC filtered FTIR spectra of crude leaf extracts from C. arabica cv. IAPAR 59 plants. Leaves were collected from (●) inferior and (■) superior layers of the irrigated plants.
Bands around 1699 and 1657 cm-1 may be attributed to the presence of caffeine,3030 Wang, N.; Fu, Y.; Lim, L. T.; J. Agric. Food Chem. 2011, 59, 3220.,7171 Ohnsmann, J.; Quintás, G.; Garrigues, S.; de la Guardia, M.; Anal. Bioanal. Chem. 2002, 374, 561. suggesting that the compound may be the one responsible for the superior layer discrimination. This result agrees with the findings by Delaroza et al.,99 Delaroza, F.; Rakocevic, M.; Malta, G. B.; Bruns, R. E.; Scarminio, I. S.; J. Braz. Chem. Soc. 2014, 25, 1929. in which sun exposed crude extracts from coffee leaves show higher amount of caffeine, suggesting greater stress suffered by leaves exposed to direct sunlight. Moreover, the superior layers were also discriminated by the 1075 and 1031 cm-1 bands indicating the presence of carbohydrates, as discussed before.4848 Zavoi, S.; Fetea, F.; Ranga, F.; Pop, R. M.; Baciu, A.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca 2011, 39, 82.,5151 Pop, R. M.; Buzoianu, A. D.; Rati, I. V.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca 2014, 42, 545.,5454 Kozarski, M.; Klaus, A.; Vunduk, J.; Zizak, Z.; Niksic, M.; Jakovljevic, D.; Vrvic, M. M.; Van Griensven, L. J. L. D.; Food Funct. 2015, 6, 1875.,6363 Černá, M.; Barros, A. S.; Nunes, A.; Rocha, S. M.; Delgadillo, I.; Čopíková, J.; Coimbra, M. A.; Carbohydr. Polym. 2003, 51, 383.,6464 Kačuráková, M.; Mathlouthi, M.; Carbohydr. Res. 1996, 284, 145. The band at 1679 cm-1 is consistent with C=C stretching of lipids,66 Wang, N.; Lim, L. T.; J. Agric. Food Chem. 2012, 60, 5446.,5555 Schulz, H.; Baranska, M.; Vib. Spectrosc. 2006, 43, 13. suggesting that these compounds may be responsible for the inferior layer discrimination. This agrees with that reported by Delaroza et al.,99 Delaroza, F.; Rakocevic, M.; Malta, G. B.; Bruns, R. E.; Scarminio, I. S.; J. Braz. Chem. Soc. 2014, 25, 1929. in which the lipid contents were higher for the self-shaded leaves than those exposed to direct sunlight, indicating the efficient lipid protective role of secondary metabolism products in the shade.
Conclusions
PCA showed the effects of mixture design solvents for metabolite extraction in coffee leaves, indicating that extractor solvents containing ethanol are most suitable for FTIR analysis. This occurs because the crude extracts obtained, using solvents containing ethanol, presented various functional groups besides those of bands indicated in the extracts without the presence of ethanol. OSC-PLS-DA applied to FTIR spectra discriminated and classified coffee leaves subject to different environmental conditions. This tool indicated discrimination between leaf extracts, under different irrigation conditions, by bands related to carbohydrates. In addition, lipid bands, such as carotenoids, and bands relative to amino acids were also important for discrimination. The modification of such compound levels may indicate an osmotic adjustment strategy developed by coffee plants of IAPAR 59 cultivar. Samples harvested at different heights, i.e. receiving different intensities of solar radiation, have bands suggesting that discrimination may be attributed to caffeine, carbohydrates and lipids. So, these metabolites in coffee leaves appear to be quite sensitive to environmental stress. Furthermore, FTIR spectroscopy coupled with chemometric tools showed to be a useful, powerful and simple technique to discriminate and classify coffee leaves subjected to different environmental conditions.
Acknowledgments
The authors acknowledge financial support of CNPq, CAPES, Consórcio Brasileiro de Pesquisa e Desenvolvimento de Cafés, INCT - Energia e Ambiente and Araucária Foundation as well as for the scholarships granted. The authors also thank the Agronomic Institute of Paraná (IAPAR) for supplying the coffee leaf samples used in this work.
References
-
1Davis, A. P.; Govaerts, R.; Bridson, D. M.; Stoffelen, P.; Bot. J. Linn. Soc. 2006, 152, 465.
-
2de Lima, R. B.; Santos, T. B.; Vieira, L. G. E.; Ferrarese, M. L. L.; Ferrarese-Filho, O.; Donatti, L.; Boeger, M. R. T.; Petkowicz, C. L. O.; Carbohydr. Polym. 2013, 93, 135.
-
3Komes, D.; Vojvodić, A. In Processing and Impact on Antioxidants in Beverages; Preedy, V. R., ed.; Academic Press, Elsevier Inc.: Waltham, USA, 2014, ch. 9.
-
4Barbin, D. F.; Felicio, A. L. S. M.; Sun, D. W.; Nixdorf, S. L.; Hirooka, E. Y.; Food Res. Int. 2014, 61, 23.
-
5Marquetti, I.; Link, J. V.; Lemes, A. L. G.; Scholz, M. B. S.; Valderrama, P.; Bona, E.; Comput. Electron. Agric. 2016, 121, 313.
-
6Wang, N.; Lim, L. T.; J. Agric. Food Chem. 2012, 60, 5446.
-
7Campa, C.; Mondolot, L.; Rakotondravao, A.; Bidel, L. P. R.; Gargadennec, A.; Couturon, E.; La Fisca, P.; Rakotomalala, J. J.; Jay-Allemand, C.; Davis, A. P.; Ann. Bot. 2012, 110, 595.
-
8Clifford, M. N.; Kirkpatrick, J.; Kuhnert, N.; Roozendaal, H.; Salgado, P. R.; Food Chem. 2008, 106, 379.
-
9Delaroza, F.; Rakocevic, M.; Malta, G. B.; Bruns, R. E.; Scarminio, I. S.; J. Braz. Chem. Soc. 2014, 25, 1929.
-
10Domingues Júnior, A. P.; Shimizu, M. M.; Moura, J. C. M. S.; Catharino, R. R.; Ramos, R. A.; Ribeiro, R. V.; Mazzafera, P.; Photochem. Photobiol. 2012, 88, 928.
-
11Mondolot, L.; la Fisca, P.; Buatois, B.; Talansier, E.; de Kochko, A.; Campa, C.; Ann. Bot. 2006, 98, 33.
-
12Zheng, X. Q.; Ashihara, H.; Plant Sci. 2004, 166, 807.
-
13Brum, C. N. F.; Melo, E. F.; Barquero, L. O. B.; Alves, J. D.; Chalfun-Júnior, A.; Coffee Sci. 2013, 8, 140.
-
14Lima, R. B.; Santos, T. B.; Vieira, L. G. E.; Ferrarese, M. L. L.; Farrarese-Filho, O.; Donatti, L.; Boeger, M. R. T.; Petkowicz, C. L. O.; Carbohydr. Polym. 2014, 112, 686.
-
15Zhao, J.; Hu, C.; Zeng, J.; Zhao, Y.; Zhang, J.; Chang, Y.; Li, L.; Zhao, C.; Lu, X.; Xu, G.; Metabolomics 2014, 10, 805.
-
16Chaves, A. R. M.; Ten-Caten, A.; Pinheiro, H. A.; Ribeiro, A.; DaMatta, F. M.; Trees 2008, 22, 351.
-
17Da Matta, F. M.; Ramalho, J. D. C.; Braz. J. Plant Physiol. 2006, 18, 55.
-
18Jones, M. M.; Turner, N. C.; Plant Physiol. 1978, 61, 122.
-
19Hsiao, T. C.; Acevedo, E.; Fereres, E.; Henderson, D. W.; Philos. Trans. R. Soc., B 1976, 273, 479.
-
20Praxedes, S. C.; DaMatta, F. M.; Loureiro, M. E.; Ferrão, M. A. G.; Cordeiro, A. T.; Environ. Exp. Bot. 2006, 56, 263.
-
21Gobbo-Neto, L.; Lopes, N. P.; Quim. Nova 2007, 30, 374.
-
22Dunn, W. B.; Ellis, D. I.; TrAC, Trends Anal. Chem. 2005, 24, 285.
-
23Ivanišević, J.; Thomas, O. P.; Lejeusne, C.; Chevaldonné, P.; Pérez, T.; Metabolomics 2011, 7, 289.
-
24Viaene, J.; Goodarzi, M.; Dejaegher, B.; Tistaert, C.; Hoang Le Tuan, A.; Nguyen Hoia, N.; Chau Van, M.; Quetin-Leclercq, J.; Vander Heyden, Y.; Anal. Chim. Acta 2015, 877, 41.
-
25Ellis, D. I.; Harrigan, G. G.; Goodacre, R. In Metabolic Profiling: Its Role in Biomarker Discovery and Gene Function Analysis; Harrigan, G. G.; Goodacre, R., eds.; Kluwer Academic Publishers: London, UK, 2003, ch. 7.
-
26Butler, H. J.; McAinsh, M. R.; Adams, S.; Martin, F. L.; Anal. Methods 2015, 7, 4059.
-
27Ord, J.; Butler, H. J.; McAinsh, M. R.; Martim, F. L.; Anal. Methods 2016, 141, 2896.
-
28Gallardo-Velázquez, T.; Osorio-Revilla, G.; Zuñiga-de Loa, M.; Rivera-Espinoza, Y.; Food Res. Int. 2009, 42, 313.
-
29Hendriks, M. M. W. B.; Cruz-Juarez, L.; de Bont, D.; Hall, R. D.; Anal. Chim. Acta 2005, 545, 53.
-
30Wang, N.; Fu, Y.; Lim, L. T.; J. Agric. Food Chem. 2011, 59, 3220.
-
31Wold, S.; Esbensen, K.; Geladi, P.; Chemom. Intell. Lab. Syst. 1987, 2, 37.
-
32García-Reiriz, A. G.; Olivieri, A. C.; Teixidó, E.; Ginebreda, A.; Tauler, R.; Environ. Sci.: Processes Impacts 2014, 16, 124.
-
33Ballabio, D.; Consonni, V.; Anal. Methods 2013, 5, 3790.
-
34Matsunaga, F. T.; Tosti, J. B.; Androcioli-Filho, A.; Brancher, J. D.; Costes, E.; Rakocevic, M.; SpringerPlus 2016, 5, 2015.
-
35Snyder, L. R.; Carr, P. W.; Rutan, S. C.; J. Chromatogr. A 1993, 656, 537.
-
36Barker, M.; Rayens, W. J.; J. Chemom. 2003, 17, 166.
-
37Geladi, P.; Kowalski, B. R.; Anal. Chim. Acta 1986, 185, 1.
-
38Farrés, M.; Piña, B.; Tauler, R.; Metabolomics 2015, 11, 210.
-
39Carlos, C.; Maretto, D. A.; Poppi, R. J.; Sato, M. I. Z.; Ottoboni, L. M. M.; Microchem. J. 2011, 99, 15.
-
40Wold, S.; Antti, H.; Lindgren, F.; Öhman, J.; Chemom. Intell. Lab. Syst. 1998, 44, 175.
-
41Bylesjö, M.; Rantalainen, M.; Cloarec, O.; Nicholson, J. K.; Holmes E.; Trygg, J.; J. Chemom. 2006, 20, 341.
-
42de Almeida, M. R.; Correa, D. N.; Rocha, W. F. C.; Scafi, F. J. O.; Poppi, R. J.; Microchem. J. 2013, 109, 170.
-
43Kennard, R. W.; Stone, L. A.; Technometrics 1969, 11, 137.
-
44Alves, F. C. G. B. S.; Valderrama, P.; Anal. Methods 2015, 7, 9702.
-
45de Almeida, M. R.; Fidelis, C. H. V.; Barata, L. E. S.; Poppi, R. J.; Talanta 2013, 117, 305.
-
46Pérez, N. F.; Ferré, J.; Boqué, R.; Chemom. Intell. Lab. Syst. 2009, 95, 122.
-
47Movasaghi, Z.; Rehman, S.; Rehman, I. U.; Appl. Spectrosc. Rev. 2008, 43, 134.
-
48Zavoi, S.; Fetea, F.; Ranga, F.; Pop, R. M.; Baciu, A.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca 2011, 39, 82.
-
49Kushwaha, K.; Saxena, J.; Tripathi, B. K.; Agarwal, M. K.; J. Biosci. Biotechnol. 2014,3, 253.
-
50Craig, A. P.; Franca, A. S.; Oliveira, L. S.; Food Chem. 2012, 132, 1368.
-
51Pop, R. M.; Buzoianu, A. D.; Rati, I. V.; Socaciu, C.; Not. Bot. Horti Agrobot. Cluj-Napoca 2014, 42, 545.
-
52Ribeiro, J. S.; Salva, T. J.; Ferreira, M. M. C.; J. Food Qual. 2010, 33, 212.
-
53Terrile, A. E.; Marcheafave, G. G.; Oliveira, G. S.; Rakocevic, M.; Bruns, R. E.; Scarminio, I. S.; J. Braz. Chem. Soc. 2016, 7, 1254.
-
54Kozarski, M.; Klaus, A.; Vunduk, J.; Zizak, Z.; Niksic, M.; Jakovljevic, D.; Vrvic, M. M.; Van Griensven, L. J. L. D.; Food Funct. 2015, 6, 1875.
-
55Schulz, H.; Baranska, M.; Vib. Spectrosc. 2006, 43, 13.
-
56Lóránd, T.; Deli, J.; Molnár, P.; Tóth, G.; Helv. Chim. Acta 2002, 85, 1691.
-
57Sandman, G.; Arch. Biochem. Biophys. 2009, 483, 169.
-
58Bartley, G. E.; Scolnik, P. A.; Plant Cell 1995, 7, 1027.
-
59Wollman, F. A.; EMBO J. 2001, 20, 3623.
-
60Szabò, I.; Bergantino, E.; Giacometti, G. M.; EMBO Rep. 2005, 6, 629.
-
61Lechinoski, A.; de Freitas, J. M. N.; Castro, D. S.; Lobato, A. K. S.; de Oliveira Neto, C. F.; Cunha, R. L. M.; Rev. Bras. Bioc. 2007, 5, 927.
-
62Pavia, D. L.; Lampman, G. M.; Kriz, G. S.; Introduction to Spectroscopy: A Guide of Students of Organic Chemistry, 3rd ed.; Thomson Learning, Books/Cole: Seneca, Mexico, 2011.
-
63Černá, M.; Barros, A. S.; Nunes, A.; Rocha, S. M.; Delgadillo, I.; Čopíková, J.; Coimbra, M. A.; Carbohydr. Polym. 2003, 51, 383.
-
64Kačuráková, M.; Mathlouthi, M.; Carbohydr. Res. 1996, 284, 145.
-
65Chaves, M. M.; J. Exp. Bot. 1991, 42, 1.
-
66Chaves Filho, J. T.; Stacciarini-Seraphin, E.; Rev. Bras. Bot. 2011, 24, 199.
-
67Pinheiro, C.; Chaves, M. M.; Ricardo, C. P.; J. Exp. Bot. 2001, 52, 1063.
-
68Gzik, A.; Environ. Exp. Bot. 1996, 36, 29.
-
69Lawlor, D. W.; Cornic, G.; Plant, Cell Environ. 2002, 25, 275.
-
70Pelleschi, S.; Rocher, J. P.; Prioul, J. L.; Plant, Cell Environ. 1997, 20, 493.
-
71Ohnsmann, J.; Quintás, G.; Garrigues, S.; de la Guardia, M.; Anal. Bioanal. Chem. 2002, 374, 561.
Publication Dates
-
Publication in this collection
Jan 2018
History
-
Received
14 Mar 2017 -
Accepted
04 July 2017