Next Article in Journal
Identification of Interleukin-8-Reducing Lead Compounds Based on SAR Studies on Dihydrochalcone-Related Compounds in Human Gingival Fibroblasts (HGF-1 cells) In Vitro
Previous Article in Journal
Effect of Terpenoids and Flavonoids Isolated from Baccharis conferta Kunth on TPA-Induced Ear Edema in Mice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Altitude on Phytochemical Composition of Hemp Inflorescence: A Metabolomic Approach

1
Centre of Applied Studies for the Sustainable Management and Protection of Mountain Areas (CRC Ge.S.Di.Mont.), University of Milan, Via Morino 8, 25048 Edolo (BS), Italy
2
Department of Agricultural and Environmental Sciences—Production, Landscape, Agroenergy (DISAA), Via Celoria 2, 20133 Milan, Italy
*
Author to whom correspondence should be addressed.
Molecules 2020, 25(6), 1381; https://doi.org/10.3390/molecules25061381
Submission received: 3 March 2020 / Revised: 16 March 2020 / Accepted: 17 March 2020 / Published: 18 March 2020

Abstract

:
The phytochemical profiling of hemp inflorescences of clonal plants growing in different conditions related to altitude was investigated. Four strains of industrial hemp (Cannabis sativa L., family Cannabaceae) of Kompolti variety were selected and cloned to provide genetically uniform material for analyses of secondary metabolites (terpenes, cannabinoids, and flavonoids) at two different elevations: mountain (Alagna Valsesia 1200 m ASL) and plains (Vercelli Province 130 m ASL). Environmental conditions influenced by elevation have proven to be important factors inducing variations in hemp inflorescences’ secondary metabolite composition. In fact, all plants grown at altitude exhibited a higher total amount of terpenes when compared with plains counterparts, with β-Myrcene, trans-Caryophyllene and α-Humulene as the main contributors. A metabolomic, un-targeted approach performed by HPLC-Q-Exactive-Orbitrap®-MS platform with subsequent data processing performed by Compound Discoverer™ software, was crucial for the appropriate recognition of many metabolites, clearly distinguishing mountain from plains specimens. Cannabidiolic acid CBDA was the most abundant phytocannabinoid, with significantly higher concentrations in the mountain samples. The metabolic pathway of CBGA (considered as the progenitor/precursor of all cannabinoids) was also activated towards the production of CBCA, which occurs in considerably 3 times higher quantities than in the clones grown at high altitude. Isoprenoid flavones (Cannaflavins A, B, and C) were correspondingly upregulated in mountain samples, while apigenin turned out to be more abundant in plains samples. The possibility to use hemp inflorescences in pharmaceutical/nutraceutical applications opens new challenges to understand how hemp crops respond in terms of secondary metabolite production in various environments. In this regard, our results with the applied analytical strategy may constitute an effective way of phytochemical profiling hemp inflorescences.

1. Introduction

Hemp (cannabis, Cannabis sativa L., family: Cannabaceae) was frequently cultivated in the past, but its agronomic expansion was interrupted in the beginning of the 1950s for many reasons, one of them being the presence of psychoactive substance Δ-9-tetrahydrocannabinol (THC) that is produced by some hemp varieties. Nowadays, this has been partly obliterated and the European Union permits the cultivation of hemp with THC content of less than 0.20% [1]. For example, in Italy, regulation n°242/2017 [2], delineates the conditions for hemp production, its commercialization, and utilization for specific industrial purposes. Therefore, different genotypes have been registered, along with standardized cultivation procedures [3]. Many ecologically, agronomically, and pharmaceutically positive properties qualify this multifunctional crop as an opportune raw material for various traditional (fiber, food, oil, medicine) or innovative industrial applications (new biomaterials and biofuels) [4,5].
The European Union has regulated commercial production and distribution of more than 70 hemp varieties [6], among which the Kompolti variety is one of the oldest. This hemp variety is frequently cultivated in continental European countries. This well-known dioecious variety was developed in Hungary as the first hybrid breed [7] in order to produce seeds for oil production, but also its stalks have been largely exploited. The particularity of this strain is its growing dynamics: it needs a whopping 6-month flowering period that usually finishes in October. Its macro/microscopical botanical aspects, geographical distribution, and agricultural status have been comprehensively studied, but surprisingly, no information is currently available on the content of the high added-value bioactive substances that are characteristic of its flowers.
The C. sativa inflorescence contains a number of chemically active compounds, such as cannabinoids, terpenoids, and flavonoids. The most important secondary metabolites are phytocannabinoids whose acidic forms are exclusively biosynthesized in the glandular trichomes, which are abundant on the surface of the female flowers [8]. Inflorescences of industrial hemp varieties, Kompolti included, are particularly rich in cannabidiolic acid (CBDA) that is disposed to the spontaneous decarboxylation to cannabidiol (CBD) under favorable environmental/conservational circumstances such as warming and light [9]. CBD is responsible for a variety of pharmacological actions that may have some remarkable applications, but unlike THC, CBD does not possess any psychoactive effects [10]. CBD has been studied intensively over the past decade due to its biomedical relevance. Several studies suggest that CBD can be effective in treating epilepsy and other neuropsychiatric disorders, including anxiety and schizophrenia [11]. This is the reason why CBD dietary supplements obtained from different industrial C. sativa chemotypes have become particularly widespread [12]. On the other hand, it is important to evaluate the main factors that determine CBD production in hemp inflorescence—if it depends on the genetic predisposition of hemp chemovar, or its production is conditioned by the environment, in particularly the geographical position where plant is bred [13].
Although CBD and THC are the crucial phytocannabinoids, hemp trichomes themselves are capable of generating a whole series of acidic/decarboxylated phytocannabinoids: about 120 have been isolated to date [8]. Based on their appearance in the metabolic pathway that involves their formation, all phytocannabinoids are categorized into 11 subclasses [9], where the central position belongs to cannabigerolic acid that is synthesized from geranyl diphosphate and olivetolic acid [14]. CBGA further provides tetrahydrocannabinolic acid (THCA), CBDA, and cannabichromenic acid (CBCA). Cannabigerol (CBG), THC, CBD, and cannabichromene (CBC) are corresponding neutral equivalents. Other phytocannabinoids detected in hemp inflorescences comprise main oxidation products of THCA and CBDA: cannabinol (CBN) and cannabinolic acid (CBNA) obtained from THC(A), and cannabielsoin (CBE) and cannabielsoinic acid that derive from CBDA. The “cannabivarin” group, commonly following the above-mentioned ones, is produced from condensation of geranyl diphosphate with divarinic acid, which results in a propyl instead of the pentyl side-chain.
Even though the attention of the scientific community has been focused on major phytocannabinoids, the phytochemical characterization of C. sativa highlights the presence of various non-cannabinoids constituents including flavonoids [15]. Their characterization is scarce and random, especially when the inflorescences of industrial hemp are concerned [13,16,17]. In any case, the characterization of the cannaflavones—compounds isolated exclusively from hemp—needs further elucidation, mainly when the inflorescences are concerned.
One non-phytocannabinoid category of bioactive secondary metabolites that is studied in much more detail is the terpene family. They represent volatile components that has been claimed to have a synergic action with cannabinoids [18]. Many different monoterpenes and sesquiterpenes are important components of C. sativa essential oils [19,20], as they define some of the unique organoleptic properties and may also influence nutraceutical potentials of different hemp strains and varieties [21].
Selecting a genotype appropriate for a particular end-use application that is adaptable to an environment is of principal importance for successful hemp cultivation. Hemp is a plant adaptable to various growing and ecological conditions, but there are no data available in the literature that concern the differences that may arise from cultivating the same chemovar contemporarily in the plains and mountain habitats. Different hemp varieties cultivated at high altitudes showed a characteristic phytochemical and ecological behavior [13].
In this research, our aim was to study whether two very different ecological environments (mountain and plains) would have a significant impact on phytocannabinoids qualitative and quantitative content, flavonoids presence, and terpenoids profile, in order to study the plant phytochemical behavior and its potential to provide nutraceutical substances.
Four strains of industrial hemp (Kompolti) were selected and cloned to provide genetically uniform material for analyses of secondary metabolites (cannabinoids, terpenes, and flavonoids) of clones of the same plant grown at different elevation in two sites representative of lowland (Vercelli Province 130 m ASL) and mountain (Alagna Valsesia 1200 m ASL) during the growing season 2018.

2. Results and Discussion

Multi-targeted applications of industrial hemp with the environmental benefits related to its cultivation, have raised interest in its production. Special attention has been paid to C. sativa inflorescences that represent a promising added-value product with remarkable pharmacological and nutraceutical effects [13,16,22,23]. The phytochemical composition of inflorescences has been studied intensively, but there is not substantial information that regards the differences that may rise due to geographical/microclimate factors. In this experiment, raw inflorescences material obtained from plants cultivated in the Italian Alps at two different elevations was evaluated. Mountain (M) region was located in the commonality of Alagna Valsesia (1200 m ASL), whereas plains (P) cultivation was performed in the Province of Vercelli (130 m ASL).

2.1. Terpenoids Profile Estimated by HS-SPME- GC-MS Analytical Procedure

Monoterpenes, diterpenes, triterpenes, and sesquiterpenes are important components of the C. sativa resin responsible for its unique aromatic properties [20,21]. Considering the terpenoids fraction characterized by high-volatile features, the headspace solid-phase microextraction-gas chromatographic mass-spectrometric (HS-SPME-GC-MS) analytical approach presents the best methodology for their comprehensive profiling, as was confirmed by our recently published studies [12,24,25]. SPME is a simple and fast modern tool used to characterize the volatile fraction of secondary metabolites of different parts of plants [26,27] and animals [28].
Complete data concerning the terpenes fingerprint in the two groups of Kompolti inflorescences are summarized and reported in Table 1 It was possible to define 20 compounds that belong to the mono/di/tri terpenes and 21 sesquiterpenes. Our results are qualitatively comparable with those reported by others [19,20,22,29]. The most remarkable aspect is that four M clones expressed significantly higher amounts of both terpenoids subgroups, with high variation in the individual quantitative profile. This is the reason why the results are elaborated by statistical approach that consisted of the comparation of relative amounts of each compound between mountain and corresponding plain clone. Generally, the predominant monoterpene was β-myrcene, followed by both α-/β-pinene and limonene, although without uniformity between four plants (Table 1). For example, the plant M2 was particularly rich in β-myrcene, followed by β-pinene and limonene, but not α-pinene. On the other hand, the plant 1 (both M and P samples) expressed its specificity in the accumulation of β-ocimene, while others were very poor in its presence. Also, plant 1 contained the oxygenated terpene 4,8-epoxy-p-menth-1-ene that was completely absent from other inflorescences.
Geographic origin, accompanied by environmental conditions, turned out to be an important variable that determined the sesquiterpenes’ quantitative characteristics: mountain plants exhibited higher total amounts than plains counterparts, with trans-caryophyllene and α-humulene as the main contributors (Table 1). Those two compounds are typical constituents of C. sativa essential oil [30], and also here showed a quite stable relevance. Selina-3,7(11)-diene was detected in a moderate amount in all mountain samples, but its presence was not detected in two plains clones. α-ylangene, α-bergamotene and β-farnesene expressed highly inconsistent trends. For example, P1 samples were particularly abundant in β-farnesene, while in plant 3 (both P and M samples) it was absent. All samples contained longicyclene, the sesquiterpene rarely identified in C. sativa inflorescences, with the exception of one drug-type chemovar [31].
In all cases, the qualitative and quantitative differences observed in the chemical profile of terpene fractions were conditioned by many factors such as, hemp variety, cultivation and environmental conditions, harvest time and post-harvest conditions, storage and drying of raw plant, and extraction procedure applied. Despite the fact that within each group, plants were grown under identical conditions and treated in the same way, it remains an open question why each of them had its specific terpenoids fingerprint.

2.2. Quantification of Cannabinoids by HPLC-Q-Exactive-Orbitrap®-MS Analysis

Based on the available commercial cannabinoids standards, quantification of inflorescence extracts was performed by applying our validated method [13,24,25], as explained in detail in the material and methods section. Quantitative data related to the analysis of the content of phytocannabinoids in the two inflorescences groups performed by means of the HPLC high-resolution mass spectrometry (HRMS, Q-Exactive-Orbitrap®-MS) are shown in Table 2. The quantification was performed for CBD (cannabidiol), Δ9-THC (delta-9- tetrahydrocannabinol), CBN (cannabinol), CBC (cannabichromene), CBG (cannabigerol), CBDV (cannabidivarin), Δ9-THCV (delta-9-tetrahydrocannabivarin) and the acid forms CBDA (cannabidiolic acid), Δ9-THCA (delta-9- tetrahydrocannabinolic acid), CBNA (cannabinolic acid), CBCA (cannabichromenic acid), CBGA (cannabigerolic acid), CBDVA (cannabidivarin acid) and Δ9-THCVA (delta-9-tetrahydrocannabivarin acid). In contrast to what was found for volatile terpenoids, the results obtained for cannabinoids were more uniform with respect to the cultivation site. It was therefore feasible to perform paired statistical evaluations combining the specimens from both locations. Since the Kompolti chemovar belongs to the fiber-type hemp, it is not surprising that CBDA was the most abundant phytocannabinoid with significantly higher concentration in the mountain than in the plains. Those values (both for M and P group) are higher than those recently reported by [32,33]. The THC and THCA content were below legal limits, but the occurrence of traces of CBNA —THCA’s non-enzymatic, oxidative product—should be taken into consideration when so-called “total THC” amount is concerned. Furthermore, it is noted that the metabolic pathway of CBGA (considered as the progenitor/precursor of all cannabinoids) was activated towards the production of CBCA, which occurred in considerably 3 times higher quantities in the samples cultivated at high altitude.
The presence of similar quantities of CBD in the mountain and plains specimens indicated that geographical location did not significantly influence decarboxylation of CBDA. This process naturally occurs under the action of heat and light, but here it is more probable that it was caused by Kompolti’s predisposition to have prolonged flowering in the late summer/beginning autumn season, when the average daily temperatures are moderately higher in both locations. Therefore, a genetically predisposed flowering season may have led to a partial conversion of the parent CBDA into its neutral counterpart, as was demonstrated for Futura and Finola 75 varieties [13].

2.3. HPLC-Q-Exactive-Orbitrap®-MS Untargeted Metabolomics Approach: Phytocannabinoids Profiling and Identification of Polyphenolic Structures

Metabolomic fingerprinting of the inflorescences that can be used for pharmacological/nutraceutical purposes is important to evaluate a plant metabolite quality and variability. Chromatographic/HRMS fingerprints have been used in the modelling and prediction of pharmacological activities of many medicinal plants [34], but only sporadically for cannabis. Also in this study, as performed in our recently published work [13] the compounds that characterize hemp inflorescences were identified by HPLC-Q-Exactive-Orbitrap®-MS untargeted metabolomics approach that consist of chromatographic separation, HRMS acquisition, and post-analysis data elaboration applying the Compound Discoverer software. Unlike in the above-mentioned paper, here we analyzed the samples separately, using two types of ionization: both positive and negative (Table 3) that enabled a more profound approach. The negative mode revealed the presence of about 80 compounds, whereas the positive mode individuated more than 190. By introducing the detection in negative polarity, it was possible to identify compounds not detected previously [13]. For example, CBGA methyl ester (CBGMA) was not detected in positive ionization, but its presence was clearly confirmed by fragmentation pattern obtained in negative mode (Figure 1). This strategy consequently allowed for an in-depth cluster analysis that clearly demonstrated the differences between the samples coming from the mountains and those cultivated in the plains (Figure 2). Two acquisition modes showed to be complementary in the statistical evaluation of differences between two samples group. As regards the data obtained, it was possible to detect other minor secondary metabolites already identified in other cannabis species. The CBD-family remains to be most abundant (Table 3), enriched by the presence of two sesquiterpene CBDA esters. This CBDA-terpene (inter)reaction needs further elucidation, especially whether it depends on climate/environmental conditions [8,9,13].
In any case, special attention must be paid to cannaflavins that belong to the class of prenylflavonoids. This flavins are secondary metabolites exclusive to the Cannabis genus and were detected in both environments. Their notable presence in mountain samples points towards alternations in their synthesis [16]. The higher levels detected in the mountain samples may possibly be a consequence of the lower average temperature, combined with high solar radiation experienced at the beginning of plant flowering. In fact, it was reported that different classes of flavonoids are involved in plant protection mechanisms, specifically for their radical scavenger activity and screening ability against short wavelength UV-B light [35].
Furthermore, in the Kompolti inflorescence, the remarkable presence of a signal with m/z value at 433.14931 was identified and leads towards the recognition of a particular flavonoid: 3-methoxinobiletin (3,3′,4′,5,6,7,8-heptamethoxyflavone) (Figure 3). This compound has not been reported for any variety of hemp so far. Furthermore, it belongs to the class of polymethoxyflavones that have anti-inflammatory and anti-carcinogenic activities and occur as “novel nutraceutical compounds” [36]. In order to perform its absolute identification, it is necessary to isolate it from the inflorescence and define it with a more detailed analytical approach, also including NMR analysis.
Our metabolomics mapping identified two important phytohormones: salicylic and abscisic acid that have not been reported for cannabinoids inflorescences so far. Their presence was confirmed in our mountain samples, but their amount turned to be more than 10 times higher in lowland specimens. These findings are important to highlight because of the fact that salicylic acid was identified as ‘calorigen’, the plant hormone that induces heat-production in some inflorescences [37]. Also, salicylic acid plays a critical role in the defense against biotrophic pathogens and in the response of plants to abiotic stress, predominantly drought, temperature, heavy metals and, osmotic stress [38]. As far as abscisic acid is concerned, it is well known that stress conditions affect its endogenous production and catabolism rates, while exogenously applied abscisic acid influenced the content and biosynthesis of terpenoids in C. sativa [39].

3. Materials and Methods

3.1. Experimental Fields and Samples Collection

3.1.1. Clones

Seedlings of C. sativa were grown starting from commercial certified seeds of the Kompolti variety. In a 120-plant seedling tray, the most vital plants were chosen for cloning to provide material of genetic homogeneity. Four pairs of clone plants were chosen to grow at altitude versus in the plains, assigning them a number indicating a single plant strain (1,2,3,4) and a letter for the growing site (M for ones destined to be grown at altitude and P for ones destined to be grown in the plains). These strains included only low Δ 9-THC fiber strains. Plants were grown potted in a loam-vermiculite-sand mixture (6:2:1) under ambient greenhouse conditions. Cuttings were taken from the parent pistillate plant of each strain, treated with Rootone, and rooted in perlite. The standard soil permitted to avoid pedotrophic variability. Then, the three-week rooted plants potted in the same substrate were transported in the mountain location (municipality: Alagna Valsesia; elevation 1200 m ASL; latitude 45°51′ N; longitude 7°56′ E) and a lowland location (municipality: Vercelli 130 m ASL; latitude 45°19′ N; longitude 8°22′ E). The geographical area is located in Piedmonte, north of Italy, in the Western Alps ecoregional section for what concerns the mountain location, with prevailing temperate semi-continental bioclimates and in the Po Plain ecoregional section, with prevailing temperate subcontinental bioclimate for what concerns the lowland experimental station [40].

3.1.2. Plant Parts Sampled

Harvest of inflorescences was carried out at flowering, corresponding to the phenological codes 2202 [41]. It was considered “inflorescence” only the 15 cm upper part of the stem. The sectioned parts of the inflorescences were left to air-dry, protected from light in open containers at room temperature (25 °C) for 2 weeks [42]. They were subsequently preserved in plastic bags under vacuum stored in a cool room until analysis. The low temperature avoided changes in metabolites, cannabinoids, and terpenes. The inflorescences material of each clone was sampled five times to realize the analyses. Then, the extracts were injected three times in the analytical instruments.

3.2. Chemical and Reagents

For head-space (HS) analysis, the SPME coating fiber (DVB/CAR/PDMS, 50/30 µm) was obtained from Supelco (Bellefonte, PA, USA) while acetonitrile, 2-propanol, formic acid (all LC-MS grade) were purchased from Carlo Erba (Milan, Italy). Ultrapure water was obtained through a Milli-Q system (Millipore, Merck KGaA, Darmstadt, Germany). All cannabinoids were analytical standards at concentration 1mg/mL (methanolic solution) were purchased from Sigma Aldrich, Round Rock, Texas).

3.3. Superfine Grinding (SFG) Sample Preparation

Samples (1.0 g each) were transformed in fine powder in a high intensity planetary mill at a frequency of 25 Hz for 1 min, using two 50 mL jars (precooled with liquid nitrogen) with 20 mm stainless steel balls.

3.4. Accelerated Solvent Extraction (ASE) for Cannabinoids Profiling

The extraction procedure was done according to the our already-published procedure [12,24,25]. In brief, all extractions were performed by accelerated solvent extraction apparatus using an ASE 350 (Thermo-Fisher Scientific, Waltham, MA, USA) with 34-mL stain steel cells. Inflorescence powder (100 mg) was mixed with an equal weight of diatomaceous earth and transferred into the cell. One-hundred μL of solution containing the IS (diazepam 1 mg/mL) was added and cell was filled with diatomaceous earth. ASE operation parameters were as following: room temperature of 25 °C, pressure (1500 psi), number of static cycles (2 cycles, 5 min each), purging time (60 s with nitrogen) and rinse volume (90%). Extracts (25 mL) obtained using pure methanol and were dried under vacuum; the residue was dissolved in 1 mL of acetonitrile. The resultant solution was diluted (1:10) in starting mobile phase, 2 μL were submitted to analysis by HPLC-Q-Exactive-Orbitrap-MS. Commercially available officinal plants mixture previously analyzed for the absences of cannabinoids served as blank samples and were used to obtain the matrix-matched calibration curves. Matrix-matched calibration curves were obtained by spiking the standard solutions of 14 commercially available cannabinoids that covered the two-concentration range: 0.1 to 10 μg/g and 10–1000 μg g−1.

3.5. Cannabinoids HPLC-Q-Exactive-Orbitrap-MS Evaluation

The cannabinoids profile was assessed employing the method recently published by us [13,24]. HPLC-Q-Exactive-Orbitrap®-MS analysis was achieved on an HPLC Surveyor MS quaternary pump, a Surveyor AS autosampler with a column oven, and a Rheodyne valve with a 20-μL loop system (Thermo Fisher Scientific, San Jose, CA, USA) using a reverse-phase HPLC column 150 × 2 mm i.d., 4 μm, Synergi Hydro RP, with a 4 × 3 mm i.d. C18 guard column (Phenomenex, Torrance, CA, USA). The mobile phase consisted of water and acetonitrile gradient both acidified with 0.1% formic acid. The gradient (flow 0.3 mL/min) started with 95% of 0.1% aqueous formic acid with a linear decrease up to 5% in 30 min. The mobile phase was returned to initial conditions at 35 min, followed by a 5-min re-equilibration period. The column and sample temperatures were 30 °C and 5 °C, respectively. The mass spectrometer Thermo Q-Exactive Plus (Thermo Scientific, San Jose, CA, USA) was equipped with a heated electrospray ionization (HESI) source. Capillary temperature and vaporizer temperature were set at 330 and 380 °C, respectively, while the electrospray voltage was set at 3.30 kV. Sheath and auxiliary gas were 35 and 15 arbitrary units, with S lens RF level of 60. The mass spectrometer was controlled by Xcalibur 3.0 software (Thermo Fisher Scientific, San Jose, CA, USA. Qual Browser in Xcalibur 3.0 software) was used for the exact mass and isotopic pattern determination. The FS-dd-MS2 (full scan data-dependent acquisition) in positive and negative mode was used for both screening and quantification purposes. Resolving power of FS adjusted on 70,000 FWHM at m/z 200, with scan range of m/z 100–900. Automatic gain control (AGC) was set at 3e6, with an injection time of 200 ms. A targeted MS/MS (dd-MS2) analysis operated in both positive and negative mode at 35,000 FWHM (m/z 200). The AGC target was programmed at 2e5, and maximum injection time was set at 100 ms. Fragmentation of precursors was optimized as three-stepped normalized collision energy (NCE) (20, 40 and 40 eV). Detection was based on retention time and on calculated exact mass of the protonated/deprotonated molecular ions, accompanied with fragmentation pattern [13].

3.6. HPLC-Q-Exactive-Orbitrap-MS Untargeted Metabolomics Approach

Raw data from high resolution mass spectrometry were elaborated with Compound Discoverer™ (Thermo Scientific), that facilitated the peak recognition, retention times arrangement, profile assignment, and isotope pattern [43]. Metabolite identification was based on accurate mass and mass fragmentation pattern spectra against MS-MS spectra of compounds available on mzCloud database (HighChem LLC, Slovakia, https://www.mzcloud.org) and in the literature [44]. The ChemSpider Web services (https://www.chemspider.com) and Human Metabolome platform (https://hmdb.ca/) was used as supplementary confirmation tools. If mass fragmentation pattern did not correspond to any of databases annotated by Compound Discoverer™ software, manual confirmation of their fragments using program ChemDrow was completed.

3.7. HS-SPME and GC-MS Analysis for Terpenes Examination

Complete analytical technique was provided in detail in our recently published article [12,24,25]. In brief, 100 mg of inflorescence powder was put into 20 mL glass vials along with 100 μL of the IS (4-metil-2-pentanone, 20 mg/mL in 2-propanol). A cap with a silicon/PTFE septum (Supelco, Bellefonte, PA, USA) was used to close the vial, which was then kept in the temperature block (37 °C), (CTC Analytics, Zwingen, Switzerland). At the end of the sample equilibration time (30 min), a conditioned (60 min at 280 °C) SPME fiber was subjected to the sample for 120 min using a CombiPAL system injector autosampler (CTC Analytics, Zwingen, Switzerland).
Analyses were performed with a Trace GC Ultra coupled to a Trace DSQII quadrupole mass spectrometer (MS) (Thermo-Fisher Scientific, Waltham, MA, USA) equipped with an Rtx-Wax column (30 m × 0.25 mm i.d. × 0.25 µm film thickness) (Restek, Bellefonte, PA, USA). The oven temperature program was: from 35 °C, held for 8 min, to 60 °C at 4 °C/min, then from 60 to 160 °C at 6 °C/min and finally from 160 to 200 at 20 °C/min. Helium was the carrier gas, at a flow rate of 1 mL/min. The MS was operated in electron impact (EI) ionization mode at 70 eV, m/z range of 35–350. An alkanes mixture (C8-C22, Sigma R 8769, Saint Louis, MO, USA) was run under the same chromatographic conditions as the samples to calculate the Kovats Retention Indices (RI) of the detected compounds [27,45]. Compounds were recognized by comparing with authentic standards or by using the Kovats retention indices in combination with the literature and via the National Institute of Standards and Technology (NIST) MS spectral database. The semi-quantitative evaluation was achieved using the internal standard procedure and the results were expressed as µg/g IS equivalents.

3.8. Statistical Analysis

The relative intensity of chromatographic peaks was processed by Compound Discoverer platform that enabled Hierarchical Cluster Analysis. Differences between two groups were evaluated using a two-tailed Student’s t-test from the BioVinci statistical program (Version 1.1.4., BioTuring, Inc. 2018 California, CA, USA). A p-value of less than 0.05 was deemed statistically significant.

4. Conclusions

The quantity and quality of secondary plant metabolites are often attributed to a combination of genetic and environmental factors. Eliminating genetic and pedotrophic factors, the results accomplished in this study indicate qualitatively and quantitatively intraspecific variations in secondary metabolites, other than a major effect attributable to the ecological conditions related to the elevation of the location. A mountain environment, with condition of UV length exposure and critical conditions, deeply influences the quantity of the inflorescence compounds, favoring the production of CBDA and cannaflavins. Information regarding the differences in industrial hemp inflorescences phytochemical profile supports hemp cultivation in mountain areas as a source of pharmacologically active cannabinoids, terpenes and cannaflavones that are considered also as promising nutraceuticals. Metabolomics approaches delineated this crop as resourceful and highly adaptable to the variation of climate/geographical conditions.

Author Contributions

A.G. and V.L. conceptualized the research activities; V.L. completed the GC-MS methodology and formal analysis; R.P. performed LC Orbitrap analysis; L.G. performed statistical analyses and botanical evaluations; V.L., L.G., and R.P., analyzed the data and wrote the article. The first three authors equally contributed to the realization of the research. A.G. was responsible for funding acquisition and team coordination. All authors have read and agreed to the published version of the manuscript.

Funding

The present paper was funded and realized within the “Italian mountain Lab” project and by DARA-CRC Ge.S.Di.Mont. agreement.

Conflicts of Interest

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

References

  1. Regulation (EU) No 1307/2013 of the European Parliament and of the Council of 17 December 2013 Establishing Rules for Direct Payments to Farmers under Support Schemes within the Framework of the Common Agricultural Policy and Repealing Council Regulation (EC) No 637/2008 and Council Regulation (EC) No 73/2009. Available online: https://eur-lex.europa.eu/ (accessed on 27 November 2019).
  2. LEGGE 2 dicembre 2016, n. 242. Disposizioni per la promozione della coltivazione e della filiera agroindustriale della canapa. (16G00258). Available online: https://www.gazzettaufficiale.it/eli/gu/2016/12/30/304/sg/pdf (accessed on 27 November 2019).
  3. Salentijn, E.M.J.; Zhang, Q.; Amaducci, S.; Yang, M.; Trindade, L. New developments in fiber hemp (Cannabis sativa L.) breeding. Ind. Crops Prod 2015, 68, 32–41. [Google Scholar] [CrossRef]
  4. Amaducci, S.; Scordia, D.; Liu, F.H.; Zhang, Q.; Guo, H.; Testa, G.; Cosentino, S.L. Key cultivation techniques for hemp in Europe and China. Ind. Crops Prod 2015, 68, 2–16. [Google Scholar] [CrossRef]
  5. Bonini, S.A.; Premoli, M.; Tambaro, S.; Kumar, A.; Maccarinelli, G.; Memo, M.; Mastinu, A. Cannabis sativa: A comprehensive ethnopharmacological review of a medicinal plant with a long history. J. Ethnopharmacol. 2018, 227, 300–315. [Google Scholar] [CrossRef] [PubMed]
  6. Plant Variety Catalogues, Databases & Information Systems. Available online: https://ec.europa.eu/food/plant/plant_propagation_material/plant_variety_catalogues_databases_en (accessed on 18 January 2019).
  7. de Meijer, E.P.M. Fibre hemp cultivars: A survey of origin, ancestry, availability and brief agronomic characteristics. J. Int. Hemp Assoc. 1995, 2, 66–73. [Google Scholar]
  8. ElSohly, M.; Radwan, M.; Gul, W.; Chandra, S.; Galal, A. Phytochemistry of Cannabis sativa. In Phytocannabinoids: Progress in the Chemistry of Organic Natural Products 103; Kinghorn, A.D., Falk, H., Gibbons, S., Kobayashi, J., Eds.; Springer: Basel, Switzerland, 2017. [Google Scholar]
  9. Hanuš, L.O.; Meyer, S.M.; Muñoz, E.; Taglialatela-Scafati, O.; Appendino, G. Phytocannabinoids: A unified critical inventory. Nat. Prod. Rep. 2016, 33, 1357–1392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Citti, C.; Braghiroli, D.M.; Vandelli, A.; Cannazza, G. Pharmaceutical and biomedical analysis of cannabinoids: A critical review. J. Pharm. Biomed. Anal. 2018, 147, 566–579. [Google Scholar] [CrossRef]
  11. Whiting, P.F.; Wolff, R.F.; Deshpande, S.; Di Nisio, M.; Duffy, S.; Hernandez, A.V.; Keurentjes, J.C.; Lang, S.; Misso, K.; Ryder, S.; et al. Cannabinoids for medical use: A systematic review and meta-analysis. JAMA 2015, 313, 2456–2473. [Google Scholar] [CrossRef]
  12. Pavlovic, R.; Nenna, G.; Calvi, L.; Panseri, S.; Borgonovo, G.; Giupponi, L.; Cannazza, G.; Giorgi, A. Quality traits of “cannabidiol oils”: Cannabinoids content, terpene fingerprint and oxidation stability of european commercially available preparations. Molecules 2018, 23, 1920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Pavlovic, R.; Panseri, S.; Giupponi, L.; Leoni, V.; Citti, C.; Cattaneo, C.; Cavaletto, M.; Giorgi, A. Phytochemical and Ecological Analysis of Two Varieties of Hemp (Cannabis sativa L.) Grown in a Mountain Environment of Italian Alps. Front. Plant Sci. 2019, 10, 1265. [Google Scholar] [CrossRef]
  14. Brenneisen, R. Chemistry and Analysis of Phytocannabinoids and Other Cannabis Constituents. In Marijuana and the Cannabinoids; ElSohly, M.A., Ed.; Humana Press: Clifton, NJ, USA, 2007; pp. 17–49. [Google Scholar]
  15. Pollastro, F.; Minassi, A.; Fresu, L.G. Cannabis Phenolics and their Bioactivities. Curr. Med. Chem. 2018, 25, 1160–1185. [Google Scholar] [CrossRef]
  16. Calzolari, D.; Magagnini, G.; Lucini, L.; Grassi, G.; Appendino, G.B.; Amaducci, S. High added-value compounds from Cannabis threshing residues. Ind. Crop. Prod. 2017, 108, 558–563. [Google Scholar] [CrossRef]
  17. Frassinetti, S.; Moccia, E.; Caltavuturo, L.; Gabriele, M.; Longo, V.; Bellani, L.; Giorgi, G.; Giorgetti, L. Nutraceutical potential of hemp (Cannabis sativa L.) seeds and sprouts. Food Chem. 2018, 262, 56–66. [Google Scholar] [CrossRef] [PubMed]
  18. Russo, E.B. The Case for the Entourage Effect and Conventional Breeding of Clinical Cannabis: No “Strain,” No Gain. Front. Plant Sci. 2019, 9, 1969. [Google Scholar]
  19. Elzinga, S.; Fischedick, R.; Podkolinski, J.; Raber, C. Cannabinoids and terpenes as chemotaxonomic markers in cannabis. Nat. Prod. Chem. Res. 2015, 3, 81. [Google Scholar]
  20. Aizpurua-Olaizola, O.; Soydaner, U.; Öztürk, E.; Schibano, D.; Simsir, Y.; Navarro, P.; Etxebarria, N.; Usobiaga, A. Evolution of the Cannabinoid and Terpene Content during the Growth of Cannabis sativa Plants from Different Chemotypes. J. Nat. Prod. 2016, 79, 324–331. [Google Scholar] [CrossRef] [PubMed]
  21. Lewis, M.A.; Russo, E.B.; Smith, K.M. Pharmacological Foundations of Cannabis Chemovars. Planta Med. 2018, 84, 225–233. [Google Scholar] [CrossRef]
  22. Bertoli, A.; Tozzi, S.; Pistelli, L.; Angelini, L.G. Fibre hemp inflorescences: From crop-residues to essential oil production. Ind. Crop. Prod. 2010, 32, 329–337. [Google Scholar] [CrossRef]
  23. Schluttenhofer, C.; Yuan, L. Challenges towards Revitalizing Hemp: A Multifaceted Crop. Trends Plant Sci. 2017, 11, 917–929. [Google Scholar] [CrossRef] [Green Version]
  24. Calvi, L.; Pentimalli, D.; Panseri, S.; Giupponi, L.; Gelmini, F.; Beretta, G.; Vitali, D.; Bruno, M.; Zilio, E.; Pavlovic, R.; et al. Comprehensive quality evaluation of medical Cannabis sativa L. inflorescence and macerated oils based on HS-SPME coupled to GC–MS and LC-HRMS (q-exactive orbitrap®) approach. J. Pharm. Biomed. Anal. 2018, 150, 208–219. [Google Scholar] [CrossRef]
  25. Calvi, L.; Pavlovic, R.; Panseri, S.; Giupponi, L.; Leoni, V.; Giorgi, A. Quality traits of Medical Cannabis sativa L. inflorescences and derived products based on comprehensive analytical investigation. In Recent Advcance in Cannabinoid Research; IntechOpen: London, UK, 2018. [Google Scholar] [CrossRef] [Green Version]
  26. Giupponi, L.; Pentimalli, D.; Manzo, A.; Panseri, S.; Giorgi, A. Effectiveness of fine root fingerprinting as a tool to identify plants of the Alps: Results of preliminary study. Plant Biosyst. 2017, 152, 464–473. [Google Scholar] [CrossRef]
  27. Giorgi, A.; Manzo, A.; Nanayakkara, N.N.; Giupponi, L.; Cocucci, M.; Panseri, S. Effect of biotic and abiotic stresses on volatile emission of Achillea collina Becker ex Rchb. Nat. Prod. Res. 2015, 29, 1695–1702. [Google Scholar] [CrossRef] [PubMed]
  28. Pavlovic, R.; Borgonovo, G.; Leoni, V.; Giupponi, L.; Ceciliani, G.; Sala, S.; Bassoli, A.; Giorgi, A. Effectiveness of Different Analytical Methods for the Characterization of Propolis: A Case of Study in Northern Italy. Molecules 2020, 25, 504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Namdar, D.; Moran, M.; Ion, A.; Koltai, H. Variation in the compositions of cannabinoid and terpenoids in Cannabis sativa derived from inflorescence position along the stem and extraction methods. Ind. Crops Prod. 2018, 113, 376–382. [Google Scholar] [CrossRef]
  30. Ascrizzi, R.; Ceccarini, L.; Tavarini, S.; Flaminia, G.; Angelini, L.G. Valorisation of hemp inflorescence after seed harvest: Cultivation site and harvest time influence agronomic characteristics and essential oil yield and composition. Ind. Crops Prod. 2019, 139, 111541. [Google Scholar] [CrossRef]
  31. Marchini, L.M.; Charvoz, C.; Dujourdy, L.; Baldovini, N.; Filippi, J.J. Multidimensional analysis of cannabis volatile constituents: Identification of 5,5-dimethyl-1-vinylbicyclo [2 .1.1] hexane as a volatile marker of hashish, the resin of Cannabis Sativa. J. Chromatogr. A 2014, 1370, 200–215. [Google Scholar] [CrossRef]
  32. Glivar, T.; Eržen, J.; Kreft, S.; Zagožen, M.; Čerenak, A.; Čeh, B.; Benković, E.T. Cannabinoid content in industrial hemp (Cannabis sativa L.) varieties grown in Slovenia. Ind. Crops Prod. 2020, 145, 112082. [Google Scholar] [CrossRef]
  33. Zengin, G.; Menghini, L.; Di Sotto, A.; Mancinelli, R.; Sisto, F.; Carradori, S.; Cesa, S.; Fraschetti, C.; Filippi, A.; Angiolella, L.; et al. Chromatographic analyses, in vitro biological activities, and cytotoxicity of Cannabis sativa L. essential oil: A multidisciplinary study. Molecules 2018, 23, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Kharbach, M.; Marmouzic, I.; El Jemlid, M.; Bouklouzeb, A.; Heydena, Y.V. Recent advances in untargeted and targeted approaches applied in herbal-extracts and essential-oils fingerprinting—A review. J. Pharm. Biomed. Anal. 2020, 177, 112849. [Google Scholar] [CrossRef]
  35. Agati, G.; Tattini, M. Multiple functional roles of flavonoids in photoprotection. New Phytol. 2010, 186, 786–793. [Google Scholar] [CrossRef]
  36. Cimaglia, P.; Vieceli Dalla Sega, F.; Vitali, F.; Lodolini, V.; Bernucci, D.; Passarini, G.; Fortini, F.; Marracino, L.; Aquila, G.; Rizzo, P.; et al. Effectiveness of a Novel Nutraceutical Compound Containing Red Yeast Rice, Polymethoxyflavones and Antioxidants in the Modulation of Cholesterol Levels in Subjects With Hypercholesterolemia and Low-Moderate Cardiovascular Risk: The NIRVANA Study. Front. Physiol. 2019, 10, 217. [Google Scholar] [CrossRef]
  37. Raskin, I.; Turner, I.; Melander, W.R. Regulation of Heat Production in the Inflorescences of an Arum Lily by Endogenous Salicylic Acid. PNAS 1989, 86, 2214–2218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Rivas-San Vicente, M.; Plasencia, J. Salicylic acid beyond defence: Its role in plant growth and development. J. Exp. Bot. 1991, 62, 3321–3338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Mansouri, H.; Asrar, Z. Effects of abscisic acid on content and biosynthesis of terpenoids in Cannabis sativa at vegetative stage. Biol. Plant. 2012, 56, 153–156. [Google Scholar] [CrossRef]
  40. Blasi, C.; Capotorti, G.; Copiz, R.; Guida, D.; Mollo, B.; Smiraglia, D.; Zavattero, L. Classification and mapping of the ecoregions of Italy. Plant Biosyst. 2014, 148, 1255–1345. [Google Scholar] [CrossRef]
  41. Mediavilla, V.; Jonquera, M.; Schmid-Slembrouck, I.; Soldati, A. Decimal code for growth stages of hemp (Cannabis sativa L.). J. Int. Hemp Assoc. 1998, 5, 68–74. [Google Scholar]
  42. Hillig, K.W. A chemotaxonomic analysis of terpenoid variation in Cannabis. Biochem. Syst. Ecol. 2004, 32, 875–891. [Google Scholar] [CrossRef]
  43. Li, L.; Li, R.; Zhou, J.; Zuniga, A.; Stanislaus, A.E.; Wu, Y.; Huan, T.; Zheng, J.; Shi, Y.; Wishart, D.S.; et al. MyCompoundID: Using an evidence-based metabolome library for metabolite identification. Anal. Chem. 2013, 85, 3401–3408. [Google Scholar] [CrossRef]
  44. Berman, P.; Futoran, K.; Lewitus, G.M.; Mukha, D.; Benami, M.; Shlomi, T.; Meiri, D. A new ESI-LC/MS approach for comprehensive metabolic profiling of phytocannabinoids in Cannabis. Sci. Rep. 2018, 8, 14280. [Google Scholar] [CrossRef] [Green Version]
  45. Giorgi, A.; Panseri, S.; Mattara, M.S.; Andreis, C.; Chiesa, L.M. Secondary metabolites and antioxidant capacities of waldheimia glabra (decne.) regel from Nepal. J. Sci. Food Agric. 2013, 93, 1026–1034. [Google Scholar] [CrossRef]
Sample Availability: Samples of the compounds are not available from the authors.
Figure 1. Extracted ion chromatogram (m/z = 373.2384) of CBGMA (cannabigerolic acid methyl ester) and corresponding MS/MS spectra identified exclusively in negative mode.
Figure 1. Extracted ion chromatogram (m/z = 373.2384) of CBGMA (cannabigerolic acid methyl ester) and corresponding MS/MS spectra identified exclusively in negative mode.
Molecules 25 01381 g001
Figure 2. Hierarchical cluster analysis: heat-map reflecting the differences between compounds revealed in Kompolti inflorescences in respect to different geographical/climatic conditions.
Figure 2. Hierarchical cluster analysis: heat-map reflecting the differences between compounds revealed in Kompolti inflorescences in respect to different geographical/climatic conditions.
Molecules 25 01381 g002
Figure 3. Tentative identification of 3-methoxynobiletin. Extracted ion chromatogram (m/z = 433.1493) and corresponding MS/MS spectra.
Figure 3. Tentative identification of 3-methoxynobiletin. Extracted ion chromatogram (m/z = 433.1493) and corresponding MS/MS spectra.
Molecules 25 01381 g003
Table 1. Mono/di/triterpenes and sesquiterpenes extracted and identified by headspace solid-phase microextraction-gas chromatographic mass-spectrometry (HS-SPME-GC/MS) in mountain (M) and plains (P) inflorescences.
Table 1. Mono/di/triterpenes and sesquiterpenes extracted and identified by headspace solid-phase microextraction-gas chromatographic mass-spectrometry (HS-SPME-GC/MS) in mountain (M) and plains (P) inflorescences.
MONO/DI/TRI TerpenesM1 P1 p
Value b
M2 P2 p
Value
M3 P3 p
Value
M4 P4 p
Value
mean a± SDmean± SDmean± SDmean± SDmean± SDmean± SDmean± SDmean± SD
α-Pinene3823.981.12705.8152.30.0142401.356.34659.2394.30.0134870.93.71944.865.0<0.0011788.6190.11560.389.1n.s.
α-Fenchenen.d.n.d.26.33.0<0.00144.427.6n.d.n.d.<0.001n.d.n.d.13.11.1<0.001n.d.n.d.12.42.6<0.001
Camphene202.64.7139.15.20.008991.66.8270.210.3<0.001292.55.391.01.1<0.001162.020.1121.60.20.004
β-Pinene2074.0152.82067.3149.80.05012965.655.62979.8233.7<0.0012258.620.71058.594.80.00151456.0336.71230.83.0n.s.
β-Myrcene26,294.5450.323745.11070.40.02076,993.14716.516,893.5519.50.000823,597.9716.75949.924.10.000226,723.62988.911,660.51778.30.004
Limonene2603.1112.53023.2150.3n.s.11472.7276.04184.1134.1<0.0015565.0846.0962.21.40.0115538.993.83883.58.0<0.001
β-Phellandrene642.86.2628.540.6n.s.2017.810.1515.934.60.0012662.0348.3188.52.20.0015681.654.9417.520.40.026
Cis-ocimene352.701.6321.560.3n.s.132.02.238.84.5<0.00144.61.08.60.7<0.00134.49.834.53.2n.s.
γ-Terpinene37.05.923.37.60.00286.51.077.60.40.00469.010.949.91.40.05418.55.315.60.8n.s.
β-Ocimene6571.325.55555.340.10.001135.205.1138.02.0n.s.493.716.328.928.9<0.001139.66.7110.46.30.001
α-Terpinolene136.415.4143.126.3n.s.486.148.4190.241.20.0002274.37.745.40.9<0.001192.06.5144.610.30.002
Terpene6.00.14.80.10.00712.80.95.30.7<0.0014.50.21.20.1<0.00111.40.40.80.0<0.001
α-Fenchone47.52.620.07.00.00984.228.274.624.60.04574.314.39.40.60.01741.10.536.61.60.020
Alloocimene80.90.191.717.1n.s.37.64.124.21.20.01512.14.12.21.2<0.00123.38.317.26.20.036
Linalyl oxide19.60.716.22.10.0708.50.17.25.30.04018.71.217.41.7n.s.21.10.816.61.6n.s.
4,8-Epoxy-p-menth-1-ene147.816.8167.913.00.012n.d.n.d.n.d.n.d.-n.d.n.d.n.d.n.d.-n.d.n.d.n.d.n.d.-
Pinalol117.01.6130.020.0n.s.271.22.1269.611.7n.s.248.231.2229.820.20.1139.812.4131.25.90.1
β-Linalool149.67.7163.049.5n.s.693.63.4576.431.30.028617.049.854.92.4<0.0011118.594.2880.78.60.010
α-Fenchol68.043.029.010.0n.s.164.648.2130.835.9n.s.368.736.612.22.30.0014171.971.612.72.70.060
Verbenol89.836.437.718.90.035204.5155.165.10.6n.s.212.1159.92.40.5n.s.6.86.511.71.3n.s.
tot43,464.5 39,038.8 109,203.3 31,100.5 40,141.9 11,038.4 38,269.1 20,299.2
SesquiterpenesM1 P1 p
Value
M2 P2 p
Value
M3 P3 p
Value
M4 P4 p
Value
mean± SDmean± SDmean± SDmean± SDmean± SDmean± SDmean± SDmean± SD
α-Ylangene130.87.3201190.04447.113.760.76.60.0515.66.53.40.90.00713.36.59.20.2n.s.
α-Copaene41.71.250.74n.s.21.90.148.814.6n.s.9.70.61.90.20.00222.12.320.80.80.027
Zingiberene62.66.4106.810.50.0268.36.841.711.2<0.001260.63512.30.081122.711.7111.84.9n.s.
Longicyclene314.280.9417.559.90.0944.20.72517.2n.s.430.92.50.8<0.001302.663.2207.111.20.08
α-Bergamotene976.958.11440.9115.3n.s.831.8310.9381.849.2n.s.326.554.177.913.90.0043140.1204.72395.886.20.008
Trans-Caryophyllene6487.1113.73668.75020.0089206.9144.93797.741.80.0073345.5456.2829.228.2<0.0017017.4922.61684.638.2<0.001
Aristolene55.83.61059.80.00518.40.72810.3n.s.2.632.630.10.91n.s.n.d.819114.40<0.001
Isoledene33.26.967.310.90.00431.911.931.810n.s.8.44.82.691.8n.s.13.41.112.70.4n.s.
β-Santalene12.32.223.93.70.0056.90.272.5n.s.4.13.150.380.21n.s.8.21.310.72.5n.s.
Aromadendrene15240.276.39.3n.s.195.71.47.60.1<0.00119.91.414.32.8n.s.22.41.42.80.30.001
α-Humulene3190152099183.20.0073107.1177.41284.29.90.0023206.1204.5306.215.80.0143431.487.22660.4100.2<0.001
β-Farnesene51.941.31495.44.6<0.0011105.714.3249.659.60.002n.d.n.d.n.d.n.d.-1494.217.41244.813.40.006
β-Selinene416.429.7632.449.20.003245.547.2502.850.4<0.001123.21.9211.710.60.003112.628.5271.57.60.004
α-Selinene260.811.9379.946.20.02761.86.4145.435.10.0344.90.423.36.10.03374.54.271.28.7n.s.
β-Bisabolene739.630.4293.7156.90.051515.629.5581.221.5<0.00155751.2129.98.80.0031278.28.91024.214.30.003
α-Farnesene213.829.6357.268.40.045661.221.3314.214.40.008160.826.637.921.20.031500.218.2583.286.2n.s.
δ-Cadinene128.220.8215.138.5n.s.63.54.274.24.60.013124.40.223.93.70.00866.85.233.420.2n.s.
β-Maaliene657.225.91048224.20.008297.812.4197.210.20.007863.310.559.614.1n.s. 354.225.13622.7n.s.
Selina-3,7(11)-diene1821134.4n.d.n.d.<0.0011623.3157.2n.d.n.d.<0.001581.681.6132.512.60.0071196.225.2992.3102n.s.
Caryophyllene oxide59.26.363.232n.s.68.26.930.14.30.02418.94.160.70.00667.22.473.64.6n.s.
Guaiol114.520.9144.848.2n.s.326.228.1280.470.7n.s.214.456.633720.2n.s.289.215.4227.98.6n.s.
10-Epi-γ-Eudesmol329.626.7170.747.20.003425.242.4331.23.10.004243.251.442.50.10.005181.121295.23.1n.s.
tot16,248.8 13,057.5 19,974.2 8420.6 9099.2 2277.9 19,708.0 12,409.6
a Data are given as mean ± SD (standard deviation), n = 3 (expressed as µg/g SI equivalents). b p-value—t-test with 95% two-tailed confidence interval for difference of means.
Table 2. Results regarding the bioaccumulation of the main cannabinoids in the inflorescences of mountain Kompolti samples and corresponding plains clones (µg/g, mean of four biological samples ± SD).
Table 2. Results regarding the bioaccumulation of the main cannabinoids in the inflorescences of mountain Kompolti samples and corresponding plains clones (µg/g, mean of four biological samples ± SD).
Mountain PlainsStatistical Significance
MeanSD (±)Mean SD (±)
Neutral forms
CBD5300350060003800ns
Δ9-THC<LOQ/<LOQ//
CBN<LOQ/<LOQ//
CBC460120120500.005
CBG1101018080<0.001
CBDV25040045040<0.001
Δ9-THCV<LOQ/<LOQ//
Acid forms
CBDA99,60024,80068,22015,0000.01
Δ9-THCA8402001010400ns
CBNA4045010ns
CBCA1570200570300.008
CBGA741090045104000.015
CBDVA3107024020ns
Δ9- THCVA <LOQ/<LOQ//
LOQ–limit of quantification 1 µg/g for all phytocannabinoids. ns: not significative
Table 3. Results regarding metabolomic identification in Kompolti inflorescences.
Table 3. Results regarding metabolomic identification in Kompolti inflorescences.
ClassCompoundFormula(M + H)+/Main Fragment(M − H)/Main FragmentRegulationMountain vs. Pains
Phytocannabinoids
CBG cannabigerol classCBGC21H32O2317.2475/193.1223315.2329/191.1078Upregulated in mountain
Sesqui-CBGC26H40O2385.3173/193.1223n.i.
6,7-epoxy-CBGC21H32O3333.2424/315.1867n.i.
CBGVAC20H27O4333.2060/173.0962331.1915/313.1809
6,7-epoxy-CBGAC22H32O5377.2323/341.2113375.2185/257.3077
CBGAC22H31O4361.2375/219.1017359.2228/191.1078
CBGMAC23H34O4n.i.373.2384/355.2293
Sesqui-CBGAC27H40O4n.i.427.2854/409.2748
CBD (cannabidiol) classCBDVC19H26O2287.2006/165.0914285.1860/217.1234Upregulated in mountain
Nor-CBDC20H28O2301.2162/179.1070299.2017/231.1391
CBDC21H30O2315.2319/193.1223313.2173/191.1078
CBDMC22H32O2329.2475/229.0812327.2329/205.1234
CBDVAC20H26O4331.1904/313.1801329.1758/217.1123
Nor-CBDAC20H28O4345.2060/327.1956343.1915/231.1391
CBDAC22H30O4359.2219/341.2114357.2017/245.1547
CBDMAC32H46O4n.i.371.2228/259.1704
Sesquiterpene-CBDA esterC32H46O4495.3469/341.2114493.3323/357.2017
γ-Eudesmyl-CBDA esterC37H54O4562.4017/341.2114561.3949/357.2017
Δ9-THC tetrahydrocannabinol class THCVC19H26O2287.2006/165.0914285.1860/217.1234ns
THCC21H30O2315.2319/193.1223313.2173/n.i.
THCVAC20H26O4331.1904/313.1801329.1758/189.0921
THCAC22H30O4359.2219/341.2114357.2071/245.1547
CBC cannabichromene classCBCVC19H26O2287.2006/165.0914285.1860/163.0765Upregulated in mountain
CBCC21H30O2315.2319/193.1223313.2173/n.d.
CBCVAC20H26O4331.1904/313.1801329.1758/189.0921
CBCAC22H30O4359.2219/341.2114357.2071/313.2179
CBN cannabinol classCBNC21H26O2311.2007/223.1118309.1860/n.i.ns
CBNAC22H26O4355.1904337.1800
Cannaflavin AC26H28O6437.1964/313.0709435.1813/309.0413
Isoprenoid flavonesCannaflavin BC21H20O6369.1333/313.0706367.1195/309.0499Upregulated in mountain
Cannaflavin CC26H28O6437.1964/313.0709435.1813/309.0414
Polymethoxyflavones3-MethoxynobiletinC22H24O9433.14980/403.10296n.i.Upregulated in mountain
FlavonesapigeninC15H10O5271.0601/nd269.0455/117.0348Upregulated in plains
Phenolic acid Salicylic acidC7H6O3n.i.137.0426/95.8554Upregulated in plains
Abscisic acidC15H20O4n.i.263.1289/219.1391
(M + H)+: exact mass of pseudomolecular ion acquired in full scan positive ionization mod; (M − H): exact mass of pseudomolecular ion acquired in full scan negative ionization mode; main fragment: the base fragment in MS/MS spectrum; ns: not significant differences between two chemovars; A: acid; V: C3 chain length; Nor: C4 side chain length; M: methyl ester; n.i.: not identified; Regulation: hierarchical cluster analysis (Figure 3).

Share and Cite

MDPI and ACS Style

Giupponi, L.; Leoni, V.; Pavlovic, R.; Giorgi, A. Influence of Altitude on Phytochemical Composition of Hemp Inflorescence: A Metabolomic Approach. Molecules 2020, 25, 1381. https://doi.org/10.3390/molecules25061381

AMA Style

Giupponi L, Leoni V, Pavlovic R, Giorgi A. Influence of Altitude on Phytochemical Composition of Hemp Inflorescence: A Metabolomic Approach. Molecules. 2020; 25(6):1381. https://doi.org/10.3390/molecules25061381

Chicago/Turabian Style

Giupponi, Luca, Valeria Leoni, Radmila Pavlovic, and Annamaria Giorgi. 2020. "Influence of Altitude on Phytochemical Composition of Hemp Inflorescence: A Metabolomic Approach" Molecules 25, no. 6: 1381. https://doi.org/10.3390/molecules25061381

APA Style

Giupponi, L., Leoni, V., Pavlovic, R., & Giorgi, A. (2020). Influence of Altitude on Phytochemical Composition of Hemp Inflorescence: A Metabolomic Approach. Molecules, 25(6), 1381. https://doi.org/10.3390/molecules25061381

Article Metrics

Back to TopTop