Abstract

An NMR-based metabolomic approach has been applied to analyse seven aromatic Mediterranean plant species used in traditional cuisine. Based on the ethnobotanical use of these plants, the approach has been employed in order to study the metabolic changes during different seasons. Primary and secondary metabolites have been detected and quantified. Flavonoids (apigenin, quercetin, and kaempferol derivatives) and phenylpropanoid derivatives (e.g., chlorogenic and rosmarinic acid) are the main identified polyphenols. The richness in these metabolites could explain the biological properties ascribed to these plant species.

1. Introduction

Aromatic plants are widespread throughout the world and they are extensively added to different food preparations. The use of these plants is very popular and has a long tradition in Mediterranean area [1].

Plants in general have been shown to produce a wide range of chemicals, traditionally categorized into primary and secondary metabolites. For the sake of simplicity, primary metabolites can be thought of to serve nutritional purposes, while secondary metabolites are required by plants as weapons against competitors, herbivores or pathogens, and so forth [2]. However, both classes of metabolites are important for the plant itself but also for their actions on plant consumers and mainly in case of edible plants.

Many aromatic plants are added to foods and eaten. The whole plants or one or some of their components are used as, for example, food preservatives, flavour, and additives. Nevertheless, it has been shown that chemical composition of plants is highly variable along the year. Several analytical techniques are available for studying plant metabolites’ content. Most of them are targeted techniques as an a priori knowledge on the metabolites to be analysed is required [2, 3]. Furthermore, for aromatic plants, a great effort has been devoted to the study of essential oils [2, 4].

Given this background, a wider knowledge about their whole metabolite content is needed. To this end, a very powerful approach is metabolomics, the comprehensive analysis of the set of low molecular weight compounds of a biological system under a given condition [5]. Analogously, other related approaches, like metabolic profiling [6], could be used.

In particular, NMR-based metabolomics has been shown to be very useful due to its untargeted and unbiased features [7]. Furthermore, it is highly reproducible, it allows the contemporary identification and quantification of a large number of compounds and needs short times of analysis (including the extraction procedures) [6]. The only limitation of NMR is its low sensitivity when compared to mass spectrometry, although sensitivity has been drastically increased with recent advances like higher magnetic fields and the introduction of microcryoprobes [5, 8]. On the other hand, NMR allows the identification of unknown compounds in the analysed mixtures, as it gives important structural information [8].

In order to demonstrate the potentiality of this approach, it has been applied to seven aromatic plant species characteristic of Mediterranean garrigue: the metabolites’ content of these plants has been determined and the seasonality of their accumulation has been studied.

2. Materials and Methods

2.1. Plant Material Sampling and Processing

Seven plant species (Table 1) were collected in a garrigue on the calcareous hills of Durazzano, (41°3′N, 14°27′E; southern Italy) in winter (February 2012), spring (May 2012), summer (July 2012), and autumn (October 2012). The plants were selected based on their occurrence in the study site. Origanum vulgare samples were not available in autumn. Plant leaf samples were collected in the field always at the same moment of the day, in order to minimize differences due to metabolites changing based on circadian clock.

Three leaf samples (biological replicates) of each plant species were harvested and immediately frozen in liquid in order to avoid unwanted enzymatic reactions and stored at −80°C up to the freeze drying process. Once freeze dried they were powdered in liquid nitrogen and stored at −20°C. Each sample was extracted and analysed by NMR.

Voucher specimens for all the plant species were deposited at the herbarium of the Second University of Naples (Table 1).

2.2. Metabolomics Analysis

Freeze-dried plant material (50 mg) was transferred to a 2 mL microtube. NMR samples were prepared in a mixture of phosphate buffer (Fluka Chemika; 90 mM; pH 6.0) in D2O (Cambridge Isotope Laboratories) containing 0.1% w/w trimethylsilylpropionic-2,2,3,3- acid sodium salt (TMSP, Sigma-Aldrich) and methanol- (Sigma-Aldrich). A volume of 1.5 mL of phosphate buffer and methanol- (1 : 1) was added to the plant samples. The mixture was vortexed at room temperature for 1 min, ultrasonicated (Elma Transonic Digitals) for 40 min, and centrifuged (Beckman Allegra 64R) at 13000 rpm for 10 min. An aliquot of 0.6 mL was transferred to an NMR tube and analysed by NMR [9]. NMR spectra were recorded at 25°C on a 300.03 MHz for 1H and 75.45 MHz for 13C on a Varian Mercury Plus 300 Fourier transform NMR. CD3OD was used as the internal lock. Each 1H NMR spectrum consisted of 256 scans with the following parameters: 0.16 Hz/point, acquisition time (AQ) = 1.0 s, relaxation delay (RD) = 1.5 s, 90° pulse width (PW) = 13.8 μs. A presaturation sequence was used to suppress the residual H2O signal. Free induction decays (FIDs) were Fourier transformed with LB = 0.3 Hz and the resulting spectra were manually phased and baseline-corrected and calibrated to TMSP at 0.0 ppm, using 1H NMR processor (MestReNova, version 8.0.2).

2.3. Quantitative Analysis

The main metabolites identified in plant extracts were analyzed by quantitative analysis. 1H-NMR spectra were bucketed, reducing it to integral segments with a width of 0.02 ppm with ACDLABS 12.0 1H-NMR processor (ACDLABS 12.0, Toronto, Canada). Spectra were scaled to the internal standard (whose area, from −0.01 to 0.01 ppm, was set equal to 1). For each metabolite, buckets corresponding to nonoverlapping signals were used to calculate the relative amount as follows: where is the metabolite signal area, but it is also equal to the signal area/standard area ratio, as standard area is equal to 1; is a constant equal to 9 (the number of protons responsible for the signal between −0.01 and 0.01 ppm) and is the number of protons of the metabolite signal area [10].

3. Results and Discussion

Recent research has shown culinary herbs and spices as a source of bioactive compounds [11]. Although most of them have been extensively studied for their essential oil composition, far less information is available on their polar and semipolar chemical composition.

Herewith, seven Mediterranean plants (Calamintha nepeta, Helichrysum italicum, Foeniculum vulgare, Micromeria graeca, Origanum vulgare, Satureja montana, and Thymus longicaulis) have been studied for their metabolite content by NMR. The identification of metabolites was carried out by comparing NMR data with an in-house library, with databases [12], and with some literature data [10, 1315]. 1H-NMR data and extract composition are given in Table 2 and spectra are shown in Figure 1.

Primary metabolites were easily identified based on data extensively reported in literature of spectra acquired in the same solvent mixture [10, 1215]. Among free amino acids, alanine was observed in all of the plants, while threonine was only detected in Micromeria graeca and Foeniculum vulgare.

The sugar content was highly variable, with glucose and sucrose as the main free carbohydrates detected.

Finally, some organic acids were identified. Quinic acid was present in all of the plants but Helichrysum italicum and Satureja montana, while malic acid was clearly detected in all of the plants.

Concerning the secondary metabolite content, the analysed Lamiaceae plants were all characterized by the presence of high amounts of rosmarinic acid (with the exception of C. nepeta), along with analogous compounds. Caffeic acid was identified based on comparison of NMR data with the literature [10, 15] and confirmed by comparison with NMR spectra of an in-house library. Rosmarinic acid was identified based on the comparison with already reported data [10] and the structure was confirmed by 2D NMR analysis. Indeed, the olefinic proton at δ 7.50 (H7) (showing HSQC correlation with the carbon at δ 145.9) and that at δ 6.30 (showing HSQC correlation with the carbon at δ 114.3) showed long range correlations with a carbon at δ 168.4 (C9). This carbon was in turn correlated with the proton at δ 5.02 (H8′), confirming the linkage between a caffeoyl moiety and the 3,4-dihydroxyphenyl lactic acid moiety. Furthermore, the former was identified based on the long range correlation of the H7 olefinic proton with the aromatic carbon at δ 126.4 (C1), showing further correlations with the signals belonging to an ortho/para trisubstituted aromatic ring (Table 2). The latter was identified as follows: the proton H8′ showed correlations with a carboxylic carbon at δ 176.5 and with a methylene carbon at δ 36.9 (C7′), showing HSQC correlations with the diastereotopic protons H7′ (Table 2). The proton H8′ also showed long range correlation with a quaternary aromatic carbon at δ 130.0 (C1′), in turn correlated with the signals belonging to a second ortho/para trisubstituted aromatic ring (Table 2).

Some phenylpropanoids in the extracts were not definitively characterized, inasmuch as, based on their scarce abundance and/or strong signal overlapping, they did not show clear correlations in 2D NMR spectra. However, the characteristic signals and correlations of the trans-propenylic chain suggested their presence. Indeed, correlations were observed in the HSQC, among the olefinic signals with carbons at 140–145 ppm (for the proton at lower fields) and at 114–120 ppm (for the proton at higher fields) and long range correlations were shown with carbon resonances attributable to ester carboxyl carbons and with quaternary aromatic carbons.

Calamintha nepeta extracts were also rich in several flavonoids and phenylpropanoids. Unfortunately, it was not possible to definitely characterize these compounds, but all of the flavonoids were identified as apigenin derivatives (Table 2). Indeed, several sets of resonances attributable to meta coupled protons H7/H8 (ring A), to proton H3 and to B ring ortho coupled protons were detected. Interestingly, the compounds, probably characterized by a different degree of glycosylation, showed a peculiar distribution along the seasons. Two apigenin derivatives were detected in spring and autumn (apigenin derivatives 1 and 2), while two different couples of these compounds were detected in summer (apigenin derivatives 3 and 4) and winter (apigenin derivatives 5 and 6) samples. Moreover, apigenin derivatives 5 and 6 were detected only in winter also in Satureja montana.

Analogously, as shown in Table 2, the presence of some phenylpropanoids was strongly dependent on the collection season (Table 2): phenylpropanoid 2 was detected only in winter, phenylpropanoid 5 only in spring, phenylpropanoid 6 only in summer, and phenylpropanoid 7 only in autumn samples.

The most stable metabolome along the seasons was detected for Thymus longicaulis while Micromeria graeca and Origanum vulgare only changed for some metabolites. However, differences in the amounts of the compounds were observed. Indeed, for all the Lamiaceae plants, a higher amount of aromatic compounds (Table 2) was observed in spring and summer samples compared to autumn and winter samples. Helichrysum italicum and Foeniculum vulgare extracts showed an analogous behaviour, with changes of metabolites mainly on the quantitative point of view.

Helichrysum italicum extracts, besides chlorogenic acids, also showed signals attributable to a 3-hydroxybenzofuran and an isobenzofuranone derivative. Chlorogenic, neochlorogenic, and dicaffeoylquinic acids were identified based on comparison of 1H-NMR data with the literature [10, 15] and with the in-house library. The caffeoyl moiety was clearly identified based on 1D and 2D NMR data and the linkage(s) with the quinic acid moiety was confirmed by the correlation observed in the long range spectrum. The 3-hydroxybenzofuran and isobenzofuranone derivatives were identified based on comparison with the NMR spectrum of the compound previously isolated [16].

Finally, Foeniculum vulgare was characterized by chlorogenic acids and flavonoids, identified, based on 1H-NMR data as kaempferol and quercetin [17]. Chlorogenic acid was reported for the first time from this species, to the best of our knowledge.

The identification of water soluble compounds in these plants is very important as most of them are added to dishes; hence they might be eaten or however they could release bioactive compounds into food. In this framework, it is worth to underlining that first of all the health promoting capacity of bioactive compounds could be dependent on synergisms. Secondly, these plants could also contain potential toxic compounds, as reported for several Lamiaceae [17]. Nevertheless, the role in nutrition of primary metabolites is often disregarded [18]. These considerations raise the attention to the need for a comprehensive profiling of their metabolites.

Furthermore, the NMR-based metabolomic approach here proposed, due to the short time of analyses and to lower costs compared to other analytical methods, was very useful for the study of seasonal variation of metabolites.

Indeed, although it is clear that secondary metabolites show peculiar trends, only fragmentary information is available, mainly because of the used approaches.

The most common and easiest procedure was to perform the phytochemical study on samples collected in a specific time of the year and then compare the other months (or seasons) by setting up a series of target analyses, often by HPLC [19, 20] or by less time and resource consuming approaches based on colorimetric assays [21].

The improvement of analytical methods, and especially the availability of a high-throughput approach like metabolomics [5], gives the chance to further explore the issue of seasonality [22] and to thoroughly study these changes.

Concerning the metabolites identified within the extracts, it is important to underline their biological activity.

Phenols are well known for their antioxidant activity [11]. As the studied plants are all very rich in phenolics, they have a great antioxidant potential. This could also explain the use of some of these herbs as natural food preservatives (Table 1).

Among the detected compounds, rosmarinic acid is the most widespread and the most abundant. A plethora of biological activities has been attributed to this compound, among them: adstringent, antioxidative, anti-inflammatory, antimutagen, antibacterial, and antiviral [23].

Many other phenylpropanoids and flavonoids were also detected. Several properties have been reported for these compounds, such as anti-inflammatory, antimicrobial, and antitumor activity [24].

The richness in these compounds of the studied plants supports their traditional uses.

However, the study evidenced that qualitative and quantitative variations of metabolites are observed along the year.

This is, to the best of our knowledge, the first report of metabolite content and seasonal qualitative and quantitative variation of this set of food spices. Furthermore, it is evidenced that a higher content of compounds known for their health promoting capacities can be found in spring and summer samples of all the analysed species, although season-specific compounds were also detected.

4. Conclusions

NMR-based metabolomics has been applied to the study of chemical composition of selected aromatic plants of Mediterranean vegetation.

The method allowed determining the chemical composition of plant extracts in terms of primary and secondary metabolites. The abundance of aromatic secondary metabolites suggested that the traditional uses of these plants might be supported by their chemical composition.

Furthermore, the seasonality of the accumulation of these metabolites was studied.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors would like to thank the anonymous reviewers for their valuable suggestions.