Next Article in Journal
Extensive Metabolic Profiles of Leaves and Stems from the Medicinal Plant Dendrobium officinale Kimura et Migo
Previous Article in Journal
Reactive Carbonyl Species as Potential Pro-Oxidant Factors Involved in Lichen Planus Pathogenesis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 9
Issue 10
10.3390/metabo9100214
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolomic Markers for the Early Selection of Coffea canephora Plants with Desirable Cup Quality Traits

by
Roberto Gamboa-Becerra
1,2,*,
María Cecilia Hernández-Hernández
3,
Óscar González-Ríos
3,
Mirna L. Suárez-Quiroz
3,
Eligio Gálvez-Ponce
4,
José Juan Ordaz-Ortiz
5 and
Robert Winkler
1,*
1
Department of Biotechnology and Biochemistry, Cinvestav Unidad Irapuato, Irapuato, Km 9.6 Libramiento Norte Carr. Irapuato-León, Guanajuato 36824, Mexico
2
Red de Biodiversidad y Sistemática, Instituto de Ecología A.C. Xalapa, Veracruz 91070, Mexico
3
Laboratorio de Tecnología del Café, Unidad de Investigación y Desarrollo en Alimentos, TNM/Instituto Tecnológico de Veracruz, Veracruz 91897, Mexico
4
Agroindustrias Unidas de México S.A. de C.V. (AMSA), Bosque de Alisos 45-A 2do Piso, Bosques de las Lomas, Ciudad de México 05120, Mexico
5
Laboratorio de Metabolómica y Espectrometría de Masas, Unidad de Genómica Avanzada, CINVESTAV-IPN, Km. 9.6 Libramiento Norte Carr. Irapuato-León, Irapuato 36824, Mexico
*
Authors to whom correspondence should be addressed.
Metabolites 2019, 9(10), 214; https://0-doi-org.brum.beds.ac.uk/10.3390/metabo9100214
Submission received: 23 August 2019 / Revised: 19 September 2019 / Accepted: 20 September 2019 / Published: 4 October 2019
Browse Figures
Review Reports Versions Notes

Abstract

:
Genetic improvement of coffee plants represents a great challenge for breeders. Conventional breeding takes a too long time for responding timely to market demands, climatic variations and new biological threads. The correlation of genetic markers with the plant phenotype and final product quality is usually poor. Additionally, the creation and use of genetically modified organisms (GMOs) are often legally restricted and rejected by customers that demand natural products. Therefore, we developed a non-targeted metabolomics approach to accelerate conventional breeding. Our main idea was to identify highly heritable metabolites in Coffea canephora seedlings, which are linked to coffee cup quality. We employed a maternal half-sibs approach to estimate the metabolites heritability in open-pollinated plants in both leaves and fruits at an early plant development stage. We evaluated the cup quality of roasted beans and correlated highly heritable metabolites with sensory quality traits of the coffee beverage. Our results provide new insights about the heritability of metabolites of C. canephora plants. Furthermore, we found strong correlations between highly heritable metabolites and sensory traits of coffee beverage. We revealed metabolites that serve as predictive metabolite markers at an early development stage of coffee plants. Informed decisions can be made on plants of six months old, compared to 3.5 to 5 years using conventional selection methods. The metabolome-wide association study (MWAS) drastically accelerates the selection of C. canephora plants with desirable characteristics and represents a novel approach for the focused breeding of crops.

Graphical Abstract

1. Introduction

Coffee is one of the most important crops in agriculture around the world since coffee beverages prepared from roasted and ground beans are the second most popular beverage only after tea [1,2]. More than 100 species have been reported belonging to Coffea genus, but Coffea arabica and Coffea canephora are the most commercialized. Because of their economic importance, these two species have been the main aim of the genetic breeding programs of coffee [3,4]. Traditional genetic breeding strategies for coffee mainly include traditional clones selection, generational selection, selection of pure lines following generational selection, generation of F1 hybrids, interspecific and intraspecific hybridization and selection by backcrossing. The best strategy is chosen according to the behavior of the coffee species, the aim of breeding (disease resistance, quality, yield) and the kind of variety (clone, pure line or variety, hybrid) [5,6,7]. Coffee is a perennial plant, and the first harvest is possible after 3.5 to 5 years. Thus, the selection and breeding of new coffee varieties require at least 25 years, with vast investments of research and development resources.
Last years, priority was given to the improvement of the cup quality of C. canephora [8]. This species is considered to have a lower cup quality compared to C. arabica. On the other side, C. canephora is more resistant against diseases and more tolerant against climatic variations [9,10,11]. Furthermore, the content of caffeine and soluble solids extraction yields are significantly higher, which makes C. canephora the principal raw material in the production of instant coffee [4,12,13]. Bean size, soluble solids content, cup quality and cup quality-related compounds (the content of sugars, lipids, caffeine, trigonelline and chlorogenic acids) has been the aim of breeding programs in C. canephora. Despite the economic importance of coffee growing and breeding, technological progress is limited compared to other crops [14,15,16,17,18,19]. Therefore, new methodologies for the selection of coffee varieties with desirable agronomical and cup quality characteristics are urgently needed.
Plant phenotypes are tightly related to metabolites. Consequently, the employment of metabolites as predictive markers for complex traits promises to be of high potential. Metabolomics is a powerful tool for the generation of targeted and non-targeted metabolic profiles. The use of advanced statistical analysis provides new possibilities for discovering metabolites, which are predictive for phenotypical characteristics [20,21]. In several crops, such as barley, maize, pine, potato, rice, tomato, wheat and grape [22,23,24,25,26,27,28,29], metabolomics approaches permitted the identification of metabolic markers as accurate predictors of key agronomical traits. The chemical composition of food and beverages also has a direct relation with sensory traits and product quality, such as its nutritional value. Therefore, the raw material used in food productions has a crucial influence on the satisfaction of the consumer. For this reason, we developed new breeding strategies to enable the selection of plant varieties with desirable agronomic and sensory traits.
Metabolome-wide association studies (MWAS) represent a potential tool for the improvement of plants, regardless of the availability of genomic information. Therefore, we applied a metabolomics strategy to identify metabolites predicting sensory traits of coffee beverage in the early development stages of C. canephora plants.

2. Results and Discussion

2.1. Variability of Metabolic Phenotypes and the Relationship Between C. Canephora Parent Plants

Forty families of C. canephora were included in this study, comprising a total of 120 coffee plants. Three individuals integrated every family: one parental (mother plant) and two mother half-sibs for every parental (Table 1). The identity of one of the parental was unknown since the coffee plants sibs were reproduced by open-pollinated seeds. To evaluate the natural variation in metabolite levels present in plants of C. canephora families, we performed an untargeted metabolomics approach, Ultra Performance Liquid Chromatography-Ion Trap-Mass Spectrometry (UPLC-IT-MS) of aqueous methanolic extracts of fruit and leaf tissues. First, we conducted a hierarchical clustering analysis (HCA) to determine the metabolic variability and relationship between coffee mother plants. We generated two dendrograms depicting the clustering of samples employing the abundances of 755 metabolites detected in fruits and 527 metabolites for leaf tissue. Then, we generated a tanglegram plot in order to compare both phylogenies (Figure 1). Our analyses of phylogeny congruence revealed a high degree (0.81) of entanglement, and comparisons indicated that this was mainly due to the change of position of regions or entire clades in the tree. Naturally, fruit and leaf tissues show very different metabolic profiles. The co-clustering of genetically-related mother plants indicates their metabolic similarity and supports the strong correlation between genotype and metabolic phenotype.

2.2. Selection of Highly Heritable Metabolites in Coffee Fruits and Leaves of C. Canephora Families

In coffee fruits, we detected 755 metabolites, which were present in all 120 samples, whereas a minor quantity of metabolites was detected in leaves (527 metabolites). Metabolites in plants are tightly associated with phenotypic characteristics and represent a high diversity at a quantitative and qualitative level within plant species. The application of untargeted metabolomics analysis enables the detection of a wide range of metabolite classes in plants and support the study of aspects regulating biochemical variation such as environmental and genetic control and how these compounds influence complex traits [20,30,31].
To evaluate the variability in metabolite levels and determine the contribution of the genetic factor, we determined the abundance of every metabolite and estimated the metabolic heritability in the narrow sense by offspring-parents regression. Heritability levels for traits range from 0 to 1. Scores less than 0.2 are categorized as low heritability, scores between 0.2 and 0.39 are considered as moderate heritability and high heritability corresponds to values above 0.4. In fruits, our results revealed a high percentage of metabolites (mass peaks) with high heritability (59.6%, 450 metabolites), while 97 and 208 metabolites exhibited low and moderate heritability levels, respectively (Figure 2A). In contrast, for metabolites heritability estimated in leaf tissues, we detected 22.01% (116 metabolites) with high heritability values, 125 metabolites of moderate heritability (23.72%) and 286 metabolites (54.27%) were classified as low heritability metabolites (Figure 2B). Highly heritable metabolites identified from mother plants to F1 progeny: 450 in fruits and 116 in leaves were considered for the posterior analysis of correlation with sensory traits of coffee beverage. Our results demonstrated that for a considerable number of metabolites, a high percentage of the observed quantitative variation could be explained by genetic factors, as demonstrated by their high heritability levels. A study of metabolic diversity by using metabolomics approaches permits to evaluate natural variation occurring within and between species [32]. Our findings are relevant for the knowledge of the variability naturally occurring in the levels of metabolites of coffee tissues and the contribution of genetic factors. Studies of heritability of qualitative and quantitative characters of the coffee crop have been the basis for the genetic improvement of quantitative traits such as yielding, disease resistance, adaptability to climatic and soil conditions, bean physical quality and sensory traits of coffee beverage [19,33,34,35,36,37]. Our highlights, together with association studies between metabolites and coffee traits, could be of great importance to support the classical breeding of coffee plants with desirable traits, especially those that are associated with the natural variation of metabolite composition and levels.

2.3. Metabolites of Roasted Coffee Beans Linked to Sensory Traits of Coffee Beverage

2.3.1. Tasting and Correlation Between Sensory Attributes of Coffee Beverage

For investigating the relationships between sensory traits of coffee beverage, a panel of 11 trained cup-testers evaluated 21 roasted coffee bean samples. Eleven samples randomly selected corresponded to C. canephora mother plants (M480–M495) and eight C. arabica and C. canephora samples were used as reference samples. Ten sensory attributes of coffee beverage were tested, including aromatic quality, aromatic intensity, astringency, body, flavor, after taste, bitterness, acidity, sourness and global preference, using a scale ranging of 0 to 5 points (Table 2). First, we performed a correlation analysis between sensory attributes (Figure 3). Based on a Spearman correlation analysis, the aromatic intensity attribute only showed a positive correlation with aromatic quality, with a Spearman correlation coefficient (SCC) of 0.48 and statistical significance (p < 0.05). Aromatic quality was positively correlated with eight out of nine attributes with the exception of acidity (SCC: −0.55; p < 0.01) and sourness (SCC: −0.62; p < 0.01) for which negative correlation was observed. We obtained a similar correlation pattern for flavor, astringency and aftertaste. The latter demonstrated positive correlation with statistical significance for most of the attributes except for acidity and sourness. Interestingly, for acidity, we detected negative correlation for six attributes (aromatic quality, flavor, aftertaste, bitterness, body and global preference) and only for the sourness attribute, we found a high Spearman correlation coefficient of 0.87 (p < 0.001). Global preference represents the rating of the overall perception of all the organoleptic characters of coffee beverage. This attribute revealed positive correlation for aromatic quality, flavor, astringency, aftertaste, bitterness and body. However, we found a high negative correlation for acidity and sourness attributes (SCC: −0.85 and −0.8; p < 0.001, respectively). The close relationship between acidity and sourness and the correlation between all the sensory traits assessed is illustrated in a dendrogram of correlation distances (Figure 4). Altogether, our results demonstrated the existence of relationships between the majority of the sensory traits, and suggest that attributes do not only affect coffee quality individually.

2.3.2. Selection of Marker Metabolites in Roasted Coffee Beans Linked to Sensory Traits

Coffee cup quality is determined by diverse factors, such as the coffee plant genotype, environmental factors (geographical origin, climate, altitude, temperature, percentage of crop shade, soil nutrient capacity, fertilization), agricultural practices, post-harvest processing and grain storage, as well as the basic chemical composition of coffee bean. A wide range of compounds contributes to the coffee cup quality, including volatile and non-volatile compounds, which determine the organoleptic quality of coffee beverage [38,39,40,41,42].
Non-volatile compounds of coffee roasted beans include carbohydrates, sucrose, free sugars, cell wall polysaccharides, lignin, nitrogen compounds (proteins, free amino acids, caffeine, trigonelline), lipids, minerals (potassium, phosphorus), organic acids, esters (chlorogenic acids, aliphatic acids, quinic acid derivates) and melanoidins [19,35]. To identify metabolite markers of coffee cup quality present in roasted coffee beans, we analyzed aqueous extracts of roasted and ground beans by UPLC-Electrospray (ESI)-MS. We detected 737 metabolites that were present in all 21 samples of C. canephora and C. arabica samples included for this analysis. We performed a Spearman correlation analysis employing the levels of detected metabolites versus the scores of each of the 10 sensory traits of coffee beverage before mentioned. We detected 392 correlations that fulfilled those characteristics, involving 174 metabolites and eight sensory attributes of coffee beverage, with SCC values ranging from 0.5 to 0.813 for positive correlations and −0.5 to −0.849 for negative correlations (see Table S1). The results indicated no metabolites (with statistical significance) linked to aromatic intensity and astringency attributes, whereas for all the other sensory attributes, at least 21 correlations were found. To identify the metabolites, we compared high-resolution accurate mass and MS data by using the programs XCMS, MSDIAL and SpiderMass. We putatively identified 91 (52.3%) of the metabolites (Table S1), which were classified in 11 class compounds, comprising alkaloids, carbohydrates, carotenoids, chlorogenic acids, fatty acids, flavonoids, hydrocarbons, lipids, organic acids and terpenoids (Figure 5). Several metabolites of these class compounds have been identified as biochemical compounds playing a role in roasting chemistry. For example, carbohydrates and free sugars suffer pyrolysis, thermal and Maillard (with proteins and amino acids) reactions generate aroma precursors such as furanes, carboxylic acids, pyridines, pyrroles, pyrazines, thiazoles and ketones. Ferulic, caffeic and quinic acids from chlorogenic acids were degraded to produce phenols, and lipids through oxidation were transformed to alcohols, ketones and aldehydes [2,43,44,45,46].
Interestingly, we detected a high number of lipids influencing the sensory attributes of coffee beverage. Fifty-two metabolites (114 associations) were annotated as lipids comprising four classes: glycerolipids (GL), glycerophospholipids (GP), prenol lipids (PR), sphingolipids (SP) and sterol lipids (ST) (Figure S1). Lipids of roasted coffee beans enhance the organoleptic properties of coffee beverage such as flavor, color and foam. Lipid composition of beans is made mainly of triacylglycerides (75%) and terpene esters (14%), free fatty acids (1%), free sterols (1.5%), sterol esters (1%) and phospholipids (PLs) which are present in lower levels [47]. We detected a phospholipids group containing SP and GP classes. SPs were associated with acidity, aftertaste, bitterness, global preference and sourness, whereas GPs were negatively or positively linked to all the evaluated attributes (Figure S1). GLs were mainly associated with acidity and body attributes, for which the majority of lipids were negatively associated (SCC = −0.50 to −0.584). Although there is much information about compounds present in coffee roasted beans or coffee beverage, limited information exists about the individual or specific chemical components that are directly linked to organoleptic characteristics of coffee beverage. Identification of metabolites linked to several sensory traits is crucial because of their value as marker metabolites for specific coffee sensory traits. In our study, we successfully identified some metabolites that were negatively or positively correlated to eight out of 10 assessed sensory attributes.
Figure 6 shows association levels of sensory traits with metabolites putatively identified as quercetin 4’-glucoside and trigonelline representing two examples of marker metabolites presented as critical components which confer specific coffee sensory properties. Quercetin 4´-glucoside is an organic compound that belongs to flavonoid o-glycosides, which contain a carbohydrate moiety linked to the 2-phenylchromen-4-one flavonoid backbone by O-glycosidic bond. Quercetin 4´glucoside was found to have a positive correlation with acidity and sourness (SCC: 0.56 and 0.64; p:0.0089 and 0.0017, respectively), whereas a negative correlation was observed for bitterness, flavor, aromatic quality, aftertaste, body and global preference attributes (Figure 6). Some glycosides have been determined as marker metabolites influencing the sensory profile of the coffee beverage. Mozambioside was reported as an Arabica-specific bitter-tasting furokauran glucoside detected in roasted coffee beans and beverage [48], whereas 3-methylbutanoyl glycosides were identified as markers for the quality of coffee flavors [49]. Trigonelline was linked to seven sensory traits, but a positive correlation was observed only for acidity, whereas for aromatic quality, bitterness, body, aftertaste, flavor and global preference, we detected a negative association. It is well-known that 50–80% of trigonelline is degraded by pyrolysis after the roasting process and contributes to the formation of other compounds such as nicotinic acid, pyridine, 3-methyl-pyridine, the methyl ester of nicotinic acid and 1-methylpyridinium [50,51]. All these products of trigonelline degradation finally influence organoleptic features of coffee beverage. Thus, the trigonelline content in roasted coffee beans represents a predictive marker for more than one sensory trait.

2.4. Highly Heritable Metabolites in Fruits Linked to Coffee Cup Quality

Some aspects, such as ploidy incompatibility between C. arabica and C. canephora, as well as the time that elapses from seed germination until the first harvest (3.5 to 5 years), represent a disadvantage for traditional breeding methods. For these reasons, obtaining an improved coffee variety might take up to 30 years [52,53]. Responding to the aforementioned, we performed a correlation analysis, which included the 450 metabolites mentioned, identified as highly heritable metabolites and scores of 10 sensory traits of coffee beverage. We detected 228 correlations involving 158 metabolites that were associated to at least one sensory trait and fulfilled the following characteristics: values of SCC > 0.5 for positive correlation and < −0.5 for negative correlation, and statistical significance (p < 0.05). 141 associations (61.8%) were negative correlations (SCC: −0.666 to −0.96) whereas 87 associations were detected as positive correlations (SCC: 0.66 to 0.945) (See Table S2). Aromatic quality, aromatic intensity and bitterness displayed the highest number of associations to highly heritable metabolites, 39, 29 and 29, respectively. On average, we detected 1.5 significant associations per metabolite feature. In Table 3, we list the top metabolites of coffee fruits with optimal characteristics, which explicated high predictive value for sensory traits of coffee beverage and for selection of plants producing coffee beans with desirable sensory characteristics. In the top metabolites list, we included 2,4-decadienoyl-CoA (m/z 940.2271), a metabolite constituted by unsaturated fatty acids with trans double bonds and which is degraded into acetyl-CoA by a beta-oxidation. High levels of this fatty acyl thioester of coffee fruits were positively associated with aftertaste (heritability: 0.881; SCC: 0.754).
The mass peak m/z 791.4983, putatively identified as nonaprenyl diphosphate, demonstrated a high heritability level (0.984) and positive correlation with global preference (SCC = 0.728). We identified two compounds that belong to polyphenols/flavonoids class: 2,3-trans-proanthocyanidin with high predictive value for astringency (heritability: 0.957; SCC: 0.833) and delphinidin 3-O-glucosyl-5-O-caffeoylglucoside which presented a positive correlation with body sensory trait (heritability: 0.690; SCC: 0.833). Six unidentified mass peaks (m/z 963.4134, m/z 860.5079, m/z 536.0019, m/z 734.4907, m/z 598.5406, m/z 521.5618) were associated to acidity, aromatic intensity, aromatic quality, bitterness, flavor and sourness, respectively. Improvement of organoleptic traits in coffee requires knowledge about the chemical composition of tissues of coffee plants, which contribute to the final product: roasted beans and coffee beverage [53]. The presented markers assist in the improvement and manipulation of sensory traits of coffee beverages. In the past, quality improvement of C. canephora has been performed for bean size, extractable soluble solids, sugars, caffeine, trigonelline, total lipids content and chlorogenic acids [8]. Here, we reported the first study that combined analysis of high throughput metabolite heritability in combination with a correlation analysis of sensory traits of coffee beverage. In other plant species, such as Arabidopsis [30], maize [28], tomato [54,55,56] and raspberry [57], the application of metabolomics together with complementary tools, such as genome-wide selection, has permitted the study of metabolite heritability, the elucidation of factors controlling biochemical variability and the purpose for modification of chemical composition by classical breeding. We reported a considerable number of metabolites linked to sensory traits, which supports the early selection of C. canephora plants.

2.5. Highly Heritable Metabolites in Leaf Linked to Coffee Cup Quality

To evaluate the link between highly heritable metabolites that have been detected before in coffee leaves and the values of 10 sensory traits of coffee beverage assessed by the trained cup testers, we performed a correlation analysis. We found 75 associations between 45 metabolites, which were linked at least with one sensory attribute of coffee beverage (Table S2). The average of significant associations per mass peak was 1.7, slightly higher than values found for fruits (1.5). Forty-six associations, corresponding to 61.3%, were negative correlations (SCC: −0.663 to −0.9746), and the remaining 38.7% of associations were positive correlations with values of SCC between 0.666 and 0.866) (Table S2). Mass peaks m/z 251.4295, m/z 387.3897, m/z 473.6946, m/z 517.7342, m/z 587.1469, m/z 695.3553 and m/z 805.7835 are interesting because all presented that they were associated at least with three coffee sensory traits of coffee beverage. This fact is of great importance, because we can evaluate multiple sensory traits simultaneously by screening the levels of a single metabolite. In Table 3, we listed the top metabolites purposed as excellent candidates for early selection of plants with desirable coffee sensory traits due to their outstanding characteristics of high heritability and strong positive correlation level to cup quality traits. 6-deoxoteasterone (m/z 473.3339), involved in brassinosteroid biosynthesis I in coffee plants, was linked to three sensory traits: acidity, sourness and bitterness. For the last one, a negative correlation (SCC = −0.857) was detected. It was included in the list, because for bitterness, only negative associations were observed. We found a similar situation for the ‘body’ attribute: 4,4-dimethylzymosterol (m/z 463.3451) was negatively linked to this trait. Two lipids, PA(16:0/18:2) (a glycerophosphate: GP) and MGDG(23:0/26:0) (a glycosyldiradylglycerol: GL05) were found in coffee leaves and presented high predictive values for aftertaste and astringency. Levels of mirystic acid, a fatty acid (FA) resulted as a suitable predictor for aromatic quality and flavor (heritability: 0.994; SCC: 0.766 and 0.672, respectively). Global preference was strongly associated with levels of gibberellin A28, with heritability levels of 0.996 and an SCC of 0.677. A level distribution analysis of all that metabolites in the 40 C. canephora parent plants (Figure 7) highlights individuals which can be subject for the breeding of plants with specific sensory characteristics. Therefore, we provided new knowledge about the inheritance of coffee leaf metabolites, but also information about the link between metabolite levels in the leaf and complex agronomic traits such as coffee cup quality. Remarkably, the use of highly heritable metabolite markers for coffee breeding in 6-month-old plants, instead of waiting for 3.5–5 years until the first coffee fruit harvest, opens perspectives to reduce the time for generation of coffee varieties with desirable sensory traits. Nevertheless, it is essential to emphasize that the chemical composition and sensory properties of coffee are not only determined by genetic factors, but also by cultivation conditions, plant condition, cultivation practices, harvest and post-harvest processes [19,39,58,59]. Therefore, a strategy for improving traits of interest must integrate all those non-genetic factors as well. Selecting plants with the proposed strategy would comprise: a) screening of a parent plant population using marker metabolites, b) selection of plants with desirable characteristics based in metabolite levels and c) plant propagation.

3. Materials and Methods

3.1. Heritability Analysis

3.1.1. Leaves and Fruits Sampling

Samples for heritability analysis were collected from an experimental farm in Tapachula, Chiapas, Mexico. We designed a maternal half-sibs approach in order to estimate the heritability of metabolites in both tissues, leaves and fruits. In this study, 120 C. canephora plants were included, consisting of 40 mother plants and their F1 progeny (two sibs for every mother plant), which was reproduced by open-pollinated seeds. Leaves of all plants (mothers and sibs) were sampled at 6 months of age for F1 plants and, three years later in a second sampling (at first fruit harvest of F1 plants), completely mature fruits of both mother plants and sibs were obtained.

3.1.2. Leaves and Fruits Tissue Metabolite Extraction

Approximately two grams of freeze-dried tissue were ground in a Mixer Mill MM 301 (Retsch®) at 30 Hz for 15 s and sieved with a mesh for <300 µm particle size. Then, 35 mg of every sample was extracted twice in 450 µL of solvent (75% methanol/24.85% water, 0.15% formic acid), 25 min in an ultrasound bath at room temperature. After centrifugation at 15,870× g for 5 min at room temperature, the supernatants of first and second extraction were collected in the same 1.5 mL tube and filtered using a 0.2 µm membrane. These fractions were subjected to UPLC-MS analysis.

3.1.3. Estimation of Metabolites Heritability

Narrow sense heritability of metabolites was calculated through offspring-parents regression, which is suitable for cross-fertilizing crops, and when the estimation of heritability is performed for half-sib families regressed on a single parent. Measures corresponded to intensity levels of every metabolite detected. The regression coefficient (Yi), which measures the rate between covariance of parents to offspring covariance, with its respective standard error was calculated as follows: Yi = β0 + β1Xi + εi, where: Yi = Average measures of offspring of the X families; β0 = intercept; β1 = regression coefficient; Xi = Measurement of only one parent (mother plant) per family; εi = Standard error. The slope of the regression line (β1) represents an estimate of narrow-sense heritability (h2), which depends mainly on additive genetic values. Narrow-sense heritability (h2) is very important for plant breeding programs because artificial selection mainly depends on additive genetic variance.

3.2. UPLC-ESI-MS Analysis

Metabolite analysis was performed employing an UPLC-ESI-MS (Ultra Performance Liquid Chromatography Electrospray Mass Spectrometry) system consisting of an Accela UPLC coupled to LCQ Fleet Ion Trap (Thermo Finnigan, San José, CA, USA). Separation of compounds was performed on a Hypersil Gold C18 column (50 × 2.1 mm, 1.9 µm particle size). In the UPLC system, mobile phase A consisted of 0.1% FA in water and solvent B was MetOH with 0.1% FA. For leaves analysis, the solvent gradient was as follows: 8% B, 0–0.5 min; 8–65% B, 0.5–6.5 min; 65%–100% B, 6.5–20 min; 100% B during 1 min (20–21 min); 0.5 min for comeback to re-equilibration conditions and finally, 8% B during 3.5 min (21.5–25 min). The gradient standardized for fruit tissues was: 3% B during 1 min; 3–50%, 1–5 min; 50–73% B, 5–5.5 min; 73–100% B, 5.5–19.3 min; 100% B during 1min; 100% B–3% B, 20.3–21.5; and 3.5 min 3% B for column re-equilibration. For roasted coffee analysis, the solvent gradient was: 0–1 min, 2% B; 1–6 min, 2–66% B; 6–12 min, 66–88% B; 12–19 min, 88–93% B; 93–97% B in 1 min; 97–100% B, 20–25 min; from 25 to 27 min was maintained at 100% B; 27–27.5 min, 100% B–98% A and finally, column re-equilibration for 3.5 min with 2% B. The column oven temperature was always maintained at 40 °C, and the flow rate was 550 µL/min. Spectra were acquired in a range of 50–1000 m/z. The analysis was made for separated in both positive and negative mode. The scan time was 500 ms (3 micro-scans). The ESI source parameters employed for positive mode were: capillary temperature 330 °C; capillary voltage 35 V; spray voltage 4 kV; Tube lens 80 V; nitrogen sheath gas 30 arbitrary units (AU); auxiliary gas 15 AU. For negative mode, the capillary voltage was −45 V and −123.38 V for tube lens. For the MS analyses, the collision energy (CID) was normalized to 35.

3.3. Spectra Data Analysis and Putative Identification of Metabolites

Raw mass spectrometry data processing for feature detection was performed in MZmine 2.37 software [60]. A base-line correction filter was applied and peak detection was performed with GridMass, with the following parameters: 0.32 m/z tolerance, 0.1 to 0.8 min width and a minimum peak height of 500. Isotopic peaks grouper was applied, and then, RANSAC (Random Sample Consensus) alignment was conducted with 0.35 m/z tolerance and 1.3 min of retention time. Finally, Gap filling (same RT and m/z) and duplicate peak filter were applied with 0.35 m/z tolerance. Intensities of features (m/z values) were employed for all statistical data analysis. Additionally, pools of samples were analyzed by high resolution-MS comprising a QTOF Synapt G1 (Waters, UK) in order to get accurate mass data. Three strategies were used for the putative identification of metabolites: (1) analysis of high-resolution MS data in XCMS [61] by matching with the METLIN database; (2) matching of accurate masses list by using a Coffea sp. in-house generated database with SpiderMass software [62]; (3) comparison of tandem MS data from LCQ Fleet Ion Trap with METLIN database, which was conducted in XCMS Online.

3.4. Cup Quality-Related Metabolite Analysis

3.4.1. Coffee Beans Processing and Sensory Analysis of Coffee Beverage

About five kilograms of fully ripe fruits were harvested for every plant included in the analysis. Fruit samples were subjected to post-harvest drying processing to obtain green coffee beans free of peel and parchment with moisture content about 12%. Green coffee beans of C. canephora and C. arabica (as reference samples) were medium roasted during 7–8 min according to SCAA specifications (Specialty Coffee Association of America 2015). The degree of roasting was standardized to a value of 55 (medium roast) in a colorimetric Agtron scale developed by SCAA. The coffee beverage was prepared using 50 g of roasted beans ground to medium size, and 1000 mL of boiled water (97 °C). After cool down to 50 °C, the beverage of all roasted coffee samples were tested by a trained panel of 11 trained cup-testers, who evaluated 10 attributes related to coffee cup quality: aromatic quality, aromatic intensity, astringency, body, flavor, after taste, bitterness, acidity, sourness and global preference, employing a scale ranging of six points (0 to 5) where 0 = no presence and 5 represents the higher intensity of the attribute.

3.4.2. Roasted Beans Metabolite Extraction

35 mg of ground roasted beans were extracted twice in 450 µL of 95 °C water in an ultrasound bath for 25 min. Samples were centrifuged at 15,870× g for 5 min at room temperature, and the supernatant was filtered through a 0.2 µm membrane and employed for UPLC-MS analysis.

4. Conclusions

Our study represents a comprehensive analysis of the heritability of metabolites in Coffea canephora plants and reports predictive metabolic markers for the coffee cup quality markers. We proved the feasibility to use metabolic profiles for the early selection of coffee plants in breeding projects and to speed up the development of new coffee varieties drastically. Compared to conventional breeding strategies, a ‘metabolome-wide association study’ (MWAS) permits a shorter time-to-market, and therefore, enables a more agile development of commercial products. Moreover, the creation of more resistant plant varieties with adequate cup quality is accelerated. The latter might become critical for responding to environmental challenges due to the climatic change.
MWAS is applicable to plant collections of any origin: natural varieties, conventional breeding and genetically modified organisms (GMOs). The analytical and statistical methods are highly scalable. Thus, MWAS has high potential in the breeding of crops of plants of agricultural or industrial interest. In further studies, we will test metabolic markers for sensory quality with C. arabica and study other complex traits, such as fungal resistance.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2218-1989/9/10/214/s1, Figure S1: Lipid classes of coffee roasted beans linked to sensory traits, Table S1. Metabolites of roasted coffee beans linked to coffee cup quality. Table S2. Highly heritable metabolites of coffee fruits and coffee leaf linked to sensory traits of coffee beverage.

Author Contributions

Conceptualization, R.G.-B. and R.W.; data curation, R.G.-B.; formal analysis, R.G.-B., Ó.G.-R. and J.J.O.-O.; funding acquisition, R.W.; investigation, R.G.-B., M.C.H.-H., Ó.G.-R., M.L.S.-Q., E.G.-P. and J.J.O.-O.; methodology, R.G.-B. and Ó.G.-R.; project administration, R.W.; resources, M.L.S.-Q., E.G.-P. and R.W.; supervision, R.W.; writing—original draft, R.G.-B.; writing—review and editing, R.W.

Funding

The project was funded by the CONACyT Fronteras project 2015-2/814, FIT S0017-2015-01, CONACyT-DFG 2016/277850, and SAGARPA-CONACyT 292546. RGB and MCHH acknowledge their CONACyT scholarships.

Acknowledgments

We thank Christophe Montagnon for fruitful scientific discussions, as well as Maria Teresa Carrillo Rayas, Maria Isabel Cristina Elizarraraz Anaya and Thermo Mexico for excellent technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Esquivel, P.; Jiménez, V.M. Functional properties of coffee and coffee by-products. Food Res. Int. 2012, 46, 488–495. [Google Scholar] [CrossRef]
  2. Farah, A. 15 - Coffee as a speciality and functional beverage A2 - Paquin, Paul. In Functional and Speciality Beverage Technology; Woodhead Publishing Series in Food Science, Technology and Nutrition; Woodhead Publishing: Cambridge, UK, 2009; pp. 370–395. ISBN 978-1-84569-342-8. [Google Scholar]
  3. Charrier, A.; Eskes, A.B. Botany and Genetics of Coffee. In Coffee: Growing, Processing, Sustainable Production; Wiley-VCH: Weinheim, Germany, 2008; pp. 25–56. ISBN 978-3-527-30731-9. [Google Scholar]
  4. Leroy, T.; Ribeyre, F.; Bertrand, B.; Charmetant, P.; Dufour, M.; Montagnon, C.; Marraccini, P.; Pot, D. Genetics of coffee quality. Braz. J. Plant Physiol. 2006, 18, 229–242. [Google Scholar] [CrossRef]
  5. Montagnon, C.; Marraccini, P.; Bertrand, B. Breeding for coffee quality. In Specialty Coffee: Managing Quality, Oberthür, T., Läderach, P., Pohlan, H.A.J., Cock, J., Eds.; IPNI: Southeast Asia, 2012; pp. 89–117. [Google Scholar]
  6. Philippe, L.; Benoít, B.; Hervé, E. Breeding coffee (Coffea arabica) for sustainable production. In Breeding Plantation Tree Crops: Tropical Species; Springer: Berlin/Heidelberg, Germany, 2009; pp. 525–543. [Google Scholar]
  7. Wintgens, J.N. Coffee: Growing, Processing, Sustainable Production. Available online: http://0-onlinelibrary-wiley-com.brum.beds.ac.uk/book/10.1002/9783527619627 (accessed on 23 March 2016).
  8. Leroy, T.; De Bellis, F.; Legnate, H.; Kananura, E.; Gonzales, G.; Pereira, L.F.; Andrade, A.C.; Charmetant, P.; Montagnon, C.; Cubry, P.; et al. Improving the quality of African robustas: QTLs for yield- and quality-related traits in Coffea canephora. Tree Genet. Genomes 2011, 7, 781–798. [Google Scholar] [CrossRef]
  9. Davis, A.P.; Gole, T.W.; Baena, S.; Moat, J. The Impact of Climate Change on Indigenous Arabica Coffee (Coffea arabica): Predicting Future Trends and Identifying Priorities. PLoS ONE 2012, 7, e47981. [Google Scholar] [CrossRef] [PubMed]
  10. Silva, M.D.C.; Várzea, V.; Guerra-Guimarães, L.; Azinheira, H.G.; Fernandez, D.; Petitot, A.-S.; Bertrand, B.; Lashermes, P.; Nicole, M. Coffee resistance to the main diseases: leaf rust and coffee berry disease. Braz. J. Plant Physiol. 2006, 18, 119–147. [Google Scholar] [CrossRef] [Green Version]
  11. Santos, C.; Leitao, A.; Pais, I.; Lidon, F.; Ramalho, J. Perspectives on the potential impacts of climate changes on coffee plant and bean quality. Emir. J. Food Agric. 2015, 27, 152. [Google Scholar] [CrossRef]
  12. Corso, M.P.; Vignoli, J.A.; Benassi, M. de T. Development of an instant coffee enriched with chlorogenic acids. J. Food Sci. Technol. 2016, 53, 1380–1388. [Google Scholar] [CrossRef]
  13. Farah, A. Coffee Constituents. In Coffee; Chu, Y.-F., Ed.; Wiley-Blackwell: Hoboken, NJ, USA, 2012; pp. 21–58. ISBN 978-1-119-94989-3. [Google Scholar]
  14. Akaffou, D.S.; Hamon, P.; Doulbeau, S.; Keli, J.; Legnate, H.; Campa, C.; Hamon, S.; de Kochko, A.; Zoro, B.I.A. Inheritance and relationship between key agronomic and quality traits in an interspecific cross between Coffea pseudozanguebariae Bridson and C. canephora Pierre. Tree Genet. Genomes 2012, 8, 1149–1162. [Google Scholar] [CrossRef]
  15. Ky, C.-L.; Louarn, J.; Guyot, B.; Charrier, A.; Hamon, S.; Noirot, M. Relations between and inheritance of chlorogenic acid contents in an interspecific cross between Coffea pseudozanguebariae and Coffea liberica var “dewevrei”. Theor. Appl. Genet. 1999, 98, 628–637. [Google Scholar] [CrossRef]
  16. Ky, C.-L.; Guyot, B.; Louarn, J.; Hamon, S.; Noirot, M. Trigonelline inheritance in the interspecific Coffea pseudozanguebariae×C. liberica var. dewevrei cross. Theor. Appl. Genet. 2001, 102, 630–634. [Google Scholar] [CrossRef]
  17. Mérot-L’Anthoëne, V.; Mangin, B.; Lefebvre-Pautigny, F.; Jasson, S.; Rigoreau, M.; Husson, J.; Lambot, C.; Crouzillat, D. Comparison of three QTL detection models on biochemical, sensory, and yield characters in Coffea canephora. Tree Genet. Genomes 2014, 10, 1541–1553. [Google Scholar] [CrossRef]
  18. Montagnon, C.; Guyot, B.; Cilas, C.; Leroy, T. Genetic parameters of several biochemical compounds from green coffee, Coffea canephora. Plant Breed. 1998, 117, 576–578. [Google Scholar] [CrossRef]
  19. Tran, H.T.; Lee, L.S.; Furtado, A.; Smyth, H.; Henry, R.J. Advances in genomics for the improvement of quality in coffee: Genomics for improving coffee quality. J. Sci. Food Agric. 2016, 96, 3300–3312. [Google Scholar] [CrossRef] [PubMed]
  20. Fernandez, O.; Urrutia, M.; Bernillon, S.; Giauffret, C.; Tardieu, F.; Le Gouis, J.; Langlade, N.; Charcosset, A.; Moing, A.; Gibon, Y. Fortune telling: metabolic markers of plant performance. Metabolomics 2016, 12. [Google Scholar] [CrossRef] [PubMed]
  21. Kumar, R.; Bohra, A.; Pandey, A.K.; Pandey, M.K.; Kumar, A. Metabolomics for Plant Improvement: Status and Prospects. Front. Plant Sci. 2017, 8, 1302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Cuperlovic-Culf, M.; Wang, L.; Forseille, L.; Boyle, K.; Merkley, N.; Burton, I.; Fobert, P.R. Metabolic Biomarker Panels of Response to Fusarium Head Blight Infection in Different Wheat Varieties. PloS ONE 2016, 11, e0153642. [Google Scholar] [CrossRef] [PubMed]
  23. Heuberger, A.L.; Broeckling, C.D.; Kirkpatrick, K.R.; Prenni, J.E. Application of nontargeted metabolite profiling to discover novel markers of quality traits in an advanced population of malting barley. Plant Biotechnol. J. 2014, 12, 147–160. [Google Scholar] [CrossRef] [PubMed]
  24. Kang, J.W.; Kim, H.-T.; Lee, W.Y.; Choi, M.N.; Park, E.-J. Identification of a potential metabolic marker, inositol, for the inherently fast growth trait by stems of Pinus densiflora via a retrospective approach. Can. J. For. Res. 2015, 45, 770–775. [Google Scholar] [CrossRef]
  25. Lima, M.R.M.; Felgueiras, M.L.; Graça, G.; Rodrigues, J.E.A.; Barros, A.; Gil, A.M.; Dias, A.C.P. NMR metabolomics of esca disease-affected Vitis vinifera cv. Alvarinho leaves. J. Exp. Bot. 2010, 61, 4033–4042. [Google Scholar] [CrossRef]
  26. López, M.G.; Zanor, M.I.; Pratta, G.R.; Stegmayer, G.; Boggio, S.B.; Conte, M.; Bermúdez, L.; Coluccio Leskow, C.; Rodríguez, G.R.; Picardi, L.A.; et al. Metabolic analyses of interspecific tomato recombinant inbred lines for fruit quality improvement. Metabolomics 2015, 11, 1416–1431. [Google Scholar] [CrossRef]
  27. Riedelsheimer, C.; Lisec, J.; Czedik-Eysenberg, A.; Sulpice, R.; Flis, A.; Grieder, C.; Altmann, T.; Stitt, M.; Willmitzer, L.; Melchinger, A.E. Genome-wide association mapping of leaf metabolic profiles for dissecting complex traits in maize. Proc. Natl. Acad. Sci. 2012, 109, 8872–8877. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Riedelsheimer, C.; Czedik-Eysenberg, A.; Grieder, C.; Lisec, J.; Technow, F.; Sulpice, R.; Altmann, T.; Stitt, M.; Willmitzer, L.; Melchinger, A.E. Genomic and metabolic prediction of complex heterotic traits in hybrid maize. Nat. Genet. 2012, 44, 217–220. [Google Scholar] [CrossRef] [PubMed]
  29. Steinfath, M.; Strehmel, N.; Peters, R.; Schauer, N.; Groth, D.; Hummel, J.; Steup, M.; Selbig, J.; Kopka, J.; Geigenberger, P.; et al. Discovering plant metabolic biomarkers for phenotype prediction using an untargeted approach. Plant Biotechnol. J. 2010, 8, 900–911. [Google Scholar] [CrossRef] [PubMed]
  30. Keurentjes, J.J.B.; Fu, J.; de Vos, C.H.R.; Lommen, A.; Hall, R.D.; Bino, R.J.; van der Plas, L.H.W.; Jansen, R.C.; Vreugdenhil, D.; Koornneef, M. The genetics of plant metabolism. Nat. Genet. 2006, 38, 842–849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Alseekh, S.; Bermudez, L.; de Haro, L.A.; Fernie, A.R.; Carrari, F. Crop metabolomics: from diagnostics to assisted breeding. Metabolomics 2018, 14, 148. [Google Scholar] [CrossRef] [Green Version]
  32. Davies, H.V.; Shepherd, L.V.T.; Stewart, D.; Frank, T.; Röhlig, R.M.; Engel, K.-H. Metabolome variability in crop plant species – When, where, how much and so what? Regul. Toxicol. Pharmacol. 2010, 58, S54–S61. [Google Scholar] [CrossRef] [PubMed]
  33. Eskes, A.B.; Leroy, T. Coffee Selection and Breeding. In Coffee: Growing, Processing, Sustainable Production: A Guidebook for Growers, Processors, Traders and Researchers; John Wiley & Sons: Hoboken, NJ, USA, 2008; pp. 57–86. ISBN 978-3-527-61962-7. [Google Scholar]
  34. Singh, R.J. Genetic Resources, Chromosome Engineering, and Crop Improvement: Medicinal Plants; CRC Press: Boca Raton, FL, USA, 2012; ISBN 978-1-4200-7384-3. [Google Scholar]
  35. Sunarharum, W.B.; Williams, D.J.; Smyth, H.E. Complexity of coffee flavor: A compositional and sensory perspective. Food Res. Int. 2014, 62, 315–325. [Google Scholar] [CrossRef]
  36. Vvan der Vossen, H.; Bertrand, B.; Charrier, A. Next generation variety development for sustainable production of arabica coffee (Coffea arabica L.): a review. Euphytica 2015, 204, 243–256. [Google Scholar] [CrossRef]
  37. Bertrand, B.; Etienne, H.; Cilas, C.; Charrier, A.; Baradat, P. Coffea arabica hybrid performance for yield, fertility and bean weight. Euphytica 2005, 141, 255–262. [Google Scholar] [CrossRef]
  38. Costa Freitas, A.M.; Mosca, A.I. Coffee geographic origin — an aid to coffee differentiation. Food Res. Int. 1999, 32, 565–573. [Google Scholar] [CrossRef]
  39. Bertrand, B.; Vaast, P.; Alpizar, E.; Etienne, H.; Davrieux, F.; Charmetant, P. Comparison of bean biochemical composition and beverage quality of Arabica hybrids involving Sudanese-Ethiopian origins with traditional varieties at various elevations in Central America. Tree Physiol. 2006, 26, 1239–1248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Bosselmann, A.S.; Dons, K.; Oberthur, T.; Olsen, C.S.; Ræbild, A.; Usma, H. The influence of shade trees on coffee quality in small holder coffee agroforestry systems in Southern Colombia. Agric. Ecosyst. Environ. 2009, 129, 253–260. [Google Scholar] [CrossRef]
  41. Bhumiratana, N.; Adhikari, K.; Chambers IV, E. Evolution of sensory aroma attributes from coffee beans to brewed coffee. LWT - Food Sci. Technol. 2011, 44, 2185–2192. [Google Scholar] [CrossRef] [Green Version]
  42. Poltronieri, Y.; Martinez, H.E.P.; Cecon, P.R. Effect of zinc and its form of supply on production and quality of coffee beans. J. Sci. Food Agric. 2011, 91, 2431–2436. [Google Scholar] [CrossRef]
  43. Buffo, R.A.; Cardelli-Freire, C. Coffee flavour: an overview. Flavour Fragr. J. 2004, 19, 99–104. [Google Scholar] [CrossRef]
  44. Viani, R.; Illy, A. Espresso Coffee: The Science of Quality, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2005; ISBN 978-0-12-370371-2. [Google Scholar]
  45. Baggenstoss, J.; Poisson, L.; Kaegi, R.; Perren, R.; Escher, F. Coffee Roasting and Aroma Formation: Application of Different Time−Temperature Conditions. J. Agric. Food Chem. 2008, 56, 5836–5846. [Google Scholar] [CrossRef]
  46. Ribeiro, J.S.; Augusto, F.; Salva, T.J.G.; Thomaziello, R.A.; Ferreira, M.M.C. Prediction of sensory properties of Brazilian Arabica roasted coffees by headspace solid phase microextraction-gas chromatography and partial least squares. Anal. Chim. Acta 2009, 634, 172–179. [Google Scholar] [CrossRef]
  47. Zhou, L.; Khalil, A.; Bindler, F.; Zhao, M.; Marcic, C.; Ennahar, S.; Marchioni, E. Effect of heat treatment on the content of individual phospholipids in coffee beans. Food Chem. 2013, 141, 3846–3850. [Google Scholar] [CrossRef]
  48. Lang, R.; Klade, S.; Beusch, A.; Dunkel, A.; Hofmann, T. Mozambioside Is an Arabica-Specific Bitter-Tasting Furokaurane Glucoside in Coffee Beans. J. Agric. Food Chem. 2015, 63, 10492–10499. [Google Scholar] [CrossRef]
  49. Iwasa, K.; Setoyama, D.; Shimizu, H.; Seta, H.; Fujimura, Y.; Miura, D.; Wariishi, H.; Nagai, C.; Nakahara, K. Identification of 3-Methylbutanoyl Glycosides in Green Coffea arabica Beans as Causative Determinants for the Quality of Coffee Flavors. J. Agric. Food Chem. 2015, 63, 3742–3751. [Google Scholar] [CrossRef]
  50. Stadler, R.H.; Varga, N.; Hau, J.; Vera, F.A.; Welti, D.H. Alkylpyridiniums. 1. Formation in Model Systems via Thermal Degradation of Trigonelline. J. Agric. Food Chem. 2002, 50, 1192–1199. [Google Scholar] [CrossRef]
  51. Bicho, N.C.; Leitão, A.E.; Ramalho, J.C.; Lidon, F.C. Identification of chemical clusters discriminators of the roast degree in Arabica and Robusta coffee beans. Eur. Food Res. Technol. 2011, 233, 303–311. [Google Scholar] [CrossRef]
  52. Los Santos-Briones, D.; Hernández-Sotomayor, S.M. Coffee biotechnology. Braz. J. Plant Physiol. 2006, 18, 217–227. [Google Scholar] [CrossRef]
  53. Mishra, M.K.; Slater, A. Recent Advances in the Genetic Transformation of Coffee. Biotechnol. Res. Int. 2012, 2012, 1–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Schauer, N.; Semel, Y.; Balbo, I.; Steinfath, M.; Repsilber, D.; Selbig, J.; Pleban, T.; Zamir, D.; Fernie, A.R. Mode of Inheritance of Primary Metabolic Traits in Tomato. Plant Cell Online 2008, 20, 509–523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Toubiana, D.; Semel, Y.; Tohge, T.; Beleggia, R.; Cattivelli, L.; Rosental, L.; Nikoloski, Z.; Zamir, D.; Fernie, A.R.; Fait, A. Metabolic Profiling of a Mapping Population Exposes New Insights in the Regulation of Seed Metabolism and Seed, Fruit, and Plant Relations. PLoS Genet. 2012, 8, e1002612. [Google Scholar] [CrossRef] [PubMed]
  56. Alseekh, S.; Tohge, T.; Wendenberg, R.; Scossa, F.; Omranian, N.; Li, J.; Kleessen, S.; Giavalisco, P.; Pleban, T.; Mueller-Roeber, B.; et al. Identification and Mode of Inheritance of Quantitative Trait Loci for Secondary Metabolite Abundance in Tomato. Plant Cell 2015, 27, 485–512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Stewart, D.; McDougall, G.J.; Sungurtas, J.; Verrall, S.; Graham, J.; Martinussen, I. Metabolomic approach to identifying bioactive compounds in berries: advances toward fruit nutritional enhancement. Mol. Nutr. Food Res. 2007, 51, 645–651. [Google Scholar] [CrossRef]
  58. Barbosa, J.N.; Borem, F.M.; Cirillo, M.A.; Malta, M.R.; Alvarenga, A.A.; Alves, H.M.R. Coffee Quality and Its Interactions with Environmental Factors in Minas Gerais, Brazil. J. Agric. Sci. 2012, 4, 181. [Google Scholar] [CrossRef]
  59. Bertrand, B.; Boulanger, R.; Dussert, S.; Ribeyre, F.; Berthiot, L.; Descroix, F.; Joët, T. Climatic factors directly impact the volatile organic compound fingerprint in green Arabica coffee bean as well as coffee beverage quality. Food Chem. 2012, 135, 2575–2583. [Google Scholar] [CrossRef]
  60. Pluskal, T.; Castillo, S.; Villar-Briones, A.; Orešič, M. MZmine 2: Modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinformatics 2010, 11, 395. [Google Scholar] [CrossRef] [PubMed]
  61. Gowda, H.; Ivanisevic, J.; Johnson, C.H.; Kurczy, M.E.; Benton, H.P.; Rinehart, D.; Nguyen, T.; Ray, J.; Kuehl, J.; Arevalo, B.; et al. Interactive XCMS Online: Simplifying Advanced Metabolomic Data Processing and Subsequent Statistical Analyses. Anal. Chem. 2014, 86, 6931–6939. [Google Scholar] [CrossRef] [PubMed]
  62. Winkler, R. SpiderMass: Semantic database creation and tripartite metabolite identification strategy. J. Mass Spectrom. 2015, 50, 538–541. [Google Scholar] [CrossRef] [PubMed]
Metabolites 09 00214 g001 550
Figure 1. The tanglegram compares the hierarchical clustering of mother plants employing fruits versus leaves metabolite profiles. Individual dendrograms were calculated using hierarchical clustering with the ward.D2 method. The signal intensity of every peak mass used for generation of dendrograms consisted of the average of three technical replicates.
Figure 1. The tanglegram compares the hierarchical clustering of mother plants employing fruits versus leaves metabolite profiles. Individual dendrograms were calculated using hierarchical clustering with the ward.D2 method. The signal intensity of every peak mass used for generation of dendrograms consisted of the average of three technical replicates.
Metabolites 09 00214 g001
Metabolites 09 00214 g002 550
Figure 2. Frequency distribution of heritability of mass peaks detected in fruits (a) and leaves (b). The histograms demonstrate the distribution of the heritability level. Colors depict values according to low, moderate and high heritability levels.
Figure 2. Frequency distribution of heritability of mass peaks detected in fruits (a) and leaves (b). The histograms demonstrate the distribution of the heritability level. Colors depict values according to low, moderate and high heritability levels.
Metabolites 09 00214 g002
Metabolites 09 00214 g003 550
Figure 3. The correlation plot shows histograms and density distribution of cup quality attributes of coffee beverage. The Spearman correlation coefficient (SCC) between sensory traits of the coffee beverage is depicted. Data with statistical significance levels are indicated with asterisks: p-values ≤ 0.05 *, ≤0.01 ** and ≤0.001 ***.
Figure 3. The correlation plot shows histograms and density distribution of cup quality attributes of coffee beverage. The Spearman correlation coefficient (SCC) between sensory traits of the coffee beverage is depicted. Data with statistical significance levels are indicated with asterisks: p-values ≤ 0.05 *, ≤0.01 ** and ≤0.001 ***.
Metabolites 09 00214 g003
Metabolites 09 00214 g004 550
Figure 4. Cluster dendrogram of correlation between cup quality attributes of coffee beverage. Hierarchical clustering was calculated with the ward.D2 method; correlation distances were calculated by employing the scores of every sensory trait determined by the trained cup-tester panel.
Figure 4. Cluster dendrogram of correlation between cup quality attributes of coffee beverage. Hierarchical clustering was calculated with the ward.D2 method; correlation distances were calculated by employing the scores of every sensory trait determined by the trained cup-tester panel.
Metabolites 09 00214 g004
Metabolites 09 00214 g005 550
Figure 5. Classification of roasted bean metabolites linked to sensory traits of coffee beverage. Only putatively identified metabolites were included for compound classification by using KEGG and Pubchem tools.
Figure 5. Classification of roasted bean metabolites linked to sensory traits of coffee beverage. Only putatively identified metabolites were included for compound classification by using KEGG and Pubchem tools.
Metabolites 09 00214 g005
Metabolites 09 00214 g006 550
Figure 6. Metabolites of roasted coffee beans linked to multiple sensory traits. The figure shows a significant Pearson correlation (positive and negative) of quercetin 4´glucoside and trigonelline with multiple sensory traits of coffee beverage.
Figure 6. Metabolites of roasted coffee beans linked to multiple sensory traits. The figure shows a significant Pearson correlation (positive and negative) of quercetin 4´glucoside and trigonelline with multiple sensory traits of coffee beverage.
Metabolites 09 00214 g006
Metabolites 09 00214 g007 550
Figure 7. Distribution of abundances of top metabolites in parent C. canephora plant tissues. Level distribution across the 40 C. canephora parent plants are depicted for the metabolites listed in Table 3. Yellow-colored boxplots indicate metabolites found in fruits and blue color represents metabolites from the coffee leaf.
Figure 7. Distribution of abundances of top metabolites in parent C. canephora plant tissues. Level distribution across the 40 C. canephora parent plants are depicted for the metabolites listed in Table 3. Yellow-colored boxplots indicate metabolites found in fruits and blue color represents metabolites from the coffee leaf.
Metabolites 09 00214 g007
Table
Table 1. C. Canephora families included in heritability analysis.
Table 1. C. Canephora families included in heritability analysis.
NumberMother Plant IDSib 1Sib 2NumberMother Plant IDSib 1Sib 2
1MP196196-1196-321MP349349-2349-3
2MP203203-1203-222MP360360-3360-4
3MP205205-1205-223MP390390-1390-3
4MP212212-1212-224MP408408-1408-4
5MP292292-3292-425MP416416-1416-3
6MP295295-2295-326MP418418-3418-4
7MP312312-3312-427MP419419-2419-4
8MP315315-1315-228MP420420-1420-4
9MP316316-1316-229MP430430-2430-3
10MP318318-2318-330MP432432-1432-2
11MP322322-3322-431MP433433-2433-3
12MP323323-1323-432MP434434-2434-4
13MP328328-3328-433MP442442-2442-3
14MP333333-1333-434MP443443-3443-4
15MP336336-1336-235MP449449-1449-3
16MP337337-3337-436MP450450-3450-4
17MP339339-1339-237MP451451-1451-4
18MP340340-1340-338MP452452-2452-3
19MP343343-1343-339MP453453-2453-3
20MP346346-1346-440MP454454-2454-4
Table
Table 2. Sensory profile of coffee beverages.
Table 2. Sensory profile of coffee beverages.
SampleAromatic
Intensity		Aromatic
Quality		FlavorAstringencyAftertasteAcidityBitternessBodySournessGlobal Preference
MP3023.27 ± 0.52 3.09 ± 0.332.73 ± 0.392.45 ± 0.493.27 ± 0.390.64 ± 0.463.45 ± 0.492.73 ± 0.520.27 ± 0.392.82 ± 0.29
MP3072.14 ± 0.331.86 ± 0.472.23 ± 0.412.05 ± 0.332.27 ± 0.380.55 ± 0.333.50 ± 0.542.23 ± 0.610.36 ± 0.462.14 ± 0.24
MP3093.18 ± 0.443.45 ± 0.493.09 ± 0.332.82 ± 0.293.00 ± 0.180.45 ± 0.492.82 ± 0.292.82 ± 0.290.55 ± 0.503.09 ± 0.16
MP3242.55 ± 0.492.64 ± 0.462.18 ± 0.292.45 ± 0.492.55 ± 0.670.55 ± 0.493.18 ± 0.472.36 ± 0.460.45 ± 0.492.36 ± 0.46
MP3162.64 ± 0.462.45 ± 0.492.91 ± 0.332.36 ± 0.462.64 ± 0.570.64 ± 0.463.73 ± 0.392.55 ± 0.490.45 ± 0.492.91 ± 0.16
MP3182.82 ± 0.292.64 ± 0.462.36 ± 0.463.09 ± 0.163.00 ± 0.180.55 ± 0.493.55 ± 0.492.73 ± 0.390.45 ± 0.492.55 ± 0.49
MP3233.00 ± 0.183.09 ± 0.163.09 ± 0.331.73 ± 0.392.64 ± 0.460.36 ± 0.463.55 ± 0.492.36 ± 0.570.18 ± 0.293.09 ± 0.33
MP3462.27 ± 0.392.09 ± 0.162.00 ± 0.362.09 ± 0.162.73 ± 0.390.45 ± 0.503.45 ± 0.492.73 ± 0.391.00 ± 0.182.00 ± 0.18
MP3602.64 ± 0.522.36 ± 0.461.82 ± 0.292.18 ± 0.292.55 ± 0.491.27 ± 1.072.36 ± 0.822.45 ± 0.771.27 ± 0.391.82 ± 0.29
MP4083.27 ± 0.392.82 ± 0.292.36 ± 0.462.82 ± 0.293.00 ± 0.360.27 ± 0.393.36 ± 0.462.73 ± 0.390.36 ± 0.362.82 ± 0.29
MP4202.73 ± 0.393.00 ± 0.183.00 ± 0.362.55 ± 0.492.91 ± 0.331.18 ± 0.293.36 ± 0.462.36 ± 0.460.55 ± 0.492.45 ± 0.49
MP4532.45 ± 0.491.82 ± 0.291.82 ± 0.472.45 ± 0.492.73 ± 0.391.00 ± 0.18 3.09 ± 0.162.64 ± 0.460.55 ± 0.491.91 ± 0.16
MP4542.91 ± 0.332.00 ± 0.181.91 ± 0.172.73 ± 0.392.55 ± 0.490.73 ± 0.42.91 ± 0.332.82 ± 0.290.64 ± 0.461.91 ± 0.16
M4952.36 ± 0.461.45 ± 0.491.45 ± 0.491.55 ± 0.591.73± 0.390.29 ± 0.742.27 ± 0.392.27 ± 0.390.91 ± 0.331.73 ± 0.52
M4982.73 ± 0.391.91 ± 0.161.36 ± 0.462.27 ± 0.392.09 ± 0.163.00 ± 0.181.91 ± 0.332.27 ± 0.402.00 ± 0.360.45 ± 0.49
M4993.09 ± 0.331.82 ± 0.291.55 ± 0.492.18 ± 0.291.73 ± 0.392.09 ± 0.492.18 ± 0.292.55 ± 0.491.18 ± 0.301.18 ± 0.29
M5003.00 ± 0.252.38 ± 0.462.25 ± 0.382.50 ± 0.502.50 ± 0.503.38 ± 0.462.25 ± 0.372.63 ± 0.461.63 ± 0.471.75 ± 0.37
M5023.43 ± 0.482.29 ± 0.401.71 ± 0.401.43 ± 0.482.43 ± 0.483.14 ± 0.242.29 ± 0.401.86 ± 0.241.14 ± 0.241.43 ± 0.49
M5062.88 ± 0.221.25 ± 0.381.50 ± 0.502.00 ± 0.252.50 ± 0.503.88 ± 0.441.75 ± 0.382.00 ± 0.253.13 ± 0.440.75 ± 0.38
M5072.57 ± 0.491.29 ± 0.411.57 ± 0.492.29 ± 0.691.57 ± 0.494.00 ± 0.571.43 ± 0.611.71 ± 0.411.43 ± 0.490.43 ± 0.49
M5082.71 ± 0.262.18 ± 0.292.00 ± 0.181.71 ± 0.401.43 ± 0.493.86 ± 0.241.29 ± 0.402.14 ± 0.241.14 ± 0.241.00 ± 0.28
The values represent the average value with the standard deviation of 11 trained cup-testers.
Table
Table 3. Top metabolites of coffee fruits and leaf with optimal characteristics (high heritability and strong correlation) linked to sensory traits.
Table 3. Top metabolites of coffee fruits and leaf with optimal characteristics (high heritability and strong correlation) linked to sensory traits.
Sensory AttributeHR m/zAdduct IonNameTissueHeritabilitySCCp-Value
Acidity963.4134n.a.UnknownFruits0.960.6674.98 × 10-2
Aftertaste940.2271[M+Na]+2,4-Decadienoyl-CoAFruits0.8810.7541.88 × 10-2
Aromatic intensity860.5079n.a.UnknownFruits0.9990.8127.89 × 10-3
Aromatic quality536.0019n.a.UnknownFruits0.6560.8009.63 × 10-3
Astringency601.1263[M+Na]+2,3-trans-proanthocyanidinFruits0.9570.8335.27 × 10-3
Bitterness734.4907n.a.UnknownFruits0.8170.7402.27 × 10-2
Body789.1658[M+H]+Delphinidin 3-O-glucosyl-5-O-caffeoylglucosideFruits0.9490.6903.98 × 10-2
Flavor598.5406n.a.UnknownFruits0.9980.7402.27 × 10-2
Global preference791.4983[M+H]+Nonaprenyl diphosphateFruits0.9840.7282.61 × 10-2
Sourness521.5618n.a.UnknownFruits0.8420.6894.00 × 10-2
Acidity473.3339[M+H]+6-DeoxoteasteroneLeaf0.6430.7831.13 × 10-2
Aftertaste695.4571[M+Na]+PA(16:0/18:2)Leaf0.5350.7372.40 × 10-2
Aromatic intensity713.317n.a.UnknownLeaf0.9580.8039.01 × 10-3
Aromatic quality251.2095[M+Na]+Myristic acidLeaf0.9940.7661.62 × 10-2
Astringency967.7852[M]−MGDG(23:0/26:0)Leaf0.6120.6664.90 × 10-2
*Bitterness473.3339[M+H]+6-DeoxoteasteroneLeaf0.643−0.8573.14 × 10-3
*Body463.3451[M+Na]+4,4-DimethylzymosterolLeaf0.891−0.6893.98 × 10-2
Flavor251.2095[M+Na]+Myristic acidLeaf0.9940.6724.89 × 10-2
Global preference433.1626[M+K]+Gibberellin A28Leaf0.9960.6774.48 × 10-2
Sourness473.3339[M+H]+6-DeoxoteasteroneLeaf0.6430.7731.47 × 10-2
Sensory attributes are presented by alphabetic order and metabolites are showed in decreasing order of heritability level for positive correlation. In some cases, only a negative correlation was detected. HR: High resolution; SCC: Spearman correlation coefficient; n.a.: not applicable; * only negative correlation was detected for this sensory attributes.

Share and Cite

MDPI and ACS Style

Gamboa-Becerra, R.; Hernández-Hernández, M.C.; González-Ríos, Ó.; Suárez-Quiroz, M.L.; Gálvez-Ponce, E.; Ordaz-Ortiz, J.J.; Winkler, R. Metabolomic Markers for the Early Selection of Coffea canephora Plants with Desirable Cup Quality Traits. Metabolites 2019, 9, 214. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo9100214

AMA Style

Gamboa-Becerra R, Hernández-Hernández MC, González-Ríos Ó, Suárez-Quiroz ML, Gálvez-Ponce E, Ordaz-Ortiz JJ, Winkler R. Metabolomic Markers for the Early Selection of Coffea canephora Plants with Desirable Cup Quality Traits. Metabolites. 2019; 9(10):214. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo9100214

Chicago/Turabian Style

Gamboa-Becerra, Roberto, María Cecilia Hernández-Hernández, Óscar González-Ríos, Mirna L. Suárez-Quiroz, Eligio Gálvez-Ponce, José Juan Ordaz-Ortiz, and Robert Winkler. 2019. "Metabolomic Markers for the Early Selection of Coffea canephora Plants with Desirable Cup Quality Traits" Metabolites 9, no. 10: 214. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo9100214

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

Article Metrics

Back to TopTop