Next Article in Journal
Molecular, Metabolic, and Nutritional Changes after Metabolic Surgery in Obese Diabetic Patients (MoMen): A Protocol for a Multicenter Prospective Cohort Study
Previous Article in Journal
Analysis of Metabolite Distribution in Rat Liver of High-Fat Model by Mass Spectrometry Imaging
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Medium Roasting and Brewing Methods Differentially Modulate Global Metabolites, Lipids, Biogenic Amines, Minerals, and Antioxidant Capacity of Hawai‘i-Grown Coffee (Coffea arabica)

1
Laboratory of Metabolic Disorders and Alternative Medicine, Department of Molecular Biosciences and Bioengineering (MBBE), College of Tropical Agriculture and Human Resources (CTAHR), University of Hawai‘i at Manoa (UHM), Honolulu, HI 96822, USA
2
Department of Tropical Plant and Soil Sciences (TPSS), CTAHR, UHM, Honolulu, HI 96822, USA
3
Bioinformatics Core, Departmentt of Quantitative Health Sciences, University of Hawai‘i Cancer Center (UHCC), John A. Burns School of Medicine (JABSOM), UHM, Honolulu, HI 96813, USA
*
Author to whom correspondence should be addressed.
Metabolites 2023, 13(3), 412; https://doi.org/10.3390/metabo13030412
Submission received: 6 January 2023 / Revised: 1 March 2023 / Accepted: 6 March 2023 / Published: 10 March 2023
(This article belongs to the Section Food Metabolomics)

Abstract

:
In the United States, besides the US territory Puerto Rico, Hawai‘i is the only state that grows commercial coffee. In Hawai’i, coffee is the second most valuable agricultural commodity. Health benefits associated with moderate coffee consumption, including its antioxidant capacity, have been correlated to its bioactive components. Post-harvest techniques, coffee variety, degree of roasting, and brewing methods significantly impact the metabolites, lipids, minerals, and/or antioxidant capacity of brewed coffees. The goal of our study was to understand the impact of roasting and brewing methods on metabolites, lipids, biogenic amines, minerals, and antioxidant capacity of two Hawai‘i-grown coffee (Coffea arabica) varieties, “Kona Typica” and “Yellow Catuai”. Our results indicated that both roasting and coffee variety significantly modulated several metabolites, lipids, and biogenic amines of the coffee brews. Furthermore, regardless of coffee variety, the antioxidant capacity of roasted coffee brews was higher in cold brews. Similarly, total minerals were higher in “Kona Typica” cold brews followed by “Yellow Catuai” cold brews. Hawai‘i-grown coffees are considered “specialty coffees” since they are grown in unique volcanic soils and tropical microclimates with unique flavors. Our studies indicate that both Hawai‘i-grown coffees contain several health-promoting components. However, future studies are warranted to compare Hawai‘i-grown coffees with other popular brand coffees and their health benefits in vivo.

1. Introduction

Coffee is one of the most popular drinks worldwide and is a widely consumed beverage. During 2022–2023, global coffee consumption is expected to reach over 168.7 million 60 kg (132 lb.) bags [1], and it is forecasted that the global production will be 172.8 million bags in the 2022–2023 growing season [2]. In the United States, it is estimated that more than 157 million cups of coffee are consumed per day, making the United States the world’s leading coffee consumer [3]. The health benefits of coffee, including its antioxidant properties, are associated with its complex array of bioactive chemicals, the most influential being alkaloids (caffeine and trigonelline), phenolic compounds (chlorogenic acids), and diterpenes (cafestol and kahweol). Epidemiological studies have reported several health benefits of coffee such as reducing the risk of metabolic syndrome (MetS), type 2 diabetes (T2D), cardiovascular disease (CVD), various types of cancer, kidney stones, liver disease, Parkinson’s disease, gout, and neurological disorders [4,5,6,7]. Moderate coffee consumption may also confer protective effects against overall mortality [8,9,10,11,12,13]. However, a few studies indicate that drinking more than four to nine cups of coffee per day may have some detrimental health effects such as increasing plasma cholesterol or triglyceride levels, and decreasing bone density in women [6,14,15,16]. Consuming large amounts of caffeinated coffee can also have negative effects on health such as increased blood pressure, anxiety, or difficulty falling asleep [17]. Abrupt cessation of caffeine consumption may induce withdrawal symptoms such as headache, fatigue, and/or depression. Coffee consumption during pregnancy is associated with adverse birth outcomes and neonatal health including lower infant birth weight [18,19,20,21,22,23,24,25,26,27,28]. Similarly, biogenic amines identified in coffee brews, resulting from amino acid decarboxylation, can be toxic to humans at high concentrations [29]. Inconsistencies observed in the health benefits of coffee may possibly be associated with differences in the chemical compositions of coffee [30,31,32,33,34,35,36].
It has been demonstrated that differences in coffee metabolites and caffeine levels [36,37] are influenced by their geographic origins or environment [38,39,40,41,42,43,44], post-harvest processing [30,33,45,46,47,48], instant versus fresh ground coffees [49,50], degrees of roasting [33,51,52,53], and/or types of brewing methods [31,33,36,49,54,55,56]. Different coffee metabolites not only impart flavors [38,57,58] or aromas [59,60,61,62,63], but also determine therapeutic impact by influencing antioxidant potential [34,64,65,66,67]. Prospective epidemiological studies conducted among populations of Hawai‘i indicate that a moderate level of coffee consumption (one to three cups per day) was associated with a lower risk of T2D, chronic liver diseases, dementia, and Parkinson’s disease [68,69,70,71]. In contrast, drinking more than nine cups of coffee per day was positively associated with increased serum cholesterol among several population-based studies, including Japanese men in Hawai‘i [14,15].
Besides the US territory Puerto Rico, Hawai‘i is the only state in the United States that grows commercial coffee. In Hawai’i, coffee is the second most valuable agriculture commodity after seed crops [72,73]. In 2021, unroasted coffee was valued at USD 102.91 million, and its roasted value was more than USD 148.48 million [74]. In 2022, annual revenue for the US coffee industry was estimated at USD 90.27 billion [75]. Hawai’i-grown coffees account for less than 1% of total global coffee production. However, they are considered specialty coffees since they are grown in unique volcanic soils and tropical microclimates with unique flavors [76]. The main varieties of coffee grown in Hawai‘i are “Kona Typica” and “Yellow Catuai”, which represent two major Coffea arabica botanical varieties of “Typica” and “Bourbon” used worldwide.
Metabolomic profiles of coffees brewed by different methods have been reported [33,36]. However, studies on the effect of roasting conditions on coffee metabolites are limited [33]. The antioxidant capacities and mineral contents of several coffee varieties and the impact of roasting conditions have also been reported [64,77,78,79,80,81,82,83]. Green coffee beans (grounds) are generally not brewed for consumption. Only roasted coffee grounds are brewed for consumption worldwide. Roasting generates the flavor of coffee. Green coffee brews would not be considered “coffee” since they do not taste like the “coffee” brewed from roasted grounds. Overall, no studies have addressed the lipidomic profiles of brewed coffee from roasted beans, while only one study evaluated selected biogenic amines in brewed coffee [29]. To date, one study has identified metabolites from spent grounds of Hawai‘ian Kona coffee [84]. However, variability in metabolites, lipids, biogenic amines, mineral contents, or antioxidant capacity of Hawai‘i-grown coffee (Coffea arabica) has yet to be elucidated.
The primary objective of our study was to identify the impact of roasting and brewing methods on global metabolites, lipids, and biogenic amines in the two Hawai‘i-grown coffee varieties, “Kona Typica” and “Yellow Catuai”. Drip filter paper, drip filter mesh and French press are the most widely used methods for brewing roasted coffee. For our initial investigation, to evaluate the differences in metabolites, lipids, and biogenic amines among brews of roasted and green coffee beans, we used the above three brewing methods of coffee preparation. However, to understand the effects of different methods of brewing roasted coffee, we also included the cold brewing method, which is gaining popularity worldwide and has been reported to provide more beneficial compounds than other brewing methods [85,86]. The secondary objective of our study was to evaluate the influence of brewing methods on the antioxidant capacity and mineral contents of these two Hawai‘i-grown coffees. Brewing green coffee bean grounds is also not common. Therefore, we analyzed the minerals and the antioxidant capacity of brews prepared only from roasted coffee grounds.

2. Materials and Methods

Green coffee beans, 50 lbs. each, of “Kona Typica” (purchased from Waialua Coffee and Cacao Estate, Oahu, Honolulu, HI, USA) and “Yellow Catuai” (purchased from Kaua’i Coffee Company, Kaua’i, Kalaheo, HI, USA) were stored at room temperature. “Kona Typica” is also grown in Kaua’i, and “Yellow Catuai” is grown in Kona and on Oahu. Although metabolite differences are noted based on genetic variety, geographic location can also make a difference [38,39,40,41,42,43,44,46]. For ease of readership and future study reference, we will refer to “Kona Typica” as Waialua and “Yellow Catuai” as Kaua’i coffee.

2.1. Green Coffee Grounds

Nylon bags filled with green coffee were dunked in liquid nitrogen for three to five seconds in order to freeze the green coffee beans. The frozen, dried green coffee beans were ground using a “Santos burr grinder” (Sao Paulo, Brazil). Grounds were collected in a jar, and larger grounds that did not pass completely through were discarded. Green coffee grounds were brewed using the drip filter paper, drip filter mesh and French press methods described below.

2.2. Coffee Roasting

Prime-grade green beans were further cleaned by hand to remove any defective green beans and roasted using the electric, programmable rotating drum type “Has Garanti roaster” (Model HSR 1 kg (92.2 lbs.), Turkey), at 230 °C temperature for about 13 min, until coffee reached a medium roast color according to the “Roast Color Classification System” (Agtron 55-65, SCAA, Long Beach, CA, USA), and cooled at room temperature. Coffee beans were ground using the “Santos coffee grinder” (Burr type, Santa Fe Springs, CA, USA) to a size for paper filter brewing. For each type of brewing method, three roasted batches were mixed thoroughly and pooled as one sample. To understand inter-roasting variability, a total of six batches were roasted, providing two pooled samples for preparing the different brews.

2.3. Drip Filter Paper (FP) Method

Two tablespoons (Tbsp) of coffee grounds with six fluid ounces (fl oz) of room-temperature tap water was brewed in the “Toastmaster” coffee machine (Star International Holdings group of brands, Star manufacturing, Smithville, TN, USA) using the “Total Home #4 Cone Style paper Filter” (CVS Pharmacy, Honolulu, HI, USA). After brewing was completed, the hot plate was turned off, and the coffee was cooled in the pot for 15 min. Cooled coffee brews were aliquoted into 50 mL tubes and stored at −80 °C until lyophilization.

2.4. Drip Metallic Filter Mesh (FM) Method

Six fl oz of room-temperature tap water with two Tbsp of coffee grounds was brewed in the “Black and Decker Brew ‘N Go” (San Diego, CA, USA) coffee machine, which contains a metallic (steel) mesh filter. After brewing was completed, coffee was cooled in “Brew ‘N Go cup” for 15 min, aliquoted into 50 mL tubes, and stored at −80 °C until lyophilization.

2.5. French Press (FrP) Method

Six fl oz of boiling water with two Tbsp of coffee grounds was set to steep for 15 min in the “Bodum Brazil 3 cup French Press Coffee Maker 12 oz” (Bodum incorporation, Triengen, Switzerland). The ground coffee beans were pressed with a perforated plunger plate and then cooled for 20 min in the brewing container. The coffee brew was then aliquoted into 50 mL tubes and stored at −80 °C until lyophilization.

2.6. Cold Brew (CB) Method

A damp, reusable “Toddy Filter” (fabric-like, compostable filter made from tree-free specialized material by Toddy, LLC (Loveland, CO, USA) ) was placed in the bottom of the “Toddy Cold Brew System” (Toddy, LLC, Loveland, CO, USA). One fl oz of room-temperature tap water and 4 Tbsp of roasted coffee grounds were added to the system. An additional three fl oz of water was gently added to the mixture and left to sit for 5 min. Another batch of four Tbsp of grounds and three fl oz of water was carefully added. The grounds were lightly pressed down to ensure all grounds were wet. The system was covered with foil and kept to steep for 24 h. After steeping, the coffee was measured and 2 times the amount of water was added to dilute the concentrate. Ten-milliliter aliquots were distributed, frozen, and stored in a −80 °C freezer until lyophilization. Samples were lyophilized in a Martin Christ Alpha 2-4 LD plus (Christ, Osterode am Harz, Germany) for 24 h, pooled, and then stored in a −80 °C freezer until analysis.

2.7. Brewed Coffee Omics (Global Metabolites, Lipids, and Biogenic Amines)

Roasted coffee grounds were brewed by drip filter paper, drip filter mesh, French press, and cold brew methods. Green coffee bean grounds were brewed by drip filter paper, drip filter mesh, and French press methods. Metabolomics, lipidomics, and analysis of biogenic amines were conducted at the Fiehn Laboratory, NIH West Coast Metabolomics Center. Global metabolites (targeted and untargeted) were analyzed using an automated liner exchange cold injection system gas chromatography time of flight mass spectrometer (ALEX-CIS GCTOF MS) as described previously [87,88,89,90,91]. In brief, 10 mg brewed coffee samples were extracted with 1 mL of 3:3:2 acetonitrile (ACN):isopropanol (IPA):water by vortexing for 10 s and shaking for 6 min at 4 °C. After centrifugation at 14,000 RCF (relative centrifugal force) for 2 min, the supernatant was aliquoted into 475 μL aliquots, dried, and stored until further analysis. Half of the dried sample was derivatized with 10 μL of 40 mg·mL−1 of methoxyamine in pyridine and shaken at 30 °C for 1.5 h. Ninety-one microliters of N-methyl-N-(trimethylsilyl) trifluoroacetamide (MSTFA) + fatty acid methyl esters (FAMEs) was added to each sample and further shaken at 37 °C for 0.5 h to complete derivatization. Then, 0.5 μL derivatized samples were injected on a 7890A gas chromatogram (GC) coupled with a time of flight mass spectrometer (TOF; LECO Corporation, St. Joseph, MI, USA) using a splitless method onto a RESTEK RTX-5SIL MS column with an Intergra-Guard at 275 °C with a helium flow of 1 mL.min−1. The GC oven was set at 50 °C for 1 min and then ramped to 330 °C at the rate of 20 °C.min−1 and held for 5 min. The transfer line was set to 280 °C and the EI ion source was set to 250 °C. The MS parameters collect data from 85 m/z to 500 m/z at an acquisition rate of 17 spectra/s.
Lipids were determined using the charged surface hybrid column electrospray method using a quadrupole time of flight mass spectrometer and tandem mass spectrometry (CSH-ESI QTOF MS/MS) as described previously [92,93]. In brief, 10 mg brewed coffee samples were vortexed with LCMS grade methanol (225 μL) and methyl tert-butyl ether (MTBE, 750 μL) and extracted by shaking for 6 min at 4 °C as previously described [90,93]. Samples were then vortexed with LCMS-grade water (188 μL). After centrifugation at 14,000 RCF (relative centrifugal force) for 2 min, the polar and non-polar layers were separated, dried, and stored until further analysis. Free fatty acids (FFAs); mono-, di-, and triglycerides (TGs); cholesteryl ester (CE); phospholipids (PLs); and sphingolipids (SLs) were analyzed by CSH-ESI QTOF MS/MS. Data were collected in both positive electrospray ionization (ESI) and negative ESI mode. Peaks were annotated by comparing MS/MS spectra and accurate masses of the precursor ion to spectra provided in the Fiehn Laboratory’s LipidBlast spectral library [94,95]. To profile biogenic amines, samples were analyzed using the hydrophilic interaction chromatography electrospray method, using a quadrupole time of flight mass spectrometer and tandem mass spectrometry (HILIC-ESI QTOF MS/MS) as described previously [90,93]. Data were collected in both positive and negative ESI modes and processed using MS-DIAL and MS-FLO programs as described [96].

2.8. Mineral and Metal Analysis of Brewed Coffee

The brewed coffees were lyophilized and sent to the Agricultural Diagnostic Service Center (ADSC), College of Tropical Agriculture and Human Resources (CTAHR), for analysis. The minerals analyzed included boron (B), calcium (Ca), copper (Cu), iron (Fe), magnesium (Mg), manganese (Mn), molybdenum (Mo), phosphorus (P), potassium (K), sodium (Na), and Zinc (Zn). The metals that were analyzed included arsenic (As), cadmium (Cd), chromium (Cr), nickel (Ni), lead (Pb), selenium (Se), and vanadium (Va). Three independent brews for each type of coffee and each brewing method were analyzed in triplicate.
Minerals and metals were analyzed by the standard EPA 3050B method (https://www.epa.gov/sites/default/files/2015-06/documents/epa-3050b.pdf, accessed on 5 January 2023) with slight modifications. In brief, 0.5 g lyophilized brewed coffee samples were used for analysis. Samples were digested with 3. 5 mL of concentrated nitric acid (15.8 N) for about 10 min at 110 °C. The samples were then mixed with 100 mL of ddH2O, incubated on a shaker, filtered through Whatman 42 filter paper, and analyzed on an Avio 560 Max Inductively Coupled Plasma Optical Emission Spectrometer (ICP-OES, Perkin Elmer, Waltham, MA, USA).

2.9. Antioxidant Capacity of Brewed Coffee

The antioxidant activity of brews prepared from roasted coffee grounds was analyzed using a commercial ORAC Antioxidant Assay Kit, Cat# AOX-2-RB(Zen-Bio, Inc., Research Triangle Park, North Carolina). The ORAC is a kinetic assay that measures fluorescein decay and antioxidant protection over time. It is an assay that is regularly used to measure the total antioxidant capacity of biological fluids, cells, and tissues as well as dietary supplements, therapeutics [97], and food extracts [98]. In this assay, the substrate 2,2′-azobis-2-methyl-propanimidamine dihydrochloride (AAPH) generates a peroxyl radical (ROO-) which is formed upon thermal homolysis at 37 °C. This peroxyl radical oxidizes fluorescein, which thereby produces a non-fluorescent product. The principle of the assay relies on the hydrogen atom transfer (HAT) mechanism of antioxidants to suppress oxidative degradation of the signal and measures the subsequent fluorescent activity vs. reaction time [99]. ORAC measures the inhibition of the peroxyl radical-induced oxidation as carried out by antioxidant compounds in coffee. The fluorescent measurements were expressed relative to the initial reading (excitation/emission at 485 nm/528–538 nm). Three independent brews for each type of coffee and each brewing method were analyzed in triplicate at four different concentrations of 1:100, 1:200, 1:500, and 1:750. Time course reactions of fluorescein decay from AAPH were measured using a 96-well plate reader (Perkin Elmer Wallac 1420 Victor2) and plate reading software (Wallac 1420).
Raw data from the ORAC assays were exported from the Wallac 1420 software to Excel and were normalized according to the following equation: AUC = 0.5 + (F1/F0) + (F2/F0) + … + 0.5 (F30/F0), where F0 = initial fluorescence reading at time = 0 min and F1, F2, etc., are the fluorescence readings at different time points. The area under the curve (AUC) values were generated using GraphPad Prism 7.0 for Mac OS X, version 7.0e. The net AUC was determined by subtracting the AUC of the blank from that of the coffee sample. A Trolox standard curve was generated using the net AUC vs. µM Trolox, in accordance with the one site-specific binding model (GraphPad Prism 7, GraphPad Software, CA, USA). Trolox equivalencies were determined by the one-site specific equation (Michaelis–Menten), as determined by the following equation: Y= Bmax * X/(Kd + X), where X is the equivalent Trolox concentration for coffee, Y is the net AUC for each coffee dilution, and Vmax and Kd are values from the Trolox standard curve. The ORAC data are represented as Trolox equivalences (TEs) of coffee. Three independent experiments were performed in replicates of six (n = 18).

2.10. Statistical Data Analysis

All coffee metabolites, lipids, and biogenic amines (known and unknown) were analyzed by univariate analysis using the parametric t-test to compare the mean differences of two groups. The one-way analysis of variance (ANOVA) was followed by Holm–Sidak correction for multiple comparisons. The known metabolites were further subjected to multivariate analysis using MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/docs/Publications.xhtml, accessed between 6 September and 20 October 2022) [31]. Data were normalized to sample median, log10 transformed and Pareto scaled. Principal component analysis (PCA) was used to observe clustering trends and exclude outliers. A discriminant model was created by further subjecting the data to partial least squares discriminant analysis (PLS-DA). One thousand permutation tests were performed to check model validity and potential over-fitting of the PLS-DA model and visualized using a validation plot. Known metabolites, lipids, and biogenic amines in each category were used to build the PLS-DA models. After validation of the PLS-DA model, data were further analyzed by orthogonal partial least squares discriminant analysis (OPLS-DA) to identify discriminant metabolites, lipids, and biogenic amines and distinguish between categories at false discovery rate (FDR) < 0.05.
Data for mineral and metal contents and ORAC assays are expressed as a mean ± standard deviation of triplicate or quadruplet values. Statistical significance was assessed using GraphPad Prism 7.0. Data were analyzed using either a two-tailed unpaired t-test or one-way ANOVA followed by Tukey’s test. p values < 0.05 were considered to be significant.

3. Results

A total of 442 metabolites tentatively assigned as global metabolites (139 known and 303 unknown), 1617 metabolites tentatively assigned as positive ESI lipids (80 known and 1537 unknown), 2862 metabolites tentatively assigned as negative ESI lipids (40 known and 2822 unknown), and 1747 metabolites tentatively assigned as biogenic amines (47 known and 1700 unknown) were detected in both coffee varieties.

3.1. Differences in Global Metabolites, Lipids, and Biogenic Amines among the Green and Roasted Coffee Varieties of “Kona Typica” (Waialua) and “Yellow Catuai” (Kaua’i) Brews

Univariate analysis indicated that 92 metabolites tentatively assigned as global metabolites (36 known and 54 unknown) were significantly different between the green varieties of Kaua’i and Waialua coffee (p < 0.05, supplement file 1). Similarly, 32 metabolites tentatively assigned as lipids (positive ESI; 4 known and 28 unknown), 103 metabolites tentatively assigned as negative ESI lipids (10 known and 93 unknown), and 132 metabolites tentatively assigned as biogenic amines (9 known and 123 unknown) were significantly different between the green varieties of Kaua’i and Waialua coffee (p < 0.05, supplement file 1).
Kaua’i green coffee brews contained higher amounts of 2-hydroxybutanoic acid (14.64-fold), LPC (18:0; 37.97-fold), chlorogenic acid (18:34-fold), and N6,N6,N6-Trimethyllysine (12.53-fold) as compared to Waialua green coffee brews (p < 0.05, supplement file 1). Orthogonal projections to latent structures discriminant analysis (OPLS-DA) score plots for green varieties of Kaua’i and Waialua coffee are represented in Figure 1 as determined by multivariate analysis.
Brews of the roasted coffees also demonstrated differences in 134 metabolites tentatively assigned as global metabolites (44 known and 89 unknown), 145 metabolites tentatively assigned positive ESI lipids (5 known and 140 unknown), 146 metabolites tentatively assigned as negative ESI lipids (12 known and 134 unknown), and 433 metabolites tentatively assigned as biogenic amines (15 known and 418 unknown) which were significantly different among Kaua’i and Waialua coffees (p < 0.05, supplement file 2).
Among the known compounds, galactinol, pyrogallol, uracil, phosphatidyl choline (PC) (34:1), PC (34:2), PC (36:2), PC (16:0), PE (38:2), PE (36:1), PC (38:2), PC (36:4), PC (36:2), PC (36:1), PC (35:4), PC (34:1), and DL-Indole-3-lactic acid were more than 2-fold higher in roasted Kaua’i coffee brews as compared to the roasted Waialua coffee brews (p < 0.05, supplement file 2). The OPLS-DA score plots for roasted Kaua’i and Waialua coffee brews are represented in Figure 2.

3.2. Influence of Roasting on Metabolites, Lipids, and Biogenic Amines of “Kona Typica” (Waialua) and “Yellow Catuai” (Kaua’i) Coffee Brews

As expected, roasting significantly increased caffeine and reduced chlorogenic acid in both coffee varieties (Figure 3A and Figure 3B, respectively). Similar effects of roasting on caffeine and chlorogenic contents of coffee have been reported by others [100]. Roasted Kaua’i coffee contained significantly high amounts of pyrogallol (Figure 3C, p < 0.05) compared to Waialua coffee. Overall differences in Kaua’i and Waialua coffee varieties as well as the effect of medium roasting conditions are noted in Figure 3D–G which depict hierarchal clustering (heatmaps) of all identified and unidentified metabolites for the green and roasted coffee brews: (C) global metabolites, (D) positive ESI lipid profiles, (E) negative ESI lipid profiles, and (G) biogenic amines.
Roasted Kaua’i coffee brews demonstrated significant differences in 102 known and 209 unknown metabolites tentatively assigned as global metabolites, 119 metabolites tentatively assigned as positive ESI mode lipids (5 known and 114 unknown), 678 metabolites tentatively assigned as negative ESI mode lipids (8 known and 670 unknown), and 1132 metabolites tentatively assigned as biogenic amines (38 known and 1094 unknown) as compared to the green Kaua’i coffee brews (p < 0.05, supplement file 3). Figure 4 compares relative changes in selected metabolites, lipids, and biogenic amines of green and roasted Kaua’i coffee brews (adjusted p value (FDR) < 0.05).
Similarly, 312 metabolites tentatively assigned as global metabolites (103 known and 209 unknown), 105 metabolites tentatively assigned as positive ESI mode lipids (5 known and 100 unknown), 697 metabolites tentatively assigned as negative ESI mode lipids (7 known and 690 unknown), and 1129 metabolites tentatively assigned as biogenic amines (38 known and 1091 unknown) were significantly different in roasted Waialua coffee brews as compared to the green Waialua coffee brews (p < 0.05, supplement file 4). Figure 4 depicts comparative selected box plots of known metabolites (Figure 4A–E), positive ESI lipids (Figure 4F–J), negative ESI lipids (Figure 4K–O), and biogenic amines (Figure 4P–T) as analyzed using MetaboAnalyst 5.0 (supplement file 5).

3.3. Influence of Brewing Methods on Metabolites, Lipids, and Biogenic Amines of “Kona Typica” (Waialua) and “Yellow Catuai” (Kaua’i) Coffee Brews

A total of six batches were roasted from each coffee variety. Three roasted batches were pooled into one sample, and the procedure was repeated twice to understand inter-roasting variability. Each variety of roasted coffee was brewed twice by four different brewing methods (n = 2 for each method). Known metabolites, lipids, and biogenic amines for both coffee types were compared for each brewing method by volcano plots using GraphPad Prism 7.0. For example, the cold brew of roasted Kaua’i coffee was compared with the cold brew of roasted Waialua coffee. Unpaired data were analyzed using parametric t-tests. The threshold for p value comparison was set at p < 0.05 and corrected for multiple comparisons using the Holm–Sidak method. No significant differences were noted among both coffee varieties brewed by the same brewing method. Therefore, data from both coffee varieties were combined for each brewing method (n = 4) and were analyzed using MetaboAnalyst 5.0.
Twenty-two metabolites tentatively assigned as global metabolites (11 known and 11 unknown) were significantly different in the four types of coffee brews as determined by one-way ANOVA followed by Fisher’s least significant difference method (Fisher’s LSD, adjusted p value (FDR) p < 0.05; supplement file 6). Interestingly, brewing methods significantly affected levels of several lipids. About 891 (49 known and 842 unknown) metabolites tentatively assigned as positive ESI mode lipids and 1447 (12 known and 1435 unknown) metabolites tentatively assigned as negative ESI mode lipids were significantly different in each type of brew (p < 0.05, supplement file 6). Selected significant metabolites (Figure 5A–D), positive ESI mode lipids (Figure 5E–H), and negative ESI mode lipids (Figure 5I–L) are depicted in Figure 5 below. Among the metabolites tentatively assigned as biogenic amines, only nine unknown amines were significantly different in the four types of brews (p < 0.05, supplement file 6).

3.4. Mineral and Metal Analysis of “Kona Typica” (Waialua) and “Yellow Catuai” (Kaua’i) Coffee Brews

Brews of roasted coffee grounds are most widely consumed worldwide as compared to brews of green beans. Hence, the mineral and metal profiles of only the roasted coffee bean brews were analyzed. As indicated in Table 1, the amounts of most minerals in cold brews were significantly higher compared to those brewed by the French press, filter paper, and filter mesh brewing methods (Table 1, p < 0.05). Although not statistically significant, total minerals were higher in W-CB > K-CB > K-FrP > W-FrP > W-FM > K-FM > W-FP > K-FP (Table 1).
An 8 oz cup of cold brew Kaua’i coffee contained 5.19%, 13.73%, 1.66%, 13.79%, 1.79%, 0.44%, 22.85%, 0.45%, and 3.15% of the RDA for P, K, Ca, Mg, Na, Fe, Mn, Zn, and Cu, respectively. All RDA values were calculated based on recommended values (National Academy of Sciences, 2019). The amounts of minerals in cold brew Kaua’i coffee were in the order of K > Mg > P > Na > Ca > Mn > B > Fe > Zn > Cu. For all other brews of Kaua’i coffee, the order of mineral abundance was K > Mg > Na > P > Ca > Mn > B > Fe > Zn = Cu. In general, except for boron, Kaua’i coffee brewed by filter paper method had the lowest mineral contents, followed by filter mesh < French press < cold brew.
Mineral contents in brewed Waialua coffee also showed similar trends to those observed in Kaua’i coffee (Table 1). An 8 oz cup of cold brew Waialua coffee contained 6.03%, 17.16%, 1.99%, 15.29%, 1.99%, 0.75%, 37.05%, 0.74%, and 11.83% of the RDAs for P, K, Ca, Mg, Na, Fe, Mn, Zn, and Cu, respectively. For Waialua coffee, most of the mineral contents were similar in those brewed by filter paper and filter mesh except for K, while the amounts of minerals in French press coffee were less than those in cold brews.
The Food and Nutrition Board has not established an RDA or AI for boron, which also does not have a DV. Total median boron intakes from dietary supplements and foods are about 1.0 to 1.5 mg.day−1 for adults (NIH.gov, accessed on 10 November 2022).

3.5. Antioxidant Capacity of “Kona Typica” (Waialua) and “Yellow Catuai” (Kaua’i) Coffee Brews

The antioxidant capacity, as measured by ORAC values, was influenced by the methods of brewing rather than the coffee variety. The highest ORAC values were noted for the cold brewing method and were about 2.6- to 2.8-fold higher compared to the filter paper method (Figure 6, p < 0.05). For both coffee varieties, the antioxidant capacities for different brewing methods were in the order of cold brew > filter mesh > French press > filter paper (Figure 6). Overall, there were no significant differences between the ORAC values for Kaua’i and Waialua coffee brews prepared by the same brewing method (Figure 6, p < 0.05).

4. Discussion

To understand the effects of health benefits and/or disease risks of coffee consumption, several studies have identified the metabolite composition of commercially brewed coffees based on the brewing methods, roast levels, coffee varieties, and/or caffeine contents [31,33,36,49,51,52,53,54,55,56]. The goal of our study was to identify differences in metabolites, lipids, and biogenic amines between the two types of Hawai‘i-grown coffees and the effects of roasting and brewing methods. Previous studies have used coffee beans to identify varietal differences in metabolites, while we used coffee brews [41,101]. The metabolomic approach has been previously used to identify specialty coffees and characterize their quality and exotic coffee tastes. For example, higher levels of sucrose and lower levels of γ-aminobutyric acid (GABA), quinic acid, choline, acetic acid, and fatty acids were observed in specialty or high-grade green coffees [102]. Similarly, arachidic acid and stearic acid were identified as markers for the Bourbon genealogical group; myristic and linoleic acids and sucrose, for the exotic genotype coffees; and lauric, palmitoleic, and oleic acids, for the Timor Hybrid group [103]. Metabolomic profiling has been also used to identify differences in Philippine coffee to distinguish between Coffea arabica (Arabica) and Coffea canephora var. Robusta coffees [41], sensory values [104], geographic diversity of Indonesian arabica coffees [105], defective coffee seeds of Brazilian coffee [101], and degrees of coffee adulteration in civet coffee blends [39]; discriminate Arabica and Robusta blends [106,107,108]; compare caffeinated and decaffeinated coffee [37]; and compare fermented coffees [109]. In keeping with published studies, the two Hawai‘i-grown coffee varieties also demonstrate distinct profiles of several metabolites in green beans as well as roasted beans. As previously noted by others [33], roasting significantly reduced chlorogenic acid and increased caffeine in both coffee varieties. Among the several coffee polyphenols, chlorogenic acid is most affected by roasting levels [110]. Compared to Waialua coffee, Kaua’i coffee had higher concentrations of pyrogallol, which is known to inhibit cellular glutathione (GSH) and induce apoptosis (cell death) in human platelets [111].
Lipid content in coffee beans accounts for about 10–17% of their dry weight, and most of these lipids are triacylglycerols (TAGs) and small quantities of phospholipids (PLs). TAGs contain both saturated and unsaturated fatty acids (FAs), with unsaturated FAs being oleic (18:1(n-9)), linoleic (18:2(n-6)), and linolenic (18:3(n-3)) [112]. Coffee lipids are major contributors to organoleptic properties, quality, and formation and stabilization of coffee foam and emulsion and influence flavor and aroma, specifically in espresso coffees [113,114,115]. Lipids in coffee beans are influenced not only by growing conditions such as altitude, shade, and temperature [115,116,117], but also by genotype [113,116,117,118]. For example, Arabica coffees generally have lipid contents of 15%, while Robusta has about 10% lipids [113]. Palmitic (16:0), arachidic (20:0), stearic (18:0), and linolenic (18:3) acid contents are higher and oleic acid (18:1) content is lower in Arabica compared to Robusta [113,119,120]. Considering the fact that lipids play an important role in coffee bean development, coffee brew properties, and the effects of coffee on human health, few studies have focused on lipidomic profiles in coffee beans [112,113,121]). To our knowledge, our study is the first to investigate the effect of coffee varieties, roasting, and brewing on the lipidomic profiles of coffee brews. In our study, palmitic (16:0), stearic (18:0), and oleic (18:1) acids were higher in the Kaua’i compared to Waialua roasted coffee brews. No major differences in these FAs were noted in green coffee beans. Both green and roasted Waialua and Kaua’i coffee brews also differed in a few phospholipids (PC; 34:1, 36:2, 34:2), phosphatidyl ethanolamine (1-PE. 17:0/17:0), and several unannotated lipids that warrant further investigation into their identities.
Biogenic amines (BAs) in food and beverages are mainly formed due to the decarboxylation or amination of proteins and/or free amino acids via microbial or natural enzyme activity. Putrescine, spermidine, spermine, and serotonin are the most abundant BAs in coffee beans and coffee beverages, while cadaverine and tyramine are present in smaller amounts [122,123]. Processing methods of unripe coffee beans can affect the final levels of some BAs. Histamine, tryptamine, and cadaverine were detected in low-quality and defective coffee beans. The presence of BAs is an indicator of undesired microbial activity and at high concentrations can be toxic to humans [29,124,125]. BA concentrations have been studied to identify coffee origins, e.g., Asia, South America, and Africa [126].
Putrescine is the most abundant amine in both Robusta and Arabica coffee, followed by spermidine, spermine, and serotonin, while cadaverine and tyramine are generally present in smaller amounts [127,128,129,130]. Putrescine is used to discriminate coffee species, while tyramine is considered a chemical marker for Angolan robusta, and low levels of histamine are present in low-quality or immature coffee beans [129,131]. In our study, although brews of green beans or the roasted coffee of both varieties demonstrated differences in several biogenic amines, putrescine, spermidine, spermine, and serotonin, cadaverine and tyramine were not detected. The effect of roasting on BAs is still controversial since some studies report a reduction in BAs while other studies indicate high BA levels after roasting [29,122]. Overall, brewed coffee contains a very low level of BAs compared to green or roasted ground coffee beans [29].
Metabolites, lipids, and BA composition are also affected by methods of brewing coffee, as noted in 76 commercial coffees [36]. Studies have indicated that different brewing methods affect not only the chemical composition but also the aroma [36,63]. Metabolite variations are also influenced by temperature as in hot or cold brews [86], brewing time, or the size of ground solids [45,58]. In contrast to studies by Rothwell et al. [36], studies by Kim et al. [33] indicated that the composition of bioactive compounds was dependent upon roasting rather than brewing methods. Kim et al. [33] demonstrated significantly different levels of resveratrol, eugenol, ferulic acid, and vanillin between hot and cold brews, but only for dark roasted coffee. In our studies, we did not detect resveratrol, eugenol, or vanillin. However, ferulic acid was highest in our coffees brewed by filter mesh (hot brewing method) and French press (cold brewing temperatures) as compared to cold brew and filter paper (hot brew) methods. Differences between our study and that of Kim et al. [33] could be attributed to the type of coffee, roasting level, or temperatures of brewing. Similar to Kim et al.’s study [33], caffeine levels were unaffected by brewing methods in our study. Ferulic acid is known for its antioxidant, anti-inflammatory, and antimicrobial properties [33]. In our study, chlorogenic acid was the lowest in cold-brewed coffee compared to other brewing methods. Caffeine levels were unaffected by brewing methods, which is similar to the result noted by Kim et al. [33].
Among several biological properties, the health benefits of coffee are attributed to its antioxidant capacity [43,64,80,81,83,132,133,134,135,136]. A large variation in antioxidant capacity has been observed among commercially brewed coffees [135] as well as variety and origin of coffee [137], degree of roasting [138], type of roast and blend [139], and brewing methods [43,64,81,136,140]. Among the 12 varieties of Arabica and 1 variety of Robusta studied by Priftis et al. [132], roasting increased the antioxidant capacity in 1 coffee variety and reduced it in 5 varieties of Arabica and resulted in no difference in the other varieties. On the other hand, slower roasting speeds and dark roasts reduced antioxidant capacity, and lightly roasted coffee had more antioxidants [83,133]. Besides roasting conditions, the preparation of coffee brews with different coffeemakers also influenced their antioxidant capacity. Mocha coffee had the highest antioxidant capacity compared to filter, plunger, and espresso coffees [134]. Another study indicated that cappuccino, a milk-based coffee drink, had the highest antioxidant activity as compared to instant coffee and Turkish coffee [141]. Wolska et al. [80] further demonstrated that brewing methods did not affect the antioxidant capacity of Robusta coffee brews, but simple infusion had the highest activity in Arabica brews prepared by French press, espresso maker, overflow espresso, and Turkish coffee. In contrast to the results of Wolska et al. [80], antioxidant capacity was unaffected by coffee types in our study but was influenced by brewing methods. Similar to the study by Perez-Martinez et al. [134], we also demonstrate that filter coffee had lower antioxidant capacity compared to other brews.
Studies have also indicated that besides antioxidant capacity, brewing methods also influence the aroma and mineral contents of coffee [64,142]. To date, one study has demonstrated that among the five brews from a coffee shop in Szczecin, Poland, coffee prepared by espresso machine had the lowest antioxidant capacity, followed by French press < drip = simple infusion coffee < Aeropress [64]. Furthermore, coffees prepared by simple infusion and Aeropress had higher levels of magnesium, manganese, chromium, cobalt, and potassium, while the drip brew had higher silicon levels [64]. Similar to Janda et al. [64], our study also indicated that French press coffee had a lower antioxidant capacity. We also observed an effect of brewing methods on the mineral contents of the coffee. Although the results are incomparable between our study and that of Janda et al. [64] due to differences in brewing methods, the individual mineral values in French press-brewed coffees are different in the two studies possibly due to the coffee variety or degree of roasting.

5. Conclusions

Our studies demonstrate that levels of global metabolites, lipids, and biogenic amines were significantly influenced by roasting and differed between the two Hawai‘i-grown coffee varieties: “Kona Typica” and “Yellow Catuai”. Interestingly, the mineral contents and the antioxidant capacity of both these coffee varieties were influenced only by the brewing methods. Total minerals were higher in cold brews compared to other brewing methods. Similarly, regardless of the coffee variety, cold brew coffees had the highest antioxidant capacity, followed by coffees brewed by the French press method. Future studies are warranted to understand and extrapolate the influence of brewing methods on antioxidant capacity in vivo.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo13030412/s1.

Author Contributions

Project conceptualization, experimental design, data analysis (MetaboAnalyst), data interpretation, manuscript writing, and funding, P.V.N.; partial manuscript writing, analyzing mineral and antioxidant data, and reviewing manuscript, J.Y.; roasting, grinding, brewing, lyophilizing coffee brews, and reviewing manuscript, K.I.; grinding, brewing, lyophilizing coffee brews, and reviewing manuscript, S.K.; brewing coffee, lyophilizing coffee brews, conducting antioxidant ORAC assays, and reviewing manuscript, J.W.; purchasing green coffee, roasting coffee, grinding coffee, editing and reviewing manuscript, H.C.B.; conducting univariate analysis, editing and reviewing manuscript, Y.D. All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by the Agriculture and Food Research Initiative (AFRI), National Institute of Food and Agriculture (NIFA; HAW05023-R, HAW00598-H, W3122, HAW00526-H, 2004-34135-15182), to P.V.N. The project was also supported in part by the Undergraduate Research Opportunities Program (UROP) to J.Y., K.I. and S.K. and by National Institutes of Health (NIH) grants P30GM114737, P20GM103466, U54MD007601, P20GM139753, and U54GM138062 to Y.D.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data supporting reported results can be found in this article and supplementary materials.

Acknowledgments

We would like to acknowledge Shawn Steiman (Coffea Consulting LLC, Honolulu, HI, USA) for helping us with the coffee roasting. We would also like to acknowledge the assistance of Rui Fang for partial univariate data analysis under the guidance of Youping Deng.

Conflicts of Interest

All authors affirm no conflict of interest. The funders had no role in study design; data collection, analyses, or interpretation; writing of the manuscript; or the decision to publish the results.

References

  1. News, L.A. Global Coffee Production Expected to Recover in 2022/23 with Brazil Crop. Available online: https://www.riotimesonline.com/brazil-news/rio-business/global-coffee-production-expected-to-recover-in-2022-23-with-brazil-crop/ (accessed on 5 January 2023).
  2. Halstead, T. Coffee: World Markets and Trade. Worls Production, Markets and Trade Report. Available online: https://apps.fas.usda.gov/psdonline/circulars/coffee.pdf (accessed on 5 January 2023).
  3. Gunter, M. 22 Coffee Consumption Statistics and Facts in 2022: Fascinating Data! Available online: https://coffeeaffection.com/coffee-consumption-statistics/ (accessed on 5 January 2023).
  4. Nerurkar, P.V.; Gandhi, K.; Chen, J.J. Correlations between Coffee Consumption and Metabolic Phenotypes, Plasma Folate, and Vitamin B12: NHANES 2003 to 2006. Nutrients 2021, 13, 1348. [Google Scholar] [CrossRef] [PubMed]
  5. Poole, R.; Kennedy, O.J.; Roderick, P.; Fallowfield, J.A.; Hayes, P.C.; Parkes, J. Coffee consumption and health: Umbrella review of meta-analyses of multiple health outcomes. BMJ 2017, 359, j5024. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Pourshahidi, L.K.; Navarini, L.; Petracco, M.; Strain, J.J. A Comprehensive Overview of the Risks and Benefits of Coffee Consumption. Compr. Rev. Food Sci. Food Saf. 2016, 15, 671–684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Santos, R.M.; Lima, D.R. Coffee consumption, obesity and type 2 diabetes: A mini-review. Eur. J. Nutr. 2016, 55, 1345–1358. [Google Scholar] [CrossRef]
  8. Cho, H.J.; Yoo, J.Y.; Kim, A.N.; Moon, S.; Choi, J.; Kim, I.; Ko, K.P.; Lee, J.E.; Park, S.K. Association of coffee drinking with all-cause and cause-specific mortality in over 190,000 individuals: Data from two prospective studies. Int. J. Food Sci. Nutr. 2022, 73, 513–521. [Google Scholar] [CrossRef]
  9. Liu, D.; Li, Z.H.; Shen, D.; Zhang, P.D.; Song, W.Q.; Zhang, W.T.; Huang, Q.M.; Chen, P.L.; Zhang, X.R.; Mao, C. Association of Sugar-Sweetened, Artificially Sweetened, and Unsweetened Coffee Consumption with All-Cause and Cause-Specific Mortality: A Large Prospective Cohort Study. Ann. Intern. Med. 2022, 175, 909–917. [Google Scholar] [CrossRef]
  10. Surma, S.; Banach, M. Coffee and caffeine consumption and overall mortality. Pleasure with restrictions-where do we really stand in 2022? Nutrition 2022, 102, 111747. [Google Scholar] [CrossRef]
  11. Zheng, H.; Lin, F.; Xin, N.; Yang, L.; Zhu, P. Association of Coffee, Tea, and Caffeine Consumption with All-Cause Risk and Specific Mortality for Cardiovascular Disease Patients. Front. Nutr. 2022, 9, 842856. [Google Scholar] [CrossRef]
  12. Nordestgaard, A.T. Causal relationship from coffee consumption to diseases and mortality: A review of observational and Mendelian randomization studies including cardiometabolic diseases, cancer, gallstones and other diseases. Eur. J. Nutr. 2022, 61, 573–587. [Google Scholar] [CrossRef]
  13. Sakamaki, T.; Kayaba, K.; Kotani, K.; Namekawa, M.; Hamaguchi, T.; Nakaya, N.; Ishikawa, S. Coffee consumption and mortality in Japan with 18 years of follow-up: The Jichi Medical School Cohort Study. Public Health 2021, 191, 23–30. [Google Scholar] [CrossRef]
  14. Curb, J.D.; Reed, D.M.; Kautz, J.A.; Yano, K. Coffee, caffeine, and serum cholesterol in Japanese men in Hawaii. Am. J. Epidemiol. 1986, 123, 648–655. [Google Scholar] [CrossRef] [PubMed]
  15. Davis, B.R.; Curb, J.D.; Borhani, N.O.; Prineas, R.J.; Molteni, A. Coffee consumption and serum cholesterol in the hypertension detection and follow-up program. Am. J. Epidemiol. 1988, 128, 124–136. [Google Scholar] [CrossRef] [PubMed]
  16. Hallstrom, H.; Byberg, L.; Glynn, A.; Lemming, E.W.; Wolk, A.; Michaelsson, K. Long-term coffee consumption in relation to fracture risk and bone mineral density in women. Am. J. Epidemiol. 2013, 178, 898–909. [Google Scholar] [CrossRef] [PubMed]
  17. Kujawska, A.; Kujawski, S.; Hajec, W.; Skierkowska, N.; Kwiatkowska, M.; Husejko, J.; Newton, J.L.; Simoes, J.A.; Zalewski, P.; Kedziora-Kornatowska, K. Coffee Consumption and Blood Pressure: Results of the Second Wave of the Cognition of Older People, Education, Recreational Activities, Nutrition, Comorbidities, and Functional Capacity Studies (COPERNICUS). Nutrients 2021, 13, 3372. [Google Scholar] [CrossRef]
  18. Greenop, K.R.; Miller, M.; Attia, J.; Ashton, L.J.; Cohn, R.; Armstrong, B.K.; Milne, E. Maternal consumption of coffee and tea during pregnancy and risk of childhood brain tumors: Results from an Australian case-control study. Cancer Causes Control. 2014, 25, 1321–1327. [Google Scholar] [CrossRef]
  19. Cheng, J.; Su, H.; Zhu, R.; Wang, X.; Peng, M.; Song, J.; Fan, D. Maternal coffee consumption during pregnancy and risk of childhood acute leukemia: A metaanalysis. Am. J. Obstet. Gynecol. 2014, 210, 151.e1–151.e10. [Google Scholar] [CrossRef]
  20. Bonaventure, A.; Rudant, J.; Goujon-Bellec, S.; Orsi, L.; Leverger, G.; Baruchel, A.; Bertrand, Y.; Nelken, B.; Pasquet, M.; Michel, G.; et al. Childhood acute leukemia, maternal beverage intake during pregnancy, and metabolic polymorphisms. Cancer Causes Control. 2013, 24, 783–793. [Google Scholar] [CrossRef]
  21. Thomopoulos, T.P.; Ntouvelis, E.; Diamantaras, A.A.; Tzanoudaki, M.; Baka, M.; Hatzipantelis, E.; Kourti, M.; Polychronopoulou, S.; Sidi, V.; Stiakaki, E.; et al. Maternal and childhood consumption of coffee, tea and cola beverages in association with childhood leukemia: A meta-analysis. Cancer Epidemiol. 2015, 39, 1047–1059. [Google Scholar] [CrossRef]
  22. Virtanen, S.M.; Uusitalo, L.; Kenward, M.G.; Nevalainen, J.; Uusitalo, U.; Kronberg-Kippila, C.; Ovaskainen, M.L.; Arkkola, T.; Niinisto, S.; Hakulinen, T.; et al. Maternal food consumption during pregnancy and risk of advanced beta-cell autoimmunity in the offspring. Pediatr. Diabetes 2011, 12, 95–99. [Google Scholar] [CrossRef]
  23. Pei, L.; Kang, Y.; Cheng, Y.; Yan, H. The Association of Maternal Lifestyle with Birth Defects in Shaanxi Province, Northwest China. PLoS ONE 2015, 10, e0139452. [Google Scholar] [CrossRef] [Green Version]
  24. Parazzini, F.; Chiaffarino, F.; Chatenoud, L.; Tozzi, L.; Cipriani, S.; Chiantera, V.; Fedele, L. Maternal coffee drinking in pregnancy and risk of small for gestational age birth. Eur. J. Clin. Nutr. 2005, 59, 299–301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Munoz, L.M.; Lonnerdal, B.; Keen, C.L.; Dewey, K.G. Coffee consumption as a factor in iron deficiency anemia among pregnant women and their infants in Costa Rica. Am. J. Clin. Nutr. 1998, 48, 645–651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Bech, B.H.; Nohr, E.A.; Vaeth, M.; Henriksen, T.B.; Olsen, J. Coffee and fetal death: A cohort study with prospective data. Am. J. Epidemiol. 2005, 162, 983–990. [Google Scholar] [CrossRef] [PubMed]
  27. Wisborg, K.; Kesmodel, U.; Bech, B.H.; Hedegaard, M.; Henriksen, T.B. Maternal consumption of coffee during pregnancy and stillbirth and infant death in first year of life: Prospective study. BMJ 2003, 326, 420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Leviton, A. Biases Inherent in Studies of Coffee Consumption in Early Pregnancy and the Risks of Subsequent Events. Nutrients 2018, 10, 1152. [Google Scholar] [CrossRef] [Green Version]
  29. Restuccia, D.; Spizzirri, U.G.; Parisi, O.I.; Cirillo, G.; Picci, N. Brewing effect on levels of biogenic amines in different coffee samples as determined by LC-UV. Food Chem. 2015, 175, 143–150. [Google Scholar] [CrossRef]
  30. Zhang, S.J.; De Bruyn, F.; Pothakos, V.; Torres, J.; Falconi, C.; Moccand, C.; Weckx, S.; De Vuyst, L. Following Coffee Production from Cherries to Cup: Microbiological and Metabolomic Analysis of Wet Processing of Coffea arabica. Appl. Environ. Microbiol. 2019, 85, e02635-18. [Google Scholar] [CrossRef] [Green Version]
  31. Chan, M.Z.A.; Lau, H.; Lim, S.Y.; Li, S.F.Y.; Liu, S.Q. Untargeted LC-QTOF-MS/MS based metabolomics approach for revealing bioactive components in probiotic fermented coffee brews. Food Res. Int. 2021, 149, 110656. [Google Scholar] [CrossRef]
  32. Hu, G.L.; Wang, X.; Zhang, L.; Qiu, M.H. The sources and mechanisms of bioactive ingredients in coffee. Food Funct. 2019, 10, 3113–3126. [Google Scholar] [CrossRef]
  33. Kim, C.H.; Park, S.J.; Yu, J.S.; Lee, D.Y. Interactive effect of post-harvest processing method, roasting degree, and brewing method on coffee metabolite profiles. Food Chem. 2022, 397, 133749. [Google Scholar] [CrossRef]
  34. Nosal, B.M.; Sakaki, J.R.; Kim, D.O.; Chun, O.K. Impact of coffee preparation on total phenolic content in brewed coffee extracts and their contribution to the body’s antioxidant status. Food Sci. Biotechnol. 2022, 31, 1081–1088. [Google Scholar] [CrossRef]
  35. Ciaramelli, C.; Palmioli, A.; Airoldi, C. Coffee variety, origin and extraction procedure: Implications for coffee beneficial effects on human health. Food Chem. 2019, 278, 47–55. [Google Scholar] [CrossRef]
  36. Rothwell, J.A.; Loftfield, E.; Wedekind, R.; Freedman, N.; Kambanis, C.; Scalbert, A.; Sinha, R. A Metabolomic Study of the Variability of the Chemical Composition of Commonly Consumed Coffee Brews. Metabolites 2019, 9, 17. [Google Scholar] [CrossRef]
  37. Chang, K.L.; Ho, P.C. Gas chromatography time-of-flight mass spectrometry (GC-TOF-MS)-based metabolomics for comparison of caffeinated and decaffeinated coffee and its implications for Alzheimer’s disease. PLoS ONE 2014, 9, e104621. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Amalia, F.; Aditiawati, P.; Yusianto; Putri, S.P.; Fukusaki, E. Gas chromatography/mass spectrometry-based metabolite profiling of coffee beans obtained from different altitudes and origins with various postharvest processing. Metabolomics 2021, 17, 69. [Google Scholar] [CrossRef] [PubMed]
  39. Jumhawan, U.; Putri, S.P.; Yusianto; Bamba, T.; Fukusaki, E. Quantification of coffee blends for authentication of Asian palm civet coffee (Kopi Luwak) via metabolomics: A proof of concept. J. Biosci. Bioeng. 2016, 122, 79–84. [Google Scholar] [CrossRef] [PubMed]
  40. Putri, S.P.; Irifune, T.; Yusianto; Fukusaki, E. GC/MS based metabolite profiling of Indonesian specialty coffee from different species and geographical origin. Metabolomics 2019, 15, 126. [Google Scholar] [CrossRef] [PubMed]
  41. Ongo, E.A.; Montevecchi, G.; Antonelli, A.; Sberveglieri, V.; Sevilla Iii, F. Metabolomics fingerprint of Philippine coffee by SPME-GC-MS for geographical and varietal classification. Food Res. Int. 2020, 134, 109227. [Google Scholar] [CrossRef] [PubMed]
  42. Bertrand, B.; Boulanger, R.; Dussert, S.; Ribeyre, F.; Berthiot, L.; Descroix, F.; Joet, 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]
  43. Ali, A.; Zahid, H.F.; Cottrell, J.J.; Dunshea, F.R. A Comparative Study for Nutritional and Phytochemical Profiling of Coffea arabica (C. arabica) from Different Origins and Their Antioxidant Potential and Molecular Docking. Molecules 2022, 27, 5126. [Google Scholar] [CrossRef]
  44. Villarreal, D.; Laffargue, A.; Posada, H.; Bertrand, B.; Lashermes, P.; Dussert, S. Genotypic and environmental effects on coffee (Coffea arabica L.) bean fatty acid profile: Impact on variety and origin chemometric determination. J. Agric. Food Chem. 2009, 57, 11321–11327. [Google Scholar] [CrossRef] [PubMed]
  45. Farag, M.A.; Zayed, A.; Sallam, I.E.; Abdelwareth, A.; Wessjohann, L.A. Metabolomics-Based Approach for Coffee Beverage Improvement in the Context of Processing, Brewing Methods, and Quality Attributes. Foods 2022, 11, 864. [Google Scholar] [CrossRef] [PubMed]
  46. Tassew, A.A.; Yadessa, G.B.; Bote, A.D.; Obso, T.K. Influence of location, elevation gradients, processing methods, and soil quality on the physical and cup quality of coffee in the Kafa Biosphere Reserve of SW Ethiopia. Heliyon 2021, 7, e07790. [Google Scholar] [CrossRef] [PubMed]
  47. van Mullem, J.J.; de Sousa Bueno Filho, J.S.; Dias, D.R.; Schwan, R.F. Chemical and sensory characterization of coffee from Coffea arabica cv. Mundo Novo and cv. Catuai Vermelho obtained by four different post-harvest processing methods. J. Sci. Food Agric. 2022, 102, 6687–6695. [Google Scholar] [CrossRef] [PubMed]
  48. Zhai, X.; Yang, M.; Zhang, J.; Zhang, L.; Tian, Y.; Li, C.; Bao, L.; Ma, C.; Abd El-Aty, A.M. Feasibility of Ultrasound-Assisted Extraction for Accelerated Cold Brew Coffee Processing: Characterization and Comparison with Conventional Brewing Methods. Front. Nutr. 2022, 9, 849811. [Google Scholar] [CrossRef] [PubMed]
  49. Niseteo, T.; Komes, D.; Belscak-Cvitanovic, A.; Horzic, D.; Budec, M. Bioactive composition and antioxidant potential of different commonly consumed coffee brews affected by their preparation technique and milk addition. Food Chem. 2012, 134, 1870–1877. [Google Scholar] [CrossRef] [PubMed]
  50. Villalon-Lopez, N.; Serrano-Contreras, J.I.; Tellez-Medina, D.I.; Gerardo Zepeda, L. An H NMR-based metabolomic approach to compare the chemical profiling of retail samples of ground roasted and instant coffees. Food Res. Int. 2018, 106, 263–270. [Google Scholar] [CrossRef]
  51. Perez-Miguez, R.; Castro-Puyana, M.; Sanchez-Lopez, E.; Plaza, M.; Marina, M.L. Untargeted HILIC-MS-Based Metabolomics Approach to Evaluate Coffee Roasting Process: Contributing to an Integrated Metabolomics Multiplatform. Molecules 2020, 25, 887. [Google Scholar] [CrossRef] [Green Version]
  52. Perez-Miguez, R.; Sanchez-Lopez, E.; Plaza, M.; Castro-Puyana, M.; Marina, M.L. A non-targeted metabolomic approach based on reversed-phase liquid chromatography-mass spectrometry to evaluate coffee roasting process. Anal. Bioanal. Chem. 2018, 410, 7859–7870. [Google Scholar] [CrossRef]
  53. Perez-Miguez, R.; Sanchez-Lopez, E.; Plaza, M.; Marina, M.L.; Castro-Puyana, M. Capillary electrophoresis-mass spectrometry metabolic fingerprinting of green and roasted coffee. J. Chromatogr. A 2019, 1605, 360353. [Google Scholar] [CrossRef]
  54. Angeloni, G.; Guerrini, L.; Masella, P.; Bellumori, M.; Daluiso, S.; Parenti, A.; Innocenti, M. What kind of coffee do you drink? An investigation on effects of eight different extraction methods. Food Res. Int. 2019, 116, 1327–1335. [Google Scholar] [CrossRef] [PubMed]
  55. Xu, L.; Lao, F.; Xu, Z.; Wang, X.; Chen, F.; Liao, X.; Chen, A.; Yang, S. Use of liquid chromatography quadrupole time-of-flight mass spectrometry and metabolomic approach to discriminate coffee brewed by different methods. Food Chem. 2019, 286, 106–112. [Google Scholar] [CrossRef] [PubMed]
  56. Zhang, J.; Sun, X.; Liu, P.; Zhang, T.; Jelderks, J.A.; Corke, H. Preliminary Characterization of Phytochemicals and Polysaccharides in Diverse Coffee Cascara Samples: Identification, Quantification and Discovery of Novel Compounds. Foods 2022, 11, 1710. [Google Scholar] [CrossRef] [PubMed]
  57. Aditiawati, P.; Astuti, D.I.; Kriswantoro, J.A.; Khanza, S.M.; Kamarisima; Irifune, T.; Amalia, F.; Fukusaki, E.; Putri, S.P. GC/MS-based metabolic profiling for the evaluation of solid state fermentation to improve quality of Arabica coffee beans. Metabolomics 2020, 16, 57. [Google Scholar] [CrossRef] [PubMed]
  58. Cordoba, N.; Pataquiva, L.; Osorio, C.; Moreno, F.L.M.; Ruiz, R.Y. Effect of grinding, extraction time and type of coffee on the physicochemical and flavour characteristics of cold brew coffee. Sci. Rep. 2019, 9, 8440. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Bressanello, D.; Liberto, E.; Cordero, C.; Rubiolo, P.; Pellegrino, G.; Ruosi, M.R.; Bicchi, C. Coffee aroma: Chemometric comparison of the chemical information provided by three different samplings combined with GC-MS to describe the sensory properties in cup. Food Chem. 2017, 214, 218–226. [Google Scholar] [CrossRef] [PubMed]
  60. Bressanello, D.; Liberto, E.; Cordero, C.; Sgorbini, B.; Rubiolo, P.; Pellegrino, G.; Ruosi, M.R.; Bicchi, C. Chemometric Modeling of Coffee Sensory Notes through Their Chemical Signatures: Potential and Limits in Defining an Analytical Tool for Quality Control. J. Agric. Food Chem. 2018, 66, 7096–7109. [Google Scholar] [CrossRef]
  61. Bressanello, D.; Marengo, A.; Cordero, C.; Strocchi, G.; Rubiolo, P.; Pellegrino, G.; Ruosi, M.R.; Bicchi, C.; Liberto, E. Chromatographic Fingerprinting Strategy to Delineate Chemical Patterns Correlated to Coffee Odor and Taste Attributes. J. Agric. Food Chem. 2021, 69, 4550–4560. [Google Scholar] [CrossRef]
  62. Liberto, E.; Bressanello, D.; Strocchi, G.; Cordero, C.; Ruosi, M.R.; Pellegrino, G.; Bicchi, C.; Sgorbini, B. HS-SPME-MS-Enose Coupled with Chemometrics as an Analytical Decision Maker to Predict In-Cup Coffee Sensory Quality in Routine Controls: Possibilities and Limits. Molecules 2019, 24, 4515. [Google Scholar] [CrossRef] [Green Version]
  63. Caprioli, G.; Cortese, M.; Sagratini, G.; Vittori, S. The influence of different types of preparation (espresso and brew) on coffee aroma and main bioactive constituents. Int. J. Food Sci. Nutr. 2015, 66, 505–513. [Google Scholar] [CrossRef] [PubMed]
  64. Janda, K.; Jakubczyk, K.; Baranowska-Bosiacka, I.; Kapczuk, P.; Kochman, J.; Rebacz-Maron, E.; Gutowska, I. Mineral Composition and Antioxidant Potential of Coffee Beverages Depending on the Brewing Method. Foods 2020, 9, 121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. El-Hawary, E.A.; Zayed, A.; Laub, A.; Modolo, L.V.; Wessjohann, L.; Farag, M.A. How Does LC/MS Compare to UV in Coffee Authentication and Determination of Antioxidant Effects? Brazilian and Middle Eastern Coffee as Case Studies. Antioxidants 2022, 11, 131. [Google Scholar] [CrossRef] [PubMed]
  66. Seremet, D.; Fabecic, P.; Vojvodic Cebin, A.; Mandura Jaric, A.; Pudic, R.; Komes, D. Antioxidant and Sensory Assessment of Innovative Coffee Blends of Reduced Caffeine Content. Molecules 2022, 27, 448. [Google Scholar] [CrossRef] [PubMed]
  67. Stanek, N.; Zarebska, M.; Bilos, L.; Barabosz, K.; Nowakowska-Bogdan, E.; Semeniuk, I.; Blaszkiewicz, J.; Kulesza, R.; Matejuk, R.; Szkutnik, K. Influence of coffee brewing methods on the chromatographic and spectroscopic profiles, antioxidant and sensory properties. Sci. Rep. 2021, 11, 21377. [Google Scholar] [CrossRef] [PubMed]
  68. Doo, T.; Morimoto, Y.; Steinbrecher, A.; Kolonel, L.N.; Maskarinec, G. Coffee intake and risk of type 2 diabetes: The Multiethnic Cohort. Public Health Nutr. 2014, 17, 1328–1336. [Google Scholar] [CrossRef] [Green Version]
  69. Gelber, R.P.; Petrovitch, H.; Masaki, K.H.; Ross, G.W.; White, L.R. Coffee intake in midlife and risk of dementia and its neuropathologic correlates. J. Alzheimer Dis. 2011, 23, 607–615. [Google Scholar] [CrossRef] [Green Version]
  70. Ross, G.W.; Abbott, R.D.; Petrovitch, H.; Morens, D.M.; Grandinetti, A.; Tung, K.H.; Tanner, C.M.; Masaki, K.H.; Blanchette, P.L.; Curb, J.D.; et al. Association of coffee and caffeine intake with the risk of Parkinson disease. JAMA 2000, 283, 2674–2679. [Google Scholar] [CrossRef]
  71. Setiawan, V.W.; Wilkens, L.R.; Lu, S.C.; Hernandez, B.Y.; Le Marchand, L.; Henderson, B.E. Association of coffee intake with reduced incidence of liver cancer and death from chronic liver disease in the US multiethnic cohort. Gastroenterology 2015, 148, 118–125; quiz e115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Hawaii. Available online: https://hfbf.org/mac-nuts/ (accessed on 5 January 2023).
  73. Top 20 Agricultural Commodities Produced State of Hawaii. Available online: https://hdoa.hawaii.gov/add/files/2021/12/Top-20-Commodities-2019_SOH_01.28.21.pdf (accessed on 5 January 2023).
  74. Kostenko, A. Coffee Industry Valued Higher Now. Available online: https://hawaiicoffeeassoc.org/News/10075301 (accessed on 5 January 2023).
  75. Sky, A. 33+ Coffee Industry Statistics [2022]: Cafes, Consumption, And Market Trends. Available online: https://www.zippia.com/advice/coffee-industry-statistics/ (accessed on 5 January 2023).
  76. A Guide to Hawaiian Coffee Production. Available online: https://perfectdailygrind.com/2022/01/a-guide-to-hawaiian-coffee-production/ (accessed on 5 January 2023).
  77. Grembecka, M.; Malinowska, E.; Szefer, P. Differentiation of market coffee and its infusions in view of their mineral composition. Sci. Total Environ. 2007, 383, 59–69. [Google Scholar] [CrossRef]
  78. Stelmach, E.; Pohl, P.; Szymczycha-Madeja, A. The suitability of the simplified method of the analysis of coffee infusions on the content of Ca, Cu, Fe, Mg, Mn and Zn and the study of the effect of preparation conditions on the leachability of elements into the coffee brew. Food Chem. 2013, 141, 1956–1961. [Google Scholar] [CrossRef]
  79. Stelmach, E.; Szymczycha-Madeja, A.; Pohl, P. A simplified determination of total concentrations of Ca, Fe, Mg and Mn in addition to their bioaccessible fraction in popular instant coffee brews. Food Chem. 2016, 197, 388–394. [Google Scholar] [CrossRef] [PubMed]
  80. Wolska, J.; Janda, K.; Jakubczyk, K.; Szymkowiak, M.; Chlubek, D.; Gutowska, I. Levels of Antioxidant Activity and Fluoride Content in Coffee Infusions of Arabica, Robusta and Green Coffee Beans in According to their Brewing Methods. Biol. Trace Elem. Res. 2017, 179, 327–333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Liang, N.; Kitts, D.D. Antioxidant property of coffee components: Assessment of methods that define mechanisms of action. Molecules 2014, 19, 19180–19208. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Di Girolamo, F.G.; Mazzucco, S.; Situlin, R.; Mohorko, N.; Jenko-Praznikar, Z.; Petelin, A.; Tence, M.; Pisot, R.; Navarini, L.; Biolo, G. Roasting intensity of naturally low-caffeine Laurina coffee modulates glucose metabolism and redox balance in humans. Nutrition 2016, 32, 928–936. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Jung, S.; Kim, M.H.; Park, J.H.; Jeong, Y.; Ko, K.S. Cellular Antioxidant and Anti-Inflammatory Effects of Coffee Extracts with Different Roasting Levels. J. Med. Food 2017, 20, 626–635. [Google Scholar] [CrossRef]
  84. Ho, K.V.; Schreiber, K.L.; Park, J.; Vo, P.H.; Lei, Z.; Sumner, L.W.; Brown, C.R.; Lin, C.H. Identification and Quantification of Bioactive Molecules Inhibiting Pro-inflammatory Cytokine Production in Spent Coffee Grounds Using Metabolomics Analyses. Front. Pharmacol. 2020, 11, 229. [Google Scholar] [CrossRef]
  85. Angeloni, G.; Guerrini, L.; Masella, P.; Innocenti, M.; Bellumori, M.; Parenti, A. Characterization and comparison of cold brew and cold drip coffee extraction methods. J. Sci. Food Agric. 2019, 99, 391–399. [Google Scholar] [CrossRef]
  86. Rao, N.Z.; Fuller, M.; Grim, M.D. Physiochemical Characteristics of Hot and Cold Brew Coffee Chemistry: The Effects of Roast Level and Brewing Temperature on Compound Extraction. Foods 2020, 9, 902. [Google Scholar] [CrossRef]
  87. Cajka, T.; Fiehn, O. Toward Merging Untargeted and Targeted Methods in Mass Spectrometry-Based Metabolomics and Lipidomics. Anal. Chem. 2016, 88, 524–545. [Google Scholar] [CrossRef]
  88. Fiehn, O. Metabolomics by Gas Chromatography-Mass Spectrometry: Combined Targeted and Untargeted Profiling. Curr. Protoc. Mol. Biol. 2016, 114, 30.4.1–30.4.32. [Google Scholar] [CrossRef] [Green Version]
  89. Sud, M.; Fahy, E.; Cotter, D.; Azam, K.; Vadivelu, I.; Burant, C.; Edison, A.; Fiehn, O.; Higashi, R.; Nair, K.S.; et al. Metabolomics Workbench: An international repository for metabolomics data and metadata, metabolite standards, protocols, tutorials and training, and analysis tools. Nucleic Acids Res. 2016, 44, D463–D470. [Google Scholar] [CrossRef] [Green Version]
  90. Johnson, C.M.; Rosario, R.; Brito, A.; Agrawal, K.; Fanter, R.; Lietz, G.; Haskell, M.; Engle-Stone, R.; Newman, J.W.; La Frano, M.R. Multiassay nutritional metabolomics profiling of low vitamin A status versus adequacy is characterized by reduced plasma lipid mediators among lactating women in the Philippines: A pilot study. Nutr. Res. 2022, 104, 118–127. [Google Scholar] [CrossRef] [PubMed]
  91. Granda, M.L.; Prince, D.K.; Fiehn, O.; Chen, Y.; Rajabi, T.; Yeung, C.K.; Hoofnagle, A.N.; Kestenbaum, B. Metabolomic Profiling Identifies New Endogenous Markers of Tubular Secretory Clearance. Kidney360 2023, 4, 23–31. [Google Scholar] [CrossRef] [PubMed]
  92. Miao, G.; Deen, J.; Struzeski, J.B.; Chen, M.; Zhang, Y.; Cole, S.A.; Fretts, A.M.; Lee, E.T.; Howard, B.V.; Fiehn, O.; et al. Plasma lipidomic profile of depressive symptoms: A longitudinal study in a large sample of community-dwelling American Indians in the strong heart study. Mol. Psychiatry 2023, 1–10. [Google Scholar] [CrossRef] [PubMed]
  93. Mazi, T.A.; Borkowski, K.; Newman, J.W.; Fiehn, O.; Bowlus, C.L.; Sarkar, S.; Matsukuma, K.; Ali, M.R.; Kieffer, D.A.; Wan, Y.Y.; et al. Ethnicity-specific alterations of plasma and hepatic lipidomic profiles are related to high NAFLD rate and severity in Hispanic Americans, a pilot study. Free. Radic. Biol. Med. 2021, 172, 490–502. [Google Scholar] [CrossRef]
  94. Kind, T.; Liu, K.H.; Lee, D.Y.; DeFelice, B.; Meissen, J.K.; Fiehn, O. LipidBlast in silico tandem mass spectrometry database for lipid identification. Nat. Methods 2013, 10, 755–758. [Google Scholar] [CrossRef] [Green Version]
  95. Kind, T.; Okazaki, Y.; Saito, K.; Fiehn, O. LipidBlast templates as flexible tools for creating new in-silico tandem mass spectral libraries. Anal. Chem. 2014, 86, 11024–11027. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Tsugawa, H.; Cajka, T.; Kind, T.; Ma, Y.; Higgins, B.; Ikeda, K.; Kanazawa, M.; VanderGheynst, J.; Fiehn, O.; Arita, M. MS-DIAL: Data-independent MS/MS deconvolution for comprehensive metabolome analysis. Nat. Methods 2015, 12, 523–526. [Google Scholar] [CrossRef] [PubMed]
  97. ORAC Antioxidant Assay Kit. Available online: http://www.zen-bio.com/products/kits/orac-antioxidant_screening.php (accessed on 29 May 2018).
  98. Pertuzatti, P.B.; Barcia, M.T.; Rodrigues, D.; da Cruz, P.N.; Hermosín-Gutiérrez, I.; Smith, R.; Godoy, H.T. Antioxidant activity of hydrophilic and lipophilic extracts of Brazilian blueberries. Food Chem. 2014, 164, 81–88. [Google Scholar] [CrossRef] [Green Version]
  99. Apak, R.; Özyürek, M.; Güçlü, K.; Çapanoğlu, E. Antioxidant Activity/Capacity Measurement. 2. Hydrogen Atom Transfer (HAT)-Based, Mixed-Mode (Electron Transfer (ET)/HAT), and Lipid Peroxidation Assays. J. Agric. Food Chem. 2016, 64, 1028–1045. [Google Scholar] [CrossRef]
  100. Pilipczuk, T.; Kusznierewicz, B.; Zielinska, D.; Bartoszek, A. The influence of roasting and additional processing on the content of bioactive components in special purpose coffees. J. Food Sci. Technol. 2015, 52, 5736–5744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Toci, A.T.; Farah, A. Volatile fingerprint of Brazilian defective coffee seeds: Corroboration of potential marker compounds and identification of new low quality indicators. Food Chem. 2014, 153, 298–314. [Google Scholar] [CrossRef] [PubMed]
  102. Kwon, D.-J.; Jeong, H.-J.; Moon, H.; Kim, H.-N.; Cho, J.-H.; Lee, J.-E.; Hong, K.S.; Hong, Y.-S. Assessment of green coffee beanmetabolites dependent on coffee quality using a 1H NMR-based metabolomics approach. Food Res. Int. 2015, 67, 175–182. [Google Scholar] [CrossRef]
  103. Malta, M.R.; Fassio, L.O.; Liska, G.R.; Carvalho, G.R.; Pereira, A.A.; Botelho, C.E.; Ferraz, V.P.; Silva, A.D.; Pedrosa, A.W.; Alvaro, L.N.; et al. Discrimination of genotypes coffee by chemical composition of the beans: Potential markers in natural coffees. Food Res. Int. 2020, 134, 109219. [Google Scholar] [CrossRef] [PubMed]
  104. Rocchetti, G.; Braceschi, G.P.; Odello, L.; Bertuzzi, T.; Trevisan, M.; Lucini, L. Identification of markers of sensory quality in ground coffee: An untargeted metabolomics approach. Metabolomics 2020, 16, 127. [Google Scholar] [CrossRef] [PubMed]
  105. Happyana, N.; Pratiwi, A.; Hakim, E.H. Metabolite Profiles of the Green Beans of Indonesian Arabica Coffee Varieties. Int. J. Food Sci. 2021, 2021, 5782578. [Google Scholar] [CrossRef]
  106. Wongsa, P.; Khampa, N.; Horadee, S.; Chaiwarith, J.; Rattanapanone, N. Quality and bioactive compounds of blends of Arabica and Robusta spray-dried coffee. Food Chem. 2019, 283, 579–587. [Google Scholar] [CrossRef]
  107. Colzi, I.; Taiti, C.; Marone, E.; Magnelli, S.; Gonnelli, C.; Mancuso, S. Covering the different steps of the coffee processing: Can headspace VOC emissions be exploited to successfully distinguish between Arabica and Robusta? Food Chem. 2017, 237, 257–263. [Google Scholar] [CrossRef]
  108. Monakhova, Y.B.; Ruge, W.; Kuballa, T.; Ilse, M.; Winkelmann, O.; Diehl, B.; Thomas, F.; Lachenmeier, D.W. Rapid approach to identify the presence of Arabica and Robusta species in coffee using 1H NMR spectroscopy. Food Chem. 2015, 182, 178–184. [Google Scholar] [CrossRef]
  109. Pereira, T.S.; Batista, N.N.; Santos Pimenta, L.P.; Martinez, S.J.; Ribeiro, L.S.; Oliveira Naves, J.A.; Schwan, R.F. Self-induced anaerobiosis coffee fermentation: Impact on microbial communities, chemical composition and sensory quality of coffee. Food Microbiol. 2022, 103, 103962. [Google Scholar] [CrossRef]
  110. Król, K.; Gantner, M.; Tatarak, A.; Hallmann, E. The content of polyphenols in coffee beans as roasting, origin and storage effect. Eur. Food Res. Technol. 2020, 246, 33–39. [Google Scholar] [CrossRef] [Green Version]
  111. Bruges, G.; Venturini, W.; Crespo, G.; Zambrano, M.L. Pyrogallol Induces Apoptosis in Human Platelets. Folia Biol. 2018, 64, 23–30. [Google Scholar]
  112. Kanai, N.; Yoshihara, N.; Kawamura, I. Solid-state NMR characterization of triacylglycerol and polysaccharides in coffee beans. Biosci. Biotechnol. Biochem. 2019, 83, 803–809. [Google Scholar] [CrossRef] [Green Version]
  113. Anagbogu, C.F.; Zhou, J.; Olasupo, F.O.; Baba Nitsa, M.; Beckles, D.M. Lipidomic and metabolomic profiles of Coffea canephora L. beans cultivated in Southwestern Nigeria. PLoS ONE 2021, 16, e0234758. [Google Scholar] [CrossRef] [PubMed]
  114. Figueiredo, L.P.; Borem, F.M.; Ribeiro, F.C.; Giomo, G.S.; Henrique, d.S.J.; Malta, M.R. Fatty acid profiles and parameters of quality of specialty coffees produced in different Brazilian regions. Afr. J. Agric. Res. 2015, 10, 3484–3493. [Google Scholar] [CrossRef]
  115. Martin, M.J.; Pablos, F.; Gonzalez, A.G.; Valdenebro, M.S.; Leon-Camacho, M. Fatty acid profiles as discriminant parameters for coffee varieties differentiation. Talanta 2001, 54, 291–297. [Google Scholar] [CrossRef]
  116. Joet, T.; Laffargue, A.; Descroix, F.; Doulbeau, S.; Bertrand, B.; de Kochko, A.; Dussert, S. Influence of environmental factors, wet processing and their interactions on the biochemical composition of green Arabica coffee beans. Food Chem. 2010, 118, 693–701. [Google Scholar] [CrossRef]
  117. Mehari, B.; Redi-Abshiro, M.; Chandravanshi, B.S.; Combrinck, S.; McCrindle, R.; Atlabachew, M. GC-MS profiling of fatty acids in green coffee (Coffea arabica L.) beans and chemometric modeling for tracing geographical origins from Ethiopia. J. Sci. Food Agric. 2019, 99, 3811–3823. [Google Scholar] [CrossRef] [PubMed]
  118. Silva, A.C.R.; da Silva, C.C.; Garrett, R.; Rezende, C.M. Comprehensive lipid analysis of green Arabica coffee beans by LC-HRMS/MS. Food Res. Int. 2020, 137, 109727–109736. [Google Scholar] [CrossRef]
  119. Cossignani, L.; Montesano, D.; Simonetti, M.S.; Blasi, F. Authentication of Coffea arabica according to Triacylglycerol Stereospecific Composition. J. Anal. Methods Chem. 2016, 2016, 7482620. [Google Scholar] [CrossRef] [Green Version]
  120. Romano, R.; Santini, A.; Le Grottaglie, L.; Manzo, N.; Visconti, A.; Ritien, i.A. Identification markers based on fatty acid composition to differentiate between roasted Arabica and Canephora (Robusta) coffee varieties in mixtures. J. Food Compos. Anal. 2014, 35, 1–9. [Google Scholar] [CrossRef] [Green Version]
  121. Silva, F.L.F.; Nascimento, G.O.; Lopes, G.S.; Matos, W.O.; Cunha, R.L.; Malta, M.R.; Liska, G.R.; Owen, R.W.; Trevisan, M.T.S. The concentration of polyphenolic compounds and trace elements in the Coffea arabica leaves: Potential chemometric pattern recognition of coffee leaf rust resistance. Food Res. Int. 2020, 134, 109221. [Google Scholar] [CrossRef] [PubMed]
  122. Visciano, P.; Schirone, M. Update on Biogenic Amines in Fermented and Non-Fermented Beverages. Foods 2022, 11, 353. [Google Scholar] [CrossRef] [PubMed]
  123. Papageorgiou, M.; Lambropoulou, D.; Morrison, C.; Kłodzínska, E.; Namièsnik, J.; Płotka-Wasylka, J. Literature update of analytical methods for biogenic amines determination in food and beverages. Trends Anal. Chem. 2018, 98, 128–142. [Google Scholar] [CrossRef] [Green Version]
  124. Kalac, P. Recent advances in the research on biological roles of dietary polyamines in man. J. Appl. Biomed. 2009, 7, 65–74. [Google Scholar] [CrossRef] [Green Version]
  125. Shalaby, A.R. Significance of biogenic amines to food safety and human health. Food Res. Int. 1996, 29, 675–690. [Google Scholar] [CrossRef]
  126. Choi, M.-Y.; Choi, W.; Park, J.H.; Lim, J.; Kwon, S.W. Determination of coffee origins by integrated metabolomic approach of combining multiple analytical data. Food Chem. 2010, 121, 1260–1268. [Google Scholar] [CrossRef]
  127. Amorim, H.V.; Basso, L.C.; Crocomo, O.J.; Teixeira, A.A. Polyamines in green and roasted coffee. J. Agric. Food Chem. 1977, 25, 957–958. [Google Scholar] [CrossRef]
  128. Casal, S.; Mendes, E.; Oliveira, M.B.P.P.; Ferreira, M.A. Determination of biogenic amines in coffee by an optimized liquid chromatographic method. J. Liq. Chromatogr. Relat. Technol. 2002, 25, 2535–2549. [Google Scholar] [CrossRef]
  129. Casal, S.; Mendes, E.; Alves, M.R.; Alves, R.C.; Beatriz, M.; Oliveira, P.P.; Ferreira, M.A. Free and conjugated biogenic amines in green and roasted coffee beans. J. Agric. Food Chem. 2004, 52, 6188–6192. [Google Scholar] [CrossRef]
  130. Cirilo, M.P.G.; Coelho, A.F.S.; Araujo, C.M.; Goncalves, F.R.B.; Nogueira, F.D.; Gloria, M.B.A. Profiles and levels of bioactive amines in green and roasted coffee. Food Chem. 2003, 82, 397–402. [Google Scholar] [CrossRef]
  131. Dias, E.C.; Pereira, R.G.; Borem, F.M.; Mendes, E.; de Lima, R.R.; Fernandes, J.O.; Casal, S. Biogenic amine profile in unripe Arabica coffee beans processed according to dry and wet methods. J. Agric. Food Chem. 2012, 60, 4120–4125. [Google Scholar] [CrossRef] [PubMed]
  132. Priftis, A.; Stagos, D.; Konstantinopoulos, K.; Tsitsimpikou, C.; Spandidos, D.A.; Tsatsakis, A.M.; Tzatzarakis, M.N.; Kouretas, D. Comparison of antioxidant activity between green and roasted coffee beans using molecular methods. Mol. Med. Rep. 2015, 12, 7293–7302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  133. Opitz, S.E.; Goodman, B.A.; Keller, M.; Smrke, S.; Wellinger, M.; Schenker, S.; Yeretzian, C. Understanding the Effects of Roasting on Antioxidant Components of Coffee Brews by Coupling On-line ABTS Assay to High Performance Size Exclusion Chromatography. Phytochem. Anal. 2017, 28, 106–114. [Google Scholar] [CrossRef] [PubMed]
  134. Perez-Martinez, M.; Caemmerer, B.; De Pena, M.P.; Cid, C.; Kroh, L.W. Influence of brewing method and acidity regulators on the antioxidant capacity of coffee brews. J. Agric. Food Chem. 2010, 58, 2958–2965. [Google Scholar] [CrossRef] [PubMed]
  135. Contreras-Calderon, J.; Mejia-Diaz, D.; Martinez-Castano, M.; Bedoya-Ramirez, D.; Lopez-Rojas, N.; Gomez-Narvaez, F.; Medina-Pineda, Y.; Vega-Castro, O. Evaluation of antioxidant capacity in coffees marketed in Colombia: Relationship with the extent of non-enzymatic browning. Food Chem. 2016, 209, 162–170. [Google Scholar] [CrossRef]
  136. Trandafir, I.; Nour, V.; Ionica, M.E. Antioxidant capacity, phenolic acids and caffeine contents of some commercial coffees available on the Romanian market. Arch. Latinoam. Nutr. 2013, 63, 87–94. [Google Scholar]
  137. Parras, P.; Martínez-Tomè, M.; Jimènez, A.M.; Murcia, M.A. Antioxidant capacity of coffees of several origins brewed following three different procedures. Food Chem. 2007, 102, 582–592. [Google Scholar] [CrossRef]
  138. Nicoli, M.C.; Anese, M.; Manzocco, L.; Lerici, C.R. Antioxidant properties of coffee brews in relation to the roasting degree. Food Sci. Technol.-Lebensm. Wiss. Technol. 1997, 30, 292–297. [Google Scholar] [CrossRef]
  139. Lopez-Galilea, I.; de Peña, M.P.; Cid, C. Application of multivariate analysis to investigate potential antioxidants in conventional and torrefacto roasted coffee. Eur. Food Res. Technol. 2008, 227, 141–149. [Google Scholar] [CrossRef]
  140. Yashin, A.; Yashin, Y.; Wang, J.Y.; Nemzer, B. Antioxidant and Antiradical Activity of Coffee. Antioxidants 2013, 2, 230–245. [Google Scholar] [CrossRef] [PubMed]
  141. Ramadan-Hassanien, M.F. Total antioxidant potential of juices, beverages and hot drinks consumed in Egypt screened by DPPH in vitro assay. Grasas Aceites 2008, 59, 254–259. [Google Scholar] [CrossRef] [Green Version]
  142. Abdelwareth, A.; Zayed, A.; Farag, M.A. Chemometrics-based aroma profiling for revealing origin, roasting indices, and brewing method in coffee seeds and its commercial blends in the Middle East. Food Chem. 2021, 349, 129162. [Google Scholar] [CrossRef] [PubMed]
Figure 1. OPLS-DA score plots for brews of Kaua’i green (KG, n = 3) and Waialua green (WG, n = 3) coffees for known (A) global metabolites, (B) positive ESI mode lipids, (C) negative ESI mode lipids, and (D) biogenic amines.
Figure 1. OPLS-DA score plots for brews of Kaua’i green (KG, n = 3) and Waialua green (WG, n = 3) coffees for known (A) global metabolites, (B) positive ESI mode lipids, (C) negative ESI mode lipids, and (D) biogenic amines.
Metabolites 13 00412 g001
Figure 2. OPLS-DA score plots for brews of roasted Kaua’i coffee (KR, n = 3) and roasted Waialua coffee (WR, n = 3) for known (A) global metabolites, (B) positive ESI mode lipids, (C) negative ESI mode lipids, and (D) biogenic amines. Metabolites 13 00412 i001 = Kaua’i roasted (KR) and Metabolites 13 00412 i002 = Waialua roasted (WR).
Figure 2. OPLS-DA score plots for brews of roasted Kaua’i coffee (KR, n = 3) and roasted Waialua coffee (WR, n = 3) for known (A) global metabolites, (B) positive ESI mode lipids, (C) negative ESI mode lipids, and (D) biogenic amines. Metabolites 13 00412 i001 = Kaua’i roasted (KR) and Metabolites 13 00412 i002 = Waialua roasted (WR).
Metabolites 13 00412 g002
Figure 3. Effect of roasting on caffeine (A), chlorogenic acid (B), and pyrogallol (C), analyzed using GraphPad Prism 7.0. Values are mean ± SD (n = 8). a,b Mean values with common letters do not differ (p < 0.05). Heatmaps of Kaua’i and Waialua coffee for known global metabolites (D), positive ESI mode lipids (E), negative ESI mode lipids (F), and biogenic amines (G).
Figure 3. Effect of roasting on caffeine (A), chlorogenic acid (B), and pyrogallol (C), analyzed using GraphPad Prism 7.0. Values are mean ± SD (n = 8). a,b Mean values with common letters do not differ (p < 0.05). Heatmaps of Kaua’i and Waialua coffee for known global metabolites (D), positive ESI mode lipids (E), negative ESI mode lipids (F), and biogenic amines (G).
Metabolites 13 00412 g003aMetabolites 13 00412 g003b
Figure 4. Relative changes of selected known metabolites (AE), positive ESI lipids (FJ), negative ESI lipids (KO), and amines (PT) in green (n = 3) vs. roasted (n = 8) coffees (p < 0.05). Metabolites 13 00412 i003 KG = Kaua’i green, Metabolites 13 00412 i004 KR = Kaua’i roasted, Metabolites 13 00412 i005 WG = Waialua green, and Metabolites 13 00412 i006 WR = Waialua roasted.
Figure 4. Relative changes of selected known metabolites (AE), positive ESI lipids (FJ), negative ESI lipids (KO), and amines (PT) in green (n = 3) vs. roasted (n = 8) coffees (p < 0.05). Metabolites 13 00412 i003 KG = Kaua’i green, Metabolites 13 00412 i004 KR = Kaua’i roasted, Metabolites 13 00412 i005 WG = Waialua green, and Metabolites 13 00412 i006 WR = Waialua roasted.
Metabolites 13 00412 g004aMetabolites 13 00412 g004b
Figure 5. Relative changes of selected known metabolites (AD) and positive (FH) and negative (IL) ESI mode lipids in different coffee brews are depicted (n = 4, p < 0.05). TG, triacylglycerol; LPC, lysophosphatidylcholine; FA, fatty acid; PE, phosphatidylethanolamine. Metabolites 13 00412 i007 cold brew; Metabolites 13 00412 i008 filter mesh; Metabolites 13 00412 i009 filter paper; Metabolites 13 00412 i010 French press.
Figure 5. Relative changes of selected known metabolites (AD) and positive (FH) and negative (IL) ESI mode lipids in different coffee brews are depicted (n = 4, p < 0.05). TG, triacylglycerol; LPC, lysophosphatidylcholine; FA, fatty acid; PE, phosphatidylethanolamine. Metabolites 13 00412 i007 cold brew; Metabolites 13 00412 i008 filter mesh; Metabolites 13 00412 i009 filter paper; Metabolites 13 00412 i010 French press.
Metabolites 13 00412 g005
Figure 6. Antioxidant capacity of Kaua’i and Waialua coffee. FrP, French press; FP, filter paper; FM, filter mesh; CB, cold brew; K, Kailua; W, Waialua. Data are represented as mean ± SD (n = 18). a,b Mean values with common letters do not differ (p < 0.05).
Figure 6. Antioxidant capacity of Kaua’i and Waialua coffee. FrP, French press; FP, filter paper; FM, filter mesh; CB, cold brew; K, Kailua; W, Waialua. Data are represented as mean ± SD (n = 18). a,b Mean values with common letters do not differ (p < 0.05).
Metabolites 13 00412 g006
Table 1. Mineral contents of Kaua’i and Waialua coffee prepared by different brewing methods. a,b,c Mean values with common letters do not differ (p < 0.05).
Table 1. Mineral contents of Kaua’i and Waialua coffee prepared by different brewing methods. a,b,c Mean values with common letters do not differ (p < 0.05).
MineralAverage RDA, AI, and UL in Males and Females in mg/day (19-70 years)Minerals in 8 oz of Kaua’i Coffee (mg)Minerals in 8 oz of Waialua Coffee (mg)
French Press (K-FrP)Filter Paper (K-FP)Filter Mesh (K-FM)Cold Brew (K-CB)French Press (W-FrP)Filter Paper (W-FP)Filter Mesh (W-FM)Cold Brew (W-CB)
P (mg.day−1) 1700 *
1250 **
5.721 ± 0.26 a3.352 ± 0.42 a4.524 ± 1.46 a12.128 ± 1.30 b5.511 ± 0.004 a4.315 ± 0.13 a4.572 ± 0.45 a14.074 ± 1.36 b
K (mg.day−1) 22600 *
4700 **
110.395 ± 9.70 a,b79.83 ± 7.16 a108.60 ± 33.17 a,b137.26 ± 16.88 b112.26 ± 2.57 a,b101.65 ± 3.42 a,b112.74 ± 8.60 a,b171.61 ± 17.77 c
Ca (mg.day−1) 11000 *
1300 **
2.883 ± 0.19 a2.116 ± 0.23 a2.363 ± 0.80 a5.711 ± 0.56 b2.299 ± 0.10 a2.032 ± 0.18 a1.946 ± 0.12 a6.871 ± 0.36 c
Mg (mg.day−1) 1310–320 *
400–420 *
11.17 ± 1.28 a7.31 ± 0.14 b8.69 ± 2.26 a,b16.78 ± 0.97 c8.92 ± 0.18 a,b7.45 ± 0.42 b7.562± 0.65 b18.61 ± 1.73 c
Na (mg.day−1) 22300 *7.00 ± 0.008 a5.384 ± 0.20 a5.977 ± 1.63 a8.773 ± 0.23 b,c6.654 ± 0.05 a,b5.92 ± 0.13 a6.309 ± 0.86 a9.99 ± 1.26 c
Fe (mg.day−1) 118 *
8 *
0.0084 ± 0.0011 a,d0.0073 ± 0.0007 d0.01 ± 0.002 a,d0.017 ± 0.001 b0.011 ± 0.0008 a0.011 ± 0.00001 a,d0.012 ± 0.001 a0.029 ± 0.003 c
Mn (mg.day−1) 22.3 **0.056 ± 0.004 a0.036 ± 0.008 a0.05 ± 0.02 a0.156 ± 0.03 b0.057 ± 0.0004 a0.046 ± 0.003 a0.05 ± 0.003 a0.252 ± 0.009 c
Zn (mg.day−1) 18 *
11 *
0.006 ± 0.0002 a0.004 ± 0.001 a0.006 ± 0.001 a0.014 ± 0.003 b0.006 ± 0.001 a0.004 ± 0.00002 a0.007 ± 0.003 a,b0.023 ± 0.006 c
Cu (μg.day−1) 10.9 **0.005 ± 0.0001 a0.004 ± 0.0004 a0.006 ± 0.001 a0.009 ± 0.001 a0.005 ± 0.0001 a0.006 ± 0.001 a0.006 ± 0.001 a0.034 ± 0.023 b
B (mg.day−1) 31.50.042 ± 0.003 a,c0.035 ± 0.003 a,b0.048± 0.01 a.c0.025 ± 0.003 b0.044 ± 0.0008 a,c0.048 ± 0.003 a,c0.052± 0.003 c0.039 ± 0.004 a
1 Recommended dietary allowance (RDA); 2 adequate intake (AI); 3 tolerable upper intake level (UL); K-, Kaua’i; W-, Waialua; FrP, French press; FP, filter paper; FM, filter mesh; CB, cold brew. Sourced from * the Dietary Guidelines for Americans 2020–2025, ** FDA.gov (accessed on 10 November 2022).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nerurkar, P.V.; Yokoyama, J.; Ichimura, K.; Kutscher, S.; Wong, J.; Bittenbender, H.C.; Deng, Y. Medium Roasting and Brewing Methods Differentially Modulate Global Metabolites, Lipids, Biogenic Amines, Minerals, and Antioxidant Capacity of Hawai‘i-Grown Coffee (Coffea arabica). Metabolites 2023, 13, 412. https://doi.org/10.3390/metabo13030412

AMA Style

Nerurkar PV, Yokoyama J, Ichimura K, Kutscher S, Wong J, Bittenbender HC, Deng Y. Medium Roasting and Brewing Methods Differentially Modulate Global Metabolites, Lipids, Biogenic Amines, Minerals, and Antioxidant Capacity of Hawai‘i-Grown Coffee (Coffea arabica). Metabolites. 2023; 13(3):412. https://doi.org/10.3390/metabo13030412

Chicago/Turabian Style

Nerurkar, Pratibha V., Jennifer Yokoyama, Kramer Ichimura, Shannon Kutscher, Jamie Wong, Harry C. Bittenbender, and Youping Deng. 2023. "Medium Roasting and Brewing Methods Differentially Modulate Global Metabolites, Lipids, Biogenic Amines, Minerals, and Antioxidant Capacity of Hawai‘i-Grown Coffee (Coffea arabica)" Metabolites 13, no. 3: 412. https://doi.org/10.3390/metabo13030412

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