Mining the chemical diversity of the hemp seed (Cannabis sativa L.) metabolome: discovery of a new molecular family widely distributed across hemp
1 Enhanced Partnerships Department Royal Botanic Gardens, Kew, London, United Kingdom
2 Department of Botany and Plant Biotechnology, University of Johannesburg (Auckland Park Campus), Auckland Park, Johannesburg, South Africa
*Correspondence: Guillermo Federico Padilla-González, federico.padillagonzalez@wur.nlAbstract
Hemp (Cannabis sativa L.) is a widely researched industrial crop with a variety of applications in the pharmaceutical, nutraceutical, food, cosmetic, textile, and materials industries. Although many of these applications are related to its chemical composition, the chemical diversity of the hemp metabolome has not been explored in detail and new metabolites with unknown properties are likely to be discovered. In the current study, we explored the chemical diversity of the hemp seed metabolome through an untargeted metabolomic study of 52 germplasm accessions to 1) identify new metabolites and 2) link the presence of biologically important molecules to specific accessions on which to focus on in future studies. Multivariate analysis of mass spectral data demonstrated large variability of the polar chemistry profile between accessions. Five main groups were annotated based on their similar metabolic fingerprints. The investigation also led to the discovery of a new compound and four structural analogues, belonging to a previously unknown chemical class in hemp seeds: cinnamic acid glycosyl sulphates. Although variability in the fatty acid profiles was not as marked as the polar components, some accessions had a higher yield of fatty acids, and variation in the ratio of linoleic acid to α-linolenic acid was also observed, with some varieties closer to 3:1 (reported as optimal for human nutrition). We found that that cinnamic acid amides and lignanamides, the main chemical classes of bioactive metabolites in hemp seed, were more concentrated in the Spanish accession Kongo Hanf (CAN58) and the French accession CAN37, while the Italian cultivar Eletta Campana (CAN48) demonstrated the greatest yield of fatty acids. Our results indicate that the high variability of bioactive and novel metabolites across the studied hemp seed accessions may influence claims associated with their commercialization and inform breeding programs in cultivar development.
1Introduction
Cannabis sativa L., specifically the subspecies, varieties and cultivars known as hemp, are considered one of the most versatile herbaceous plants in agriculture, because of its commercial value to the pharmaceutical, food, nutraceutical, cosmetics, textile, papermaking, and construction industries (Schluttenhofer and Yuan, 2017; Farinon et al., 2020). In the context of its nutraceutical use, the focus has been on its seeds, which are rich in polyunsaturated fatty acids, proteins, minerals, and specialized metabolites with noteworthy implications for human and animal nutrition (Callaway, 2004; Russo and Reggiani, 2015; Galasso et al., 2016; Farinon et al., 2020). In this regard, hemp seeds are processed to make edible oil, cake flour and protein powder (Farinon et al., 2020; Leonard et al., 2021b).
Considering that some varieties of C. sativa express psychoactive cannabinoids that are occasionally detected in Cannabis-based products, a distinction between hemp and the recreational variety has been made (Schluttenhofer and Yuan, 2017). The European Industrial Hemp Association (EIHA) defines hemp as “the plant Cannabis sativa L., or any part of it, with a delta-9 tetrahydrocannabinol (THC) concentration up to 0.3% on a dry weight basis” to differentiate it from psychoactive C. sativa in which the THC concentration exceeds 0.3% and can be up to 20% (European Industrial Hemp Association, https://eiha.org/). Subsequently, several countries mainly from Europe, North America and Asia legalized the cultivation of low THC hemp cultivars (Leonard et al., 2021b). By 2018 more than 50,000 hectares were allocated to the cultivation of this crop in Europe, one of the main producers of industrial hemp. Within Europe, France is currently the largest producer, followed by Italy and the Netherlands (European Industrial Hemp Association, https://eiha.org/). Over the years, varieties of hemp were developed to select for specific traits and seeds from many of these varieties are currently stored in international germplasm collections and botanical gardens. There are currently 75 varieties registered in the EU Catalogue, all of them with a THC content below 0.3% and a diverse profile of other metabolites (European Industrial Hemp Association, https://eiha.org/).
Beyond the well-known cannabinoids, C. sativa contains a diverse metabolome of bioactive metabolites implicated in human health and nutrition. Hemp seeds are low in cannabinoids, and rich in phenylpropionamides (cinnamic acid amides and lignanamides) and unsaturated fatty acids, attracting nutraceutical and commercial interest (Crescente et al., 2018; Leonard et al., 2021a; Leonard et al., 2021b). The principal nutraceutical value of hemp seed oil is in its fatty acid composition, dominated by >90% of polyunsaturated fatty acids (Schluttenhofer and Yuan, 2017). It contains two dietary essential fatty acids, linoleic acid and α-linolenic acid in the ratio of 2.5–3:1. This is allegedly ideal for human nutrition and cardiovascular health (Leizer et al., 2000; Simopoulos, 2008; Galasso et al., 2016). However, the presence of high quantities of phenylpropionamides in hemp seeds have also been linked to some of their biological properties (Chen et al., 2012; Moccia et al., 2020; Leonard et al., 2021a). Hemp seeds accumulate the highest structural diversity of lignanamides among other lignanamides-producing taxa, with more than 80 different compounds (Leonard et al., 2021a).
Cinnamic acid amides and lignanamides display potent anti-inflammatory, antioxidant and anti-cancer activities in both in-vitro and in-vivo studies (Chen et al., 2012; Farinon et al., 2020; Moccia et al., 2020; Leonard et al., 2021a). Previous studies have promoted these hemp polyphenols as protective agents against human chronic diseases (Leonard et al., 2021a). For example, phenylpropionamides such as cannabisin B and N-trans-caffeoyltyramine are significantly stronger antioxidants than the standard soy isoflavones (Chen et al., 2012). Subsequent studies demonstrated that cannabisin B also enacts antiproliferative activity by inducing autophagic cell death in liver hepatocarcinoma (HepG2) cells (Chen et al., 2013). Other hemp lignanamides, such as cannabisin F, have demonstrated a potential neuro-protective effect by reducing mRNA levels of pro-inflammatory mediators, intracellular reactive oxygen species (ROS) and tumor necrosis factor α (TNF-α) in lipopolysaccharide-stimulated BV2 microglia cells (Wang et al., 2019).
Thus, the health and nutritional properties of hemp seeds are dependent upon their chemical composition that differs mainly according to the variety used. With such a high number of hemp varieties currently available in international germplasm collections, a comprehensive chemical assessment of these accessions is especially relevant, specifically to discover new metabolites and link the presence of biologically important molecules to certain accessions, ultimately to exploit the hemp germplasm by designing commercial cultivars with specific attributes. However, to date the phytochemical diversity of the hemp seed metabolome has not been explored in detail. In the present study, we focused on characterizing the chemical diversity of the hemp seed metabolome across cultivars/accessions using an untargeted metabolomics study of 52 hemp seed accessions.
2Materials and methods
2.1Plant material
A collection of 51 hemp seed accessions from fourteen countries were obtained from the Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), Gatersleben, Germany (http://gbis.ipk-gatersleben.de/) ( Table 1 ). These accessions, maintained ex-situ by the IPK, were provided as small quantities of seeds (< 1 g). An additional sample of hemp seeds (1 kg) (S1, Table 1 ) was acquired in a local supermarket for further isolation of targeted metabolites (see section 2.5), bringing the total number of samples to 52.
| Code | Scientific name | Cultivar or Accession name | Country of origin | EU catalogue status |
|---|---|---|---|---|
| CAN16 | Cannabis sativa L. | – | Slovakia | – |
| CAN17 | Cannabis sativa L. | – | Hungary | – |
| CAN18 | Cannabis sativa L. subsp. sativa | – | Germany | – |
| CAN19 | Cannabis sativa L. | – | Italy | – |
| CAN20 | Cannabis sativa L. | – | North Korea | – |
| CAN21 | Cannabis sativa L. subsp. sativa | – | Romania | – |
| CAN22 | Cannabis sativa L. subsp. spontanea Serebr. | – | Georgia | – |
| CAN23 | Cannabis sativa L. subsp. sativa | – | North Korea | – |
| CAN24 | Cannabis sativa L. | – | Italy | – |
| CAN26 | Cannabis sativa L. | – | Turkey | – |
| CAN27 | Cannabis sativa L. | – | unknown | – |
| CAN28 | Cannabis sativa L. | Forose | unknown | – |
| CAN29 | Cannabis sativa L. subsp. sativa | – | Romania | – |
| CAN30 | Cannabis sativa L. | Bernburger | Germany | – |
| CAN31 | Cannabis sativa L. subsp. sativa | Cunepa | Romania | – |
| CAN32 | Cannabis sativa L. subsp. sativa | Conopla (ukr.) | Romania | – |
| CAN33 | Cannabis sativa L. subs. sativa | Cinepa | Romania | – |
| CAN34 | Cannabis sativa L. subs. sativa | Cinepa | Romania | – |
| CAN35 | Cannabis sativa L. subs. sativa | Cinepa | Romania | – |
| CAN36 | Cannabis sativa L. subsp. sativa | – | unknown | – |
| CAN37 | Cannabis sativa L. subsp. sativa | – | France | – |
| CAN38 | Cannabis sativa L. subsp. sativa | Cinepa | Romania | – |
| CAN39 | Cannabis sativa L. | – | China | – |
| CAN40 | Cannabis sativa L. | – | Italy | – |
| CAN42 | Cannabis sativa L. subsp. sativa | – | Croatia | – |
| CAN43 | Cannabis sativa L. subsp. sativa | Hohenthurmer Gleichzeitig Reifender | Germany | – |
| CAN44 | Cannabis sativa L. subsp. sativa | Fibrimon | Germany | – |
| CAN45 | Cannabis sativa L. subsp. sativa | Krasnodarskaya | Russia* | – |
| CAN46 | Cannabis sativa L. | Fibrimon 56 | France* | Deleted |
| CAN47 | Cannabis sativa L. subsp. sativa | – | Turkey | – |
| CAN48 | Cannabis sativa L. | Eletta campana | Italy | Registered |
| CAN50 | Cannabis sativa L. subsp. sativa | Fibrimon | France* | – |
| CAN51 | Cannabis sativa L. subsp. sativa | – | Argentina | – |
| CAN52 | Cannabis sativa L. subsp. sativa | Fibridia | unknown | – |
| CAN53 | Cannabis sativa L. subsp. sativa | Fibrimon | France* | – |
| CAN54 | Cannabis sativa L. subsp. sativa | Fibrimon 21 | France* | Deleted |
| CAN55 | Cannabis sativa L. subsp. sativa | Havellander | Russia* | – |
| CAN56 | Cannabis sativa L. subsp. sativa | Kompolti | Hungary* | Registered |
| CAN58 | Cannabis sativa L. subsp. sativa | Kongo Hanf | Spain | – |
| CAN59 | Cannabis sativa L. subsp. sativa | Turken P 494 | unknown | – |
| CAN60 | Cannabis sativa L. subsp. sativa | Schurig | Russia* | – |
| CAN61 | Cannabis sativa L. subsp. sativa | Carmagnola in Selezione | Italy | Registered |
| CAN62 | Cannabis sativa L. subsp. sativa | Bredemann P | Germany* | – |
| CAN63 | Cannabis sativa L. subsp. sativa | Ramo | Germany | – |
| CAN64 | Cannabis sativa L. subsp. sativa | – | Turkey | – |
| CAN65 | Cannabis sativa L. subsp. sativa | Juso 14 | Germany | Deleted |
| CAN66 | Cannabis sativa L. subsp. sativa | Lovrin 110 | Germany | Registered |
| CAN67 | Cannabis sativa L. subsp. sativa | Fasamo | Germany | Deleted |
| CAN68 | Cannabis sativa L. subsp. sativa | Kompolti | Germany | Registered |
| CAN69 | Cannabis sativa L. subsp. sativa | Futura | Germany | Deleted |
| CAN70 | Cannabis sativa L. subsp. sativa | Kompolti | Germany | Registered |
| S1 | Cannabis sativa L. | – | unknown | – |
The 51 externally sourced samples derive from sixteen recognised cultivars, including five historic cultivars whose seeds have a “deleted” status ( Table 1 ) in the EU database of registered plant varieties (https://ec.europa.eu/info/index_en). Samples that were attained in replicate include: three accessions of the cultivar “Fibrimon” (plus two related but deleted strains, Fibrimon 21 and 56), four of “Cinepa” and three of “Kompolti” ( Table 1 ).
2.2Extraction and LC-MS analyses
The extraction and metabolic profiling by LC-MS was based on the protocol reported by De Vos et al. (2007). Seeds from each of the hemp accessions were ground in liquid nitrogen with a pestle and mortar. Fifty milligrams of the resultant seed powder were then transferred to Eppendorf tubes and 2 mL of 80% MeOH (HPLC grade) were added as extraction solvent. Tubes were then placed in an ultrasonic bath for 15 min at 25°C using a frequency of 40 kHz. After extraction, samples were centrifuged at 19,975 × g for 10 min at room temperature, and the supernatant was transferred to an LC-MS vial. Samples were randomly analyzed by arbitrarily selecting vials to minimize statistical bias. Furthermore, 30 µL of each of the 52 samples were pooled into a LC-MS vial labelled as “Quality Control”. This sample was analyzed multiple times along the chromatographic sequence to check for reproducibility (De Vos et al., 2007).
The metabolic profiles of the hemp seed accessions were recorded on a Vanquish UHPLC system coupled to a 100 Hz photodiode array detector (PDA) and an Orbitrap Fusion Tri-hybrid high-resolution tandem mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). Chromatographic separation of hemp seed extracts (5 μL) was performed on a Luna C18 column (150 mm × 3 mm i.d., 3 μm, Phenomenex, Torrance, CA, USA) using a linear mobile phase gradient of 0:90:10 to 90:0:10 [MeOH (A): water (C): acetonitrile +1% formic acid (D)] over 20 min at a flow rate of 400 μL min−1. Ultraviolet data were recorded between 210 nm and 550 nm.
Mass spectrometric detection was performed in the positive and negative ionization modes using the full scan and data dependent MS2 and MS3 acquisition modes. Total ion current (TIC) chromatograms were obtained over the range of 125 – 1800 m/z using a spray voltage of +3.5 and −2.5 kV for the positive and negative ionization modes, respectively. Four different scan events were recorded for each ionization mode as follows: (1) full scan, (2) MS2 of the most intense ion in full scan, (3) MS3 of the most intense ion in MS2, and (4) MS3 of the second most intense ion in MS2. Additional parameters for the MS included the following: full scan resolution set to 60,000 (full-width at half-maximum, FWHM), capillary temperature set to 350°C, ion transfer tube temperature set to 325°C, RF lens set to 50%, automatic gain control target set to 4.0 × 105 (full scan) or 1.0 × 104 (MS2 and MS3), intensity threshold set to 1.0 × 104, collision-induced dissociation energy set to 35 eV, activation Q set to 0.25, and isolation window set to 4 m/z. Nitrogen was used as the drying, nebulizer and fragmentation gas.
2.3LC-MS data processing and multivariate analysis
The LC-MS raw data was sliced into two sets according to the ionization mode (positive and negative) and transformed to mzXML format using the MSConvert package from the software ProteoWizard 3.0.9798 (Proteowizard Software Foundation, Palo Alto, CA, USA). Each data set was then processed by MZmine 2.53 (Pluskal et al., 2010) following the protocol previously described (Padilla-González et al., 2020a; Padilla-González et al., 2020b). After MZmine pre-processing, the results were exported as.csv and.mgf files which were then uploaded to the GNPS platform (https://gnps.ucsd.edu) for Feature Based Molecular Networking (FBMN) analysis (Wang et al., 2016; Nothias et al., 2020). The.csv file containing quantitative information related to ion abundances in each hemp seed sample was also submitted to multivariate statistical analysis by principal component analysis (PCA), hierarchical clustering analysis (HCA) and heatmaps in the software R 3.0.3 (R Foundation for Statistical Computing, Austria). Prior to the multivariate analyses the data matrix was centered and scaled by unit variance to avoid statistical bias influenced by differences between high and low abundant metabolites. LC-MS files, R scripts and detailed parameters used in MZmine (available as an MZmine batch file) are freely available in the MassIVE repository (study identifier MSV000090725).
2.4Feature-based molecular networking and annotation tools
Feature-based molecular networking (FBMN) was performed following the workflow by (Nothias et al., 2020) on the GNPS platform (https://gnps.ucsd.edu). For this analysis, the data were first filtered by removing all MS/MS fragment ions within a window of +/−17 Da of the precursor ion value and by choosing only the top 6 fragment ions in the +/− 50 Da window throughout the spectrum. The precursor and fragment ion mass tolerance were set to 0.05 Da and 0.99 Da respectively, in accordance with the MS settings used. A cosine score above 0.65 and more than six matched peaks were considered for creating the molecular network. Furthermore, edges between two nodes were kept in the network only if each of the nodes appeared in each other’s respective top 10 most similar nodes. Experimental MS2 spectra were then searched against the GNPS spectral libraries (Wang et al., 2016) using the same filters as for the input data. To enhance chemical structural information within the molecular network and provide a more comprehensive overview of the hemp seed chemical space, information from in silico structure annotations from the GNPS library search were incorporated into the network using the GNPS MolNetEnhancer and MS2LDA workflows (Wandy et al., 2018; Ernst et al., 2019). MolNetEnhancer and MS2LDA were performed using the default values reported in GNPS, except for bin width in MS2LDA, which was set at 0.1. Chemical class annotations were performed using the ClassyFire chemical ontology (Feunang et al., 2016). Final molecular networks including quantitative data, substructure information and chemical classifications were then visualized using the software Cytoscape (Ono et al., 2014).
Lastly, to confirm and expand the spectral library annotation made by molecular networking and its in-silico match against the GNPS spectral libraries, accurate mass values, MS2 and MS3 data of all the annotated metabolites were manually inspected and compared with information available in the literature to achieve different confidence values following the Metabolomics Standards Initiative (Sumner et al., 2007).
2.6Extraction and derivatization of fatty acids
Powdered seeds (100 mg) from each hemp accession were extracted with 2 mL of n-hexane in an ultrasonic bath for 15 min at 25°C using a frequency of 40 kHz. After extraction, samples were evaporated under vacuum to obtain an oily residue corresponding to hemp seed fixed oil. Five microliters of the oil were added to a 2 mL glass vial containing 500 µL of boron trifluoride in methanol (20%). The vial was filled with nitrogen, tightly sealed and gently mixed. The sample vial was subsequently heated at 100 °C for 15 min in a heating block. After cooling to room temperature, 500 µL of water and 750 µL of pentane were added to the vial and vigorously shaken. The organic phase (upper) was then transferred to another vial and evaporated under nitrogen. This partition was repeated for a second time and the organic phases were mixed and evaporated under nitrogen. The concentrated sample containing fatty acids methyl esters was re-dissolved in 200 µL of n-hexane and placed in a glass vial with insert for GC-MS analysis.
2.7GC-MS analyses of fatty acids methyl esters
GC-MS analyses of fatty acids methyl esters were performed in an Agilent 7890A GC chromatograph coupled to a single quadrupole MS analyser (Agilent 5975C). Chromatographic separation was performed on a ZB-WAX column (30 m x 0.25 mm i.d. x 0.25 µm, Phenomenex), using helium as carried gas at a flow rate of 1 mL per minute. The temperature program was configured from 70 – 250 °C at a linear increment of 3 °C per minute. The injection temperature was set to 220 °C and the source temperature to 180 °C. 70 eV was used as ionization energy and the mass spectrometer was configured to record data from 38 to 650 m/z. Compounds were identified by comparing their retention indices (calculated against a series of n-alkanes analysed under the same experimental conditions) and by comparing their mass spectra with the NIST library.
The percentage composition of the fatty acid methyl esters detected in each hemp cultivar/accession was semi-quantitatively obtained by manual integration of GC-MS data. This information was then centered and scaled by unit variance and submitted to statistical analyses in the software R.
3Results
3.2Structural characterization of new cinnamic acid glycosyl sulphates
Considering that the analysis by FBMN and our dereplication approach suggested the presence of several potentially new metabolites belonging to a previously unreported molecular family in hemp (cinnamic acid glycosyl sulphates), we performed a target isolation of the most abundant metabolite from this chemical class, detected at 6.84 min with a mass of 435.06061 m/z [M−H]− ( Table 2 ). The process started with a dried 20% methanolic extract of hemp seeds, which was submitted to the classic isolation processes using Sephadex LH-20 column chromatography followed by semipreparative HPLC (see Materials and Methods section). Approximately 4 mg of ferulic acid 4-O-glucosyl-6’-O-sulphate (compound 5) was afforded. Furthermore, the tentative structures of four additional cinnamic acid glycosyl sulphates (compounds 6 – 9; Figure 4 ) in the crude extracts are tentatively reported based on similarity to 5 according to HRMS and MS2 data. All the compounds in this chemical class (compounds 5 – 9, Figure 4 ) showed a characteristic base peak ion in their MS2 spectra at 241 m/z [M−H]−, representing a deprotonated glucopyranosyl-sulphate unit formed after the neutral loss of the phenylpropanoid moiety. A subsequent neutral loss of a water molecule generated the ion at 223 m/z ( Figure 4 ). In addition to these diagnostic ions, minor peaks representing a deprotonated ferulic acid (193 m/z), caffeic acid (179 m/z), cinnamic acid (163 m/z) and dimethoxycinnamic acid unit (223 m/z), were observed in the MS2 spectra of compounds 5, 6, 7 and 8, respectively. The suggested fragmentation mechanism leading to the formation of these diagnostic ions is provided in Figure 4 .
Compounds 5 and 9 were assigned as two isomers of the previously undescribed ferulic acid glucosyl sulphate metabolite, based on chromatographic and spectrometric data (HRMS and MS2). Both compounds showed a deprotonated molecule [M−H]− at m/z 435.0606 ( Table 2 ), consistent with the molecular formula of C16H20SO12 (calculated for C16H19SO12, 435.0603). Their online UV spectra displayed absorbance maxima at 289 and 315 nm, characteristic of phenylpropanoid derivatives.
The 1H NMR spectrum of compound 5 included a proton spin system consistent with a six membered meta-para-disubstituted aromatic ring ( Table 3 ), and two adjacent exocyclic olefinic monoprotonated carbons with a coupling constant (15.8 Hz) consistent with the trans configuration. The 13C spectrum is consistent with that of a substituted free ferulic acid moiety. Positions 3 and 4 of the ferulic acid moiety were assigned by HMBC couplings with other protons on the ring, and by an online NMR shifts simulator (nmrdb.org). The position of the methoxy singlet was established by a HMBC coupling to the carbon at position 3, and an NOE coupling to the proton doublet (1.2 Hz) at position 2. The attachment of the glucosyl moiety was established by a HMBC coupling from the anomeric proton to the fully substituted carbon at position 4.
| Feruloyl | Glucose | ||||||
|---|---|---|---|---|---|---|---|
| No. | 13C | 1H | HMBC, H→C | No. | 13C | 1H | |
| 1 | 130.73 | – | – | 1’* | 102.34 | 4.94, d (3.5 Hz) | |
| 2 | 112.37 | 7.22, d (1.2 Hz) | 1, 3 | 2’ | 74.8 | 3.52, dd (3.5, 8.4 Hz) | |
| 3 | 150.98 | – | – | 3’ | 71.27 | 3.42, dd (7.6, 8.4 Hz) | |
| OCH3 | 56.8 | 3.82, 3H s | 3 | 4’ | 77.58 | 3.48, dd (7.5, 7.6 Hz) | |
| 4 | 149.8 | – | – | 5’ | 71.27 | 3.68, dd (5.9, 7.5 Hz) | |
| 5 | 116.3 | 6.69, d (7.9 Hz) | 6 | 6’ | 68.15 | 4.13, dd (5.9, 11.1 Hz), | |
| 6 | 130.95 | 7.03, dd (1.2, 7.9 Hz) | 1 | 4.36, d (11.1 Hz) | |||
| 7 | 145.13 | 7.56, d (15.8 Hz) | 2, 8, 9 | ||||
| 8 | 119.3 | 6.40, d (15.8 Hz) | 1 | ||||
| 9 | 171.2 | – | – | ||||
The glucosyl moiety was assigned using COESY couplings, coupling constants and HSQC. The position of attachment of the sulphate was confirmed by the re-shielding of the carbon at position 6’, which is normally measured at 61 – 62 ppm in similar glycosyl-arene structures (Schuster et al., 1986).
4Discussion
Hemp is an industrial crop that has been widely investigated due to its diversity of applications in the pharmaceutical, nutraceutical, food, cosmetic, textile, and material industries (Schluttenhofer and Yuan, 2017; Farinon et al., 2020). However, the chemical diversity of the hemp seed metabolome has not been explored in detail and new metabolites with unknown properties are likely to be discovered. In the current study, following an untargeted metabolomics approach by LC-MS and GC-MS, we discovered a previously unreported molecular family in hemp seeds “cinnamic acid glycosyl sulphates” and linked the presence of biologically important molecules to certain accessions. Chemical knowledge of these cultivars and accessions could be used as a base in future to optimize commercial cultivars in the context of nutritional and other functional uses.
Metabolic fingerprinting by LC-MS and data mining tools revealed that the 52 hemp seed accessions can be categorized into five main groups based on the similarities of their metabolic fingerprints. Although it is expected that when growing the same seeds under homogeneous conditions there will be some differences in the metabolic grouping, we believe that we can only have an accurate representation of the true metabolic phenotype of an organism by studying the chemical composition of plants/organs obtained from individuals collected in their natural environment. Germplasm and herbarium collections provide valuable samples for this. In the current study, among all the studied accessions and cultivars, the Spanish accession Kongo Hanf (CAN58) and the French accession CAN37 showed a remarkable structural diversity of bioactive and novel metabolites, some of them accumulated in high amounts. These two varieties showed the highest accumulation of phenylpropionamides (including cinnamic acid amides and their oxidative coupling products lignanamides), as well as new cinnamic acid glycosyl sulphates.
Phenylpropionamides (cinnamic acid amides and lignanamides) represents the main chemical class of bioactive metabolites in hemp seeds (Moccia et al., 2020; Leonard et al., 2021a; Leonard et al., 2021b). These compounds display potent antioxidant, anti-inflammatory, cytotoxic and acetylcholinesterase inhibitory activities, both in vitro and in vivo studies (Yan et al., 2015; Farinon et al., 2020; Moccia et al., 2020; Leonard et al., 2021a). For instance, previous studies have recognized the strong potential of hemp seed phenylpropionamides as protective agents against human chronic diseases (Leonard et al., 2021a), although results from clinical trials or cohort studies are not available to support this. Our results indicate that N-trans-caffeoyltyramine (compound 11) and N-trans-feruloyltyramine (compound 12) constitute the main phenylpropionamides accumulated in hemp seeds, especially in the Spanish accession Kongo Hanf (CAN58) and the French accession CAN37. These two compounds have demonstrated strong antioxidant activities (Chen et al., 2012; Moccia et al., 2020) and in vivo anti-neuroinflammatory action (Maiolo et al., 2018; Zhou et al., 2018). Furthermore, the lignanamides fraction of hemp seeds are also highly active. Like compounds 11 and 12, lignanamides (compounds 16 – 38) are also accumulated in high amounts in the accession CAN58 and CAN37.
Recent studies by Chen et al. (2012) observed stronger DPPH scavenging activity by cannabisin B (compound 17) isolated from hemp seed than common soy (Glycine max (L.) Merr.) isoflavones and secoisolariciresinol diglucoside from flax seed (Linum usitatissimum L.). Cannabisins M, A, C, D, and 3,3’-demethylgrossamide (compounds 21, 23, 27, 32 and 20, respectively) from hemp seed also displayed potent antioxidant activities in vitro (Yan et al., 2015). These compounds are especially accumulated in accessions CAN37, CAN58, CAN68 and CAN42.
The anti-inflammatory properties that these and other hemp seed lignanamides display have also been widely recognized. For example, grossamide, (compound 35) accumulated in higher amounts in accession CAN37 ( Figure 3 ), showed a significant reduction in the production of the proinflammatory cytokines interleukin-6 (IL-6) and tumor necrosis factor-a (TNF-a) at a concentration of 10 – 20 µM in a lipopolysaccharide (LPS) induced model (Luo et al., 2017). These cytokines have been linked to the stimulation of inflammatory and auto-immune processes in many diseases such as rheumatoid arthritis, multiple sclerosis, psoriasis as well as different types of cancer (Bradley, 2008; Blauvelt, 2017; Shang et al., 2017). Like grossamide, cannabisin F (compound 34) has also been reported to reduce the production of IL-6 and TNF-a in LPS-induced models, although through a different mechanism (Wang et al., 2019).
Regarding the cytotoxic effects of hemp seed lignanamides, previous studies have reported that it varies in accordance with their chemical structures. For example, cannabisin B has been shown to induce autophagic cell death on HepG2, a human liver cancer cell line (Chen et al., 2013), while cannabisin F is especially active against SH-SY5Y, a human neuroblastoma cell line (Li et al., 2019). Nevertheless, other hemp seed lignanamides, such as grossamide, cannabisin D and cannabisin G, have displayed weak cytotoxic activity against human prostate cancer LNCaP cells (Ma et al., 2002). Considering that CAN37, CAN58, CAN68 and CAN42 have a higher concentration of bioactive cinnamic acid amides and lignanamides, these accessions could be commercially prioritized as potential sources of antioxidant, anti-inflammatory and cytotoxic metabolites or in selective breeding programs to improve the concentration of specific metabolites in other commercialized cultivars. Interestingly, none of the five historic cultivars with deleted status in the EU database ( Table 1 ) showed a particularly high concentration of bioactive metabolites compared to the newer accessions.
In addition to the 29 phenylpropionamide peaks annotated in the current study ( Table 2 ), our results suggest that other lignanamides might also be present. Manual inspection of the FBMN results indicate that there are at least eleven additional nodes representing potentially new phenylpropionamides, some of them showing m/z values previously unreported in hemp seeds. For example, two nodes at m/z 600.270 and m/z 614.286 [M+H]+ appear to be structurally related to tricoumaroyl spermidine (compound 14). Manual interpretation of MS data suggests that the node at m/z 600.270 [M+H]+ likely represents a hydroxylated analogue, while the node at m/z 614.286 [M+H]+ is likely its methoxylated version. Several isomers of m/z 597.33 and m/z 611.23 [M+H]+ were also observed. The isolation and complete structural characterization of these potentially new metabolites is still necessary. Recent studies have proved that despite being exhaustively commercialized, hemp seed continues to be a promising source of new lignanamides. For example, Yan et al. (2015) identified the first trimeric feruloyltyramine in hempseeds, of which its coupling resulted in the formation of an 8-O-4 linkage. Considering that hemp hulls are the main site of accumulation of lignanamides, future agricultural waste valorisation programs are highly encouraged. Hemp hulls are a by-product of the hemp seed oil industry (Pojić et al., 2014; Mattila et al., 2018), which represents an opportunity. This is further supported by the fact that the total lignanamides in hemp hull can be up to four times higher than in whole hemp seed (3.36 mg/g vs 0.77 mg/g, respectively) (Mattila et al., 2018).
Aside from phenylpropionamides, hemp seed oil is also a commodity held in high regard for its nutraceutical properties due to a perfectly balanced content of omega 3 and omega 6 polyunsaturated fatty acids (Crescente et al., 2018; Citti et al., 2019). Mean concentration values for the six main fatty acids detected in the current study: palmitic acid oleic acid, linoleic acid, γ-linolenic acid, α-linolenic acid and stearidonic acid ( Table 4 ), are within ranges comparable to previous studies (Galasso et al., 2016; Pavlovic et al., 2019). The Italian cultivar Eletta Campana (CAN48) showed the highest accumulation of fatty acids among the investigated hemp seed accessions. However, CAN32, CAN50 and CAN53 were the accessions with the ratio of linoleic acid: α-linolenic acid closest to 3:1, reported as optimal for human nutrition (Leizer et al., 2000; Simopoulos, 2008; Galasso et al., 2016). These results agree and further expand the studies of Galasso et al. (2016). Similar to previous studies, the cultivar with the lowest proportion of α-linolenic acid was Eletta Campana (CAN48) (Galasso et al., 2016). However, while the mean concentration values for the six fatty acids are also within comparable ranges to previous reports (Galasso et al., 2016; Pavlovic et al., 2019), there are some discrepancies in the concentration of specific fatty acids in certain accessions. For example, we found a higher proportion of oleic acid in CAN48 than in CAN39 and the higher percentage of α-linolenic acid was found in CAN40 instead of CAN24 as reported by Galasso et al. (2016). These discrepancies can be due to different extraction methods and the fact that Galasso et al. (2016) studied the oil composition of seeds coming from individuals grown under the same conditions while our study directly analysed the seeds stored in the IPK genebank.
In addition to the main fatty acids previously mentioned, other fatty acids have also been reported in minor concentrations in hemp seed oil (Galasso et al., 2016; Pavlovic et al., 2019). These metabolites were not reported in the current paper due to identification and quantification discrepancies in these low concentrated metabolites. However, considering that the biological properties of hemp seed oil are related to the presence of the major fatty acids, especially the ratio between linoleic acid and α-linolenic acid, our results are biologically relevant. Complete reports of minor metabolites in hemp seed oil can be found in other studies (Galasso et al., 2016; Pavlovic et al., 2019).
A recent study reported the cannabinoid profile of ten commercially available hemp seed oils (Citti et al., 2019). Besides tetrahydrocannabinol and cannabidiol, another 30 cannabinoids were also identified (Citti et al., 2019). In the current study, we found traces of only three cannabinoids using a non-selective LC-MS method. Our results demonstrated that Eletta Campana (CAN48) was the cultivar with the highest concentration of cannabinoids, although their concentration was around 200 times lower than that of the most abundant metabolite (N-trans-caffeoyltyramine). Previous studies on the same accession (CAN48) obtained from the same genebank (IPK) reported an absolute quantification value of THC of 0.08%, which is well below the 0.3% threshold set by international regulations (Citti et al., 2019). Based on the low THC content of this cultivar, we might assume that all the hemp seed accessions included in the current study meet the international regulatory standards of < 0.3% THC content.
Recent studies have proved that despite being heavily studied, hemp seed continues to be a promising source of new metabolites. In the current study we discovered a previously unknown molecular family in hemp seeds “cinnamic acid glycosyl sulphates” widely distributed across all hemp seed accessions in similar concentrations. The identity of one novel compound from this chemical class was confirmed by isolation and interpretation of spectroscopic data, and four structural analogues, three of them new, is also suggested. The presence of these compounds in all hemp accessions suggest that cinnamic acid glycosyl sulphates might have important implications as quality control markers for the authentication of commercial hemp seed products extracted with highly aqueous solvents. Considering that these highly polar extracts, commonly used in the manufacture of cosmetics, are usually devoid of characteristic hemp seed phenylpropionamides due to insolubility issues, the need of more suitable quality control markers is justified. Further studies are still needed to identify the biological properties of these new metabolites.
In conclusion, our results confirm that some hemp accessions present in the IPK germplasm collection appear to contain more interesting profiles of metabolites than some of the cultivated hemp cultivars, as first suggested by Galasso et al. (2016). We suggest that accessions CAN37 and CAN58, should be scientifically prioritized for the discovery of new metabolites as well as in selective breeding programs for the development of new cultivars with high contents of bioactive phenylpropionamides. However, despite the functional and nutraceutical properties of this crop, the seeds of hemp are not totally free of antinutritional compounds, such as phytic acid (Russo, 2013; Galasso et al., 2016). According to previous reports, including ten of the accessions investigated here, CAN48 presents the highest levels of phytate, although other cultivars not included here like Futura75 and Felina32 have even higher values. The presence of high phytate contents in some cultivars limit their use in human nutrition and other monogastric animals since a high level of phytic acid may lead to mineral deficiencies of macro-and microelements, protein digestibility and poor organoleptic properties (Russo, 2013; Russo and Reggiani, 2015; Galasso et al., 2016). Therefore, an improvement for this trait might be necessary if the high yielding oil cultivar CAN48 is prioritized. According to previous studies (Galasso et al., 2016), the Italian accession CAN40 might be a good candidate to reduce the content of phytate through hybridization and selective breeding. Other accessions such as CAN32, or the Fibrimon cultivars CAN50 and CAN53, could also be prioritized if a 3:1 ratio of linoleic acid/α-linolenic acid is desired.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://massive.ucsd.edu/ProteoSAFe/static/massive.jsp, MSV000090725.
Acknowledgments
The authors would like to thank the Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), Gatersleben, Germany, for kindly proving the plant material. Li Hsun-Wen and Lin Che-Hsien (China Medical University, Taichung, Republic of China) are also acknowledged for their technical assistance preparing the LC-MS data for the heatmap analysis and to Rui Fang and Christopher Wallis (Royal Botanic Gardens, Kew) for the enlightening discussions.
Conflict of interest
The authors declare that we have received funding from companies to assess the quality of commercialized hemp seed extracts including “The Procter & Gamble Company”. However, no company contributed scientifically or financially to this research or was involved in the decision to publish it.
The authors declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1114398/full#supplementary-material