# Faba Bean (Vicia faba L. minor) Bitterness: An Untargeted Metabolomic Approach to Highlight the Impact of the Non-Volatile Fraction

Adeline Karolkowski, Emmanuelle Meudec, Antoine Bruguière, Anne-Claire

Mitaine-Offer, Emilie Bouzidi, Loïc Levavasseur, Nicolas Sommerer, Loïc

Briand, Christian Salles

► To cite this version:

Adeline Karolkowski, Emmanuelle Meudec, Antoine Bruguière, Anne-Claire Mitaine-Offer, Emilie Bouzidi, et al. Faba Bean (Vicia faba L. minor) Bitterness: An Untargeted Metabolomic Approach to Highlight the Impact of the Non-Volatile Fraction. Metabolites, 2023, 13 (8), pp.964. 10.3390/metabo13080964. hal-04213983

HAL Id: hal-04213983

https://hal.inrae.fr/hal-04213983

Submitted on 21 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Faba Bean (*Vicia faba L. minor*) Bitterness: An Untargeted Metabolomic Approach to Highlight the Impact of the Non-Volatile Fraction

Adeline Karolkowski <fnref n="1" />,<fnref n="2" />, Emmanuelle Meudec <fnref n="3" />,<fnref n="4" /> ID, Antoine Bruguière <fnref n="1" />, Anne-Claire Mitaine-Offer <fnref n="1" />, Emilie Bouzidi <fnref n="5" />, Loïc Levavasseur <fnref n="2" />, Nicolas Sommerer <fnref n="3" />,<fnref n="4" /> ID, Loïc Briand <fnref n="1" />,* and Christian Salles <fnref n="1" />,* ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Karolkowski, A.; Meudec, E.; Bruguière, A.; Mitaine-Offer, A.-C.; Bouzidi, E.; Levavasseur, L.; Sommerer, N.; Briand, L.; Salles, C. Faba Bean (*Vicia faba L. minor*)

**Bitterness:** An Untargeted Metabolomic Approach to Highlight the Impact of the Non-Volatile Fraction. *Metabolites* **2023**, *13*, 964. https://doi.org/10.3390/metabo13080964

Received: 2 August 2023

Revised: 17 August 2023

Accepted: 18 August 2023

Published: 21 August 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** In the context of climate change, faba beans are an interesting alternative to animal proteins but are characterised by off-notes and bitterness that decrease consumer acceptability. However, research on pulse bitterness is often limited to soybeans and peas. This study aimed to highlight potential bitter non-volatile compounds in faba beans. First, the bitterness of flours and air-classified fractions (starch and protein) of three faba bean cultivars was evaluated by a trained panel. The fractions from the high-alkaloid cultivars and the protein fractions exhibited higher bitter intensity. Second, an untargeted metabolomic approach using ultra-high-performance liquid chromatography–diode array detector–tandem–high resolution mass spectrometry (UHPLC–DAD–HRMS) was correlated with the bitter perception of the fractions. Third, 42 tentatively identified non-volatile compounds were associated with faba bean bitterness by correlated sensory and metabolomic data. These compounds mainly belonged to different chemical classes such as alkaloids, amino acids, phenolic compounds, organic acids, and terpenoids. This research provided a better understanding of the molecules responsible for bitterness in faba beans and the impact of cultivar and air-classification on the bitter content. The bitter character of these highlighted compounds needs to be confirmed by sensory and/or cellular analyses to identify removal or masking strategies.

**Keywords:** pulses; faba beans; off-flavours; bitterness; sensory analysis; UHPLC-HRMS; metabolomic approach

## 1. Introduction

In the context of climate change, it seems appropriate to reconcile the economic, social, and environmental impacts of human activities. The transition to a more plant-based diet, in particular in pulses, thus appears to be one of the main levers for improving human health and the sustainability of food systems [1]. For example, faba bean (*Vicia faba L. minor*) ingredients are increasingly used in the formulation of food products due to agronomic, environmental, nutritional, and functional interests [2,3]. Despite the great interest in pulses, consumer acceptability remains low due to off-flavours, such as off-notes, bitterness, and astringency [4]. However, faba beans are still of sensory interest. They are not astringent and have a lower intensity of off-notes than other pulses, such as peas or soybeans [5,6]. Research on the off-flavours of faba beans has focused on the involvement of volatile compounds in off-notes [7,8]. In addition, a few studies have been carried out on the

<footnote n="1">¹ Centre des Sciences du Goût et de L'Alimentation, CNRS, INRAE, Institut Agro, Université de Bourgogne, F-21000 Dijon, France; adeline.karolkowski@inrae.fr (A.K.); antoine.bruguiere@u-bourgogne.fr (A.B.); anne-claireoffert@u-bourgogne.fr (A.-C.M.-O.)</footnote>

<footnote n="2">² Groupe Soufflet-Invivo, F-10400 Nogent-sur-Seine, France; llevavasseur@invivo-group.com</footnote>

<footnote n="3">³ SPO, Université de Montpellier, INRAE, Institut Agro, F-34000 Montpellier, France; emmanuelle.meudec@inrae.fr (E.M.); nicolas.sommerer@inrae.fr (N.S.)</footnote>

<footnote n="4">⁴ INRAE, PROBE Research Infrastructure, PFP Polyphenol Analysis Facility, F-34060 Montpellier, France</footnote>

<footnote n="5">5 Vivien Paille (Groupe Avril), F-59300 Valenciennes, France; ebouzidi@vivienpaille.fr</footnote>

<footnote>* Correspondence: loic.briand@inrae.fr (L.B.); christian.salles@inrae.fr (C.S.)</footnote>
molecules responsible for the bitterness of pulses, but they were mainly conducted on peas and soybeans [9].

Many compounds from plant defence metabolisms could be responsible for pulse bitterness [9]. Recent studies have highlighted the role of saponins, phenolic compounds, peptides, lipids, and lipid oxidation products in the bitterness of peas [10–13]. Concerning lupins, bitter perception should be more related to specific alkaloids and tannins [14]. Such research is similarly scarce for the bitter compounds in faba beans. In general, only the profile of non-volatile compounds has been studied in faba beans without making a link to their sensory perception. The main identified metabolites belonged to alkaloids, phenolic compounds, saponins, jasmonates, organic acids, and amino acids [15–22]. It has been suggested that the presence of free phenolic compounds, vicine and convicine was linked to bitterness in faba bean ingredients (flour, protein concentrate, and isolate). Moreover, saponins and tannins have been proposed to play a minor role in bitter perception [8]. This observation was partly confirmed by a recent study that demonstrated that the concentration of saponins in faba beans was too low to contribute to their bitterness. However, vicine (one of the major alkaloids) activated one of the 25 human bitter taste receptors (TAS2Rs), and its concentration in faba beans should partly account for bitterness [23].

The aim of this study was to determine the non-volatile compounds that could contribute to faba bean bitterness by linking chemical and sensory data. To this end, three cultivars and three fractions (flour (F), starch fraction (S), and protein fraction (P) (concentrate)) of faba beans were studied. First, the sensory bitterness of the fractions was evaluated by a trained panel. Second, the non-volatile profiles of the different samples were characterised by a metabolomic approach using ultra-high-performance liquid chromatography–diode array detector–tandem–high resolution mass spectrometry (UHPLC–DAS–HRMS). Third, the metabolomic and sensory data were correlated, and the compounds positively correlated with bitterness were tentatively identified. All the selected sensory and physico-chemical methods have already been validated and used to make a link between the chemical composition of mushrooms or peas with their sensory attributes [10,24].

## 2. Materials and Methods

### 2.1. Faba Bean Fractions

Three cultivars (C1, C2, and C3) of faba beans (*Vicia faba L. minor*) were studied. Mature seeds were harvested in 2020 and provided by Soufflet Agriculture (Groupe Soufflet-Invivo, Nogent-sur-Seine, France). C1 was cultivated under organic conditions. C1 and C2 were richer in alkaloids (vicine and convicine) than C3. The seeds were dehulled, ground, and air-classified by an external laboratory (SAS IMPROVE, Amiens, France). The F, P, and S fractions obtained for each cultivar were vacuum-packed in a glass container and stored at –20 °C before analysis.

### 2.2. Sensory Analysis

Faba bean fractions were incorporated into a gel made of xanthan gum and water. Twenty-one trained panellists (35–73 years old, 12 females and 9 males) evaluated the bitterness of the 9 gels along a linear scale (from 0 (absence) to 10 (very intense)). A solution of caffeine (0.5 g/L) was used as an external reference for bitterness. Additional details on the gel formulation and sensory profiling can be found elsewhere [23,25].

### 2.3. Non-Volatile Compound Extraction

Extraction of the non-volatile compounds was adapted from previous protocols [10,26]. For each fraction, 100 mg of sample was extracted twice with 5 mL of a methanol (Bio-solve Chimie, Dieuze, France) and formic acid (Biosolve Chimie, Dieuze, France) mixture (99/1 *v*/*v*) under a 1 h constant agitation at room temperature. The suspensions were kept at –20 °C for 20 min and centrifuged (4500 rpm, 20 min, 4 °C). Then, the supernatant was evaporated at 35 °C using a centrifugal evaporator (EZ-2 Plus Evaporating System, Genevac, Ipswich, UK). After evaporation, the samples were solubilised in 200 µL of a methanol,
Milli-Q® ultrapure water (Merck Millipore, USA), and formic acid mixture (49.5/49.5/1 v/v) before being centrifuged (15,000 rpm, 15 min, 4 °C). The supernatants were collected through a 0.22 µM PTFE (polytetrafluoroethylene) filter, placed in HPLC vials and stored at -80 °C in darkness until analysis. The samples were prepared in triplicate.

## 2.4. Non Volatile Compound Untargeted Analysis

The extracts were analysed by UHPLC (Vanquish, Thermo Scientific, Waltham, MA, USA) using an Acquity UPLC HSST3 C18 column (100 mm × 1 mm i.d., 1.7 µm; Waters, Saint-Quentin-en-Yvelines, France). The mobile phase consisted of (A) water–formic acid (99/1 v/v) and (B) acetonitrile–water–formic acid (79.5/19.5/1 v/v/v). The gradient program with a flow rate of 0.22 mL/min was set as follows: 1–1.5 min A/B, 98/2%; 4.5–7 min A/B, 88/12%; 12 min A/B, 76/24%; 15 min A/B, 72/28%; and 16 min A/B, 40/60%. The column and injector temperatures were maintained at 35 and 10 °C, respectively. The injection volume was 0.5 µL. The UHPLC system was coupled with a DAD (UV–visible diode array detector) covering the full range of acquisition (190–600 nm) and an HRMS (Orbitrap Exploris™ 480, Thermo Scientific, Waltham, MA, USA) equipped with a heated electrospray ionisation probe. HRMS was operated in both negative and positive ion modes. The parameters for the ion source were as follows: ion transfer tube temperature—280 °C; voltage (+)—3500 V; voltage (−)—2500 V; sheath gas—40 a.u.; auxiliary gas—10 a.u.; sweep gas—2 a.u.; vaporiser temperature—300 °C; mass range—100–1800 Th, resolution—at *m/z* 240,000 and 480,000. The sample sequence was adapted from a previous method [27]. It included extracts from S, F, and P fractions, a quality control (QC, which was a mix of all samples) sample, and blank. Samples were injected in biological triplicate and in a random order. After equilibration, the sequence started with 3 injections of a blank, followed by 5 injections of QC samples. The blank and QC sample were then injected every 7 real samples. The sequence ended with 1 injection of QC sample and 1 injection of blank. Xcalibur™ 4.4 (Thermo Scientific, Waltham, MA, USA) was used for instrument control, data acquisition, and data analysis.

## 2.5. Metabolomic Analysis

Raw data from HRMS were processed with Compound Discoverer™ (v 3.2.2.421, Thermo Scientific, Waltham, MA, USA), which facilitated peak recognition and dereplication of raw data (retention time alignment, adducts, and isotopic peak clustering). Details about the workflow are available in the Supplementary Materials. Then, these data were correlated (Pearson correlation) with the bitter intensity of the faba bean fractions to identify those that were positively correlated. The tentative identification of the positively correlated compounds was based on ultra-high-resolution mass spectrum, mass accuracy below 1 ppm, characteristic MS/MS fragmentation, UV spectrum, and retention time available in the literature. In particular, the isotope-ratio at ultra-high resolution of ¹⁵N vs. ¹³C and ¹⁸O vs. 2x¹³C and ³⁴S was used to determine the chemical formula of the studied compounds (FreeStyle 1.5, Thermo Scientific, Waltham, MA, USA).

## 2.6. Statistical Analysis

Analyses were performed using XLSTAT (Addinsoft, Paris, France). For the bitter intensity, a one-way analysis of variance (ANOVA) was performed, and significant differences were evaluated by Tukey's honest significant difference (HSD) post hoc test (*p* < 0.05). Linear and logarithmic Pearson correlations (α = 2.5%) were used to determine the relationships between the non-volatile compound areas, in both negative and positive modes, and the bitter intensity of the 9 fractions. Only compounds positively correlated with bitterness were then tentatively identified. In addition, a principal component analysis (PCA, Pearson correlation) was also carried out to visually explore the differences in the non-volatile compound profile of the studied samples.
## 3. Results and Discussion

### 3.1. Bitterness of the Faba Bean Fractions

The bitter taste of the nine gels was evaluated by a trained panel (Figure 1). No significant difference in the bitter intensity was observed for the three flours. However, the P fractions, except for P3, were perceived as more bitter than the F and S fractions. P1 was also perceived as more bitter than P2 and P3. These sensory data were correlated with those from the untargeted metabolomic approach to highlight potential bitter compounds in faba bean fractions.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Bitter intensity (over 10) of the 9 fractions evaluated by a trained panel. Significant differences are indicated by different letters (Tukey's HSD test, $\alpha$ = 5.0%). S: starch fraction; F: flour; P: protein fraction—the number after the fraction corresponds to the cultivar (1, 2, or 3).</figcaption></figure>

### 3.2. Tentative Identification of Non-Volatile Compounds Correlated with Faba Bean Bitterness

A list of 355 and 155 potential compounds was established in negative and positive modes, respectively. In the negative mode, 14 non-volatile compounds were positively correlated with bitterness using linear and/or logarithmic correlations, whereas in the positive mode, 31 compounds were positively correlated with bitterness. Only three compounds were detected and positively correlated with bitterness in both positive and negative modes. Figure 2 shows the different non-volatile compounds resulting from the metabolomic analysis of the fractions and those that were positively correlated with bitterness (variables in blue for linear model, variables in dark green for logarithmic model and variables in pink for both models). Few adducts or fragments of these compounds were also positively correlated with bitterness. Those adducts should have been dereplicated by the workflow used but were not. An expert overlooked of the features after dereplication was carried out, and they were easily identified. As there were few mis-dereplications, the workflow was not modified. Those mis-dereplicated compounds were not further considered for this analysis; however, they were represented by black variables on the PCA.

Then, the positively correlated compounds were tentatively identified by comparing the UV spectrum and the MS/MS spectra from the literature. Moreover, the occurrence of tentatively identified compounds in faba beans and/or other pulses/plants was also specified. For all the proposed identifications, the use of high-mass accuracy (<1 ppm) and ultra-high-resolution MS spectra (R = 480 k at *m/z* 200) of the isotopic profile allowed unique and unambiguous molecular formulas. The correlated compounds are presented in Table 1. They are distributed into the following chemical classes: alkaloids, amino acids, phenolic compounds, organic acids, other compounds, and unidentified compounds.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Biplot representation of the PCA (centred reduced variables, Pearson correlation, $\alpha = 5.0\%$) of the detected compound areas in the negative (A) or positive (B) modes and the perceived bitterness (as a supplementary variable in light green) of the 9 faba bean fractions. The compounds positively correlated with bitterness are related to variables in blue for the linear model, in dark green for the logarithmic model and in pink for both the linear and logarithmic models, whereas the black highlighted compounds correspond to the mis-dereplicated data of the positively correlated compounds. S: starch fraction; F: flour; P: protein fraction—the number after the fraction corresponds to the cultivar (1, 2, or 3).</figcaption></figure>

### 3.2.1. Alkaloids

A total of six compounds belonging to the alkaloid family were tentatively identified. Vicine (1) had one absorption band at 274 nm and displayed a major molecular ion at *m/z* 305.1093 in the positive mode. The fragment at *m/z* 143.0564 in the *MS*² spectra should indicate the presence of divicine ([C₄H₆N₄O₂ + H]⁺) and the loss of a hexoside residue (−162). Vicine and convicine are the main alkaloids widely spread in *Vicia faba* [16,28]. Although convicine had been identified in both positive and negative modes, it was not positively correlated with bitterness (linear and logarithmic models). The compounds (2)–(4) and (6) had one absorption band between 270 and 274 nm and exhibited two characteristic fragments of vicine in their *MS*² spectra: *m/z* 305.1093 ([vicine + H]⁺) and 143.0564 *m/z*; they were tentatively identified as vicine derivatives. Only the derivative of vicine (4), *m/z* 387.1521 (negative) and 389.1667 (positive)) has already been identified in faba beans [19,21]. It also demonstrated the loss of the C₅H₈O moiety in both positive and negative modes that could correspond to a valeric or isovaleric residue according to Kowalczyk et al. (2021) [21]. Moreover, a compound derived from convicine (5) was detected at 7.49 min and exhibited a mass of *m/z* 388.1361 and 390.1507 in negative and positive modes, respectively. It was characterised by two fragments: *m/z* 304.0791 (negative mode) corresponding to convicine and *m/z* 142.0264 (negative mode) and 144.0405 (positive mode) corresponding to isouramil [23,28]. The mass difference between the alkaloid derivative and the ion concerning the respective alkaloid was similar (−84) for the compounds (4) and (5), which also suggested the presence of a valeric or isovaleric residue for convicine derivative (5). Many derivatives of alkaloids have been identified in faba beans, even if they did not all correspond to those tentatively identified in this article [16,19,21].

### 3.2.2. Amino Acids

Many amino acids and their derivatives were detected in positive mode and positively correlated with faba bean bitterness. The compound (7), with a mass of *m/z* 175.1192 ([M+H]⁺), was also found under a potassium adduct *m/z* 213.0747 ([M+2K-H]⁺) and showed a neutral molecular formula of C₆H₁₄N₄O₂. The fragments *m/z* 158.0926, 130.0975, 116.0707, 112.0870, and 60.0557 were characteristic of L-arginine [29]. The L-arginine ion and the associated fragments were found in the fragmentation of the three following compounds (8, 9, and 17), suggesting that they were derived from this amino acid. Only
N-formyl-L-arginine (9) showing the same molecular formula and its fragmentation has already been identified in another plant, the black cohosh [30]. At 1.00 min, the mass of *m/z* 198.0762 associated with the fragments *m/z* 181.0496, 152.0707, 139.0390, and 135.0441 was related to L-DOPA (12) [22]. An isomer of L-DOPA hexoside (13) was observed in its protonated (*m/z* 360.1286) and stacked (*m/z* 719.2506) forms. The fragmentation patterns revealed the loss of a hexoside residue (*m/z* 198.0762; -162) and an ammonia moiety (*m/z* 181.0495; -17); these observations were consistent with previous results [16,18]. Few isomers have already been identified in faba beans [16,18], but it was not possible to determine which isomer was positively correlated with bitterness. Moreover, three other L-DOPA derivatives were tentatively identified. They exhibited the same characteristic fragments related to L-DOPA, but only the fragmentation of the compounds (11) and (19) revealed the neutral loss of a hexoside compared to the compound (20). However, it was not possible to determine a molecular formula by the isotope-ratio method due to the presence of other signals at the same retention time for the compound (19). L-phenylalanine (16) and L-tryptophan (18) were also tentatively identified, and their fragmentations were in agreement with a previous study [29]. The compound (15) exhibited the same fragment *m/z* 146.0601 as L-tryptophan, suggesting that it was derived from this amino acid. Concerning the compounds (10) and (14), their fragments (*m/z* 224.0918, 178.0864, and 165.0547) were similar to those of N-acetyl-L-tyrosine [31]. Moreover, their fragmentation patterns revealed the loss of a hexoside residue (*m/z* 224.0918; -162); which suggested that these compounds were two isomers of N-acetyl-L-tyrosine hexoside. N-acetyl-L-tyrosine was produced in common bean plants during fungal pathogen attacks [32]. Finally, it is important to note that L-arginine, L-DOPA, L-tryptophan, L-phenylalanine, and L-tyrosine have already been identified in faba beans [16,18,33], which suggested that the tentatively identified compounds were indeed derived from these amino acids.

### 3.2.3. Phenolic Compounds

Three phenolic compounds, only detected in the negative mode, were positively correlated with bitterness. p-Hydroxybenzoic acid (21) was tentatively identified by comparing the accurate MS¹ and MS² spectra from the literature (obtained on a standard) [10]. The compound (22) had two absorption bands at 257 and 293 nm characteristic of phenolic acids [10,34]. The molecular ion at *m/z* 315.0721 and the fragment *m/z* 153.0194 should indicate the presence of a protocatechuic acid moiety and the loss of a hexoside residue [34]. Concerning the compound (23), the band at 282 nm and the fragments from MS² spectra suggested the presence of methylfukiic acid [16,35]. This compound was identified by NMR in *Piscidia Erythrina* L., belonging to the *Fabaceae* family [35]. p-Hydroxybenzoic acid has already been identified in peas [10], whereas protocatechuic acid hexoside and methylfukiic acid were detected out in faba beans [16,17,19–21].

The compound (24) had two absorption bands at 270 and 340 nm corresponding to the phenolic core and the conjugated system of a flavonoid, respectively [36]. Moreover, its fragment ion at *m/z* 287.0551 indicated the presence of a kaempferol moiety. Many kaempferol derivatives have already been observed in faba beans (leaves, pods, and mature seeds) and peas (leaves and protein isolates) [10,15,19–21].

### 3.2.4. Organic Acids

A total of three organic acids were tentatively characterised. The fragmentation patterns in the positive mode of the compound (25) were consistent with those of hydroxy aspergillic acid, and this compound has already been identified in brown alga [37]. The compound (26) exhibited the fragments *m/z* 190.0499, 172.0394, and 144.0444, which should indicate that it was derived from kynurenic acid, which has been identified in olive fruit [38–40]. The MS¹/MS² spectra of the compound (27) were characteristic of pantothenic acid hexoside [41,42]. Pantothenic acid has been found in faba beans [16,18], whereas pantothenic acid hexoside has also been detected in chickpeas and tomatoes [41,42].
### 3.2.5. Terpenoids

A total of two terpenoids were detected. The compound (28) was tentatively identified as 8-β-D-glucopyranosyloxy-2,7-dimethyl-2,4-decadiene-1,10-dioic acid by comparing the fragmentation patterns with the faba bean literature [16]. Another terpenoid (29) was detected in both negative and positive modes and exhibited the same fragments as dihydrophaseic acid 4'-O-β-D-glucopyranoside; dihydrophaseic acid and its derivatives were detected in faba beans [16].

### 3.2.6. Other Non-Volatile Compounds

Two compounds were also tentatively identified, but they did not belong to a specific class. First, the compound (30) detected at 0.74 min exhibited fragments at *m/z* 104.1071, 60.0808, and 58.0657, indicating that it was derived from choline [43], which has already been detected in faba beans [44]. Choline was also tentatively identified in samples with the same *MS*² spectrum but was not positively correlated with bitterness. Second, the fragments *m/z* 315.1812, 191.0564, 161.0456, 149.0456, 143.0347, and 131.0349 indicated that the compound (31) should be a geraniol pentoxide hexoside [45]. However, it was not possible to determine the corresponding isomer. Geraniol has already been detected in the volatile content of faba beans [7].

### 3.2.7. Unidentified Compounds

A total of 11 compounds were not identified among the 45 metabolites positively correlated with bitterness; however, a few comments could be made. The fragmentation patterns of the compound (35) revealed the loss of a hexoside moiety (*m/z* 321.1010; –162). It was not possible to tentatively identify the compounds (40) and (42) because their fragmentation data were not available in the literature. The compound (40) could be a diterpene glycoside called 19-hydroxycinnzeylanol 19-glucoside, which has already been identified in the leaf of a dicotyledon plant [46]. Concerning the compound (42), the *MS*¹ spectrum was consistent with those of 4-chloro-oxoindole-acetic acid, which have been detected in faba bean seeds [21]. However, these hypotheses must be verified.
<table><caption>**Table 1.** Tentative identification of the non-volatile compounds tentatively identified and positively correlated with faba bean bitterness (linear and logarithmic models; Pearson correlation, $\alpha$ = 2.5%). The *p*-value is shown in bold when a positive correlation is observed (*p*-value < 0.025).</caption><thead><tr><th rowspan="2">No.</th><th rowspan="2">RT (min)</th><th colspan="2">Linear Model</th><th colspan="2">Logarithmic Model</th><th rowspan="2">UV (nm)</th><th rowspan="2">Mode</th><th rowspan="2">Experimental m/z</th><th rowspan="2">Formula (Neutral)</th><th rowspan="2">Expected m/z</th><th rowspan="2">Error (ppm)</th><th rowspan="2">Main MS/MS Fragment Ions</th><th rowspan="2">Compound</th><th rowspan="2">RI</th><th rowspan="2">Ref.</th></tr><tr><th>R</th><th>p-Value</th><th>R</th><th>p-Value</th></tr></thead><tbody><tr><td colspan="14" style="text-align:center; font-weight:bold;">ALKALOIDS</td></tr><tr><td>1</td><td>0.86</td><td>0.737</td><td><strong>0.024</strong></td><td>0.549</td><td>0.126</td><td>274</td><td>POS</td><td>305.1093</td><td>C<sub>10</sub>H<sub>14</sub>N<sub>4</sub>O<sub>7</sub></td><td>305.1092</td><td>0.10</td><td>143.0564 (100)</td><td>Vicine</td><td>1,2,3,4</td><td>[23,28]</td></tr><tr><td>2</td><td>1.27</td><td>0.796</td><td><strong>0.010</strong></td><td>0.598</td><td>0.089</td><td>274</td><td>POS</td><td>391.1097</td><td>C<sub>13</sub>H<sub>18</sub>N<sub>4</sub>O<sub>10</sub></td><td>391.1107</td><td>1.00</td><td>305.1099 (1); 143.0564 (100)</td><td>Vicine derivative</td><td>1,2</td><td></td></tr><tr><td>3</td><td>6.61</td><td>0.83</td><td><strong>0.006</strong></td><td>0.628</td><td>0.070</td><td>278</td><td>POS</td><td>613.1988</td><td>C<sub>25</sub>H<sub>32</sub>N<sub>4</sub>O<sub>14</sub></td><td>613.1988</td><td>0.00</td><td>305.1093 (1); 147.0442 (60); 143.0564 (100)</td><td>Vicine derivative</td><td>1,2</td><td></td></tr><tr><td rowspan="2">4</td><td rowspan="2">7.61</td><td>0.785</td><td><strong>0.012</strong></td><td>0.609</td><td>0.082</td><td>274</td><td>NEG</td><td>775.3110</td><td rowspan="2">C<sub>15</sub>H<sub>24</sub>N<sub>4</sub>O<sub>8</sub></td><td>387.1521</td><td>0.90</td><td>387.1512 (50); 303.0947 (5); 141.0418 (100)</td><td rowspan="2">Vicine derivative (ester with valeric/isovaleric acid)</td><td rowspan="2">1,2,4</td><td rowspan="2">[19,21]</td></tr><tr><td>0.871</td><td><strong>0.002</strong></td><td>0.701</td><td>0.035</td><td>278</td><td>POS</td><td>389.1670</td><td>389.1667</td><td>-0.30</td><td>305.1092 (5); 143.0564 (100)</td></tr><tr><td rowspan="2">5</td><td rowspan="2">7.68</td><td>0.829</td><td><strong>0.006</strong></td><td>0.645</td><td>0.060</td><td>274</td><td>NEG</td><td>777.2793</td><td rowspan="2">C<sub>15</sub>H<sub>23</sub>N<sub>3</sub>O<sub>9</sub></td><td>388.1361</td><td>0.90</td><td>388.1352 (100); 304.0791 (5); 142.0264 (5)</td><td rowspan="2">Convicine derivative (ester with valeric/isovaleric acid)</td><td rowspan="2">1,2,4</td><td rowspan="2">[19,21]</td></tr><tr><td>0.883</td><td><strong>0.002</strong></td><td>0.648</td><td>0.059</td><td>274</td><td>POS</td><td>390.1509</td><td>390.1507</td><td>-0.20</td><td>229.1070 (20); 144.0405 (100); 127.0390 (20); 85.0648 (30); 57.0699 (20)</td></tr><tr><td>6</td><td>11.56</td><td>0.755</td><td><strong>0.019</strong></td><td>0.579</td><td>0.102</td><td>274</td><td>POS</td><td>815.3052</td><td>C<sub>32</sub>H<sub>46</sub>N<sub>8</sub>O<sub>17</sub></td><td>815.3054</td><td>0.20</td><td>305.1093 (5); 143.0565 (100)</td><td>Vicine derivative</td><td>1,2</td><td></td></tr><tr><td colspan="14" style="text-align:center; font-weight:bold;">AMINO ACIDS</td></tr><tr><td>7</td><td>0.71</td><td>0.736</td><td><strong>0.024</strong></td><td>0.695</td><td>0.038</td><td></td><td>POS</td><td>213.0747 (175.1192)</td><td>C<sub>6</sub>H<sub>14</sub>N<sub>4</sub>O<sub>2</sub></td><td>175.1201</td><td>0.90</td><td>175.1192 (100); 158.0925 (20); 130.0975 (10); 116.0707 (30); 112.0869 (5); 60.0557 (20)</td><td>L-arginine</td><td>1,4</td><td>[29,33]</td></tr><tr><td>8</td><td>0.74</td><td>0.885</td><td><strong>0.002</strong></td><td>0.792</td><td><strong>0.011</strong></td><td></td><td>POS</td><td>292.1979</td><td>C<sub>11</sub>H<sub>25</sub>N<sub>5</sub>O<sub>4</sub></td><td>292.1979</td><td>0.00</td><td>175.1190 (70); 158.0926 (10); 118.0863 (100); 116.0705 (1); 60.0556 (1)</td><td>L-arginine derivative</td><td>1</td><td></td></tr><tr><td rowspan="2">9</td><td rowspan="2">0.85</td><td rowspan="2">0.846</td><td rowspan="2"><strong>0.004</strong></td><td rowspan="2">0.811</td><td rowspan="2"><strong>0.008</strong></td><td></td><td>POS</td><td>203.1138</td><td>C<sub>7</sub>H<sub>15</sub>N<sub>4</sub>O<sub>3</sub></td><td>203.1150</td><td>1.20</td><td>203.1143 (100); 186.0876 (10); 175.1191 (20); 158.0928 (10); 144.0657 (20); 130.0974 (1); 116.0707 (10); 112.0870 (5); 88.0870 (5); 70.0651 (1)</td><td>N-formyl-L-arginine</td><td>1,6</td><td>[30]</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>N-acetyl-L-tyrosine hexoside (unknow isomer)</td><td>N-acetyl-L-tyrosine hexoside (unknow isomer)</td><td>1,4</td><td>[16,18]</td></tr><tr><td>10</td><td>0.93</td><td>0.787</td><td><strong>0.012</strong></td><td>0.765</td><td><strong>0.016</strong></td><td></td><td>POS</td><td>178.0863 (386.1447)</td><td>C<sub>17</sub>H<sub>23</sub>NO<sub>9</sub></td><td>386.1456</td><td>0.90</td><td>224.0917 (100); 178.0863 (40); 85.0284 (15)</td><td>L-DOPA hexoside derivative</td><td>1</td><td></td></tr><tr><td>11</td><td>0.96</td><td>0.753</td><td><strong>0.019</strong></td><td>0.706</td><td>0.033</td><td></td><td>POS</td><td>568.1873 (730.2401)</td><td>C<sub>25</sub>H<sub>43</sub>NO<sub>21</sub></td><td>730.2411</td><td>1.00</td><td>198.0761 (100); 181.0496 (40); 152.0707 (30)</td><td>L-DOPA hexoside derivative</td><td>1</td><td></td></tr><tr><td>12</td><td>1.00</td><td>0.888</td><td><strong>0.001</strong></td><td>0.815</td><td><strong>0.007</strong></td><td></td><td>POS</td><td>198.0762</td><td>C<sub>9</sub>H<sub>11</sub>NO<sub>4</sub></td><td>198.9761</td><td>-0.10</td><td>181.0496 (40); 152.0707 (100); 139.0390 (40); 135.0441 (20)</td><td>L-DOPA</td><td>1,4</td><td>[22]</td></tr><tr><td>13</td><td>1.02</td><td>0.797</td><td><strong>0.010</strong></td><td>0.658</td><td>0.054</td><td></td><td>POS</td><td>719.2506 (360.1286)</td><td>C<sub>15</sub>H<sub>21</sub>NO<sub>9</sub></td><td>360.1289</td><td>0.30</td><td>360.1286 (5); 198.0760 (100); 181.0495 (10); 152.0704 (10); 139.0390 (20); 85.0284 (5)</td><td>L-DOPA hexoside (unknow isomer)</td><td>1,4</td><td>[16,18]</td></tr><tr><td>14</td><td>1.11</td><td>0.770</td><td><strong>0.015</strong></td><td>0.736</td><td><strong>0.024</strong></td><td></td><td>POS</td><td>386.1446</td><td>C<sub>17</sub>H<sub>23</sub>NO<sub>9</sub></td><td>386.1456</td><td>1.00</td><td>224.0918 (100); 178.0864 (40); 165.0547 (1); 85.0284 (10)</td><td>N-acetyl-L-tyrosine hexoside (unknow isomer)</td><td>1</td><td></td></tr><tr><td>15</td><td>1.86</td><td>0.844</td><td><strong>0.004</strong></td><td>0.728</td><td><strong>0.026</strong></td><td></td><td>POS</td><td>323.0874</td><td>C<sub>13</sub>H<sub>14</sub>N<sub>2</sub>O<sub>7</sub></td><td>323.0874</td><td>0.00</td><td>146.0601 (100)</td><td>L-tryptophan derivative</td><td>1</td><td></td></tr><tr><td>16</td><td>2.19</td><td>0.813</td><td><strong>0.008</strong></td><td>0.788</td><td><strong>0.012</strong></td><td></td><td>POS</td><td>166.0863</td><td>C<sub>9</sub>H<sub>12</sub>NO<sub>2</sub></td><td>166.0863</td><td>0.00</td><td>149.0597 (5); 120.0808 (100); 103.0542 (5)</td><td>L-phenylalanine</td><td>1,4</td><td>[29,33]</td></tr><tr><td rowspan="2">17</td><td rowspan="2">3.60</td><td rowspan="2">0.900</td><td rowspan="2"><strong>0.001</strong></td><td rowspan="2">0.717</td><td rowspan="2"><strong>0.030</strong></td><td></td><td>POS</td><td>259.1764</td><td>C<sub>11</sub>H<sub>22</sub>N<sub>4</sub>O<sub>3</sub></td><td>259.1765</td><td>0.10</td><td>259.1766 (100); 242.1499 (20); 200.1279 (10); 175.1191 (10); 158.0925 (20); 112.0871 (1); 116.0705 (1); 70.0651 (5)</td><td>L-arginine derivative</td><td>1</td><td></td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>L-tryptophan</td><td>L-tryptophan derivative</td><td>1,4</td><td>[29,33]</td></tr><tr><td>18</td><td>5.30</td><td>0.838</td><td><strong>0.005</strong></td><td>0.794</td><td><strong>0.011</strong></td><td></td><td>POS</td><td>188.0706 (205.0969)</td><td>C<sub>11</sub>H<sub>12</sub>N<sub>2</sub>O<sub>2</sub></td><td>205.0971</td><td>0.20</td><td>205.0969; 146.0602 (100); 118.0652 (10)</td><td>L-DOPA hexoside derivative</td><td>1</td><td></td></tr><tr><td>19</td><td>9.68</td><td>0.734</td><td><strong>0.024</strong></td><td>0.689</td><td>0.040</td><td></td><td>POS</td><td>222.0648 (443.1226)</td><td>ND</td><td>-</td><td>-</td><td>281.0705 (30); 252.0441 (20); 237.0800 (80); 198.0761 (50); 181.0497 (80); 152.0708 (100); 139.0392 (30); 135.0442 (20); 85.0650 (30)</td><td>L-DOPA hexoside derivative</td><td>1</td><td></td></tr><tr><td>20</td><td>11.28</td><td>0.734</td><td><strong>0.024</strong></td><td>0.652</td><td>0.057</td><td></td><td>POS</td><td>444.1867</td><td>C<sub>20</sub>H<sub>29</sub>NO<sub>10</sub></td><td>444.1864</td><td>-0.30</td><td>152.0706 (60); 139.0390 (10); 135.0440 (10); 85.0648 (30)</td><td>L-DOPA derivative</td><td>1</td><td></td></tr><tr><td colspan="14" style="text-align:center; font-weight:bold;">PHENOLIC COMPOUNDS</td></tr><tr><td>21</td><td>1.57</td><td>0.851</td><td><strong>0.004</strong></td><td>0.781</td><td><strong>0.013</strong></td><td>257;</td><td>NEG</td><td>299.0773</td><td>C<sub>13</sub>H<sub>16</sub>O<sub>8</sub></td><td>299.0772</td><td>-0.10</td><td>137.0244 (100); 93.0345 (30)</td><td>p-Hydroxybenzoic hexoside</td><td>1,4,5</td><td>[10,42,47]</td></tr><tr><td>22</td><td>2.63</td><td>0.808</td><td><strong>0.008</strong></td><td>0.730</td><td><strong>0.025</strong></td><td>293</td><td>NEG</td><td>315.0721</td><td>C<sub>13</sub>H<sub>16</sub>O<sub>9</sub></td><td>315.0722</td><td>0.10</td><td>153.0194 (50); 152.0114 (100); 109.0295 (30); 108.0217 (40)</td><td>Protocatechuic acid Hexoside</td><td>1,2,4,5</td><td>[16,34,42,47]</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th rowspan="2">No.</th><th rowspan="2">RT (min)</th><th colspan="2">Linear Model</th><th colspan="2">Logarithmic Model</th><th rowspan="2">UV (nm)</th><th rowspan="2">Mode</th><th rowspan="2">Experimental m/z</th><th rowspan="2">Formula (Neutral)</th><th rowspan="2">Expected m/z</th><th rowspan="2">Error (ppm)</th><th rowspan="2">Main MS/MS Fragment Ions</th><th rowspan="2">Compound</th><th rowspan="2">RI</th><th rowspan="2">Ref.</th></tr><tr><th>R</th><th>p-Value</th><th>R</th><th>p-Value</th></tr></thead><tbody><tr><td colspan="15" style="text-align:center; font-weight:bold;">PHENOLIC COMPOUNDS</td></tr><tr><td>23</td><td>4.94</td><td>0.860</td><td><strong>0.003</strong></td><td>0.842</td><td><strong>0.004</strong></td><td>282</td><td>NEG</td><td>285.0616</td><td>C<sub>12</sub>H<sub>14</sub>O<sub>8</sub></td><td>285.0616</td><td>0.00</td><td>223.0613 (20); 209.0456 (40); 195.0663 (100); 137.0608 (30)</td><td>3'-O-Methylfukiic acid (3-O-methyl (3',4'-dihydroxybenzyl tartaric acid))</td><td>1,2,3,4,5</td><td>[16,35]</td></tr><tr><td>24</td><td>12.80</td><td>0.824</td><td><strong>0.006</strong></td><td>0.742</td><td><strong>0.022</strong></td><td>270; 340</td><td>POS</td><td>595.1658</td><td>C<sub>27</sub>H<sub>30</sub>O<sub>15</sub></td><td>595.1657</td><td>-0.10</td><td>287.0551</td><td>Kaempferol derivative</td><td>1,2,4,5</td><td>[10,15,48]</td></tr><tr><td colspan="15" style="text-align:center; font-weight:bold;">ORGANIC ACIDS</td></tr><tr><td>25</td><td>1.07</td><td>0.732</td><td>0.025</td><td>0.735</td><td><strong>0.024</strong></td><td></td><td>POS</td><td>241.1546</td><td>C<sub>12</sub>H<sub>20</sub>N<sub>2</sub>O<sub>3</sub></td><td>241.1547</td><td>0.10</td><td>241.1546 (100); 242.1585 (10); 196.0965 (1); 168.0365 (1); 128.1069 (5); 84.0444 (5)</td><td>Hydroxy aspergillic acid</td><td>1,6</td><td></td></tr><tr><td>26</td><td>2.87</td><td>0.907</td><td><strong>0.001</strong></td><td>0.816</td><td><strong>0.007</strong></td><td></td><td>POS</td><td>305.0768</td><td>C<sub>14</sub>H<sub>12</sub>N<sub>2</sub>O<sub>6</sub></td><td>305.0779</td><td>1.10</td><td>215.0814 (20); 190.0499 (50); 172.0394 (100); 144.0444 (10)</td><td>Kynurenic acid derivative</td><td>1</td><td>[38,39]</td></tr><tr><td>27</td><td>4.35</td><td>0.902</td><td><strong>0.001</strong></td><td>0.843</td><td><strong>0.004</strong></td><td></td><td>NEG</td><td>380.1552</td><td>C<sub>15</sub>H<sub>27</sub>O<sub>10</sub>N</td><td>380.1562</td><td>1.00</td><td>362.1441 (10); 308.1351 (20); 218.1036 (10); 146.0822(90)</td><td>Pantothenic acid hexoside</td><td>1,5,6</td><td>[41,42]</td></tr><tr><td colspan="15" style="text-align:center; font-weight:bold;">TERPENOIDS</td></tr><tr><td>28</td><td>6.14</td><td>0.778</td><td><strong>0.013</strong></td><td>0.769</td><td><strong>0.016</strong></td><td></td><td>NEG</td><td>403.1602</td><td>C<sub>18</sub>H<sub>28</sub>O<sub>10</sub></td><td>403.1599</td><td>-0.30</td><td>403.1602 (100); 223.0976 (5); 179.1077 (5); 161.0455 (10); 59.0138 (50)</td><td>8-β-D-glucopyranosyloxy-2,7-dimethyl-2,4-deadiene-1,10-dioic acid</td><td>1,4</td><td>[16]</td></tr><tr><td rowspan="2">29</td><td rowspan="2">6.25</td><td>0.780</td><td><strong>0.013</strong></td><td>0.751</td><td><strong>0.020</strong></td><td></td><td>NEG</td><td>887.3914 (443.1916)</td><td rowspan="2">C<sub>21</sub>H<sub>32</sub>O<sub>10</sub></td><td>443.1923</td><td>0.70</td><td>281.1314 (1); 237.1497 (5); 219.1391 (10); 161.0454 (10); 101.0244 (70)</td><td>Dihydrophaseic acid</td><td rowspan="2">1,4</td><td rowspan="2">[16,19,21]</td></tr><tr><td>0.834</td><td><strong>0.005</strong></td><td>0.781</td><td><strong>0.013</strong></td><td></td><td>POS</td><td>467.1887 (445.2068)</td><td>445.2079</td><td>1.10</td><td>284.0916 (50); 143.0563 (100)</td><td>4'-O-β-D-glucopyranoside</td></tr><tr><td colspan="15" style="text-align:center; font-weight:bold;">OTHERS</td></tr><tr><td>30</td><td>0.74</td><td>0.921</td><td><strong>0.000</strong></td><td>0.877</td><td><strong>0.002</strong></td><td></td><td>POS</td><td>221.1859</td><td>C<sub>10</sub>H<sub>24</sub>N<sub>2</sub>O<sub>3</sub></td><td>221.1860</td><td>0.10</td><td>104.1071 (100); 60.0808 (1); 58.0657 (1)</td><td>Choline derivative</td><td>1</td><td></td></tr><tr><td>31</td><td>15.98</td><td>0.765</td><td><strong>0.016</strong></td><td>0.654</td><td>0.056</td><td></td><td>NEG</td><td>447.2228</td><td>C<sub>21</sub>H<sub>36</sub>O<sub>10</sub></td><td>447.2236</td><td>0.80</td><td>315.1812 (20); 191.0564 (10); 161.0456 (70); 149.0456 (5); 143.0347 (10); 131.0349 (5); 113.0244 (70); 101.0244 (100)</td><td>Geraniol pentoside hexoside (unknown isomer)</td><td>1,6</td><td>[45]</td></tr><tr><td colspan="15" style="text-align:center; font-weight:bold;">UNIDENTIFIED COMPOUNDS</td></tr><tr><td>32</td><td>0.87</td><td>0.788</td><td><strong>0.012</strong></td><td>0.683</td><td>0.042</td><td></td><td>NEG</td><td>545.1620</td><td>ND</td><td>-</td><td>-</td><td>201.0709 (20); 196.0614 (30); 142.0509 (100); 100.0404 (40)</td><td>Unknown</td><td></td><td></td></tr><tr><td>33</td><td>1.56</td><td>0.877</td><td><strong>0.002</strong></td><td>0.636</td><td>0.066</td><td></td><td>POS</td><td>317.1092</td><td>C<sub>11</sub>H<sub>16</sub>N<sub>4</sub>O<sub>7</sub></td><td>317.1103</td><td>1.10</td><td>155.0564 (100)</td><td>Unknown</td><td></td><td></td></tr><tr><td>34</td><td>2.16</td><td>0.902</td><td><strong>0.001</strong></td><td>0.747</td><td><strong>0.021</strong></td><td></td><td>POS</td><td>374.1446</td><td>C<sub>16</sub>H<sub>23</sub>NO<sub>9</sub></td><td>374.1457</td><td>1.10</td><td>212.0918 (90); 195.0652 (90); 153.0547 (20); 152.0707 (100); 85.0284 (40); 69.0335 (10)</td><td>Unknown</td><td></td><td></td></tr><tr><td>35</td><td>6.46</td><td>0.792</td><td><strong>0.011</strong></td><td>0.738</td><td><strong>0.023</strong></td><td></td><td>POS</td><td>242.0803 (483.1538)</td><td>ND</td><td>-</td><td>-</td><td>363.1116 (30); 339.1116 (30); 321.1010 (100); 303.0902 (40)</td><td>Unknown (+ hexoside)</td><td></td><td></td></tr><tr><td>36</td><td>6.46</td><td>0.846</td><td><strong>0.004</strong></td><td>0.787</td><td><strong>0.012</strong></td><td></td><td>POS</td><td>490.2646</td><td>C<sub>23</sub>H<sub>39</sub>NO<sub>10</sub></td><td>490.2647</td><td>0.10</td><td>462.0538 (20); 320.0827 (90); 311.0769 (100); 265.1437 (40); 247.1310 (20)</td><td>Unknown</td><td></td><td></td></tr><tr><td>37</td><td>6.46</td><td>0.854</td><td><strong>0.003</strong></td><td>0.779</td><td><strong>0.013</strong></td><td></td><td>POS</td><td>942.3170</td><td>ND</td><td>-</td><td>-</td><td>499.1261 (10); 378.0767 (80); 320.0706 (100)</td><td>Unknown</td><td></td><td></td></tr><tr><td>38</td><td>6.46</td><td>0.863</td><td><strong>0.003</strong></td><td>0.772</td><td><strong>0.015</strong></td><td></td><td>POS</td><td>927.3515</td><td>ND</td><td>-</td><td>-</td><td>483.1540 (100); 363.1117 (40)</td><td>Unknown</td><td></td><td></td></tr><tr><td>39</td><td>7.86</td><td>0.741</td><td><strong>0.022</strong></td><td>0.708</td><td>0.033</td><td></td><td>NEG</td><td>161.0819</td><td>C<sub>7</sub>H<sub>14</sub>O<sub>4</sub></td><td>161.0819</td><td>0.00</td><td>117.0557 (50); 99.0451 (50); 71.0502 (20)</td><td>Unknown</td><td></td><td></td></tr><tr><td rowspan="2">40</td><td rowspan="2">11.82</td><td>0.748</td><td><strong>0.021</strong></td><td>0.683</td><td>0.043</td><td></td><td>NEG</td><td>561.2550</td><td>C<sub>26</sub>H<sub>42</sub>O<sub>13</sub></td><td>561.2553</td><td>0.30</td><td>519.2444 (80); 387.2013 (100); 207.1386 (40); 191.0561 (40); 161.0454 (40); 113.0244 (30); 101.0244 (50); 99.0087 (70); 89.0244 (50); 71.0138 (40)</td><td>Unknown</td><td></td><td></td></tr><tr><td>12.38</td><td>0.824</td><td><strong>0.006</strong></td><td>0.777</td><td><strong>0.014</strong></td><td></td><td>NEG</td><td>529.2652</td><td>C<sub>26</sub>H<sub>42</sub>O<sub>11</sub></td><td>529.2654</td><td>0.20</td><td>180.0222 (70)</td><td>Unknown</td><td></td><td></td></tr><tr><td>41</td><td>14.21</td><td>0.888</td><td><strong>0.001</strong></td><td>0.820</td><td><strong>0.007</strong></td><td></td><td>NEG</td><td>224.0120</td><td>C<sub>10</sub>H<sub>8</sub>CINO<sub>3</sub></td><td>224.0120</td><td>0.00</td><td></td><td>Unknown</td><td>Unknown</td><td></td><td></td></tr></tbody></table>

RT: retention time; RI: reliability of the attempted identification (1: main MS/MS fragment ions, 2: UV, 3: NMR, 4: identified in faba beans; 5: identified in other pulses; 6: identified in other plants); ND: not determined.
### 3.3. *Non-Volatile Compounds Potentially Responsible for Faba Bean Bitterness*

The difference in the bitter intensities has been partly explained by the presence of alkaloids, including vicine and convicine, in faba bean fractions [23]. However, the aim of this research was to gain a better understanding of the non-volatile compounds responsible for bitterness. Thus, linear and logarithmic Pearson correlations (α = 2.5%) were used to determine the relationships between the non-volatile compound areas, in both negative and positive modes, and the bitter intensity of the nine fractions. Indeed, linear regression should be used when the compound concentration is above the threshold, whereas non-linear models are more suitable when its concentration is below the threshold [10,49].

A total of 6 alkaloids, 4 phenolic compounds, 14 amino acids, 3 organic acids, 2 terpenoids, 2 other compounds, and 11 unidentified compounds positively correlated with bitterness were tentatively identified (Table 1). Then, PCA was performed to explore differences in the non-volatile content of the fractions related to bitterness (Figure 2). In the negative mode (Figure 2A), component 1 (47.44%) was more related to the cultivars, whereas component 2 (29.22%) separated the samples among the fractions. Most of the compounds were positively correlated with bitterness according to the linear model alone (6/14) or to both models (8/14). Concerning the positive mode (Figure 2B), component 1 (51.27%) was related to the type of fractions, whereas the cultivars were much more discriminated for component 2 (29.30%). Only one compound (25) was positively correlated with bitterness according to the logarithmic model, whereas the majority was highlighted by linear (12/31) or both models (18/31). Unlike the PCA from the negative mode, there were two groups of positively correlated compounds with bitterness, each of which appeared to be dependent on the P1 (bottom right) or P2 fractions (top right). N-acetyl-L-tyrosine hexoside (10 and 14) and L-DOPA hexoside derivative (19) were more characteristic of P2, whereas P1 was more related to the presence of vicine (1) and alkaloid derivatives (2, 3, 4, 5, and 6), L-arginine (7) and a derivative (17), L-tryptophan (15), and L-DOPA (13) and its derivatives (11, 12, and 15).

Among the alkaloids, vicine (1) is known to activate TAS2R16 and should be partly responsible for faba bean bitterness [23]. This suggested that vicine derivatives (2, 3, 4, and 6) also activated this receptor, notably if their β-glucopyranoside moiety was free to bind to the active site of TAS2R16 [50]. Concerning amino acids, L-phenylalanine (16) has been shown to activate TAS2R1 and 8, whereas L-tryptophan (18) only activates TAS2R4 and 39 receptors [51,52]. These results were consistent with the literature, which confirms that this multidisciplinary approach was relevant for identifying potential bitter compounds. In addition, the compound (15) was a derivative of L-tryptophan and could therefore activate bitter taste receptors. L-tyrosine, also a hydrophobic amino acid, exhibits a bitter taste [53], which suggested that N-acetyl-tyrosine hexoside (10 and 14) was characterised by a bitter taste. L-arginine (7) did not activate TAS2Rs, but this amino acid enhanced the bitter intensity of 6-n-propylthiouracil (PROP) and caffeine [54]. This suggested that the presence of free L-arginine in faba bean fractions could increase the bitter perception of the bitter compounds. There was no information on the bitterness of p-hydroxybenzoic hexoside (21), protocatechuic acid hexose (22), kaempferol derivative (24), and kynurenic acid derivative (26), but the bitter characteristics of the molecules from which they were derived have been determined. p-Hydroxybenzoic acid and protocatechuic acid were perceived as slightly strong bitter and moderately bitter, respectively, at a concentration of 2 g/L in water [55]. TAS2R14 and 39 receptors have been shown to be sensitive to kaempferol at a concentration in the µM range [56]. The bitter detection threshold of kynurenic acid in distilled water was 78.1 ppm [57]. In addition, it would be interesting to identify the compounds that were negatively correlated with bitterness that would potentially mask bitter perception and behave as potential TAS2R blockers.

These results were consistent with our previous research focused on the role of saponins (DDMP soyasaponin and soyasaponin βb) and alkaloids (vicine and convicine) in faba bean bitterness. Indeed, alkaloids should be partly responsible for bitterness, whereas the soyasaponin concentration was too low to activate TAS2R receptors [23]. In this current
study, soyasaponins were detected but were not positively correlated with bitterness. In addition, vicine was also correlated with the bitter taste, which suggested anew the role of this molecule in faba bean off-flavours. Unlike peas, saponins did not contribute to the bitterness of faba beans [10,12].

Most of the tentatively identified compounds were derived from plant defence metabolism. Alkaloids limit the mollusc repellency of lupin plants [58]. It has been shown that water stress promotes the generation of alkaloids in faba beans [59]. Moreover, vicine exhibits fungicidal and insecticidal properties [60]. Amino acids serve mainly as intermediates in the generation of metabolites, but sometimes they have direct defensive functions against herbivores. For example, L-DOPA, a derivative of L-tyrosine, is a precursor or mimic of neurotransmitters that have been associated with insecticidal properties [61]. There was no information on the defence properties of N-acetyl-L-tyrosine hexoside (10 and 14), but N-acetyl-L-tyrosine was produced in common beans during fungal pathogen attacks [32]. Phenolic compounds are known to protect plant tissues against UV irradiation and attacks from herbivores, fungi, and viruses [62]. Concerning organic acids, kynurenic acid, which comes from tryptophan metabolism, is involved in plant defence mechanisms [38,40]. In pine, the catabolism of abscisic acid was activated during infection by *Fusarium circinate*, leading to the production of dihydrophaseic acid [63]. There were two hypotheses to explain the bitterness of a compound. The first hypothesis supposed that bitter molecules in plants would repel aggressors (insects, herbivores) to allow their development for survival. The second hypothesis suggests that bitter compounds are toxic to animals and prevent them from ingesting harmful doses [64]. In this article, the first hypothesis could explain the bitter taste of these compounds, which ensures the plant's healthy development and survival. However, vicine and convicine are dangerous for people suffering from favism [65], and the second hypothesis thus could be plausible for specific compounds to prevent them from being ingested. It should be noted that these molecules, which were potentially responsible for off-flavours, also ensured good agricultural yields and highlighted that cultivars must be carefully selected. Thus, it would be interesting to also identify compounds that were not positively correlated with bitterness and verify if they are also involved in plant defence mechanisms to better select cultivars.

Finally, combining sensory and untargeted metabolomic approaches allowed us to target potential bitter compounds among a high number of metabolites. However, it was possible that some of these compounds were specific to the P1 and P2 fractions without activating the bitter taste receptors. For most compounds, only the chemical class or main moiety from which they were derived was tentatively used to characterise them. It would be interesting to confirm the identification of these compounds by standards and then determine their bitter characteristics by sensory and/or cellular approaches [9]. However, this approach was probably not complete for identifying all the bitter compounds. Indeed, some metabolites could be present in all the samples but insufficient concentrations to be detected. If several molecules were insufficient to activate the same bitter taste receptor alone, it was possible that all the molecules together, in the same sample, activated a TAS2R at the origin of perceived bitterness.

4. Conclusions

In this study, a combination of sensory and untargeted metabolomic analyses was used to tentatively identify the potential non-volatile compounds responsible for bitterness in air-classified faba bean fractions. Bitterness was more related to high-alkaloid cultivars and P fractions. A total of 42 compounds were positively correlated with bitterness, including 6 alkaloids, 4 phenolic compounds, 14 amino acids, 3 organic acids, 2 terpenoids, and 2 other compounds. However, the tentative identification of 11 compounds remains to be determined. To our knowledge, this work is the first one illustrating the variety of compounds from faba beans that are strong markers of bitterness. It should be interesting to apply our experimental design and this untargeted metabolomic approach to other faba bean cultivars. These results showed that bitter compounds in faba beans were different
from those in peas, whose bitterness was more related to flavonoids and saponins. Then, the majority of these compounds were defence metabolites that were generated in the plants or seeds when they were exposed to abiotic or biotic stresses. This suggested that faba beans produced bitter compounds to prevent them from being ingested by pests to ensure their survival. One strategy for reducing the bitterness of faba beans is to select cultivars with low levels of these compounds while ensuring the presence of other defence metabolites to guarantee sensory, agronomic, and economic benefits. It should also be interesting to confirm the identification of the highlighted compounds to verify their ability to activate TAS2Rs by in vitro cellular-based assays or to be sensory detected by panellists. Finally, this approach allowed us to select potential compounds responsible for faba bean bitterness, but it did not consider the interactions between the compounds in the fractions, which could modulate their bitter intensity. However, this approach was relevant for reducing the number of studied compounds and can be applied to other food matrices or to analyse negative perceptions related to volatile compounds in faba beans such as green, rancid, and metallic notes, for example.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo13080964/s1. Supplementary Materials about the workflow.

**Author Contributions:** Conceptualization, A.K., N.S. and C.S.; methodology, A.K.; formal analysis, A.K.; investigation, A.K. and E.M.; resources, E.M., A.B., A.-C.M.-O. and E.B.; writing—original draft preparation, A.K.; writing—review and editing, A.K., E.M., N.S., L.B. and C.S.; supervision, N.S., L.B. and C.S.; project administration, L.L.; funding acquisition, L.B. and C.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** The PhD of Adeline Karolkowski was funded by Groupe Soufflet-Invivo (France), Groupe AVRIL (France), and ANRT (National Agency for Research and Technology, France) (ANRT-CIFRE 2019/1744).

**Institutional Review Board Statement:** The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of Inserm (French National Institute for Health and Medical Research) (No. 21-773, approved in February 2021).

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** The data are not publicly available due to confidential reasons.

**Acknowledgments:** The authors are indebted to Christophe Martin and David Pertuit (CSGA, INRAE, Dijon, France) for their helpful advice in this study. They are also thankful to Damien Florès (SPO/PFP, INRAE, Montpellier, France) for the UHPLC-DAD-HRMS experiments. INRAE, Regional Council of Bourgogne Franche-Comté and the European Regional Development Fund (FEDER) are thanked for the equipment funding.

**Conflicts of Interest:** Adeline Karolkowski and Loïc Levavasseur are employees of Groupe Soufflet-Invivo. Emilie Bouzidi is employee of Vivien Paille (Groupe Avril). The paper reflects the views of the scientists, and not the companies.

## References

1. Willett, W.; Rockström, J.; Loken, B.; Springmann, M.; Lang, T.; Vermeulen, S.; Garnett, T.; Tilman, D.; DeClerck, F.; Wood, A.; et al. Food in the Anthropocene: The EAT-Lancet Commission on Healthy Diets from Sustainable Food Systems. *Lancet* **2019**, *393*, 447–492. [CrossRef]

2. Sharan, S.; Zanghelini, G.; Zotzel, J.; Bonerz, D.; Aschoff, J.; Saint-Eve, A.; Maillard, M.-N. Fava Bean (*Vicia faba* L.) for Food Applications: From Seed to Ingredient Processing and Its Effect on Functional Properties, Antinutritional Factors, Flavor, and Color. *Compr. Rev. Food Sci. Food Saf.* **2021**, *20*, 401–428. [CrossRef] [PubMed]

3. Mayer Labba, I.-C.; Frøkiær, H.; Sandberg, A.-S. Nutritional and Antinutritional Composition of Fava Bean (*Vicia faba* L., Var. Minor) Cultivars. *Food Res. Int.* **2021**, *140*, 110038. [CrossRef] [PubMed]

4. Roland, W.S.U.; Pouvreau, L.; Curran, J.; van de Velde, F.; de Kok, P.M.T. Flavor Aspects of Pulse Ingredients. *Cereal Chem.* **2017**, *94*, 58–65. [CrossRef]

5. Amarowicz, R.; Troszyńska, A.; Baryłko-Pikielna, N.; Shahidi, F. Polyphenolics Extracts from Legume Seeds: Correlations Between Total Antioxidant Activity, Total Phenolics Content, Tannins Content and Astringency. *J. Food Lipids* **2004**, *11*, 278–286. [CrossRef]
6. Khrisanapant, P.; Kebede, B.; Leong, S.Y.; Oey, I. A Comprehensive Characterisation of Volatile and Fatty Acid Profiles of Legume Seeds. *Foods* **2019**, *8*, 651. [CrossRef] [PubMed]

7. Karolkowski, A.; Gourrat, K.; Bouzidi, E.; Albouy, J.-F.; Levavasseur, L.; Briand, L.; Guichard, E.; Salles, C. Origins of Volatile Compounds and Identification of Odour-Active Compounds in Air-Classified Fractions of Faba Bean (*Vicia faba* L. *Minor*). *Food Res. Int.* **2023**, *163*, 112260. [CrossRef]

8. Tuccillo, F.; Kantanen, K.; Wang, Y.; Martin Ramos Diaz, J.; Pulkkinen, M.; Edelmann, M.; Knaapila, A.; Jouppila, K.; Piironen, V.; Lampi, A.-M.; et al. The Flavor of Faba Bean Ingredients and Extrudates: Chemical and Sensory Properties. *Food Res. Int.* **2022**, *162*, 112036. [CrossRef]

9. Karolkowski, A.; Belloir, C.; Briand, L.; Salles, C. Non-Volatile Compounds Involved in Bitterness and Astringency of Pulses: A Review. *Molecules* **2023**, *28*, 3298. [CrossRef]

10. Cosson, A.; Meudec, E.; Ginies, C.; Danel, A.; Lieben, P.; Descamps, N.; Cheynier, V.; Saint-Eve, A.; Souchon, I. Identification and Quantification of Key Phytochemicals in Peas—Linking Compounds with Sensory Attributes. *Food Chem.* **2022**, *385*, 132615. [CrossRef]

11. Cosson, A.; Oliveira Correia, L.; Descamps, N.; Saint-Eve, A.; Souchon, I. Identification and Characterization of the Main Peptides in Pea Protein Isolates Using Ultra High-Performance Liquid Chromatography Coupled with Mass Spectrometry and Bioinformatics Tools. *Food Chem.* **2022**, *367*, 130747. [CrossRef] [PubMed]

12. Heng, L.; Vincken, J.-P.; van Koningsveld, G.; Legger, A.; Gruppen, H.; van Boekel, T.; Roozen, J.; Voragen, F. Bitterness of Saponins and Their Content in Dry Peas. *J. Sci. Food Agric.* **2006**, *86*, 1225-1231. [CrossRef]

13. Gläser, P.; Mittermeier-Kleßinger, V.K.; Spaccasassi, A.; Hofmann, T.; Dawid, C. Quantification and Bitter Taste Contribution of Lipids and Their Oxidation Products in Pea-Protein Isolates (*Pisum sativum* L.). *J. Agric. Food Chem.* **2021**, *69*, 8768–8776. [CrossRef] [PubMed]

14. Dupont, M.S.; Muzquiz, M.; Estrella, I.; Fenwick, G.R.; Price, K.R. Relationship between the Sensory Properties of Lupin Seed with Alkaloid and Tannin Content. *J. Sci. Food Agric.* **1994**, *65*, 95–100. [CrossRef]

15. Neugart, S.; Rohn, S.; Schreiner, M. Identification of Complex, Naturally Occurring Flavonoid Glycosides in *Vicia faba* and *Pisum sativum* Leaves by HPLC-DAD-ESI-MSn and the Genotypic Effect on Their Flavonoid Profile. *Food Res. Int.* **2015**, *76*, 114–121. [CrossRef]

16. Abu-Reidah, I.M.; Contreras, M.d.M.; Arráez-Román, D.; Fernández-Gutiérrez, A.; Segura Carretero, A. UHPLC-ESI-QTOF-MS Based Metabolic Profiling of *Vicia faba* L. (*Fabaceae*) Seeds as a Key Strategy for Characterization in Foodomics. *Electrophoresis* **2014**, *35*, 1571–1581. [CrossRef]

17. Abu-Reidah, I.M.; Arráez-Román, D.; Warad, I.; Fernández-Gutiérrez, A.; Segura-Carretero, A. UHPLC/MS2-Based Approach for the Comprehensive Metabolite Profiling of Bean (*Vicia faba* L.) by-Products: A Promising Source of Bioactive Constituents. *Food Res. Int.* **2017**, *93*, 87–96. [CrossRef]

18. Mekky, R.H.; Thabet, M.M.; Rodríguez-Pérez, C.; Elnaggar, D.M.Y.; Mahrous, E.A.; Segura-Carretero, A.; Abdel-Sattar, E. Comparative Metabolite Profiling and Antioxidant Potentials of Seeds and Sprouts of Three Egyptian Cultivars of *Vicia faba* L. *Food Res. Int.* **2020**, *136*, 109537. [CrossRef] [PubMed]

19. Valente, I.M.; Cabrita, A.R.J.; Malushi, N.; Oliveira, H.M.; Papa, L.; Rodrigues, J.A.; Fonseca, A.J.M.; Maia, M.R.G. Unravelling the Phytonutrients and Antioxidant Properties of European *Vicia faba* L. Seeds. *Food Res. Int.* **2019**, *116*, 888–896. [CrossRef]

20. Valente, I.M.; Maia, M.R.G.; Malushi, N.; Oliveira, H.M.; Papa, L.; Rodrigues, J.A.; Fonseca, A.J.M.; Cabrita, A.R.J. Profiling of Phenolic Compounds and Antioxidant Properties of European Varieties and Cultivars of *Vicia faba* L. Pods. *Phytochemistry* **2018**, *152*, 223–229. [CrossRef]

21. Kowalczyk, M.; Rolnik, A.; Adach, W.; Kluska, M.; Juszczak, M.; Grabarczyk, Ł.; Wozniak, K.; Olas, B.; Stochmal, A. Multifunctional Compounds in the Extract from Mature Seeds of *Vicia faba* Var. *Minor*: Phytochemical Profiling, Antioxidant Activity and Cellular Safety in Human Selected Blood Cells in in Vitro Trials. *Biomed. Pharmacother.* **2021**, *139*, 111718. [CrossRef] [PubMed]

22. Purves, R.W.; Zhang, H.; Khazaei, H.; Vandenberg, A. Rapid Analysis of Medically Relevant Compounds in Faba Bean Seeds Using FAIMS and Mass Spectrometry. *Int. J. Ion Mobil. Spectrom.* **2017**, *20*, 125–135. [CrossRef]

23. Karolkowski, A.; Belloir, C.; Lucchi, G.; Martin, C.; Bouzidi, E.; Levavasseur, L.; Salles, C.; Briand, L. Activation of Bitter Taste Receptors by Saponins and Alkaloids Identified in Faba Beans (*Vicia faba* L. *Minor*). *Food Chem.* **2023**, *426*, 136548. [CrossRef]

24. Aisala, H.; Manninen, H.; Laaksonen, T.; Linderborg, K.M.; Myoda, T.; Hopia, A.; Sandell, M. Linking Volatile and Non-Volatile Compounds to Sensory Profiles and Consumer Liking of Wild Edible Nordic Mushrooms. *Food Chem.* **2020**, *304*, 125403. [CrossRef] [PubMed]

25. Karolkowski, A.; Martin, C.; Bouzidi, E.; Albouy, J.-F.; Levavasseur, L.; Briand, L.; Salles, C. Heat Treatment, Cultivar and Formulation Modify the Sensory Properties and Consumer Acceptability of Gels Containing Faba Bean (*Vicia faba* L. *Minor*) Protein Concentrates. *Foods* **2022**, *11*, 3018. [CrossRef]

26. Mané, C.; Souquet, J.M.; Ollé, D.; Verriés, C.; Véran, F.; Mazerolles, G.; Cheynier, V.; Fulcrand, H. Optimization of Simultaneous Flavanol, Phenolic Acid, and Anthocyanin Extraction from Grapes Using an Experimental Design: Application to the Characterization of Champagne Grape Varieties. *J. Agric. Food Chem.* **2007**, *55*, 7224–7233. [CrossRef] [PubMed]

27. Arapitsas, P.; Mattivi, F. LC-MS Untargeted Protocol for the Analysis of Wine. In *Metabolic Profiling: Methods and Protocols*; Theodoridis, G.A., Gika, H.G., Wilson, I.D., Eds.; Methods in Molecular Biology; Springer: New York, NY, USA, 2018; pp. 225–235. ISBN 978-1-4939-7643-0.
28. Pulkkinen, M.; Gautam, M.; Lampi, A.-M.; Ollilainen, V.; Stoddard, F.; Sontag-Strohm, T.; Salovaara, H.; Piironen, V. Determination of Vicine and Convicine from Faba Bean with an Optimized High-Performance Liquid Chromatographic Method. *Food Res. Int.* **2015**, *76*, 168–177. [CrossRef]

29. Piraud, M.; Vianey-Saban, C.; Petritis, K.; Elfakir, C.; Steghens, J.-P.; Morla, A.; Bouchu, D. ESI-MS/MS Analysis of Underivatised Amino Acids: A New Tool for the Diagnosis of Inherited Disorders of Amino Acid Metabolism. Fragmentation Study of 79 Molecules of Biological Interest in Positive and Negative Ionisation Mode. *Rapid Commun. Mass Spectrom.* **2003**, *17*, 1297–1311. [CrossRef]

30. Nikolić, D.; Gödecke, T.; Chen, S.-N.; White, J.; Lankin, D.C.; Pauli, G.F.; van Breemen, R.B. Mass Spectrometric Dereplication of Nitrogen-Containing Constituents of Black Cohosh (*Cimicifuga racemosa* L.). *Fitoterapia* **2012**, *83*, 441–460. [CrossRef]

31. Li, H.; Qin, Q.; Shi, X.; He, J.; Xu, G. Modified Metabolites Mapping by Liquid Chromatography-High Resolution Mass Spectrometry Using Full Scan/All Ion Fragmentation/Neutral Loss Acquisition. *J. Chromatogr. A* **2019**, *1583*, 80–87. [CrossRef]

32. Chen, L.; Wu, Q.; He, T.; Lan, J.; Ding, L.; Liu, T.; Wu, Q.; Pan, Y.; Chen, T. Transcriptomic and Metabolomic Changes Triggered by *Fusarium solani* in Common Bean (*Phaseolus vulgaris* L.). *Genes* **2020**, *11*, 177. [CrossRef] [PubMed]

33. Alghamdi, S.S. Chemical Composition of Faba Bean (*Vicia faba* L.) Genotypes under Various Water Regimes. *Pak. J. Nutr.* **2009**, *8*, 447–482.

34. Abu-Reidah, I.M.; Arráez-Román, D.; Segura-Carretero, A.; Fernández-Gutiérrez, A. Extensive Characterisation of Bioactive Phenolic Constituents from Globe Artichoke (*Cynara scolymus* L.) by HPLC-DAD-ESI-QTOF-MS. *Food Chem.* **2013**, *141*, 2269–2277. [CrossRef]

35. Heller, W.; Tamm, C. Fukiinsäure Und 3'-O-Methyl-Fukiinsäure, Zwei Phenolische Hydroxycarbonsäuren Aus *Piscidia Erythrina*. *Helv. Chim. Acta* **1975**, *58*, 974–979. [CrossRef]

36. Mabry, T.; Markham, K.R.; Thomas, M.B. *The Systematic Identification of Flavonoids*; Springer Science & Business Media: New York, NY, USA, 1970; ISBN 978-3-642-88458-0.

37. Carvalho, L.G.; McDougall, G.; Allwood, W.; Sungurtas, J.; Verrall, S.; Stewart, D.; Neilson, K.; O'Byrne, M. Tissues and Industrial Co-Products Formed during Alginate Extraction from Laminaria Hyperborea Provide Different Metabolite Profiles Depending on Harvest Season. *J. Appl. Phycol.* **2023**, *35*, 849–865. [CrossRef]

38. Wang, Y.; Sun, W.; Zheng, J.; Xu, C.; Wang, X.; Li, T.; Tang, Y.; Li, Z. Urinary Metabonomic Study of Patients with Acute Coronary Syndrome Using UPLC-QTOF/MS. *J. Chromatogr. B* **2018**, *1100–1101*, 122–130. [CrossRef]

39. Zhang, Z.-H.; Zhao, Y.-Y.; Cheng, X.-L.; Lin, R.-C.; Dai, Z.; Zhou, C. Metabonomic Study of Biochemical Changes in the Rat Urine Induced by *Pinellia ternata* (Thunb.) Berit. *J. Pharm. Biomed. Anal.* **2013**, *85*, 186–193. [CrossRef]

40. Novelli, S.; Gismondi, A.; Di Marco, G.; Canuti, L.; Nanni, V.; Canini, A. Plant Defense Factors Involved in Olea Europaea Resistance against *Xylella Fastidiosa* Infection. *J. Plant Res.* **2019**, *132*, 439–455. [CrossRef] [PubMed]

41. Mintz-Oron, S.; Mandel, T.; Rogachev, I.; Feldberg, L.; Lotan, O.; Yativ, M.; Wang, Z.; Jetter, R.; Venger, I.; Adato, A.; et al. Gene Expression and Metabolism in Tomato Fruit Surface Tissues. *Plant Physiol.* **2008**, *147*, 823–851. [CrossRef]

42. Mekky, R.H.; Contreras, M.d.M.; El-Gindi, M.R.; Abdel-Monem, A.R.; Abdel-Sattar, E.; Segura-Carretero, A. Profiling of Phenolic and Other Compounds from Egyptian Cultivars of Chickpea (*Cicer arietinum* L.) and Antioxidant Activity: A Comparative Study. *RSC Adv.* **2015**, *5*, 17751–17767. [CrossRef]

43. Min, M.-H.; Park, J.-H.; Hur, J.-H.; Shin, H.-C.; Cho, Y.; Kim, D.-D. Formulation and Bioequivalence Studies of Choline Alfoscerate Tablet Comparing with Soft Gelatin Capsule in Healthy Male Volunteers. *Drug Des. Dev. Ther.* **2019**, *13*, 1049–1058. [CrossRef] [PubMed]

44. Hefni, M.; McEntyre, C.; Lever, M.; Slow, S. A Simple HPLC Method with Fluorescence Detection for Choline Quantification in Foods. *Food Anal. Methods* **2015**, *8*, 2401–2408. [CrossRef]

45. Dong, L.; Miettinen, K.; Goedbloed, M.; Verstappen, F.W.A.; Voster, A.; Jongsma, M.A.; Memelink, J.; van der Krol, S.; Bouwmeester, H.J. Characterization of Two Geraniol Synthases from *Valeriana officinalis* and *Lippia dulcis*: Similar Activity but Difference in Subcellular Localization. *Metab. Eng.* **2013**, *20*, 198–211. [CrossRef]

46. Jayeoye, T.J.; Eze, F.N.; Olatunji, O.J.; Tyopine, A.A. Synthesis of Biocompatible Konjac Glucomannan Stabilized Silver Nanoparticles, with Asystasia Gangetica Phenolic Extract for Colorimetric Detection of Mercury (II) Ion. *Sci. Rep.* **2022**, *12*, 9176. [CrossRef] [PubMed]

47. Saied, D.B.; Farag, M.A. How Does Maturity Stage Affect Seeds Metabolome via UPLC/MS Based Molecular Networking and Chemometrics and in Relation to Antioxidant Effect? A Case Study in 4 Major Cereals and Legumes. *Food Chem.* **2023**, *426*, 136491. [CrossRef]

48. García-Cayuela, T.; Gómez-Maqueo, A.; Guajardo-Flores, D.; Welti-Chanes, J.; Cano, M.P. Characterization and Quantification of Individual Betalain and Phenolic Compounds in Mexican and Spanish Prickly Pear (*Opuntia ficus-Indica* L. Mill) Tissues: A Comparative Study. *J. Food Compos. Anal.* **2019**, *76*, 1–13. [CrossRef]

49. Chambers, E.; Koppel, K. Associations of Volatile Compounds with Sensory Aroma and Flavor: The Complex Nature of Flavor. *Molecules* **2013**, *18*, 4887–4905. [CrossRef]

50. Bufe, B.; Hofmann, T.; Krautwurst, D.; Raguse, J.-D.; Meyerhof, W. The Human TAS2R16 Receptor Mediates Bitter Taste in Response to β-Glucopyranosides. *Nat. Genet.* **2002**, *32*, 397–401. [CrossRef]

51. Bassoli, A.; Borgonovo, G.; Caremoli, F.; Mancuso, G. The Taste of D- and L-Amino Acids: In Vitro Binding Assays with Cloned Human Bitter (TAS2Rs) and Sweet (TAS1R2/TAS1R3) Receptors. *Food Chem.* **2014**, *150*, 27–33. [CrossRef]
52. Kohl, S.; Behrens, M.; Dunkel, A.; Hofmann, T.; Meyerhof, W. Amino Acids and Peptides Activate at Least Five Members of the Human Bitter Taste Receptor Family. *J. Agric. Food Chem.* **2013**, *61*, 53–60. [CrossRef]

53. Ishibashi, N.; Sadamori, K.; Yamamoto, O.; Kanehisa, H.; Rouge, K.; Kikuchi, E.; Okai, H.; Fukui, S. Bitterness of Phenylalanine- and Tyrosine-Containing Peptides. *Agric. Biol. Chem.* **1987**, *51*, 3309–3313. [CrossRef]

54. Melis, M.; Arca, M.; Aragoni, M.C.; Cabras, T.; Caltagirone, C.; Castagnola, M.; Crnjar, R.; Messana, I.; Tepper, B.J.; Barbarossa, I.T. Dose-Dependent Effects of L-Arginine on PROP Bitterness Intensity and Latency and Characteristics of the Chemical Interaction between PROP and L-Arginine. *PLoS ONE* **2015**, *10*, e0131104. [CrossRef] [PubMed]

55. Ferrer-Gallego, R.; Hernández-Hierro, J.M.; Rivas-Gonzalo, J.C.; Escribano-Bailón, M.T. Sensory Evaluation of Bitterness and Astringency Sub-Qualities of Wine Phenolic Compounds: Synergistic Effect and Modulation by Aromas. *Food Res. Int.* **2014**, *62*, 1100–1107. [CrossRef]

56. Roland, W.S.U.; van Buren, L.; Gruppen, H.; Driesse, M.; Kouka, R.J.; Smit, G.; Vincken, J.-P. Bitter Taste Receptor Activation by Flavonoids and Isoflavonoids: Modeled Structural Requirements for Activation of HTAS2R14 and HTAS2R39. *J. Agric. Food Chem.* **2013**, *61*, 10454–10466. [CrossRef]

57. Takahashi, K.; Tadenuma, M.; Kitamoto, K.; Sato, S. L-Prolyl-l-Leucine Anhydride A Bitter Compound Formed in Aged Sake. *Agric. Biol. Chem.* **1974**, *38*, 927–932. [CrossRef]

58. Wink, M. Chemical Defense of Lupins. Mollusc-Repellent Properties of Quinolizidine Alkaloids. *Z. Für Naturforschung C* **1984**, *39*, 553–558. [CrossRef]

59. Ulyanych, O.; Poltoretskyi, S.; Liubych, V.; Yatsenko, A.; Yatsenko, V.; Lazariev, O.; Kravchenko, V. Effect of Surface Drip Irrigation and Cultivars on Physiological State and Productivity of Faba Bean Crop. *J. Agric. Sci.* **2021**, *32*, 139–149. [CrossRef]

60. Rimando, A.M.; Duke, S.O. Natural Products for Pest Management. In *Natural Products for Pest Management*; ACS Symposium Series; American Chemical Society: Washington, DC, USA, 2006; Volume 927, pp. 2–21. ISBN 978-0-8412-3933-3.

61. Huang, T.; Jander, G.; de Vos, M. Non-Protein Amino Acids in Plant Defense against Insect Herbivores: Representative Cases and Opportunities for Further Functional Analysis. *Phytochemistry* **2011**, *72*, 1531–1537. [CrossRef] [PubMed]

62. Cheynier, V. Phenolic Compounds: From Plants to Foods. *Phytochem. Rev.* **2012**, *11*, 153–177. [CrossRef]

63. Amaral, J.; Correia, B.; Escandón, M.; Jesus, C.; Seródio, J.; Valledor, L.; Hancock, R.D.; Dinis, L.-T.; Gomez-Cadenas, A.; Alves, A.; et al. Temporal Physiological Response of Pine to *Fusarium circinatum* Infection Is Dependent on Host Susceptibility Level: The Role of ABA Catabolism. *Tree Physiol.* **2021**, *41*, 801–816. [CrossRef]

64. Meyerhof, W.; Batram, C.; Kuhn, C.; Brockhoff, A.; Chudoba, E.; Bufe, B.; Appendino, G.; Behrens, M. The Molecular Receptive Ranges of Human TAS2R Bitter Taste Receptors. *Chem. Senses* **2010**, *35*, 157–170. [CrossRef] [PubMed]

65. Yannai, S.; Marquardt, R.R. Induction of Favism-like Symptoms in the Rat: Effects of Vicine and Divicine in Normal and Buthionine Sulphoximine-Treated Rats. *J. Sci. Food Agric.* **1985**, *36*, 1161–1168. [CrossRef]

**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.