ORIGINAL RESEARCH article

Front. Nutr., 04 August 2022

Sec. Food Chemistry

Volume 9 - 2022 | https://doi.org/10.3389/fnut.2022.968783

Impact of agro-forestry systems on the aroma generation of coffee beans

  • 1. Food and Pharmaceutical Engineering Institute, Guiyang University, Guiyang, China

  • 2. Division of Food Sciences, University of Nottingham, Nottingham, United Kingdom

  • 3. Centro Agronómico Tropical de Investigación y Enseñanza, CATIE, Turrialba, Costa Rica

Abstract

A long experiment has been established since 2000 at CATIE (Tropical Agricultural Research and Higher Education Center), Turrialba, Costa Rica. Twenty agro-forestry systems with different shade types and managements (organic and non-organic) consisting of an incomplete randomized block-design with shade tree as main effect and subplots represented by management were set up. The effects of different managements and shade types on the aroma and color generation of roasted coffee beans were investigated. The total protein content was significantly higher (P < 0.05) under the intensive conventional (IC) (168 g/Kg) and intensive organic (IO) (167 g/Kg) managements than under the moderate conventional (MC) (153 g/Kg in IC vs. MC group, 157 g/Kg in MC vs. IO group). Comparing with the moderate conventional (MC) management, the intensive organic (IO) management had a stronger ability to generate more flavor and color. The total protein content was significantly higher (P < 0.05) under the full sun system (172 g/Kg) than under the shaded (159 g/Kg) and Erythrina system (155 g/Kg), under the service system (165 g/Kg) than under the timber system (146 g/Kg), under the legume timber system (170 g/Kg) than under the non-legume timber system (152 g/Kg). The full sun system had a greater flavor generation and color after roasting. Comparing with the timber system, the service system produced roasted beans with the more flavor and color. Comparing with the non-legume shade tree, the legume shade tree improved the performance of flavor and color in the roasted coffee beans.

Introduction

Coffee is widely appreciated and consumed due to its pleasant flavor generated during the roasting process (1). In addition, expert coffee tasters assess the quality of coffee, to a greater extent, by its aroma and taste, and the highest quality coffee beans can obtain a considerable price premium (2). A large number of investigations, therefore, have been carried out to evaluate the types and quantities of volatile compounds as well as its formation principals in the roasted coffee beans. Most volatile compounds formed during roasting are derived from a range of chemical precursors in the green coffee beans, these include sugars, amino acids, organic acids and phenolic compounds, and these react through many complex reactions including the Maillard reaction, degradation of sugars and lipid oxidation (37). Nine specific pathways to generate the volatile compounds during roasting of coffee have been outlined by Dart and Nursten (2) in 1985. There were the Maillard reaction, Strecker degradation, trigonelline degradation, phenolic acids degradation, lipid degradation, sugar degradation, breakdown of sulphur amino acid, breakdown of hydroxyl amino acids and proline, and hydroxyproline degradation. Moreover, the dominating volatile compounds, according to their functional groups, have been divided into ten classes by Dart and Nursten (2). They included sulphur compounds, pyrazines, pyridines, pyrroles, oxazoles, furans, aldehydes, ketones, phenols, and other heterocyclic compounds.

There are many sourcess affecting the aroma generation of roasted coffee beans: the origin of coffee beans (i.e., Arabica, Robusta) may contribute to the differences (8); seasonal variation may also play an essential role in differences (9, 10); additionally, differences in geographical locations (11), differences in post-harvest processing (wet or dry processing) and aging before roasting also affect the final aroma generation (1). Furthermore, roasting profile (time-temperature) and the roaster models are aslo significant contributors to the differences of final coffee flavor (12, 13).

Coffee roasting plays an essential role in the formation of organoleptic properties including flavor and color. During the roasting process, moisture content decreases and many chemical reactions occur, along with crucial variations in terms of flavor, color, volume, weight, bean pop, pH, density and presence of volatile compounds (14).

Gas chromatograph is widely used to analyse the coffee volatile composition followed by mass spectrometer (15) or other specific detectors [flame ionization detectors (16), nitrogen-phosphorous detectors, photo-ionization detectors] due to its sensitivity and operability. It can separate and identify a complex mixture of aromas in one operation as well as quantify the volatiles at extremely low concentrations (2).

Numerous researches have been carried to identify the volatile compounds of coffee, explore the formation mechanisms of these pleasant aromas and investigate the relationship between the composition of green coffee beans and aroma generation of roasted coffee beans (1721), but only a limited number of studies have actually connected these to the growth conditions, levels and types of fertilizers, the species of shade trees. As a result, it is necessary for coffee scientists to explore the interactions and detail the specific effects of fertilizer levels and shade types on aroma generation of roasted coffee beans. The aim of this research was to investigate the impact of agro-forestry systems including fertilizer levels and shade types on the aroma generation of roasted coffee bean.

Materials and methods

Agro-forestry methodologies

The experiment was established in 2000 at CATIE (Tropical Agricultural Research and Higher Education Center), Turrialba, Costa Rica (9°53'44” N, 83°40'7” W, CATIE, Turrialba, Costa Rica), which is defined as a low altitude (600 m above sea level), wet coffee zone without a marked dry season. Average annual rainfall, temperature, relative humidity and solar radiation was 2,915 mm/year, 22°C, 90.2%, and 15.9 MJ/m2/year (2000–2013).

Twenty agroforestry systems with different shade types and managements consisting of an incomplete randomized block-design with shade tree as main effect and subplots represented by management were set up (Table 1). For each system, three replicates were established. Shade type [initially 417 trees per ha-1 (6 × 4 m2 spacing)] consisted of timber and service tree species with contrasting characteristics (Table 2). Trees were progressively thinned to maintain a reasonable shade environment for coffee production (Table 3).

Table 1

Shade types*ETCC+TE+TC+EFull Sun
Managements**ICICICIC
MCMCMCMCMCMCMC
IOIOIOIOIOIO
LOLOLO

Agroforestry systems with main plot (Shade type) and subplot (Management) treatments.

*

E, Erythrina poepiggiana; C, Chloroleucon eurycyclum; T, Terminalia amazonia.

**

IC, Intensive conventional; MC, Moderate conventional; IO, Intensive organic, LO, Low organic; (n = 3).

Table 2

SpeciesPhenologyCanopyN-fixerUse
Erythrina poepiggiana (E)EvergreenLow compactYesService
Chloroleucon eurycyclum (C)Deciduous*High spreadingYesTimber
Terminalia amazonia (T)Deciduous*High compactNoTimber

Characteristics of shade trees (22).

*

Deciduous for about 20–30 days per year.

Table 3

Agroforestry systemTree density per ha−1
SystemTree species200820112013
Monocultures
EE360285241
CC38115465
TT31716773
Polycultures
C+EC18310045
E181134115
C+TC1667739
T1707734
E+TE147143109
T1588134

Mean shade tree density after thinning.

E, Erythrina poepiggiana; C, Chloroleucon eurycyclum; T, Terminalia Amazonia.

Intensive conventional (IC) Erythrina trees were biannually pollarded to a 1.8–2.0 m main trunk. Whilst this is normal practice in Costa Rica, Muschler (24) found that coffee quality benefited from increased Erythrina shade levels, therefore, for all the other treatments with Erythrina, trees were pollarded to 4 m leaving three branches for partial shade. Temporary shade was planted in the form of Ricinus in organic treatments. This took place 1 year after the coffee plants, to improve coffee plant survival and impede weed growth. Lower branches of the timber trees were pruned annually (year 1–7) to improve stem quality. In all pruning scenarios, pruning residuals from coffee trees and shade trees were left on the ground (trunks were removed). Management consisted of fertilization, weed, disease and pest control, detailed in Table 4.

Table 4

ManagementFertilization N:P:K **Weed controlDisease/Pest control
IC287:20:150Six*
herbicides
3–4*
Fungicides/
insecticides
MC150:10:75Five
herbicides
four
manual
1–4
Fungicides/
insecticides
as required
IO248:205:326Four
Manual
Organic substances as required
LO66:2:44Four
Manual
No

Mean input levels of fertilizers (Kg ha−1year−1) and weed/disease control since 2006, adapted from Haggar et al. (23) and Noponen et al. (25).

IC, Intensive conventional, MC, Moderate conventional, IO, Intensive organic, LO, Low organic. IO fertilization: Chicken manure 10 t ha−1year−1and K-Mag (fertilizer)100 kg ha−1year−1; LO fertilization: Coffee pulp 5 t ha−1year−1.

*

Number of treatments applied per year.

**

Fertilization levels (Kg ha−1year−1) are 7 years means (2003–2009), from the second to fourth year LO systems received the same fertilization as IO ones, due to the site limitations that did not allow organic coffee to establish effectively with lower inputs.

Coffea Arabica L. var. Caturra, was planted at 5000 holes ha−1 with dead plants replaced each year. Two plants per planting hole were planted (local practice) but were treated as one plant in every analysis. The distance between rows and holes were 2 m and 1 m. Sub-plots were 500–600 m2 of which the central 225–300 m2 was studied (100 coffee plants and 24 shade trees).

Coffee plants were manually pruned from 2004 leaving 1–4 resprouts per stump, according to the productive potential of each coffee resprout in the next harvest. Every coffee planting hole thus comprised 1–2 stumps and a total of 1–4 resprouts per stump.

Coffee roasting

All green coffee beans harvested in different agro-forestry systems from CATIE, Costa Rica in both 2013 and 2014 were sent to Edgehill Coffee Company (Warwick, UK) for roasting. The coffee roaster (Roastilino, Fracino) with a small fluidized bed and a nominal capacity of about 200 g of green coffee beans was used and 50 g of green coffee beans were prepared and roasted. The air outlet temperature and roasting time were 166°C and 1 min 50 s, respectively. All coffee beans were taken past first crack and roasted to medium roasting. After roasting, 30 g of different roasted coffee beans were ground into powders using a mini chopper (CH180 mini chopper, 300 W, Kenwood) and a 710 μm sieve was used to obtain even coffee powder samples. All coffee beans and powders were stored at −20°C.

Measurement of color

The colorimeter (ColorQuest XE, HunterLab Inc.) was used to measure the color of the coffee bean (25). The colorimeter was calibrated using the white plate and packed coffee powder was added to a transparent preservative film and then placed on the detection port to obtain data. Results were calculated by the equipment using the Hunter Lab color scale. In this scale, L ranges from 0 (black) to 100 (white), a indicates degree of greenness (for negative a values) and degree of redness (for positive a results), b axis also ranges from negative to positive values indicating, respectively, degree of blueness to yellowness.

Measurement of total protein content

A bicinchoninic acid (BCA) assay kit (23225/23227, Thermo Scientific) was used to measure the total protein content (26).

0.1 g of coffee powders (wet basis, 11% water content) was weighed and 2 mL of chloroform was used to extract lipid with three replicates. After lipid extraction, the solid phase was dried at 60°C for 30 min and placed in a 15 mL centrifuge tube. After that, 1 mL of 2% sodium dodecyl sulfate (SDS) solution was added into the tube to eliminate interference from lipids and the tube was heated at 60°C in a water bath for 30 min. The sample was then vortexed for 1 min and centrifuged at 13,000 g for 3 min. And then, samples were filtered and diluted 100 times using 2% SDS.

According to the instructions of BCA protein assay kit, 0.1 mL of each standard and sample replicate was added into test tubes. 2.0 mL of the working reagent containing 50 parts of BCA Reagent A (sodium carbonate, sodium bicarbonate, bicinchoninic acid and sodium tartrate in 0.1 mol/L sodium hydroxide) and 1 part of BCA Reagent B (4% cupric sulfate) was added to each tube and mixed well. All tubes were covered and incubated at 37°C for 30 min, after incubation, tubes were cooled to room temperature. The absorbance of the samples was measured at 562 nm by the ultraviolet spectrophotometer (Evolution 350, Thermo Scientific) and subtracted the average 562 nm absorbance measurement of the Blank standard replicates. A standard curve was prepared by plotting the average blank-corrected 562 nm measurements for each bovin serum albumin (BSA) standard vs. its concentration in μg/mL. Finally, the standard curve was used to determine the protein concentration of each sample.

Measurement of sucrose

The sucrose standard was prepared at a concentration of 10, 20, 40, and 50 mg/L. 0.1 g of coffee powder was placed in a 50 mL centrifuge tube with 20 mL of boiling water and vortexed for 5 min. Samples then were centrifuged at 1,600 g for 10 min at room temperature. After centrifugation, the liquid phase was transferred into a new glass vial. The above processes were repeated three times for extracting sucrose completely. The mixture was cooled to room temperature and then filtered using a syringe filter (0.45 μm, hydrophilic nylon syringe filter, Millipore Corporation). The final extract was diluted with water (1:1) prior to Liquid chromatography-mass spectrometry (LC/MS) analysis [the method was modified from Ky et al. (27), Mccusker et al. (28), Perrone et al. (29)].

The LC equipment (1,100 Series, Agilent) consisted of a degasser (G1322A, Agilent), a pump (G1312A, Agilent), an auto-sampler (G1313A, Agilent). This LC system was interfaced with a Quattro Ultima mass spectrometer (Micromass, UK Ltd.) fitted with an electrospray ion source.

The Kromasil 5 ODS (C18) column (250 × 3.20 mm, 5 μm, Phenomenex) was used to separate sucrose and the column was maintained at room temperature. The mobile phase consisted of 0.3% aqueous formic acid (eluent A) and methanol (eluent B) which was delivered at a flow rate of 0.4 ml/min. Initially the solvent was 25% B. After the injection, this proportion was changed immediately to 60% B until the end of the run at 8 min. Between injections, the column was re-equilibrated with 25% B for 4 min.

The electrospray ionization source was operated in the negative mode (3.5 KV) from 0.0 to 4.0 min to generate the formic acid adduct of sucrose (M+HCOO) and in the positive mode (3.5 KV) from 4.0 to 8.0 min to generate caffeine (M+H)+ ions. Desolvation gas (N2) flow rate was 500 L/h and desolvation temperature was 400°C. The mass spectrometer was operated in the selected ion recording (SIR) model to measure sucrose. Identification of compounds was performed by comparing retention time and molecular weight of the respective standards. Quantification was achieved by comparing peak areas of samples and standards using Masslynx software (version 4.0, Waters) (the method was modified from 24 and 25).

Measurement of volatile compounds

2 g of each coffee sample powders and 2.5 μL of 1% heptanone (internal standard) were placed GC vials for analysis (30).

A GC-ISQ (Thermo Electron Corporation Inc.) was used to detect the volatile compounds of coffee in full scan mode over the mass range 20–200 m/z (ion source 200°C). The sample was incubated at 55°C for 5 min with shaking using the autosampler. A 50/30 μm DVB/CAR/PDMS (Divinylbenzene/Carboxen/Polydimethylsiloxane) was the fiber sampled the volatiles for 3.5 min followed by a 1.5 min desorption (inlet temperature 200°C). A splitless mode was used and the constant carrier pressure was at 103 kPa.

The GC oven was held at 40°C for 5 min, and then increased from 40 to 240°C at a rate of 5°C/min, and held at 240°C for 2 min. Separation was carried out on a ZB-WAX Capillary GC Column (Length 30 m, inner diameter 0.25 mm, and film thickness 0.25 μm; Phenomenex Inc., Macclesfield, UK). The identity and quantity of volatile compounds was determined using Thermo Xcalibur software. The mass spectrum was used to identify the volatile compounds by comparing with the library (NIST/EPA/NIH Mass Spectral Library. Version 2.0, Faircom Corporation, U.S.) in Xcalibur software. The peak areas of compounds and the internal standard were determined using Thermo Xcalibur software (30). The following formula could be used to calculate the concentration of volatile compounds:

Statistical analysis

The experimental design consisted of shade types as main plot and subplots represented by managements but with an unbalanced structure since not all managements were represented under all Shade types (Table 1). Therefore, the specific pre-planned contrast models (Modified from 19) (Table 5) were used to compare the differences between different managements and shade types.

Table 5

ContrastTreatments compared
Management
IC vs. MCIC(FS, E, T, CE) vs. MC(FS, E, T, CE)
MC vs. IOMC(E, T, C, CE, CT, ET) vs. IO(E, T, C, CE, CT, ET)
IO vs. LOIO(E, CE) vs. LO(E, CE)
IC vs. IOIC(E, T, CE) vs. IO(E, T, CE)
Shade type
FS vs. shaded

Erythrina vs. FS*
FS(IC, MC) vs. E(IC, MC) + T(IC, MC) + CE(IC, MC)
E(IC, MC) vs. FS(IC, MC)
Service vs. timber treesE(MC, IO) vs. T(MC, IO) + C(MC, IO) + TC(MC, IO)
Legume timber vs. non-legume timberC(MC, IO) vs. T(MC, IO)

Principal contrasts used in the analysis of shade type and management effects.

IC, Intensive conventional; MC, Moderate conventional; IO, Intensive organic; LO, Low organic; FS, Full sun. E, Erythrina poepiggiana; C, Chloroleucon eurycyclum; T, Terminalia Amazonia; CE, Chloroleucon eurycyclum and Erythrina poepiggiana; CT, Chloroleucon eurycyclum and Terminalia Amazonia; ET, Erythrina poepiggiana and Terminalia Amazonia.

*

Erythrina was regarded as a low canopy tree with low shade cover and compared with full sun.

Data was analyzed using One-Way ANOVA contrasts model with SPSS statistics software (Version 22nd, IBM), treatment as contrast factor to analyze eight contrast groups as shown in Table 5. Significant value (P < 0.05) in contrast tests was used to decide if there was a significant difference in each contrast model and the standard error of each contrast could be obtained from contrast tests table. In addition, the mean value of every management and shade type under each contrast model was calculated.

Results and discussion

Impact of fertilizer levels and shade types on the total protein and sucrose contents of green coffee bean

During the coffee roasting progress, carbohydrates especially sucrose are degraded, dehydrated and interacted with protein to generate a large number of volatile compounds (31, 32), which contributes to the formation of organoleptic properties of roasted coffee beans including flavor, aroma and color (14). As a result, the sucrose and protein were analyzed due to the main contribution to the generation of volatile compounds.

The results of total protein content were shown in Table 6. The total protein content was significantly lower (P < 0.05) under the MC (153, 11.3 g/Kg) than under the IC (168, 12.2 g/Kg). Furthermore, the total protein content was significantly higher (P < 0.05) under the IO (167, 12.4 g/Kg) than under the MC (157, 11.8 g/Kg). Nevertheless, no significant differences (P > 0.05) were detected in the total protein content between the IO and LO, and the IC and IO.

Table 6

ContrastTotal protein
(g/Kg)
Sucrose
(g/Kg)
ManagementsMeanS.E.Dp-valueMeanS.E.Dp-value
IC vs. MC167.6217.250.00191.529.120.544
152.5390.09
MC vs. IO157.1121.310.0188.1211.250.053
166.5192.25
IO vs. LO162.3912.180.9486.976.910.915
161.9287.33
IC vs. IO166.9214.730.41489.898.210.328
162.7987.22
Shade types
FS vs. S171.8829.870.01191.7316.050.647
159.0590.49
E vs. FS155.3612.11<0.00190.367.360.68
171.7891.73
Ser. vs. Tim.165.3330.24<0.00188.5016.180.984
146.2988.56
LT vs. NLT170.3312.21<0.00184.216.870.077
152.2690.14

Contrast results for total protein and sucrose contents (g/Kg) by dry weight of green coffee beans.

Values were presented as mean, standard error of the contrast difference (S.E.D), and significance of the difference (P-value), P-values <0.05 were shown in bold. IC, Intensive conventional; MC, Moderate conventional; IO, Intensive organic; LO, Low organic; FS, Full sun; S, Shaded; E, Erythrina poepiggiana; Ser., Service, Tim, Timber; LT, Legume timber; NLT, Non-legume timber.

For the different shade types, the total protein content of green coffee beans was significantly higher (P < 0.05) under the full sun system (172 g/Kg) than under the shaded system (159 g/Kg) and the Erythrina system (155 g/Kg). In addition, the total protein content of green coffee beans was significantly higher (P < 0.05) under the service system (165 g/Kg) than under the timber system (146 g/Kg), while it was significantly lower (P < 0.05) under the non-legume timber system (152 g/Kg) than under the legume timber system (170 g/Kg).

The sucrose content was shown in the Table 6. There were no significant differences (P > 0.05) detected between the different managements and shade types.

The total protein content was significantly higher (P < 0.05) under the IC and IO managements than under the MC due to the higher input amount of fertilizers especially the nitrogen. After a long period management (about 13 years), the IO management had a similar total protein content in the green coffee beans to the IC management and even a significantly higher (P < 0.05) value than the MC management. Nonetheless, there was no significant difference (P > 0.05) in the total protein content between the IO and LO system. When compared with the IO management, the LO management was combined with the legume species Erythrina and Chloroleucon. On the one hand, the legume species can replenish the lack of N2 due to the low input level by N2 fixation (33); on the other hand, the high pruning intensity like Erythrina may compensate the lower external inputs, as the leaves and branches pruned from the Erythrina shade tree were left on the field (22, 34).

The result of the effect of different shade types on the total protein content suggested that light intensity was the main factor to drive the protein generation in the green coffee beans in our experiment, as the protein biosynthesis in the ribosome is adenosine triphosphate (ATP) dependent (35). As a result, the total protein content was significantly higher (P < 0.05) under the full sun system than under the shaded and Erythrina system due to the higher shade cover and lower light intensity in the shaded and Erythrina systems, moreover, it was also significantly higher (P < 0.05) under the service system than under the timber system due to the lower shade cover and higher light intensity in the service system (the Erythrina was the only service tree in our experiment). Additionally, the total protein content in the green coffee beans was significantly higher (P < 0.05) under the legume timber system than under the non-legume timber system due to the N2 fixation of the legume species to supply the organic N2 for the coffee plant to synthesize the amino acids and protein.

It is well-known that the inorganic carbon is converted to the organic carbon in the plant mainly through photosynthesis, and also, the plant can also absorb the organic carbon fixed by fungi as well as from organic matters on the field (36, 37). In addition, photosynthesis is affected by light intensity, carbon dioxide concentration, water availability, temperature (38, 39). For each contrast group of different managements in our experiment, the shade types have been fixed. Thus, it could be assumed that light intensity and temperature were the same when comparing the different managements in each contrast group. Different managements had no effect on the carbon dioxide concentration due to the same atmospheric conditions. Furthermore, Costa Rica site was considered as the wet coffee zone, which had the adequate water resource. Therefore, no significant differences (P > 0.05) were detected in the sucrose content between the different managements.

In our experiment, there were no significant differences (P > 0.05) found in the sucrose content between the full sun and shaded system, service and timber, and the legume timber and non-legume timber system, which demonstrated that the light intensity, shaded cover and the N2 fixation of shade tree were not the main factors to affect the sucrose content of green coffee beans.

Impact of different fertilizer levels and shade types on the aroma and the color generation of roasted coffee beans

Sixty-four volatile aroma compounds, which had previously been identified in coffee, were detected in the roasted coffee beans (30, 40). These aroma compounds were classified into thirteen groups based on their chemical attributes (three alcohols, four aldehydes, 11 ketones, three acids, one ether, two esters, 14 pyrazines, two pyridines, seven pyrroles, two sulfides, eight furans, five phenolic compounds, and 2 N-containing heterocyclic compounds), and their odors and functional group were described (Table 7). Figure 1 showed the gas chromatogram of roasted coffee beans and identified some main peaks of volatile aroma compounds based on the library (NIST/EPA/NIH Mass Spectral Library. Version 2.0, Faircom Corporation, U.S.). The quantitative data sets of the volatile aroma compounds measured in the current work can be found in Supporting Information.

Table 7

Aroma compoundIDKILiterature KIOdor descriptionFunctional group
12,5-Dimethylfuranm, r721706EtherealFuran
22-Methylfuranm, r893888ChocolateFuran
32-Methylbutanalm, r918930MaltyAldehyde
43-Methylbutanalm, r950934MaltyAldehyde
52,3-Butanedionem, r9831,000Buttery, cheesyKetone
6Dimethyl-Disulphidem, r1,0531,050OnionSulphide
7Hexanalm, r1,0631,064Grassy, green oilyAldehyde
82,3-Pentanedionem, r1,0781,072Oily butteryKetone
92-Vinylfuranm, r1,1121,096Ethereal, rum, cocoa noteFuran
102,3-Hexanedionem, r1,1541,138Buttery, cheesy, sweet, creamyKetone
111-Methylpyrrolem, r1,1711,168White bread, woodyPyrrole
122,4,5-Trimethyloxazolem, r1,2191,200Green, woody, mustyHeterocyclic N
13Pyridinem, r1,2251,213Bitter, astringent, roasted, burntPyridine
142-Pentylfuranm, r1,2391,235Fruity, greenFruan
15Furfuryl methyl etherm, r1,2621,243Nutty, rich, phenolicEther
16Dihydro-2-methyl-3-furanonem, r1,2671,246Sweet, roastedKetone
172-Methylpyrazinem, r1,2821,297Nutty, roasted, chocolatePyrazine
184-Methylthiazolem, r1,3351,312Nutty, greenHeterocyclic N
19Acetoinm, r1,348-Buttery, creamyKetone
202,5-Dimethylpyrazinem, r1,3701,357Nutty, roasted, grassy, cornPyrazine
212,6-Dimethylpyrazinem, r1,3751,362Nutty, sweet, friedPyrazine
222-Ethyl-3-methylpyrazinem, r1,3821,363NuttyPyrazine
232-Methyl-2-cyclopentenonem, r1,3881,366FruityKetone
242-Ethylpyrazinem, r1,3911,370Nutty, roastedPyrazine
252,3-Dimethylpyrazinem,r1,3951,383Nutty, roasted, greenPyrazine
262,3,5-Trimethylpyrazinem, r1,4101,395Nutty, roastedPyrazine
272-Ethyl-5-methylpyrazinem, r1,4161,397Nutty, roastedPyrazine
282-Ethyl-6-methylpyrazinem, r1,4261,420Roasted, hazelnut-likePyrazine
29Propyl pyrazinem, r1,4371,428GreenPyrazine
30Vinyl pyrazinem, r1,4461,434Nutty, greenPyrazine
312,6-Diethylpyrazinem, r1,4581,444Nutty, roastedPyrazine
32Acetic acidm, r1,4711,454SourAcid
33Furfuralm, r1,4791,462Bread, almond, sweetAldehyde
342-Ethyl-3,5-dimethylpyrazinem, r1,4831,464NuttyPyrazine
35Acetoxy acetonem, r1,4911,469ButteryKetone
362-Ethyl-3,6-dimethylpyrazinem, r1,4981,480Nutty, roastedPyrazine
372-Fufurylmethylsulfidem, r1,506-Alliaceous, sulfurousSulphide
382-Acetylfuranm, r1,5141,500Balsamic-sweetFuran
39Pyrrolem, r1,5181,512Nutty, hay-like, herbaceousPyrrole
402,3-Dimethyl-2-cyclopentenonem, r1,5231,524Grassy, bitterKetone
41Propionic acidm, r1,5361,527SourAcid
42Acetoxy-2-butanonem, r1,546-SourKetone
432-Furfurylacetatem, r1,5521,566Fruity, greenEster
445-Methylfurfuralm, r1,5621,566Sweet, caramel, breadyFuran
453-Methylpyrrolem, r1,5741,569WoodyPyrrole
462-Acetylpyridinem, r1,5881,602Popcorn type, corn typePyridine
471-Methyl-2-formylpyrrolem, r1,5931,626Bread, burnt, caramelPyrrole
482-Furfuryl-5-methylfuranm, r1,6471,659ChocolateFuran
49g-Butyrolactonem, r1,658-Creamy, milkyEster
50Furfuryl alcoholm, r1,6711,662BurntAlcohol
51Isovaleric acidm, r1,6811,682CheesyAcid
522,5-Dihydrofuranonem, r1,7491,767CaramelKetone
531-Furfurylpyrrolem, r1,8141,822VegetablePyrrole
54Guaiacolm, r1,8431,855Phenolic, woodyPhenolic
552-Thiophenemethanolm,r1,8721,890savoryAlcohol
56Phenylethyl alcoholm, r1,9021,896FloralAlcohol
572-Acetylpyrrolem, r1,9531,949White breadPyrrole
58Difurfuryl etherm,r1,9911,977Coffee-like, toasted odourEther
59Phenolm, r2,011-SmokyPhenolic
604-Ethylguaiacolm,r2,0422,034PhenolicPhenolic
612-Formylpyrrolem, r2,0482,036White bread, jasmine ricePyrrole
62Furaneolm, r2,0532,037FruityKetone
63p-Cresolm, r2,0912,078Phenolic-typePhenolic
642-Methoxy-4-vinylguaiacolm, r2,112-Woody, smokyPhenolic

Detection of sixty-four volatile aroma compounds in the roasted coffee beans in the different agricultural treatments.

ID, Identification method; m, mass spectrum based on Mass library - NIST/EPA/NIH Mass Spectral Library. Version 2.0, Faircom Corporation, U.S.; r, identified by retention index; KI, Kovats index calculated by n-alkanes on the ZB-WAX Capillary GC Column (Length 30 m, inner diameter 0.25 mm, and film thickness 0.25 μm; Phenomenex Inc., Macclesfield, UK); Literature KI was collected from NIST database: webbook.nist.gov; Odor Description: data was collected from odor database: www.flavornet.org; origin-scifinder.cas.org).

Figure 1

For each contrast model, the fertilizer level with the relatively lower amount of volatile compounds generation was assumed as the control group and the other one with the relatively higher amount of volatile compounds generation was supposed to be the non-control group. The values of these volatile compounds in the control group were defined as 100. The ratio of the non-control group to the control group were calculated, and then, the final values of these volatile compounds in the non-control group were expressed as the ratio multiplied by 100. For example, in the Figure 2A, the MC with the relatively lower volatile compounds generation was assumed as the control group and showed as the circle with round markers at 100%, while the IC with the relatively higher abundance of volatile compounds was defined as the non-control group and showed as the spider diagram with diamond markers. Although sixty-four volatile compounds were detected in our experiment, only sixty volatile compounds were displayed in Figures 2, 3, as there were sno significant differences (P > 0.05) in the amounts of these four volatile compounds (2-Methylbutanal, Isovaleraldehyde, Phenol and Furaneol) between the different fertilizer levels and shade types.

Figure 2

Figure 3

Figure 2 showed the effect of different managements on the aroma generation of roasted coffee beans. The relative abundance of twenty volatile compounds were significantly higher (P < 0.05) in the roasted coffee beans under the IC management than in those roasted coffee beans under the MC management (Figure 2A, more details see Table 1 in supplementary). The amounts of twenty-five volatile compounds were significantly higher (P < 0.05) in the roasted coffee beans under the IO management than in those roasted coffee beans under the MC management (Figure 2B, more details see Table 2 in supplementary file). Nevertheless, only five volatile compounds had significantly higher (P < 0.05) amounts in the roasted coffee beans under the IO management than in those roasted coffee beans under the LO management (Figure 2C, more details see Table 3 in supplementary file). Moreover, there were eight volatile compounds with a significantly higher (P < 0.05) amount in the roasted coffee beans under the IC management than in those roasted coffee beans under the IO management (Figure 2D, more details see Table 4 in supplementary file).

Similarly, for each contrast model, the shade type with the relatively lower volatile compounds generation was assumed as the control group and the other one with the relatively higher volatile compounds generation was supposed to be the non-control group. The values of these volatile compounds in the control group were defined as 100. The ratio of the non-control group to the control group was calculated, and then, the final values of these volatile compounds in the non-control group were expressed as the ratio multiplied by 100.

Figure 3 showed the effect of different shade types on the aroma generation of roasted coffee beans. Thirty-six volatile compounds had a significantly higher (P < 0.05) amounts in the roasted coffee beans under the full sun system than in those roasted coffee beans under the shaded system (Figure 3A, more details see Table 5 in supplementary file). The amounts of forty volatile compounds were significantly higher (P < 0.05) in the roasted coffee beans under the full sun system than in those roasted coffee beans under the Erythrina system (Figure 3B, more details see Table 6 in supplementary file). Additionally, there were forty-three volatile compounds with a significantly higher (P < 0.05) amount in the roasted coffee beans under the service shade tree system than in those roasted coffee beans under the timber shade tree system (Figure 3C, more details see Table 7 in supplementary file). The amounts of fifty-one volatile compounds were significantly higher (P < 0.05) in the roasted coffee beans under the legume timber shade tree system than in those roasted coffee beans under the non-legume timber shade tree system (Figure 3D, more details see Table 8 in supplementary file).

As mentioned in the methods, L value ranges from 0 to 100, which represents the relevant color is from black (0) to white (100). As a result, the color turns darker with the decreasing of L value. Table 8 showed the effects of fertilizer levels and shade tree types on the color of roasted coffee beans.

Table 8

ContrastL-value
MeanS.E.Dp-value
Managements
IC vs. MC50.772.010.042
53.96
MC vs. IO54.022.210.006
50.26
IO vs. LO53.041.330.085
54.52
IC vs. IO53.162.210.102
54.03
Shade types
FS vs. S52.973.050.016
56.82
E vs. FS57.701.21<0.001
52.97
Ser.v. Tim.50.533.11<0.001
56.98
LT vs. NLT52.091.27<0.001
59.48

Contrast results for the color of roasted coffee beans.

Values were presented as mean, standard error of the contrast difference (S.E.D) and significance of the difference (P-value), P-values <0.05 were shown in bold. IC, Intensive conventional; MC, Moderate conventional; IO, Intensive organic; LO, Low Organic; FS, Full sun; S, Shaded; E, Erythrina poepiggiana; Ser., Service; Tim, Timber; LT, Legume timber; NLT, Non-legume timber.

With regards to the different fertilizer levels, L-values of roasted coffee beans were significantly lower (P < 0.05) under the IC (51) and the IO (50) management than under the MC (54) management. As a result, the color of roasted coffee beans was significantly darker (P < 0.05) under the IC and the IO management than under the MC management. Nevertheless, no significant differences (P > 0.05) were detected in the L values between the IO and LO, and the IC and IO.

In terms of shade types, the roasted coffee beans under the full sun system (53) had a significantly lower (P < 0.05) L-value than under the shaded (57) and Erythrina (58) systems. Thus, the color of roasted coffee beans was significantly darker (P < 0.05) under the full sun system than under the shaded and Erythrina systems. Furthermore, the L value of roasted coffee beans was significantly lower (P < 0.05) under the service shade tree system (51) than under the timber shade tree system (57), while there was a significantly higher (P < 0.05) L value detected in the roasted coffee beans under the non-legume timber shade tree system (59) than under the legume timber shade tree system (52). Therefore, the color of roasted coffee beans was significantly darker (P < 0.05) under the service shade tree system than under the timber shade tree system, whilst the roasted coffee beans under the non-legume timber shade tree system had a significantly lighter (P < 0.05) color than under the legume shade tree system.

Although many complex chemical reactions such as the degradation of sugar and acids occur during coffee bean roasting progress, Maillard reaction is obviously one of the most important chemical reactions to promote the formation of coffee flavor and color (2, 41). Maillard reaction is related to the interactions between reducing sugars and amino acids (42). Even though there is a large number of polysaccharides present in the green coffee beans (43), most carbohydrates are relatively stable to roasting (44, 45) and therefore it is the relative abundance of free or available sugars drives the thermal flavor generation. Moreover, on roasting, sucrose is rapidly degraded and its content dramatically decreases by about 97% (46). As a result, during coffee bean roasting, the reducing sugars and amino acids, as the main precursors of Maillard reaction, are mainly derived from the decomposition of sucrose and protein in the green coffee beans, respectively. Furthermore, the melanoidins, which contributes to the color formation of roasted coffee beans, are derived from Maillard reaction and dehydration of sucrose (47). Therefore, it is proposed that the contents of sucrose and total protein in the green coffee beans are responsible for the formation of coffee flavor and color.

There were no significant differences (P > 0.05) detected in the sucrose content of green coffee beans among the different fertilizer levels and shade types (see Table 6). Nonetheless, significant differences (P < 0.05) were found in the total protein content of green coffee beans among the different fertilizer levels and shade types (see Table 6).

In terms of the different fertilizer levels, the total protein content of green coffee beans was significantly lower (P < 0.05) under the MC management than under the IC and IO management. This explains why the roasted coffee beans under MC management presented a significantly lower (P < 0.05) number of volatile compounds and lighter color than those roasted coffee beans under the IC and IO management.

With regards to the different shade types, the amounts of volatile compounds in the roasted coffee beans were significantly higher (P < 0.05) and the color of roasted coffee beans was significantly darker (P < 0.05) under the full sun system than under the shaded and Erythrina system due to a significantly higher (P < 0.05) total protein content of green coffee beans in the full sun system. Similarly, due to the significantly higher (P < 0.05) total protein content of green coffee beans under the service and legume timber shade tree system, the roasted coffee beans under the service and legume timber shade tree system had a significantly higher (P < 0.05) number of volatile compounds and significantly darker (P < 0.05) color than those roasted coffee beans under the timber and non-legume timber shade tree system, respectively.

Conclusion

To sum up, in our experiment, the total protein content of green coffee beans was significantly affected by the fertilizer levels and shade types, while there was no significant effect on the sucrose content of green coffee beans. As a result, the total protein content was the main factor to drive the formation of aroma and color of roasted coffee beans. As has been discussed, the green coffee beans under the intensive conventional (IC) and intensive organic (IO) management had a significantly higher (P < 0.05) total protein content than those under the moderate conventional (MC) management due to the higher fertilizer input amount, especially the nitrogen. Thus, the higher N input amount could increase the protein generation in the green coffee beans and further improve the formation of final flavor and color of roasted coffee beans. Furthermore, the legume shade tree could improve the protein generation in the green coffee beans by fixing more N2 and further promote the formation of final flavor and color of roasted coffee beans. Compared with shaded system, the light intensity of full sun system was significant higher, which improved the protein biosynthesis in the ribosome, and finally, promoted the formation of flavor and color of roasted coffee beans.

Funding

This research was financially supported by Science and Technology Fund Project of Guizhou (Qian Ke He Zhicheng [2022] Zhongdian No. 015), Youth Science and Technology Talent Development Project of Education Department of Guizhou Province (Qian Jiao He KY Zi [2018]303), the Special Funding of Guiyang Science, and Technology Bureau and Guiyang University (GYU-KY-[2021]).

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

SX, IF, and EV conceived the experiments. SX, YL, and FM drafted the manuscript. SX conducted all the experiments. NY helped to discuss the results and perfected the language assisted with the structure elucidation and manuscript revision. IF and EV designed and supervised the research and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.968783/full#supplementary-material

References

  • 1.

    FiskIDKettleAHofmeisterSVirdieAKennyJS. Discrimination of roast and ground coffee aroma. Flavour. (2012) 1:14. 10.1186/2044-7248-1-14

  • 2.

    DartSNurstenH. Volatile Components Coffee. Dordrecht: Springer (1985).

  • 3.

    TresslR. Formation of flavor components in roasted coffee. ACS Publications. (1989) 9:27. 10.1021./bk-1989-0409.ch027

  • 4.

    HoCHwangHYuTZhangJ. An overview of the Maillard reactions related to aroma generation in coffee. COLLOQUE Scientifique International sur le Café. (1993) 15:611.

  • 5.

    HolscherWSteinhartH. Formation pathways for primary roasted coffee aroma compounds. ACS Publications. (1994) 3:16. 10.1021./bk-1994-0543.ch016

  • 6.

    HolscherWSteinhartH. Aroma compounds in green coffee. Develop Food Sci. (1995) 37:785803.

  • 7.

    ReinecciusG. The Maillard reaction and coffee flavor. 16th International Scientific Colloquium on Coffee, Kyoto (1995) 914.

  • 8.

    SemmelrochPGroschW. Studies on character impact odorants of coffee brews. J Agric Food Chem. (1996) 44:53743.

  • 9.

    AgrestiPDMFrancaASOliveiraLSAugustiR. Discrimination between defective and non-defective Brazilian coffee beans by their volatile profile. Food Chem. (2008) 106:78796.10.1016/j.foodchem.06019

  • 10.

    SilvaEADMazzaferaPBruniniOSakaiEArrudaFBMattosoLHCet al. The influence of water management and environmental conditions on the chemical composition and beverage quality of coffee beans. Br J Plant Physiol. (2005) 17:22938. 10.1590/S1677-04202005000200006

  • 11.

    RisticevicSCarasekEPawliszynJ. Headspace solid-phase microextraction–gas chromatographic–time-of-flight mass spectrometric methodology for geographical origin verification of coffee. Analytica Chimica Acta. (2008) 617:7284. 10.1016/j.aca.2008.04.009

  • 12.

    FrancaASOliveiraLSOliveiraRCAgrestiPCMAugustiR. A preliminary evaluation of the effect of processing temperature on coffee roasting degree assessment. J Food Enginee. (2009) 92:34552. 10.1016/j.jfoodeng.2008.12.012

  • 13.

    BhumiratanaNAdhikariKChambersE. Evolution of sensory aroma attributes from coffee beans to brewed coffee. LWT-Food Sci Technol. (2011) 44:218592. 10.1016/j.lwt.2011.07.001

  • 14.

    BottazziDFarinaSMilaniMMontorsiL. A numerical approach for the analysis of the coffee roasting process. J Food Enginee. (2012) 112:24352. 10.1016/j.jfoodeng.2012.04.009

  • 15.

    YeretzianCJordanALindingerW. Analysing the headspace of coffee by proton-transfer-reaction mass-spectrometry. Int J Mass Spectromet. (2003) 223:11539. 10.1016/S1387-3806(02)00785-6

  • 16.

    FreitasACMoscaA. Coffee geographic origin—an aid to coffee differentiation. Food Res Int. (1999) 32:56573.

  • 17.

    AeschbacherHWollebULöligerJSpadoneJLiardonR. Contribution of coffee aroma constituents to the mutagenicity of coffee. Food Chem Toxicol. (1989) 27:22732. 10.1016/0278-6915(89)90160-9

  • 18.

    Bade-WegnerHHolscherWVitzthumO. Quantification of 2-methylisoborneol in roasted coffee by, GCMSCOLLOQUE. Scientifique International Sur Le, Cafe Asic Association Scientifique Internationale. (1993). 56:537.

  • 19.

    Amrani-HemaimiMCernyCFayLB. Mechanisms of formation of alkylpyrazines in the Maillard reaction. J Agri Food Chem. (1995) 43:281822. 10.1021/jf00059a009

  • 20.

    BlankIDevaudSFayLB. Study on the formation and decomposition of the amadori compound. Maillard Rea Foods Med. (1998) 223:43.

  • 21.

    CzernyMGroschW. Potent odorants of raw Arabica coffee. Their changes during roasting. J Agricul Food Chem. (2000) 48:86872.

  • 22.

    HaggarJBarriosMBolañosMMerloMMoragaPMunguiaRet al. Coffee agroecosystem performance under full sun, shade, conventional, and organic management regimes in Central America. Agrofor Sys. (2011) 82:285301. 10.1007/s10457-011-9392-5

  • 23.

    NoponenMREdwards-JonesGHaggarJPSotoGAttarzadehNHealey. Greenhouse gas emissions in coffee grown with differing input levels under conventional and organic management. Agri Ecosys Environ. (2012) 151:615. 01, 019.

  • 24.

    MuschlerRG. Shade improves coffee quality in a sub-optimal coffee-zone of Costa Rica. Agrofor Sys. (2001) 51:1319.

  • 25.

    BaggenstossJPoissonLLuethiRPerrenREscherF. Influence of water quench cooling on degassing and aroma stability of roasted coffee. J Agri Food Chem. (2007) 55:668591. 10.1021/JF070338D

  • 26.

    FiskIDGkatzionisKLadMDoddCERGrayDA. Gamma-irradiation as a method of microbiological control, and its impact on the oxidative labile lipid component of Cannabis sativa and Helianthus annus. Eu Food Res Technol. (2009) 228:61321. 10.1007/s00217-008-0970-3

  • 27.

    KyCLDoulbeauSGuyotBAkaffouSCharrierAHamonSet al. Inheritance of coffee bean sucrose content in the interspecific cross Coffea pseudozanguebariaex Coffea liberica 'dewevrei'. Plant Breeding. (2000) 119:1658. 10.1046/j.1439-0523.2000.00464.x

  • 28.

    MccuskerRRGoldbergerBAConeEJ. Caffeine content of specialty coffees. J Anal Toxicol. (2003) 27, 520522. 10.1093/jat/27.7.520

  • 29.

    PerroneDDonangeloCMFarahA. Fast simultaneous analysis of caffeine, trigonelline, nicotinic acid and sucrose in coffee by liquid chromatography–mass spectrometry. Food Chem. (2008) 110:10305. 10.1016/j.foodchem.2008.03.012

  • 30.

    YangNLiuCLiuXDegnTKMunchowMFiskI. Determination of volatile marker compounds of common coffee roast defects. Food Chem. (2016) 211:20614. 10.1016/j.foodchem.2016.04.124

  • 31.

    TresslRWondrakGTGarbeLAKrügerRPRewickiD. Pentoses and hexoses as sources of new melanoidin-like Maillard polymers. J Agri Food Chem. (1998) 46:17656.

  • 32.

    GinzMBalzerHHBradburyAGMaierHG. Formation of aliphatic acids by carbohydrate degradation during roasting of coffee. Eu Food Res Technol. (2000) 211:40410.

  • 33.

    NygrenPFernándezMPHarmandJMLeblancHA. Symbiotic dinitrogen fixation by trees: an underestimated resource in agroforestry systems?Nutrient Cycl Agroecosys. (2012) 94:123160. 10.1007/s10705-012-9542-9

  • 34.

    BoreuxVVaastPMadappaLPCheppudiraKGGarciaCGhazoulJ. Agroforestry coffee production increased by native shade trees, irrigation, and liming. Agronomy Sustain Develop. Paris (2016) 36:19. 10.1007/s13593-016-0377-7

  • 35.

    BoulterD. Protein synthesis in plants. Annu Rev Plant Physiol. (1970) 21:91114.

  • 36.

    KuceyRPaulEA. Carbon flow in plant microbial associations. Science. (1981) 6:473. 10.1126./science.213.4506.473

  • 37.

    StumpfWConnPMPreissJ. The Biochemistry of Plants: Carbohydrates: Academic Press (2012). 10.1016/C2009-0-02820-7

  • 38.

    BoardmanNT. Comparative photosynthesis of sun and shade plants. Annu Rev Plant Physiol. (1977) 28:35577.

  • 39.

    RobartsRDZoharyT. Temperature effects on photosynthetic capacity, respiration, and growth rates of bloom-forming cyanobacteria. N Z J Mar Freshwater Res. (1987) 21:3919.

  • 40.

    LiuYangNYangQAyedCLinforthRFiskIDet al. Enhancing Robusta coffee aroma by modifying flavour precursors in the green coffee bean. Food Chem. (2019). 281:817.

  • 41.

    GroschW. (2001). Chemistry III: volatile compounds. Coffee: recent Develop6889. 10.1002./9780470690499.ch3

  • 42.

    MaillardLC. Action des acidesamines sur les sucres: formation des melanoidines par voie methodique. CR Acad Sci(Paris).Paris (1912) 154:668.

  • 43.

    BradburyAGWHallidayDJ. Polysaccharides in Green Coffee Beans. Paris: ASIC (1987).

  • 44.

    BradburyAGWHallidayDJ. Chemical structures of green coffee bean polysaccharides. J Agric Food Chem. (1990) 38:38992.

  • 45.

    WilsonAPetraccoMIllyE. Some preliminary investigations of oil biosynthesis in the coffee fruit and its subsequent re-distribution within green and roasted beans. Proceedings of the 17th International Conference on Coffee, Science. (1997) 92.

  • 46.

    TrugoLCMacraeR. The Use of the Mass Detector for Sugar Analysis of Coffee Products. Paris, France: ASIC (1985).

  • 47.

    Bradbury AGW. LA Carbohydrates. In: Clarke RJ, Vitzthum OG (eds.) Coffee: Recent Development. Cornwall: Blackwell Science Ltd (2001).

Summary

Keywords

agro-forestry systems, shade trees, organic management, coffee color, aroma generation

Citation

Xu S, Liu Y, Ma F, Yang N, Virginio Filho EM and Fisk ID (2022) Impact of agro-forestry systems on the aroma generation of coffee beans. Front. Nutr. 9:968783. doi: 10.3389/fnut.2022.968783

Received

14 June 2022

Accepted

12 July 2022

Published

04 August 2022

Volume

9 - 2022

Edited by

Mingquan Huang, Beijing Technology and Business University, China

Reviewed by

Ye Liu, Beijing Technology and Business University, China; Meigui Huang, Nanjing Forestry University, China

Updates

Copyright

*Correspondence: Su Xu

This article was submitted to Food Chemistry, a section of the journal Frontiers in Nutrition

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics