Skip to main content
Canna~Fangled Abstracts

Cultivar-dependent phenotypic and chemotypic responses of drug-type Cannabis sativa L. to polyploidization

By August 11, 2023September 7th, 2023No Comments


Abstract

Cannabis sativa L. is a plant with a wide range of potential medicinal applications. In recent years, polyploidy has gained attention as a potential strategy for rapidly improving C. sativa, which, unlike other modern crops, has not yet benefitted from this established biotechnological application. Currently, no reports on high THCA and CBDA drug-type polyploid cultivars have been published. Moreover, it still needs to be clarified if different cultivars react similarly to polyploidization. For these reasons, we set out to evaluate and compare the phenotype and chemotype of three high Δ9-tetrahydrocannabinolic acid (THCA) and one high cannabidiolic acid (CBDA) drug-type cultivars in their diploid, triploid and tetraploid state through agronomic and metabolomic approaches. Our observations on plant morphology revealed a significant increase in plant height and leaf size with increasing ploidy levels in a cultivar-dependent manner. In contrast, cannabinoids were negatively affected by polyploidization, with the concentration of total cannabinoids, THCA, CBDA and cannabigerolic acid (CBGA) decreasing significantly in higher ploidy levels across all four cultivars. Headspace analysis of volatiles revealed that total volatile content decreased in triploids. On the other hand, tetraploids reacted differently depending on the cultivars. Two THCA dominant cultivars showed an increase in concentrations, while in the other two cultivars, concentrations decreased. Additionally, several rare compounds not present in diploids appeared in higher ploidy levels. Moreover, in one high THCA cultivar, a couple of elite tetraploid genotypes for cannabinoid and volatile production were identified, highlighting the role of cultivar and genotypic variability as an important factor in Cannabis sativa L. polyploids. Overall, our observations on plant morphology align with the giga phenotype observed in polyploids of other plant species. The decrease in cannabinoids and volatiles production in triploids have relevant implications regarding their commercial use. On the other hand, this study found that tetraploidization is a suitable approach to improve Cannabis sativa L. medicinal potential, although the response is cultivar and genotype-dependent. This work lays the ground for further improving, evaluating and harnessing Cannabis sativa L. chemical diversity by the breeding, biotechnological and pharmaceutical sectors.

Keywords: Cannabis sativa L., polyploidy, cannabinoids, chemical diversity, terpenes, metabolomics, biotechnology, medicinal & aromatic plants

Introduction

Polyploidization is a widespread phenomenon in the plant kingdom and it plays a significant role in plant evolution and speciation. Indeed, genomic doubling provides additional genetic material upon which natural selection can act. At the same time, it establishes reproductive barriers, which can limit the capacity of polyploid organisms to mate with their diploid ancestors (). Polyploidization can lead to a multitude of effects on the plants’ phenology, morphology, physiology and metabolic processes. One of the most common effects of polyploidization is the “giga” effect. Cell and nucleus volume increases compared to diploid ancestors, with cell areas being 1.5 times larger (). This can lead to an increase in organ size and, ultimately, biomass (). Other frequently observed effects in different species are changes in size and number of stomata, shoots, roots, leaves, flowers, seeds and plant height (). Because of these effects, there are many examples of plant species that have improved agronomic qualities due to polyploidization. However, polyploidy does not necessarily work for the better, as documented by  and , who found that for Centella asiatica L. and Echinacea purpurea L. the number of organs was reduced, and biomass did not increase. Nevertheless, polyploidization has been a powerful breeding tool for crop improvement ().  estimated that approximately 40% of the cultivated species are polyploids, and today, common crops such as Vitis spp. L. (), Actinidia chinensis. (), Solanum tuberosum L. (Prunus domestica L. (), and Fragaria (), among many others, are polyploids.

Polyploidization has also been shown to affect secondary metabolite production in several medicinal and aromatic plants (MAPs) (), including plant species of the Cannabaceae family such as Humulus lupulus L. (). Duplication of the genome can induce a series of changes, the most obvious being the increase in gene copy number. Genome duplications can also affect the genomic organization, gene regulation, transposon activation and can induce epigenetic reprogramming by modifying histones and chromatin structures (). The final outcome is a change in gene expression which can directly influence the metabolome, either through increased, decreased, or silenced transcription ().

The general trend in polyploid MAPs is an increase in secondary metabolite production (), with some exceptions (). Various studies have reported quantitative differences of diverse secondary metabolites in synthetic autotetraploids, which could be attributed to changes in gene expression. For example, genes involved in the artemisinin biosynthetic pathway were upregulated, resulting in an increased production of artemisinin in Artemisia annua plants (), and a similar change was reported for morphine and related biosynthetic genes in Papaver somniferum L. ().

In light of these findings, Cannabis sativa L. (Cannabaceae) appears to be an interesting candidate for polyploidization as C. sativa has been widely cultivated due to its industrial (), ornamental (), nutritional (), and broad medicinal potentials () and C. sativa is a prolific producer of secondary metabolites with at least 348 well-characterized compounds () classified as cannabinoids (150), terpenoids (120), phenolics (42), flavonoids (34) and alkaloids (2).

For these reasons, lately, polyploidization has been at the center of attention in the C. sativa industry. In the past, due to legal restrictions, the generation and testing of C. sativa polyploids posed considerable challenges as permits for handling C. sativa material and hazardous chemicals as well as access to specialized equipment and technical knowledge were required. These hurdles delayed basic research needed to elucidate the changes brought about by polyploidization of C. sativa and already documented in other plant species.

To date, there are only a handful of studies assessing the polyploid state of C. sativa.  were the first to compare tetraploid to diploid C. sativa. They observed that the tetraploids had more total proteins, total flavonoids and starch while having reduced cellulose content compared to diploid plants. Unfortunately, samples were not decarboxylated before cannabinoid analysis and only Δ9-tetrahydrocannabinol (THC) and Cannabidiol (CBD) were measured therefore missing THCA and CBDA.  expanded their measurements in a second experiment. They found significant differences in leaf morphology, with tetraploids having shorter and broader leaves and a decrease in glandular trichome density, fiber content, height, and soluble sugars. At the same time, the amount of total proteins doubled compared to diploids. However, cannabinoids were measured as in their first publication, and no specifications were made regarding the starting plant material.  made observations on a single balanced CBDA/THCA drug-type tetraploid genotype and its diploid background. Measurements were conducted on ten replicates of a clone and found an 8% increase in CBDA but no increase in THCA, total cannabinoids and terpenes, while limonene decreased significantly compared to diploids. Like , they also observed that the central leaflet width was significantly bigger in tetraploids.  published a safe and efficient method using pregerminated seeds and colchicine for developing tetraploids in C. sativa. They successfully generated tetraploids of 5 drug-type CBDA dominant cultivars and were able to generate triploids via embryo rescue. However, they only presented results on stomata for which they observed an increased size and reduced density with increasing ploidy levels. Lastly,  developed triploids and tetraploids of a CBGA dominant drug-type cultivar. Although not significant (p<0.1), they observed an increase in total cannabinoid concentrations, being 8.66%, 10.18%, and 12.38% of the dry flower weight, for diploid, triploid and tetraploid, respectively. Additionally, they also observed an increase in flower weight in triploids of 23% and 26% compared to diploids in outdoor and indoor trials, respectively.

C. sativa is a highly variable crop () and , as well as , showed that the response of C. sativa to polyploidization induction is genotype and cultivar dependent. Additionally, observations of  and  on cannabinoid concentrations hinted that polyploidization might affect secondary metabolite quantity differently, likely dependent on the starting plant material. However, up to now, no tests on THCA and CBDA dominant drug-type varieties have been published, and it still needs to be clarified how different cultivars are affected by polyploidization. Therefore, we use a targeted metabolomics approach to investigate changes in cannabinoids and volatile content in triploids and tetraploids of 4 drug-type cultivars compared to their diploid ancestors. The cultivars used in the study, to fill the gaps in the field, were THCA and CBDA dominant and varied in secondary metabolite composition. Additionally, as leaf morphology has been proposed as proxy for polyploidy identification, we also measured leaf characteristics to check if previous observations apply to different cultivars. Here we set out to answer the following research questions: 1) Does ploidy level affect leaf morphology, cannabinoid and volatile quantity and diversity in C. sativa? 2) Is the response to polyploidization concerning leaf morphology, cannabinoid and volatile quantity and diversity cultivar-dependent in C. sativa?

Material and methods

Plant material

C. sativa seeds of four feminized commercial drug-type cultivars, in their respective diploid, triploid and tetraploid states, were obtained from F1 SeedTech. A pilot experiment was conducted, and among 29 cultivars tested, four cultivars were selected for their productivity, relatively low phenotypic variability, differential genetic background and chemical composition to evaluate their polyploid status. Cultivars used in the experiment were THCA dominant (cultivars A, B and C) and CBDA dominant (cultivar D).

Plant growth

Plant growth was carried out following standard agricultural practices for the production of C. sativa inflorescences, as previously described by . To minimize possible variation in the treatment of individual plants, the same nutrient solution concentration and regime of fertigation was used for all plants. Additionally, plants were placed in the growth chamber following a completely randomized design. Plants were grown in an indoor growth chamber fitted with 12, 400w HPS horticultural lights (Osram, Munich, Germany), providing a light intensity of 200 ± 50 µmol/m²/s PAR at canopy level. Environmental parameters were kept constant with a temperature of 25°C during the day and 21°C during the night, while RH was set to 70%. Diploid, triploid and tetraploid seeds were sterilized with 0.05% NaClO solution for 3 minutes and were subsequently rinsed three times with sterile deionized water and placed in Petri dishes lined with a moist paper filter. After 72h germinated seeds were then transferred to Rockwool plugs (Rockwool, Roermond, NL) moistened with hydroponic solution (Dutchpro Nutrients, Almere, NL). The plant nutrient solution was prepared by adding A and B of the vegetative hydroponic stock solution (Dutchpro) to deionized water to reach an electric conductivity (EC) of 0.4 mS/cm, and then pH was adjusted to 5.8 with NaOH. After a week from germination, plantlets were transplanted into 1l Rockwool cubes (Grodan) fitted with an automated irrigation system. During the vegetative phase, plants were fertigated three times a day with circa 150 ml vegetative hydroponic solution A and B (Dutchpro), supplemented with MgSO4 and CaCl2 to a final concentration of 0.127g/l and 0.05g/l, respectively. The EC was increased gradually as the plants grew to reach a value of 1.0 mS/cm while the surplus solution released by the Rockwool cubes was drained.

After three weeks from germination, flowering was induced by a change of photoperiod to short daylight conditions (12h/12h light/darkness). Plants then received the flowering solution A and B (Dutchpro), supplemented with MgSO4 and CaCl2, as stated above, until harvest. The EC value was increased to 1.6 mS/cm as plants entered the reproductive phase between week 7 to week 12 from germination. The nutrient solution concentration was then gradually lowered to EC 0.3 mS/cm during the last 3 weeks of flowering to promote maturation and senescence. Once all stigmas turned brown, flowers ceased to grow, and leaves progressed into senescence, plants were deemed mature. At week 15 (total of 12 weeks of flowering), flowers were harvested and dried in darkness in a climate chamber at 20°C and 30% RH for 7 days. For the number of C. sativa plants grown per cultivar and ploidy level, see Table 1. Samples for each plant were taken in equal proportion from the apical, middle and lower flowers and pooled to account for possible variability in secondary metabolite accumulation throughout the plant height. As a result, a representative 1g of dried flower material per plant was obtained. Samples were stored at -20°C in sealed 10ml Falcon tubes, and before chemical analysis, samples were grounded with a mortar and pestle until a fine powder was obtained.

Table 1

Number of Cannabis sativa L. plants grown for the experiment.

Cultivar Diploid Triploid Tetraploid
A 10 4 5
B 10 7 4
C 9 5 12
D 9 5 9

Ploidy analysis

The method for ploidy analysis was carried out following a modified protocol by . Nuclei isolation buffer containing 4.15g/l MOPS, 9.15 g/l MgCl2, 8.8 g/l TriSodium citrate, 1.55 g/l DTT, 25g/l Polyvinylpyrrolidone and 1ml/l Triton-x was prepared before the analysis and stored at 4°C. Subsequently, 100-200 mg of fresh plant material was collected from newly emerged leaves and kept on ice until analysis. The samples were prepared by chopping the leaves with a razor blade in a Petri dish containing 1 ml of buffer. The buffer containing the free nuclei was then filtered through a 20 µm Minisart nylon filter (Sartorius, Goettingen, Germany) into a 1.5 ml Eppendorf and kept on ice while the other samples were processed. After all samples were prepared, Propidium iodide was added to a final concentration of 50 µg/ml. The samples were mixed and incubated in darkness on ice for 10 minutes prior to analysis. Ploidy analysis on a flow cytometer (Milliliter Guava Easycyte, Merck, Darmstadt, Germany) was performed following analytical parameters as stated by  and . The laser at 488nm was used for excitation, and fluorescence was recorded with Yellow H linear channel at 583/26 nm. Diploid C. sativa plants were used as standards. Samples were measured in triplicate and all plants part of the experiment were tested twice throughout their development (during the vegetative and flowering phase) to ensure ploidy stability, including control diploids.

Morphological measurements

Leaf and plant height measurements were carried out two weeks before the final harvest. All plants were measured on the same day. For the leaflet measurements, the central leaflet of the leaf attached to the third node counting from the apical meristem of the plant was used. The length of the central leaflet and the width at the minor axis of the leaflet ellipse was measured. Plant height was measured from the base of the stalk to the apical meristem.

Cannabinoid analysis

Cannabinoids were extracted with 1ml of methanol from 10 mg of homogenized flower material. Sonication was carried out for 20 minutes. Subsequently, the tube was centrifuged at 13000g to pellet plant material. The supernatant was filtered through a 20 µm RC Minisart filter (Sartorius). The extract was kept in a sealed dark glass vial and stored at -20°C until analysis. Cannabinoids were quantified following a modified protocol by  with a reversed‐phase HPLC (Agilent 1200 chromatographic system, Agilent, Folsom, CA, USA) equipped with a UV-photodiode array detector (UV-DAD) and an auto-sampler. The separation was achieved on a Luna-C18 column (Phenomenex, Utrecht, NL). The mobile phase consisted of solvent A 0.1% (v/v) formic acid in water and solvent B 0.1% (v/v) formic acid in acetonitrile. The flow rate was set at 1.2 ml/minute. Solvent B was initially set at 30%, and after 3 minutes, it was increased to 75%. At 15 minutes from the start of the program, solvent B was increased to 100% and held for 3 minutes, after which it was decreased to 30% in the final two minutes of the sample run. The auto-sampler automatically injected 5μl of the sample. Light absorption was detected by the DAD at 228, 230, 280, 320 and 360 nm. The wavelength chosen for quantification was 320 nm. Quantification was performed with standards, THCA, CBDA, CBGA, THC, CBD and Cannabinol (CBN) (Merck). A five-fold serial dilution of the standards was performed on the range of quantification, and linear regression in RStudio software (version 4.3.0) was performed. R² for all calibration regressions were above 0.99.

Headspace analysis

Thirty milligrams of dried plant material were placed into a 20 mL headspace vial and sealed. The samples were incubated at 150°C for 30 minutes at 500 rpm before analyses. Before capping the vial, an aliquot of 15 µL of methyl palmitate (5mg/mL), used as internal standard, was added.

The GC-MS analysis was performed following a modified protocol () using a 7890A gas chromatograph equipped with a 7693 automatic sampler and a 5975C single-quadrupole mass detector (Agilent). Compounds were separated on a DB-5 column (30 m x 0.25 mm, 0.25 µm film, J&W Science, Folsom, CS, USA). Helium (99.9% purity) was used as a carrier gas at a flow rate of 1.6 mL/min. The oven temperature was held at 80°C for 1 minute, then increased at 3°C/min to 150°C, and then increased to 250°C at 7°C/min and held for 3 min. The samples were injected with a 20:1 split and split flow of 8 mL/min. The injection volume was 500 µL. The ionization energy in EI mode was 70 eV, and a mass scan range of 50–550 amu. Data was processed using MassHunter (B.07, Agilent). Peak identification was done by comparing the ion spectra obtained from samples to ion spectra in the NIST V.2008 library.

Statistical analysis

Two-way ANOVAs were performed with cultivar and ploidy levels as fixed factors and plant height, leaflet length, leaflet width and cannabinoid levels as dependent variables. When the normality of residuals and homoscedasticity assumptions were not met, data were Box-Cox transformed. Differences between all possible pairs of means were obtained by performing a Tukey’s Honest Significance Difference (HSD) test or a Dunnett post hoc test. Boxplots were made to visualize differences between treatments and/or cultivars.

The GC headspace chromatograms’ total peak areas (TIC) of compounds relative to that of the internal standard (methyl palmitate) were used for multivariate analysis. SIMCA-P (16.0) software (Umetrics, Umeå, Sweden) was used to perform multivariate analysis, applying univariate scaling to the data. Loading plots, as well as VIP, were set up to visualize results and predict possible biomarkers. The unsupervised method PCA, not requiring any a priori knowledge, was applied first to explore the dataset. Subsequently, a supervised method OPLS-DA was applied to identify differential components among the samples. The VIP (Variable Influence on Projection) is a parameter used for calculating the cumulative measure of the influence of individual X-variables on the model. VIP values larger than 1, point to the most relevant variables, and generally, VIP values below 0.5 are considered irrelevant variables (Galindo-Prieto et al., 2013). Thus, in the present analysis, we considered only VIP values > 1. A two-way ANOVA was performed with cultivar and ploidy level as fixed factors and total peak area as dependent variable. For the number of peaks and peak area of 33 compounds with the highest VIP score, a non-parametric Kruskal-Wallis test with cultivar and ploidy as factors was performed. Following a Kruskal-Wallis test per each individual cultivar with ploidy as factor was performed. Post hoc analysis to detect differences between all possible pairs of means was performed with a non-parametric Dunn test. The software that was used for these analyses is R version 4.3.0. Boxplots were made to visualize differences between means using R version 4.3.0.

Results

Morphological measurements

 

Plant height

Measurements of plant height, as expected, revealed significant differences between cultivars. Indeed, cultivar C did have a significantly larger height than cultivars A, B and D (Figure 1A). Also, ploidy level did affect plant height significantly, with tetraploid C. sativa plants being significantly taller (71.38 ± 15.42 cm) than diploids (61.54 ± 12.99 cm). At the same time, triploids were found in-between (65.3 ± 20.77 cm). Further analysis of separate cultivars revealed a differential response to polyploidization, with cultivar A showing a decreasing trend with increasing ploidy levels. On the other hand, triploids of cultivars C and D increased compared to diploids, with D being significantly taller than its diploid ancestors. In contrast, in cultivar B a slightly significant (0.070) difference was found between the triploids, which decreased and the tetraploids, which increased compared to diploids (Figure 1A).

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g001.jpg

Boxplots of the (A) plant height (cm), (B) length of the middle leaflet (cm) and (C) width of the middle leaflet of a leaf (cm) of diploid (2n), triploid (3n) and tetraploid (4n) plants of four Cannabis sativa L. cultivars. (A): Kruskal-Wallis rank sum test: Cultivar; χ2 = 38.269, p = 2.479e-08. Ploidy; χ2 = 7.7563, p = 0.02069. Different letters above bars indicate significant differences at p<0.05 with an all-pairs Dunn post-hoc analysis. The inserted graph displays the length of the middle leaflet length of the ploidy levels. (B) Logarithmic transformation was performed as indicated by a Box-Cox test. ANOVA: Cultivar; F3,80 = 46.12, p < 2e-16. Ploidy; F2,80 = 11.72, p = 3.42e-05. The inserted graph displays the width of the middle leaflet length of the ploidy levels. (C) Two-way ANOVA: Cultivar; F3,77 = 8.629, p = 4.98e-05. Ploidy; F2,77 = 27.244, p = 9.47e-10. The insert graphs show plant height between ploidy levels. Different letters above bars indicate significant differences at p<0.05 with a Tukey post-hoc analysis. For two samples no measurements were taken.

 

Leaf morphological traits

The middle leaflet length measurements also revealed significant differences between cultivars. Cultivars A and B had the shortest leaflet length, cultivar C had the longest leaflet length, while cultivar D was found in between. As expected, ploidy level also affected the leaflet length. Regardless of the cultivar, an increase was observed with increasing ploidy levels, with diploids (11.6 ± 2.0 cm) showing a significantly shorter leaflet length compared to triploids (13.04 ± 3.0 cm) and tetraploids (13.21 ± 2.6 cm) (Figure 1B). Both cultivar and ploidy had a significant effect on middle leaflet width. C. sativa cultivars B and C had significantly broader middle leaflets than cultivars A and D (Figure 1C). Regardless of the cultivars, the middle leaflet width increased with ploidy level, following a similar pattern observed with leaflet length. Diploids (3.17 ± 0.4 cm) showed a significantly shorter leaflet width compared to triploids (3.73 ± 0.47 cm) and tetraploids (4.09± 0.77 cm) (Figure 1C). Overall, leaf dimensions increased with ploidy level ( Figure S1 ).

Cannabinoids

As expected, for all measured cannabinoid concentrations, significant differences were found between cultivars except for CBN. In all cases, CBN levels were below the threshold of detection. These findings, in combination with the detection of very low levels of THC and CBD, indicate that degradation of the cannabinoid acids during the harvesting, sample preparation and storage was minimal. The total cannabinoid levels differed significantly between cultivars, with cultivar D, the only CBDA dominant cultivar, having a significantly lower total cannabinoid content (16.47 ± 3.98%) than cultivars A (21.43 ± 2.74%), B (21.72 ± 2.80%) and C (20.64 ± 3.62%) (Figure 2). Total cannabinoid content was influenced by ploidy level. A decreasing trend with increasing ploidy level was observed, with total cannabinoids of diploids (21.59 ± 3.25%) being significantly higher than triploids (19.19 ± 3.73%) and tetraploids (18.51 ± 4.26%) (Figure 2). The precursor cannabinoid CBGA displayed differences for both cultivars and ploidy (Figure 3). Cultivar A had the highest level (0.34 ± 0.12%), cultivar D had the lowest level (0.18 ± 0.11%), while cultivars B (0.28 ± 0.10%) and C (0.22± 0.10%) had intermediate levels. Significant differences were detected primarily between cultivar D and the three high THCA cultivars. Similar to total cannabinoids, CBGA content was significantly affected by ploidy level, with lower concentrations associated with increasing ploidy level. Diploids had significantly higher CBGA concentration (0.31 ± 0.12%) than triploids (0.24 ± 0.11%) and tetraploids (0.18 ± 0.10%) (Figure 3). Moreover, the interaction was also found significant. This significant difference was driven by the very low levels in triploids and tetraploids of cultivar D and by the differential response of cultivar A, B and C. Indeed, for cultivar D, a significant -30.8% and a -71.9% decrease was detected in triploids and tetraploids, respectively. On the other hand, only a -10.5% decrease in triploids and a slight increase of +3.2% in tetraploids was detected for cultivar B. While for cultivars C and A, although a decrease was observed, it was not found to be significant, with a -37.0% and -32.2% for cultivar C, while -13.8% and -29.5% for cultivar A in the triploid and tetraploids, respectively.

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g002.jpg

Boxplot of concentration total cannabinoids (expressed as % of dry flower weight) of a HPLC analysis of diploid (2n), triploid (3n) and tetraploid (4n) plants of four Cannabis sativa L. cultivars. Two-way ANOVA: Cultivar; F3,82 = 12.248, p = 1.06e-06. Ploidy; F2,82 = 8.672, p = 0.000383. The insert graph displays total cannabinoids of the ploidy levels: diploid, triploid and tetraploid. Different letters above bars indicate significant differences at p<0.05 with a Tukey post-hoc analysis.

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g003.jpg

Boxplot of the concentration of CBGA (expressed as % of dry flower weight) of a HPLC analysis of flowers of diploid (2n), triploid (3n) and tetraploid (4n) plants of four Cannabis sativa L. cultivars. A logarithmic transformation was performed as indicated by a Box-Cox test. Two-way ANOVA: Cultivar; F3,76 = 15.885, p = 4.12e-08. Ploidy; F2,76 = 21.644, p = 3.63e-08. Cultivar x Ploidy; F6,76 = 4.727, p = 0.000386. Different letters above bars indicate significant differences at p<0.05 with a Tukey post-hoc analysis.

As expected, the individual cannabinoids CBDA, CBD, THCA and THC significantly differed between cultivars since cultivar D is high in CBDA and very low in THCA. For this reason, a statistical analysis of THCA, THC, CBDA and CBD was performed by separating the cultivars according to the chemotype.

In THCA dominant cultivars (A, B and C), THCA concentrations were influenced only by ploidy levels. For all three cultivars a significant decrease in the triploids (18.8 ± 2.8%) compared to diploids (21.0 ± 2.8%) was found. Tetraploids (19.5 ± 3.0%) were not significantly different from diploids and triploids (Figure 4A). In the three THCA dominant cultivars, the measured CBDA content was minimal, yet an effect of polyploidization could be detected. A significant decrease with increasing ploidy level was observed. Diploids displayed the highest levels of CBDA (0.21 ± 0.09%), followed by triploids (0.18 ± 0.08%) and tetraploids (0.14 ± 0.05%) (Figure 4B). Analysis of THC did not reveal any differences, while CBD could not be detected for the three high THCA cultivars (data not shown).

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g004.jpg

Boxplot of the concentration of THCA (A) and CBDA (B) (expressed as % of dry flower weight) of a HPLC analysis of flowers of diploid (2n), triploid (3n) and tetraploid (4n) plants of three high THCA Cannabis sativa L. cultivars. (A) One-way ANOVA: Ploidy; F2,62 = 3.452, p = 0.0379. (B) An inverted square root transformation was performed as indicated by a Box-Cox test. Two-way ANOVA: Cultivar; F2,60 = 11.127, p = 0.000196. Ploidy; F2,60 = 9.879, p = 0.000196. Different letters above bars indicate significant differences at p<0.05 with a Tukey post-hoc analysis. The insert graphs display CBDA and THCA concentration of the ploidy levels: diploid, triploid and tetraploid.

CBDA concentration in the CBDA dominant cultivar (D) was also affected significantly by ploidy level. In this case, the highest ploidy level presented the lowest CBDA concentration (12.14 ± 1.19%), while triploids had intermediate concentrations (15.02 ± 4.76%) in comparison to diploids (17.59 ± 3.20%) (Figure 5). For the high CBDA cultivar no differences between ploidy levels for CBD and THCA, could be found, while THC could not be detected in cultivar D.

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g005.jpg

Boxplot of the concentration of CBDA (expressed as % of dry flower weight) of a HPLC analysis of flowers of diploid (2n), triploid (3n) and tetraploid (4n) plants of the high CBDA Cannabis sativa L. cultivar D. One-way ANOVA: Ploidy; F2,20 = 7.281, p = 0.004. Different letters above bars indicate significant differences at p<0.05 with a Tukey post-hoc analysis.

Headspace-GC-MS analysis

In total, 248 compounds were annotated (with an 85% NIST library score) from the flowers of the four C. sativa cultivars. As expected, most of the compounds annotated were monoterpenes (57) and sesquiterpenes (126). At the same time, also esters (2), alcohols (4), fatty acids (2), and cannabinoids such as Δ9-Tetrahydrocannabinol, Δ1-Tetrahydrocannabinol, Cannabichromene were also detected, while for the remaining 54 compounds the chemical family remained unassigned ( Table S1 ). A comparison of the volatile compounds showed that cultivar C was the most diverse, containing the largest number of compounds detected (81 ± 14 peaks) compared to A (72 ± 10 peaks) and B (71 ± 14 peaks), while cultivar D was found in-between (77 ± 13 peaks). Interestingly no effect of ploidy or a ploidy x cultivar interaction was detected for number of peaks (Figure 6A). Cultivar B showed the highest concentration of volatiles followed by C, A and D (Figure 6B). The interaction between ploidy and cultivar was found significant while the effect of ploidy was detected as almost significant (p = 0.071). A decrease of -21.89% was found in the triploids compared to diploids (Figure 6B). While analysis of the interaction revealed a differential response of the cultivars to polyploidization regarding volatile concentrations. The main differences were detected between tetraploids of cultivar B and C, presenting the highest overall concentrations compared to tetraploids of cultivar A and tetraploids and triploids of cultivar D, which displayed the lowest overall concentrations. Additionally, a significant decrease was also detected between the diploids and tetraploids of cultivar D (-62.7%), while although not significant in cultivar B, tetraploids displayed a +75.7% increase in volatile concentrations compared to B diploids (Figure 6B).

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g006.jpg

(A) Boxplot of total number of peaks of a headspace GC-MS analysis of diploid (2n), triploid (3n) and tetraploid (4n) plants of four Cannabis sativa L. cultivars. Kruskal-Wallis rank sum test: Cultivar; χ2 = 12.624, p = 0.005. Ploidy; χ2 = 0.8091, p = 0.667. Different letters above bars indicate significant differences at p<0.05 with a Kruskal-Wallis and a post-hoc all-pairs Dunn test. (B) Boxplot of total peak area per mg dry weight of a headspace GC-MS analysis of diploid, triploid and tetraploid plants of four Cannabis sativa L. cultivars. A logarithmic transformation was performed as indicated by a Box-Cox test. Two-way ANOVA: Cultivar; F3,76 = 6.558, p = 0.00053. Ploidy; F2,76 = 2.741, p = 0.071. Cultivar x Ploidy; F6,76 = 4.082, p = 0.001. Different letters above bars indicate significant differences at p<0.05 with a Tukey post-hoc analysis.

To investigate the changes in volatiles, a principal component analysis (PCA) was carried out. The multivariate analysis showed a trend in the sample separation driven by C. sativa cultivars ( Figure S2A ). The main separation along the PC1 was due mainly to cultivars C and D, while the PCA plot did not show separation among the three ploidy levels ( Figure S2B ). To better discriminate cultivars and ploidy level effects, a supervised discriminant analysis, Orthogonal partial least square analysis-discriminant analysis (OPLS-DA), was used. The OPLS-DA model showed separation between the cultivars C and D while A and B clustered together (Figure 7A). The model was highly validated (R2X= 0.92, Q2 = 0.75 and p-value in CV-ANOVA= 0). Using the same OPLS-DA model, but differentiating the samples based on ploidy levels, a gradient separation was achieved, though the model was not validated R2X= 0.45, Q2 = 0.07 and p-value in CV-ANOVA = 1 (Figure 7B).

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g007.jpg

Orthogonal partial least square analysis-discriminant analysis of GC-MS analysis of 4 Cannabis sativa L. cultivars A, B, C and D (A). The model was highly validated (R2X= 0.92, Q2 = 0.75 and p value in CV-ANOVA= 0). OPLS-DA on ploidy levels (B). The model was not validated (R2X= 0.45, Q2 = 0.07 and p value in CV-ANOVA= 1).

Given the different volatile compositions of the cultivars, we therefore also analyzed the effect of ploidy level on each cultivar separately with OPLS-DA (Figure 8). Although this strategy increased the resolution of ploidy separation in all four cultivars, only for cultivar B the model was validated with a Q2 value larger than 0.4 (Figure 8B).

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g008.jpg

Orthogonal partial least square analysis-discriminant analysis of GC-MS analysis of four Cannabis sativa L. cultivars based on ploidy levels (A) cultivar A: Q2 0.24, (B) cultivar B: Q2 0.44, (C) cultivar C: Q2 0.04 and (D) cultivar D: Q2 0.27.

To further investigate the changes in volatiles caused by the different ploidy levels on individual compounds, 25 compounds were extracted from the model separating the cultivars. In addition, 29 compounds with a VIP score larger than one were extracted from the model separating the ploidy levels, although this model was not validated. Checking for overlapping compounds between models revealed that of the 29 extracted from the ploidy model, 21 were already present in the cultivar model. This resulted in a pool of 33 individual compounds, which except for one case (Butanoic acid, hexyl ester), concerned monoterpenes and sesquiterpenes (Figures 9,10).10). Further investigation revealed that changes in ploidy levels greatly influenced these individual compounds and often, cultivars were affected differently (Tables 23).

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g009.jpg

Boxplots of the total peak area per mg dry weight of 25 compounds with the highest VIP score in an OPLS-DA analysis separating cultivars of a headspace GC-MS analysis of diploid (2n), triploid (3n) and tetraploid (4n) plants of four Cannabis sativa L. cultivars. Kruskal-Wallis rank sum tests were performed for all compounds, for results of the tests see corresponding Table 2. Different letters above bars indicate significant differences at p<0.05 with a Kruskal-Wallis within the cultivar and between cultivars and post hoc all pairs Dunn test. (A) Fenchone, (B) Linalool, (C) Germacrene B, (D) L-α-Terpineol, (E) Fenchol, (F) γ-selinene, (G) Selina-3,7(11)-diene, (H) β-Maaliene, (I) cis-2-Norbornanol, (J) 3-Eudesmen-11-ol, (K) Butanoic acid, hexyl ester, (L) Caryophyllene, (M) Bicyclogermacrene, (N) β-Selinene, (O) Guaiol, (P) 1,4,7,-Cycloundecatriene,1,5,9,9-tetramethyl-,Z,Z,Z-, (Q) Ylangene, (R) 3-Methylcamphenilol, (S) β-Bisabolene, (T) D-Limonene, (U) Selina-5,11-diene, (V) Eudesm-7(11)-en-4-ol, (W) cis-α-Bergamotene, (X) γ-Eudesmol, (Y) (-)-Guaia-6,9-diene.

Table 2

Results of Kruskal-Wallis rank sum tests of the total peak area per mg flower dry weight of 25 compounds with the highest VIP score in an OPLS-DA analysis of a headspace GC-MS analysis of diploid, triploid and tetraploid plants of four different Cannabis sativa L. cultivars.

Fig Common name Class Cultivar Ploidy Ploidy levels per cultivar
A B C D
χ2 p χ2 p χ2 p χ2 p χ2 p χ2 p
A Fenchone Monoterpene 36.7 0.000 4.9 0.087 0.8 0.661 0.1 0.961 2.0 0.362 14.2 0.001
B Linalool Monoterpene 48.1 0.000 0.9 0.649 4.4 0.112 0.8 0.673 1.8 0.400 1.9 0.384
C Germacrene B Sesquiterpene 48.1 0.000 2.8 0.251 8.2 0.016 4.4 0.111 3.5 0.175 3.4 0.187
D L-α-Terpineol Monoterpene 47.2 0.000 1.6 0.457 2.5 0.283 9.3 0.009 5.4 0.067 4.1 0.128
E Fenchol Monoterpene 37.0 0.000 8.7 0.013 6.0 0.050 8.1 0.018 5.7 0.059 12.0 0.002
F γ-selinene Sesquiterpene 37.4 0.000 0.7 0.718 7.0 0.030 0.7 0.704 2.1 0.356 2.6 0.275
G Selina-3,7(11)-diene Sesquiterpene 34.7 0.000 2.6 0.268 7.8 0.020 4.8 0.093 2.7 0.264 5.3 0.070
H β-Maaliene Sesquiterpene 35.0 0.000 0.1 0.973 6.7 0.035 8.1 0.018 3.4 0.186 4.3 0.115
I cis-2-Norbornanol Monoterpene 30.6 0.000 6.4 0.039 0.9 0.633 3.7 0.155 4.8 0.089 3.5 0.171
J 3-Eudesmen-11-ol Sesquiterpene 26.2 0.000 3.3 0.195 1.2 0.543 5.5 0.063 1.3 0.530 16.0 0.000
K Butanoic acid, hexyl ester Ester 26.9 0.000 0.3 0.841 2.7 0.256 1.9 0.395 NA NA 1.2 0.541
L Caryophyllene Sesquiterpene 43.8 0.000 2.0 0.368 3.7 0.150 7.5 0.023 9.8 0.007 16.9 0.000
M Bicyclogermacrene Sesquiterpene 39.8 0.000 2.1 0.348 0.3 0.830 4.2 0.120 2.2 0.338 1.5 0.463
N β-selinene Sesquiterpene 41.6 0.000 1.0 0.589 6.2 0.044 2.6 0.275 1.5 0.475 1.4 0.495
O Guaiol Sesquiterpene 20.4 0.000 2.7 0.255 2.9 0.237 3.5 0.171 1.3 0.534 14.9 0.001
P 1,4,7,-Cycloundecatriene, 1,5,9,9-tetramethyl-, Z,Z,Z- Sesquiterpene 15.9 0.001 1.3 0.509 3.2 0.204 1.7 0.437 1.8 0.399 3.1 0.217
Q Ylangene Sesquiterpene 39.9 0.000 1.7 0.419 6.3 0.043 0.2 0.884 8.3 0.016 1.0 0.613
R 3-Methylcamphenilol Monoterpene 25.5 0.000 4.2 0.121 1.4 0.491 2.3 0.313 2.4 0.302 3.7 0.158
S β-Bisabolene Sesquiterpene 24.0 0.000 4.2 0.121 0.9 0.638 NA NA NA NA 7.2 0.027
T D-Limonene Monoterpene 26.2 0.000 7.4 0.024 5.0 0.082 3.1 0.218 4.0 0.135 12.2 0.002
U Selina-5,11-diene Sesquiterpene 18.9 0.000 4.1 0.127 4.2 0.120 0.7 0.714 2.3 0.317 5.4 0.066
V Eudesm-7(11)-en-4-ol Sesquiterpene 19.6 0.000 4.4 0.108 1.9 0.389 3.0 0.224 2.8 0.248 9.2 0.010
W cis-α-Bergamotene Sesquiterpene 19.5 0.000 0.6 0.705 3.8 0.150 2.4 0.299 1.0 0.618 1.7 0.417
X γ-Eudesmol Sesquiterpene 10.3 0.016 4.4 0.108 1.9 0.389 3.0 0.224 2.8 0.248 9.2 0.010
Y (-)-guaia-6,9-diene Sesquiterpene 31.9 0.000 4.2 0.117 1.9 0.387 1.9 0.387 3.4 0.179 1.6 0.459

Compounds are ranked with the first compound having the highest VIP score. Letters in the first row refer to corresponding boxplots in Figure 9. NA indicates that the compounds were not detected in the cultivar and therefore a statistical test could not be performed.

Table 3

Results of Kruskal-Wallis rank sum tests of the total peak area per mg flower dry weight of 8 compounds with the highest VIP score in an OPLS-DA analysis of a headspace GC-MS analysis of diploid, triploid and tetraploid plants.

Fig Common name Class Cultivar Ploidy Ploidy levels per cultivar
A B C D
χ2 p-value χ2 p-value χ2 p-value χ2 p-value χ2 p-value χ2 p-value
A α-Costol Sesquiterpene 10.5 0.015 3.7 0.154 1.900 0.387 1.222 0.543 3.192 0.203 NA NA
B α-Amorphene Sesquiterpene 5.6 0.131 3.6 0.160 1.608 0.448 1.807 0.405 4.987 0.083 1.556 0.459
C γ-Cadinene Sesquiterpene 16.3 0.001 2.2 0.324 0.916 0.633 0.586 0.746 1.245 0.537 3.600 0.165
D cis-Caryophyllene Sesquiterpene 26.3 0.000 4.9 0.084 NA NA NA NA 4.146 0.126 3.253 0.197
E Citral Monoterpene 11.7 0.008 6.2 0.045 NA NA NA NA NA NA 7.196 0.027
F Guai-1(10)-en-11-ol Sesquiterpene 8.1 0.045 4.0 0.134 2.735 0.255 2.385 0.303 2.557 0.279 8.273 0.016
G β-Cymene Monoterpene 11.7 0.008 3.2 0.199 NA NA NA NA NA NA 3.613 0.164
H γ-Muurolene Sesquiterpene 21.0 0.000 4.3 0.115 1.943 0.379 4.079 0.130 5.416 0.067 1.553 0.460

Compounds are ranked with the first compound having the highest VIP score. Letters in the first row refer to corresponding boxplots in Figure 10. NA indicates that the compounds were not detected in the cultivar and therefore a statistical test could not be performed.

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g010.jpg

Boxplots of the total peak area per mg dry weight of 8 compounds with the highest VIP score in an OPLS-DA analysis separating ploidy levels of a headspace GC-MS analysis of diploid (2n), triploid (3n) and tetraploid (4n) plants of four Cannabis sativa L. cultivars. Kruskal-Wallis rank sum tests were performed for all metabolites, for their results see Table 3(A) α-Costol, (B) α-Amorphene, (C) γ-Cardinene, (D) cis-Caryoplhyllene, (E) Citral, (F) Guai-1(10)-en-11-ol, (G) β-Cymene, (H) γ-Muurolene. Different letters above bars indicate significant differences at p<0.05 with a Kruskal-Wallis and an all-pairs Dunn test within the cultivar.

Analysis of the 33 compounds in triploids revealed the same pattern observed in total volatile content but on a deeper level. In triploid plants, the majority of compounds decreased with cultivar A, showing a decrease in the concentration of 22 compounds while 8 increased compared to diploids. In cultivar B, 21 compounds decreased while 7 increased. In cultivar C, 18 compounds decreased while 11 increased and in cultivar D the concentration of 27 compounds decreased while only 5 increased compared to diploids ( Table S2 ). Overall, this data indicates that although a small proportion of compounds increased in concentration, the majority decreased, with cultivar D being the most affected and cultivar C the least ( Table S2 ).

An interesting trend was observed in tetraploids, with cultivars B and C showing an important difference compared to cultivars A and D. In cultivar A the concentrations of 24 compounds decreased while 6 increased compared to diploids. Cultivar D also showed a decrease in concentrations of 26 compounds and an increase in 5 compared to diploids. In contrast, in cultivar B, the concentration of only 8 compounds decreased while 18 increased. Similarly in cultivar C the concentration showed a reduction in 9 compounds while 20 increased compared to diploids ( Table S2 ).

Additionally, several compounds present in diploids were no longer detected in triploids and tetraploids ( Table S2 ). In cultivar A, β-bisabolene, Guaia-6,9-diene, α-Costol were not detected in both triploids and tetraploids. The compounds γ-Cadinene, Eudesm-7(11)-en-4-ol and 1,4,7,-Cycloundecatriene, 1,5,9,9-tetramethyl-, Z,Z,Z- were not detected in triploids and Selina-5,11-diene was not detected in tetraploids A. In cultivar B α-Costol was not detected in both triploids and tetraploids and Eudesm-7(11)-en-4-ol and α-Amorphene were not detected in triploids while 3-Eudesmen-11-ol, Guaiol, cis-α-Bergamotene, Guai-1(10)-en-11-ol and γ-Muurolene were not detected in tetraploids. In cultivar C only α-Costol disappeared in triploids. On the other hand, several compounds were no longer detected in both triploids and tetraploids of cultivar D, Eudesm-7(11)-en-4-ol, Guaia-6,9-diene, α-Amorphene and cis-Caryophyllene while γ-Muurolene and Butanoic acid, hexyl ester were not present in triploids and Selina-5,11-diene was not present in tetraploids ( Table S2 ).

Analysis across cultivars revealed that Selina-5,11-diene was not present in both tetraploids of cultivars A and D, Eudesm-7(11)-en-4-ol was not present in triploids of A, B and D, Guaia-6,9-diene was not present in both triploids and tetraploids of A and D. α-Costol was not present in triploids of A, B, C and tetraploids A and B while α-Amorphene was not present in triploids of A, B and D.

In cultivars C and D, some compounds that were not present in diploids were detected in the higher ploidy levels. In cultivar C, α-Amorphene was present in triploids and tetraploids, in concentrations that were more than doubled in tetraploids compared to triploids. Interestingly in cultivar D Citral and β-Cymene were present in both triploids and tetraploids while γ-Cadinene only was present in triploids.

Statistical analysis revealed that, for four compounds, significant differences in ploidy level were detected when pooling the four cultivars. Three compounds (cis-2-Norbornanol, Fenchol, and D-Limonene) showed the same pattern with strikingly similar numbers. Indeed, a significant decrease in triploids (cis-2-Norbornanol -55.31%; Fenchol -43.86%; D-Limonene -47.22%) and a decrease in tetraploids (cis-2-Norbornanol -20.94%; Fenchol -18.91%; D-Limonene -27.54%) was observed (Figure 9). On the other hand, Eudesm-7(11)-en-4-ol showed an increase in tetraploids (+85.39%) while there was a decrease in triploids (-81.77%), leading to a significant difference between triploids and tetraploids plants. Further analysis on Eudesm-7(11)-en-4-ol revealed that in triploids of cultivar A, B and D, the compound was undetectable, while the increase in the tetraploids was driven by cultivar B (+111.23%) and C (+112.92%) (Figure 9).

Fenchone decreased in both ploidy levels across the four cultivars, apart from tetraploid C, where a slight increase of +5.45% was detected. Moreover, the decrease in cultivar D was found to be significant, with -81.52% in triploids and -91.73% in tetraploids. β-Bisabolene was undetectable in cultivars B and C, while in cultivar D, it significantly decreased in both triploids (-88.30%) and tetraploids (-94.43%) and decreased under detection level in both triploids and tetraploids of cultivar A. The remaining compounds all showed a differential response to polyploidization depending on the cultivars under analysis, particularly comparing tetraploids of cultivar A and D versus B and C ( Table S2 ).

Some of the most striking differences were detected in β-Maaliene, where a significant decrease was detected in tetraploids of cultivar A (-57.67%), while for cultivar B a significant increase was detected in tetraploids (+164.82%). Ylangene also showed a differential response of the cultivars in the tetraploid state, with cultivar A showing a significant decrease of -72.75%. In contrast, a significant increase of +182.13% for cultivar C was observed. Similarly, for β-Selinene, cultivar A showed a significant decrease in tetraploids (-56.66%) while cultivar B, although not significantly, increased by +113.96% in tetraploids. For γ-Eudesmol, cultivar D showed a significant -69.06% decrease in tetraploids, while cultivar A showed a +105.39% increase. Analysis of cis-α-Bergamotene revealed that concentrations in tetraploids of cultivar A increased by +78.04%, in cultivar C by +83.09%, while in tetraploids of cultivar B, the compound was not detected anymore. Interestingly, the rare compound Guaia-6,9-diene increased in tetraploids of cultivar C by +145.45% and +3378.74% in tetraploids of cultivar B. On the other hand, the compound could not be detected anymore in tetraploids of cultivars A and D.

Since cultivar B showed for several traits a positive reaction to polyploidization and the OPLS-DA model was validated, the analysis was also performed on the VIPs compounds of cultivar B separately (Table 4). The analysis was performed on 15 compounds, of which 8 showed a significant increase in cultivar B tetraploids compared to diploids (β-Maaliene +164.83%; Endo-Borneol +120.50%; Caryophyllene +118.4%; 3-Carene +1874.68%; L-α-Terpineol +150.87%; Fenchol +59.39%; β-Thujene +692.22%; Trans-2-Pinanol +309.3%) while two increased although the statistics could not detect a significant effect of ploidy (D-Guaiene +298.38%; Bicyclogermacrene +274.77%) (Figure 11).

An external file that holds a picture, illustration, etc.
Object name is fpls-14-1233191-g011.jpg

Boxplots of the total peak area per mg flower dry weight of 8 compounds with the highest VIP score in an OPLS-DA analysis separating ploidy levels of a headspace GC-MS analysis of diploid (2n), triploid (3n) and tetraploid (4n) plants of Cannabis sativa L. cultivar B. ANOVAs were performed to assess the differences between ploidy levels for this cultivar in the case of non-parametric data the non-parametric alternative Kruskal-Wallis rank sum tests was performed, for results see Table 4. Different letters above bars indicate significant differences at p<0.05 of an ANOVA and Tukey post-hoc analysis or a Kruskal-Wallis and an all-pairs Dunn test. (A) β-Maaliene, (B) endo-Borneol, (C) Caryophyllene, (D) 3-Carene, (E) L-α-Terpineol, (F) Fenchol, (G) β-Thujene, (H) trans-2-Pinanol.

Table 4

Results of ANOVAs or Kruskal-Wallis rank sum tests of the total peak area per mg dry flower weight of 8 compounds with the highest VIP score in an OPLS-DA analysis of a headspace GC-MS analysis of diploid, triploid and tetraploid plants of cultivar B. Letters in the first row refer to corresponding boxplots in Figure 11.

Fig. Common name Class Ploidy Ploidy
F-value p-value χ2 p-value
A β -Maaliene Sesquiterpene 8.565 0.003
B endo-Borneol Monoterpene 10.040 0.001
C Caryophyllene Sesquiterpene 8.460 0.003
D 3-Carene Monoterpene 10.723 0.005
E L- α-Terpineol Monoterpene 10.980 0.001
F Fenchol Monoterpene 6.562 0.008
G β -Thujene Monoterpene 9.505 0.009
H trans-2-Pinanol Monoterpene 5.776 0.056

Additionally, analysis of outliers detected an interesting effect of polyploidization in cultivar C. Several genotypes in this cultivar consistently presented very high levels of rare compounds compared to diploid ancestors, identifying them as elite tetraploid genotypes. Indeed, plant ID 82 presented a large increase in 8 compounds compared to diploid (β-Maaliene +336.89%; Caryophyllene +354.78%; 1,4,7,-Cycloundecatriene, 1,5,9,9-tetramethyl-, Z,Z,Z- +504.81%; Ylangene +830.39%; β-Selinene +145.86%; Guaia-6,9-diene +1024.29%; α-Amorphene +inf%; α-Costol +932.78%) while plant ID 91 showed a marked increase in 3 rare compounds (Guaiol + 569.55%; cis-Caryophyllene +913.95%; 3-Eudesmen-11-ol +527.32%). Moreover, further investigation revealed that plant ID 82 was also a top performer for total volatile content being sixth out of 89 plants. Given these results, cannabinoid data was also screened, and C tetraploid genotype ID 85 was found to be the best performer for total cannabinoids (26.84%) and linalool (+292% compared to average diploid).

Discussion

In recent years, polyploidy has gained attention as a potential strategy for rapidly improving C. sativa. Compared to other modern crops, C. sativa has yet to benefit from this established biotechnological application. Indeed, the observed effects on plant morphology and secondary metabolite production in other medicinal and aromatic plants make it a very promising field of research.

Here we found that contrary to previous observations (), plant height significantly increased with increasing ploidy levels. However, analysis of separate cultivars also highlighted a differential response, with cultivar D showing significantly taller triploids than diploids, while B triploids were shorter than their diploids counterpart. On the other hand, leaf morphology observations agreed with previous observations (). Leaf size increased significantly, with both leaflet length and width increasing with increasing ploidy level. Overall, our observations also align with the “giga” phenotype observed in other plant species after polyploidization ().

Investigation of secondary metabolites revealed several surprising changes. First, cannabinoids appear to be negatively affected by polyploidization, as the concentration of total cannabinoids decreased in higher ploidy levels. Moreover, major cannabinoids THCA, CBDA, as well as CBGA followed the same trend as observed for total cannabinoids. The significant decrease in concentration observed in total cannabinoids, -13.27% in triploids and – 28.89% in tetraploids of the CBDA dominant cultivar, and -11.76% in triploids and -6.89% in tetraploids of THCA dominant cultivars, has relevant implications regarding the uses of polyploids for producing major cannabinoids. These results are in contrast with the observation by , which observed an increase in total cannabinoids in polyploids of a drug-type CBGA dominant cultivar, highlighting once more the cultivar-dependent response to polyploidization. Regarding the concentration of CBGA, the results agree with the observation of  who also detected a circa -30% decrease. Moreover, considering that the limitation on the availability of the precursor cannabinoid CBGA is a strong limiting factor for the accumulation of downstream cannabinoids (). These data point to the probable downregulation of the cannabinoid pathway in the higher ploidy levels. From these results, it appears that in modern drug-type C. sativa cultivars, where cannabinoid production has already been pushed to extremely high levels (>20% of dry flower weight) by artificial selection, in an event such as polyploidization, where plants undergo a highly stressful genomic reorganization, changes might favor other metabolic pathways as observed by .

A much more complex suite of changes was detected regarding the presence and concentration of volatiles. First of all, triploids showed a clear decrease in total volatile content. On the other hand, tetraploids depending on the cultivars, reacted differently, with cultivars B and C affected positively, while in cultivars A and D, concentrations decreased.

The VIPs extracted from the multivariate analysis yielded the 33 individual compounds having the most influence on the separation between cultivars and ploidy levels. Analysis of these compounds highlighted to a deeper level several interesting findings. Generally, triploids’ volatiles were negatively affected across cultivars, apart from specific cases. Plants in the tetraploid state showed a differential response, with many compounds in C and B increasing while at the same time in A and D decreasing. Interestingly, several rare compounds appeared in higher ploidy levels. Cultivar C presented a couple of elite tetraploid genotypes, highlighting the role of within cultivar genotypic variability as an important factor in C. sativa polyploids. In cultivars A and D, the concentrations of most of the compounds decreased, and several were not detected in higher ploidy levels, pointing to the downregulation of the terpene pathway in these cultivars. On the other hand, in cultivar B, although an increase in concentration was detected for most of the compounds at higher ploidy levels, some also disappeared. These changes point to the probable gene-silencing effects caused by polyploidization in this specific cultivar ().

This complex suite of changes have important implications for the medicinal value of the plants. Indeed, terpenes play an essential role in the bioactivity of C. sativa extracts and can contribute to pharmacological activity via the entourage effect (). Some compounds that were found to increase significantly in cultivars B and C display important medicinal properties. “Endo-Borneol can increase drug delivery across various physiological barriers (), 3-Carene has sleep-enhancing effects by targeting the GABAA-benzodiazepine receptors (), L-α-Terpineol displays anticancer, anticonvulsant, antiulcer, antihypertensive, anti-nociceptive properties () and β-Caryophyllene might become effective for the treatment of diabetes and associated complications (). Moreover, Guaia-6,9-diene, a rare compound with promising anti-SARS-CoV-2 properties, also greatly increased ().

Likewise, a decrease in cultivars A and D of compounds with relevant medicinal properties would probably affect the cultivars’ medicinal potential. The concentration of some compounds which decreased in higher ploidy levels show promising antitumor activity, such as Guaiol and β-bisabolene () and the ability to counteract the accumulation of Amyloid-beta in Alzheimer’s disease, Fenchol ()

To the best of our knowledge, this is the first report applying a metabolomics approach to investigate C. sativa polyploids in their diploid, triploid and tetraploid status across different cultivars. Overall, the response to polyploidization on morphological traits aligns with the observations on other plant species, showing larger organs with increasing ploidy levels. Moreover, analysis of secondary metabolites revealed for cannabinoids a marked decrease across cultivars in the higher ploidy levels. On the other hand, volatiles were negatively affected in triploids, while a differential response based on the cultivars’ genetic background was observed in tetraploids.

Our observations revealed that the effect of gene doubling appears to be more evident on minor compounds, making polyploidy a very promising tool to increase their concentrations. C. sativa presents at least 150 minor cannabinoids apart from THCA and CBDA () and many other secondary metabolites () with medicinal potential () that were not considered in this study and could be affected in the tetraploid state. Further research on C. sativa tetraploids applying metabolomic approaches could reveal changes also in minor cannabinoids and other metabolite classes with potential biological activity and medicinal properties. Research on preclinical and animal models pointed out that C. sativa botanical preparations appear more effective than pure THC () or CBD (). Yet it is not clear what compounds or combination of compounds are responsible for such increased activity. Given the shift observed in secondary metabolite quantity and diversity, follow-up tests on tetraploid C. sativa botanical preparations could be a first step to reveal what kind of changes polyploidization bring about in the plant’s pharmacological potential. Moreover, the only CBDA dominant cultivar of the study showed the most severe changes in cannabinoid and volatile reduction in the higher ploidy levels. Follow-up studies with more cultivars might elucidate what kind of relationship, if present, there might be between the effects of polyploidization and the genetic background of the plants.

The research on C. sativa polyploids is still in its infancy, and with this study, we have just scratched the surface. More research applying genomics, transcriptomics, metabolomics and phytocannabinomics () could reveal significant changes induced by genome doubling and elucidate mechanisms by which secondary metabolites are influenced in C. sativa polyploids.

Overall, this study found that polyploidization of C. sativa is a suitable approach to improve its medicinal potential, while the response is cultivar and genotype-dependent. The findings and methodology of this research lay the ground for further improving, evaluating and harnessing C. sativa chemical diversity by the breeding, biotechnological and pharmaceutical sectors.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

HPF: Headspace analysis, data acquisition, multivariate analysis, manuscript writing, review & editing. YHC: Supervising chemical analysis & manuscript review. KV: Conceptualization, project administration, formal analysis, supervision, manuscript writing, review & editing. MB: Conceptualization. BS: Plant growth, cannabinoid analysis, ploidy analysis, data acquisition, statistical analysis, formal analysis, manuscript writing, review & editing. FT: Conceptualization, supervision, methodology, plant growth, ploidy analysis, cannabinoid analysis, data acquisition, statistical analysis, formal analysis, investigation, manuscript writing, review & editing. All authors contributed to the article and approved the submitted version.

Acknowledgments

We thank Lina Bayona Maldonado and Hye Kyong Kim for help with the method development of the cannabinoids analysis and F1 SeedTech for the financial support of the project.

Funding Statement

The authors declare that this study received funding from F1 SeedTech. The funder provided the seeds of the cultivars studied. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.

Conflict of interest

Author MB is employed by the company F1 SeedTech.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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.1233191/full#supplementary-material

Supplementary Figure 1

The picture shows a side branch of a diploid (left side) and tetraploid (right side) genotypes of cultivar A.

Supplementary Figure 2

Principal component analysis PC1 against PC2, using UV scaling preprocessing data, of 4 Cannabis sativa L. cultivars A–D (A) and three ploidy levels 2n, 3n and 4n (B). The PCA is based on 248 identified compounds. The explained variation is R2X= 0.28.

References

  • Amparo T. R., Seibert J. B., Silveira B. M., Costa F. S. F., Almeida T. C., Braga S. F. P., et al.. (2021). Brazilian essential oils as source for the discovery of new anti-COVID-19 drug: a review guided by in silico studyPhytochem. Rev. 20, 1013–1032. doi:  10.1007/S11101-021-09754-4 [PMC free article] [PubMed] [CrossRef[]
  • Andre C. M., Hausman J. F., Guerriero G. (2016). Cannabis sativa: The plant of the thousand and one moleculesFront Plant Sci. 7, 19. doi:  10.3389/fpls.2016.00019 [PMC free article] [PubMed] [CrossRef[]
  • Bennett M. D., Leitch I. J. (2005). Nuclear DNA amounts in angiosperms: progress, problems and prospectsAnn. Bot. 95, 45–90. doi:  10.1093/AOB/MCI003 [PMC free article] [PubMed] [CrossRef[]
  • Blasco-Benito S., Seijo-Vila M., Caro-Villalobos M., Tundidor I., Andradas C., García-Taboada E., et al.. (2018). Appraising the “entourage effect”: Antitumor action of a pure cannabinoid versus a botanical drug preparation in preclinical models of breast cancerBiochem. Pharmacol. 157, 285–293. doi:  10.1016/J.BCP.2018.06.025 [PubMed] [CrossRef[]
  • Carputo D., Frusciante L., Peloquin S. J. (2003). The role of 2n gametes and endosperm balance number in the origin and evolution of polyploids in the tuber-bearing SolanumsGenetics 163, 287–294. doi:  10.1093/GENETICS/163.1.287 [PMC free article] [PubMed] [CrossRef[]
  • Cerrato A., Citti C., Cannazza G., Capriotti A. L., Cavaliere C., Grassi G., et al.. (2021). Phytocannabinomics: Untargeted metabolomics as a tool for cannabis chemovar differentiationTalanta 230, 122313. doi:  10.1016/J.TALANTA.2021.122313 [PubMed] [CrossRef[]
  • Clarke R. C., Merlin M. D. (2016). Cannabis domestication, breeding history, present-day genetic diversity, and future prospectsCRC Crit. Rev. Plant Sci. 35, 293–327. doi:  10.1080/07352689.2016.1267498 [CrossRef[]
  • Crawford S., Rojas B. M., Crawford E., Otten M., Schoenenberger T. A., Garfinkel A. R., et al.. (2021)Characteristics of the diploid, triploid, and tetraploid versions of a cannabigerol-dominant F1 hybrid industrial hemp cultivar, Cannabis sativa “Stem cell CBG.” Genes (Basel) 12, 923 doi:  10.3390/GENES12060923 [PMC free article] [PubMed] [CrossRef[]
  • Eckardt N. A. (2010). A double lock on polyploidy-associated epigenetic gene silencingPlant Cell 22, 3. doi:  10.1105/TPC.110.220112 [PMC free article] [PubMed] [CrossRef[]
  • Edger P. P., McKain M. R., Yocca A. E., Knapp S. J., Qiao Q., Zhang T. (2020). Reply to: Revisiting the origin of octoploid strawberryNat. Genet. 52, 1 52, 5–7. doi:  10.1038/s41588-019-0544-2 [PMC free article] [PubMed] [CrossRef[]
  • Gallily R., Yekhtin Z., Hanuš L. (2015). Overcoming the bell-shaped dose-response of cannabidiol by using Cannabis extract enriched in cannabidiolPharmacol. Pharm. 6, 75–85. doi:  10.4236/pp.2015.62010 [CrossRef[]
  • Gul W., Gul S. W., Radwan M. M., Wanas A. S., Mehmedic Z., Khan I. I., et al.. (2015). Determination of 11 cannabinoids in biomass and extracts of different varieties of cannabis using high-performance liquid chromatographyJ. AOAC Int. 98 (6), 1523–1528. doi:  10.5740/jaoacint.15-095 [PubMed] [CrossRef[]
  • Hanuš L. O., Meyer S. M., Muñoz E., Taglialatela-Scafati O., Appendino G. (2016). Phytocannabinoids: a unified critical inventoryNat. Prod Rep. 33 (12), 1357–1392. doi:  10.1039/c6np00074f [PubMed] [CrossRef[]
  • Hashiesh H. M., Nagoor Meeran M. F., Sharma C., Sadek B., Kaabi J., Al, Ojha S. K. (2020). Therapeutic potential of β-caryophyllene: A dietary cannabinoid in diabetes and associated complicationsNutrients 12, 1–30. doi:  10.3390/NU12102963 [PMC free article] [PubMed] [CrossRef[]
  • Hesami M., Pepe M., Baiton A., Salami S. A., Jones A. M. P. (2022). New insight into ornamental applications of cannabis: perspectives and challengesPlants (Basel). 11 (18), 2383. doi:  10.3390/plants11182383 [PMC free article] [PubMed] [CrossRef[]
  • Hopping M. E. (1994). Flow cytometric analysis of Actinidia speciesNew Z. J. Bot. 32 (1), 85–93. doi:  10.1080/0028825X.1994.10410410 [CrossRef[]
  • Iannicelli J., Guariniello J., Tossi V. E., Regalado J. J., Di Ciaccio L., van Baren C. M., et al.. (2020). The “polyploid effect” in the breeding of aromatic and medicinal speciesScientia Hortic. 260, 108854. doi: 10.1016/j.scienta.2019.108854 [CrossRef[]
  • Kaensaksiri T., Soontornchainaksaeng P., Soonthornchareonnon N., Prathanturarug S. (2011). In vitro induction of polyploidy in Centella asiatica (L.) UrbanPlant Cell Tissue Organ Cult. 107, 187–194. doi:  10.1007/s11240-011-9969-8 [CrossRef[]
  • Karche T., Singh M. R. (2019). The application of hemp (Cannabis sativa L.) for a green economy: a reviewTurkish J. Bot. 43, 6. doi:  10.3906/bot-1907-15 [CrossRef[]
  • Khaleel C., Tabanca N., Buchbauer G. (2018). α-Terpineol, a natural monoterpene: A review of its biological propertiesOpen Chem. 16, 349–361. doi:  10.1515/CHEM-2018-0040/ASSET/GRAPHIC/J_CHEM-2018-0040_FIG_002.JPG [CrossRef[]
  • Krüger M., van Eeden T., Beswa D. (2022). Cannabis sativa cannabinoids as functional ingredients in snack foods-historical and developmental aspectsPlants (Basel). 11 (23), 3330. doi:  10.3390/plants11233330 [PMC free article] [PubMed] [CrossRef[]
  • Kulkarni M., Sawant N., Kolapkar A., Huprikar A., Desai N. (2021). Borneol: a promising monoterpenoid in enhancing drug delivery across various physiological barriersAAPS PharmSciTech 22, 145. doi:  10.1208/S12249-021-01999-8 [PubMed] [CrossRef[]
  • Kurtz L. E., Brand M. H., Lubell-Brand J. D. (2020). Production of tetraploid and triploid hempHortScience 55, 1703–1707. doi:  10.21273/HORTSCI15303-20 [CrossRef[]
  • Lavania U. C. (2005). Genomic and ploidy manipulation for enhanced production of phyto-pharmaceuticalsPlant Genet. Resour. 3, 170–177. doi:  10.1079/PGR200576 [CrossRef[]
  • Lavania U. C. (2013). Polyploidy, body size, and opportunities for genetic enhancement and fixation of heterozygosity in plantsNucleus (India) 56, 1–6. doi:  10.1007/S13237-013-0075-7 [CrossRef[]
  • Laverty K. U., Stout J. M., Sullivan M. J., Shah H., Gill N., Holbrook L., et al.. (2019). A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase lociGenome Res. 29, 146–156. doi:  10.1101/GR.242594.118 [PMC free article] [PubMed] [CrossRef[]
  • Lin X., Zhou Y., Zhang J., Lu X., Zhang F., Shen Q., et al.. (2011). Enhancement of artemisinin content in tetraploid Artemisia annua plants by modulating the expression of genes in artemisinin biosynthetic pathwayBiotechnol. Appl. Biochem. 58, 50–57. doi:  10.1002/BAB.13 [PubMed] [CrossRef[]
  • Madani H., Escrich A., Hosseini B., Sanchez-Muñoz R., Khojasteh A., Palazon J. (2021). Effect of polyploidy induction on natural metabolite production in medicinal plantsBiomolecules 11, 899. doi:  10.3390/BIOM11060899 [PMC free article] [PubMed] [CrossRef[]
  • Magagnini G., Grassi G., Kotiranta S. (2018). The effect of light spectrum on the morphology and cannabinoid content of Cannabis sativa LMed. Cannabis Cannabinoids. 1 (1), 19–27. doi:  10.1159/000489030 [PMC free article] [PubMed] [CrossRef[]
  • Mansouri H., Bagheri M. (2017). “Induction of polyploidy and its effect on Cannabis sativa L.” in Cannabis sativa L. – Botany and biotechnology. Eds. (Cham: Springer; ), 365–383. doi:  10.1007/978-3-319-54564-6_17 [CrossRef[]
  • Mishra B. K., Pathak S., Sharma A., Trivedi P. K., Shukla S. (2010). Modulated gene expression in newly synthesized auto-tetraploid of Papaver somniferum LSouth Afr. J. Bot. 76, 447–452. doi:  10.1016/J.SAJB.2010.02.090 [CrossRef[]
  • Motosugi H., Okudo K., Kataoka D., Naruo T. (2002). Comparison of growth characteristics between diploid and colchicine-induced tetraploid grape rootstocksJ. Japanese Soc. Hortic. Sci. 71, 335–341. doi:  10.2503/JJSHS.71.335 [CrossRef[]
  • Parsons J. L., Martin S. L., James T., Golenia G., Boudko E. A., Hepworth S. R. (2019). Polyploidization for the genetic improvement of Cannabis sativa Front. Plant Sci. 10. doi:  10.3389/fpls.2019.00476 [PMC free article] [PubMed] [CrossRef[]
  • Pollastro F., Minassi A., Fresu L. G. (2018). Cannabis phenolics and their bioactivitiesCurr. Med. Chem. 25 (10), 1160–1185. doi:  10.2174/0929867324666170810164636 [PubMed] [CrossRef[]
  • Radwan M. M., Chandra S., Gul S., Elsohly M. A. (2021). Cannabinoids, phenolics, terpenes and alkaloids of Cannabis Molecules 26, 2774. doi:  10.3390/MOLECULES26092774 [PMC free article] [PubMed] [CrossRef[]
  • Razazan A., Karunakar P., Mishra S. P., Sharma S., Miller B., Jain S., et al.. (2021). Activation of microbiota sensing – free fatty acid receptor 2 signaling ameliorates amyloid-β Induced neurotoxicity by modulating proteolysis-senescence axisFront. Aging Neurosci. 13. doi:  10.3389/FNAGI.2021.735933/BIBTEX [PMC free article] [PubMed] [CrossRef[]
  • Russo E. B. (2011). Taming THC: potential cannabis synergy and phytocannabinoid-terpenoid entourage effectsBr. J. Pharmacol. 163, 1344. doi:  10.1111/J.1476-5381.2011.01238.X [PMC free article] [PubMed] [CrossRef[]
  • Sattler M. C., Carvalho C. R., Clarindo W. R. (2016). The polyploidy and its key role in plant breedingPlanta 243, 281–296. doi:  10.1007/S00425-015-2450-X [PubMed] [CrossRef[]
  • Shimadzu Scientific Instruments (2016) Simplified cannabis terpene profiling by GCMS. Available at: https://www.ssi.shimadzu.com/sites/ssi.shimadzu.com/files/pim/pim_document_file/ssi/others/14371/GCMS-1604-TerpeneProfilingCannabis.pdf.
  • Silvarolla M. B., Mazzafera P., Lima M. M. A., de, Medina Filho H. P., Fazuoli L. C. (1999). Ploidy level and caffeine content in leaves of CoffeaSci. Agric. 56, 661–663. doi:  10.1590/S0103-90161999000300021 [CrossRef[]
  • Simmonds N. (1980). Polyploidy in plant breedingSpan 23 (2), 73–75. []
  • Small E. (2018). Dwarf germplasm: the key to giant Cannabis hempseed and cannabinoid cropsGenet. Resour Crop Evol. 65, 1071–1107. doi:  10.1007/s10722-017-0597-y [CrossRef[]
  • Song Q., Chen J. Z. (2015). Epigenetic and developmental regulation in plant polyploidsCurr. Opin. Plant Biol. 24, 101–109. doi:  10.1016/J.PBI.2015.02.007 [PMC free article] [PubMed] [CrossRef[]
  • Tan F. Q., Tu H., Wang R., Wu X. M., Xie K. D., Chen J. J., et al.. (2017). Metabolic adaptation following genome doubling in citrus doubled diploids revealed by non-targeted metabolomicsMetabolomics 13, 143. doi:  10.1007/S11306-017-1276-X [CrossRef[]
  • Trojak-Goluch A., Kawka-Lipińska M., Wielgusz K., Praczyk M. (2021). Polyploidy in industrial crops: applications and perspectives in plant breedingAgron. 2021 11, 2574. doi:  10.3390/AGRONOMY11122574 [CrossRef[]
  • Trojak-Goluch A., Skomra U. (2013). Artificially induced polyploidization in Humulus lupulus L. and its effect on morphological and chemical traitsBreed Sci. 63 (4), 393–399. doi:  10.1270/jsbbs.63.393 [PMC free article] [PubMed] [CrossRef[]
  • Urwin N. A. R. (2014). Generation and characterisation of colchicine-induced polyploid Lavandula x intermedia Euphytica 197, 331–339. doi:  10.1007/S10681-014-1069-5 [CrossRef[]
  • Wendel J. F., Lisch D., Hu G., Mason A. S. (2018). The long and short of doubling down: polyploidy, epigenetics, and the temporal dynamics of genome fractionationCurr. Opin. Genet. Dev. 49, 1–7. doi:  10.1016/J.GDE.2018.01.004 [PubMed] [CrossRef[]
  • Wohlmuth H., Leach D. N., Smith M. K., Myers S. P. (2005). Gingerol content of diploid and tetraploid clones of ginger (Zingiber officinale Roscoe)J. Agric. Food Chem. 53, 5772–5778. doi:  10.1021/JF050435B [PubMed] [CrossRef[]
  • Woo J., Yang H., Yoon M., Gadhe C. G., Pae A. N., Cho S., et al.. (2019). 3-carene, a phytoncide from pine tree has a sleep-enhancing effect by targeting the GABAA-benzodiazepine receptorsExp. Neurobiol. 28, 593. doi:  10.5607/EN.2019.28.5.593 [PMC free article] [PubMed] [CrossRef[]
  • Wood T. E., Takebayashi N., Barker M. S., Mayrose I., Greenspoon P. B., Rieseberg L. H. (2009). The frequency of polyploid speciation in vascular plantsProc. Natl. Acad. Sci. U.S.A. 106, 13875–13879. doi:  10.1073/PNAS.0811575106 [PMC free article] [PubMed] [CrossRef[]
  • Xu C., Tang Y., Chen R., Liang C., Liu X., Wu C., et al.. (2014). A comparative study of bioactive secondary metabolite production in diploid and tetraploid Echinacea purpurea (L.) MoenchPlant Cell Tissue Organ Cult. 116, 323–332. doi:  10.1007/s11240-013-0406-z [CrossRef[]
  • Yang X., Yang J., Gu X., Tao Y., Ji H., Miao X., et al.. (2023). (-)-Guaiol triggers immunogenic cell death and inhibits tumor growth in non-small cell lung cancerMol. Cell Biochem. 478, 1611–1620. doi:  10.1007/S11010-022-04613-Y [PMC free article] [PubMed] [CrossRef[]
  • Yeo S. K., Ali A. Y., Hayward O. A., Turnham D., Jackson T., Bowen I. D., et al.. (2016). β-bisabolene, a sesquiterpene from the essential oil extract of opoponax (Commiphora guidottii), exhibits cytotoxicity in breast cancer cell linesPhytother. Res. 30, 418–425. doi:  10.1002/PTR.5543 [PubMed] [CrossRef[]
  • Zonneveld B. J. M., Van Iren F. (2001). Genome size and pollen viability as taxonomic criteria: Application to the genus Hosta Plant Biol. 3, 176–185. doi:  10.1055/S-2001-12900 [CrossRef[]

Articles from Frontiers in Plant Science are provided here courtesy of Frontiers Media SA

Leave a Reply