Next Article in Journal
Zeolite Properties, Methods of Synthesis, and Selected Applications
Next Article in Special Issue
Phytochemical Profiling by UHPLC–Q-TOF/MS and Chemopreventive Effect of Aqueous Extract of Moringa oleifera Leaves and Benzyl Isothiocyanate on Murine Mammary Carcinogenesis
Previous Article in Journal
Effect of Sulfuric Acid Corrosion on Flotation Performance of Calcite by Changing Surface Roughness
Previous Article in Special Issue
Ultrasound-Assisted Extraction of Paeonol from Moutan Cortex: Purification and Component Identification of Extract
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploration of the Diversity of Vicine and Convicine Derivatives in Faba Bean (Vicia faba L.) Cultivars: Insights from LC-MS/MS Spectra

1
Biotechnologies and Environmental Analytics Platform (BEAP), Environmental Research and Innovation Department (ERIN), Luxembourg Institute of Science and Technology (LIST), 5, Rue Bommel, L-4940 Hautcharage, Luxembourg
2
NPZ Innovation GmbH, Hohenlieth-Hof 1, 24363 Holtsee, Germany
3
Université de Reims Champagne-Ardenne, INRAE, RIBP USC 1488, 51100 Reims, France
4
Institut fir Biologësch Landwirtschaft an Agrarkultur Luxemburg a.s.b.l (IBLA), 1 Wantergaass, L-7664 Medernach, Luxembourg
5
Louvain Institute of Biomolecular Science and Technology (LIBST), UCLouvain, Croix du 11 Sud 4-5/L7.07.03, B-1348 Louvain-la-Neuve, Belgium
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(5), 1065; https://doi.org/10.3390/molecules29051065
Submission received: 6 December 2023 / Revised: 16 February 2024 / Accepted: 19 February 2024 / Published: 29 February 2024

Abstract

:
While numerous Fabaceae seeds are a good nutritional source of high-quality protein, the use of some species is hampered by toxic effects caused by exposure to metabolites that accumulate in the seeds. One such species is the faba or broad bean (Vicia faba L.), which accumulates vicine and convicine. These two glycoalkaloids cause favism, the breakdown of red blood cells in persons with a glucose-6-phosphate dehydrogenase deficiency. Because this is the most common enzyme deficiency worldwide, faba bean breeding efforts have focused on developing cultivars with low levels of these alkaloids. Consequently, quantification methods have been developed; however, they quantify vicine and convicine only and not the derivatives of these compounds that potentially generate the same bio-active molecules. Based on the recognition of previously unknown (con)vicine-containing compounds, we screened the fragmentation spectra of LC-MS/MS data from five faba bean cultivars using the characteristic fragments generated by (con)vicine. This resulted in the recognition of more than a hundred derivatives, of which 89 were tentatively identified. (Con)vicine was mainly derivatized through the addition of sugars, hydroxycinnamic acids, and dicarboxylic acids, with a group of compounds composed of two (con)vicine residues linked by dicarboxyl fatty acids. In general, the abundance profiles of the different derivatives in the five cultivars mimicked that of vicine and convicine, but some showed a derivative-specific profile. The description of the (con)vicine diversity will impact the interpretation of future studies on the biosynthesis of (con)vicine, and the content in potentially bio-active alkaloids in faba beans may be higher than that represented by the quantification of vicine and convicine alone.

1. Introduction

The Fabaceae or Leguminosae family, the third largest family of flowering plants, contains numerous species whose protein-rich seeds are of interest as food and feed or as a source of compounds/extracts for medical use [1,2]. The application range of legumes will furthermore likely increase due to the continued exploration of poorly used species [3] and the discovery of new/better activities in species that are already valorized. In human nutrition, either the entire pods or the seeds are consumed, with the latter making up the bulk of the consumption. Seeds from species such as soybean, peas, peanuts, lentils, and chickpeas are consumed worldwide and are an indispensable part of many local cuisines, providing protein, dietary fibers, and micronutrients. Furthermore, the capacity to enter into symbiosis with nitrogen-fixing bacteria makes the legumes low-input crops [4] that, depending on the species, also provide specialty phytochemicals such as isoflavones [5]. Globally, there is a growing interest in the cultivation of grain legumes, stemming not only from the demand for protein-rich food in vegan/vegetarian diets but also from supra-national efforts to promote low-input protein production such as the EU Protein Strategy (https://www.europarl.europa.eu/RegData/etudes/BRIE/2023/751426 (accessed on 23 September 2023)).
While legumes are thus attractive for cultivation and human consumption, some species also contain compounds with an “anti-nutritional” effect [6,7]. The protein-rich seeds are attractive for herbivores; therefore, legume seeds are stocked with a plethora of chemicals and proteins that inhibit digestibility to limit this biotic stress. These include saponins, the glycosides of steroids or triterpenes, which are studied for their antifungal activities [8]. Contrary to the saponins produced in numerous Fabaceae, some species contain specific metabolites that are toxic when consumed. Excessive or prolonged consumption of seeds from the grass pea (Lathyrus sativus L.) can cause neurolathyrism, a disease characterized by paralysis of the lower body and degenerative changes to the spinal cord. Neurolathyrism is caused by the neurotoxin l-β-N-oxalyl-α,β-diaminopropionic acid (β-ODAP) present in Lathyrus seeds [9], although the toxicology has most likely been oversimplified [10]. Another legume with seeds containing potentially toxic compounds is Vicia faba L., known as faba beans, fava beans, or broad beans.
Faba beans contain the fungicidal glycoalkaloids vicine and convicine [11], distributed over all plant parts [12] but with the highest concentrations in seeds, reaching 5 and 2 mg/g dry weight. In planta, vicine and convicine are part of the chemical defense against biotic stress, showing, for instance, an antifungal activity [11]. The hydrolysis of the glucose unit results in the release of the aminopyrimidine aglycones divicine and isouramil. After ingesting faba beans, individuals with a low-activity variant of glucose 6-phosphate dehydrogenase in erythrocytes suffer from favism, a life-threatening decrease in viable red blood cells [13,14]. The use of faba beans as animal feed is likewise limited, as this is known to impact productivity negatively, for instance, in poultry farming [15]. Contrary to this, the abundance of vicine and convicine seems to have a negligible impact on weaner pig performance, which is influenced by other antinutritional factors present in the faba beans [16]. The presence of these compounds thus limits the cultivation and use of this otherwise excellent source of protein, carbohydrates, fiber, and interesting phytochemicals, which, furthermore, has highly favorable agricultural and ecological properties [17,18]. Due to the limitations imposed by these compounds, plant breeders use vicine/convicine content as a selection target [19,20], with a low vicine–convicine faba bean line that has been available for 40 years. As an alternative, enzymatic methods have been developed to eliminate vicine and convicine from food matrices that contain faba beans [13,21].
Numerous analytical approaches have been developed for the analysis of glycoalkaloids [22] and specifically for the quantification of (con)vicine. Initially, these methods were based on spectrophotometry [23] and HPLC with UV detection [24]. As for the analysis of natural products in general [25], most current methods for (con)vicine analysis use mass spectrometry hyphenated to HPLC [26]. The HPLC operated is mostly reversed-phase [13], while the hydrophilic nature of vicine and convicine also allows the use of normal-phase separations, for instance, with HILIC [26,27]. However, the highest throughput was obtained using flow injection MS/MS, resulting in a 60 s analysis time [28].
(Con)vicine derivatives were described for the first time in 2021 [29] but had already been annotated as derivatives without further identification in 2019 [30]. Kowalczyk et al. (2021) describe 11 compounds, of which 7 were identified as (con)vicine derivatives based on the characteristic signal of divicine and isouramil in positive mode. During the metabolite profiling of faba bean seedlings with UPLC-MS/MS, we recognized a number of compounds that present the characteristic fragments of vicine and convicine [28]. We realized that the diversity of (con)vicine-based compounds is higher than that described in the literature and that this may have an impact on studies on their biosynthesis; therefore, the data were systematically screened for MS/MS spectra with m/z 141.018 and 141.040 base peaks. This resulted in more than 130 hits, of which the majority were identified based on observed mass shifts; however, the exact connectivity between the different groups could not be established. From the dataset, it is observed that the relative MS signal intensity over the five varieties of most, but not all, derivatives follows that of the mother compounds. We present here the identification of a collection of derivatives of these highly studied compounds, including compounds composed of two (con)vicine residues connected with different dicarboxylic acids.

2. Results

While the negative ion of isouramil is expected to generate an m/z value of 142.026 ([M − H]), the m/z observed experimentally is 141.018. This has been reported previously for the negative mode LC-MS/MS identification of these compounds [30]. Purves et al. developed a Multiple Reaction Monitoring (MRM)-based method for the quantification of vicine/convicine and used the transitions 303 to 141 for vicine and 304 to 141 for convicine [28]. The calculated mass for C4H3N3O3, the [M − 2H] ion of isouramil, is 141.017, which corresponds to the m/z of 141.018 that we observed. For vicine, the observed base peak fragment ion is 141.042, corresponding to C4H5N4O2 or the [M − H] ion of divicine.
After the identification of vicine and convicine (compounds 1 and 2, Table 1), hereafter called the mother compounds, based on these characteristic fragment ions, and a comparison to published MS/MS spectra, it was observed that other compounds showed the same fragments. To extract potential (con)vicine derivatives, the data were filtered for the presence of the fragment at m/z 141.03 ± 0.015 (covering both vicine and convicine) using the PeakView fragment filter tool (version 2.2.0.11391, AB Sciex). The extracted fragment ion chromatogram generated is shown in Figure 1. A manual inspection of the MS/MS spectra, at which the fragment at 141.03 was observed, resulted in the identification of 163 features containing either vicine or convicine. The convicine-based compound m/z 720.322 at RT 27.08 is the derivative with the highest retention time. Compounds based on (con)vicine span the retention time range from 0.5 to 27.08 min and an m/z range from 303 to 1156.
At a retention time of 19.58 min, the m/z 141.03 fragment is observed in the MS/MS spectra of the compounds at m/z 480.124 (compound 38, Table 1) and 961.256 (dimer of 480.124). The exact m/z of 141.017 and 304.077 allows us to identify convicine, while the remainder of the highest fragment peaks (193.05, 175.04, 160.015, 134.036, and 132.02, Figure 2) agree with the MS/MS spectra of a ferulic acid standard and published data on ferulic acid conjugates [31]. Compound 38 is thus feruloyl convicine (calculated m/z 480.126; Figure 2, upper spectrum), one of the convicine derivatives previously identified [29]. The vicine homologue (compound 39; m/z 479.140; calculated m/z 479.142) was subsequently identified at a retention time of 19.95 min (Figure 2, lower spectrum). These identifications were confirmed in positive mode.
There is no fragment with an m/z corresponding to the ferulic acid attached to the aminopyrimidine moiety. However, in both MS/MS spectra in Figure 2, m/z 337.091, corresponding to [feruloyl hexoside − H2O], is observed. This shows that the ferulic acid is attached to the glucose of (con)vicine. The fact that derivatization occurs via the sugar moiety is also confirmed by the observation of the neutral loss of 124.039 amu, corresponding to [divicine − H2O], from most vicine-based compounds. For instance, for compound 50 (Table 1), identified as vicine-diferuloyl hexoside, an in-source fragment with m/z 693.197 (Δm = 1.2 ppm from the calculated m/z of the compound 50 with the in-source loss of divicine) is observed. This confirms that no derivatization takes place on the aminopyrimidine moiety of vicine. The analogy between vicine and convicine and the fact that the same derivatives are found for both mother compounds, allows us to hypothesize that convicine is also derivatized on the glucose moiety.
The gross composition of more than half of the recognized LC-MS/MS features containing (con)vicine could be identified (that is, the groups attached to the mother compounds could be determined), while limited data were obtained on how the groups are linked. Therefore, the identifications reported are tentative, and confirmation by using either standards or NMR is needed. The elimination of those features corresponding to neutral losses and dimer formation resulted in the identification of 89 vicine-/convicine-based compounds, represented in Table 1. In Supplemental Table S1, 42 compounds are represented that contain either vicine or convicine but for which no further identification of the derivatives was possible.
The compounds in Table 1 can be separated into three different groups. The first group (compounds 1 to 25) comprises derivatives with small molecules, with a maximum of six carbons and at least one carboxyl function. The second group of (con)vicine-based compounds are derivatives with phenylpropanoids or guaiacylglycerol (compounds 26 to 52). Most striking in this group are the derivatives containing multiple ferulic acid and hexose residues (compounds 45 to 52). The MS/MS spectra of the compounds resemble the spectra shown in Figure 2, with little to no MS/MS fragments observed above m/z 479/480 (vicine/convicine ferulic acid). This lack of ions precludes the discernment of connectivity between the different groups. Only in the MS/MS spectra of m/z 1155 (compound 52) is some further structural indication observed. Low-intensity fragment ions are observed at 479.146 (vicine ferulic acid) + 162 and 479.146 + 176, which may indicate the formation of a branched rather than a linear ferulic acid–glucose chain. No derivatives were found in which a combination of ferulic acid and hexoses connects two aminopyrimidines.
The third group contains (con)vicine derivatized with groups with more than six carbons, or with such group of more than six carbons linking two (con)vicine residues. The first identification in this group were compounds wherein the derivatizing group is attached to a single (con)vicine, and the formula reported corresponds to the “(con)vicine + the formula of the derivatizing group − H2O”. To illustrate that some of the compounds are formed by a bifunctional linker and two (con)vicine residues, this formula is maintained, and identification is thus reported as “(con)vicine + the formula of the derivatizing group − H2O (con)vicine”, although the linker obviously loses an additional H2O to form the link with the second (con)vicine.
This group of compounds is illustrated by the identification of vicine + C10H14O3, convicine + C10H14O3, and both alkaloids linked with this derivative. In Figure 3, the MS/MS spectrum of compound 84 is shown. The base peak is formed by the fragment at m/z 141.018, and a minor peak is observed at 304.074, indicating that this is a convicine-based compound. The mass difference between the experimental mass at m/z 486.172 and the calculated mass for convicine is 182.094 amu; this corresponds to the formula C10H14O3 (mass error of 1.6 ppm). When the water lost during the formation of the bond with vicine is added, this becomes C10H16O4, a group observed in the MS/MS spectrum at m/z 199.098 (mass error of 3.9 ppm). From the fragment at 199, there is a loss of 43.994 amu, corresponding to CO2. A similar MS/MS is reported for decenedioic acid [32], with a matching chemical formula. The fragment peaks shared between compound 84 and compound 85, and not present in the MS/MS of vicine or convicine, were used in a library search with MZCloud™. The highest-scoring compound from this library search is 2-decenedioic acid, with a formula that matches the observed added mass. A further comparison of the MS/MS spectra of these compounds in both ion modes with available predicted spectra (for instance HMDB0000603) confirms that fragmentation spectra from m/z 486.172 are indeed as expected for convicine linked to decenedioic acid on the glucose moiety.
The derivatizing molecule resulting in the compound shown in Figure 3 is confirmed to be a dicarboxylic acid through the identification of compounds wherein two (con)vicine residues are linked. At RT 21.94 and 22.17, m/z 772.265 and 771.279 are fragmented (compounds 76 and 80, respectively, Figure 4A,B). In Figure 4A, the peak at m/z 686.260, loss of 86.03 amu, indicates that convicine is present, while the MS/MS base peak at 141.04 indicates vicine, as does m/z 303.09, the negative ion of vicine. These two latter fragments are also observed in Figure 4B, while no peak due to the loss of 86.03 is observed.
Details of these spectra confirm that compound 76 indeed contains both vicine and convicine. A detail of spectra Figure 4A shows that both 141.042 and 141.019, characteristic for vicine and convicine, respectively, are observed (Figure 4C, upper spectrum), while in the lower spectrum, a detail of Figure 4B, only the fragment indicating the presence of vicine is observed. Similar to this is the loss of both vicine and convicine, resulting in fragments at m/z 486.177 and 485.195 (Figure 4D upper spectrum), observed in the detail of Figure 4A. In the detail of spectrum 4B, only the loss of vicine to m/z 485.193 is observed (Figure 4D lower spectrum), confirming that compound 80 contains only vicine. Figure 4E shows the lower mass range of Figure 4A, illustrating that the same fragments that were used for the tentative identification of decenedioic acid as the derivatizing molecule of the compound 84 (Figure 3) are present. Taking all this into account, compound 86 is tentatively identified as a vicine linked to a convicine by decenedioic acid. Compound 80 is then two molecules of vicine linked by decenedioic acid.
In a similar way, the other compounds of the third group of (con)vicine derivatives in Table 1, dicarboxylic fatty acids bound to a vicine or convicine or forming a linker between different combinations of these, were identified. To make it clear that these different forms exist, the molecular formula of the dicarboxylic acid with the loss of a single water molecule was used as the annotation.
The 89 compounds represented in Table 1 were used for quantitative analysis. The principal component analysis with these compounds resulted in the separation of varieties D and E from the others on the first principal component, accounting for 58% of the variance (Figure 5). Varieties D and E are separated on the second principal component, accounting for 18% of the variance. Varieties A to C show a higher variance among the replicates and are less well separated on the PCA plot (Figure 5). Most of the compounds follow the profile of the mother compounds, that is, high in varieties A, B, and C and low varieties in D and E. For further quantitative assessment, more analyses are needed, as well as the use of standards for absolute quantification. However, the high MS signal intensities recorded for some of the derivatives reported here indicate that they are present in higher quantities than trace amounts.

3. Discussion

In total, 163 LC-MS features containing one or more convicine or vicine groups were recognized. After eliminating dimers and features due to neutral losses, 131 LC-MS features remained. Of these, the gross structure could be established for 89 compounds, indicated in Table 1 with quantitative data, with the main fragment ions indicated in Supplemental Table S2. Based on the characteristic fragments, the presence of vicine and/or convicine could be established for the 42 remaining compounds (Supplemental Table S1). Although details of the connectivity between the different groups in these compounds could not be established, for all compounds for which this could be observed, derivatizing groups were bound on the sugar moiety of (con)vicine. The relative quantification of the compounds in the different varieties was performed based on the MS signal intensity. Although this is only a proxy for the abundance of a compound, it allows a relative cross-variety comparison of the abundance of a compound. The MS signal intensity data presented in Supplemental Table S2 can therefore only be interpreted as such.
The first group of derivatives in Table 1 (compounds 1 to 25) are derivatives found on numerous other groups of biomolecules. They include malonic, (iso)valeric, and benzoic acid, which are relatively common on biomolecules, including on glycosylated flavonoids isolated from Fabaceae [33]; some of these (con)vicine derivatives were previously tentatively identified in faba bean seeds [29,30]. The identification of these compounds is supported by well-described observations such as the loss of CO2 from malonyl derivatives in negative mode [34] or known masses for derivatives (among others, those described in maize [35]).
Compounds 26 to 52 are derivatives with hydroxycinnamic acids. The highest MS signal intensity was observed for derivatives with ferulic acid (Figure 2). As for the first group of derivatives, derivatives of biomolecules with hydroxycinnamic acids are commonly found in plants. Relatively complex combinations of hydroxycinnamic acids and sugars have been identified on faba flavonoids [36]. Flavonoids carrying derivatives as complex as compounds 45 to 52 were identified in kale [37], and Arabidopsis anthocyanins are known to contain large chains composed of hydroxycinnamic acids and sugars. In the latter case, the addition of hydroxycinnamic acid groups influences the absorption spectrum of the anthocyanins, modulating their protective activity [38].
The third group of compounds (compounds 53 to 89) are tentatively identified as (con)vicine with dicarboxylic acids or two molecules of (con)vicine linked by dicarboxylic acids. While derivatives with short-chain dicarboxylic acids such as malonic, succinic, and (hydroxy)adipic acids are known for a range of biomolecules and are found among the (con)vicine derivatives in Table 1 (compounds 5 to 9), derivatives of biomolecules with longer dicarboxylic acids are rare. One group of such compounds is dicarboxylic acids conjugated to amino acids and amines, such as carnitine [39]. In humans, these acylcarnitines are involved in the transport of fatty acids to the mitochondria, and a transport, stabilizing, and/or protective function can be hypothesized for these (con)vicine derivatives.
The PCA analysis, shown in Figure 5, illustrates that quantifying these 89 compounds in shoots of 14-day-old seedlings allows us to differentiate varieties classified as low-VC from high-VC varieties (Table 2). As is known and illustrated in our dataset, this can also be performed using standardized methods for the absolute quantification of vicine and convicine [26,28]. A correlation analysis of the 89 derivatives in Table 1 resulted in 5 clusters. In cluster 1, 20 derivatives are found that are high in A and B and low in C–E. Vicine (compound 1) is in this cluster, and the 19 other compounds in this cluster all contain vicine, with only compound 55 being convicine linked to vicine. The abundance of compounds in this cluster is thus driven by the abundance of vicine. Contrary to this, cluster 2 contains compounds high in varieties A, C, and to a lesser degree, E. The compounds in this cluster are convicine linked to convicine (compounds 56, 70, and 73) or vicine (compound 79). Surprisingly, these are higher in the three varieties with the lowest convicine/vicine ratio. Cluster 4, with 6 compounds that have a minimum MS signal intensity in shoots from variety B, includes both hydroxybenzoic acid and both dihydroxy-tetradecadienedioic acid derivatives (compounds 17, 19, 62, and 66). The remaining derivatives follow the description of the varieties (Table 2) and the abundance profile of the mother compounds (Supplemental Table S2), with the majority found in cluster 5. The general observation is thus that the MS signal intensity of derivatives follows that of the mother compounds with only eleven derivatives (clusters 2 and 4) that have a clearly different profile.
Based on the tentative identifications presented in Table 1, some general observations can be made. Convicine derivatives have a lower retention time than the corresponding vicine derivative. For most derivatives, the MS signal for convicine derivatives is higher than that of the vicine derivative. The absolute quantification of the mother compounds in different tissues shows that vicine is higher than convicine in the seeds and stems, a ratio that is reversed in the roots [40,41]. The higher MS signal intensity of convicine observed here may be explained by the better ionization of this convicine derivatives.
To validate the presence of these compounds, different samples were created, and their metabolomic datasets were searched with our in-house database containing the 89 compounds in Table 1. In the roots of the five varieties used here, the same (con)vicine derivative were recognized (Supplemental Figure S1A,B), albeit with different ratios. Finally, the identifications were checked against a dataset of faba bean seedlings grown in-house. Extracts were made from the seedlings of 18 varieties of faba beans involved in an annual variety trial of grain-Leguminosae, germinated on pleated paper. The extracts were pooled to obtain an overall view of the different (con)vicine derivatives present, and an analysis was carried out as described in the Section 4. The derivatives presented here were also identified in this dataset (Supplemental Figure S1C,D). Contrary to the samples from NPZ on which this study is based, the samples used for verification were grown and prepared at LIST, thus avoiding artefacts that may be caused during transport. The identification of these derivatives in completely independent samples shows that these derivatives are not artefacts and can be expected to be present in all extracts from faba bean seedlings. Finally, datasets from extracts of seedlings of other Fabaceae (Pisum sativum, Glycine max, Lupinus spp., Cicer arietinum, and Medicago sativa) were explored and, as expected, based on the absence of the mother compounds in these species, no identification of (con)vicine or their derivatives was obtained.
The origin of these derivatives is not known. The key enzyme in the biosynthesis of (con)vicine is only expressed in the seeds of faba plants [12]. However, in the Supplemental Data of this publication, the characteristic fragment ions for vicine (RT 1.065 m/z 305 and 143) are found in all tissues, although at a low MS signal intensity [12]. Other studies quantified vicine and convicine in germinating seeds, stems, and roots at different developmental stages. Vicine and convicine were found in most tested tissues including vegetative excluding vegetative roots and stems [41]. Given the dry weight of vegetative-stage plants compared to seeds, it is unlikely that the redistribution of (con)vicine from seeds is the only source of these compounds [41]. Older studies calculated the total amount of (con)vicine present in seeds and in vegetative tissues and discovered that not all (con)vicine found in vegetative tissues can originate from the amounts present in the seed, indicating de novo synthesis in V. narbonensis [40,42]. For V. faba, a first study indicated that the sum of vicine and convicine present in the seeds accounted for the amounts found in the roots, but would require the conversion of vicine to convicine [40]. This conclusion was later revoked by the same group, showing that both Vicia species must have the capacity to synthesize vicine and convicine in young plants [42]. With the addition of the derivates described here, the amount of divicine and isouramil in seedlings is much higher than represented by only vicine and convicine. The absolute quantification of the -(con)vicine derivatives described here in all tissues including seeds potentially resolves these contradictions.
So far, few (con)vicine derivatives have been identified in faba bean seeds. Valente et al. reported five derivatives and an unidentified compound that is likely also a vicine derivative [30]. Compound 33, reported by Valente et al., m/z 387 and its [2M − H] ion at m/z 775, is identified here as (iso)valeryl vicine (compound 24, Table 1), also identified in [29]. Compound 59 in Valente et al., m/z 527, is identified here as vicine + C12H16O4 (compound 77, Table 1). For two of the (con)vicine derivatives reported by Valente (compound 75 (m/z 553) and compound 37 (m/z 525)), no matching feature could be found in our data, as is also the case for the vicine pentoside reported by Kowalczyk [29]. While this may be due to the presence of seed-specific derivatives or the specific variety used in the different studies, it indicates that the overall diversity in (con)vicine derivatives is higher than described here. The analysis of more Vicia tissues from more varieties will likely result in the discovery of other derivatives.
Whether the (con)vicine derivatives reported in this study have an impact on the quality of food containing faba beans cannot currently be excluded and requires further study. First, although some (con)vicine derivatives have been identified in seeds, the total diversity of derivatives in seeds is not known. Furthermore, the anti-nutritional effects of vicine and convicine are caused by the free aminopyrimidines. Further research is needed to establish whether these derivatives contribute to the pool of bio-active divicine and isouramil in food or after digestion of faba beans containing food or feed.
In conclusion, in this paper, we described the diversity in compounds containing vicine and convicine, which is found reproducibly in the shoots and roots of Vicia faba seedlings. Three groups of compounds can be discerned, of which the third one, derivatives of (con)vicine with bifunctional molecules forming a linker between two (con)vicine residues, is a type of compound that has not been described previously. In general, for the five varieties analyzed, the MS signal intensity of the derivatives follows that of the mother compounds. To establish the impact of these derivatives on current knowledge of faba beans requires more research; the identifications need to be confirmed and refined. Secondly, detailed studies are needed to determine whether, and if so how, these (con)vicine compounds impact the amount of divicine and isouramil present in food containing faba beans. Because vicine and convicine relate to biotic stress [11,43,44], detailed in planta studies are needed to understand the in planta functions of the compounds reported here.

4. Materials and Methods

4.1. Plants and Sampling

In this trial, five different cultivars from the breeder Norddeutsche Pflanzenzucht KG (NPZ) (Holtsee, Germany) were analyzed. Four cultivars or varieties are officially registered in Germany: Taifun (2011), Trumpet (2017), Augusta (2018), and Allison (2019). The fifth variety, “E”, is a non-public breeding line. All varieties are spring types, except for Augusta, which is a winter type. The Taifun and “E” varieties are white-flowering (Zt) and thus tannin-free, and the other three varieties, Augusta, Trumpet, and Allison, are color-flowering (Bb) and tannin-rich. Based on the breeder-internal vicine threshold, the Taifun, Augusta, and Trumpet varieties are classified as high in vicine (HVC) and the Allison and “E” varieties as low in vicine (LVC).
From each cultivar, 30 seeds were randomly selected from a field-produced seed lot and seedlings were grown in the NPZ seed quality laboratory. The standard method for ISTA germination was utilized. Thirty seeds were sown in sand substrate in a plastic box and put in a climate chamber at 20 °C with 16 h under light and 8 h under dark conditions. The boxes were covered with black foil for seven days, before the foil was removed, and the seedlings were grown for another seven days, corresponding to the end of the sprout stage [41]. Three replicates of ten seedlings were sampled; sand was washed off the root system, and the plants were transferred to a closed plastic bag on ice packs. Upon sampling, the fresh seedlings in plastic bags were packaged and shipped per express delivery to LIST in Hautcharage, Luxembourg.

4.2. Metabolite Extraction

Upon arriving at LIST, the samples were lyophilized for 48 h using a Lyocube (Christ, Osterode am Harz, Germany). The seed parts still attached to the shoot were removed, and the lyophilized vegetative tissue was separated into shoots and roots. The respective organs of three individual plants were pooled to form one biological replicate. For each variety, three replicates were formed, the plants were ground in liquid nitrogen with a mortar and pestle, and the powder obtained was stored at −20 °C prior to extraction.
Approximately 100 mg of powder was weighted and 80% MeOH/20% MQ (Milli-Q water, Merck Millipore, Darmstadt, Germany) was added in a 1/12 ratio w/v; this was incubated at 22 °C with 1450 rpm agitation. After 4 h, the samples were centrifuged at 20,000× g for 30 min at 4 °C. The supernatant was recovered and dried under vacuum, after which the dried extract was stored at −20°C prior to analysis.

4.3. UPLC-MS/MS and Identification

The dried extracts were resolubilized in 5% MeOH/MQ in the same volume used for extraction and filtered through a 0.22μm PTFE syringe filter (Millex-LG, Merck KGaA, Darmstadt, Germany). The filtered extracts were analyzed using an Acquity UPLC I-class UPLC system equipped with a diode array detector (Waters, Milford, MA, USA) hyphenated to a TripleTOF 6600 (SCIEX, Framingham, MA, USA) mass spectrometer as previously described [45]. Five μL of the extracts was injected and separated on a reversed-phase column (Acquity UPLC BEH C18 column 2.1 × 100 mm, 1.7 μm particle size; Waters) with the column temperature set at 50 °C and a flow rate of 0.5 mL/min. The eluents were 0.1% (v/v) formic acid in water (A) and 0.1% (v/v) formic acid in acetonitrile (B) with the following gradient: 0 min, 1% B; 4 min, 1% B; 16 min, 5% B; 35 min, 40% B; 45 min, 100% B; 50 min, 100% B; 54 min, 1% B; 60 min, 1% B.
Electrospray ionization (ESI) was performed using the following parameters: source temperature 650 °C; ion spray voltage of 4.5 and −4.5 kV; curtain gas (nitrogen) of 30 psi; nebulizer gas (air) of 55 psi; and turbine gas (air) of 50 psi. The declustering potential was set up at 60 V and −60 V in positive and negative modes, respectively. The precursor charge state selection was set at 1. For information-dependent acquisition, MS scans (from m/z 100 to 2000) were acquired for 175 ms and the 10 most abundant m/z values were selected for MS/MS scans (from m/z 50 to 2000) of 200 ms. A sweeping collision energy of 15 V below and above 15 and −15 V, for the positive and negative modes, respectively, was applied to all precursor ions. After three occurrences, an m/z value was excluded for MS/MS for 2 s. Data for all replicates were acquired in negative mode, while positive mode analysis was run on one replicate of each couple variety × organ, and positive mode data were used to aid in the identification of the compounds.
For the identification and relative quantification, the raw data files were uploaded in Progenesis QI (v2.3, Nonlinear Dynamics, Newcastle upon Tyne, UK) in order to align all runs, normalize the data, and carry out a quantitative analysis based on the groups of samples (varieties). A separate work file was created for the roots and shoots. LC-MS features (retention time × m/z couples) without fragmentation spectrum were filtered out, as were peaks corresponding to known or suspected contaminants; the remaining features were used for identification. Initially, msp files downloaded from NIST, Massbank, and from our in-house database were used, and hits were manually validated. Compounds not found in these databases were identified through a literature search and MS/MS comparison with the following databases: LipidMaps (https://www.lipidmaps.org/); HMBD (https://hmdb.ca/spectra/ms/search) GNPS (https://gnps.ucsd.edu/ProteoSAFe/libraries.jsp); MZCloud™ (https://beta.mzcloud.org/); PubChem (https://pubchem.ncbi.nlm.nih.gov); and Metlin (https://metlin.scripps.edu).
For the semi-targeted identification reported here, a hypothetical identification was proposed after manual spectral interpretation; the structure was generated in Chemdraw (v: 19.0.0.22) and the theoretical MS/MS spectrum of this structure, generated by MetFrag (https://ipb-halle.github.io/MetFrag/) and CFM-ID (https://cfmid.wishartlab.com/) [46,47], were compared to the experimental spectrum. Accepted identifications were subsequently added to our in-house database for dereplication. All identifications reported here are at level 2 and are based on the m/z of the precursor and the comparison of the experimental MS/MS with experimental or predicted MS/MS spectra from the sources cited above [48].

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules29051065/s1. Figure S1: Extracted fragment ion chromatograms of the fragment at 141.03 ± 0.015, covering the base peak fragment of vicine and convicine from the root samples and from the pool of extracts from the seedlings of 18 faba bean varieties. Table S1: Overview of 42 compounds showing the characteristic MS/MS signature of (con)vicine derivatives for which no identification was obtained. Table S2: Quantitative data on the 89 (con)vicine derivatives represented in Table 1, normalized quantitative data of the three replicates of the five varieties, and statistical data on the abundance of these compounds over the varieties; the m/z values of these fragments higher than 10% of the base peak were ordered from the most-intense to the least-intense are likewise provided.

Author Contributions

Conceptualization, K.S., S.Z., S.G. and J.R.; methodology, C.J., K.S., H.T., S.Z. and S.G.; software, K.S.; validation, S.H., H.H., S.Z. and M.M.; formal analysis, K.S., S.H., M.M., H.H. and S.G.; data curation, K.S.; writing—original draft preparation, K.S.; writing—review and editing, all co-authors; supervision, K.S. and C.J.; project administration, S.Z., K.S., J.R. and C.J.; funding acquisition, K.S. and C.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research is part of the TASSILI project, supported by the Luxembourg National Research Fund (FNR), grant number O22/17207949/TASSILI.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data resulting in the here presented study can be found within the article or the Supplementary Materials. The identification database containing MS and MS/MS data as well as meta data concerning the identifications can be obtained on request to the corresponding author.

Acknowledgments

The authors acknowledge Laurent Solinhac and Diego Rios Salgado for technical support.

Conflicts of Interest

Author Simon Goertz and Hanna Tietgen were employed by the NPZ Innovation GmbH. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Usman, M.; Khan, W.R.; Yousaf, N.; Akram, S.; Murtaza, G.; Kudus, K.A.; Ditta, A.; Rosli, Z.; Rajpar, M.N.; Nazre, M. Exploring the Phytochemicals and Anti-Cancer Potential of the Members of Fabaceae Family: A Comprehensive Review. Molecules 2022, 27, 3863. [Google Scholar] [CrossRef]
  2. Moreno-Valdespino, C.A.; Luna-Vital, D.; Camacho-Ruiz, R.M.; Mojica, L. Bioactive proteins and phytochemicals from legumes: Mechanisms of action preventing obesity and type-2 diabetes. Food Res. Int. 2020, 130, 108905. [Google Scholar] [CrossRef]
  3. Grygier, A.; Chakradhari, S.; Ratusz, K.; Rudzińska, M.; Patel, K.S.; Lazdiņa, D.; Górnaś, P. Seven underutilized species of the Fabaceae family with high potential for industrial application as alternative sources of oil and lipophilic bioactive compounds. Ind. Crops Prod. 2022, 186, 115251. [Google Scholar] [CrossRef]
  4. Ritchie, H.; Rosado, P.; Roser, M. Environmental impacts of food production. Our World Data 2021. Available online: https://ourworldindata.org/environmental-impacts-of-food (accessed on 3 October 2023).
  5. Wang, Q.; Ge, X.; Tian, X.; Zhang, Y.; Zhang, J.; Zhang, P. Soy isoflavone: The multipurpose phytochemical (Review). Biomed. Rep. 2013, 1, 697–701. [Google Scholar] [CrossRef] [PubMed]
  6. Gupta, Y.P. Anti-nutritional and toxic factors in food legumes: A review. Plant Foods Hum. Nutr. 1987, 37, 201–228. [Google Scholar] [CrossRef]
  7. Banti, M.; Bajo, W. Review on Nutritional Importance and Anti-nutritional Factors of Legumes. Int. J. Nutr. Food Sci. 2020, 9, 138–149. [Google Scholar] [CrossRef]
  8. Zaynab, M.; Sharif, Y.; Abbas, S.; Afzal, M.Z.; Qasim, M.; Khalofah, A.; Ansari, M.J.; Khan, K.A.; Tao, L.; Li, S. Saponin toxicity as key player in plant defense against pathogens. Toxicon 2021, 193, 21–27. [Google Scholar] [CrossRef] [PubMed]
  9. Hussien, F.M.; Dagnaw, M.M.; Ahmed, A.Y.; Hassen, H.Y. Lathyrism and socioeconomic disparities: A neglected public health problem in Northeast Ethiopia. Am. J. Trop. Med. Hyg. 2021, 104, 1889–1894. [Google Scholar] [CrossRef]
  10. Lambein, F.; Travella, S.; Kuo, Y.H.; Van Montagu, M.; Heijde, M. Grass pea (Lathyrus sativus L.): Orphan crop, nutraceutical or just plain food? Planta 2019, 250, 821–838. [Google Scholar] [CrossRef]
  11. Pavlík, M.; Váňová, M.; Laudová, V.; Harmatha, J. Fungitoxicity of natural heterocycle glucoside vicine obtained from Vicia faba L. against selected microscopic filamentous fungi. Rostl. Vyroba 2002, 48, 543–547. [Google Scholar] [CrossRef]
  12. Björnsdotter, E.; Nadzieja, M.; Chang, W.; Escobar-Herrera, L.; Mancinotti, D.; Angra, D.; Xia, X.; Tacke, R.; Khazaei, H.; Crocoll, C.; et al. VC1 catalyses a key step in the biosynthesis of vicine in faba bean. Nat. Plants 2021, 7, 923–931. [Google Scholar] [CrossRef] [PubMed]
  13. Rizzello, C.G.; Losito, I.; Facchini, L.; Katina, K.; Palmisano, F.; Gobbetti, M.; Coda, R. Degradation of vicine, convicine and their aglycones during fermentation of faba bean flour. Sci. Rep. 2016, 6, 32452. [Google Scholar] [CrossRef]
  14. McMillan, D.C.; Bolchoz, L.J.C.; Jollow, D.J. Favism: Effect of divicine on rat erythrocyte sulfhydryl status, hexose monophosphate shunt activity, morphology, and membrane skeletal proteins. Toxicol. Sci. 2001, 62, 353–359. [Google Scholar] [CrossRef]
  15. Naber, E.C.; Vogt, H.; Harnish, S.; Krieg, R.; Ueberschaer, K.H.; Rauch, H.W. Reproductive performance of hens fed field beans and potential relationships to vicine metabolism. Poult. Sci. 1988, 67, 455–462. [Google Scholar] [CrossRef]
  16. Crépon, K.; Marget, P.; Peyronnet, C.; Carrouée, B.; Arese, P.; Duc, G. Nutritional value of faba bean (Vicia faba L.) seeds for feed and food. F. Crop. Res. 2010, 115, 329–339. [Google Scholar] [CrossRef]
  17. Mejri, F.; Selmi, S.; Martins, A.; Benkhoud, H.; Baati, T.; Chaabane, H.; Njim, L.; Serralheiro, M.L.M.; Rauter, A.P.; Hosni, K. Broad bean (Vicia faba L.) pods: A rich source of bioactive ingredients with antimicrobial, antioxidant, enzyme inhibitory, anti-diabetic and health-promoting properties. Food Funct. 2018, 9, 2051–2069. [Google Scholar] [CrossRef] [PubMed]
  18. Köpke, U.; Nemecek, T. Ecological services of faba bean. F. Crop. Res. 2010, 115, 217–233. [Google Scholar] [CrossRef]
  19. Multari, S.; Stewart, D.; Russell, W.R. Potential of Fava Bean as Future Protein Supply to Partially Replace Meat Intake in the Human Diet. Compr. Rev. Food Sci. Food Saf. 2015, 14, 511–522. [Google Scholar] [CrossRef]
  20. Khazaei, H.; Purves, R.W.; Hughes, J.; Link, W.; O’Sullivan, D.M.; Schulman, A.H.; Björnsdotter, E.; Geu-Flores, F.; Nadzieja, M.; Andersen, S.U.; et al. Eliminating vicine and convicine, the main anti-nutritional factors restricting faba bean usage. Trends Food Sci. Technol. 2019, 91, 549–556. [Google Scholar] [CrossRef]
  21. Pulkkinen, M.; Coda, R.; Lampi, A.M.; Varis, J.; Katina, K.; Piironen, V. Possibilities of reducing amounts of vicine and convicine in faba bean suspensions and sourdoughs. Eur. Food Res. Technol. 2019, 245, 1507–1518. [Google Scholar] [CrossRef]
  22. Larcher, R.; Nardin, T. Suspect screening of glycoalkaloids in plant extracts using neutral loss—High resolution mass spectrometry. J. Chromatogr. A 2019, 1596, 59–68. [Google Scholar] [CrossRef]
  23. Collier, H. The Estimation of Vicine in Fababeans by an Ultraviolet Spectrophotometric Method. Can. Inst. Food Sci. Technol. J. 1976, 9, 155–159. [Google Scholar] [CrossRef]
  24. Marquardt, R.R.; Frohlich, A.A. Rapid reversed-phase high-performance liquid chromatography method for the quantitation of vicine, convicine and related compounds. J. Chromatogr. A 1981, 208, 373–379. [Google Scholar] [CrossRef]
  25. Aydoğan, C. Recent advances and applications in LC-HRMS for food and plant natural products: A critical review. Anal. Bioanal. Chem. 2020, 412, 1973–1991. [Google Scholar] [CrossRef]
  26. Elkins, A.C.; Rochfort, S.J.; Maharjan, P.; Panozzo, J. A Simple High-Throughput Method for the Analysis of Vicine and Convicine in Faba Bean. Molecules 2022, 27, 6288. [Google Scholar] [CrossRef]
  27. Purves, R.W.; Khazaei, H.; Vandenberg, A. Quantification of vicine and convicine in faba bean seeds using hydrophilic interaction liquid chromatography. Food Chem. 2018, 240, 1137–1145. [Google Scholar] [CrossRef] [PubMed]
  28. Purves, R.W.; Khazaei, H.; Vandenberg, A. Toward a high-throughput method for determining vicine and convicine levels in faba bean seeds using flow injection analysis combined with tandem mass spectrometry. Food Chem. 2018, 256, 219–227. [Google Scholar] [CrossRef]
  29. 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. [Google Scholar] [CrossRef] [PubMed]
  30. 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. [Google Scholar] [CrossRef]
  31. Quéméner, B.; Ralet, M. Evidence for linkage position determination in known feruloylated mono- and disaccharides using electrospray ion trap mass spectrometry. J. Mass Spectrom. 2004, 39, 1153–1160. [Google Scholar] [CrossRef] [PubMed]
  32. Kokotou, M.G.; Mantzourani, C.; Babaiti, R.; Kokotos, G. Study of the royal jelly free fatty acids by liquid chromatography-high resolution mass spectrometry (LC-HRMS). Metabolites 2020, 10, 40. [Google Scholar] [CrossRef]
  33. Wojakowska, A.; Piasecka, A.; García-lópez, P.M.; Zamora-natera, F.; Krajewski, P.; Marczak, Ł.; Kachlicki, P.; Stobiecki, M. Phytochemistry Structural analysis and profiling of phenolic secondary metabolites of Mexican lupine species using LC—MS techniques. Phytochemistry 2013, 92, 71–86. [Google Scholar] [CrossRef]
  34. Shi, X.; Yang, W.; Huang, Y.; Hou, J.; Qiu, S.; Yao, C.; Feng, Z.; Wei, W.; Wu, W.; Guo, D. Direct screening of malonylginsenosides from nine Ginseng extracts by an untargeted profiling strategy incorporating in-source collision-induced dissociation, mass tag, and neutral loss scan on a hybrid linear ion-trap/Orbitrap mass spectrometer coupled t. J. Chromatogr. A 2018, 1571, 213–222. [Google Scholar] [CrossRef]
  35. Desmet, S.; Saeys, Y.; Verstaen, K.; Dauwe, R.; Kim, H.; Niculaes, C.; Fukushima, A.; Goeminne, G.; Vanholme, R.; Ralph, J.; et al. Maize specialized metabolome networks reveal organ-preferential mixed glycosides. Comput. Struct. Biotechnol. J. 2021, 19, 1127–1144. [Google Scholar] [CrossRef]
  36. 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. [Google Scholar] [CrossRef]
  37. Schmidt, S.; Zietz, M.; Schreiner, M.; Rohn, S.; Kroh, L.W.; Krumbein, A. Identification of complex, naturally occurring flavonoid glycosides in kale (Brassica oleracea var. sabellica) by high-performance liquid chromatography diode-array detection/electrospray ionization multi-stage mass spectrometry. Rapid Commun. Mass Spectrom. 2022, 24, 2009–2022. [Google Scholar] [CrossRef]
  38. Gould, K.S.; Jay-Allemand, C.; Logan, B.A.; Baissac, Y.; Bidel, L.P.R. When are foliar anthocyanins useful to plants? Re-evaluation of the photoprotection hypothesis using Arabidopsis thaliana mutants that differ in anthocyanin accumulation. Environ. Exp. Bot. 2018, 154, 11–22. [Google Scholar] [CrossRef]
  39. Dambrova, M.; Makrecka-Kuka, M.; Kuka, J.; Vilskersts, R.; Nordberg, D.; Attwood, M.M.; Smesny, S.; Sen, Z.D.; Guo, A.C.; Oler, E.; et al. Acylcarnitines: Nomenclature, Biomarkers, Therapeutic Potential, Drug Targets, and Clinical Trials. Pharmacol. Rev. 2022, 74, 506–551. [Google Scholar] [CrossRef] [PubMed]
  40. Ramsay, G.; Griffiths, D.W. Accumulation of vicine and convicine in Vicia faba and V. narbonensis. Phytochemistry 1996, 42, 63–67. [Google Scholar] [CrossRef]
  41. Duan, S.; Kwon, S.J.; Lim, Y.J.; Gil, C.S.; Jin, C.; Eom, S.H. L-3,4-dihydroxyphenylalanine accumulation in faba bean (Vicia faba L.) tissues during different growth stages. Agronomy 2021, 11, 502. [Google Scholar] [CrossRef]
  42. Griffiths, D.W.; Ramsay, G. The distribution of pyrimidinone glucosides in developing seedlings of Vicia faba and Vicia narbonensis. J. Sci. Food Agric. 1996, 72, 469–475. [Google Scholar] [CrossRef]
  43. Desroches, P.; El Shazly, E.; Mandon, N.; Duc, G.; Huignard, J. Development of Callosobruchus chinensis (L.) and C. maculatus (F.) (Coleoptera: Bruchidae) in seeds of Vicia faba L. differing in their tannin, vicine and convicine contents. J. Stored Prod. Res. 1995, 31, 83–89. [Google Scholar] [CrossRef]
  44. Bjerg, B.; Heide, M.; Nørgaard Knudsen, J.C.; Sørensen, H. Inhibitory effects of convicine, vicine and dopa from Vicia faba on the in vitro growth rates of fungal pathogens. J. Plant Dis. Prot. 1984, 91, 483–487. [Google Scholar]
  45. Backes, A.; Charton, S.; Planchon, S.; Esmaeel, Q.; Sergeant, K.; Hausman, J.-F.; Renaut, J.; Barka, E.A.; Jacquard, C.; Guerriero, G. Gene expression and metabolite analysis in barley inoculated with net blotch fungus and plant growth-promoting rhizobacteria. Plant Physiol. Biochem. 2021, 168, 488–500. [Google Scholar] [CrossRef] [PubMed]
  46. Ruttkies, C.; Schymanski, E.L.; Wolf, S.; Hollender, J.; Neumann, S. MetFrag relaunched: Incorporating strategies beyond in silico fragmentation. J. Cheminform. 2016, 8, 1–16. [Google Scholar] [CrossRef] [PubMed]
  47. Wang, F.; Allen, D.; Tian, S.; Oler, E.; Gautam, V.; Greiner, R.; Metz, T.O.; Wishart, D.S.; Cfm-id, T. CFM-ID 4. 0—A web server for accurate MS-based metabolite identification. Nucleic Acids Res. 2022, 50, 165–174. [Google Scholar] [CrossRef] [PubMed]
  48. Salek, R.M.; Steinbeck, C.; Viant, M.R.; Goodacre, R.; Dunn, W.B. The role of reporting standards for metabolite annotation and identification in metabolomic studies. Gigascience 2013, 2, 13. [Google Scholar] [CrossRef]
Figure 1. Extracted fragment ion chromatogram of the fragment at 141.03 ± 0.015, thus covering the base peak fragment of vicine and convicine. The inset shows the RT range from 5 to 30 min; the peak at 22.38 is cut to illustrate the number of MS/MS spectra containing the fragment at 141.03 ± 0.015. Output generated during the data analysis with Progenesis QI v2.3.
Figure 1. Extracted fragment ion chromatogram of the fragment at 141.03 ± 0.015, thus covering the base peak fragment of vicine and convicine. The inset shows the RT range from 5 to 30 min; the peak at 22.38 is cut to illustrate the number of MS/MS spectra containing the fragment at 141.03 ± 0.015. Output generated during the data analysis with Progenesis QI v2.3.
Molecules 29 01065 g001
Figure 2. The MS/MS spectra of feruloyl convicine and feruloyl vicine indicated are the characteristic fragment ions for the two gluco-aminopyrimidines and feruloyl glucoside showing the attachment of ferulic acid to the glucose. Other major fragments at m/z 193, 175, 160, and 134 correspond to ferulic acid fragments.
Figure 2. The MS/MS spectra of feruloyl convicine and feruloyl vicine indicated are the characteristic fragment ions for the two gluco-aminopyrimidines and feruloyl glucoside showing the attachment of ferulic acid to the glucose. Other major fragments at m/z 193, 175, 160, and 134 correspond to ferulic acid fragments.
Molecules 29 01065 g002
Figure 3. MS/MS spectrum of compound 84 (m/z 486.172 RT 22.49), tentatively identified as convicine + C10H14O3; characteristic losses for convicine are indicated: the loss of 86.03 (C2H2N2O2) from isouramil, the loss of isouramil (−143.03 C4H5N3O3), convicine at m/z 304.074, and isouramil − 2H at m/z 141.016. The highest intensity m/z values that were used in the MZCloud library search and matched to decenedioic acid are encircled in red.
Figure 3. MS/MS spectrum of compound 84 (m/z 486.172 RT 22.49), tentatively identified as convicine + C10H14O3; characteristic losses for convicine are indicated: the loss of 86.03 (C2H2N2O2) from isouramil, the loss of isouramil (−143.03 C4H5N3O3), convicine at m/z 304.074, and isouramil − 2H at m/z 141.016. The highest intensity m/z values that were used in the MZCloud library search and matched to decenedioic acid are encircled in red.
Molecules 29 01065 g003
Figure 4. Spectra resulting in the identification of the MS features (21.94, 772.26) and (22.17, 771.28) as a convicine and a vicine linked by decenedioic acid and two vicines linked by decenedioic acid, respectively. Spectra (A,B) show the total fragmentation spectrum, while spectra (CE) show details of these spectra.
Figure 4. Spectra resulting in the identification of the MS features (21.94, 772.26) and (22.17, 771.28) as a convicine and a vicine linked by decenedioic acid and two vicines linked by decenedioic acid, respectively. Spectra (A,B) show the total fragmentation spectrum, while spectra (CE) show details of these spectra.
Molecules 29 01065 g004
Figure 5. PCA plot calculated based on MS signal intensities of the compounds reported in Table 1; Output generated during the data analysis with Progenesis QI v2.3. The capitals A to E represent the five varieties used in this study for which details are provided in Table 2.
Figure 5. PCA plot calculated based on MS signal intensities of the compounds reported in Table 1; Output generated during the data analysis with Progenesis QI v2.3. The capitals A to E represent the five varieties used in this study for which details are provided in Table 2.
Molecules 29 01065 g005
Table 1. Overview of the 89 (con)vicine derivatives identified in negative ion mode in the five faba bean varieties.
Table 1. Overview of the 89 (con)vicine derivatives identified in negative ion mode in the five faba bean varieties.
#RT am/z bFormula FC cMax dMin dIDΔM eDerivative f
mother compounds and short-chain acid (# C from 3 to 6) derivatives
10.71303.094C10H16N4O77.5B E VC−0.8
20.71304.079C10H15N3O85.7B E CV2.3
31.10466.131C16H25N3O139.6C D CV Hex−1.9
41.14465.145C16H26N4O127.6C D VC Hex−5.8
51.49390.078C13H17N3O11104.3B E malonyl CV−1.7
61.69389.093C13H18N4O1054.5B E malonyl VC−4.6
73.64448.120C16H23N3O12254.0B E hydroxyadipoyl CV−2.9
84.19447.136C16H24N4O1196.6B E hydroxyadipoyl VC−2.1
94.62448.118C16H23N3O12487.5B E hydroxyadipoyl CV−5.9
105.64409.097C16H18N4O93.2B D CV + C6H3NO−7.9nicotinic acid
116.19404.130C15H23N3O104.1A E CV + C5H8O2−3.8hydroxy-valerate
126.42360.104C13H19N3O93.4A D CV + C3H4O−3.2propionic acid
136.54408.114C16H19N5O85.4B E VC + C6H3NO−5.1nicotinic acid
147.01403.145C15H24N4O94.0A E VC + C5H8O2−5.0hydroxy-valerate
158.30403.145C15H24N4O94.9B E VC + C5H8O2−4.5hydroxy-valerate
168.50585.170C23H30N4O149.6A E VC + hydroxybenzoic acid Hex1.6
1711.04424.099C17H19N3O104.3A B CV + hydroxybenzoic acid−2.4
1811.63374.120C14H21N3O93.8C E butyryl CV−0.5
1912.02423.112C17H20N4O92.5A E VC + hydroxybenzoic acid−8.1
2012.76373.136C14H22N4O84.0C E butyryl VC−1.1
2115.77386.119C15H21N3O93.9A D CV + C5H6O−4.2pentenoic acid
2217.44385.134C15H22N4O82.7B E VC + C5H6O−5.2pentenoic acid
2317.21388.136C15H23N3O92.8C D (iso)valeryl CV−0.4
2418.23387.152C15H24N4O83.1C E (iso)valeryl VC0.7
2519.48387.150C15H24N4O83.5B C (iso)valeryl VC−4.6
derivatives with phenylpropanoids and guaiacylglycerol
2614.05612.167C25H31N3O158.8A D CV coumaroyl Hex−2.0
2714.59611.181C25H32N4O1420.1B D VC coumaroyl Hex−4.7
2815.57466.109C19H21N3O114.1A D CV caffeic acid−2.4
2916.43465.124C19H22N4O104.4B E VC caffeic acid−4.1
3016.70642.178C26H33N3O168.7A D CV feruloyl Hex −1.2
3117.36641.193C26H34N4O159.9A D VC feruloyl Hex −3.5
3217.99672.189C27H35N3O173.9C D CV sinapoyl Hex−1.3
3318.42671.203C27H36N4O162.4B E VC sinapoyl Hex−4.2
3418.42450.113C19H21N3O108.0C E CV coumaric acid−4.8
3518.89554.162C23H29N3O133.9A D CV + C13H14O5−1.0glyceryl ferulate
3619.25553.177C23H30N4O123.3A D VC + C13H14O5−3.0glyceryl ferulate
3719.48500.151C20H27N3O126.8E D CV + C10H12O4−1.7guaiacylglycerol
3819.60480.126C20H23N3O112.9A E CV ferulic acid0.1
3919.99479.141C20H24N4O107.4A E VC ferulic acid−2.9
4020.03499.166C20H28N4O116.7B D VC + C10H12O4−4.2guaiacylglycerol
4120.14510.135C21H25N3O123.5C E CV sinapinic acid−3.9
4220.53509.150C21H26N4O115.4B E VC sinapinic acid−4.0
4320.97676.197C30H35N3O154.2A E feruloyl CV + C10H12O4−3.1guaiacylglycerol
4421.32675.213C30H36N4O147.8A E feruloyl VC + C10H12O4−3.3guaiacylglycerol
4522.571142.329C48H61N3O299.1A E CV diferuloyl triHex−2.4
4622.60980.278C42H51N3O248.1A E CV diferuloyl diHex−1.3
4722.761141.345C48H62N4O285.1A E VC diferuloyl triHex−2.3
4822.84979.292C42H52N4O2312.9A E VC diferuloyl diHex−2.9
4922.96818.224C36H41N3O1916.9A E CV diferuloyl Hex−2.8
5023.66817.236C36H42N4O182.8B E VC diferuloyl Hex−7.5
5125.041156.323C52H59N3O273.3A D CV triferuloyl diHex−2.8
5225.231155.337C52H60N4O263.7A E VC triferuloyl diHex−4.2
Medium-chain acid (# C above 6) derivatives
5315.61446.140C17H25N3O113.6A D CV + C7H10O3−3.7pimelic acid
5416.62445.156C17H26N4O105.1B E VC + C7H10O3−4.7pimelic acid
5518.07732.234C27H39N7O1710.5B E CV + C7H10O3 VC1.7pimelic acid
5619.32789.247C30H42N6O194.0C D CV + C10H14O4 CV4.3hydroxy decenedioic acid
5719.52746.248C28H41N7O172.0A D VC + C8H12O3 CV−1.5suberic acid
5819.83788.261C30H43N7O1831.2C D CV + C10H14O4 VC2.9hydroxy decenedioic acid
5919.95502.168C20H29N3O126.6C D CV + C10H14O40.3hydroxy decenedioic acid
6020.30504.183C20H31N3O1216.1C D CV + C10H16O4−1.2hydroxy sebacic acid
6120.50501.184C20H30N4O119.4C D VC + C10H14O4−0.4hydroxy decenedioic acid
6220.61572.211C24H35N3O133.2A B CV + C14H20O51.8Dihydroxy-tetra-decadienedioic acid
6320.77503.199C20H32N4O1115.9C D VC + C10H16O4−0.7hydroxy sebacic acid
6420.85472.156C19H27N3O114.9A D CV + C9H12O3−2.1nonenedioic acid
6520.97814.274C32H45N7O1824.4A E VC + C12H16O4 CV−1.5hydroxy-dodecadienedioic acid
6621.08571.224C24H36N4O122.8A B VC + C14H20O5−2.7Dihydroxy-tetra-decadienedioic acid
6721.24474.173C19H29N3O1149.6C E CV + C9H14O30.3azelaic acid
6821.28813.288C32H46N8O1735.8A E VC + C12H16O4 VC−3.7hydroxy-dodecadienedioic acid
6921.43471.171C19H28N4O1013.8A D VC + C9H12O3−5.0nonenedioic acid
7021.43773.244C30H42N6O186.5A D CV + C10H14O3 CV −5.3decenedioic acid
7121.47528.180C22H31N3O125.6C E CV + C12H16O4−6.8hydroxy-dodecadienedioic acid
7221.63859.330C34H52N8O1810.7B E VC + C14H22O5 VC−3.1dihydroxy-tetra-decenedioic acid
7321.67775.266C30H44N6O184.1A B CV C10H16O3 CV2.7sebacic acid
7421.71530.198C22H33N3O122.5C E CV + C12H18O4−2.0hydroxy-dodecenedioic acid
7521.71473.188C19H30N4O102.5C E VC + C9H14O3−3.5azelaic acid
7621.94772.265C30H43N7O172.7A D VC + C10H14O3 CV 0.4decenedioic acid
7722.02527.198C22H32N4O113.6A E VC + C12H16O4−3.3hydroxy-dodecadienedioic acid
7822.14574.224C24H37N3O132.9A E CV + C14H22O5−2.8dihydroxy-tetra-decenedioic acid
7922.17774.281C30H45N7O173.2A D VC + C10H16O3 CV1.0sebacic acid
8022.17771.279C30H44N8O163.0A E VC + C10H14O3 VC −2.1decenedioic acid
8122.17529.213C22H34N4O113.3C E VC + C12H18O4−3.8hydroxy-dodecenedioic acid
8222.37573.239C24H38N4O123.5A E VC + C14H22O5−3.7dihydroxy-tetra-decenedioic acid
8322.41773.295C30H46N8O164.1A D VC + C10H16O3 VC−1.2sebacic acid
8422.49486.172C20H29N3O113.0A B CV + C10H14O3−1.3decenedioic acid
8522.92485.186C20H30N4O102.1A E VC + C10H14O3−5.5decenedioic acid
8623.07488.187C20H31N3O115.3A D CV + C10H16O3−2.2sebacic acid
8723.27485.187C20H30N4O103.3A C VC + C10H14O3−3.2decenedioic acid
8823.58487.203C20H32N4O105.1A D VC + C10H16O3−3.2sebacic acid
a: retention time in minutes; b: based on the output of progenesis QI, the given m/z value corresponds to the experimental m/z of the [M − H] ion; c: maximum fold change based on the average normalized intensity of the different varieties; d: variety with the highest or lowest average MS signal intensity, details on the varieties A to E are given in Table 2; e: mass error in ppm between the experimental m/z and the calculated m/z for the chemical formula in negative mode; f: proposed name for the added group, determined based on mass/formula, the literature, and MS/MS fragments; VC: vicine; CV: convicine; Hex: hexose.
Table 2. The five faba bean varieties used in this study.
Table 2. The five faba bean varieties used in this study.
CodeGenotypeRegistration GERColorVicineType
ATaifun2011ZtHVCSpring
BAugusta2018BbHVCWinter
CTrumpet2017BbHVCSpring
DAllison2019BbLVCSpring
EBreeding line ZtLVCSpring
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

Sergeant, K.; Goertz, S.; Halime, S.; Tietgen, H.; Heidt, H.; Minestrini, M.; Jacquard, C.; Zimmer, S.; Renaut, J. Exploration of the Diversity of Vicine and Convicine Derivatives in Faba Bean (Vicia faba L.) Cultivars: Insights from LC-MS/MS Spectra. Molecules 2024, 29, 1065. https://doi.org/10.3390/molecules29051065

AMA Style

Sergeant K, Goertz S, Halime S, Tietgen H, Heidt H, Minestrini M, Jacquard C, Zimmer S, Renaut J. Exploration of the Diversity of Vicine and Convicine Derivatives in Faba Bean (Vicia faba L.) Cultivars: Insights from LC-MS/MS Spectra. Molecules. 2024; 29(5):1065. https://doi.org/10.3390/molecules29051065

Chicago/Turabian Style

Sergeant, Kjell, Simon Goertz, Salma Halime, Hanna Tietgen, Hanna Heidt, Martina Minestrini, Cédric Jacquard, Stephanie Zimmer, and Jenny Renaut. 2024. "Exploration of the Diversity of Vicine and Convicine Derivatives in Faba Bean (Vicia faba L.) Cultivars: Insights from LC-MS/MS Spectra" Molecules 29, no. 5: 1065. https://doi.org/10.3390/molecules29051065

APA Style

Sergeant, K., Goertz, S., Halime, S., Tietgen, H., Heidt, H., Minestrini, M., Jacquard, C., Zimmer, S., & Renaut, J. (2024). Exploration of the Diversity of Vicine and Convicine Derivatives in Faba Bean (Vicia faba L.) Cultivars: Insights from LC-MS/MS Spectra. Molecules, 29(5), 1065. https://doi.org/10.3390/molecules29051065

Article Metrics

Back to TopTop